Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CaSe”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Intervening alcohol marketing to reduce harmful alcohol use and lessons learned from the theory of changes: Case studies in Thailand

Objectives Globally, the burden of disease caused by alcohol use has been steadily increasing, including in Thailand. In this study, we aim to test the effectiveness of Anderson et al.‘s suggested three approaches to change the collective social norms, which comprise of: (1) providing information and an understanding about alcohol use behaviour, its causes and distribution; (2) focusing strategies on groups rather than individuals; and (3) strengthening supportive laws, regulations and approaches. Study design We employed a mixed-methods approach. Evidence was gathered from literature review and in-depth interviews with key individuals who are responsible for community-based interventions to alcohol marketing strategies in Thailand. Methods We chose to focus on two case studies in Nan and Surin provinces, where hospital-based longitudinal data (8 years) were available. Changes in casualties related to the harmful use of alcohol, resulting from interactions between community-based interventions and alcohol marketing during the time of annual festivals were investigated. We employed the theory of change (ToC) defined by Vogel to guide the data collection and analysis. We reviewed literature from online databases and grey literature to generate causal-loop diagrams. Results We created a causal-loop diagram to describe the complexity of harmful alcohol use, its related factors, context, interventions and outcomes. Over the decade between 2006 and 2015, community-based strategies led to a substantial reduction of casualties (initially a 50% reduction, rising to an 80–90% reduction by the end of the study period) during the time of the festivals. Conclusions The reduction in injuries and fatalities could be a result of the concerted actions, including legal sanctions of alcohol beverage sales and advertisement, and public education to raise awareness and impart knowledge of the harmful use of alcohol. The actions were organised by a coalition of civil society, health professionals, public authorities and community leaders using hospital-based data on the adverse effects of harmful alcohol use to mobilise political support at the provincial level. The availability of long-term financial support as a catalytic source of funds and the presence of a comprehensive alcohol control act enabled framing and mobilisation of local resources and political support.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The Circular Economy Life Cycle Assessment and Visualization Framework: A Multistate Case Study of Wind Blade Circularity in United States

A circular economy (CE) aims to decouple human activities from economic activities and resource use, and its overall goal is reducing or avoiding negative environmental externalities. The newly developed Circular Economy Lifecycle Assessment and Visualization (CELAVI) framework simulates changes in supply chain environmental impacts as it transitions toward circularity. This study expands CELAVI by incorporating detailed spatial resolution and real-world road routes connecting all facilities within the system. The case study on end-of-life decision making of wind turbine blades in the states of Iowa and Missouri explores how supply chain circularity and environmental impacts are affected by pathway costs and level of wind turbine installations. It demonstrates how high circularity costs might be beneficial for circularity transitions given revenue generated from circular pathways. Finally, impacts have important contributions to the supply chain design and thus show the importance of including detailed spatial resolution in CELAVI and CE studies in general.

circular economy↗
DOE OSTI2025

Bayesian Optimized Deep Ensemble for Uncertainty Quantification of Deep Neural Networks: a System Safety Case Study on Sodium Fast Reactor Thermal Stratification Modeling

Deep neural networks (DNNs) are increasingly important to scientific computing and engineering system simulations. Accurate uncertainty quantification (UQ) for DNNs is critical in safety-sensitive engineering domains. Traditional Deep Ensemble (DE) methods, while easy to implement, frequently suffer from poorly calibrated uncertainty estimates and limited predictive accuracy due to reliance on fixed architectures with varied weight initializations. To address these issues, we introduce a workflow that combines Bayesian Optimization (BO) and DE. The workflow is modular, scalable, and integrates parallel BO initialized with Sobol sequences to individually optimize the hyperparameters of each ensemble member. This method enhances ensemble diversity, improves predictive accuracy, and provides reliable uncertainty estimates. We evaluate the proposed BODE approach in a sodium fast reactor thermal stratification modeling case study, where we used a densely connected convolutional neural network to predict turbulent viscosity during the reactor transient with consideration of data noise. We benchmark its performance against several optimization approaches, including baseline deep ensemble, evolutionary algorithm-optimized ensemble, ensemble formed via random search combined with greedy selection, and a BO ensemble using random initialization. Here, our results demonstrate superior performance of the developed BODE approach. In noise-free scenarios, BODE notably reduces incorrect aleatoric uncertainty and significantly enhances predictive accuracy. Under conditions of 5% and 10% Gaussian noise, BODE adaptively quantifies uncertainty proportional to data noise, achieving up to an 80% reduction in root mean square error compared to baseline methods and producing well-calibrated prediction intervals.

