Abstract: Dynamics of magnetic skyrmions in hybrid ferromagnetic films harbors interesting physical phenomena and holds promise for technological applications. In this work, we discuss the behavior of magnetic skyrmions when coupled to superconducting vortices in a ferromagnet-superconductor heterostructure. We use numerical simulations and analytic arguments within London and Thiele formalisms to reveal broader possibilities for manipulating the skyrmion-vortex dynamic correlations in the hybrid system, that are not possible in its separated constituents. We explore the thresholds of particular dynamic phases, and quantify the phase diagram as a function of the relevant material parameters, applied current, and induced magnetic torques. Finally, we demonstrate the broad and precise tunability of the skyrmion Hall angle in the presence of vortices, with respect to currents applied to either or both the superconductor and the ferromagnet within the heterostructure.

Abstract: The electronic and transport properties of a three-dimensional (3D) Dirac system are investigated theoretically, which is motivated by recent experimental measurements on quantum and transport mobilities in the 3D Dirac semimetal Na3Bi by J. Xiong et al. [Science 350, 413 (2015); Europhys. Lett. 114, 27002 (2016)]. The electron Hamiltonian is taken from a simplified k center dot p approach. From the obtained electronic band structure and the Fermi energy, we explain why the anomalous effect induced by the chiral anomaly and the Berry curvature in the energy band can be observed experimentally in magnetotransport coefficients in both low-and high-density samples. Moreover, the quantum and transport mobilities are calculated on the basis of the momentum-balance equation derived from a semiclassical Boltzmann equation with the electron-impurity interaction. The quantum and transport mobilities obtained from this study agree both qualitatively and quantitatively with those measured experimentally. We also examine the electron mobilities along different crystal directions in Na3Bi and find them largely anisotropic. The theoretical findings from this work can be helpful in gaining an in-depth understanding of the experimental results and of the basic electronic and transport properties of newly developed 3D Dirac systems.

Abstract: Within a tight-binding approach, we investigate the effect of twisting angle on the energy levels of circular bilayer graphene (BLG) quantum dots (QDs) in both the absence and presence of a perpendicular magnetic field. The QDs are defined by an infinite-mass potential, so that the specific edge effects are not present. In the absence of magnetic field (or when the magnetic length is larger than the moire length), we show that the low-energy states in twisted BLG QDs are completely affected by the formation of moire patterns, with a strong localization at AA-stacked regions. When magnetic field increases, the energy gap of an untwisted BLG QD closes with the edge states, localized at the boundaries between the AA- and AB-stacked spots in a twisted BLG QD. Our observation of the spatial localization of the electrons in twisted BLG QDs can be experimentally probed by low-bias scanning tunneling microscopy measurements.

Abstract: The hunt for high-temperature superfluidity has received new impetus from the discovery of atomically thin stable materials. Electron-hole superfluidity in coupled MoSe2-WSe2 monolayers is investigated using a mean-field multiband model that includes band splitting caused by strong spin-orbit coupling. This splitting leads to a large energy misalignment of the electron and hole bands which is strongly modified by interchanging the doping of the monolayers. The choice of doping determines if the superfluidity is tunable from one to two components.

Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: Atomic scale defects in semiconductors enable their technological applications and realization of different quantum states. Using scanning tunneling microscopy and spectroscopy complemented by ab initio calculations we determine the nature of defects in the anisotropic van der Waals layered semiconductor ReS2. We demonstrate the in-plane anisotropy of the lattice by directly visualizing chains of rhenium atoms forming diamond-shaped clusters. Using scanning tunneling spectroscopy we measure the semiconducting gap in the density of states. We reveal the presence of lattice defects and by comparison of their topographic and spectroscopic signatures with ab initio calculations we determine their origin as oxygen atoms absorbed at lattice point defect sites. These results provide an atomic-scale view into the semiconducting transition metal dichalcogenides, paving the way toward understanding and engineering their properties.

