Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular processes”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Random field optimization

Herein we present a new modeling paradigm for optimization that we call random field optimization. Random fields are a powerful modeling abstraction that aims to capture the behavior of random variables that live on infinite-dimensional spaces (e.g., space and time) such as stochastic processes (e.g., time series, Gaussian processes, and Markov processes), random matrices, and random spatial fields. This paradigm involves sophisticated mathematical objects (e.g., stochastic differential equations and space-time kernel functions) and has been widely used in neuroscience, geoscience, physics, civil engineering, and computer graphics. Despite of this, however, random fields have seen limited use in optimization; specifically, existing optimization paradigms that involve uncertainty (e.g., stochastic programming and robust optimization) mostly focus on the use of finite random variables. This trend is rapidly changing with the advent of statistical optimization (e.g., Bayesian optimization) and multi-scale optimization (e.g., integration of molecular sciences and process engineering). Our work extends a recently-proposed abstraction for infinite-dimensional optimization problems by capturing more general uncertainty representations. Moreover, we discuss solution paradigms for this new class of problems based on finite transformations and sampling, and identify open questions and challenges.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Gas Phase Synthesis of Benzene via the Propargyl Radical Self-Reaction

Polycyclic Aromatic Hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs - benzene - via the self-recombination of two resonantly stabilized propargyl (C 3 H 3 ) radicals in dilute environments using isomer selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, 2-ethynyl-1,3-butadiene) and o-benzyne are detected and their branching ratios quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. Furthermore, these results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments, but also in combustion systems gaining us a better understanding of the hydrocarbon chemistry of our universe.

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

Spectroscopic Investigations of Molecular Symmetry Breakdown

This project addressed precise measurements of the spectra of molecules to address fundamental questions in chemical physics in the realm of nuclear spin symmetries and their breakdown in molecules, and parity non-conservation caused by the weak nuclear force. Nuclear spin symmetry is intimately connected to parity in symmetrical molecules and plays a central role in the physical and chemical properties of symmetric molecules, their collisional relaxation, and the establishment of thermodynamic equilibrium in an ensemble. Such effects are present in every molecular system or process, but their small size means they are often only directly observable under special conditions. However, their measurement is of fundamental, and often practical importance, and strongly relates to the themes of coherence in light and matter in the recent Basic Energy Sciences Advisory Committee (BESAC) report “Challenges at the Frontiers of Matter and Energy: Transformative Opportunities for Discovery Science”.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

A single-cell atlas of the bobtail squid visual and nervous system highlights molecular principles of convergent evolution

Abstract The cephalopod and vertebrate visual systems are a textbook example of convergent evolution with unknown molecular underpinnings. Here we characterize 98,537 single-cell transcriptomes in the bobtail squidEuprymna berryito understand how the cephalopod retina and optic lobes relate to the vertebrate retina. We confirm the overall relative simplicity of the cephalopod retina but identify two related photoreceptor cell subtypes expressing distinct r-opsins. By contrast, the adult optic lobe contains a diverse repertoire of neuronal and glial cell types, with a predominance of dopaminergic neurons. We show that cephalopod-specific gene duplicates probably contributed to this cell type diversification. Comparing neuronal cell population in the optic lobes of hatchlings and adults, we reveal a switch towards dopaminergic neurotransmitter usage with age, indicative of a maturation process. We further identify an FMRF-amide-based retrograde signal from the optic lobe towards the retina that supports the functional analogy of the cephalopod optic lobe cortex and the vertebrate inner retina in visual signal processing from a molecular standpoint. Finally, comparative analyses with vertebrate and arthropod cells suggest a scenario in which two photoreceptor types and two neuronal populations may have already been present in the eye of the bilaterian ancestor.

Environmental Sciences & Ecology↗
DOE OSTI2022

Molecular-scale mechanisms of CO 2 mineralization in nanoscale interfacial water films

