Abstract: Recently, rhenium disulfide (ReS2) monolayers were experimentally extracted by conventional mechanical exfoliation technique from as-grown ReS2 crystals. Unlike the well-known members of transition metal dichalcogenides (TMDs), ReS2 crystallizes in a stable distorted-1T structure and lacks an indirect to direct gap crossover. Here we present an experimental and theoretical study of the formation, energetics, and stability of the most prominent lattice defects in monolayer ReS2. Experimentally, irradiation with 3-MeV He+2 ions was used to break the strong covalent bonds in ReS2 flakes. Photoluminescence measurements showed that the luminescence from monolayers is mostly unchanged after highly energetic a particle irradiation. In order to understand the energetics of possible vacancies in ReS2 we performed systematic first-principles calculations. Our calculations revealed that the formation of a single sulfur vacancy has the lowest formation energy in both Re and S rich conditions and a random distribution of such defects are energetically more preferable. Sulfur point defects do not result in any spin polarization whereas the creation of Re-containing point defects induce magnetization with a net magnetic moment of 1-3 mu B. Experimentally observed easy formation of sulfur vacancies is in good agreement with first-principles calculations.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: The influence of silicene's strong spin-orbit interaction and of an external electric field E-z on the transport coefficients are investigated in the presence of a perpendicular magnetic field B. For finite E-z the spin and valley degeneracy of the Landau levels is lifted and leads to additional plateaus in the Hall conductivity, at half-integer values of 4e(2)/h, due to spin intra-Landau-level transitions that are absent in graphene. These plateaus are more sensitive to disorder and thermal broadening than the main plateaus, occurring at integral values of 4e(2)/h, when the Fermi level passes through the Landau levels. We also evaluate the Hall and longitudinal resistivities and critically contrast the results with those for graphene on a substrate.

Abstract: We numerically investigate the spin-orbit (SO) coupling effect on the magnetic flux evolution of energy and supercurrent in mesoscopic d-wave superconducting loops by solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. It is found that the energy spectrum splits when the SO interaction is involved and the Majorana zero mode can be realized in the [100] edges of square systems for an appropriate SO coupling strength. Superconducting phase transitions appear when the energy gap closes, accompanied by energy jumps between different energy parabolas in the ground state, which provides a possible mechanism to support fractional flux periodicity of supercurrent. Moreover, in the case of rectangular loops with SO coupling, the jumps of the ground-state energy gradually disappear by increasing the ratio of length to height of the sample, and a paramagnetic response with opposite direction of the screening current around zero flux value can occur in such systems.

Abstract: The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various nonlinear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the fidelity susceptibility of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time scales as long as the ground-state fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating phase, a dissipative regime, or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

Abstract: Using advanced three-dimensional simulations, we show that an Abrikosov vortex, trapped inside a cavity perpendicular to an artificial Josephson junction, can serve as a very efficient source for generation of Josephson vortex-antivortex pairs in the presence of the applied electric current. In such a case, the nucleation rate of the pairs can be tuned in a broad range by an out-of-plane ac magnetic field in a broad range of frequencies. This parametrically amplified vortex-antivortex nucleation can be considered as a macroscopic analog of the dynamic Casimir effect, where fluxon pairs mimic the photons and the ac magnetic field plays the role of the oscillating mirrors. The emerging vortex pairs in our system can be detected by the pronounced features in the measured voltage characteristics, or through the emitted electromagnetic radiation, and exhibit resonant dynamics with respect to the frequency of the applied magnetic field. Reported tunability of the Josephson oscillations can be useful for developing high-frequency emission devices.

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