Abstract: Twisted graphene layers produce a moire pattern (MP) structure with a predetermined wavelength for a given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory, we define the MP structure as the minimum-energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB-and AA-stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy results across multiple substrates is reported as a function of twist angle.

Abstract: In honeycomb Dirac systems with broken inversion symmetry, orbital magnetic moments coupled to the valley degree of freedom arise due to the topology of the band structure, leading to valley-selective optical dichroism. On the other hand, in Dirac systems with prominent spin-orbit coupling, similar orbital magnetic moments emerge as well. These moments are coupled to spin, but otherwise have the same functional form as the moments stemming from spatial inversion breaking. After reviewing the basic properties of these moments, which are relevant for a whole set of newly discovered materials, such as silicene and germanene, we study the particular impact that these moments have on graphene nanoengineered barriers with artificially enhanced spin-orbit coupling. We examine transmission properties of such barriers in the presence of a magnetic field. The orbital moments are found to manifest in transport characteristics through spin-dependent transmission and conductance, making them directly accessible in experiments. Moreover, the Zeeman-type effects appear without explicitly incorporating the Zeeman term in the models, i.e., by using minimal coupling and Peierls substitution in continuum and the tight-binding methods, respectively. We find that a quasiclassical view is able to explain all the observed phenomena.

Abstract: Using the quasiclassical formalism, we determine the low-temperature phase diagram of a quasi-one-dimensional superconductor, taking into account the interchain Josephson coupling and the paramagnetic spin splitting. We show that the anisotropy of the onset of superconductivity changes in the FFLO state as compared with the conventional superconducting phase. It can result in anomalous peaks in the field-direction dependence of the upper critical field when the magnetic field length equals to the FFLO period. This regime is characterized by the lock-in effect of the FFLO modulation wave vector, which is governed by the magnetic length. Furthermore, in the FFLO phase, the anisotropy of the upper critical field is inverted at T-1(**) = 0.5T(c0), where the orbital anisotropy disappears. We suggest that an experimental study of the anisotropy of the upper critical field can provide very reach information about the parameters of the FFLO phase in quasi-1D samples.

Abstract: We present a theoretical study of the plasmon and coupled plasmon-phonon modes induced by intraband electron-electron interaction in graphene in the presence of driving dc electric field. We find that the electric field dependence of these collective excitation modes in graphene differs significantly from that in a conventional two-dimensional electron gas with a parabolic energy spectrum. This is due mainly to the fact that graphene has a linear energy spectrum and the Fermi velocity of electrons in graphene is much larger than the drift velocity of electrons. The obtained results demonstrate that the plasmon and coupled plasmon-phonon modes in graphene can be tuned by applying not only the gate voltage but also the source-to-drain field. The manipulation of plasmon and coupled plasmon-phonon modes by source-to-drain voltage can let graphene be more conveniently applied as an advanced plasmonic material.

Abstract: Quantum thermal transport in armchair and zigzag graphene nanoribbons is investigated in the presence of single atomic vacancies and subject to different boundary conditions. We start with a full comparison of the phonon polarizations and energy dispersions as given by a fifth-nearest-neighbor force-constant model (5NNFCM) and by elasticity theory of continuum membranes (ETCM). For free-edge ribbons, we discuss the behavior of an additional acoustic edge-localized flexural mode, known as fourth acoustic branch (4ZA), which has a small gap when it is obtained by the 5NNFCM. Then, we show that ribbons with supported edges have a sample-size dependent energy gap in the phonon spectrum which is particularly large for in-plane modes. Irrespective to the calculation method and the boundary condition, the dependence of the energy gap for the low-energy optical phonon modes against the ribbon width W is found to be proportional to 1/W for in-plane, and 1/W-2 for out-of-plane phonon modes. Using the 5NNFCM, the ballistic thermal conductance and its contributions from every single phonon mode are then obtained by the nonequilibrium Green's function technique. We found that, while edge and central localized single atomic vacancies do not affect the low-energy transmission function of in-plane phonon modes, they reduce considerably the contributions of the flexural modes. On the other hand, in-plane modes contributions are strongly dependent on the boundary conditions and at low temperatures can be highly reduced in supported-edge samples. These findings could open a route to engineer graphene based devices where it is possible to discriminate the relative contribution of polarized phonons and to tune the thermal transport on the nanoscale.

Abstract: The electronic properties of the two-dimensional material silicene are strongly influenced by the application of a perpendicular electric field E-z and of an exchange field M due to adatoms positioned on the surface or a ferromagnetic substrate. Within the random phase approximation, we investigate how electron-electron interactions are affected by these fields and present analytical and numerical results for the dispersion of plasmons, their lifetime, and their oscillator strength. We find that the combination of the fields E-z and M brings a spin and valley texture to the particle-hole excitation spectrum and allows the formation of spin-and valley-polarized plasmons. When the Fermi level lies in the gap of one spin in one valley, the intraband region of the corresponding spectrum disappears. For zero E-z and finite M the spin symmetry is broken and spin polarization is possible. The lifetime and oscillator strength of the plasmons are shown to depend strongly on the number of spin and valley type electrons that form the electron-hole pairs.

Abstract: We revisit the problem of the dynamic response of a superconducting bridge after abruptly switching on the supercritical current. In contrast to previous theoretical works we take into account spatial gradients and use both the local temperature approach and the kinetic equation for the distribution function of quasiparticles. We find that the temperature dependence of the finite delay time t(d) in the voltage response is model dependent and relatively large t(d) is connected with temporary cooling of quasiparticles during decay of superconducting order parameter vertical bar Delta vertical bar in time. It turns out that the presence of even small inhomogeneities in the bridge or finite length of the homogenous bridge favors a local suppression of vertical bar Delta vertical bar during the dynamic response. It results in a decrease of the delay time, in comparison with the spatially uniform model, due to the diffusion of nonequilibrium quasiparticles from the region with locally suppressed vertical bar Delta vertical bar. In the case when the current density is maximal near the edge of a not very wide bridge the delay time is mainly connected with the time needed for the nucleation (entrance) of the first vortex and t(d) could be tuned by a weak external magnetic field. We also find that a short alternating current pulse (sinusoidlike) with zero time average may result in a nonzero time- averaged voltage response where its sign depends on the phase of the ac current.

Abstract: Using first principles calculations we showed that the electronic and optical properties of single-layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example, while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G(0)W(0) scheme. Applied tensile strain was shown to significantly enhance electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.

Abstract: We calculate the spin-polarized electronic transport through a molecular bilayer spin valve from first principles, and establish the link between the magnetoresistance and the spin-dependent interactions at the metal-molecule interfaces. The magnetoresistance of a Fe vertical bar bilayer-C-70 vertical bar Fe spin valve attains a high value of 70% in the linearresponse regime, but it drops sharply as a function of the applied bias. The current polarization has a value of 80% in linear response and also decreases as a function of bias. Both these trends can be modeled in terms of prominent spin-dependent Fe vertical bar C-70 interface states close to the Fermi level, unfolding the potential of spinterface science to control and optimize spin currents.

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.