Abstract: By using density functional theory (DFT) calculations, we predict a puckered, dynamically stable Janus single-layer black arsenic-phosphorus (b-AsP), which is composed of two different atomic sublayers, arsenic and phosphorus atoms. The calculated phonon spectrum reveals that Janus single-layer b-AsP is dynamically stable with either pure or coupled optical phonon branches arising from As and P atoms. The calculated Raman spectrum indicates that due to the relatively strong P-P bonds, As atoms have no contribution to the highfrequency optical vibrations. In addition, the orientation-dependent isovolume heat capacity reveals anisotropic contributions of LA and TA phonon branches to the low-temperature thermal properties. Unlike pristine single layers of b-As and b-P, Janus single-layer b-AsP exhibits additional out-of-plane asymmetry which leads to important consequences for its electronic, optical, and elastic properties. In contrast to single-layer b-As, Janus single-layer b-AsP is found to possess a direct band gap dominated by the P atoms. Moreover, real and imaginary parts of the dynamical dielectric function, including excitonic effects, reveal the highly anisotropic optical feature of the Janus single-layer. A tight-binding (TB) model is also presented for Janus single-layer b-AsP, and it is shown that, with up to seven nearest hoppings, the TB model reproduces well the DFT band structure in the low-energy region around the band gap. This TB model can be used in combination with the Green's function approach to study, e.g., quantum transport in finite systems based on Janus single-layer b-AsP. Furthermore, the linear-elastic properties of Janus single-layer b-AsP are investigated, and the orientation-dependent in-plane stiffness and Poisson ratio are calculated. It is found that the Janus single layer exhibits strong in-plane anisotropy in its Poisson ratio much larger than that of single-layer b-P. This Janus single layer is relevant for promising applications in optical dichroism and anisotropic nanoelasticity.

Abstract: High tunability of two-dimensional magnetic materials (by strain, gating, heterostructuring, or otherwise) provides unique conditions for studying versatile magnetic properties and controlling emergent magnetic phases. Expanding the scope of achievable magnetic phenomena in such materials is important for both fundamental and technological advances. Here we perform atomistic spin-dynamics simulations to explore the (chiral) magnetic phases of atomic monolayers in the limit of suppressed first-neighbors exchange interaction. We report the rich phase diagram of exotic magnetic configurations, obtained for both square and honeycomb lattice symmetries, comprising coexistence of ferromagnetic and antiferromagnetic spin cycloids, as well as multiple types of magnetic skyrmions. We perform a minimum-energy path analysis for the skyrmion collapse to evaluate the stability of such topological objects and reveal that magnetic monolayers could be good candidates to host the antiferromagnetic skyrmions that are experimentally evasive to date.