The calamitous impacts of untethered carbon emission from fossil-fuel-burning energy infrastructure calls for accelerated development of large-scale CO 2 capture, utilization, and storage technologies that are underpinned by fundamental understanding of molecular-level chemical processes. In the subsurface, rocks rich in divalent metals can react with CO 2 , permanently sequestering it in the form of stable metal carbonate minerals, with the CO 2 -H 2 O composition of the post-injection pore fluid acting as a primary control variable. Herein, we compare mechanistic reaction pathways for aqueous-mediated carbonation with carbon mineralization occurring in nanoscale adsorbed water films. In the extreme of pores filled with a CO 2 -dominant fluid, carbonation reactions are confined to nanometer-thick water films coating mineral surfaces, which enable metal cation release, transport, nucleation, and crystallization of metal carbonate minerals. Though seemingly counterintuitive, laboratory studies have demonstrated facile carbonation rates in these low water environments, which in recent years has begun to be better understood in mechanistic detail. The overarching objective of this review is to delineate the unique underlying molecular-scale reaction mechanisms that govern CO 2 mineralization in these reactive and dynamic quasi-2D interfaces. Here. we highlight the importance of understanding unique properties in thin water films, such as how water dielectric properties, and consequently ion solvation/hydration behavior, can change under nanoconfinement. We conclude by identifying important frontiers for future work and opportunities to exploit these fundamental chemical insights for decarbonization technologies in the 21st century.

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

Molecular Transformation of Dissolved Organic Matter during the Processes of Wildfire, Alum Coagulation, and Disinfection Using ESI(-) and ESI(+) FT-ICR MS

Forest fires alter the quality of terrestrial dissolved organic matter (DOM) flushed into source water, which could be further changed during alum coagulation and serve as potential disinfection byproduct (DBP) precursors in water facilities. To unveil this molecular transformation, we applied 12T Fourier transform ion cyclotron resonance mass spectrometry on a controlled study which used the leachates of black and white ashes (collected immediately after the 2013 Rim Fire, CA) along with unburned white fir (Abies concolor) and ponderosa pine (Pinus ponderosa). Here, results showed that the assigned features in DOM from burned ashes had a smaller molecular weight (MW), contained more unsaturated, aromatic, and oxidized and less fraction of labile compounds, and were enriched with nitrogen, relative to unburned litters. After alum coagulation, we observed the disappeared features with higher MW and S/C and N/C atom ratios than the resistant features. During the following disinfections, the disappeared features serving as the potential DBP precursors had higher MW and unsaturation degree, contained more S- and N-containing functionalities, and were more aromatic, compared to the resistant features. Our study suggested organic compounds with small MW and high aromaticity as the targets for water facilities to improve post-fire drinking water quality.

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

Approaching Angstrom-Scale Resolution in Lithography Using Low-Molecular-Mass Resists (<500 Da)

Resists that enable high-throughput and high-resolution patterning are essential in driving the semiconductor technology forward. The ultimate patterning performance of a resist in lithography is limited because of the trade-off between resolution, line-width roughness, and sensitivity; improving one or two of these parameters typically leads to a loss in the third. As the patterned feature sizes approach angstrom scale, the trade-off between these three metrics becomes increasingly hard to resolve and calls for a fundamental rethinking of the resist chemistry. Low-molecular-mass monodispersed metal-containing resists of high atom economy can provide not only very high resolution but also very low line-width roughness without sacrificing sensitivity. Here we describe a modular metal-containing resist platform (molecular mass <500 Da) where a molecular resist consists of just two components: a metal and a radical initiator bonded to it. This simple system not only is amenable to high-resolution electron beam lithography (EBL) and extreme ultraviolet lithography (EUVL) but also unites them mechanistically, giving a consolidated perspective of molecular and chemical processes happening during exposure. Irradiation of the resist leads to the production of secondary electrons that generate radicals in the initiator bonded to metal. This brings about an intramolecular rearrangement and causes solubility switch in the exposed resist. We demonstrate record 1.9–2.0 nm isolated patterns and 7 nm half-pitch dense line-space features over a large area using EBL. With EUVL, 12 nm half-pitch line-space features are shown at a dose of 68 mJ/cm 2 . In both of these patterning techniques, the line-width roughness was found to be ≤2 nm, a record low value for any resist platform, also leading to a low-performance trade-off metric, Z factor, of 0.6 × 10 –8 mJ·nm 3 . With the ultimate resolution limited by instrumental factors, potential patterning at the level of a unit cell can be envisaged, making low-molecular-mass resists best poised for angstrom-scale lithography.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Spectroscopic Investigations of Molecular Symmetry Breakdown

