Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local coordinate transform”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fossil Fuel Transitions Framework: Case Studies of the Decision-Making Process for Energy and Economic Development Pathways

Each community is unique, and transitions must be tailored to local and national contexts, with recognition that access to financing is a key differentiator between wealthy and low-to-middle income countries. This analysis examines six community fossil fuel transitions in two different country contexts to provide insights on inclusive processes and decision-making criteria that bring forth positive outcomes for the communities. In both the United States and Chile, a coordination mechanism is established across government agencies and the importance of stakeholder engagement is recognized. In many cases, the government, local community groups, and electric utility or industry players have a critical role in transition planning and support.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Unravelling the mechanism of pressure induced polyamorphic transition in an inorganic molecular glass

The atomic structure of a germanium doped phosphorous selenide glass of composition Ge 2.8 P 57.7 Se 39.5 is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P 4 Se 3 molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe 4 tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P 3 ring at the base of the P 4 Se 3 molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

42 ENGINEERING↗
DOE OSTI2022

Strain-Induced asymmetry and on-site dynamics of silicon defects in graphene

In the last decade, the atomically-focused electron beams utilized in scanning transmission electron microscopes (STEMs) have been shown to induce a broad set of local structural transformations in materials, opening pathways for directing material synthesis and modification atom-by-atom. The mechanisms underlying these transformations remain largely unknown, due to the intractability of modeling the myriad of reaction pathways that can be accessed through high-energy electron scattering. The information on materials’ structure and dynamics that can be extracted from STEM images is similarly left underexplored. Here, we report the observation of anomalous on-site dynamics of individual silicon impurity atoms in graphene during STEM imaging. Density functional theory-based structural optimizations of anisotropically-strained molecular nanographenes reveal two distinct (but nearly degenerate) stable structures for four-fold coordinated silicon impurities, where interconversion between the two structures manifests slight changes of the silicon position within the lattice site. Implications for defect-based strain engineering in graphene are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the Feasibility of Market Manipulation and Energy Storage Arbitrage via Load-Altering Attacks

Around the globe, electric power networks are transforming into complex cyber–physical energy systems (CPES) due to the accelerating integration of both information and communication technologies (ICT) and distributed energy resources. While this integration improves power grid operations, the growing number of Internet-of-Things (IoT) controllers and high-wattage appliances being connected to the electric grid is creating new attack vectors, largely inherited from the IoT ecosystem, that could lead to disruptions and potentially energy market manipulation via coordinated load-altering attacks (LAAs). In this article, we explore the feasibility and effects of a realistic LAA targeted at IoT high-wattage loads connected at the distribution system level, designed to manipulate local energy markets and perform energy storage (ES) arbitrage. Realistic integrated transmission and distribution (T&D) systems are used to demonstrate the effects that LAAs have on locational marginal prices at the transmission level and in distribution systems adjacent to the targeted network.

25 ENERGY STORAGE↗
DOE OSTI2020

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF–Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born–Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. In this study, using tetrahydrofuran (THF)–water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D 8 THF to examine the molecular origins of the isotope-dependent loss of the THF–water miscibility gap.

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

Structural features of solid-solid phase transitions and lattice dynamics in U 3 O 8

Triuranium octoxide (U 3 O 8 ) undergoes an orthorhombic to hexagonal structural phase transition near Ts=305°C, and a separate nonstructural phase transition at Tc=210°C. The later transition has previously been associated with temperature-induced fluctuations in the uranium oxidation state. A discontinuity in the slope of electrical conductivity versus temperature measurement at 210°C has supported this idea. The orthorhombic phase has three crystallographic sites in two distinct oxidation configurations [2 U(V) and 1 U(VI)], whereas the hexagonal phase has one distinct uranium site. High-resolution x-ray diffraction measurements eliminate the possibility of superlattice Bragg reflections to less than 0.2 e - scattering power and U 3 O 8 is not metallic; consequently, the presence of oxidation fluctuations is required for charge balancing. Interestingly, the order-to-disorder transition occurs at a much lower temperature than the structural transition. Using temperature-dependent x-ray diffraction and Raman spectroscopy, we show anisotropic lattice expansion in the in-plane b and c lattice constants. A specific discontinuity in the temperature derivatives of the b and c lattice constants are the first reported structural signatures of the order-to-disorder transition, suggestive of a change in local U–O coordination. Phonon frequencies of U 3 O 8 measured by Raman spectroscopy show significant temperature-dependent dynamics. Redshifting of several modes between 40 and 300°C cannot be explained by unit cell expansion alone because the unit cell volume decreases in this region. Instead, we show that phonon frequencies are highly correlated with the anisotropic lattice expansion/contraction along specific crystallographic directions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Envisioning U.S. Climate Predictions and Projections to Meet New Challenges

