Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results.

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Induction and Ferroelectric Switching of Flux Closure Domains in Strained PbTiO3 with Neural Network Quantum Molecular Dynamics.

We have developed an extension of the Neural Network Quantum Molecular Dynamics (NNQMD) simulation method to incorporate electric-field dynamics based on Born effective charge (BEC), called NNQMD-BEC. We first validate NNQMD-BEC for the switching mechanisms of archetypal ferroelectric PbTiO3 bulk crystal and 180° domain walls (DWs). NNQMD-BEC simulations correctly describe the nucleation-and-growth mechanism during DW switching. In triaxially strained PbTiO3 with strain conditions commonly seen in many superlattice configurations, we find that flux-closure texture can be induced with application of an electric field perpendicular to the original polarization direction. Upon field reversal, the flux-closure texture switches via a pair of transient vortices as the intermediate state, indicating an energy-efficient switching pathway. Our NNQMD-BEC method provides a theoretical guidance to study electro-mechano effects with existing machine learning force fields using a simple BEC extension, which will be relevant for engineering applications such as field-controlled switching in mechanically strained ferroelectric devices.

Dielectric Constant of Liquid Water Determined with Neural Network Quantum Molecular Dynamics.

The static dielectric constant ϵ_{0} and its temperature dependence for liquid water is investigated using neural network quantum molecular dynamics (NNQMD). We compute the exact dielectric constant in canonical ensemble from NNQMD trajectories using fluctuations in macroscopic polarization computed from maximally localized Wannier functions (MLWF). Two deep neural networks are constructed. The first, NNQMD, is trained on QMD configurations for liquid water under a variety of temperature and density conditions to learn potential energy surface and forces and then perform molecular dynamics simulations. The second network, NNMLWF, is trained to predict locations of MLWF of individual molecules using the atomic configurations from NNQMD. Training data for both the neural networks is produced using a highly accurate quantum-mechanical method, DFT-SCAN that yields an excellent description of liquid water. We produce 280×10^{6} configurations of water at 7 temperatures using NNQMD and predict MLWF centers using NNMLWF to compute the polarization fluctuations. The length of trajectories needed for a converged value of the dielectric constant at 0°C is found to be 20 ns (40×10^{6} configurations with 0.5 fs time step). The computed dielectric constants for 0, 15, 30, 45, 60, 75, and 90°C are in good agreement with experiments. Our scalable scheme to compute dielectric constants with quantum accuracy is also applicable to other polar molecular liquids.

Electric-field-induced crossover of polarization reversal mechanisms in Al1-xScxN ferroelectrics.

Scandium-doped aluminum nitride, Al1-xScxN, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al-Sc-N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al1-xScxN to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms-inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

Growth Kinetics and Atomistic Mechanisms of Native Oxidation of ZrSxSe2-x and MoS2 Crystals.

A thorough understanding of native oxides is essential for designing semiconductor devices. Here, we report a study of the rate and mechanisms of spontaneous oxidation of bulk single crystals of ZrSxSe2-x alloys and MoS2. ZrSxSe2-x alloys oxidize rapidly, and the oxidation rate increases with Se content. Oxidation of basal surfaces is initiated by favorable O2 adsorption and proceeds by a mechanism of Zr-O bond switching, that collapses the van der Waals gaps, and is facilitated by progressive redox transitions of the chalcogen. The rate-limiting process is the formation and out-diffusion of SO2. In contrast, MoS2 basal surfaces are stable due to unfavorable oxygen adsorption. Our results provide insight and quantitative guidance for designing and processing semiconductor devices based on ZrSxSe2-x and MoS2 and identify the atomistic-scale mechanisms of bonding and phase transformations in layered materials with competing anions.

Phonon-Suppressed Auger Scattering of Charge Carriers in Defective Two-Dimensional Transition Metal Dichalcogenides.

Two-dimensional transition metal dichalcogenides (TMDs) draw strong interest in materials science, with applications in optoelectronics and many other fields. Good performance requires high carrier concentrations and long lifetimes. However, high concentrations accelerate energy exchange between charged particles by Auger-type processes, especially in TMDs where many-body interactions are strong, thus facilitating carrier trapping. We report time-resolved optical pump-THz probe measurements of carrier lifetimes as a function of carrier density. Surprisingly, the lifetime reduction with increased density is very weak. It decreases only by 20% when we increase the pump fluence 100 times. This unexpected feature of the Auger process is rationalized by our time-domain ab initio simulations. The simulations show that phonon-driven trapping competes successfully with the Auger process. On the one hand, trap states are relatively close to band edges, and phonons accommodate efficiently the electronic energy during the trapping. On the other hand, trap states localize around defects, and the overlap of trapped and free carriers is small, decreasing carrier-carrier interactions. At low carrier densities, phonons provide the main charge trapping mechanism, decreasing carrier lifetimes compared to defect-free samples. At high carrier densities, phonons suppress Auger processes and lower the dependence of the trapping rate on carrier density. Our results provide theoretical insights into the diverse roles played by phonons and Auger processes in TMDs and generate guidelines for defect engineering to improve device performance at high carrier densities.

