Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ab initio study of ${(\nu_\ell,\ell^-)}$ and ${(\overline{\nu}_\ell,\ell^+)}$ inclusive scattering in $^{12}\hbox{C}$: confronting the MiniBooNE and $\hbox{T2K CCQE}$ data

We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections, measured on $^{12}$C by the MiniBooNE and T2K collaborations in the charged-current quasielastic regime. The calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one-and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with $\Lambda_A = 1.0$ or $1.15$ GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Experimental and computational studies of Criegee intermediate reactions with NH 3 and CH 3 NH 2

Ammonia and amines are emitted into the troposphere by various natural and anthropogenic sources, where they have a significant role in aerosol formation. Here, we explore the significance of their removal by reaction with Criegee intermediates, which are produced in the troposphere by ozonolysis of alkenes. Rate coefficients for the reactions of two representative Criegee intermediates, formaldehyde oxide (CH 2 OO) and acetone oxide ((CH 3 ) 2 COO) with NH 3 and CH 3 NH 2 were measured using cavity ring-down spectroscopy. Temperature-dependent rate coefficients, k(CH 2 OO + NH 3 ) = (3.1 ± 0.5) × 10 -20 T 2 exp(1011 ± 48/T) cm 3 s -1 and k(CH 2 OO + CH 3 NH 2 ) = (5 ± 2) × 10 -19 T 2 exp(1384 ± 96/T) cm 3 s -1 were obtained in the 240 to 320 K range. Both the reactions of CH 2 OO were found to be independent of pressure in the 10 to 100 Torr (N 2 ) range, and average rate coefficients k(CH 2 OO + NH 3 ) = (8.4 ± 1.2) × 10 -14 cm 3 s -1 and k(CH 2 OO + CH 3 NH 2 ) = (5.6 ± 0.4) × 10 -12 cm 3 s -1 were deduced at 293 K. An upper limit of ≤2.7 × 10 -15 cm 3 s -1 was estimated for the rate coefficient of the (CH 3 ) 2 COO + NH3 reaction. Complementary measurements were performed with mass spectrometry using synchrotron radiation photoionization giving k(CH 2 OO + CH 3 NH 2 ) = (4.3 ± 0.5) × 10-12 cm3 s-1 at 298 K and 4 Torr (He). Photoionization mass spectra indicated production of NH 2 CH 2 OOH and CH 3 N(H)CH 2 OOH functionalized organic hydroperoxide adducts from CH 2 OO + NH 3 and CH 2 OO + CH 3 NH 2 reactions, respectively. Ab initio calculations performed at the CCSD(T)(F12*)/cc-pVQZ-F12//CCSD(T)(F12*)/cc-pVDZ-F12 level of theory predicted pre-reactive complex formation, consistent with previous studies. Master equation simulations of the experimental data using the ab initio computed structures identified submerged barrier heights of -2.1 ± 0.1 kJ mol -1 and -22.4 ± 0.2 kJ mol -1 for the CH 2 OO + NH 3 and CH 2 OO + CH 3 NH 2 reactions, respectively. The reactions of NH3 and CH3NH2 with CH 2 OO are not expected to compete with its removal by reaction with (H 2 O) 2 in the troposphere. Similarly, losses of NH 3 and CH 3 NH 2 by reaction with Criegee intermediates will be insignificant compared with reactions with OH radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Naturally-meaningful and efficient descriptors: machine learning of material properties based on robust one-shot ab initio descriptors

Establishing a data-driven pipeline for the discovery of novel materials requires the engineering of material features that can be feasibly calculated and can be applied to predict a material’s target properties. Here we propose a new class of descriptors for describing crystal structures, which we term Robust One-Shot Ab initio (ROSA) descriptors. ROSA is computationally cheap and is shown to accurately predict a range of material properties. These simple and intuitive class of descriptors are generated from the energetics of a material at a low level of theory using an incomplete ab initio calculation. We demonstrate how the incorporation of ROSA descriptors in ML-based property prediction leads to accurate predictions over a wide range of crystals, amorphized crystals, metal–organic frameworks and molecules. We believe that the low computational cost and ease of use of these descriptors will significantly improve ML-based predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Revealing Local Structures through Machine-Learning-Fused Multimodal Spectroscopy

