Abstract: The effect of strain in graphene is usually modeled by a pseudomagnetic vector potential which is, however, derived in the limit of small strain. In realistic cases deviations are expected in view of graphene's very high strain tolerance, which can be up to 25%. Here we investigate the pseudomagnetic field generated by a Gaussian bump and we show that it exhibits significant differences with numerical tight-binding results. Furthermore, we calculate the electronic states in the strained region for a hexagon shaped flake with armchair edges. We find that the sixfold symmetry of the wave functions inside the Gaussian bump is directly related to the different effects of strain along the fundamental directions of graphene: zigzag and armchair. Low energy electrons are strongly confined in the armchair directions and are localized on the carbon atoms of a single sublattice.

Abstract: The electronic tunneling properties of the two stable forms of trilayer graphene (TLG), rhombohedral ABC and Bernal ABA, are examined for p-n and p-n-p junctions as realized by using a single gate (SG) or a double gate (DG). For the rhombohedral form, due to the chirality of the electrons, the Klein paradox is found at normal incidence for SG devices, while at high-energy interband scattering between additional propagation modes can occur. The electrons in Bernal ABA TLG can have a monolayer- or bilayer-like character when incident on a SG device. Using a DG, however, both propagation modes will couple by breaking the mirror symmetry of the system, which induces intermode scattering and resonances that depend on the width of the DG p-n-p junction. For ABC TLG the DG opens up a band gap which suppresses Klein tunneling. The DG induces also an unexpected asymmetry in the tunneling angle for single-valley electrons.

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: Using multicomponent Ginzburg-Landau simulations, we show a plethora of vortex states possible in mesoscopic three-band superconductors. We find that mesoscopic confinement stabilizes chiral states, with nontrivial phase differences between the band condensates, as the ground state of the system. As a consequence, we report the broken-symmetry vortex states, the chiral states where vortex cores in different band condensates do not coincide (split-core vortices), as well as fractional-flux vortex states with broken time-reversal symmetry.

Abstract: Flux entry in type-II superconductors without prominent symmetry is a first-order phase transition, where flux enters conventionally gradual in units of a flux quantum. Here we show that neither is necessarily the case in a mesoscopic superconducting sphere with a perforation. In axially applied magnetic field, vortices initially occupy the hole, and can oppose further flux entry in the sample. As a result, multiple-quanta flux entry is found at significantly higher field, and it can manifest as a second-order transition due to suppressed geometric barrier at the equatorial belt of the sample. At high fields a new state is found, with gradually destroyed condensate from the equator inwards, the exact opposite of surface superconductivity.

Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

Abstract: The effect of a parallel magnetic field on the orbital motion of electrons in high-quality superconducting nanowires resulting in a superconductor-to-normal transition which occurs through a cascade of jumps in the order parameter as a function of the magnetic field. Such cascades originate from the transverse size quantization that splits the conduction band into a series of subbands. Here, based on a numerical solution of the Bogoliubov-de Gennes equations for a hollow nanocylinder, we investigate how the quantum-size cascades depend on the confining geometry, i.e., by changing the cylinder radius R and its thickness d we cover the range from the nanowire-like to the nanofilm-like regime. The cascades are shown to become much less pronounced when increasing R/d, i.e., when the nanofilm-like regime is approached. When the temperature is non-zero they are thermally smoothed. This includes the spin-magnetic-field interaction which reduces the critical (depairing) parallel magnetic field H-c,H-parallel to but does not have any qualitative effect on the quantum cascades. From our calculations it is seen that the paramagnetic limiting field H-par significantly exceeds H-c,H-parallel to even in extremely narrow nanocylinders, i.e., when R, d are down to a few nanometers, and H-c,H-parallel to is only about 10% larger when switching-off the spin-magnetic-field interaction in this case. Both characteristic fields, H-c,H-parallel to and H-par, exhibit pronounced quantum-size oscillations. We demonstrate that the quantum cascades and the quantum-size oscillations survive in the presence of surface roughness.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: We investigate the dependence of the structural phase transitions in an infinite quasi-one-dimensional system of repulsively interacting particles on the profile of the confining channel. Three different functional expressions for the confinement potential related to real experimental systems are used that can be tuned continuously from a parabolic to a hard-wall potential in order to find a thorough understanding of the ordering of the chainlike structure transitions. We resolve the long-standing issue why the most theories predicted a 1-2-4-3-4 sequence of chain configurations with increasing density, while some experiments found the 1-2-3-4 sequence.