Two-Dimensional Lateral Epitaxy of 2H (MoSe2)-1T' (ReSe2) Phases.

Two-dimensional (2D) transition metal dichalcogenide (TMDC) heterostructures have been proposed as potential candidates for a variety of applications like quantum computing, neuromorphic computing, solar cells, and flexible field effective transistors. The 2D TMDC heterostructures at the present stage face difficulties being implemented in these applications because of lack of large and sharp heterostructure interfaces. Herein, we address this problem via a CVD technique to grow thermodynamically stable heterostructure of 2H/1T' MoSe2-ReSe2 using conventional transition metal phase diagrams as a reference. We demonstrate how the thermodynamics of mixing in the MoReSe2 system during CVD growth dictates the formation of atomically sharp interfaces between MoSe2 and ReSe2, which can be confirmed by high-resolution scanning transmission electron microscopy imaging, revealing zigzag selenium-terminated interface between the epitaxial 2H and 1T' lattices. Our work provides useful insights for understanding the stability of 2D heterostructures and interfaces between chemically, structurally, and electronically different phases.

Optical Control of Non-Equilibrium Phonon Dynamics.

The light-induced selective population of short-lived far-from-equilibrium vibration modes is a promising approach for controlling ultrafast and irreversible structural changes in functional nanomaterials. However, this requires a detailed understanding of the dynamics and evolution of these phonon modes and their coupling to the excited-state electronic structure. Here, we combine femtosecond mega-electronvolt electron diffraction experiments on a prototypical layered material, MoTe2, with non-adiabatic quantum molecular dynamics simulations and ab initio electronic structure calculations to show how non-radiative energy relaxation pathways for excited electrons can be tuned by controlling the optical excitation energy. We show how the dominant intravalley and intervalley scattering mechanisms for hot and band-edge electrons leads to markedly different transient phonon populations evident in electron diffraction patterns. This understanding of how tuning optical excitations affect phonon populations and atomic motion is critical for efficiently controlling light-induced structural transitions of optoelectronic devices.

Electronic Origin of Optically-Induced Sub-Picosecond Lattice Dynamics in MoSe2 Monolayer.

Atomically thin layers of transition metal dichalcogenide (TMDC) semiconductors exhibit outstanding electronic and optical properties, with numerous applications such as valleytronics. While unusually rapid and efficient transfer of photoexcitation energy to atomic vibrations was found in recent experiments, its electronic origin remains unknown. Here, we study the lattice dynamics induced by electronic excitation in a model TMDC monolayer, MoSe2, using nonadiabatic quantum molecular dynamics simulations. Simulation results show sub-picosecond disordering of the lattice upon photoexcitation, as measured by the Debye-Waller factor, as well as increasing disorder for higher densities of photogenerated electron-hole pairs. Detailed analysis shows that the rapid, photoinduced lattice dynamics are due to phonon-mode softening, which in turn arises from electronic Fermi surface nesting. Such mechanistic understanding can help guide optical control of material properties for functionalizing TMDC layers, enabling emerging applications such as phase change memories and neuromorphic computing.

Computational Synthesis of MoS2 Layers by Reactive Molecular Dynamics Simulations: Initial Sulfidation of MoO3 Surfaces.

Transition metal dichalcogenides (TMDC) like MoS2 are promising candidates for next-generation electric and optoelectronic devices. These TMDC monolayers are typically synthesized by chemical vapor deposition (CVD). However, despite significant amount of empirical work on this CVD growth of monolayered crystals, neither experiment nor theory has been able to decipher mechanisms of selection rules for different growth scenarios, or make predictions of optimized environmental parameters and growth factors. Here, we present an atomic-scale mechanistic analysis of the initial sulfidation process on MoO3 surfaces using first-principles-informed ReaxFF reactive molecular dynamics (RMD) simulations. We identify a three-step reaction process associated with synthesis of the MoS2 samples from MoO3 and S2 precursors: O2 evolution and self-reduction of the MoO3 surface; SO/SO2 formation and S2-assisted reduction; and sulfidation of the reduced surface and Mo-S bond formation. These atomic processes occurring during early stage MoS2 synthesis, which are consistent with experimental observations and existing theoretical literature, provide valuable input for guided rational synthesis of MoS2 and other TMDC crystals by the CVD process.

Quantifying the origin of inter-adsorbate interactions on reactive surfaces for catalyst screening and design.