Abstract: The effect of proximity-induced interactions such as Rashba spin-orbit coupling (SOC) and exchange interaction on the electronic and optical properties of n-type monolayer (ML) MoS2 is investigated. We predict and demonstrate that the Rashba SOC can induce an in-plane spin splitting with terahertz (THz) energy, while the exchange interaction lifts the energy degeneracy in different valleys. Thus, spin polarization can be achieved in an n-type ML MoS2 and valley Hall or optical Hall effect can be observed using linearly polarized THz radiation. In such a case, the transverse optical conductivity sigma(xy) (omega) results from spin-flip transition within spin-split conduction bands and from the fact that contributions from electrons with different spin orientations in different valleys can no longer be canceled out. Interestingly, we find that for fixed effective Zeeman field (or exchange interaction) the lowest spin-split conduction band in ML MoS2 can be tuned from one in the K valley to another one in the K' valley by varying the Rashba parameter lambda(R). Therefore, by changing lambda(R) we can turn the sign of the spin polarization and Im sigma(xy) (omega) from positive to negative. Moreover, we find that the dominant contribution of the selection rules to sigma(xx)(omega) is from electrons in the K valley and to sigma(xy) (omega) is from electrons in the K' valley. These important and interesting theoretical findings can be helpful to experimental observation of the optical Hall effect in valleytronic systems using linearly polarized THz radiation fields.

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: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: We calculate the Curie temperature of layered ferromagnets, chromium tri-iodide (CrI3), chromium tri-bromide (CrBr3), chromium germanium tri-telluride (CrGeTe3), and the Ned temperature of a layered antiferromagnet iron di-chloride (FeCl2), using first-principles density functional theory calculations and Monte Carlo simulations. We develop a computational method to model the magnetic interactions in layered magnetic materials and calculate their critical temperature. We provide a unified method to obtain the magnetic exchange parameters (J) for an effective Heisenberg Hamiltonian from first principles, taking into account both the magnetic ansiotropy as well as the out-of-plane interactions. We obtain the magnetic phase change behavior, in particular the critical temperature, from the susceptibility and the specific-heat, calculated using the three-dimensional Monte Carlo (METROPOLIS) algorithm. The calculated Curie temperatures for ferromagnetic materials (CrI3, CrBr3, and CrGeTe3), match well with experimental values. We show that the interlayer interaction in bulk CrI3 with R (3) over bar stacking is significantly stronger than the C2/m stacking, in line with experimental observations. We show that the strong interlayer interaction in R (3) over bar CrI3 results in a competition between the in-plane and the out-of-plane magnetic easy axes. Finally, we calculate the Ned temperature of FeCl2 to be 47 +/- 8 K and show that the magnetic phase transition in FeCl2 occurs in two steps with a high-temperature intralayer ferromagnetic phase transition and a low-temperature interlayer antiferromagnetic phase transition.

Abstract: The optoelectronic properties of confined water form one of the most active research areas in the past few years. Here we present the multiscale methodology to discern the out-of-plane electronic and dipolar dielectric constants (epsilon(el)(perpendicular to) and epsilon(diP)(perpendicular to)) of strongly confined water. We reveal that epsilon(perpendicular to el) and epsilon(diP)(perpendicular to) become comparable for water confined in angstrom-scale channels (with a height of less than 15 angstrom) within graphene (GE) and hexagonal boron nitride (hBN) bilayers. Channel height (h) associated with a minimum in both epsilon(e)(l)(perpendicular to) and epsilon(dip)(perpendicular to) is linked to the formation of the ordered structure of ice for h approximate to (7 -7.5) angstrom. The recently measured total dielectric constant epsilon(T)(perpendicular to) of nanoconfined water [L. Fumagalli et al., Science 360, 1339 (2018)] is corroborated by our results. Furthermore, we evaluate the contribution from the encapsulating membranes to the dielectric properties, as a function of the interlayer spacing, i.e., the height of the confining channel for water. Finally, we conduct analysis of the optical properties of both confined water and GE membranes, and show that the electron energy loss function of confined water strongly differs from that of bulk water.

Abstract: Atomic collapse in graphene nanoribbons behaves in a fundamentally different way as compared to monolayer graphene due to the presence of multiple energy bands and the effect of edges. For armchair nanoribbons we find that bound states gradually transform into atomic collapse states with increasing impurity charge. This is very different in zigzag nanoribbons where multiple quasi-one-dimensional bound states are found that originates from the zero-energy zigzag edge states. They are a consequence of the flat band and the electron distribution of these bound states exhibits two peaks. The lowest-energy edge state transforms from a bound state into an atomic collapse resonance and shows a distinct relocalization from the edge to the impurity position with increasing impurity charge.

Abstract: The modification of the quantum states in a Dice lattice due to a Coulomb impurity are investigated. The energy-band structure of a pristine Dice lattice consists of a Dirac cone and a flat band at the Dirac point. We use the tight-binding formalism and find that the flat band states transform into a set of discrete bound states whose electron density is localized on a ring around the impurity mainly on two of the three sublattices. Its energy is proportional to the strength of the Coulomb impurity. Beyond a critical strength of the Coulomb potential atomic collapse states appear that have some similarity with those found in graphene with the difference that the flat band states contribute with an additional ringlike electron density that is spatially decoupled from the atomic collapse part. At large value of the strength of the Coulomb impurity the flat band bound states anticross with the atomic collapse states.