In the face of a changing climate, the understanding, predictions, and projections of natural and human systems are increasingly crucial to prepare and cope with extremes and cascading hazards, determine unexpected feedbacks and potential tipping points, inform long-term adaptation strategies, and guide mitigation approaches. Increasingly complex socio-economic systems require enhanced predictive information to support advanced practices. Such new predictive challenges drive the need to fully capitalize on ambitious scientific and technological opportunities. These include the unrealized potential for very high-resolution modeling of global-to-local Earth system processes across timescales, reduction of model biases, enhanced integration of human systems and the Earth Systems, better quantification of predictability and uncertainties; expedited science-to-service pathways, and co-production of actionable information with stakeholders. Enabling technological opportunities include exascale computing, advanced data storage, novel observations and powerful data analytics, including artificial intelligence and machine learning. Looking to generate community discussions on how to accelerate progress on U.S. climate predictions and projections, representatives of Federally-funded U.S. modeling groups outline here perspectives on a six-pillar national approach grounded in climate science that builds on the strengths of the U.S. modeling community and agency goals. This calls for an unprecedented level of coordination to capitalize on transformative opportunities, augmenting and complementing current modeling center capabilities and plans to support agency missions. Tangible outcomes include projections with horizontal spatial resolutions finer than 10 km, representing extremes and associated risks in greater detail, reduced model errors, better predictability estimates, and more customized projections to support next generation climate services.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Ions Tune Interfacial Water Structure and Modulate Hydrophobic Interactions at Silica Surfaces

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy–Chapman–Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface–solvent–ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica–aqueous interface by shifting the structural makeup in the Stern layer from dominant water–surface interactions to water–water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air–water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T 1 ) of excited interfacial O–H stretching modes from T 1 ~ 600 fs to T 1 ~ 250 fs. Hence, we propose that the T 1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

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

A classical molecular dynamics simulation method for the formation of “dry” gels from boro-aluminosilicate glass structures

In contact with water, glass transforms into amorphous and porous structures called gels. A simulation method based on classical molecular dynamics is proposed here to mimic “dry” gels forming from initial oxide glass structures. Six glass compositions were investigated. Two behaviours were evidenced depending on the initial glass composition, and in particular on the quantity of elements removed. If a large quantity of soluble elements (B, Na) was removed, it induced an increase in the average pore size within the gels, and the time needed to stabilise the gel structure increased because more local atomic rearrangements occurred. The gel network displayed a higher proportion of Si-Q4 at the expense of Si-Q3 and a lower average ring size compared to the glass network, irrespective of the glass composition. Surface effects were also highlighted in the dry gels, such as the presence of 3-coordinated Al and a decrease in the average angle Si-O-Si and Al-O-Al. Our findings will be compared to both wet gels and experimental data in further studies, to help find the best procedure to simulate such structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pro‐ L * ‐ A Probabilistic L * Mapping Tool for Ground Observations

