Abstract: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: We have calculated the plasmon modes in graphene double layer structures at finite temperatures, taking into account the inhomogeneity of the dielectric background of the system. The effective dielectric function is obtained from the solution of the Poisson equation of a three-layer dielectric medium with graphene sheets located at the interfaces, separating the different materials. Due to the momentum dispersion of the effective dielectric function, the intra- and interlayer bare Coulomb interactions in the graphene double layer system acquires an additional momentum dependence-an effect that is of the order of the interlayer interaction itself. We show that the energies of the in-phase and out-of-phase plasmon modes are determined largely by different values of the spatially dependent effective dielectric function. The effect of the dielectric inhomogeneity increases with temperature, and even at high temperatures the energy shift induced by the dielectric inhomogeneity and temperature itself remains larger than the broadening of the plasmon energy dispersions due to the Landau damping. The obtained new features of the plasmon dispersions can be observed in frictional drag measurements and in inelastic light scattering and electron energy-loss spectroscopies.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: Quantum confinement of the perpendicular motion of electrons in single-crystalline metallic superconducting nanofilms splits the conduction band into a series of single-electron subbands. A distinctive feature of such a nanoscale multi-band superconductor is that the energetic position of each subband can vary significantly with changing nanofilm thickness, substrate material, protective cover and other details of the fabrication process. It can occur that the bottom of one of the available subbands is situated in the vicinity of the Fermi level. We demonstrate that the character of the superconducting pairing in such a subband changes dramatically and exhibits a clear molecule-like trend, which is very similar to the well-known crossover from the Bardeen-Cooper-Schrieffer regime to Bose-Einstein condensation (BCS-BEC) observed in trapped ultracold fermions. For Pb nanofilms with thicknesses of 4 and 5 monolayers (MLs) this will lead to a spectacular scenario: up to half of all the Cooper pairs nearly collapse, shrinking in the lateral size (parallel to the nanofilm) down to a few nanometers. As a result, the superconducting condensate will be a coherent mixture of almost molecule-like fermionic pairs with ordinary, extended Cooper pairs.

Abstract: Using highly efficient simulations of the tight-binding Bogoliubov-de-Gennes model, we solved self-consistently for the pair correlation and the Josephson current in a superconducting-bilayer graphene-superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short- and long-junction regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the long-junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if the junction length is larger than the Fermi length.

Abstract: Using the valence force field model of Perebeinos and Tersoff (2009 Phys. Rev. B 79 241409(R)), different energy modes of suspended graphene subjected to tensile or compressive strain are studied. By carrying out Monte Carlo simulations it is found that: (i) only for small strains (vertical bar epsilon vertical bar (sic) 0.02) is the total energy symmetrical in the strain, while it behaves completely differently beyond this threshold; (ii) the important energy contributions in stretching experiments are stretching, angle bending, an out-of-plane term, and a term that provides repulsion against pi-pi misalignment; (iii) in compressing experiments the two latter terms increase rapidly, and beyond the buckling transition stretching and bending energies are found to be constant; (iv) from stretching-compressing simulations we calculated the Young's modulus at room temperature 350 +/- 3.15 N m(-1), which is in good agreement with experimental results (340 +/- 50 N m(-1)) and with ab initio results (322-353) N m(-1); (v) molar heat capacity is estimated to be 24.64 J mol(-1) K-1 which is comparable with the Dulong-Petit value, i. e. 24.94 J mol(-1) K-1, and is almost independent of the strain; (vi) nonlinear scaling properties are obtained from height-height correlations at finite temperature; (vii) the used valence force field model results in a temperature independent bending modulus for graphene, and (viii) the Gruneisen parameter is estimated to be 0.64.

Abstract: The effective Hamiltonian for an electron in a quasi-one-dimensional quantum ring in the presence of spin-orbit interactions is derived. We demonstrate that, when both coupling types are simultaneously present, the effective Hamiltonian derived by the lowest-radial-state approximation produces energy spectra and charge densities which deviate strongly from the exact ones. For equal Rashba and Dresselhaus coupling constants the lowest-radial-state approximation opens artifactal avoided crossings in the energy spectra and deforms the circular symmetry of the confined charge densities. In this case, there does not exist a ring thin enough to justify the restriction to the lowest radially quantized energy state. We derive the effective Hamiltonian accounting for both the lowest and the first excited radial states, and show that the inclusion of the latter restores the correct features of the exact solution. Relation of this result to the states of a quantum wire is also discussed.

