SEARCH · Engineering Papers
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.
Cation and anion vacancy formation energies are calculated for HgTe, ZnTe, CdTe, and their dilute alloys. Harrison's tight-bonding theory is extended to a cluster embedded in an extended crystal. Only neutral vacancies have been considered, and two final states for the removed atom have been addressed: a free atom in vacuum and an atom on an ideal (111) surface.
Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 – ) or dimeric (Al 2 O(OH) 6 2– ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 – and Al 2 O(OH) 6 2– anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2– anions, while the third contains the monomeric Al(OH) 4 – anion. Here, comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.
A synergistic anion photoelectron spectroscopic and ab initio computational study of photodetachment of UF – is reported. The measurement determined a vertical detachment energy of 0.63(03) eV, which is consistent with a spinor-based relativistic coupled-cluster CCSD(T) value of 0.61 eV. The complex spectral features due to excited electronic states and vibrational progressions of UF are analyzed and assigned with the help of spin–orbit-coupled multireference perturbation theory and spinor-based relativistic coupled-cluster calculations. UF and UF – are confirmed to be dominated by ionic bonding. Furthermore, the usefulness of the spinor CCSD(T) approach is demonstrated.
Abstract The current study of the small lithium molecules Li 2 +,0,− and Li 3 +,0,− focuses on the nature of the bonding in these molecules as well as their structures and energetics (bond energies, ionization energies, and electron affinities). Valence CASSCF (2s,2p) calculations incorporate nondynamical electron correlation in the calculations, while the corresponding multireference configuration interaction and coupled cluster calculations incorporate dynamical electron correlation. Treatment of nondynamical correlation is critical for properly describing the Li 2,3 +,0,− molecules as well as the Li − anion with dynamical correlation, in general, only fine‐tuning the predictions. All lithium molecules and ions are bound, with the Li 3 + and Li 2 + ions being the most strongly bound, followed by Li 3 − , Li 2 , Li 2 − and Li 3 . The minimum energy structures of Li 3 +,0,− are, respectively, an equilateral triangle, an isosceles triangle, and a linear structure. The results of SCGVB calculations are analyzed to obtain insights into the nature of the bonding in these molecules. An important finding of this work is that interstitial orbitals, a concept first put forward by McAdon and Goddard in 1985, play an essential role in the bonding of all lithium molecules considered here except for Li 2 . The interstitial orbitals found in the Li 3 +,0 molecules likely give rise to the non‐nuclear attractors/maxima observed in these molecules.
Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.
Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.
Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.
Sub-nanometer metal oxide clusters are very important materials that are widely used, for example, in catalysis or electronic devices such as sensors. Hence, it is critical to understand the atomic structures and properties of sub-nanometer metal oxide clusters under a reactive gas environment, such as O 2 . We consider here experimentally accessible precise-size Cu clusters (Cu 4 ) supported on partially hydroxylated amorphous alumina and show that such clusters can access, in catalytic conditions at high temperature under a pressure of O 2 , a large ensemble of oxidized structures, representing a large variety of oxygen content, of geometries in link with the support, and of catalytic activities for oxidation reactions, as seen from their reducibilities. A grand canonical basin hopping method based on first-principles energy reveals an ensemble of 24 configurations for the Cu 4 O x cluster of low free energy, less than 0.8 eV above the global minimum. The low free energy ensemble consists of clusters of different stoichiometries, which are mainly Cu 4 O 3 and Cu 4 O 4 in the temperature range of 200–400 °C and under a pressure of 0.5 bar of O 2 . The presence of several competitive isomers at each composition implies that cluster fluxionality impacts the phase diagram, which should be ensemble-averaged. In terms of catalytic oxidation activity, Cu 4 O 3 isomers present highly variable O abstraction energy: the most stable isomers are inactive for alkane oxidative dehydrogenation, but isomerization to metastable isomers, that proceed with low barrier, enable to create active configurations with low O abstraction energy. O atoms with the lowest anionic character, and thus of more electrophilic nature, present the best oxidation capability. In contrast, all Cu 4 O 4 isomers show a low O abstraction energy and a high potential catalytic activity. Here, this manuscript demonstrates the unique structural and electronic properties of sub-nano Cu oxide clusters and illustrates the critical roles of configuration ensembles and rearrangement to highly reactive metastable cluster isomers in nanocatalysis.
Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.
Heterometallic clusters have attracted broad interests in the synthetic chemistry due to their various coordination modes and potential applications in heterogeneous catalysis. Here we report the synthesis, experimental, and theoretical characterizations of four ternary clusters ([M 2 (CO) 6 Sn 2 Sb 5 ] 3- (M = Cr, Mo), and [(MSn 2 Sb 5 ) 2 ] 4- , (M = Cu, Ag)) in the process of capturing the hypho - [Sn 2 Sb 5 ] 3- in ethylenediamine (en) solution. We show that the coordination of the binary anion to transition-metal ions or fragments provides additional stabilization due to the formation of locally σ-aromatic units, producing a spherical aromatic shielding region in the cages. While in the case of [Mo 2 (CO) 6 Sn 2 Sb 5 ] 3- stabilization arises from locally σ-aromatic three-centre and five-centre two-electron bonds, aromatic islands in [(AgSn 2 Sb 5 ) 2 ] 4- and [(CuSn 2 Sb 5 ) 2 ] 4- render them globally antiaromatic. This work describes the coordination chemistry of the versatile building block [Sn 2 Sb 5 ] 3- , thus providing conceptual advances in the field of metal-metal bonding in clusters.
The design of cathode materials that remain chemically and structurally stable during repetitive ion insertion and extraction poses a significant challenge in developing multivalent batteries. The cycling stability of traditional metal oxide-based cathode is challenged by sluggish diffusion of multivalent cations and parasitic reactivity at interfacial regimes, including the cathode electrolyte interphase layer (CEI). Understanding the reactions at the cathode-electrolyte interface, particularly those induced by non-stoichiometric surface layers, is a crucial design parameter for both cathode materials and electrolytes. Here, in this study, we employed multimodal analysis, including in situ and ex situ X-ray photoelectron spectroscopy (XPS), high resolution transmission electron microscopy (TEM) and electrochemical impedance spectroscopy (EIS) to examine the surface reactions and subsequent structural and chemical evolutions of the CEI on high voltage magnesium vanadium oxide (MgV 2 O 4 ) spinel cathode during the Mg 2+ insertion/extraction processes. The results revealed that the presence of non-stoichiometric surface layers in the magnesium vanadium oxide cathode drive the decomposition of bis(trifluoromethanesulfonyl)imide (TFSI - ) anion, leading to the formation of the CEI layer. The CEI layer could inhibit the Mg 2+ ion transfer processes. Accompanying this reactivity-driven degradation, the magnesium vanadium oxide cathode undergoes pulverization, forming clusters of nanosized particles. This process likely improves cycling ability by creating new intercalation sites and shortening the diffusion pathway for the Mg 2+ cations. This study demonstrates that controlling surface stoichiometry and engineering morphological properties are critical design parameters for high performance cathodes for multivalent batteries.
Photolysis reaction pathways of [Au(III)Cl 4 ] - in aqueous solution have been investigated by time-resolved X-ray absorption spectroscopy. Ultraviolet excitation directly breaks the Au-Cl bond in [Au(III)Cl 4 ] - to form [Au(II)Cl 3 ] - that becomes highly reactive within 79 ps. Disproportionation of [Au(II)Cl 3 ] - generates [Au(I)Cl 2 ] - , which is stable for ≤ 10 mu s. In contrast, intense near-infrared lasers photolyze water to generate hydrated electrons, which then reduce [Au(III)Cl 4 ] - to [Au(II)Cl 3 ] - at 5 ns. Hydrated electrons further induce a chain reaction from [Au(II)Cl 3 ] - to [Au(0)Cl] - by successively removing one Cl - . The zero-valency Au anions quickly polymerize and condense to form Au nanoparticles, which become the dominating product after 400 s. Here our results reveal that the condensation of zero-valency Au starts with dimerization of gold clusters coordinated with chloride ions rather than direct condensation of pristine Au atoms.
Quinone-based sorbents offer a tunable, energy-efficient route to electrochemical CO2 capture, but systematic guidance for molecular design is lacking. Here, we report a high-throughput computational workflow that combines density functional theory (DFT) screening with machine-learning (ML) modeling to evaluate CO2 binding thermodynamics across several quinone derivatives, spanning benzoquinones, naphthoquinones, and anthraquinones. In addition to using solvents to stabilize the quinone anion and dianion, we studied the effect of ion-pairing on the reduction potentials and the CO2 binding energy. Automated Python scripts handled geometry optimizations and adduct-formation energies on an HPC cluster, reducing manual effort significantly. This integrated platform can uncover structure–property relationships and enables rapid in silico evaluation of untested candidates. We present one example from our workflow to showcase the capability of using quinones with ion-pairing to effectively capture CO2. Our approach paves the way for the rational selection of optimal quinone sorbents and can be extended with experimental thermochemical and kinetic data, alternative redox cycles, and stability assessments to accelerate development of next-generation electrochemical CO2 capture materials.
Abstract. In the lower troposphere, rapid collisions between ions and trace gases result in the transfer of positive charge to the highest proton affinity species and negative charge to the lowest proton affinity species. Measurements of the chemical composition of ambient ions thus provide direct insight into the most acidic and basic trace gases and their ion–molecule clusters – compounds thought to be important for new particle formation and growth. We deployed an atmospheric pressure interface time-of-flight mass spectrometer (APi-ToF) to measure ambient ion chemical composition during the 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land Ecosystems (HI-SCALE) campaign at the United States Department of Energy Atmospheric Radiation Measurement facility in the Southern Great Plains (SGP), an agricultural region. Cations and anions were measured for alternating periods of ∼ 24 h over 1 month. We use binned positive matrix factorization (binPMF) and generalized Kendrick analysis (GKA) to obtain information about the chemical formulas and temporal variation in ionic composition without the need for averaging over a long timescale or a priori high-resolution peak fitting. Negative ions consist of strong acids including sulfuric and nitric acid, organosulfates, and clusters of NO3- with highly oxygenated organic molecules (HOMs) derived from monoterpene (MT) and sesquiterpene (SQT) oxidation. Organonitrates derived from SQTs account for most of the HOM signal. Combined with the diel profiles and back trajectory analysis, these results suggest that NO3 radical chemistry is active at this site. SQT oxidation products likely contribute to particle growth at the SGP site. The positive ions consist of bases including alkyl pyridines and amines and a series of high-mass species. Nearly all the positive ions contained only one nitrogen atom and in general support ammonia and amines as being the dominant bases that could participate in new particle formation. Overall, this work demonstrates how APi-ToF measurements combined with binPMF analysis can provide insight into the temporal evolution of compounds important for new particle formation and growth.
The oxygen-evolving complex (OEC) of photosystem II (PSII) is an oxomanganese cluster composed of four redox-active Mn ions, and one redox-inactive Ca 2+ ion, with two nearby bound Cl – ions. Sodium is a common counter-ion of both chloride and hydroxide anions, and a sodium-specific binding site has not been identified near the OEC. Furthermore, we find that the oxygen-evolution activity of spinach PSII increases with [Na + ], particularly at high pH. A Na + -specific binding site next to the OEC, upon deprotonation of the D1-H337 amino-acid residue, is suggested by the analysis of two recently published PSII cryo-EM maps in combination with quantum mechanical calculations and molecular dynamics simulations.
This study provides a comprehensive structural characterization of commercially available alkaline anion-exchange polymers (Fumasep FAA-3 and IRA 900) used in moisture-driven direct air capture (DAC) of carbon dioxide. Using X-ray diffraction, SAXS/WAXS, atomic force microscopy, FIB-SEM, and transmission electron microscopy, the authors identify nanoscale clustering, porosity, swelling behavior, and humidity-dependent structural changes that influence CO₂ adsorption and release. These findings establish structure–function relationships critical for designing more durable and energy-efficient DAC polymer materials.
Iron-sulfur clusters play important roles in biology as parts of electron transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we report results for the ground and low-lying electronically excited states of the complex [Fe(SCH3)4]2-/1-/2+/3+, using Multi-Reference (CASSCF, MRCISD), and Coupled Cluster [CCSD(T)] methodology in order to provide accurate adiabatic reduction energies, dissociation energies and insights into the bonding analysis. The nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH3)4] complexes offer a physical rationale for the relative stabilization, structure and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH3)4]2-/1- complexes and covalent in [Fe(SCH3)4]2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH3)4]2-®[Fe(SCH3)4]- (i.e., FeII®FeIII) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH3)4]2+®[Fe(SCH3)4]3+ of 16.6 eV, almost half of the IP of Fe2+. The ionic vs covalent bond character influences the Fe-S bond strength and length, i.e., ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for FeII) and 0.04 Å (for FeII). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (FeII) and 14.6 (FeIII) eV at the RCCSD(T) level of theory.
Anion photoelectron spectra of US – and US 2 – were recorded using the third (355 nm) and fourth (266 nm) harmonics of an Nd:YAG laser, which yielded vertical detachment energies (VDEs) of 1.71 and 2.02 eV, respectively. The experimental results are supported by extensive relativistic ab initio calculations, primarily at the coupled cluster level of theory, with systematic sequences of correlation consistent basis sets. Calculations include the closely related SUO and SUO – molecules, as well as the oxide congeners UO/UO – and UO 2 /UO 2 – which are well-known experimentally and provide benchmark systems for the sulfide calculations. Adiabatic electron detachment energies (ADEs) are computed for UO – , UO 2 – , US – , SUO – , and US 2 – using the Feller–Peterson–Dixon (FPD) composite approach. Additionally, ADEs are determined for UO – and US – using a spinor-based coupled cluster approach where spin–orbit coupling is included at the orbital level. VDEs are derived from the ab initio results from Franck–Condon simulations of the photoelectron spectra.