Bayesian optimization↗
DOE OSTI2022

Carbon footprint prediction considering the evolution of alternative fuels and cargo: A case study of Yangtze river ships

For reducing carbon emissions in the shipping sector, application of alternative low-carbon and zero-carbon fuels is the consensus. However, requirements of economic development cannot be ignored while focusing on emission reduction. Cargo development is an important prerequisite. Considering the cargo growth as a novelty, a method to calculate the annual carbon emissions of different alternative fuel-powered ships in different cargo growth prospects and power scenarios is proposed. With respect to the time dimension, the relationship between alternative fuels, cargo, and carbon reduction is revealed. The Yangtze River bulk carrier case was studied and the life cycle carbon emission analysis of diesel, LNG hybrid, LNG, hydrogen, methanol, and ammonia were carried out. The annual carbon curves of the high, steady, and low cargo growth prospects were obtained for different scenarios of the power system structure. The results show that LNG hybrid, LNG and methanol fuels are currently the suitable choices. Their life cycle carbon emissions were reduced by 31.5–38.1% compared with those of diesel power. The carbon emissions of green hydrogen and green ammonia were reduced by 78.8% and 91.3%, respectively, compared with those of hydrogen and ammonia. The carbon reduction results of green fuels in their initial application were explored which are at 0.6–10.8%. In addition, there is a balanced annual growth rate, whose annual carbon curve will be steady under its corresponding scenario of power system structure, which could help the ship owners plan their future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Defining weather scenarios for simulation-based assessment of thermal resilience of buildings under current and future climates: A case study in Brazil

In response to increasingly severe weather conditions, optimization of building performance and investment provides an opportunity to consider co-benefits of thermal resilience during energy efficiency retrofits. This work aims to assess thermal resilience of buildings using building performance simulation to evaluate the indoor overheating risk under nine weather scenarios, considering historical (2010s), mid-term future (2050s), and long-term future (2090s) typical meteorological years, and heat wave years. Such an analysis is based on resilience profiles that combine six integrated indicators. A case study with a district of 92 buildings in Brazil was conducted, and a combination of strategies to improve thermal resilience was identified. Results reflect the necessity of planning for resilience in the context of climate change. This is because strategies recommended under current conditions might not be ideal in the future. Therefore, an adaptable design should be prioritized. Cooling energy consumption could increase by 48 % by the 2050s, while excessive overheating issues could reach 37 % of the buildings. Simple passive strategies can significantly reduce the heat stress. A comprehensive thermal resilience analysis should ultimately be accompanied by a thorough reflection on the stakeholders’ objectives, available resources, and planning horizon, as well as the risks assumed for not being resilient.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A machine learning method of modern urban building energy modeling: A case study of Chicago