This project addressed precise measurements of the spectra of molecules to address fundamental questions in chemical physics in the realm of nuclear spin symmetries and their breakdown in molecules, and parity non-conservation caused by the weak nuclear force. Nuclear spin symmetry is intimately connected to parity in symmetrical molecules and plays a central role in the physical and chemical properties of symmetric molecules, their collisional relaxation, and the establishment of thermodynamic equilibrium in an ensemble. Such effects are present in every molecular system or process, but their small size means they are often only directly observable under special conditions. However, their measurement is of fundamental, and often practical importance, and strongly relates to the themes of coherence in light and matter in the recent Basic Energy Sciences Advisory Committee (BESAC) report “Challenges at the Frontiers of Matter and Energy: Transformative Opportunities for Discovery Science”. We begin by addressing lifetime broadening of spectroscopic lines caused by self-collisions in acetylene, where efficient resonant rotational energy transfer between a pair of molecules separated in energy by rotational quantum number ΔJ = ± 2, can result in a different pressure broadening for ortho- and para- nuclear spin state molecules. This work was prompted by a controversial report [Iwakuni et al. Phys. Rev Letts. 117 (2016) 143902] of the observation of such an effect in the v 1 +v 3 combination band transition of acetylene discussed in the original proposal. After delays caused by the COVID-19 pandemic, we have completed these measurements and a paper was published in the Journal of Chemical Physics describing the work. The measurements pointed out that the use of an approximate model in the analysis of the original work was the cause of the erroneous conclusion that ortho and para acetylene molecules had substantially different collisional lifetime broadening. In the process, much improved broadening measurements have been made that will permit more reliable modeling of the effects in the future. Subsequently, we have completed a new set of measurements of self- and nitrogen-collisional broadening at a series of low and intermediate temperatures to characterize the temperature dependence of the effects in the molecule. This work is currently being analyzed and will be the subject of a paper to be published in the relatively near future. In addition, some very preliminary measurements of CH 4 as described in the proposal, have been completed. If facilities and time permit, these will continue as an undergraduate research project.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Non-Covalent Interactions and Helical Packing in Thiophene-Phenylene Copolymers: Tuning Solid-State Ordering and Charge Transport for Organic Field-Effect Transistors

In this study, we introduce two thiophene-phenylene-thiophene (TPT) polymers designed to leverage noncovalent intramolecular interactions to regulate main-chain conformation and enhance solid-state ordering. By incorporating unsubstituted thiophene (T) or bithiophene (2T) units, we reveal striking divergence in the thermal, morphological, and optoelectronic properties of the resulting films, facilitated by these noncovalent interactions. Using a combination of computational and experimental approaches, we show that annealing yields remarkably different polymer conformations and, consequently, charge transport properties. TPT-T undergoes a significant structural transformation, adopting a more planar backbone conformation and a highly crystalline, edge-on molecular orientation. In contrast, the introduction of a single additional thiophene unit in TPT-2T leads to a more isotropic molecular orientation with a slight preference for face-on alignment, resulting in a heterogeneous film structure that hinders charge transport despite achieving tighter molecular packing. Remarkably, despite being composed of achiral components, TPT-2T develops chirality upon annealing, indicating the formation of a helical conformation. Organic field-effect transistor measurements reveal that the well-ordered alignment in annealed TPT-T films results in higher charge carrier mobility and a narrower distribution of mobility values than in TPT-2T. These findings provide critical insights into the structure−property relationships of conjugated polymers, offering guidance for optimizing molecular design and processing strategies for highperformance organic electronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spontaneous assembly of hierarchical phases

