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: 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: We investigate the electron states and optical absorption in square- and hexagonal-shaped two-dimensional (2D) HgTe quantum dots and quantum rings in the presence of a perpendicular magnetic field. The electronic structure is modeled by means of the sp(3)d(5)s* tight-binding method within the nearest-neighbor approximation. Both bulklike and edge states appear in the energy spectrum. The bulklike states in quantum rings exhibit Aharonov-Bohm oscillations in magnetic field, whereas no such oscillations are found in quantum dots, which is ascribed to the different topology of the two systems. When magnetic field varies, all the edge states in square quantum dots appear as quasibands composed of almost fully flat levels, whereas some edge states in quantum rings are found to oscillate with magnetic field. However, the edge states in hexagonal quantum dots are localized like in rings. The absorption spectra of all the structures consist of numerous absorption lines, which substantially overlap even for small line broadening. The absorption lines in the infrared are found to originate from transitions between edge states. It is shown that the magnetic field can be used to efficiently tune the optical absorption of HgTe 2D quantum dot and quantum ring systems.

Abstract: The densities of phonon states of thin Sn films on InSb substrates are determined during different stages of the alpha-Sn to beta-Sn phase transition using nuclear inelastic x-ray scattering. The vibrational entropy and internal energy per atom as a function of temperature are obtained by numerical integration of the phonon density of states. The free energy as a function of temperature for the nanoscale samples is compared to the free energy obtained from ab initio calculations of bulk tin in the alpha-Sn and beta-Sn phase. In thin films this phase transition is governed by the interplay between the vibrational behavior of the film (the phase transition is driven by the vibrational entropy) and the stabilizing influence of the substrate (which depends on the film thickness). This brings a deeper understanding of the role of lattice vibrations in the phase transition of nanoscale Sn.

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: By performing density functional theory based first-principles calculations, the electronic, vibrational, elastic, and piezoelectric properties of two dynamically stable crystal phases of monolayer Janus MoSTe, namely 1H-MoSTe and 1T'-MoSTe, are investigated. Vibrational frequency analysis reveals that the other possible crystal structure, 1T-MoSTe, of this Janus monolayer does not exhibit dynamical stability. The 1H-MoSTe phase is found to be an indirect band-gap semiconductor while 1T'-MoSTe is predicted as small-gap semiconductor. Notably, in contrast to the direct band-gap nature of monolayers 1H-MoS2 and 1H-MoTe2, 1H-MoSTe is found to be an indirect gap semiconductor driven by the induced surface strains on each side of the structure. The calculated Raman spectrum of each structure shows unique character enabling us to clearly distinguish the stable crystal phases via Raman measurements. The systematic piezoelectric stress and strain coefficient analysis reveals that out-of-plane piezoelectricity appears in 1H-MoSTe and the noncentral symmetric 1T'-MoSTe has large piezoelectric coefficients. Static total-energy calculations show clearly that the formation of 1T'-MoSTe is feasible by using 1T'-MoTe2 as a basis monolayer. Therefore, we propose that the Janus MoSTe structure can be fabricated in two dynamically stable phases which possess unique electronic, dynamical, and piezoelectric properties.

Abstract: Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalco-genides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory-based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber-type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectralenergy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.

Abstract: Using density functional theory-based first-principles calculations, we investigate the stacking order dependence of the electronic and vibrational properties of HfS2-Ca(OH)(2) heterobilayer structures. It is shown that while the different stacking types exhibit similar electronic and optical properties, they are distinguishable from each other in terms of their vibrational properties. Our findings on the vibrational properties are the following: (i) from the interlayer shear (SM) and layer breathing (LBM) modes we are able to deduce the AB' stacking order, (ii) in addition, the AB' stacking type can also be identified via the phonon softening of E-g(I) and A(g)(III) modes which harden in the other two stacking types, and (iii) importantly, the ultrahigh frequency regime possesses distinctive properties from which we can distinguish between all stacking types. Moreover, the differences in optical and vibrational properties of various stacking types are driven by two physical effects, induced biaxial strain on the layers and the layer-layer interaction. Our results reveal that with both the phonon frequencies and corresponding activities, the Raman spectrum possesses distinctive properties for monitoring the stacking type in novel vertical heterostructures constructed by alkaline-earth-metal hydroxides.

Abstract: Superfluidity has recently been reported in double electron-hole bilayer graphene. The multiband nature of the bilayers is important because of the very small band gaps between conduction and valence bands. The long-range nature of the superfluid pairing interaction means that screening must be fully taken into account. We have carried out a systematic mean-field investigation that includes (i) contributions to screening from both intraband and interband excitations, (ii) the low-energy band structure of bilayer graphene with its small band gap and flattened Mexican-hat-like low-energy bands, (iii) the large density of states at the bottom of the bands, (iv) electron-hole pairing in the multibands, and (v) electron-hole pair transfers between the conduction and valence band condensates. We find that the superfluidity strongly modifies the intraband contributions to the screening, but that the interband contributions are unaffected. Unexpectedly, a net effect of the screening is to suppress Josephson-like pair transfers and to confine the superfluid pairing entirely to the conduction-band condensate even for very small band gaps, making the system behave similarly to a one-band superfluid.

