Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular functionalization”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2025

One-dimensional heterocyclic carbene–Au metal–organic frameworks bridging ultra-high vacuum models and scalable liquid-phase growth

The controlled design of molecule–metal interfaces is central to the development of functional nanomaterials for catalysis, sensing, and molecular electronics. Here we show that the adsorption of a Janus-type diimidazolium precursor on gold yields one-dimensional (1D) N-heterocyclic carbene (NHC)–Au–NHC metal organic frameworks (MOFs) featuring positively charged gold nodes. Using synchrotron X-ray photoemission spectroscopy (XPS), near edge X-ray adsorption fine structure (NEXAFS) spectroscopy and scanning tunnelling microscopy (STM), we demonstrate that thermal activation promotes counterion removal and drives the formation of extended 1D arrays, characterized by ∼1.0 nm Au–Au spacing and adatom densities up to 0.6 atom nm −2 (∼4% of surface atoms). Importantly, we translate this ultra-high vacuum (UHV) benchmark into a scalable solution-phase protocol in ethanol, enabling 1D-MOF growth under mild, base-free, open-air conditions. The resulting films retain structural and electronic signatures of UHV-grown systems, bridging model studies and practical synthesis. This approach establishes NHC–metal frameworks as accessible, tunable platforms for catalysis and materials design.

Gold adatoms↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Investigating the non-classical M-H 2 bonding in OsClH 3 (PPh 3 ) 3

Dihydrogen ligation to metal centers results in a wide range of metal-hydrogen interactions. Compounds containing non-classical dihydrogen ligands in which the H-H contacts are significantly longer than the 0.74 A distance found in solid H 2 are classified as Kubas-like compounds, stretched dihydrogen ligands, or compressed dihydrides, depending on the H-H contact length. Compounds containing stretched dihydrogen and compressed dihydride ligands are somewhat rare. Developing a more complete understanding of the metal-hydrogen and hydrogen-hydrogen interactions in these species may provide insights into the mechanisms for the oxidative addition reaction and hydrogen storage. Here, we use diffraction and inelastic neutron scattering (INS) measurements, paired with both molecular and solid-state density functional theory (DFT), to fully characterize the nature of the metal-hydrogen and hydrogen-hydrogen interactions in a previously reported compound, OsClH 3 (PPh 3 ) 3 . In that report, the H-H distance was measured as 1.48(2) A by single-crystal neutron diffraction and was described as a "stretched" dihydrogen. The INS-generated phonon density of states is well described by the DFT calculations, including those normal modes dominated by the vibrational motions of the three Os-ligated hydrogen atoms. Additionally, the resulting calculated electron charge density landscape indicates that there is no remnant bonding character between the Os-ligated hydrogen atoms in the compound. Based on these findings, we update the understanding of the metal-hydrogen bonding in this complex and reclassify it as a compressed dihydride. As such, OsClH 3 (PPh 3 ) 3 represents a rare example of a non-classical compressed dihydride species.

08 HYDROGEN↗
DOE OSTI2026

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗
DOE OSTI2025

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
DOE OSTI2025

Structure-guided utilization of lignocellulose for catalysis, energy, and biomaterials

As a complex composite of cellulose, hemicellulose, and lignin, plant lignocellulose has long served as a major resource for biomass conversion, materials engineering, and bio-based product development. High-resolution structural insights enabled by solid-state nuclear magnetic resonance (ssNMR) now allow the mapping of polymer interfaces, identification of functional group accessibility, and tracking of molecular organization during processing, all of which are critical factors for optimizing catalytic strategies. These insights could drive transformative progress in lignocellulose-based applications, including selective depolymerization, improved pretreatment design, and efficient upcycling of lignin into resins, plastics, and biomedical materials. In industry-relevant contexts, such as biofuel generation and renewable material manufacturing, understanding the hydration dynamics, cross-linking patterns, and structural heterogeneity is also essential. The ability to visualize these features in native biomass presents a unique opportunity to develop new strategies for sustainability and performance. As the structural toolbox continues to expand, it is becoming a central enabler for innovations in renewable energy, green chemistry, and advanced bioproducts.

bioproduct↗
DOE OSTI2024

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2020

Engineering Electrical Conductivity in Stable Zirconium-Based PCN-222 MOFs with Permanent Mesoporosity

Electrically conductive metal–organic frameworks (MOFs) featuring zirconium-based nodes are of great interest for electrochemical and optoelectronic applications due to their exceptional thermal and chemical stability, although the number of such MOFs remains limited. Here in this paper, we demonstrate that electron deficient molecules such as transition metal(IV) bis(dicarbollide) (TM IV (C 2 B 9 H 11 ) 2 ; TM IV CB for short, where TM = Ni, Pd, and Pt) or C 60 , spatially infiltrated in the microporous channels of the zirconium-based porphyrinic MOFs PCN-222 and PCN-222-Zn, generate highly stable and electrically conductive frameworks in which the mesopores remain accessible to other guests. Solid-state density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations indicate incorporation of Ni IV CB and C 60 in the MOF micropores is energetically and structurally feasible. Interestingly, in contrast to the free species, strong donor–acceptor interactions between the MOF and Ni IV CB restrain it to a single conformation. Calculated electronic structures and charge-hopping conduction probabilities illustrate that efficient charge-transfer (CT) from photoexcited linkers to the guest molecules facilitates charge hopping in the framework, making the MOFs electrically conductive. The donor–acceptor conjugates also enhance exciton dissociation at their heterojunctions, fostering the formation of long-lived electron-trapped states with potential utility for photo- and electrochemical devices.

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

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗
DOE OSTI2021