The adsorption energy of reactant molecules and reaction intermediates is one of the key descriptors of catalytic activity of surfaces and is commonly used as a metric in screening materials for design of heterogeneous catalysts. The efficacy of such screening schemes depends on the accuracy of calculated adsorption energies under reaction conditions. These adsorption energies can depend strongly on interactions between adsorbed molecules in the adlayer. However, these interactions are typically not accounted for in screening procedures that use DFT-based zero-coverage adsorption energies. Identifying the physical mechanisms behind these interactions is essential to model realistic catalyst surfaces under reaction conditions and to understand the dependence of adsorption energies on reaction parameters like surface strain and composition. This article describes a method to quantitatively resolve the observed inter-adsorbate interactions into various direct adsorbate-adsorbate interactions (i.e. Coulombic and steric) and surface-mediated interactions (i.e. adsorbate-induced surface relaxation and change in electronic structure) by combining density functional theory and cluster-expansion calculations of coverage-dependent adsorption energies. The approach is implemented on a model catalyst surface of FeS2(100) reacting with H2S molecules. We find that the adsorption energy of H2S molecules can be affected by over 0.55 eV by the repulsive inter-adsorbate interactions caused primarily by the adsorbate-induced changes to the electronic structure of the FeS2 surface. These interactions also show a strong monotonic dependence on surface strain, being three times stronger on compressively strained surfaces than on surfaces under tensile strain. The large magnitude of inter-adsorbate interactions as well as their strong dependence on lattice strain demonstrate the need for using coverage-dependent adsorption energies for more accurate screening, for example for strained catalytic systems like core-shell and overlayer structures.

Magnetic diffusion anomaly at the Néel temperature of pyrrhotite, Fe(1-x)S.

Cation diffusion is an important rate-limiting process in the growth of pyrrhotite (Fe1-xS) in passivating films on steels exposed to sulfidic environments, and for proposed synthetic applications of Fe1-xS, for example single-phase magnetic switching devices. Above the Néel temperature TN of 315 °C, where Fe1-xS is paramagnetic and structurally disordered, iron self-diffusivity *DFe predictably follows a standard, established Arrhenius law with temperature. However, we report (57)Fe tracer diffusion measurements below TN, obtained using secondary ion mass spectrometry (SIMS), that demonstrate a 100-fold reduction in diffusion coefficient as compared to the extrapolated, paramagnetic Arrhenius trend at 150 °C. The results can be described by a magnetic diffusion anomaly, where the vacancy migration energy for the spontaneously-magnetized cation sublattice is increased by approximately 40% over the paramagnetic state. These constitute the first set of consistent diffusivity data obtained in magnetic pyrrhotite, allowing more accurate prediction of pyrrhotite growth rates and determination of magnetic properties for synthetic devices.

Surface Transfer Doping in MoO3-x/Hydrogenated Diamond Heterostructure.

Surface transfer doping is proposed to be a potential solution for doping diamond, which is hard to dope for applications in high-power electronics. While MoO3 is found to be an effective surface electron acceptor for hydrogen-terminated diamond with a negative electron affinity, the effects of commonly existing oxygen vacancies remain elusive. We have performed reactive molecular dynamics simulations to study the deposition of MoO3-x on a hydrogenated diamond (111) surface and used first-principles calculations based on density functional theory to investigate the electronic structures and charge transfer mechanisms. We find that MoO3-x is an effective surface electron acceptor and the spatial extent of doped holes in hydrogenated diamond is extended, promoting excellent transport properties. Charge transfer is found to monotonically decrease with the level of oxygen vacancy, providing guidance for engineering of the surface transfer doping process.

Inelastic Neutron Scattering Study of Phonon Density of States of Iodine Oxides and First-Principles Calculations.

Iodine oxides I2Oy (y = 4, 5, 6) crystallize into atypical structures that fall between molecular- and framework-base types and exhibit high reactivity in an ambient environment, a property highly desired in the so-called "agent defeat materials". Inelastic neutron scattering experiments were performed to determine the phonon density of states of the newly synthesized I2O5 and I2O6 samples. First-principles calculations were carried out for I2O4, I2O5, and I2O6 to predict their thermodynamic properties and phonon density of states. Comparison of the INS data with the Raman and infrared measurements as well as the first-principles calculations sheds light on their distinctive, anisotropic thermomechanical properties.

Squishing Skyrmions: Symmetry-Guided Dynamic Transformation of Polar Topologies Under Compression.

Mechanical controllability of recently discovered topological defects (e.g., skyrmions) in ferroelectric materials is of interest for the development of ultralow-power mechano-electronics that are protected against thermal noise. However, fundamental understanding is hindered by the "multiscale quantum challenge" to describe topological switching encompassing large spatiotemporal scales with quantum mechanical accuracy. Here, we overcome this challenge by developing a machine-learning-based multiscale simulation framework─a hybrid neural network quantum molecular dynamics (NNQMD) and molecular mechanics (MM) method. For nanostructures composed of SrTiO3 and PbTiO3, we find how the symmetry of mechanical loading essentially controls polar topological switching. We find under symmetry-breaking uniaxial compression a squishing-to-annihilation pathway versus formation of a topological composite named skyrmionium under symmetry-preserving isotropic compression. The distinct pathways are explained in terms of the underlying materials' elasticity and symmetry, as well as the Landau-Lifshitz-Kittel scaling law. Such rational control of ferroelectric topologies will likely facilitate exploration of the rich ferroelectric "topotronics" design space.

