Abstract: Starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs), we construct a four-band Hamiltonian describing interlayer excitons consisting of an electron in one TMD layer and a hole in the other TMD layer. An expression for the electron-hole interaction potential is derived, taking into account the effect of the dielectric environment above, below, and between the two TMD layers as well as polarization effects in the transition metal layer and in the chalcogen layers of the TMD layers. We calculate the interlayer exciton binding energy and average in-plane interparticle distance for different TMD heterostructures. The effect of different dielectric environments on the exciton binding energy is investigated and a remarkable dependence on the dielectric constant of the barrier between the two layers is found, resulting from competing effects as a function of the in-plane and out-of-plane dielectric constants of the barrier. The polarization effects in the chalcogen layers, which in general reduce the exciton binding energy, can lead to an increase in binding energy in the presence of strong substrate effects by screening the substrate. The excitonic absorbance spectrum is calculated and we show that the interlayer exciton peak depends linearly on a perpendicular electric field, which agrees with recent experimental results.

Abstract: We investigate the effects related to vortex-core deformations when vortices approach each other. As a result of these vortex-core deformations, the vortex-vortex interaction effectively acquires an attractive component leading to a variety of vortex patterns typical for systems with nonmonotonic repulsive-attractive interaction, such as stripes and labyrinths. The core deformations are anisotropic and can induce frustration in the vortex-vortex interaction. In turn, this frustration has an impact on the resulting vortex patterns, which are analyzed in the presence of additional random pinning, as a function of the pinning strength. This analysis can be applicable to vortices in multiband superconductors or to vortices in Bose-Einstein condensates.

Abstract: Magnetization reversal is a well-studied problem with obvious applicability in computer hard drives. One can accomplish a magnetization reversal in at least one of two ways: application of a magnetic field or through a spin current. The latter is more amenable to a fully quantum-mechanical analysis. We formulate and solve the problem whereby a spin current interacts with a ferromagnetic Heisenberg spin chain, to eventually reverse the magnetization of the chain. Spin flips are accomplished through both elastic and inelastic scattering. A consequence of the inelastic-scattering channel, when it is no longer energetically possible, is the occurrence of a nonequilibrium bound state, which is an emergent property of the coupled local plus itinerant spin system. For certain definite parameter values the itinerant spin lingers near the local spins for some time, before eventually leaking out as an outwardly diffusing state. This phenomenon results in spin-flip dynamics and filtering properties for this type of system.

Abstract: We report the direct measurement of the topological skyrmion energy barrier through a hysteresis of the skyrmion lattice in the chiral magnet MnSi. Measurements were made using small-angle neutron scattering with a custom-built resistive coil to allow for high-precision minor hysteresis loops. The experimental data were analyzed using an adapted Preisach model to quantify the energy barrier for skyrmion formation and corroborated by the minimum-energy path analysis based on atomistic spin simulations. We reveal that the skyrmion lattice in MnSi forms from the conical phase progressively in small domains, each of which consisting of hundreds of skyrmions, and with an activation barrier of several eV.

Abstract: The interplay between a magnetic field and the Coulomb potential from a charged vacancy on the electron states in graphene is investigated within the tight-binding model. The Coulomb potential removes locally Landau level degeneracy, while the vacancy introduces a satellite level next to the normal Landau level. These satellite levels are found throughout the positive-energy region, but in the negative-energy region, they turn into atomic collapse resonances. Crossings between Landau levels with different angular quantum number m are found. Unlike the point impurity system in which an anticrossing occurs between Landau levels of the same m, in this work anticrossing is found between the normal Landau level and the vacancy-induced level. The atomic collapse resonance hybridizes with the Landau levels. The charge at which the lowest Landau level m = -1, N = 1 crosses E = 0 increases with enhancing magnetic field. A Landau level scaling anomaly occurs when the charge is larger than the critical charge beta 0.6 and this critical charge is independent of the magnetic field.

Abstract: By using the full 3D generalized time-dependent Ginzbug-Landau equation, we study a long superconducting film of finite width and thickness under an applied transport current. We show that, for sufficiently large thickness, the vortices and the antivortices become curved before they annihilate each other. As they approach the center of the sample, their ends combine, producing a single closed vortex. We also determine the critical values of the thickness for which the closed vortex sets in for different values of the Ginzburg-Ladau parameter. Finally, we propose a model of how to detect a closed vortex experimentally.

Abstract: Starting from first principles, we show the formation and evolution of superconducting gaps in MgB2 at its ultrathin limit. Atomically thin MgB2 is distinctly different from bulk MgB2 in that surface states become comparable in electronic density to the bulklike sigma and pi bands. Combining the ab initio electron-phonon coupling with the anisotropic Eliashberg equations, we showthat monolayer MgB2 develops three distinct superconducting gaps, on completely separate parts of the Fermi surface due to the emergent surface contribution. These gaps hybridize nontrivially with every extra monolayer added to the film owing to the opening of additional coupling channels. Furthermore, we reveal that the three-gap superconductivity in monolayer MgB2 is robust over the entire temperature range that stretches up to a considerably high critical temperature of 20 K. The latter can be boosted to >50K under biaxial tensile strain of similar to 4%, which is an enhancement that is stronger than in any other graphene-related superconductor known to date.

