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: In this study, reactivity of single-layer GaAs against Ge atoms is studied by means of ab initio density functional theory calculations. Firstly, it is shown that Ge atoms interact quite strongly with the GaAs layer which allows the formation of Ge islands while it hinders the growth of detached germanene monolayers. It is also predicted that adsorption of Ge atoms on GaAs single-layer lead to formation of two novel stable single-layer crystal structures, namely 1H-GaGeAs and 1H(A)-GaGeAs. Both the total energy optimizations and the calculated vibrational spectra indicate the dynamical stability of both single layer structures. Moreover, although both structures crystallize in 1H phase, 1H-GaGeAs and 1H(A)-GaGeAs exhibit distinctive vibrational features in their Raman spectra which is quite important for distinguishing the structures. In contrast to the semiconducting nature of single-layer GaAs, both polytypes of GaGeAs exhibit metallic behavior confirmed by the electronic band dispersions. Furthermore, the linear-elastic constants, in-plane stiffness and Poisson ratio, reveal the ultrasoft nature of the GaAs and GaGeAs structures and the rigidity of GaAs is found to be slightly enhanced via Ge adsorption. With their stable, ultra-thin and metallic properties, predicted single-layer GaGeAs structures can be promising candidates for nanoscale electronic and mechanical applications.

Abstract: We present KITE, a general purpose open-source tight-binding software for accurate real-space simulations of electronic structure and quantum transport properties of large-scale molecular and condensed systems with tens of billions of atomic orbitals (N similar to 10(10)). KITE's core is written in C++, with a versatile Python-based interface, and is fully optimized for shared memory multi-node CPU architectures, thus scalable, efficient and fast. At the core of KITE is a seamless spectral expansion of lattice Green's functions, which enables large-scale calculations of generic target functions with uniform convergence and fine control over energy resolution. Several functionalities are demonstrated, ranging from simulations of local density of states and photo-emission spectroscopy of disordered materials to large-scale computations of optical conductivity tensors and real-space wave-packet propagation in the presence of magneto-static fields and spin-orbit coupling. On-the-fly calculations of real-space Green's functions are carried out with an efficient domain decomposition technique, allowing KITE to achieve nearly ideal linear scaling in its multi-threading performance. Crystalline defects and disorder, including vacancies, adsorbates and charged impurity centres, can be easily set up with KITE's intuitive interface, paving the way to user-friendly large-scale quantum simulations of equilibrium and non-equilibrium properties of molecules, disordered crystals and heterostructures subject to a variety of perturbations and external conditions.

Abstract: Integration of graphene-based two-dimensional materials is essential for nanoelectronics applications. Using density-functional theory, we systematically investigate the electronic properties of vertically stacked graphene-graphitic carbon nitrides (GE/GCN). We also studied the covalently lateral stitched graphene-graphitic carbon nitrides (GE-GCN heterojunctions). The effects of perpendicular electric field on the electronic properties of six different heterostructures, i.e., (i) one layer of GE on top of a layer of CnNm with (n, m) = (3,1), (3,4), and (4,3) and (ii) three heterostructures CnNm/Cn'Nm', where (n, m) not equal (n', m') are elucidated. The most important calculated features are (i) the systems GE/C3N4, C3N/C3N4, GE-C3N, GE-C4N3, and C3N-C3N4 exhibit semiconducting characteristics having small band gaps of Delta(0)=20, 250, 100, 100, 80 meV, respectively while (ii) the systems GE/C4N3, C3N/C4N3, and C3N-C4N3 show ferromagnetic-metallic properties. In particular, we found that, in semiconducting heterostructures, the band gap increases nontrivially with increasing the absolute value of the applied perpendicular electric field. This work is useful for designing heterojunctions and heterostructures made of graphene and other two-dimensional materials such as those proposed in recent experiments [X. Liu and M. C. Hersam Sci. Adv. 5, 6444 (2019)].

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: 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: The experimental parameter ranges needed to generate superfluidity in optical and drag experiments in GaAs double quantum wells are determined using a formalism that includes self-consistent screening of the Coulomb pairing interaction in the presence of the superfluid. The very different electron and hole masses in GaAs make this a particularly interesting system for superfluidity with exotic superfluid phases predicted in the BCS-Bose-Einstein condensation crossover regime. We find that the density and temperature ranges for superfluidity cover the range for which optical experiments have observed indications of superfluidity but that existing drag experiments lie outside the superfluid range. We also show that, for samples with low mobility with no macroscopically connected superfluidity, if the superfluidity survives in randomly distributed localized pockets, standard quantum capacitance measurements could detect these pockets.

Abstract: In the present work, the Cu-poor and Sn-rich CZTS thin films were prepared in order to study the volatility of Sn with respect to other components. Thin film compositions were kept intentionally Sn-rich to understand the behaviour of loss and segregation of Sn during sulphurization. The homogeneous composition distribution in precursor thin films turns heterogeneous with a change in morphology after sulphurization. The inability of identifying nanoscale secondary phases in CZTS thin film by conventional analytical techniques such as XRD and Raman, can be fulfilled by employing HAADF-STEM analysis. XPS and HAADF-STEM analyses provide the quantification of nanoscale secondary phases across the thin film and surface, respectively. The volatility of Sn was revealed in the form of segregation in the middle layer of CZTS cross-sectional lamella rather than loss to annealing atmosphere. It was observed that among the cations of CZTS, Sn segregates more than Cu, while Zn segregates least. The nanoscale spurious phases were observed to vary across different regions in the sulphurized CZTS sample. The reactive annealing lead to grain growth and formation of grain boundary features in the CZTS thin films, where annealing significantly modifies the potential difference and band bending at grain boundaries with respect to intra-grains.

