Abstract: In the Comment, the authors state that the separation of the relative and center of mass variables in our work is not correct. Here we point out that there is a typographical error, i.e., qi instead of -e, in two of our equations which, when corrected, makes the Comment redundant. Within the ansatzes mentioned in our paper all our results are correct, in contrast to the claims of the Comment.

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: We report characteristic vortex configurations in s + id superconductors with time-reversal symmetry breaking, exposed to magnetic field. A vortex in the s + id state tends to have an opposite phase winding between s- and d-wave condensates. We find that this peculiar feature together with the competition between s- and d-wave symmetry results in three distinct classes of vortical configurations. When either s or d condensate absolutely dominates, vortices form a conventional lattice. However, when one condensate is relatively dominant, vortices organize in chains that exhibit skyrmionic character, separating the chiral components of the s +/- id order parameter into domains within and outside the chain. Such skyrmionic chains are found stable even at high magnetic field. When s and d condensates have comparable strength, vortices split cores in two chiral components to form full-fledged skyrmions, i.e., coreless topological structures with an integer topological charge, organized in a lattice. We provide characteristic magnetic field distributions of all states, enabling their identification in, e.g., scanning Hall probe and scanning SQUID experiments. These unique vortex states are relevant for high-T-c cuprate and iron-based superconductors, where the relative strength of competing pairing symmetries is expected to be tuned by temperature and/or doping level, and can help distinguish s + is and s + id superconducting phases.

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 chiral p-wave order parameter in Sr2RuO4 would make it a special case amongst the unconventional superconductors. A consequence of this symmetry is the possible existence of superconducting domains of opposite chirality. At the boundary of such domains, the locally suppressed condensate can produce an intrinsic Josephson junction. Here, we provide evidence of such junctions using mesoscopic rings, structured from Sr2RuO4 single crystals. Our order parameter simulations predict such rings to host stable domain walls across their arms. This is verified with transport experiments on loops, with a sharp transition at 1.5 K, which show distinct critical current oscillations with periodicity corresponding to the flux quantum. In contrast, loops with broadened transitions at around 3 K are void of such junctions and show standard Little-Parks oscillations. Our analysis demonstrates the junctions are of intrinsic origin and makes a compelling case for the existence of superconducting domains.

Abstract: By means of density functional theory based first-principles calculations, the structural, vibrational, and electronic properties of 1H- and 1T-phases of single-layer CaX2 (X = F, Cl, Br, or I) structures are investigated. Our results reveal that both the 1H- and 1T-phases are dynamically stable in terms of their phonon band dispersions with the latter being the energetically favorable phase for all single-layers. In both phases of single-layer CaX2 structures, significant phonon softening occurs as the atomic radius increases. In addition, each structural phase exhibits distinctive Raman active modes that enable one to characterize either the phase or the structure via Raman spectroscopy. The electronic band dispersions of single-layer CaX2 structures reveal that all structures are indirect bandgap insulators with a decrease in bandgaps from fluorite to iodide crystals. Furthermore, the calculated linear elastic constants, in-plane stiffness, and Poisson ratio indicate the ultra-soft nature of CaX2 single-layers, which is quite important for their nanoelastic applications. Overall, our study reveals that with their dynamically stable 1T- and 1H-phases, single-layers of CaX2 crystals can be alternative ultra-thin insulators.

Abstract: We investigate experimentally the terahertz (THz) optoelectronic properties of monolayer (ML) tungsten disulfide (WS2) placed on different substrates using THz time-domain spectroscopy (TDS). We find that the THz optical response of n-type ML WS2 depends sensitively on the choice of the substrate. This dependence is found to be a consequence of substrate induced charge transfer, extra scattering centers, and electronic localization. Through fitting the experimental results with the Drude-Smith formula, we can determine the key sample parameters (e.g., the electronic relaxation time, electron density, and electronic localization factor) of ML WS2 on different substrates. The temperature dependence of these parameters is examined. Our results show that the THz TDS technique is an efficient non-contact method that can be utilized to characterize and investigate the optoelectronic properties of nano-devices based on ML WS2.

