Abstract: Superconducting correlations in an isolated metallic grain are governed by the interplay between two energy scales: the mean level spacing delta and the bulk pairing gap Delta(0), which are strongly influenced by the position of the chemical potential with respect to the closest single-electron level. In turn superconducting correlations affect the position of the chemical potential. Within the parity projected BCS model we investigate the probability distribution of the chemical potential in a superconducting grain with randomly distributed single-electron levels. Taking into account statistical fluctuations of the chemical potential due to the pairing interaction, we find that such fluctuations have a significant impact on the critical level spacing delta(c) at which the superconducting correlations cease: the critical ratio delta(c)/Delta(0) at which superconductivity disappears is found to be increased.

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: Experiments have confirmed that double monolayer graphene does not generate finite-temperature electron-hole superfluidity, because of very strong screening of the pairing attraction. The linear dispersing energy bands in monolayer graphene block any attempt to reduce the strength of the screening. We propose a hybrid device with two sheets of monolayer graphene in a modulated periodic perpendicular magnetic field. The field preserves the isotropic Dirac cones of the original monolayers but reduces the slope of the cones, making the monolayer Fermi velocity v(F) smaller. We demonstrate that with current experimental techniques, the reduction in vF can weaken the screening sufficiently to allow electron-hole superfluidity at measurable temperatures.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.

Abstract: The effect of trigonal warping on the transmission of electrons tunneling through potential barriers in graphene is investigated. We present calculations of the transmission coefficient for single and double barriers as a function of energy, incidence angle and barrier heights. The results show remarkable valley-dependent directional effects for barriers oriented parallel to the armchair or parallel to the zigzag direction. These results indicate that electrostatic gates can be used as valley filters in graphene-based devices.

Abstract: In this work we have investigated the structural and electronic properties of SnO, which is built up from layers kept together by van der Waals (vdW) forces. The combination of a vdW functional within density functional theory (DFT) and quasiparticle band structure calculations within the GW approximation provides accurate values for the lattice parameters, atomic positions, and the electronic band structure including the fundamental (indirect) and the optical (direct) band gap without the need of experimental or empirical input. A systematic comparison is made between different levels of self-consistency within the GW approach {following the scheme of Shishkin et al. [Phys. Rev. B 75, 235102 (2007)]} and the results are compared with DFT and hybrid functional results. Furthermore, the effect of the vdW-corrected functional as a starting point for the GW calculation of the band gap has been investigated. Finally, we studied the effect of the vdW functional on the electron charge density.

Abstract: The electronic and vibrational properties of germanane and fluorinated germanene are studied within density functional theory (DFT) and density functional perturbation theory frameworks. Different structural configurations of germanane and fluorinated germanene are investigated. The energy difference between the different configurations are consistently smaller than the energy of thermal fluctuations for all the analyzed DFT functionals LDA, GGA, and hybrid functionals, which implies that, in principle, it is possible to find these different configurations in different regions of the sample as minority phases or local defects. We calculate the Raman and infrared spectra for these configurations by using ab initio calculations and compare it with available experimental spectra for germanane. Our results show the presence of minority phases compatible with the configurations analyzed in this work. As these low energy configurations are metastable the present work shows that the synthesis of these energy competing phases is feasible by selectively changing the synthesis conditions, which is an opportunity to expand in this way the availability of new two-dimensional compounds.

Abstract: We investigate mesoscopic Josephson-junction arrays created by patterning superconducting disks on monolayer graphene, concentrating on the high-T/T-c regime of these devices and the phenomena which contribute to the superconducting glass state in diffusive arrays. We observe features in the magnetoconductance at rational fractions of flux quanta per array unit cell, which we attribute to the formation of flux-quantized vortices. The applied fields at which the features occur are well described by Ginzburg-Landau simulations that take into account the number of unit cells in the array. We find that the mean conductance and universal conductance fluctuations are both enhanced below the critical temperature and field of the superconductor, with greater enhancement away from the graphene Dirac point.

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”.

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.