Urban-scale building energy modeling is vital for urban planning. However, it can be challenging to assimilate reliable non-geometry building data for urban-scale modeling without extensive investment. Here, this study introduces a novel approach to developing modern urban-scale building energy stock data using geographic information systems and machine learning algorithms without necessarily requiring pre-supplied non-geometric metadata. The proposed framework integrates building footprint and height data to estimate gross floor areas, and matches each building to a pool of candidate records from ComStock or ResStock—filtered to the same county and ranked by geometric similarity—demonstrate a proof-of-concept case study in Chicago for predicting energy use intensity (EUI) using scalable datasets. The model achieved a mean bias error (MBE) of 0.08 kWh/m² and root mean square error (RMSE) of 14.84 kWh/m² under full metadata input for EUI prediction. With only location inputs, the model captured 69.2 % of EUI within predicted ranges. These results demonstrate the model’s potential to support early-stage urban planning, identify candidates for energy-efficient retrofits. By removing the dependency on detailed pre-surveys or extensive building metadata, the approach overcomes a key barrier in traditional urban-scale building energy modeling, illustrating a pathway toward broader and more cost-effective application, though further multi-city validation and improved treatment of pre-1925 buildings are needed.

Energy Use Intensity↗
DOE OSTI2023

The inclusion of uncertainty in circularity transition modeling: A case study on wind turbine blade end-of-life management

The transition to a more circular economy (CE) is complex and hard to predict, including in sustainable energy technologies. Many sources of uncertainty make it challenging to model CE scenarios and their potential benefits. As an example, the high variability in costs and revenues of different recycling options makes future wind turbine blade recycling highly uncertain. To better understand this challenge, the circular economy, life cycle assessment and visualization (CELAVI) framework - a discrete event simulation and life cycle assessment framework - is modified to incorporate uncertainty analysis capabilities. Moreover, a 3-step procedure that covers different aspects of uncertainty in CE studies and includes a Monte-Carlo analysis is proposed. The procedure is tested in a case study on wind turbine blade recycling using CELAVI. Results highlight that grinding and landfilling costs are the most influential parameters for wind power circularity. The model's output coefficients of variation (when input parameter uncertainties are propagated) are between 92% and 384% depending on the indicator. The approach developed in this study may help researchers and decision-makers who study circularity prioritize their data collection effort. Finally, our method contributes to a mounting yet critical body of research: the measurement of uncertainty in circularity transitions.

17 WIND ENERGY↗
DOE OSTI2024

A country-scale green energy-water-hydrogen nexus: Jordan as a case study

Many developing countries suffer from a shortage of clean energy-water production, causing significant dependency on imported fossil fuel to meet the local energy and freshwater demands. Furthermore, when planning a 100 % renewable energy system (RES) to match the energy demand on an hourly basis, the RES is usually oversized. Hence, unavoidable large amounts of excess energy would be generated during hours of high production, low demand, and fully charged energy storage systems (ESSs). Hence, using Jordan as a case study, this work proposes a novel integrated system of wind, solar photovoltaic (PV), and lithium-ion ESS to match 100 % of the country's energy demand while using the excess generated power to drive reverse osmosis water desalination plants to match the demand of freshwater as well. Furthermore, the remaining excess energy is used for producing green hydrogen. A techno-economic optimizer was developed to show that this concept is indeed feasible. Here, the model was used to scan Jordan for the fitting sites of the highest demand-supply matching and the lowest levelized cost of electricity (LCOE), simultaneously. Three different integration scenarios are presented, namely PV-ESS, wind-ESS, and hybrid PV-wind-ESS RESs. The hybrid system showed the best performance, especially compared to the wind-ESS RES, in terms of total installed capacity (33.37 GW), LCOE (0.0492 USD/kWh), specific water cost (0.3629 USD/m 3 ), and energy and water demand-supply fractions (99.47 and 96.16 %, respectively). Storing desalinated water in natural aquifers was found to be a favorable option to enhance the performance of the three systems, where stored water by the end of the year was sufficient to cover the freshwater requirements for several following years. Finally, the remaining excess energy is used to produce green hydrogen with an annual production of up to 1.37 million tons and a specific cost down to 1.08 USD/kg.

14 SOLAR ENERGY↗
DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

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

An evaluation of electricity system technology competitiveness metrics: The case for profitability