Atomistic structures of materials offer valuable insights into their functionality. Determining these structures remains a fundamental challenge in materials science, especially for systems with defects. While both experimental and computational methods exist, each has limitations in resolving nanoscale structures. Core-level spectroscopies, such as X-ray absorption (XAS) or electron energy-loss spectroscopies (EELS), have been used to determine the local bonding environment and structure of materials. Recently, machine learning (ML) methods have been applied to extract structural and bonding information from XAS/EELS data. However, frameworks relying solely on a single data stream, defined as characterization data derived from a single element using one technique, are often insufficient because multiple local environments can yield similar spectral features, making it challenging to differentiate between competing structural hypotheses. Here, in this work, we address this challenge by integrating multimodal ab initio simulations, experimental data acquisition, and ML techniques for structure characterization. Our goal is to determine local structures and properties using EELS and XAS data from multiple elements and edges. To showcase our approach, we use various lithium nickel manganese cobalt (NMC) oxide compounds which are used for lithium ion batteries, including those with oxygen vacancies and antisite defects, as the sample material system. We successfully inferred local element content, ranging from lithium to transition metals, with quantitative agreement with experimental data. Beyond local element inference, we find that ML model based on multimodal spectroscopic data is able to determine whether local defects such as oxygen vacancy and antisites are present, a task which is impossible for single mode spectra or other experimental techniques. Furthermore, our framework is able to provide physical interpretability, bridging spectroscopy with the local atomic and electronic structures.

battery↗
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2025

Ca X ML: Chemistry‐informed machine learning explains mutual changes between protein conformations and calcium ions in calcium‐binding proteins using structural and topological features

Proteins' flexibility is a feature in communicating changes in cell signaling instigated by binding with secondary messengers, such as calcium ions, associated with the coordination of muscle contraction, neurotransmitter release, and gene expression. When binding with the disordered parts of a protein, calcium ions must balance their charge states with the shape of calcium-binding proteins and their versatile pool of partners depending on the circumstances they transmit. Accurately determining the ionic charges of those ions is essential for understanding their role in such processes. However, it is unclear whether the limited experimental data available can be effectively used to train models to accurately predict the charges of calcium-binding protein variants. Here, we developed a chemistry-informed, machine-learning algorithm that implements a game theoretic approach to explain the output of a machine-learning model without the prerequisite of an excessively large database for high-performance prediction of atomic charges. We used the ab initio electronic structure data representing calcium ions and the structures of the disordered segments of calcium-binding peptides with surrounding water molecules to train several explainable models. Network theory was used to extract the topological features of atomic interactions in the structurally complex data dictated by the coordination chemistry of a calcium ion, a potent indicator of its charge state in protein. Our design created a computational tool of Ca X ML, which provided a framework of explainable machine learning model to annotate ionic charges of calcium ions in calcium-binding proteins in response to the chemical changes in an environment. Our framework will provide new insights into protein design for engineering functionality based on the limited size of scientific data in a genome space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Physics in the Machine: Integrating Physical Knowledge in Autonomous Phase-Mapping

Application of artificial intelligence (AI), and more specifically machine learning, to the physical sciences has expanded significantly over the past decades. In particular, science-informed AI, also known as scientific AI or inductive bias AI, has grown from a focus on data analysis to now controlling experiment design, simulation, execution and analysis in closed-loop autonomous systems. The CAMEO (closed-loop autonomous materials exploration and optimization) algorithm employs scientific AI to address two tasks: learning a material system’s composition-structure relationship and identifying materials compositions with optimal functional properties. By integrating these, accelerated materials screening across compositional phase diagrams was demonstrated, resulting in the discovery of a best-in-class phase change memory material. Key to this success is the ability to guide subsequent measurements to maximize knowledge of the composition-structure relationship, or phase map. In this work we investigate the benefits of incorporating varying levels of prior physical knowledge into CAMEO’s autonomous phase-mapping. This includes the use of ab-initio phase boundary data from the AFLOW repositories, which has been shown to optimize CAMEO’s search when used as a prior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Elucidating the Chemical Dynamics of the Elementary Reactions of the 1-Propynyl Radical (CH 3 CC; X 2 A 1 ) with Methylacetylene (H 3 CCCH; X 1 A 1 ) and Allene (H 2 CCCH 2 ; X 1 A 1 )