We report the grand promise of self-assembly is the creation of target nanostructures through spontaneous molecular organization, a process that is — in principle — simpler and more scalable than laborious top-down nanofabrication methods. Yet self-assembly faces a fundamental challenge: the structures that form thermodynamically tend to minimize surface area and maximize symmetry; the most readily-accessible structures are thus quite simple. Block copolymer (BCP) self-assembly, for example, takes advantage of the chemical incompatibility along a polymer backbone to drive local phase separation; frustrated assembly leads to a nanoscale morphology whose unit cell is dictated by chain architecture. Canonical BCPs (linear chains with two blocks) form alternating lamellae and hexagonally-packed cylinders in a straightforward manner; but more exotic phases require additional control, for example directing the assembly using external forces, such as shear or chemical templates. An elegant alternative is to use synthetic methods to design more sophisticated chain architectures, which may give rise to new morphologies. An outstanding challenge is to achieve morphological control without undue synthetic burden. Sun et al. now report in Nature Nanotechnology the design and synthesis of novel triblock copolymers that give rise to structural order at multiple length scales, thereby forming complex and hierarchical nanoscale phases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2025

Thermally Evaporated Naphthalene Diimides as Electron Transport Layers for Perovskite Solar Cells

Thermally evaporated organic electron transport layers (ETLs) have the potential to enable high-performance and scalable perovskite solar cells (PSCs). Among these, naphthalene diimide (NDI)-based ETLs are a promising family of materials that exhibit the optoelectronic properties, ambient stability and versatility required of high-performance ETLs. Here, we synthesized five NDI derivatives with varying functional groups and identified the two most promising candidates for evaluating the impact of molecular structure on processability via thermal evaporation. While phosphonic acid functionalization was shown to introduce thermal instability, leading to chemical changes during evaporation, NDI-bis N-phenyl-bromide (NDI-(PhBr) 2 ) emerged as a promising ETL candidate. NDI-(PhBr) 2 demonstrated excellent compatibility with the thermal evaporation process and enabled PSCs with power conversion efficiencies (PCEs) of 15.6%, surpassing all previously reported PSCs containing thermally evaporated NDI ETLs. Furthermore, NDI-(PhBr) 2 exhibited excellent operational stability, retaining 75% of the initial PCE after 150 h of operation under continuous illumination at 65 °C. These results highlight the potential of NDI-based ETLs for advancing the scalability and performance of PSCs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2025

Entropy-Driven Porous Liquids Allowing Gas Solubility in Solvent-Filled Imine-Based Porous Organic Cages

Porous liquids offer a promising platform for gas separation by combining fluid processability with intrinsic molecular porosity. Traditional Type II porous liquids are formed by dissolving porous molecular cages in size-excluded solvents, limiting solvent options and practical applications. In this work, we introduce a novel method of creating Type II porous liquids using common small solvents, where intrinsic porosity is achieved at elevated pressures due to the selective displacement of solvent molecules by gas molecules within the cage structures. Using molecular simulations, we investigate the behavior of CO 2 in solutions of the imine-based porous organic cage CC13 dissolved in small molecular solvents such as chloroform and 1,2-dimethoxyethane (DME). Grand canonical Monte Carlo simulations of solid-state CC13 reveal that selectivity reversal, where smaller CO 2 molecules displace larger solvent molecules inside the cage, is achievable at sufficiently high pressures. Temperature quench molecular dynamics simulations confirm that while CO 2 displacement within chloroform-filled cages is limited, DME enables entropy-driven cage CO 2 occupancy at pressures as low as ∼23 bar, setting up the foundation of an alternative way of forming Type II porous liquids.

adsorption↗
DOE OSTI2022

Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis beyond the relative scales typically employed. In principle, theoretical methods can provide such quantification but are often hindered by the innate complexity of strong electron correlation and dynamic relaxation of solvent environments. Here, we present and validate a protocol for calculating the redox potentials of ferrocene/ferrocenium redox pair in the acetonitrile. Equation-of-motion ionization potential coupled-cluster single-double (EOM-IP-CCSD) and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials (IP) as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrated the ability of modern computational methods to predict absolute redox potentials quantitatively and, more critically, quantify the correlation of dynamic effects, which underlie their origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