In this work we develop standardized functional forms for electricity technology competitiveness metrics and show how the mathematical relationship between value and cost can affect the robustness of the metric. Here, we make the case to replace established metrics with economic profitability metrics – including a new profitability-based adjustment to LCOE. These profitability metrics are shown to be more robust for evaluating technology competitiveness by keeping comparisons on an equivalent monetary basis.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Arctic shipping under global change: A case study of offshore oil exports

We explore impacts of sea ice thinning and evolutions in the energy sector on future use of the Northern Sea Route (NSR) versus the Suez Canal Route (SCR), using a case study of shipping oil extracted from the offshore Russian Arctic to China. We combine an integrated human-Earth system model with a shipping cost model to incorporate impacts on both oil production and shipping costs under internally consistent scenarios. We find that the NSR could become cost-competitive with the SCR as sea ice thickness declines, especially in an RCP8.5 scenario, due to decreasing fuel and icebreaker escort costs. In a global energy evolution scenario consistent with RCP2.6, high emissions costs on the longer SCR may outweigh the costs associated with thicker sea ice on the NSR. Our novel framework provides integrated projections of NSR shipping traffic driven by a specific commodity likely to be shipped through the Arctic.

Arctic↗
DOE OSTI2022

Life cycle assessment of utilizing freshly cut urban wood: A case study

Large amounts of urban wood resources are generated through tree removals in urban areas. Therefore, there is a growing interest to improve the environmental performance of the urban wood supply chain by quantifying the environmental impacts of creating high-value products from urban tree removals to enhance its utilization. By surveying existing urban wood utilization operations for primary data in two major cities: Baltimore, Maryland, and Milwaukee, Wisconsin, a new life cycle inventory (LCI) dataset per m 3 of dry urban hardwood lumber and live edge slabs was created. Incorporating the new LCI data into life cycle assessment method, the total global warming (GW) impact of converting urban trees to kiln-dried lumber was estimated to be 122 kg CO 2 eq and 336 kg CO 2 eq per m 3 of lumber produced for Baltimore and Milwaukee, respectively. In both cases, the total environmental impact of the product was dominated by the conversion of urban trees to kiln dried lumber at the sawmill processing stage. Using scenario analysis for the source of heat generation used in kiln (force) drying wood, substituting wood fuel from coproducts for natural gas in a boiler for Milwaukee resulted in a substantial reduction in fossil-derived GW impact. Urban and traditional hardwood lumber production follow similar GW impact trajectories over their life-cycle stages with the lumber processing stage having the highest environmental impacts for both. Finally, increasing the overall lumber yield and decreasing kiln (force) drying would substantially improve the environmental performance of urban hardwood lumber and help make it comparable to traditional lumber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of (NO)Fe(N 2 S 2 ) as a Metallodithiolate Spin Probe Ligand: A Case Study Approach