The reactions of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with two C 3 H 4 isomers, methylacetylene (H 3 CCCH; X 1 A 1 ) and allene (H 2 CCCH 2 ; X 1 A 1 ), along with their (partially) deuterated counterparts were explored at collision energies of 37 kJ mol –1 , exploiting crossed molecular beams to unravel the chemical reaction dynamics to synthesize distinct C 6 H 6 isomers under single collision conditions. The forward convolution fitting of the laboratory data along with ab initio and statistical calculations revealed that both reactions have no entrance barrier, proceed via indirect (complex-forming) reaction dynamics involving C 6 H 7 intermediates with life times longer than their rotation period(s), and are initiated by the addition of the 1-propynyl radical with its radical center to the π-electron density of the unsaturated hydrocarbon at the terminal carbon atoms of methylacetylene (C1) and allene (C1/C3). In the methylacetylene system, the initial collision complexes undergo unimolecular decomposition via tight exit transition states by atomic hydrogen loss, forming dimethyldiacetylene (CH 3 CCCCCH 3 ) and 1-propynylallene (H 3 CCCHCCCH 2 ) in overall exoergic reactions (123 and 98 kJ mol –1 ) with a branching ratio of 9.4 ± 0.1; the methyl group of the 1-propynyl reactant acts solely as a spectator. On the other hand, in the allene system, our experimental data exhibit the formation of the fulvene (c-C 5 H 4 CH 2 ) isomer via a six-step reaction sequence with two higher energy isomers—hexa-1,2-dien-4-yne (H 2 CCCHCCCH 3 ) and hexa-1,4-diyne (HCCCH 2 CCCH 3 )—also predicted to be formed based on our statistical calculations. Here, the pathway to fulvene advocates that, in the allene–1-propynyl system, the methyl group of the 1-propynyl reactant is actively engaged in the reaction mechanism to form fulvene. Because both reactions are barrierless and exoergic and all transition states are located below the energy of the separated reactants, the hydrogen-deficient C 6 H 6 isomers identified in our investigation are predicted to be synthesized in low-temperature environments, such as in hydrocarbon-rich atmospheres of planets and their moons such as Titan along with cold molecular clouds such as Taurus Molecular Cloud-1.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Competing Magnetic Interactions and the Role of Unpaired 4 f Electrons in Oxygen-Deficient Perovskites Ba 3 R Fe 2 O 7.5 ( R = Y, Dy)