Abstract Both ground and space observations are used extensively in the modeling of space weather processes within the Earth’s magnetosphere. In radiation belt physics modeling, one of the key phase‐space coordinates is L *, which indicates the location of the drift paths of energetic electrons. Global magnetic field models allow a subset of locations on the ground (mainly subauroral) to be mapped along field lines to a location in space and transformed into L *, provided that the initial ground location maps to a closed drift path. This allows observations from ground, or low‐altitude space‐based platforms to be mapped into space in order to inform radiation belt modeling. Many data‐based magnetic field models exist; however, these models can significantly disagree on mapped L * values for a single point on the ground, during both quiet times and storms. We present a state of the art probabilistic L * mapping tool, Pro‐ L *, which produces probability distributions for L * corresponding to a given ground location. Pro‐ L * has been calculated for a high resolution magnetic latitude by magnetic local time grid in the Earth’s Northern Hemisphere. We have developed the probabilistic model using 11 years of L * calculations for seven widely used magnetic field models. Usage of the tool is highlighted for both event studies and statistical models, and we demonstrate a number of potential applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Luminescent N-heterocyclic carbene Cu( I ) complexes with N^O chelating ligands exhibit microsecond lifetimes and photocatalytic activity

To replace precious noble metal-based photosensitizers in applications involving photoinduced charge separation, energy transfer, or photocatalysis, Cu(I) complexes are considered to be cost-effective, earth-abundant, and sustainable alternatives. An emerging and effective design principle in Cu(I) photosensitizers involves heteroleptic structures where the HOMO and LUMO are spatially separated over two different ligands. In the present work, we introduce a complementary class of heteroleptic, three-coordinate copper photosensitizers that pairs variable N^O chelating ligands (8-hydroxyquinoline and 10-hydroxybenzo[h]quinoline) with a bulky 2,6-diisopropylphenyl-substituted N-heterocyclic carbene (NHC). In this design, both frontier orbitals are localized on the same ligand, the N^O chelate, such that structural modulation of the electron-rich N^O-chelates can substantially tune the energy levels of the HOMO and LUMO, thereby controlling the photoluminescence properties. Detailed photophysical and electrochemical experiments as well as DFT calculations suggest charge-transfer transitions with intra-ligand charge transfer (ILCT) character, involving the N^O ligands. This strategy successfully produced long triplet excited-state lifetimes (up to 44 µs) in compounds that are strong photoreductants (E([Cu] + /*[Cu] as negative as −2.0 V vs. the ferrocenium/ferrocene couple). These properties allow these photosensitizers to be used as photocatalysts in various transformations of organic compounds, such as hydrogenation of substituted benzophenones, hydrodehalogenation of aryl/alkyl halides (including challenging C–Cl bond activation) and E/Z isomerization of (E)-stilbene (an example of triplet–triplet energy transfer).

Chakraborty, Soumi [University of Houston, TX (Uni↗
DOE OSTI2021

Individual and Collective Strategies to Limit the Impacts of Large Power Outages of Long Duration

As modern society becomes ever more dependent on the availability of electric power, the costs that could arise from individual and social vulnerability to large outages of long duration (LLD-outages) increases. During such an outage, even a small amount of power would be very valuable. Here, this article compares individual and collective strategies for providing limited amounts of electric power to residential customers in a hypothetical New England community during a large electric power outage of long duration. We develop estimates of the emergency load required for survival and assess the cost of strategies to address outages that last 5, 10, and 20 days in either winter or summer. We find that the cost of collective solutions could be as much as 10 to 40 times less than individual solutions (less than $2 per month per home). However, collective solutions would require community-wide coordination, and if local distribution system lines are destroyed, only individual back-up systems could provide contingency power until those lines are repaired. Costs might be reduced if more robust distributed generation were employed that could be operated continuously with the ability to sell power back to the grid. Our cost-effectiveness analysis only assesses what could be done, developing estimates of preparedness cost. A decision about what should be done would require additional input from a range of stakeholders as well as some form of analytical deliberative process.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Evidence of a liquid–liquid transition in a glass-forming ionic liquid

A liquid–liquid transition (LLT) is a transformation from one liquid to another through a first-order transition. The LLT is fundamental to the understanding of the liquid state and has been reported in a few materials such as silicon, phosphorus, triphenyl phosphite, and water. Furthermore, it has been suggested that the unique properties of materials such as water, which is critical for life on the planet, are linked to the existence of the LLT. However, the experimental evidence for the existence of an LLT in many molecular liquids remains controversial, due to the prevalence and high propensity of the materials to crystallize. In this work, we show evidence of an LLT in a glass-forming trihexyltetradecylphosphonium borohydride ionic liquid that shows no tendency to crystallize under normal laboratory conditions. We observe a step-like increase in the static dielectric permittivity at the transition. Furthermore, the sizes of nonpolar local domains and ion-coordination numbers deduced from wide-angle X-ray scattering also change abruptly at the LLT. We independently corroborate these changes in local organization using Raman spectroscopy. The experimental access to the evolution of local order and structural dynamics across a liquid–liquid transition opens up unprecedented possibilities to understand the nature of the liquid state.

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

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗
DOE OSTI2022

Effect of surface acidity modulation on Pt/Al 2 O 3 single atom catalyst for carbon monoxide oxidation and methanol decomposition

With ideal metal utilization efficiency and homogeneous active sites, single atom catalysts (SACs) have attracted extensive attention in heterogeneous catalysis field. The performance of SACs is highly dependent on the local coordination environment of metal single atoms. However, few works reported the effect of surface acidity modification on the performance of SACs such as Pt SACs which have been widely studied nowadays. In this work, a Pt/Al 2 O 3 SAC has been successfully prepared using a bayerite Al(OH) 3 as support through simple IWI method for Pt loading. Through post-modification of Pt/Al 2 O 3 with acidic WO 3 and basic MgO, the surface acidity of Pt/Al 2 O 3 was tuned. Here, activity of CO oxidation and methanol decomposition on the modified catalysts has been evaluated. It was found that the catalytic performance was highly related to the surface acidity of Pt SACs. Basic MgO modified Pt/Al 2 O 3 performed better in CO oxidation under the reaction condition with H 2 O. Acidic WO 3 modified Pt/Al 2 O 3 showed superior activity in methanol decomposition. Detailed characterization and in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) study suggested that the positive H 2 O effect on CO oxidation was due to the mitigated formation of bicarbonate species on Pt/Al 2 O 3 and WO 3 /Pt/Al 2 O 3 catalysts and the promoted desorption of carbonate species from MgO/Pt/Al 2 O 3 . The superior methanol decomposition performance on WO 3 /Pt/Al 2 O 3 was owing to its abundant oxygen vacancies and high resistance to carbonate accumulation. This study highlights the importance of surface acidity modulation of Pt SACs in heterogeneous catalysis and suggests that fine-tuning the surface acidity can be an effective strategy in new SACs development.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2026

Interpretable Deep Learning for Advancing Field-Enhanced Catalysis

This DOE Early Career project developed a physics-informed, interpretable AI-and-modeling framework to understand and exploit electric-field effects in heterogeneous catalysis, with ammonia cracking and synthesis as a representative pathway. The team built and validated methods to map local electric fields on metal surfaces and nanoparticles, showing that low-coordination features (tips/edges/corners) can concentrate fields by several-fold relative to flat facets. Using DFT-generated datasets, the project created physics-guided machine learning models that rapidly predict local electric fields and field-dependent adsorption energetics with near-DFT accuracy while reducing computational cost by orders of magnitude. These predictions were integrated with microkinetic modeling to quantify how field-dipole interactions reshape reaction energetics and mechanisms, enabling large increases in predicted catalytic rates and substantial reductions in operating temperature under favorable field conditions. To accelerate discovery of earth-abundant catalysts, the project combined interpretable ML screening (with electronic-structure descriptors identified as key drivers) with a generative inverse-design workflow based on diffusion models and physics constraints. The resulting closed-loop approach, linking simulation, mechanistic modeling, and AI, provides reusable tools and datasets for designing catalysts and operating conditions in field-enhanced catalysis, with broad relevance to electrostatic catalysis, plasma catalysis, electrocatalysis, and other energy-related chemical transformations.

30 DIRECT ENERGY CONVERSION↗