Influence of Cation Size on the Local Atomic Structure and Electronic Properties of Ta Perovskite Oxynitrides

Partial anion substitution in transition metal oxides provides rich opportunities to control and tune physical and chemical properties, for example, combining the merits of oxides and nitrides. In addition, the possibility of resulting anion sublattice order provides a means to target polar and chiral structures based on a wide array of accessible structural archetypes by design. Here, in this work, we investigate the local structures of a family of perovskite tantalum oxynitrides—ATaO 2 N (A = Ba, Sr, and Ca)—using a combination of experimental and theoretical approaches including neutron total scattering, density functional theory (DFT), and ab initio molecular dynamics (AIMD) simulations. We present the first experimental study of chemical short-range order (CSRO) in CaTaO 2 N, confirming local cis N ordering of the anion sub-lattice. Our systematic exploration of a local structure across the A cation size series (from the larger Ba to the smaller Ca) reveals a perovskite motif increasingly distorted with respect to long-range average structures. DFT and AIMD simulations support the observed trends. Overall, structures with cis ordering of the nitrogen anions in each TaO 4 N 2 octahedron are favored over those with trans ordering. With diminishing A cation size, local cis ordering and Ta off-centering play decreasing roles in overall lattice stability, overshadowed by the stabilizing effects of octahedral tilting. The influence of these factors on local dipole formation and frustrated dipole ordering are discussed.

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

In Situ Visualization of Local Distortions in the High- T c Molecule-Intercalated Li x (C 5 H 5 N) y Fe 2– z Se 2 Superconductor

A time-resolved synchrotron X-ray total scattering study sheds light on the evolution of the different structural length scales involved during the intercalation of the layered iron–selenide host by organic molecular donors, aiming at the formation of the expanded-lattice Li x (C 5 H 5 N) y Fe 2– z Se 2 hybrid superconductor. The intercalates are found to crystallize in the tetragonal ThCr 2 Si 2 -type structure at the average level, however, with an enhanced interlayer iron–selenide spacing ( d = 16.2 Å) that accommodates the heterocyclic molecular spacers. Quantitative atomic pair distribution function (PDF) analysis at variable times suggests distorted FeSe 4 tetrahedral local environments that appear swollen with respect to those in the parent β-FeSe. Simultaneously acquired in situ synchrotron X-ray powder diffraction data disclose that secondary phases (α-Fe and Li 2 Se) grow significantly when a higher lithium concentration is used in the solvothermal reaction or when the solution is aged. These observations are in line with the strongly reducing character of the intercalation medium’s solvated electrons that mediate the defect chemistry of the expanded-lattice superconductor. In the latter, intralayer correlated local distortions indicate electron-donating aspects that reflect in somewhat enlarged Fe–Se bonds. Additionally, they also reveal a degree of relief of chemical pressure associated with a large distance between Fe and Se sheets (“taller” anion height) and a stretched Fe–Fe square planar topology. The elongation of the latter, derived from the in situ PDF study, speaks for a plausible increase in the Fe-site vacancy concentration. The evolution of the local structural parameters suggests an optimum reaction window where kinetically stabilized phases resemble the distortions of the edge-sharing Fe–Se tetrahedra, required for a high- T c in expanded-lattice iron-chalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

NCAP: Noncanonical Amino Acid Parameterization Software for CHARMM Potentials

Noncanonical Amino Acids (NCAAs) provide numerous avenues for introduction of novel functionality to peptides and proteins. NCAAs can be incorporated through solid phase synthesis or genetic code expansion in conjugation with heterologous expression of the encoded protein modification. Due to the difficulty of synthesis, wide chemical space and lack of empirically resolved structures modeling the effects of NCAA mutation is critical for rational protein design. To evaluate the structural and functional perturbations NCAAs introduce we utilize molecular potentials that describe the forces in protein structure. Most potentials such as CHARMM are designed to model canonical residues but can be parameterized in include novel NCAAs. Here, in this work, we introduce NCAP a software package to generate CHARMM compatible parameters from quantum chemical calculation. Unlike currently available tools NCAP is designed to recognize NCAA structure and automatically bridge the gap between DFT calculations and potential parameters. For our software we discuss workflow, validation against canonical parameter sets and comparison to published NCAA-protein structures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

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

Quantum Chemical Simulations of CO 2 and N 2 Capture in Reline, a Prototypical Deep Eutectic Solvent

Deep eutectic solvents such as reline are an emerging class of low-cost, environmentally friendly solvents with tunable properties that are potentially applicable for the capture and separation of CO 2 . Experimental measurements showed that a reline-based membrane contactor can capture and separate CO 2 via physisorption through a dissolution process with 96.7% purity from a mixed gas containing CO 2 and N 2 (50:50% molar ratio). Here, we examine the nature of the interaction of CO 2 and N 2 with reline employing quantum chemical methods. We focus on explaining the mechanism by which CO 2 and N 2 bind to reline and the reason for the high selectivity for absorption of CO 2 compared to N 2 . We analyze the dynamics, energetics, and binding motifs for CO 2 and N 2 in reline employing density functional theory, density functional tight binding, and ab initio molecular dynamics. We also investigate the effect of reline on the vibrational spectra of CO 2 and reline. Our simulations indicate that the selective capture of CO 2 from the mixture of CO 2 and N 2 is due to the interplay between attractive electrostatic and charge polarization forces with opposing entropic effects, which shift the energetic balance and make the N 2 absorption unfavorable in reline.

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

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