Abstract: Using Floquet-Bloch theory, we propose to realize chiral topological phases in two-dimensional (2D) hexagonal FeX2 (X=Cl, Br, I) monolayers under irradiation of circularly polarized light. Such 2D FeX2 monolayers are predicted to be dynamically stable and exhibit both ferromagnetic and semiconducting properties. To capture the full topological physics of the magnetic semiconductor under periodic driving, we adopt ab initio Wannier-based tight-binding methods for the Floquet-Bloch bands, with the light-induced bandgap closings and openings being obtained as the light field strength increases. The calculations of slabs with open boundaries show the existence of chiral edge states. Interestingly, the topological transitions with branches of chiral edge states changing from zero to one and from one to two by tuning the light amplitude are obtained, showing that the topological Floquet phase of high Chern number can be induced in the present Floquet-Bloch systems. Published under license by AIP Publishing.

Abstract: Magnetic skyrmions are topologically stabilized nanoscale spin structures that could be of use in the development of future spintronic devices. When a skyrmion is driven by an electric current it propagates at an angle relative to the flow of current-known as the skyrmion Hall angle (SkHA)-that is a function of the drive current. This drive dependence, as well as thermal effects due to Joule heating, could be used to tailor skyrmion trajectories, but are not well understood. Here we report a study of skyrmion dynamics as a function of temperature and drive amplitude. We find that the skyrmion velocity depends strongly on temperature, while the SkHA does not and instead evolves differently in the low- and high-drive regimes. In particular, the maximum skyrmion velocity in ferromagnetic devices is limited by a mechanism based on skyrmion surface tension and deformation (where the skyrmion transitions into a stripe). Our mechanism provides a complete description of the SkHA in ferromagnetic multilayers across the full range of drive strengths, illustrating that skyrmion trajectories can be engineered for device applications. An analysis of skyrmion dynamics at different temperatures and electric drive currents is used to develop a complete description of the skyrmion Hall angle in ferromagnetic multilayers from the creep to the flow regime and illustrates that skyrmion trajectories can be engineered for device applications.

Abstract: We present the Tight-Binding Studio (TB Studio) software package that calculates the different parameters of a tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or semi-empirical band-structure theory. This will be helpful for scientists who are interested in studying electronic and optical properties of structures using Green's function theory within the tight-binding approximation. TB Studio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high performance OpenGL graphic library for real time plotting. TB Studio and its examples together with the tutorials are available for download from tight-binding.com. Program summary Program Title: Tight-Binding Studio Program Files doi:http://dx.doi.org/10.17632/j6x5mwzm2d.1 Licensing provisions: LGPL Programming language: C++ External routines: BLAS, LAPACK, LAPACKE, wxWidgets, OpenGL, MathGL Nature of problem: Obtaining Tight-Binding Hamiltonian from a set of Bloch energy bands obtained from first-principles calculations. Solution method: Starting from the simplified LCAO method, a tight-binding model in the two-center approximation is constructed. The Slater and Koster (SK) approach is used to calculate the parameters of the TB Hamiltonian. By using non-linear fitting approaches the optimal values of the SK parameters are obtained such that the TB energy eigenvalues are as close as possible to those from first-principles calculations. We obtain the expression for the Hamiltonian and the overlap matrix elements between the different orbitals of the different atoms in an orthogonal or non-orthogonal basis set. (C) 2020 Elsevier B.V. All rights reserved.

Abstract: Condensation of spatially indirect excitons, with the electrons and holes confined in two separate layers, has recently been observed in two different double layer heterostructures. High transition temperatures were reported in a double Transition Metal Dichalcogenide (TMD) monolayer system. We briefly review electron-hole double layer systems that have been proposed as candidates for this interesting phenomenon. We investigate the double TMD system WSe2/hBN/MoSe2, using a mean-field approach that includes multiband effects due to the spin-orbit coupling and self-consistent screening of the electron-hole Coulomb interaction. We demonstrate that the transition temperature observed in the double TMD monolayers, which is remarkably high relative to the other systems, is the result of (i) the large electron and hole effective masses in TMDs, (ii) the large TMD band gaps, and (iii) the presence of multiple superfluid condensates in the TMD system. The net effect is that the superfluidity is strong across a wide range of densities, which leads to high transition temperatures that extend as high as TBKT=150 K.

