SEARCH · Search NASA
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.
A novel multi-electrode glow discharge provides increased sensitivity and analytical flexibility for both atomic and molecular MS determinations.
The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.
The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.
Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.
Porous liquids (PLs) are an attractive material for gas separation and carbon sequestration due to their permanent internal porosity and high adsorption capacity. PLs that contain zeolitic imidazole frameworks (ZIFs), such as ZIF-8, form PLs through exclusion of aqueous solvents from the framework pore due to its hydrophobicity. The gas adsorption sites in ZIF-8 based PLs are historically unknown; gas molecules could be captured in the ZIF-8 pore or adsorb at the ZIF-8 interface. To address this question, ab initio molecular dynamics was used to predict CO 2 binding sites in a PL composed of a ZIF-8 particle solvated in a water, ethylene glycol, and 2-methylimidazole solvent system. Further, the results show that CO 2 energetically prefers to reside inside the ZIF-8 pore aperture due to strong van der Waals interactions with the terminal imidazoles. However, the CO 2 binding site can be blocked by larger solvent molecules that have greater adsorption interactions. CO 2 molecules were unable to diffuse into the ZIF-8 pore, with CO 2 adsorption occurring due to binding with the ZIF-8 surface. Therefore, future design of ZIF-based PLs for enhanced CO 2 adsorption should be based on the strength of gas binding at the solvated particle surface.
Oriented attachment is a critical, yet poorly understood, crystal growth pathway based on the self-assembly of nanocrystals. During oriented attachment, solvent-separated particles align and coalesce through forces that enable precise rotation and translation. While prior studies emphasized intragap forces driving crystallographic alignment, the forces enabling uniform stacking and superlattice formation remain unclear. Here, we demonstrate how macroscopic gibbsite mesocrystals emerge from nanoplates guided into staggered positions by directional sliding. Electron microscopy and X-ray scattering reveal the monoclinic superlattice structure, based on nanoplate stacking with a uniform ≈50° stagger along the gibbsite [010] direction. In situ liquid-cell TEM captures preferential sliding along the gibbsite [010] direction, decelerating with increasing particle overlap. Molecular dynamics simulations reveal that this staggered arrangement corresponds to a global free-energy minimum, rather than full alignment. The simulations also confirm that sliding along the [010] direction is energetically favored and provide insight into the role of interfacial water in achieving long-range ordered assemblies. These insights highlight the energy landscape’s role in oriented attachment, with implications for material synthesis and hierarchical structures in nature.
It is critical to find methods to control the thermodynamic driving force for photoexcited charge transfer from quantum dots (QDs) and explore how this affects charge transfer rates, since the efficiency of QD-based photovoltaic and photocatalysis technologies depends on both this rate and the associated energetic losses. In this work, we introduce a single-pot shell growth and Cu-catalyzed cation exchange method to synthesize Cd x Zn 1-x Se/Cd y Zn 1-y S QDs with tunable driving forces for electron transfer. Functionalizing them with two molecular electron acceptors—naphthalenediimide (NDI) and anthraquinone (AQ)—allowed us to probe nearly 1 eV of driving forces. For AQ, at lower driving forces, we find that higher Zn content results in a 130-fold increase of electron transfer rate constants. However, at higher driving forces electron transfer dynamics are unaltered. Here, the data are understood using an Auger-assisted electron transfer model and analyzed with computational work to determine approximate binding geometries of these electron acceptors. Our work provides a method to tune QD reducing power and produces useful metrics for optimizing QD charge transfer systems that maximize rates of electron transfer while minimizing energetic losses.
Hypothesis. Understanding the mechanisms of proton transfer on quartz surfaces in water is critical for a range of processes in geochemical, environmental, and materials sciences. The wide range of surface acidities (>9 pKa units) found on the ubiquitous mineral quartz is caused by the structural variations of surface silanol groups. Molecular scale simulations provide essential tools for elucidating the origin of site-specific surface acidities. Simulations. Here, we used density-functional tight-binding-based molecular dynamics combined with rare-event metadynamics simulations to probe the mechanisms of deprotonation reactions from ten representative surface silanol groups found on both pristine and defect-rich quartz (1 0 1) surfaces with Si vacancies. Findings. The results show that deprotonation is a highly dynamic process where both the surface hydroxyls and bridging oxygen atoms serve as the proton acceptors, in addition to water. Deprotonation of embedded silanols through intrasurface proton transfer exhibited lower pKa values with less H-bond participation and higher energy barriers, suggesting a new mechanism to explain the bimodal acidity observed on quartz surface. Defect sites, recently shown to comprise a significant portion of the quartz (1 0 1) surface, diversify the coordination and local H-bonding environments of the surface silanols, changing both the deprotonation pathways and energetics, leading to a wider range of pKa values (2.4 to 11.5) than that observed on pristine quartz surface (10.4 and 12.1).
The primary questions motivating this report are: Are there ways to increase coherence and delocalization of excitation among many molecules at moderate electronic coupling strength? Coherent delocalization of excitation in disordered molecular systems is studied using numerical calculations. The results are relevant to molecular excitons, polaritons, and make connections to classical phase oscillator synchronization. In particular, it is hypothesized that it is not only the magnitude of electronic coupling relative to the standard deviation of energetic disorder that decides the limits of coherence, but that the structure of the Hamiltonian—connections between sites (or molecules) made by electronic coupling—is a significant design parameter. Inspired by synchronization phenomena in analogous systems of phase oscillators, some properties of graphs that define the structure of different Hamiltonian matrices are explored. The report focuses on eigenvalues and ensemble density matrices of various structured, random matrices. Some reasons for the special delocalization properties and robustness of polaritons in the single-excitation subspace (the star graph) are discussed. The key result of this report is that, for some classes of Hamiltonian matrix structure, coherent delocalization is not easily defeated by energy disorder, even when the electronic coupling is small compared to disorder.
Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.
The use of pressure to obtain new materials that can be recovered under ambient conditions is a central problem in high-pressure physics. Despite decades of research, this goal has only been achieved in the laboratory for a few notable examples, such as diamond and cubic boron nitride. An area of significant interest is the transformation under compression of light-element molecular compounds to extended covalent-bonded (polymeric) solids. Among them, CO 2 has been extensively studied because of its status as a prototypical simple molecular system with a rich phase diagram and due to its fundamental role in Earth’s physics and chemistry. One of its polymeric crystalline phases, accessible at extreme pressures and temperatures, has been recently quenched to ambient pressure, but below room temperature. Here we report ab initio calculations predicting that isothermal compression of a carbon monoxide and oxygen mixture (CO+O 2 ), rather than the compound CO 2 , lowers the onset of C-polymerization at room temperature from ~ 118 GPa to ~ 7 GPa (complete by ~ 23 GPa). Moreover, it leads to the formation of an intrinsically different polymer with enhanced metastability. We predict that this dense phase is an energetic material which can potentially be recovered to ambient pressure and temperature.
Positive ion desorption following electron impact dissociative ionization of ethane adsorbed on Ar, Kr, and Xe multilayers has been studied as a function of incident electron energy from threshold to 100 eV. Based on the dependence of ion yields on the identity of the rare gas, it is likely that the majority of ethane molecules undergo indirect ionization following hole transfer from the ionized underlying rare gas. Further, this has also been corroborated by density of states calculations showing the energetic alignment of the outer valence states of ethane and the condensed rare gas ionization energies. Due to the near-resonant nature of charge transfer for single-hole states, the ethane molecular ion is excited to different final ionic states on different rare gases, which leads to differences in ion desorption yields and branching ratios. The quantitative yields increase with increasing ionization energy gap between the rare gas and ethane, in the order Ar > Kr >Xe. The large increase in yields from 25 eV onwards for all rare gases is likely due to the formation and decay of two-hole states on neighboring rare gas and ethane molecules due to interatomic and intermolecular Coulomb decay (ICD) and not electron transfer mediated decay (ETMD). The ICD and ETMD pathways become accessible when the incoming electron has sufficient energy to excite the inner valence ns level of the rare gas to a Rydberg state or ionize it. The experimental findings are supported by calculations of thresholds, density of states for the final configurations of these processes, and coupling strengths for hole transfer between ethane and rare gases. The fragment ion branching ratios vary with energy from threshold to about 35 eV, showing the fragmentation pattern changes with the mode of hole transfer and availability of excess energy. Sigma C–C bonds are more likely to break than C–H bonds in the mid-20 eV range, and this effect is most pronounced for Xe, followed by Kr, and then Ar.
Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.
Understanding the adsorption and reactivity of carboxylic acids on oxide surfaces is of great interest in catalysis for biomass upgrading via ketonization, a carbon–carbon coupling reaction. Herein, we investigate the adsorption and reaction of acetic acid on anatase TiO 2 (101) using scanning tunneling microscopy, infrared spectroscopy, temperature programmed reaction, and density functional theory calculations. We demonstrate the adsorption of acetic acid can form two intermediates: (1) dissociated, bidentate acetate with an associated bridging hydroxyl, and (2) molecular, monodentate acetic acid. The coexistence of ordered phases with increasing monolayer (ML) saturation coverages consisting of (1) pure acetate (0.5 ML), (2) mixed acetate/acetic acid (0.67 ML), (3) mixed acetate/acetic acid (1.0 ML) and (4) pure acetic acid demonstrates similar energetics for both acetate and acetic acid species. Under ultra-high vacuum conditions, the presence of both monodentate acetic acid and bidentate acetate was observed below room temperature, while solely bidentate acetate was observed up to 575 K. The deprotonation of acetic acid produces water at 280 K, while the thermal decomposition of bidentate acetate produces ketene and acetic acid at 645 K. In conclusion, this model study provides detailed insight into the stability and reactivity of carboxylic acid surface-bound intermediates, which could participate during ketonization reactions for biomass upgrading.
Concerns surround whether insensitive (or any) energetic materials are more dangerous to handle when exposed to abnormal thermal environments. This study characterizes the residual material remaining after LX-17 (92.5 % 1,3,5-triamino 2,4,6-trinitro benzene (TATB) and 7.5 % Kel-F) is exposed to various thermal environments in a sealed small-scale vessel cook-off test reactor (heated at 0.1 to 100 °C/min until the reactor opened at 3000 psi (20.7 MPa)). Previous work has shown no additional sensitivity of these residues as evaluated by small-scale safety analysis, but characterization on the molecular scale indicates the TATB is transformed to more reactive compounds as well as the residue could be precursors to toxic gases. The solids and chars were characterized by various analytical methods. Heat-flow measurements indicated exothermic release is due to a mixture of residual TATB and related decomposition products (which may be more energetic). The N/C and O/N ratios indicated a material much more degraded than TATB. Primarily, the solids were a network of amorphous C inter-dispersed with N and O. Types of bonding include C–C, C–N, N–H, N–C, N=C, N≡C, C–O=, and –OH. Solvent extracts of the solids showed TATB decomposition intermediates benzo-furazans and benzo-furoxans, substituted TATB (mono-nitroso, hydroxyl, and chlorinated) along with several unidentified smaller molecules. These results indicate thermal treatment produces an amorphous carbon residue with heteroatoms incorporated through differing functionality, varying depending upon the thermal severity of exposure. These structures also could further decompose producing toxic light gases (such as cyanide).
Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.
Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.
Simulations of optical spectra can provide key insights to aid experimental interpretation of electronic excitation phenomena. For chromophores in the condensed phase, these spectra, which incorporate the coupling between electronic excitation and molecular and solvent nuclear motions, can be simulated using excitation energies obtained from molecular dynamics simulations of the chromophore and solvent. Here, we present a hybrid scheme that exploits machine learning and physically informed spectral densities to show that as few as 25 ground and excited state energetic gradient calculations can be used to construct models that accurately predict environment-influenced vibronic coupling in optical spectra. We demonstrate our approach for the green fluorescent protein chromophore in water and the cresyl violet chromophore in methanol. We show that our hybrid approach, employing a machine learning model for the high-frequency spectral density and an ab initio parameterized Debye spectral density for the low-frequency, results in a systematic improvement of the optical absorption lineshape, leading to a simple machine learning scheme that can be used for the simulation of spectral densities and optical spectra.