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: 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: We consider two possible mechanisms for single-vortex fluctuative entry/exit through the surface barrier in ultrathin superconducting disk-shaped nanoislands made of Pb and consisting of just a few monoatomic layers, which can be fabricated using modern techniques. We estimate tunneling probabilities and establish criteria for the crossover between these two mechanisms depending on magnetic field and system sizes. For the case of vortex entry, quantum tunneling dominates on the major part of the temperature/flux phase diagram. For the case of vortex exit, thermal activation turns out to be more probable. This nontrivial result is due to the subtle balance between the barrier height and width, which determine rates of the thermal activation and quantum tunneling, respectively.

Abstract: We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k . p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extravalley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

Abstract: The adsorption characteristics of alkali, alkaline-earth, and transition metal adatoms on silicene, a graphene-like monolayer structure of silicon are analyzed by means of first-principles calculations. In contrast to graphene, interaction between the metal atoms and the silicene surface is quite strong due to its highly reactive buckled hexagonal structure. In addition to structural properties, we also calculate the electronic band dispersion, net magnetic moment, charge transfer, work function, and dipole moment of the metal adsorbed silicene sheets. Alkali metals, Li, Na, and K, adsorb to hollow sites without any lattice distortion. As a consequence of the significant charge transfer from alkalis to silicene, metalization of silicene takes place. Trends directly related to atomic size, adsorption height, work function, and dipole moment of the silicene/alkali adatom system are also revealed. We found that the adsorption of alkaline-earth metals on silicene is entirely different from their adsorption on graphene. The adsorption of Be, Mg, and Ca turns silicene into a narrow gap semiconductor. Adsorption characteristics of eight transition metals Ti, V, Cr, Mn, Fe, Co, Mo, and W are also investigated. As a result of their partially occupied d orbital, transition metals show diverse structural, electronic, and magnetic properties. Upon the adsorption of transition metals, depending on the adatom type and atomic radius, the system can exhibit metal, half-metal, and semiconducting behavior. For all metal adsorbates, the direction of the charge transfer is from adsorbate to silicene, because of its high surface reactivity. Our results indicate that the reactive crystal structure of silicene provides a rich playground for functionalization at nanoscale. DOI: 10.1103/PhysRevB.87.085423

Abstract: Typical Raman spectra of transition-metal dichalcogenides (TMDs) display two prominent peaks, E-2g and A(1g), that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak in contrast to other TMDs. As the dimensionality is lowered, the observed peak splits in two. In contrast, our ab initio calculations predict that the degeneracy is retained even for WSe2 monolayers. Interestingly, for minuscule biaxial strain, the degeneracy is preserved, but once the crystal symmetry is broken by a small uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a good match to the measured Raman spectrum, which suggests that uniaxial strain exists in WSe2 flakes, possibly induced during the sample preparation and/or as a result of the interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect-to-direct band-gap transition from bulk to monolayers, which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational and electronic properties of WSe2, but also point to effects of the interaction between the monolayer TMDs and the substrate on the vibrational and electronic properties. DOI: 10.1103/PhysRevB.87.165409

Abstract: We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improve if one adds an extra channel. This ensures that there exists a quantum analog of the Braess paradox, originating from quantum scattering and interference.

Abstract: The Green function approach to the Bardeen-Cooper-Schrieffer theory of superconductivity is used to study nanofilms. We go beyond previous models and include effects of confinement on the strength of the electron-phonon coupling as well as on the electronic spectrum and on the phonon modes. Within our approach, we find that in ultrathin films, confinement effects on the electronic screening become very important. Indeed, contrary to what has been advanced in recent years, the sudden increases of the density of states when new bands start to be occupied as the film thickness increases, tend to suppress the critical temperature rather than to enhance it. On the other hand, the increase of the number of phonon modes with increasing number of monolayers in the film leads to an increase in the critical temperature. As a consequence, the superconducting critical parameters in such nanofilms are determined by these two competing effects. Furthermore, in sufficiently thin films, the condensate consists of well-defined subcondensates associated with the occupied bands, each with a distinct coherence length. The subcondensates can interfere constructively or destructively giving rise to an interference pattern in the Cooper pair probability density.

Abstract: We investigate the influence of an external perpendicular magnetic field on the lowest-energy states of an electron bound to a donor which is located near a semiconductor-metal interface. The problem is treated within the effective mass approach and the lowest-energy states are obtained through (1) the “numerically exact” finite element method, and (2) a variational approach using a trial wave function where all image charges that emerge due to the presence of the metallic gate are taken into account. The trial wave functions are constructed such that they reduce to an exponential behavior for sufficiently small magnetic fields and become Gaussian for intermediate and large magnetic fields. The average electron-donor distance can be controlled by the external magnetic field. We find that the size of the 2p(z) state depends strongly on the magnetic field when the donor is close to the interface, showing a nonmonotonic behavior, in contrast with the ground and the other excited states. DOI: 10.1103/PhysRevB.87.075313

Abstract: A parabolic quantum dot (QD) as realized by biasing nanostructured gates on bilayer graphene is investigated in the presence of electron-electron interaction. The energy spectrum and the phase diagram reveal unexpected transitions as a function of a magnetic field. For example, in contrast to semiconductor QDs, we find a valley transition rather than only the usual singlet-triplet transition in the ground state of the interacting system. The origin of these features can be traced to the valley degree of freedom in bilayer graphene. These transitions have important consequences for cyclotron resonance experiments.

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 relation between the energy and momentum of plasmarons in bilayer graphene is investigated within the Overhauser approach, where the electron-plasmon interaction is described as a field theoretical problem. We find that the Dirac-like spectrum is shifted by Delta E(k) similar to 100 divided by 150 meV depending on the electron concentration n(e) and electron momentum. The shift increases with electron concentration as the energy of plasmons becomes larger. The dispersion of plasmarons is more pronounced than in the case of single layer graphene, which is explained by the fact that the energy dispersion of electrons is quadratic and not linear. We expect that these predictions can be verified using angle-resolved photoemission spectroscopy (ARPES).

Abstract: The effect of substitutional doping of fluorographene with boron and nitrogen atoms on its electronic and magnetic properties is investigated using first-principles calculations. It is found that boron dopants can be readily incorporated in the fluorographene crystal where they act as shallow acceptors and cause hole doping, but no changes in the magnetic properties are observed. Nitrogen dopants act as deep donors and give rise to a magnetic moment, but the resulting system becomes chemically unstable. These results are opposite to what was found for substitutional doping of graphane, i.e., hydrogenated graphene, in which case B substituents induce magnetism and N dopants do not.

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.