Abstract: The McMillan-Ginzburg-Landau (MGL) model for charge density waves (CDW) is employed in a systematic phenomenological study of the different phases that have been probed in recent experiments involving transition metal dichalcogenides. We implemented an efficient imaginary time evolution method to solve the MGL equations, which enabled us to investigate the role of different coupling parameters on the CDW patterns and to perform calculations with different energy functionals that lead to several experimentally observed singularities in the CDW phase profiles. In particular, by choosing the appropriate energy functionals, we were able to obtain phases that go beyond the well-known periodic phase slips (discommensurations), exhibiting also topological defects (i.e., vortex-antivortex pairs), domain walls where the CDW order parameter is suppressed, and even CDW with broken rotational symmetry. Finally, we briefly discuss the effect of these different CDW phases on the profile and critical temperature of the competing superconducting state.

Abstract: To investigate the impact of both the CoFe ferromagnetic layer thickness and the capping paramagnetic layer on the tunnel magnetoresistance (TMR), we performed first-principles simulations on epitaxial magnetic tunnel junctions contacted with either CoFe or Ta paramagnetic capping layers. We observed a strong oscillation of the TMR amplitude with respect to the thickness of the ferromagnetic layer. The TMR is found to be amplified whenever the MgO spin tunnel barrier is thickened. Quantization of the electronic structure of the ferromagnetic layers is found to be at the origin of this oscillatory behavior. Metals such as Ta contacting the magnetic layer are found to enhance the amplitude of the oscillations due to the occurrence of an interface dipole. The latter drives the band alignment and tunes the nature of the spin channels that are active during the tunneling process. Subsequently, the regular transmission spin channels are modulated in the magnetic tunnel junction stack and other complex ones are being activated.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: We propose a tight-binding (TB) model, that includes spin-orbit coupling (SOC), to describe the electronic properties of methyl-substituted germanane (GeCH3). This model gives an electronic spectrum in agreement with first principle results close to the Fermi level. Using the Z(2) formalism, we show that a topological phase transition from a normal insulator (NI) to a quantum spin Hall (QSH) phase occurs at 11.6% biaxial tensile strain. The sensitivity of the electronic properties of this system on strain, in particular its transition to the topological insulating phase, makes it very attractive for applications in strain sensors and other microelectronic applications.

Abstract: Using first-principles calculations, we study the dependence of the electronic and vibrational properties of multilayered PbI2 crystals on the number of layers and focus on the electronic-band structure and the Raman spectrum. Electronic-band structure calculations reveal that the direct or indirect semiconducting behavior of PbI2 is strongly influenced by the number of layers. We find that at 3L thickness there is a direct-to-indirect band gap transition (from bulk-to-monolayer). It is shown that in the Raman spectrum two prominent peaks, A(1g) and E-g, exhibit phonon hardening with an increasing number of layers due to the interlayer van der Waals interaction. Moreover, the Raman activity of the A(1g) mode significantly increases with an increasing number of layers due to the enhanced out-of-plane dielectric constant in the few-layer case. We further characterize rigid-layer vibrations of low-frequency interlayer shear (C) and breathing (LB) modes in few-layer PbI2. A reduced monatomic (linear) chain model (LCM) provides a fairly accurate picture of the number of layers dependence of the low-frequency modes and it is shown also to be a powerful tool to study the interlayer coupling strength in layered PbI2.

Abstract: Using the tight-binding (TB) approximation with inclusion of the spin-orbit interaction, we predict a topological phase transition in the electronic band structure of phosphorene in the presence of axial strains. We derive a low-energy TB Hamiltonian that includes the spin-orbit interaction for bulk phosphorene. Applying a compressive biaxial in-plane strain and perpendicular tensile strain in ranges where the structure is still stable leads to a topological phase transition. We also examine the influence of strain on zigzag phosphorene nanoribbons (zPNRs) and the formation of the corresponding protected edge states when the system is in the topological phase. For zPNRs up to a width of 100 nm the energy gap is at least three orders of magnitude larger than the thermal energy at room temperature.

Abstract: Combining first-principles calculations with Landauer-Mittiker formalism, ballistic thermoelectric transport properties of semiconducting two-dimensional transition metal dichalcogenides (TMDs) and oxides (TMOs) (namely MX2 with M = Cr, Mo, W, Ti, Zr, Hf; X = O, S, Se, Te) are investigated in their 2H and 1T phases. Having computed structural, as well as ballistic electronic and phononic transport properties for all structures, we report the thermoelectric properties of the semiconducting ones. We find that 2H phases of four of the studied structures have very promising thermoelectric properties, unlike their 1T phases. The maximum room temperature p-type thermoelectric figure of merit (ZT) of 1.57 is obtained for 2H-HfSe2, which can be as high as 3.30 at T = 800 K. Additionally, 2H-ZrSe2, 2H-ZrTe2, and 2H-HfS2 have considerable ZT values (both nand p-type), that are above 1 at room temperature. The 1T phases of Zr and Hf-based oxides possess relatively high power factors, however their high lattice thermal conductance values limit their ZT values to below 1 at room temperature.

Abstract: By making use of circularly polarized light and electrostatic gating, monolayer molybdenum disulfide (ML – MoS2) can form a platform supporting multiple types of charge carriers. They can be discriminated by their spin, valley index, or whether they are electrons or holes. We investigate the collective properties of those charge carriers and are able to identify distinct plasmon modes. We analyze the corresponding dispersion relation, lifetime, and oscillator strength, and calculate the phase relation between the oscillations in the different components of the plasmon modes. All platforms in ML-MoS2 support a long-wavelength root q plasmon branch at zero kelvins. In addition to this, for an n-component system, n-1 distinct plasmon modes appear as acoustic modes with linear dispersion in the long-wavelength limit. These modes correspond to out-of-phase oscillations in the different fermion liquids and have, although being damped, a relatively long lifetime. Additionally, we also find distinct modes at large wave vectors that are more strongly damped by intraband processes.