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 prominent goal in quantum chemistry is to solve the molecular electronic structure problem for ground state energy with high accuracy.
Optical cycling, a continuous photon scattering off atoms or molecules, plays a central role in the quantum information science. While optical cycling has been experimentally achieved for many neutral species, few molecular ions have been investigated. We present a systematic theoretical search for diatomic molecular ions suitable for optical cycling using equation-of-motion coupled-cluster methods. Inspired by the electronic structure patterns of laser-cooled neutral molecules, we establish the design principles for molecular ions and explore various possible cationic molecular frameworks. The results show that finding a perfect molecular ion for optical cycling is challenging, yet possible. Among various possible diatomic molecules we suggest several candidates, which require further attention from both theory and experiment: YF + , SiO + , PN + , SiBr + , and BO + .
Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.
This feature focuses on bridging isolated water and bulk water studies. It assembles experiments and theory on water clusters, mainly probed by vacuum ultraviolet (VUV) radiation, and summarizes what the effects are of intermolecular interactions on both the spectroscopy and the VUV-induced processes in water. In particular, it highlights studies of heterogeneous water clusters─those incorporating other molecular species such as naphthalene, methane, formic acid, and glycerol─which serve as model systems to investigate the fundamental roles of water in hydrogen bonding networks, proton transfer, and astrochemical processes. These mixed clusters provide a platform to revisit the stability of protonated water clusters and compile observations of ionization-induced structural rearrangements and fragmentation, especially in systems involving hydroxyl-rich cosolvents. This feature then explores energy transfer mechanisms in molecular clusters, following VUV photoexcitation. It concludes with potential future directions. First, we investigate excited-state dynamics in molecular clusters through direct probing, complementing the synchrotron studies discussed here. Second, exploring gas-phase molecule evaporation from confined spaces and interfaces using advanced spectroscopic techniques sheds light on these ubiquitous, yet currently debated, molecular processes.
Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.
Alkyltin clusters are exploited in nanolithography for the fabrication of microelectronics. The alkyltin Keggin family is unique among Keggin clusters across the periodic table; its members appear to favor the lower-symmetry β and γ isomers rather than the highly symmetrical α and ε isomers. Therefore, the alkyltin Keggin family may provide important fundamental information about the formation and isomerization of Keggin clusters. In this study, we have synthesized and structurally characterized a new butyltin Keggin cluster with a tetrahedral Ca 2+ center, fully formulated [(BuSn) 12 (CaO 4 )(OCH 3 ) 12 (O) 4 (OH) 8 ] 2+ (β-CaSn 12 ). The synthesis is a simple one-step process. Extensive solution characterization including electrospray ionization mass spectrometry, small-angle X-ray scattering, and multinuclear ( 1 H, 13 C, and 119 Sn) nuclear magnetic resonance shows β-CaSn 12 is essentially phase-pure and stable. This differs from the previously reported Na-centered analogues that always form a mixture of β and γ isomers, with facile interconversion. Therefore, this study has clarified prior confusion over complex spectroscopic and crystallographic characterization of the Na-centered analogues. Density functional theory calculations showed the following stability order: γ-CaSn 12 < γ-NaSn 12 < β-CaSn 12 < β-NaSn 12 . The β analogue is always more stable than the γ analogue, consistent with experiment. Notable outcomes of this study include a rare tetrahedral Ca coordination, a Na-free alkyltin cluster (important for microelectronics manufacturing), and a better understanding of Keggin families built of different metal cations.
We present a systematic ab initio study of the low-lying states in beryllium isotopes from 7 Be to 12 Be using nuclear lattice effective field theory with the N 3 LO interaction. Our calculations achieve good agreement with experimental data for energies, radii, and electromagnetic properties. We introduce a novel, model-independent method to quantify nuclear shapes, uncovering a distinct pattern in the interplay between positive and negative parity states across the isotopic chain. By combining Monte Carlo sampling of the many-body density operator with a novel nucleon-grouping algorithm, the prominent two-center cluster structures, the emergence of one-neutron halo, complex nuclear molecular dynamics such as 𝜋 orbital and 𝜎 orbital, emerge naturally.
The nitrogenase Fe protein mediates ATP-dependent electron transfer to the nitrogenase MoFe protein during nitrogen fixation, in addition to catalyzing MoFe protein-independent substrate (CO 2 ) reduction and facilitating MoFe protein metallocluster biosynthesis. The precise role(s) of the Fe protein Fe 4 S 4 cluster in some of these processes remains ill-defined. Herein, we report crystallographic data demonstrating ATP-dependent chalcogenide exchange at the Fe 4 S 4 cluster of the nitrogenase Fe protein when potassium selenocyanate is used as the selenium source, an unexpected result as the Fe protein cluster is not traditionally perceived as a site of substrate binding within nitrogenase. The observed chalcogenide exchange illustrates that this Fe 4 S 4 cluster is capable of core substitution reactions under certain conditions, adding to the Fe protein’s repertoire of unique properties.
Supercritical fluids are characterized by unique thermodynamic properties. One of these properties is the existence of two-component dynamics that is associated with distinct low-frequency and high-frequency vibrational responses of the fluid. However, the origin of this behavior remains unknown. By combining inelastic X-ray scattering and molecular dynamics simulations, we show that this behavior can be connected to density heterogeneities arising from molecular clusters. Analyses of measurements and molecular trajectories suggest that the two-component dynamics emerges due to distinct momentum fluctuations of clustered and unbound molecules. This connection between clusters and two-component dynamics highlights the importance of molecular-structural heterogeneities in supercritical fluids, colloids, and condensed-matter systems.
The superhydrides LaH 10 , EuH 9 and UH 8 are studied theoretically using chemical bonding analysis of periodic and molecular models. Structural trends are explained using a simple electron counting model based on atomic charge calculations.
Recent progress has been made on the synthesis and characterization of metal halide perovskite magic-sized clusters (PMSCs) with ABX 3 composition ( A = C H 3 N H 3 + or Cs + , B = P b 2 + , and X = C l − , Br - , or I - ). However, their mechanism of growth and structure is still not well understood. In our effort to understand their structure and growth, we discovered that a new species can be formed without the CH 3 NH 3 + component, which we name as molecular clusters (MCs). Specifically, CH 3 NH 3 PbBr 3 PMSCs, with a characteristic absorption peak at 424 nm, are synthesized using PbBr 2 and CH 3 NH 3 Br as precursors and butylamine (BTYA) and valeric acid (VA) as ligands, while MCs, with an absorption peak at 402 nm, are synthesized using solely PbBr 2 and BTYA, without CH 3 NH 3 Br. Interestingly, PMSCs are converted spontaneously overtime into MCs. An isosbestic point in their electronic absorption spectra indicates a direct interplay between the PMSCs and MCs. Therefore, we suggest that the MCs are precursors to the PMSCs. From spectroscopic and extended X-ray absorption fine structure (EXAFS) results, we propose some tentative structural models for the MCs. The discovery of the MCs is critical to understanding the growth of PMSCs as well as larger perovskite quantum dots (PQDs) or nanocrystals (PNCs).
The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.
Metal clusters featuring carbon and sulfur donors have coordination environments comparable to the active site of nitrogenase enzymes. Here, we report a series of di-iron clusters supported by the dianionic yldiide ligands, in which the Fe sites are bridged by two μ 2 -C atoms and four pendant S donors. The [L 2 Fe 2 ] (L = {[Ph 2 P(S)] 2 C} 2– ) cluster is isolable in two oxidation levels, all-ferrous Fe 2 II and mixed-valence Fe II Fe III . The mixed-valence cluster displays two peaks in the Mössbauer spectra, indicating slow electron transfer between the two sites. The addition of the Lewis base 4-dimethylaminopyridine to the Fe2 II cluster results in coordination with only one of the two Fe sites, even in the presence of an excess base. Conversely, the cluster reacts with 8 equiv of isocyanide t BuNC to give a monometallic complex featuring a new C–C bond between the ligand backbone and the isocyanide. The electronic structure descriptions of these complexes are further supported by X-ray absorption and resonant X-ray emission spectroscopies.
Clustering plays an important role in the structure of 10 Be. Exotic molecular-like configurations, such as α:2n:α, have been suggested at relatively low excitation energies. To search for the high-spin states that may belong to the molecular-like α:2n:α configuration in 10 Be. Measuring excitation functions for 6 He + α scattering, populating states in the excitation energy range from 4.5 MeV to 8 MeV in 10 Be using a 6 He rare-isotope beam and a thick helium gas target. No new excited states in 10 Be have been observed. Stringent limitation on the possible degree of a-clustering of the hypothetical yrast 6 + state has been obtained. Here, the high-spin members of the α:2n:α molecular-like rotational band configuration, that is considered to have a 0 + bandhead at 6.18 MeV, either do not exist or have small overlap with the 6 He(g.s.)+α channel.
It has been suggested by several authors that the sodium and aluminum abundance variations correlating with CN-band strength, frequently observed in CN-bimodal globular cluster giants, could be spurious manifestations of different temperature structures in the 'CN-strong' and 'CN-weak' stars, caused by different molecular line blanketing related to the C, N, and O trio. For stellar parameters generally appropriate to giants in the intermediate metallicity CN-bimodal cluster M4, we demonstrate through new model atmosphere calculations, employing opacity sampling and spherical geometry, that the observed abundance anomalies cannot be the result of atmospheric temperature structure. Our results using spherical geometry are compared to identical calculations performed with plane-parallel geometry: the effects of atmospheric extension on derived abundances for all lines considered amount to less than 0.1 dex.
Water-lean solvents are a promising technology for capturing acid gases like carbon dioxide (CO 2 ). In situ liquid time-of-flight secondary ionization mass spectroscopy (ToF-SIMS) is used to study a representative solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (2-EEMPA) with different CO 2 loadings to reveal the complex solvent structure upon CO 2 capture. Characteristic peaks of 2-EEMPA, such as m/z – 215 C 11 H 23 N 2 O 2 – (deprotonated 2-EEMPA) and m/z + 217 C 11 H 25 N 2 O 2 + (protonated 2-EEMPA), are detected due to acid gas uptake. Also, solvent molecules and carboxylate ion pairs, such as m/z – 259 C 12 H 23 N 2 O 4 – [(deprotonated 2-EEMPA∙∙∙CO 2 )] and m/z + 261 C 12 H 25 N 2 O 4 + (protonated 2-EEMPA∙∙∙CO 2 ), are observed. Interestingly, more than one CO 2 molecule can be captured per each solvent molecule as evidenced in SIMS mass spectra, for example, m/z – 321 C 13 H 25 N 2 O 7 – [(deprotonated 2-EEMPA)∙∙∙2CO 2 ∙∙∙H 2 O], m/z + 305 C 13 H 25 N 2 O 4 + [(protonated 2-EEMPA)∙∙∙2CO 2 ], m/z – 389 C 17 H 29 N 2 O 8 – [(deprotonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙3CH 2 ], and m/z + 373 C 16 H 25 N 2 O 8 + [(protonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙2C]. However, the monomer of 2-EEMPA and CO 2 seems to be most prevalent. Furthermore, solvent clusters are detected in loaded solvents, for instance m/z + 433 C 22 H 49 N 4 O 4 + [(2-EEMPA)2∙∙∙H] and m/z + 646 [(2-EEMPA) 3 -2H], while capturing CO 2 at different amounts. Relative abundance of cluster ions provides a semi-qualitative venue to assess the free energies of gas capture energetics, indicating the relative stability trend within the same solvent system, previously impossible. These observed ion clusters are verified with molecular modeling, where dimer, trimer, and cluster ions are validated for their presence either due to weak molecular interactions or hydrogen bonds. In situ molecular imaging of ionic liquids and molecular modeling reveals that the acid gas capture mechanism by ionic liquids includes both physical adsorption and chemical bonding with multiple reaction pathways, engaging cluster formation and alteration of solvent structures.
We report a joint negative ion photoelectron spectroscopy (NIPES) and computational study on the electronic structures and noncovalent interactions of a series of cyclodextrin-closo-dodecaborate dianion complexes, ?-CD·B12X122- (? = a, ß, ?; X = H, F). The measured vertical / adiabatic detachment energies (VDEs / ADEs) are 1.15/0.93, 3.55/3.20, 3.90/3.60, and 3.85/3.60 eV for B12H122- and its a-, ß-, ?-CD complexes, respectively; while the corresponding values are 1.90/1.70, 4.00/3.60, 4.33/3.95, and 4.30/3.85 eV for the X = F case. These results show that the inclusion of B12X122- into the CD cavities greatly increase the electronic stability of the dianions. The effect of electronic stabilization for ß-CD is roughly the same as for ?-CD, both being considerably stronger than that for a-CD. Density functional theory (DFT) based geometry optimization reveals that B12X122- are inserted into CDs increasingly deeper from a-CD to ?-CD. The calculated VDEs and ADEs agree with the experiments well, particularly, reproducing the electron binding energy (EBE) trends. The molecular orbital analyses indicate that the most loosely bound photodetached electrons origin from the guest B12X122- moieties. In addition to a shift of all signals to larger EBE, significant changes in the signal patterns are observed. At low EBE, this is due to the splitting of highly degenerate B12X122- orbitals, while at high EBE, photodetachment from CD oxygens contributes to the new bands. The guest B12X122- and host CD nocovalent, size-specific interaction based on the independent gradient model (IGM) and energy decomposition analysis (EDA), is dominated by electrostatic interactions. The analysis further unravels unambiguiously the existence of dihydrogen bonding and how it affects the total energy that stabilizes the host-guest complexes of CDs·B12H122- compared to the general hydrogen bonding interaction in CDs·B12F122-. This work clearly exhibits strong influences on the electronic structures of dodecaborates upon clustering with CDs, with both size (a-, ß-, ?-) and molecular (X = H or F) specificities, thus providing critical molecular-level information on the cyclodextrin-closo-dodecaborate interactions of interest to medical applications, e.g. Boron neutron capature therapy. The NIPES experiments done at PNNL were supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Bioscience (X.-B.W.) and was performed at the EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. H.S. and Z.S. acknowledge the funding support of National Natural Science Foundation of China (Nos.11727810, 61720106009 and 21603074), the Science and Technology Commission of Shanghai Municipality (Nos. 19JC1412200), and the Program of Introducing Talents of Discipline to Universities 111 project (B12024). Z. L. thanks the China Scholarship Council (CSC) for financial support. We acknowledge the ECNU Multifunctional Platform for Innovation (001) and HPC Research Computing Team for providing computational and storage resources. J.W. acknowledges support from the Alexander von Humboldt foundation (Feodor Lynen Fellowship and Rückkehrerstipendium), a Freigeist fellowship of the Volkswagenfoundation and Prof. Vladimir A. Azov for helpful discussions.
Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.