SEARCH · Engineering Papers

Results for “Anharmonic lattice dynamics”

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.

SourceResearch connectionFrom yearThrough year

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

363 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Handbook of Molecular Beam Epitaxy of Oxide Materials: Epitaxial Complex Oxide Growth on Semiconductors

This chapter covers the epitaxy of complex oxides on common inorganic semiconductors. The chapter opens with an introduction, motivating the subject and highlighting key general challenges (Section 6.1). We then proceed to describe the general steps of the growth procedure in Section 6.2, which concludes with an overall discussion about trade-offs and expectation management, with an emphasis on the oxide–semiconductor interface. From there, the text describes oxide growth on the common semiconductors, starting with silicon (Section 6.3), where surface reactions and oxidation are the most challenging aspects. Oxide epitaxy on germanium is described in Section 6.4, which is a less challenging oxide growth scheme on semiconductors. Section 6.5 describes oxide epitaxy on gallium arsenide, one of the more challenging schemes, where interface stability and surface preparation pose key challenges. Finally, in Section 6.6, oxide epitaxy on gallium nitride is described, where the lattice mismatch takes the stage as the key challenge. Altogether, this chapter aims to provide the reader with the practical knowledge, tactics and strategies for oxide epitaxy on semiconductors.

Demkov, Alexander A [The University of Texas at Au↗

DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Quantifying the dislocation content of atomically resolved grain boundary line defects using the Nye tensor

The Nye tensor, which quantifies the density of Burgers vector for a given dislocation line direction, can be effectively used to characterize dislocation content in bulk crystals from atomic-resolution transmission electron microscopy images. The Nye tensor can be calculated from these images, in part because the reference state is simply defined by the lattice of the perfect crystal. The application of the Nye tensor to interfacial line defects, for which the natural reference state is the dichromatic pattern of the two grains in their reference orientation, poses additional challenges. In this work, we present a method that employs the Nye tensor to characterize the edge dislocation content of line defects at grain boundaries from atomic-resolution images. This approach enables us to rapidly characterize all edge dislocation content along a grain boundary. Additionally, the Nye tensor provides information about line defect core structure. Finally, we demonstrate this method on two exemplar defects: a twin boundary disconnection and a facet junction in face-centered cubic Au.

Crystallographic defects↗

DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗

DOE OSTI2025

Spin-strain interactions under hydrostatic pressure in α-RuCl 3

We investigate the effects of hydrostatic pressure on 𝛼−RuCl 3 , a prototypical material for the Kitaev spin model on a honeycomb lattice with a possible spin-liquid ground state. Using ultrasound measurements at pressures up to 1.16 GPa, we reveal significant modifications of the acoustic properties and the 𝐻−𝑇 phase diagram of this material. Hydrostatic pressure suppresses the three-dimensional magnetic order and induces a dimerization transition at higher pressures. At low pressures, the sound attenuation exhibits a linear temperature dependence, while above 0.28 GPa, it becomes nearly temperature independent, suggesting a shift in the phonon scattering regime dominated by Majorana fermions. These findings provide new insights into spin-strain interactions in Kitaev magnets and deliver a detailed characterization of the 𝐻−𝑇 phase diagram of 𝛼−RuCl 3 under hydrostatic pressure.

Focused ion beam↗

DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Mechanisms of Alkali Ionic Transport in Amorphous Oxyhalides Solid State Conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity (1 mS cm –1 ), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid‐state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ machine learning interatomic potentials to construct large‐scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal‐anion tetrahedra of various lengths. By analyzing the residence time of alkali cations migrating around tetrahedrally‐coordinated metals, we find that oxygen anions limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct‐particles correlations, while alkali transport is dictated by uncorrelated self‐diffusion. By extending this analysis to chemical compositions AMX 2.5 O 0.75 , spanning different alkaline (A = Li, Na, K), metallic (M = Al, Ga, In), and halogen (X = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic isovalent chemistry.

amorphous materials↗

DOE OSTI2023

Flexible and Extensive Platinum Ion Gel Condensers for Programmable Catalysis

Catalytic condensers composed of ion gels separating a metal electrode from a platinum-on-carbon active layer were fabricated and characterized to achieve more powerful, high surface area dynamic heterogeneous catalyst surfaces. Ion gels comprised of poly(vinylidene difluoride)/1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl) imide were spin coated as a 3.8 μm film on a Au surface, after which carbon sputtering of a 1.8 nm carbon film and electron-beam evaporation of 2 nm Pt clusters created an active surface exposed to reactant gases. Electronic characterization indicated that most charge condensed within the Pt nanoclusters upon application of a potential bias, with the condenser device achieving a capacitance of ∼20 μF/cm 2 at applied frequencies of up to 120 Hz. The maximum charge of ∼10 14 |e – | cm –2 was condensed under stable device conditions at 200 °C on catalytic films with ∼10 15 sites cm –2 . Grazing incidence infrared spectroscopy measured carbon monoxide adsorption isobars, indicating a change in the CO* binding energy of ∼19 kJ mol –1 over an applied potential bias of only 1.25 V. Condensers were also fabricated on flexible, large area Kapton substrates allowing stacked or tubular form factors that facilitate high volumetric active site densities, ultimately enabling a fast and powerful catalytic condenser that can be fabricated for programmable catalysis applications.

Catalysis↗

DOE OSTI2026

Phase Stability and Electrochemical Performance of La-Site-Doped Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 High-Entropy Garnets

We investigate La-site substitution in the high-entropy garnet Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 (LLZNTH) using Ba2+, Sr2+, and Sm3+ to elucidate how dopant governs phase stability, Li-site distribution, and electrochemical behavior. X-ray diffraction shows that Sr2+ is incorporated homogeneously into the garnet lattice, whereas the larger Ba2+ and smaller Sm3+ ions partially exceed the structural tolerance, generating secondary phases. Nevertheless, the Sm-doped composition (x = 0.05) exhibits the highest room-temperature ionic conductivity (2.7 × 10–4 S cm–1). Neutron powder diffraction reveals that Sm substitution drives a redistribution of Li+ from the tetrahedral 24 d sites into the higher-mobility 96 h positions, enhancing the connectivity of the three-dimensional Li-ion migration network. A Sm-doping series (x = 0.01–0.05) further shows that only sufficiently high Sm levels induce this redistribution, whereas lower concentrations retain Li arrangements similar to the undoped garnet. Critical current density measurements demonstrate that La-site dopants also influence interfacial stability against Li metal, underscoring a trade-off between bulk transport enhancement and mechanical robustness. Collectively, these findings reveal that in high-entropy garnets improved ionic conductivity can originate not only from phase-pure structures but also from targeted modification of the Li sublattice, even when accompanied by secondary phases, offering a compositional design principle for garnet electrolytes.

Li, Chang [Mechanical Engineering, School of Scien↗

DOE OSTI2025

Should I stay or should I flow? An exploration of phase‐separated metallosupramolecular liquid crystal polymers

Abstract Dynamic liquid crystalline polymers (dLCPs) incorporate both liquid crystalline mesogens and dynamic bonds into a single polymeric material. These dual functionalities impart order‐dependent thermo‐responsive mechano‐optical properties and enhanced reprocessability/programmability enabling their use as soft actuators, adaptive adhesives, and damping materials. While many previous works studying dynamic LCPs utilize dynamic covalent bonds, metallosupramolecular bonds provide a modular platform where a series of materials can be accessed from a single polymeric feedstock through the variation of the metal ion used. A series of dLCPs were prepared by the addition of metal salts to a telechelic 2,6‐bisbenzimidazolylpyridine (Bip) ligand endcapped LCP to form metallosupramolecular liquid crystal polymers (MSLCPs). The resulting MSLCPs were found to phase separate into hard and soft phases which aids in their mechanical robustness. Variations of the metal salts used to access these materials allowed for control of the thermomechanical, viscoelastic, and adhesive properties with relaxations that can be tailored independently of the mesogenic transition. This work demonstrates that by accessing phase separation through the incorporation of metallosupramolecular moieties, highly processable yet robust MSLCP materials can be realized. This class of materials opens the door to LCPs with bulk flow behavior that can also be utilized as multi‐level adhesives.

Chemistry↗

DOE OSTI2023

Fabrication of Large-Area Metal-on-Carbon Catalytic Condensers for Programmable Catalysis

Catalytic condensers stabilize charge on either side of a high-k dielectric film to modulate the electronic states of a catalytic layer for the electronic control of surface reactions. Here, carbon sputtering provided for fast, large-scale fabrication of metal–carbon catalytic condensers required for industrial application. Carbon films were sputtered on HfO 2 dielectric/p-type Si with different thicknesses (1, 3, 6, and 10 nm), and the enhancement of conductance and capacitance of carbon films was observed upon increasing the carbon thickness following thermal treatment at 400 °C. After Pt deposition on the carbon films, the Pt catalytic condenser exhibited a high capacitance of ∼210 nF/cm 2 that was maintained at a frequency ∼1000 Hz, satisfying the requirement for a dynamic catalyst to implement catalytic resonance. Temperature-programmed desorption of carbon monoxide yielded CO desorption peaks that shifted in temperature with the varying potential applied to the condenser (−6 or +6 V), indicating a shift in the binding energy of carbon monoxide on the Pt condenser surface. A substantial increase in capacitance (∼2000 nF/cm 2 ) of the Pt-on-carbon devices was observed at elevated temperatures of 400 °C that can modulate ∼10% of charge per metal atom when 10 V potential was applied. A large catalytic condenser of 42 cm 2 area Pt/C/HfO 2 /Si exhibited a high capacitance of 9393 nF with a low leakage current/capacitive current ratio (<0.1), demonstrating the practicality and versatility of the facile, large-scale fabrication method for metal–carbon catalytic condensers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗

DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Semiconducting Electrides Derived from Sodalite: A First-Principles Study

Electrides are ionic crystals, with electrons acting as anions occupying well-defined lattice sites. These exotic materials have attracted considerable attention in recent years for potential applications in catalysis, rechargeable batteries, and display technology. Among this class of materials, electride semiconductors can further expand the horizon of potential applications due to the presence of a band gap. However, there are only limited reports on semiconducting electrides, hindering the understanding of their physical and chemical properties. In recent work, we initiated an approach to derive potential electrides via selective removal of symmetric Wyckoff sites of anions from existing complex minerals. Herein, we present a follow-up effort to design semiconducting electrides from parental complex sodalites. Among four candidate compounds, we found that a cubic Ca 4 Al 6 O 12 structure with the I-43m space group symmetry exhibits perfect electron localization at the sodalite cages, with a narrow electronic band gap of 1.8 eV, making it suitable for use in photocatalysis. Analysis of the electronic structures reveals that a lower electronegativity of the surrounding cations drives greater electron localization and promotes the formation of an electride band near the Fermi level. Our work proposes an alternative approach for designing new semiconducting electrides under ambient conditions and offers guidelines for further experimental exploration.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Tautomerism induces bending and twisting of biogenic crystals

Understanding and exploiting material flexibility through phenomena such as the bending and twisting of molecular crystals has been a subject of increased interest owing to the number of applications that benefit from these properties, such as optoelectronics, mechanophotonics, soft robotics, and smart sensors. Here, we report the growth of spontaneously bent and twisted ammonium urate crystals induced by the keto–enol tautomerism of the urate molecule. The major tautomer is native to biogenic crystals, whereas the minor tautomer functions as an effective crystal growth modifier to induce naturally bent and twisted ammonium urate crystals. We show that the degree of curvature can be tailored based on the judicious selection of growth conditions. A combination of state-of-the-art microscopy and spectroscopy techniques are used to characterize the origin of bending. Spatially resolved nano-electron diffraction and high-resolution electron microscopy of naturally bent crystals show nearly single crystallinity with local lattice deformations generated by a combination of screw and edge dislocations. These observations are consistent with photoinduced force microscopy and contact resonance atomic force microscopy, which confirmed spatially resolved changes in the intermolecular interactions and the mechanical properties throughout the cross-sectional and axial regions of bent crystals. A mechanism of bending involving the generation of regionally specific dislocations is proposed as an alternative to more commonly reported models. These findings highlight a unique characteristic of tautomeric crystals that may have broader implications for other biogenic materials.

Science & Technology - Other Topics↗

DOE OSTI2026

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