Abstract: We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of a Rashba-like spin-orbit interaction. Our results corroborate previous work performed with the momentum average approximation and with weak-coupling perturbation theory. We find that spin-orbit coupling increases the effective mass in the regime with weak electron-phonon coupling, and decreases the effective mass in the regimes of intermediate and strong electron-phonon coupling. Analytical strong-coupling perturbation theory results confirm our numerical results in the small-polaron regime. A large amount of spin-orbit coupling can lead to a significant lowering of the polaron effective mass.

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

Abstract: The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various nonlinear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the fidelity susceptibility of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time scales as long as the ground-state fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating phase, a dissipative regime, or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

Abstract: A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples, or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single-atom vacancies, we observe a strong suppression of the supercurrent, which is a consequence of strong intervalley scattering. Although lattice deformations should not induce intervalley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudomagnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e., existence of electron-hole puddles, or finite. In both cases, short-range impurities strongly affect the supercurrent, similar to the vacancies scenario.

Abstract: Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information similarly to the way spintronics uses electron spin. We show how a double-barrier structure consisting of electric and vector potentials can be used to filter massive Dirac electrons based on their valley index. We study the resonant transmission through a finite number of barriers and we obtain the energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is included, the energy bands and energy gaps at the K and K′ points are different and they can be tuned by changing the potential.

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: The Tomasch effect (TE) is due to quasiparticle interference (QPI) as induced by a nonuniform superconducting order parameter, which results in oscillations in the density of states (DOS) at energies above the superconducting gap. Quantum confinement in nanoscale superconductors leads to an inhomogenerous distribution of the Cooperpair condensate, which, as we found, triggers the manifestation of a new TE. We investigate the electronic structure of nanoscale superconductors by solving the Bogoliubov-de Gennes (BdG) equations self-consistently and describe the TE determined by two types of processes, involving two-or three-subband QPIs. Both types of QPIs result in additional BCS-like Bogoliubov-quasiparticles and BCS-like energy gaps leading to oscillations in the DOS and modulated wave patterns in the local density of states. These effects are strongly related to the symmetries of the system. A reduced 4 x 4 inter-subband BdG Hamiltonian is established in order to describe analytically the TE of two-subband QPIs. Our study is relevant to nanoscale superconductors, either nanowires or thin films, Bose-Einsten condensates, and confined systems such as two-dimensional electron gas interface superconductivity.

Abstract: Bogoliubov-de Gennes theory is used to investigate the effect of the size of a superconducting square on the vortex states in the quantum confinement regime. When the superconducting coherence length is comparable to the Fermi wavelength, the shape resonances of the superconducting order parameter have strong influence on the vortex configuration. Several unconventional vortex states, including asymmetric ones, giant-multivortex combinations, and states comprising giant antivortices, were found as ground states and their stability was found to be very sensitive on the value of k(F)xi(0), the size of the sample W, and the magnetic flux Phi. By increasing the temperature and/or enlarging the size of the sample, quantum confinement is suppressed and the conventional mesoscopic vortex states as predicted by the Ginzburg-Laudau (GL) theory are recovered. However, contrary to the GL results we found that the states containing symmetry-induced vortex-antivortex pairs are stable over the whole temperature range. It turns out that the inhomogeneous order parameter induced by quantum confinement favors vortex-antivortex molecules, as well as giant vortices with a rich structure in the vortex core-unattainable in the GL domain.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.

Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

Abstract: Two-dimensional atomic crystals can radically change their properties in response to external influences, such as substrate orientation or strain, forming materials with novel electronic structure(1-5). An example is the creation of weakly dispersive, 'flat' bands in bilayer graphene for certain 'magic' angles of twist between the orientations of the two layers(6). The quenched kinetic energy in these flat bands promotes electron-electron interactions and facilitates the emergence of strongly correlated phases, such as superconductivity and correlated insulators. However, the very accurate fine-tuning required to obtain the magic angle in twisted-bilayer graphene poses challenges to fabrication and scalability. Here we present an alternative route to creating flat bands that does not involve fine-tuning. Using scanning tunnelling microscopy and spectroscopy, together with numerical simulations, we demonstrate that graphene monolayers placed on an atomically flat substrate can be forced to undergo a buckling transition(7-9), resulting in a periodically modulated pseudo-magnetic field(10-14), which in turn creates a 'post-graphene' material with flat electronic bands. When we introduce the Fermi level into these flat bands using electrostatic doping, we observe a pseudogap-like depletion in the density of states, which signals the emergence of a correlated state(15-17). This buckling of two-dimensional crystals offers a strategy for creating other superlattice systems and, in particular, for exploring interaction phenomena characteristic of flat bands. Buckled monolayer graphene superlattices are found to provide an alternative to twisted bilayer graphene for the study of flat bands and correlated states in a carbon-based material.

Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

Abstract: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

Abstract: The interaction of an injected electron towards a quantum dot (QD) containing a single confined electron is investigated using a flexible time-dependent quantum mechanics formalism, which allows both electrons to move and undergo quantum transitions. Different scenarios combining quantum dot dimensions, dielectric constant, injected wave packet energy, and width were explored, and our main results are: (i) due to the large characteristic transitions times between the confined state in the quantum dot and the delocalized state in the continuum, it is relatively difficult to ionize the occupied QD by Coulomb interaction solely and (ii) the charging state of the quantum dot can be sensed by direct injection of charges. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4759292]

Abstract: We perform the wavepacket dynamics simulation of a graphene-based device where propagating electron trajectories are tamed by an applied magnetic field toward a normal/superconductor interface. The magnetic field controls the incidence angle of the incoming electronic wavepacket at the interface, which results in the tunable electron-hole ratio in the reflected wave function due to the angular dependence of the Andreev reflection. Here, mapped control of the quasiparticle trajectories by the external magnetic field not only defines an experimental probe for fundamental studies of the Andreev reflection in graphene but also lays the foundation for further development of magnetic focusing devices based on nanoengineered superconducting two-dimensional materials.