Abstract: We investigate the magnetic interface states of graphene quantum dots that contain p-n junctions. Within a tight-binding approach, we consider rectangular quantum dots in the presence of a perpendicular magnetic field containing p-n as well as p-n-p and n-p-n junctions. The results show the interplay between the edge states associated with the zigzag terminations of the sample and the snake states that arise at the p-n junction due to the overlap between electron and hole states at the potential interface. Remarkable localized states are found at the crossing of the p-n junction with the zigzag edge having a dumb-bell-shaped electron distribution. The results are presented as a function of the junction parameters and the applied magnetic flux. DOI: 10.1103/PhysRevB.87.035426

Abstract: At the interface of electrostatic potential kink profiles, one-dimensional chiral states are found in bilayer graphene (BLG). Such structures can be created by applying an asymmetric potential to the upper and the lower layers of BLG. We found the following: (i) due to the strong confinement by the single kink profile, the unidirectional states are only weakly affected by a magnetic field; (ii) increasing the smoothness of the kink potential results in additional bound states, which are topologically different from those chiral states; and (iii) in the presence of a kink-antikink potential, the overlap between the oppositely moving chiral states results in the appearance of crossing and anticrossing points in the energy spectrum. This leads to the opening of tunable minigaps in the spectrum of the unidirectional topological states.

Abstract: Artificial molecular states of double quantum dots defined in bilayer graphene are studied with the atomistic tight-binding method and its low-energy continuum approximation. We indicate that the extended electron wave functions have opposite parities on sublattices of the layers and that the ground-state wave-function components change from bonding to antibonding with the interdot distance. In the weak-coupling limit, the one most relevant for quantum dots defined electrostatically, the signatures of the interdot coupling include, for the two-electron ground state, formation of states with symmetric or antisymmetric spatial wave functions split by the exchange energy. In the high-energy part of the spectrum the states with both electrons in the same dot are found with the splitting of energy levels corresponding to simultaneous tunneling of the electron pair from one dot to the other.

Abstract: We numerically investigate the spin-orbit (SO) coupling effect on the magnetic flux evolution of energy and supercurrent in mesoscopic d-wave superconducting loops by solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. It is found that the energy spectrum splits when the SO interaction is involved and the Majorana zero mode can be realized in the [100] edges of square systems for an appropriate SO coupling strength. Superconducting phase transitions appear when the energy gap closes, accompanied by energy jumps between different energy parabolas in the ground state, which provides a possible mechanism to support fractional flux periodicity of supercurrent. Moreover, in the case of rectangular loops with SO coupling, the jumps of the ground-state energy gradually disappear by increasing the ratio of length to height of the sample, and a paramagnetic response with opposite direction of the screening current around zero flux value can occur in such systems.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.

Abstract: We investigated the magnetic flux dependence of the supercurrent in mesoscopic d-wave superconducting loops, containing impurities and surface defects, by numerically solving the Bogoliubovde Gennes equations self-consistently. In the presence of impurities, bound states arise close to the Fermi energy. In the case of a single impurity, the flux-induced current is found to be suppressed. This can be different when more impurities are introduced in the sample due to the quantum interference effect, which depends sensitively on the relative position between the impurities. We further analyze the effect of small surface defects at the inner or outer edge of the loop, and show that indentation and bulge defects have pronounced and different effects on the supercurrent.

Abstract: We investigate the effect of ellipticity ratio of the polarized radiation field on optoelectronic properties of a two-dimensional (2D) semi-Dirac (SD) system. The optical conductivity is calculated within the energy balance equation approach derived from the semiclassical Boltzmann equation. We find that there exists the anisotropic optical absorption induced via both the intra-and interband electronic transition channels in the perpendicular xx and yy directions. Furthermore, we examine the effects of the ellipticity ratio, the temperature, the carrier density, and the band-gap parameter on the optical conductivity of the 2D SD system placed in transverse and vertical directions, respectively. It is shown that the ellipticity ratio, temperature, carrier density, and band-gap parameter can play the important roles in tuning the strength, peak position, and shape of the optical conductivity spectrum. The results obtained from this study indicate that the 2D SD system can be a promising anisotropic and tunable optical and optoelectronic material for applications in innovative 2D optical and optoelectronic devices, which are active in the infrared and terahertz bandwidths.

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: Recent experiments show that an atomic step on the surface of atomically thin metallic films can strongly affect electronic transport. Here we reveal multiple and versatile effects that such a surface step can have on superconductivity in ultrathin films. By solving the Bogoliubov-de Gennes equations self-consistently in this regime, where quantum confinement dominates the emergent physics, we show that the electronic structure is profoundly modified on the two sides of the step, as is the spatial distribution of the superconducting order parameter and its dependence on temperature and electronic gating. Furthermore, the surface step changes nontrivially the transport properties both in the proximity-induced superconducting pair correlations and the Josephson effect, depending on the step height. These results offer a new route to tailor superconducting circuits and design atomically thin heterojunctions made of one same material.

