Abstract: Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information similarly to the way spintronics uses electron spin. We show how a double-barrier structure consisting of electric and vector potentials can be used to filter massive Dirac electrons based on their valley index. We study the resonant transmission through a finite number of barriers and we obtain the energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is included, the energy bands and energy gaps at the K and K′ points are different and they can be tuned by changing the potential.

Abstract: Iron carbides play a crucial role in steel manufacturing and processing and to a large extent determine the physical properties of steel products. The modified embedded atom method (MEAM) in combination with Lee's Fe-C potential is a good candidate for molecular dynamics simulations on larger Fe-C systems. Here, we investigate the stability and crystal structures of pure iron and binary iron carbides using MEAM and compare them with the experimental data and quantum-mechanical density functional theory calculations. The analysis shows that the Fe-C potential gives reasonable results for the relative stability of iron and iron carbides. The performance of MEAM for the prediction of the potential energy and the calculated lattice parameters at elevated temperature for pure iron phases and cementite are investigated as well. The conclusion is that Lee's MEAM Fe-C potential provides a promising basis for further molecular dynamics simulations of Fe-C alloys and steels at lower temperatures (up to 800 K).

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

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