Abstract: We perform the wavepacket dynamics simulation of a graphene-based device where propagating electron trajectories are tamed by an applied magnetic field toward a normal/superconductor interface. The magnetic field controls the incidence angle of the incoming electronic wavepacket at the interface, which results in the tunable electron-hole ratio in the reflected wave function due to the angular dependence of the Andreev reflection. Here, mapped control of the quasiparticle trajectories by the external magnetic field not only defines an experimental probe for fundamental studies of the Andreev reflection in graphene but also lays the foundation for further development of magnetic focusing devices based on nanoengineered superconducting two-dimensional materials.

Abstract: We report on the anisotropy of the London penetration depth of CaKFe4As4, discussing how it relates to its electronic structure and how it modifies under introduction of disorder, both chemically induced (by Ni substitution) and irradiation induced (by 3.5-MeV protons). Indeed, CaKFe4As4 is particularly suitable for the study of fundamental superconducting properties due to its stoichiometric composition, exhibiting clean-limit behavior in the pristine samples and having a fairly high critical temperature, T-c approximate to 35 K. The London penetration depth lambda(L) is measured with a microwave-coplanar-resonator technique that allows us to deconvolve the anisotropic contributions lambda(L,ab) and lambda(L,c) and obtain the anisotropy parameter gamma(lambda) = lambda(L,c)/lambda(L,ab). The gamma(lambda) (T) found for the undoped pristine sample is in good agreement with previous literature and is here compared to ab initio density-functional-theory and Eliashberg calculations. The dependence of gamma(lambda) (T) on both chemical and irradiation-induced disorder is discussed to highlight which method is more suitable to decrease the direction dependence of the electromagnetic properties while maintaining a high critical temperature. Lastly, the relevance of an intrinsic anisotropy such as gamma(lambda) on application-related anisotropic parameters (critical current, pinning) is discussed in light of the recent employment of CaKFe4As4 in the production of wires.

Abstract: Using density functional theory, we investigate a novel two-dimensional silicon bismotide (SiBi) that has a layered GaSe-like crystal structure.Ab initiomolecular dynamic simulations and phonon dispersion calculations suggest its good thermal and dynamical stability. The SiBi monolayer is a semiconductor with a narrow indirect bandgap of 0.4 eV. Our results show that the indirect bandgap decreases as the number of layers increases, and when the number of layers is more than six layers, direct-to-indirect bandgap switching occurs. The SiBi bilayer is found to be very sensitive to an E-field. The bandgap monotonically decreases in response to uniaxial and biaxial compressive strain, and reaches 0.2 eV at 5%, while at 6%, the semiconductor becomes a metal. For both uniaxial and biaxial tensile strains, the material remains a semiconductor and indirect-to-direct bandgap transition occurs at a strain of 3%. Compared to a SiBi monolayer with a layer thickness of 4.89 angstrom, the bandgap decreases with either increasing or decreasing layer thickness, and at a thicknesses of 4.59 to 5.01 angstrom, the semiconductor-to-metal transition happens. In addition, under pressure, the semiconducting character of the SiBi bilayer with a 0.25 eV direct bandgap is preserved. Our results demonstrate that the SiBi nanosheet is a promising candidate for designing high-speed low-dissipation devices.

Abstract: Motivated by the recent synthesis of two-dimensional VI3 [Kong et al. Adv. Mater. 31, 1808074 (2019)], we investigate the effect of V doping on the magnetic and electronic properties of monolayer VI3 by means of first-principles calculations. The dynamically stable semiconducting ferromagnetic (FM) and antiferromagnetic (AFM) phases of monolayer VI3 are found to display distinctive vibrational features that the magnetic state can be distinguished by Raman spectroscopy. In order to clarify the effect of experimentally observed excessive V atoms, the magnetic and electronic properties of the V-doped VI3 structures are analyzed. Our findings indicate that partially doped VI3 structures display FM ground state while the fully-doped structure exhibits AFM ground state. The fully-doped monolayer VI3 is found to be a semiconductor with a relatively larger band gap than its pristine structure. In addition, strain-dependent electronic and magnetic properties of fully- and partially-doped VI3 structures reveal that pristine monolayer displays a FM-to-AFM phase transition with robust semiconducting nature for 5% of compressive strain, while fully-doped monolayer VI3 structure possesses AFM-to-FM semiconducting transition at tensile strains larger than 4%. In contrast, the partially-doped VI3 monolayers are found to display robust FM ground state under biaxial strain. Its dopant and strain tunable electronic and magnetic nature makes monolayer VI3 a promising material for applications in nanoscale spintronic devices.