The ubiquity of sulfur–metal connections in nature inspires the design of bi- and multimetallic systems in synthetic inorganic chemistry. Common motifs for biocatalysts developed in evolutionary biology include the placement of metals in close proximity with flexible sulfur bridges as well as the presence of π-acidic/delocalizing ligands. This Account will delve into the development of a (NO)Fe(N 2 S 2 ) metallodithiolate ligand that harnesses these principles. The Fe(NO) unit is the centroid of a N 2 S 2 donor field, which as a whole is capable of serving as a redox-active, bidentate S-donor ligand. Its paramagnetism as well as the ν(NO) vibrational monitor can be exploited in the development of new classes of heterobimetallic complexes. We offer four examples in which the unpaired electron on the {Fe(NO)} 7 unit is spin-paired with adjacent paramagnets in proximal and distal positions. First, the exceptional stability of the (NO)Fe(N 2 S 2 )-Fe(NO) 2 platform, which permits its isolation and structural characterization at three distinct redox levels, is linked to the charge delocalization occurring on both the Fe(NO) and the Fe(NO) 2 supports. This accommodates the formation of a rare nonheme {Fe(NO)} 8 triplet state, with a linear configuration. A subsequent FeNi complex, featuring redox-active ligands on both metals (NO on iron and dithiolene on nickel), displayed unexpected physical properties. Our research showed good reversibility in two redox processes, allowing isolation in reduced and oxidized forms. Various spectroscopic and crystallographic analyses confirmed these states, and Mössbauer data supported the redox change at the iron site upon reduction. Oxidation of the complex produced a dimeric dication, revealing an intriguing magnetic behavior. The monomer appears as a spin-coupled diradical between {Fe(NO)} 7 and the nickel dithiolene monoradical, while dimerization couples the latter radical units via a Ni 2 S 2 rhomb. Magnetic data (SQUID) on the dimer dication found a singlet ground state with a thermally accessible triplet state that is responsible for magnetism. A theoretical model built on an H 4 chain explains this unexpected ferromagnetic low-energy triplet state arising from the antiferromagnetic coupling of a four-radical molecular conglomerate. For comparison, two (NO)Fe(N 2 S 2 ) were connected through diamagnetic group 10 cations producing diradical trimetallic complexes. Antiferromagnetic coupling is observed between {Fe(NO)} 7 units, with exchange coupling constants (J) of -3, -23, and -124 cm –1 for Ni II , Pd II , and Pt II , respectively. This trend is explained by the enhanced covalency and polarizability of sulfur-dense metallodithiolate ligands. A central paramagnetic trans-Cr(NO)(MeCN) receiver unit core results in a cissoid structural topology, influenced by the stereoactivity of the lone pair(s) on the sulfur donors. This {Cr(NO)} 5 radical bridge, unlike all previous cases, finds the coupling between the distal Fe(NO) radicals to be ferromagnetic (J = 24 cm –1 ). The stability and predictability of this S = 1/2 moiety and the steric/electronic properties of the bridging thiolate sulfurs suggest it to be a likely candidate for the development of novel molecular (magnetic) compounds and possibly materials. The role of synthetic inorganic chemistry in designing synthons that permit connections of the (NO)Fe(N 2 S 2 ) metalloligand is highlighted as well as the properties of the heterobi- and polymetallic complexes derived therefrom.

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

SARS-CoV-2 RNA in Wastewater Settled Solids Is Associated with COVID-19 Cases in a Large Urban Sewershed

Wastewater-based epidemiology may be useful for informing public health response to viral diseases like COVID-19 caused by SARS-CoV-2. Here, we quantified SARS-CoV-2 RNA in wastewater influent and primary settled solids in two wastewater treatment plants to inform the preanalytical and analytical approaches and to assess whether influent or solids harbored more viral targets. The primary settled solids samples resulted in higher SARS-CoV-2 detection frequencies than the corresponding influent samples. Likewise, SARS-CoV-2 RNA was more readily detected in solids using one-step digital droplet (dd)RT-PCR than with two-step RT-QPCR and two-step ddRT-PCR, likely owing to reduced inhibition with the one-step ddRT-PCR assay. We subsequently analyzed a longitudinal time series of 89 settled solids samples from a single plant for SARS-CoV-2 RNA as well as coronavirus recovery (bovine coronavirus) and fecal strength (pepper mild mottle virus) controls. SARS-CoV-2 RNA targets N1 and N2 concentrations correlated positively and significantly with COVID-19 clinically confirmed case counts in the sewershed. Together, the results demonstrate that measuring SARS-CoV-2 RNA concentrations in settled solids may be a more sensitive approach than measuring SARS-CoV-2 in influent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

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

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

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

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na 2 + ) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that for the Na 2 + molecule, solute--solvent interactions lead to two stable, chemically-distinct coordination states (Na(THF) 4 -Na(THF) 5 + and Na(THF) 5 -Na(THF) 5 + ) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Lastly, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

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

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