Oxygen-deficient perovskite compounds with the general formula Ba 3 RFe 2 O 7.5 present a good opportunity to study competing magnetic interactions between Fe 3+ 3d cations with and without the involvement of unpaired 4f electrons on R 3+ cations. From analysis of neutron powder diffraction data, complemented by ab initio density functional theory calculations, we determined the magnetic ground states when R 3+ = Y 3+ (non-magnetic) and Dy 3+ (4f 9 ). They both adopt complex long-range ordered antiferromagnetic structures below T N = 6.6 and 14.5 K, respectively, with the same magnetic space group C a 2/c (BNS #15.91). However, the dominant influence of f-electron magnetism is clear in temperature dependence and differences between the size of the ordered moments on the two crystallographically independent Fe sites, one of which is enhanced by R–O–Fe superexchange in the Dy compound, while the other is frustrated by it. We report the Dy compound also shows evidence for temperature- and field-dependent transitions with hysteresis, indicating a field-induced ferromagnetic component below TN.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Understanding Strain and Failure of a Knot in Polyethylene Using Molecular Dynamics with Machine-Learned Potentials

A neural network potential (NNP) has been developed by fitting to ab initio electronic structure data on hydrocarbons and is used to study failure of linear and knotted polyethylene (PE) chains. A linear PE chain must be highly strained before breaking as the stress is equally distributed across the chain. In contrast, the stress in a PE chain with a 31 or overhand knot, accumulates at the knot’s entrance/exit. We find the strain energy is greatest when the bond length and angle are strained simultaneously, and that the knot weakens the chain by increasing the variance of the C–C–C angle, thereby allowing rupture at lower bond strains. Here, we extend our analysis to both 51 and 52 knots and find that both break at the entrance/exit of a loop. Notably, molecular scale PE knots exhibit many of the same characteristics as knots in a macroscopic rope, with stick–slip phenomena upon tightening and similar points of failure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Isolation and electronic structures of derivatized manganocene, ferrocene and cobaltocene anions

The discovery of ferrocene nearly 70 years ago marked the genesis of metallocene chemistry. Although the ferrocenium cation was discovered soon afterwards, a derivatized ferrocenium dication was only isolated in 2016 and the monoanion of ferrocene has only been observed in low-temperature electrochemical studies. Here we report the isolation of a derivatized ferrocene anion in the solid state as part of an isostructural family of 3d metallocenates, which consist of anionic complexes of a metal centre (manganese, iron or cobalt) sandwiched between two bulky Cp ttt ligands (where Cp ttt is {1,2,4-C 5 H 2 t Bu 3 }). These thermally and air-sensitive complexes decompose rapidly above –30 °C; however, we were able to characterize all metallocenates by a wide range of physical techniques and ab initio calculations. These data have allowed us to map the electronic structures of this metallocenate family, including an unexpected high-spin S = 3/2 ground state for the 19e– derivatized ferrocene anion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Ultrafast flow chemistry for the acid-catalyzed conversion of fructose

Biomass conversion to chemicals and fuels has been epitomized by batch processing and long reaction times that limit the economic viability of future biorefineries and prevent modular manufacturing near the source-production site that is necessary due to the biomass’ large water content. Herein we introduce a continuous flow microreactor to enable process intensification and one thousand-fold device miniaturization for remote chemical processing. We demonstrate the approach in HCl-catalyzed fructose dehydration in water. We characterize and quantify the mixing using laser induced fluorescence (LIF) of fluorescent dyes and particle image velocimetry (PIV) of fluorescent microspheres. The estimated mixing times ranged from 0.03 to 4.8 s for residence times from 1 to 120 s, respectively. Moreover, a curved channel geometry induces secondary (Dean) vortices that produce a three-fold increase in mixing by diffusion through advection along the cross-sectional direction of the channel, creating dissipative chaotic mixing under laminar flow. We then obtain for the first time isothermal kinetics of fructose dehydration at short contact times and high temperatures and compare to a previously published, hybrid ab initio and data-driven kinetic model. Rigorous model optimization is carried out and an optimal HMF yield of 54% is attained in a single aqueous phase. High temperatures and low pH improve the reactor throughput by increasing the HMF space-time yield. We demonstrate the highest ever productivity (HMF yield per unit time). Importantly, this can be achieved at ~100% fructose conversion, eliminating the need for recycling to improve economics. Furthermore, optimization of the reactor dimensions and the specific energy dissipation rate results in ~60% pressure-drop related energy savings. We describe a methodology to scale-up microreactors to millimeter size and increase productivity 16-fold under optimal conditions. The energy efficiency and cost analysis for HMF production from corn stover of a typical corn farm indicate an optimal strategy for numbering-up millimeter size reactors and, importantly, they underscore the potential feasibility of producing bioproducts in a distributed manner.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin–orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3 (A) leading to the NH(X 3 Σ – , a 1 Δ) + H 2 channel

Here we simulate the photodissociation of NH 3 originating from its first excited singlet state S 1 into the NH 2 + H (radical) and NH + H 2 (molecular) channels. The states considered are the ground singlet state S 0 , the first excited singlet state S 1 and the lowest-lying triplet state T 1 , which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S 0 , S 1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T 1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H 2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH 2 + H channel. The NH(X 3 Σ_) radical produced through the intersystem crossing from S 0 to T 1 is rare (~0.2%) compared to NH(a 1 Δ) due to the process being spin forbidden.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Elasticity of two-dimensional ferroelectrics across their paraelectric phase transformation

The mechanical behavior of two-dimensional (2D) materials across 2D phase changes is unknown, and the finite temperature (T) elasticity of paradigmatic SnSe monolayers—ferroelectric 2D materials turning paraelectric as their unit cell turns from a rectangle into a square—is described here in a progressive manner. To begin with, their zero–T total energy landscape gives way to (Boltzmann-like) averages from which the elastic behavior is determined. Furthermore, these estimates are complemented with results from the strain-fluctuation method, which employs the energy landscape or ab initio molecular dynamics data. All approaches capture the coalescence of elastic moduli < C 11 (T) > = < C 22 (T) > due to the structural transformation. The broad evolution and sudden changes of elastic parameters < C 11 (T) >, < C 22 (T) >, and < C 12 (T) > of these atomically thin phase-change membranes establishes a heretofore overlooked connection among 2D materials and soft matter.

2-dimensional systems↗