Abstract: Structural phase transitions in two vertically or horizontally coupled channels of strongly interacting particles are investigated. The particles are free to move in the x direction but are confined by a parabolic potential in the y direction. They interact with each other through a screened power-law potential (r(-n)e(-r/lambda)). In vertically coupled systems, the channels are stacked above each other in the direction perpendicular to the (x, y) plane, while in horizontally coupled systems both channels are aligned in the confinement direction. Using Monte Carlo (MC) simulations we obtain the ground-state configurations and the structural transitions as a function of the linear particle density and the separation between the channels. At zero temperature, the vertically coupled system exhibits a rich phase diagram with continuous and discontinuous transitions. On the other hand, the horizontally coupled system exhibits only a very limited number of phase transitions due to its symmetry. Further, we calculated the normal modes for the Wigner crystals in both cases. From MC simulations, we found that in the case of vertically coupled systems, the zigzag transition is only possible for low densities. A Ginzburg-Landau theory for the zigzag transition is presented, which predicts correctly the behavior of this transition from which we interpret the structural phase transition of the Wigner crystal through the reduction of the Brillouin zone.

Abstract: To date, several experimental groups reported measurements of the thickness dependence of T-c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T-c oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon-mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover, and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T-c oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

Abstract: The lattice dynamics of low-frequency rigid-plane modes in metallic (graphene multilayers, GML) and in insulating (hexagonal boron-nitride multilayers, BNML) layered crystals is investigated. The frequencies of shearing and compression (stretching) modes depend on the layer number N and are presented in the form of fan diagrams. The results for GML and BNML are very similar. In both cases, only the interactions (van der Waals and Coulomb) between nearest-neighbor planes are effective, while the interactions between more distant planes are screened. A comparison with recent Raman scattering results on low-frequency shear modes in GML [Tan et al., Nat. Mater., in press, doi: 10.1038/nmat3245, (2012)] is made. Relations with the low-lying rigid-plane phonon dispersions in the bulk materials are established. Master curves, which connect the fan diagram frequencies for any given N, are derived. Static and dynamic thermal correlation functions for rigid-layer shear and compression modes are calculated. The results might be of use for the interpretation of friction force experiments on multilayer crystals.

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: 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: Vortex loops are generated by the inhomogeneous stray field of a magnetic dipole on top of a current-carrying mesoscopic superconductor. Cutting and recombination processes unfold under the applied drive, resulting in periodic voltage oscillations across the sample. We show that a direct and detectable consequence of the cutting and recombination of these vortex loops in the present setup is the onset of vortices at surfaces where they were absent prior to the application of the external current. The nonlinear dynamics of vortex loops is studied within the time-dependent Ginzburg-Landau theory to describe the profound three-dimensional features of their time evolution.

Abstract: We investigate the effect of Rashba and intrinsic spin-orbit couplings on the electronic properties and spin configurations of Dirac fermions confined in: (i) a flat graphene sheet, (ii) a graphene wire with p-n-p structure, and (iii) a superlattice of graphene wires. The interplay between the spin-orbit interaction mechanisms breaks the electron-hole symmetry and the spin configuration induced by Rashba spin-orbit coupling lacks inversion symmetry in k space. We show that the Rashba spin-orbit interaction doubles the Fabry-Perot resonant modes in the transmission spectrum of a graphene wire and opens new channels for the electron transmission. Moreover, it leads to the appearance of spin split extra Dirac cones in the energy spectrum of a graphene superlattice. It is shown that the spin of the electrons and holes confined in a flat graphene sheet is always perpendicular to their motion while this is not the case for the other nanostructures.

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