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: 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: 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: 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: 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: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Abstract: The electron states in graphene-based magnetic dot and magnetic ring structures and combinations of both are investigated. The corresponding spectra are studied as a function of the radii, the strengths of the inhomogeneous magnetic field and of a uniform background field, the strength of an electrostatic barrier and the angular momentum quantum number. In the absence of an external magnetic field we have only long-lived quasi-bound and scattering states and we assess their influence on the density of states. In addition, we consider elastic electron scattering by a magnetic dot, whose average B vanishes, and show that the Hall and longitudinal resistivities, as a function of the Fermi energy, exhibit a pronounced oscillatory structure due to the presence of quasi-bound states. Depending on the dot parameters this oscillatory structure differs substantially for energies below and above the first Landau level.

Abstract: The scattering of two-dimensional (2D) massless Dirac electrons is investigated in the presence of a random array of circular mass barriers. The inverse momentum relaxation time and the Hall factor are calculated and used to obtain parallel and perpendicular resistivity components within linear transport theory. We found a nonzero perpendicular resistivity component which has opposite sign for electrons in the different K and K′ valleys. This property can be used for valley filter purposes. The total cross section for scattering on penetrable barriers exhibits resonances due to the presence of quasibound states in the barriers that show up as sharp gaps in the cross section while for Schrödinger electrons they appear as peaks.

Abstract: The conductance, the transmission, and the reflection probabilities through rectangular potential barriers and p-n junctions are obtained for bilayer graphene taking into account the four bands of the energy spectrum. We have evaluated the importance of the skew hopping parameters gamma(3) and gamma(4) to these properties and show that for energies E > gamma(1)/100 their effect is negligible. For high energies two modes of propagation exist and we investigate scattering between these modes. For perpendicular incidence both propagation modes are decoupled, and scattering between them is forbidden. This extends the concept of pseudospin as defined within the two-band approximation to a four-band model and corresponds to the (anti) symmetry of the wave functions under in-plane mirroring. New transmission resonances are found that appear as sharp peaks in the conductance which are absent in the two-band approximation. The application of an interlayer bias to the system (1) breaks the pseudospin structure, (2) opens a band gap that results in a distinct feature of suppressed transmission in the conductance, and (3) breaks the angular symmetry with respect to normal incidence in the transmission and reflection.

Abstract: We study the effect of a hydrogen atom on the thermal stability of a trapped CH3 molecule on graphene using ReaxFF molecular dynamics simulations. Due to the hydrogen-molecule interaction, enhanced pinning of the CH3 molecule is observed when it is positioned adjacent to the graphene site with the hydrogen atom. We discuss the formation process of such a stable configuration, which originates from different adhesion and migration energies of the hydrogen atom and the CH3 molecule. We also studied the effect of the CH3-H configuration on the electronic transport properties of graphene nanoribbons using first principles density-functional calculations. We found that the formation of the CH3-H structure results in extra features in the transmission spectrum due to the formation of strongly localized states, which are absent when the CH3 molecule is trapped on pristine graphene. Our findings will be useful in exploiting gas sensing properties of graphene, especially for selective detection of individual molecules. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: We investigate how the interplay of quantum confinement and particle number-parity fluctuations affects superconducting correlations in ultra-small metallic grains. Using the number-parity projected BCS formalism we calculate the critical temperature and the excitation gap as a function of the grain size for grains with even and odd number of confined carriers. We show that the experimentally observed anomalous increase of the coupling ratio ΔE/kBTc with decreasing superconducting grain size can be attributed to an enhancement of the number-parity fluctuations in ultra-small grains.

Abstract: An analytical approach, using the Dirac-Weyl equation, is implemented to obtain the energy spectrum and optical absorption of a circular graphene quantum dot in the presence of an external magnetic field. Results are obtained for the infinite-massand zigzag boundary conditions. We found that the energy spectrum of a dot with the zigzag boundary condition exhibits a zero-energy band regardless of the value of the magnetic field, while for the infinite-mass boundary condition, the zero-energy states appear only for high magnetic fields. The analytical results are compared to those obtained from the tight-binding model: (i) we show the validity range of the continuum model and (ii) we find that the continuum model with the infinite-mass boundary condition describes rather well its tight-binding analog, which can be partially attributed to the blurring of the mixed edges by the staggered potential.

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with a ringlike shape are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum dots have volcano-like shapes. The continuum mechanical model is employed to determine the strain distribution, and the single-band envelope function approach is adopted to compute the electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the same width. Comparison with the unstrained ringlike quantum dot shows that the amplitude of the excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured in recent experiments. DOI: 10.1103/PhysRevB.87.085314

Abstract: We study Dirac plasmons and their damping in spatially separated N-layer graphene structures at finite doping and temperatures. The plasmon spectrum consists of one optical excitation with square-root dispersion and N – 1 acoustical excitations with linear dispersion, which are undamped at zero temperature and finite doping within a triangular energy region outside the electron-hole continuum. In the long-wavelength limit the energy and weight of the optical plasmon modes increase, respectively, as the square root and linearly with N in agreement with recent experimental findings. The energy and weight of the upper-lying acoustical branches also increase with N. This increase is strongest for the uppermost acoustical mode, and we find that its energy can exceed at some value of momentum the plasmon energy in an individual graphene sheet. Meanwhile, the energy of the low-lying acoustical branches decreases weakly with N as compared with the single acoustical mode in double-layer graphene structures. Our numerical calculations provide a detailed understanding of the overall behavior of the wave-vector dependence of the optical and acoustical multilayer plasmon modes and show how their dispersion and damping are modified as a function of temperature, interlayer spacing, and inlayer carrier density in (un)balanced graphene multilayer structures. DOI: 10.1103/PhysRevB.87.085401