Abstract: The Hall (R-H) and bend (R-B) resistances of a graphene Hall bar structure containing a pn-junction are calculated when in the ballistic regime. The simulations are done using the billiard model. Introducing a pn-junction-dividing the Hall bar geometry in two regions-leads to two distinct regimes exhibiting very different physics: (1) both regions are of n-type and (2) one region is n-type and the other p-type. In regime (1), a “Hall plateau”-an enhancement of the resistance-appears for R-H. On the other hand, in regime (2), we found a negative R-H, which approaches zero for large B. The bend resistance is highly asymmetric in regime (2) and the resistance increases with increasing magnetic field B in one direction while it reduces to zero in the other direction. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4704667]

Abstract: We derive the extended Ginzburg-Landau (GL) formalism for a clean s-wave two-band superconductor by employing a systematic expansion of the free-energy functional and the corresponding matrix gap equation in powers of the small deviation from the critical temperature tau = 1 – T/T-c. The two lowest orders of this expansion produce the equation for T-c and the standard GL theory. It is shown that in agreement with previous studies, this two-band GL theory maps onto the single-band GL model and thus fails to describe the difference in the spatial profiles of the two-band condensates. We prove that this difference appears already in the leading correction to the standard GL theory, which constitutes the extended GL formalism. We derive linear differential equations that determine the leading corrections to the band order parameters and magnetic field, discuss the validity of these equations, and consider examples of an important interplay between the band condensates. Finally, we present numerical results for the thermodynamic critical magnetic field and temperature-dependent band gaps for recent materials of interest, which are in very good agreement with those obtained from the full BCS approach in a wide temperature range. To this end, we emphasize the advantages of our extended GL theory in comparison with the often used two-component GL-like model based on an unreconstructed two-band generalization of the Gor'kov derivation.

Abstract: Superconductivity in nanowires made of weak coupling superconductor materials is investigated using a Green function approach. We show that these are multigap systems in which the ratio Delta(T)/k(B)T(c) is to a large extent similar to what is observed in some high-T-c two-gap systems, such as MgB2 and some of the Fe-based superconductors. On the other hand, because of confinement, the superfluid density has a temperature behavior of the form n(s)(T) = 1 – (T/T-c)(3) near T-c, thus deviating from the BCS behavior for bulk superconductors.

Abstract: The electronic structure of GaAs, InAs and InSb nanowires is studied using the six-band and the eight-band k.p models. The effect of the different Luttinger-like parameters (in the eight-band model) on the hole band structure is investigated. Although GaAs nanostructures are often treated within a six-band model because of the large bandgap, it is shown that an eight-band model is necessary for a correct description of its hole spectrum. The camel-back structure usually found in the six-band model is not always present in the eight-band model. This camel-back structure depends on the interaction between light and heavy holes, especially the ones with opposite spin. The latter effect is less pronounced in an eight-band model, but could be very sensitive to the Kane inter-band energy (E-P) value.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Spin selectivity in angle-resolved Auger photoelectron coincidence spectroscopy (AR-APECS) is used to probe electron correlation in ferromagnetic thin films. In particular, exploiting the AR-APECS capability to discriminate Auger electron emission events characterized by valence hole pairs created either in the high or in the low total spin state, a strong correlation effect in the Fe M2,3VV Auger line shape (measured in coincidence with the Fe 3p photoelectrons) of Fe/Cu(001) thin films is detected and ascribed to interactions within the majority spin subband. Such an assignment follows from a close comparison of the experimental AR-APECS line shapes with the predictions of a model based on spin polarized density functional theory and the Cini-Sawatzky approach.

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: In 1956 Cooper demonstrated (1956 Phys. Rev. 104 1189) that, no matter how weak the attraction is, two electrons in three-dimensional (3D) space just above the Fermi sea could be bound. In this work we investigate the influence of confinement on the binding energy of a Cooper pair. We show that confinement-induced modification of the Fermi sea results in a significant increase of the binding energy, when the bottom of an energy subband is very close to the Fermi surface.