Abstract: A four-band exciton Hamiltonian is constructed starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs). The angular part of the exciton wave function can be separated from the radial part, in the case of zero center of mass momentum excitons, by exploiting the eigenstates of the total exciton angular momentum operator with which the Hamiltonian commutes. We explain why this approach fails for excitons with finite center of mass momentum or in the presence of a perpendicular magnetic field and present an approximation to resolve this issue. We calculate the (binding) energy and average interparticle distance of different excited exciton states in different TMDs and compare these with results available in the literature. Remarkably, we find that the intervalley exciton ground state in the -/+ K valley has angular momentum j = +/- 1, which is due to the pseudospin of the separate particles. The exciton mass and the exciton Landau levels are calculated and we find that the degeneracy of exciton states with opposite relative angular momentum is altered by a magnetic field.

Abstract: Understanding the appearance of local electrical connectivity in liquid filled layered graphene oxide subjected to an external electric field is important to design electrically controlled smart permeable devices and also to gain insight into the physics behind electrical effects on confined water permeation. Motivated by recent experiments [K. G. Zhou et al. Nature (London) 559, 236 (2018)], we introduce a new model with random percolating paths for electrical connectivity in micron thick water filled layered graphene oxide, which mimics parallel resistors connected across the top and bottom electrodes. We find that a strong nonuniform radial electric field of the order similar to 10-50 mV/nm can be induced between layers depending on the current flow through the formed conducting paths. The maxima of the induced fields are not necessarily close to the electrodes and may be localized in the middle region of the layered material. The emergence of electrical connectivity and the associated electrical effects have a strong influence on the surrounding fluid in terms of ionization and wetting which subsequently determines the permeation properties.

Abstract: Devising magnetic nanostructures with spatially heterogeneous Dzyaloshinskii-Moriya interaction (DMI) is a promising pathway toward advanced confinement and control of magnetic skyrmions in potential devices. Here we discuss theoretically how a skyrmion interacts with a heterochiral interface using micromagnetic simulations and analytic arguments. We show that a heterochiral interface deflects the trajectory of ferromagnetic (FM) skyrmions, and that the extent of such deflection is tuned by the applied spin-polarized current and the difference in DMI across the interface. Further, we show that this deflection is characteristic of the FM skyrmion, and it is completely absent in the antiferromagnetic (AFM) case. In turn, we reveal that the AFM skyrmion achieves much higher velocities than its FM counterpart, yet experiences far stronger confinement in nanoengineered heterochiral tracks, which reinforces AFM skyrmions as a favorable choice for skyrmion-based devices.

Abstract: When slow ions incident on a surface are neutralized, the excess potential energy is passed on to an electron inside the surface, leading to emission of secondary electrons. The microscopic description of this process, as
well as the calculation of the secondary electron yield, is a challenging problem due to its complexity as well
as its sensitivity to surface properties. One of the first quantitative descriptions was articulated in the 1950s by
Hagstrum, who based his calculation on a parametrization of the density of states of the material. In this paper, we
present a model for calculating the secondary electron yield, derived from Hagstrum’s initial approach. We use
first-principles density functional theory calculations to acquire the necessary input and introduce the concept of
electron cascades to Hagstrum’s model in order to improve the calculated spectra, as well as remove its reliance
on fitting parameters. We apply our model to He+ and Ne+ ions incident on Ge(111) and Si(111) and obtain
yield spectra that match closely to the experimental results of Hagstrum.

Abstract: Interface effect in complex oxide thin-film heterostructures lies at the vanguard of current research to design technologically relevant functionality and explore emergent physical phenomena. While most of the previous works focus on the perovskite/perovskite heterostructures, the study of perovskite/brownmillerite interfaces remains in its infancy. Here, we investigate spontaneously stabilized perovskite-ferromagnet (SrCoO3-delta)/brownmillerite-antiferromagnet (SrCoO2.5) bilayer with T-N > T-C and discover an unconventional interfacial magnetic exchange bias effect. From magnetometry investigations, it is rationalized that the observed effect stems from the interfacial ferromagnet/antiferromagnet coupling. The possibility for coupled ferromagnet/spin-glass interface engendering such effect is ruled out. Strikingly, a finite coercive field persists in the paramagnetic state of SrCoO3-delta,whereas the exchange bias field vanishes at T-C . We conjecture the observed effect to be due to the effective external quenched staggered field provided by the antiferromagnetic layer for the ferromagnetic spins at the interface. Our results not only unveil a paradigm to tailor the interfacial magnetic properties in oxide heterostructures without altering the cations at the interface, but also provide a purview to delve into the fundamental aspects of exchange bias in such unusual systems, paving a big step forward in thin-film magnetism.

Abstract: The classical-map hypernetted-chain (CHNC) scheme, developed for treating fermion fluids at strong coupling and at finite temperatures, is applied to electron-electron and electron-hole double quantum wells. The pair-distribution functions and the local field factors needed in linear-response theory are determined for a range of temperatures, carrier densities, and barrier widths typical for experimental double-quantum-well systems in GaAs-GaAlAs. For electron-hole double quantum wells, a large enhancement in the pair-distribution functions is found for small carrier separations. The CHNC equations for electron-hole systems no longer hold at low densities where bound-state formation occurs.