Abstract: Zigzag graphene nanoribbons patterned on graphane are studied using spin-polarized ab initio calculations. We found that the electronic and magnetic properties of the graphene/graphane superlattice strongly depends on the degree of hydrogenation at the interfaces between the two materials. When both zigzag interfaces are fully hydrogenated, the superlattice behaves like a freestanding zigzag graphene nanoribbon, and the magnetic ground state is antiferromagnetic. When one of the interfaces is half hydrogenated, the magnetic ground state becomes ferromagnetic, and the system is very close to being a half metal with possible spintronics applications whereas the magnetic ground state of the superlattice with both interfaces half hydrogenated is again antiferromagnetic. In this last case, both edges of the graphane nanoribbon also contribute to the total magnetization of the system. All the spin-polarized ground states are semiconducting, independent of the degree of hydrogenation of the interfaces. The ab initio results are supplemented by a simple tight-binding analysis that captures the main qualitative features. Our ab initio results show that patterned hydrogenation of graphene is a promising way to obtain stable graphene nanoribbons with interesting technological applications.

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: Surface effects of a doped thin film made of a strongly correlated material are investigated both in the absence and presence of a perpendicular electric field. We use an inhomogeneous Gutzwiller approximation for a single-band Hubbard model in order to describe correlation effects. For low doping, the bulk value of the quasiparticle weight is recovered exponentially deep into the slab, but with increasing doping, additional Friedel oscillations appear near the surface. We show that the inverse correlation length has a power-law dependence on the doping level. In the presence of an electrical field, considerable changes in the quasiparticle weight can be realized throughout the system. We observe a large difference (as large as five orders of magnitude) in the quasiparticle weight near the opposite sides of the slab. This effect can be significant in switching devices that use the surface states for transport. DOI: 10.1103/PhysRevB.87.035131

Abstract: The electrical transport properties of a MoGe thin film with a honeycomb array of nanoscale holes are investigated. The critical current of the system shows nonmatching anomalies as a function of applied magnetic field, enabling us to distinguish between multiquanta vortices trapped in the holes and interstitial vortices located between the holes. The number of vortices trapped in each hole is found to be larger than the saturation number predicted for an isolated hole and shows a nonlinear field dependence, leading to the caging effect as predicted from the Ginzburg-Landau (GL) theory. Our experimental results are supplemented by numerical simulations based on the GL theory.

Abstract: The electronic structure of GaAs/AlxGa1−xAs and AlxGa1−xAs/GaAs core-shell nanowires grown in the [001] direction is studied. The k⋅p method with the 6×6 Kohn-Lüttinger Hamiltonian, taking into account the split-off band is used. The variation in the energy level dispersion, the spinor contribution to the ground state and the optical interband absorption are studied. For some range of parameters the top of the valence band exhibits a camelback structure which results in an extra peak in the optical absorption.

Abstract: Based on the Gor'kov formalism for a clean s-wave superconductor, we develop an extended version of the single-band Ginzburg-Landau (GL) theory by means of a systematic expansion in the deviation from the critical temperature T(c), i.e., tau = 1 – T/T(c). We calculate different contributions to the order parameter and the magnetic field: the leading contributions (proportional to tau(1/2) in the order parameter and. t in the magnetic field) are controlled by the standard GL theory, while the next-to-leading terms (proportional to tau(3/2) in the gap and proportional to tau(2) in the magnetic field) constitute the extended GL (EGL) approach. We derive the free-energy functional for the extended formalism and the corresponding expression for the current density. To illustrate the usefulness of our formalism, we calculate, in a semianalytical form, the temperature-dependent correction to the GL parameter at which the surface energy becomes zero, and analytically, the temperature dependence of the thermodynamic critical field. We demonstrate that the EGL formalism is not just a mathematical extension to the theory: variations of both the gap and the thermodynamic critical field with temperature calculated within the EGL theory are found in very good agreement with the full BCS results down to low temperatures, which dramatically improves the applicability of the formalism compared to its standard predecessor.