Hydrogen Bonding in Liquid Ammonia.

The nature of hydrogen bonding in condensed ammonia phases, liquid and crystalline ammonia has been a topic of much investigation. Here, we use quantum molecular dynamics simulations to investigate hydrogen bond structure and lifetimes in two ammonia phases: liquid ammonia and crystalline ammonia-I. Unlike liquid water, which has two covalently bonded hydrogen and two hydrogen bonds per oxygen atom, each nitrogen atom in liquid ammonia is found to have only one hydrogen bond at 2.24 Å. The computed lifetime of the hydrogen bond is t ≅ 0.1 ps. In contrast to crystalline water-ice, we find that hydrogen bonding is practically nonexistent in crystalline ammonia-I.

Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics.

Ferroelectric materials exhibit a rich range of complex polar topologies, but their study under far-from-equilibrium optical excitation has been largely unexplored because of the difficulty in modeling the multiple spatiotemporal scales involved quantum-mechanically. To study optical excitation at spatiotemporal scales where these topologies emerge, we have performed multiscale excited-state neural network quantum molecular dynamics simulations that integrate quantum-mechanical description of electronic excitation and billion-atom machine learning molecular dynamics to describe ultrafast polarization control in an archetypal ferroelectric oxide, lead titanate. Far-from-equilibrium quantum simulations reveal a marked photo-induced change in the electronic energy landscape and resulting cross-over from ferroelectric to octahedral tilting topological dynamics within picoseconds. The coupling and frustration of these dynamics, in turn, create topological defects in the form of polar strings. The demonstrated nexus of multiscale quantum simulation and machine learning will boost not only the emerging field of ferroelectric topotronics but also broader optoelectronic applications.

Piezoelectricity across 2D Phase Boundaries.

Piezoelectricity in low-dimensional materials and metal-semiconductor junctions has attracted recent attention. Herein, a 2D in-plane metal-semiconductor junction made of multilayer 2H and 1T' phases of molybdenum(IV) telluride (MoTe2 ) is investigated. Strong piezoelectric response is observed using piezoresponse force microscopy at the 2H-1T' junction, despite that the multilayers of each individual phase are weakly piezoelectric. The experimental results and density functional theory calculations suggest that the amplified piezoelectric response observed at the junction is due to the charge transfer across the semiconducting and metallic junctions resulting in the formation of dipoles and excess charge density, allowing the engineering of piezoelectric response in atomically thin materials.

Lattice thermal transport in two-dimensional alloys and fractal heterostructures.

Engineering thermal transport in two dimensional materials, alloys and heterostructures is critical for the design of next-generation flexible optoelectronic and energy harvesting devices. Direct experimental characterization of lattice thermal conductivity in these ultra-thin systems is challenging and the impact of dopant atoms and hetero-phase interfaces, introduced unintentionally during synthesis or as part of deliberate material design, on thermal transport properties is not understood. Here, we use non-equilibrium molecular dynamics simulations to calculate lattice thermal conductivity of [Formula: see text] monolayer crystals including [Formula: see text] alloys with substitutional point defects, periodic [Formula: see text] heterostructures with characteristic length scales and scale-free fractal [Formula: see text] heterostructures. Each of these features has a distinct effect on phonon propagation in the crystal, which can be used to design fractal and periodic alloy structures with highly tunable thermal conductivities. This control over lattice thermal conductivity will enable applications ranging from thermal barriers to thermoelectrics.

Review of strategies toward the development of alloy two-dimensional (2D) transition metal dichalcogenides.

Atomically thin two-dimensional (2D) transition metal dichalcogenides (TMDCs) have attracted significant attention owing to their prosperity in material research. The inimitable features of TMDCs triggered the emerging applications in diverse areas. In this review, we focus on the tailored and engineering of the crystal lattice of TMDCs that finally enhance the efficiency of the material properties. We highlight several preparation techniques and recent advancements in compositional engineering of TMDCs structure. We summarize different approaches for TMDCs such as doping and alloying with different materials, alloying with other 2D metals, and scrutinize the technological potential of these methods. Beyond that, we also highlight the recent significant advancement in preparing 2D quasicrystals and alloying the 2D TMDCs with MAX phases. Finally, we highlight the future perspectives for crystal engineering in TMDC materials for structure stability, machine learning concept marge with materials, and their emerging applications.