Abstract: Chiral p-wave superconductors in applied magnetic field can exhibit more complex topological defects than just conventional superconducting vortices, due to the two-component order parameter (OP) and the broken time-reversal symmetry. We investigate the electronic properties of those exotic states, some of which contain clusters of one-component vortices in chiral components of the OP and/or exhibit skyrmionic character in the relative OP space, all obtained as a self-consistent solution of the microscopic Bogoliubov-de Gennes equations. We reveal the link between the local density of states (LDOS) of the novel topological states and the behavior of the chiral domain wall between the OP components, enabling direct identification of those states in scanning tunneling microscopy. For example, a skyrmion always contains a closed chiral domain wall, which is found to be mapped exactly by zero-bias peaks in LDOS. Moreover, the LDOS exhibits electron-hole asymmetry, which is different from the LDOS of conventional vortex states with same vorticity. Finally, we present the magnetic field and temperature dependence of the properties of a skyrmion, indicating that this topological defect can be surprisingly large in size, and can be pinned by an artificially indented nonsuperconducting closed path in the sample. These features are expected to facilitate the experimental observation of skyrmionic states, thereby enabling experimental verification of chirality in emerging superconducting materials.

Abstract: We report characteristic vortex configurations in s + id superconductors with time-reversal symmetry breaking, exposed to magnetic field. A vortex in the s + id state tends to have an opposite phase winding between s- and d-wave condensates. We find that this peculiar feature together with the competition between s- and d-wave symmetry results in three distinct classes of vortical configurations. When either s or d condensate absolutely dominates, vortices form a conventional lattice. However, when one condensate is relatively dominant, vortices organize in chains that exhibit skyrmionic character, separating the chiral components of the s +/- id order parameter into domains within and outside the chain. Such skyrmionic chains are found stable even at high magnetic field. When s and d condensates have comparable strength, vortices split cores in two chiral components to form full-fledged skyrmions, i.e., coreless topological structures with an integer topological charge, organized in a lattice. We provide characteristic magnetic field distributions of all states, enabling their identification in, e.g., scanning Hall probe and scanning SQUID experiments. These unique vortex states are relevant for high-T-c cuprate and iron-based superconductors, where the relative strength of competing pairing symmetries is expected to be tuned by temperature and/or doping level, and can help distinguish s + is and s + id superconducting phases.

Abstract: We study the mobility of Dirac fermions in monolayer graphene on a GaAs substrate, limited by the combined action of the extrinsic potential of piezoelectric surface acoustical phonons of GaAs (PA) and of the intrinsic deformation potential of acoustical phonons in graphene (DA). In the high-temperature (T) regime, the momentum relaxation rate exhibits the same linear dependence on T but different dependencies on the carrier density n, corresponding to the mobility mu proportional to 1 root n and 1/n, respectively for the PA and DA scattering mechanisms. In the low-T Bloch-Gruneisen regime, the mobility shows the same square-root density dependence mu proportional to root n, but different temperature dependencies mu proportional to T-3 and T-4, respectively for PA and DA phonon scattering. DOI: 10.1103/PhysRevB.87.075443

Abstract: We present a theoretical study of the plasmon and coupled plasmon-phonon modes induced by intraband electron-electron interaction in graphene in the presence of driving dc electric field. We find that the electric field dependence of these collective excitation modes in graphene differs significantly from that in a conventional two-dimensional electron gas with a parabolic energy spectrum. This is due mainly to the fact that graphene has a linear energy spectrum and the Fermi velocity of electrons in graphene is much larger than the drift velocity of electrons. The obtained results demonstrate that the plasmon and coupled plasmon-phonon modes in graphene can be tuned by applying not only the gate voltage but also the source-to-drain field. The manipulation of plasmon and coupled plasmon-phonon modes by source-to-drain voltage can let graphene be more conveniently applied as an advanced plasmonic material.

Abstract: We analyze theoretically and experimentally vortex configurations in mesoscopic superconducting squares. Our theoretical approach is based on the analytical solution of the London equation using Green's-function method. The potential-energy landscape found for each vortex configuration is then used in Langevin-type molecular-dynamics simulations to obtain stable vortex configurations. Metastable states and transitions between them and the ground state are analyzed. We present our results of the first direct visualization of vortex patterns in micrometer-sized Nb squares, using the Bitter decoration technique. We show that the filling rules for vortices in squares with increasing applied magnetic field can be formulated, although in a different manner than in disks, in terms of formation of vortex “shells”.