Abstract: Using first-principles calculations, we study the effect of embedding various atoms into the nanopore sites of both C6N6 and C6N8 monolayers. Our results indicate that the embedded atoms significantly affect the electronic and magnetic properties of C6N6 and C6N8 monolayers and lead to extraordinary and multifarious electronic properties, such as metallic, half-metallic, spin-glass semiconductor and dilute-magnetic semiconductor behaviour. Our results reveal that the H atom concentration dramatically affects the C6N6 monolayer. On increasing the H coverage, the impurity states also increase due to H atoms around the Fermi-level. C6N6 shows metallic character when the H atom concentration reaches 6.25%. Moreover, the effect of charge on the electronic properties of both Cr@C6N6 and C@C6N8 is also studied. Cr@C6N6 is a ferromagnetic metal with a magnetic moment of 2.40 mu(B), and when 0.2 electrons are added and removed, it remains a ferromagnetic metal with a magnetic moment of 2.57 and 2.77 mu(B), respectively. Interestingly, one can observe a semi-metal, in which the VBM and CBM in both spin channels touch each other near the Fermi-level. C@C6N8 is a semiconductor with a nontrivial band gap. When 0.2 electrons are removed, it remains metallic, and under excess electronic charge, it exhibits half-metallic behaviour.

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: A periodic spatial modulation, as created by a moire pattern, has been extensively studied with the view to engineer and tune the properties of graphene. Graphene encapsulated by hexagonal boron nitride (hBN) when slightly misaligned with the top and bottom hBN layers experiences two interfering moire patterns, resulting in a so-called supermoire (SM). This leads to a lattice and electronic spectrum reconstruction. A geometrical construction of the nonrelaxed SM patterns allows us to indicate qualitatively the induced changes in the electronic properties and to locate the SM features in the density of states and in the conductivity. To emphasize the effect of lattice relaxation, we report band gaps at all Dirac-like points in the hole doped part of the reconstructed spectrum, which are expected to be enhanced when including interaction effects. Our result is able to distinguish effects due to lattice relaxation and due to the interfering SM and provides a clear picture on the origin of recently experimentally observed effects in such trilayer heterostuctures.

Abstract: Indoor environments are considered an important contributor to external human exposure to halogenated flame retardants (HFRs) due to the large amounts of chemicals currently incorporated in indoor equipment and the time humans spend every day in indoor environments. In this chapter, the presence and use of novel brominated flame retardants (NBFRs), dechlorane plus (DPs), chlorinated organophosphorus flame retardants (Cl-PFRs) and chlorinated paraffins (CPs) in indoor dust, air and consumer products collected from different indoor microenvironments (homes, public indoor spaces, and vehicles) are discussed. While data on the concentrations of HFRs in indoor dust and air are widely available, figures are still scarce for consumer products, such as textiles and foams, furnishings, flooring, electric and electronic products and building materials. This knowledge gaps still represents the biggest obstacle in linking eventual sources of contamination to the presence and chemical patterns in indoor dust and air.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: This dataset comprises a large collection of atom probe tomography datasets collected from DIN 1.4970 alloy that was irradiated with Fe ions at different conditions. The DIN 1.4970 alloy is an austenitic stainless steel with 15 wt% Cr, 15 wt% Ni, a small addition of Ti. The full composition and characterization of our material can be found published elsewhere [1,2]. Some of our material was subjected to ageing heat treatments at different temperatures for different times. Small samples of our original material and aged material was irradiated at the Michigan Ion Beam Laboratory in 2017 with 4.5 MeV Fe ions up to 40 dpa at an average dose rate of 2×10−4 dpa/s. This was done at three different temperatures: 300, 450, and 600 ºC. Atom probe samples were made of the irradiated layers (approximately 1.5 micron deep) with focused ion beam and mounted on Microtip coupons. APT measurements took place on three CAMECA LEAP-HR systems located at CAES in Idaho Falls, USA (files beginning with R33), at Montanuniversität Leoben in Leoben, Austria (R21) and at Friedrich–Alexander University in Erlangen, Germany (R56).

Abstract: Electroporation is a non-thermal electromagnetic phenomenon widely used in medical diseases treatment. Different mathematical models of electroporation have been proposed in literature to study pore evolution in biological membranes. This paper presents a nonlinear dispersive multiphysic model of electroporation in irregular shaped biological cells in which the spatial and temporal evolution of the pores size is taken into account. The model solves Maxwell and asymptotic Smoluchowski equations and it describes the dielectric dispersion of cell media using a Debye-based relationship. Furthermore, the irregular cell shape has been modeled using the Gielis superformula. Taking into account the cell in mitosis phase, the electroporation process has been studied comparing the numerical results pertaining the model with variable pore radius with those in which the pore radius is supposed constant. The numerical analysis has been performed exposing the biological cell to a rectangular electric pulse having duration of 10 μs. The obtained numerical results highlight considerable differences between the two different models underling the need to include into the numerical algorithm the differential equation modeling the spatial and time evolution of the pores size.