Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular cluster structure”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

In situ molecular imaging of ion clusters reveals the acid gas capture capacity and mechanism of water-lean ionic liquids

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.

Acid gas capture↗
DOE OSTI2020

Photoelectron spectroscopy and computational investigations of the electronic structures and noncovalent interactions of cyclodextrin-closo-dodecaborate anion complexes x-CD·B12X122- (x = a, ß, y; X = H, F)

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.

Li, Zhipeng↗
DOE OSTI2024

Large-scale atomistic model construction of subbituminous and bituminous coals for solvent extraction simulations with reactive molecular dynamics

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.

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

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Local-scale panmixia in the lichenized fungus Xanthoria parietina contrasts with substantial genetic structure in its Trebouxia photobionts

Microsatellite markers can provide valuable information about gene flow and population history. We developed and tested new microsatellites for the nitrophilic lichenized fungus Xanthoria parietina and studied its genetic diversity and structure within the urban area of Munich, Bavaria. We compared its local genetic pattern with that of its photobiont partner Trebouxia decolorans, for which existing microsatellites were applied. For comparison, a reference site with clean air was included in the sampling. We found support for three genetic clusters in the fungus X. parietina, which occurred intermingled in collecting sites. There was a high degree of admixture within fungal populations and individuals, and analysis of molecular variance revealed a lack of population structure in the mycobiont. The Trebouxia photobiont, in contrast, exhibited structured populations which grouped into two to five genetic clusters, and individuals showed less admixture than in the mycobiont. This indicates that the two lichen partners differ in their ability to move around in the landscape. The microsatellite markers we report are polymorphic and are suitable for population genetic studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Structural, Electronic, and Photophysical Insights into a Few Atom Copper-Sulfur Cluster in the Solid and Solution States

Coinage-metal chalcogenide clusters are widely studied for their attractive photoluminescence properties. Copper chalcogenides are especially promising, but are often confined to solid-state investigations due to their limited solution stability and the difficulty of synthesizing stable, well-defined clusters. Here, we investigate copper–sulfur clusters incorporating a small number of Cu atoms to elucidate fundamental atomic interactions, ground- and excited-state characteristics, and photophysical behavior in both solid and solution. We have synthesized the Cu6(4,6-dimethyl-2-mercaptopyrimidine)6 cluster in both neutral and charged states, Cu6 and Cu6-2H2+, respectively, by selective ligand protonation. The molecular structures are determined using single-crystal X-ray diffraction, while Cu K-edge X-ray absorption spectroscopy is used to probe Cu electronic structure differences arising from the ligand modification. Steady-state and pump-probe optical spectroscopy is used to investigate photophysical properties, interpreted using density functional theory methods. Both clusters exhibit good stability in the solid state and in solution and show characteristic near-infrared emission with microsecond lifetimes. Overall, the Cu6S6 clusters display favorable charge–transfer characteristics and show potential for further use in driving photochemical transformations.

Copper-sulfur clusters↗
DOE OSTI2022

Growth mechanism of glucose-based hydrochar under the effects of acid and temperature regulation

Improving the tailorability of hydrochar synthesis is an effective way to enhance its performance and utilization efficiency. In this study, the growth rate, morphology, and molecular structure of hydrochar were controlled by regulating the pH and temperature of the hydrothermal carbonization process. Growth process analysis indicates that hydrochar has three growth periods: induction, rapid, and stable growth periods. It is mainly controlled by 5-hydroxymethylfurfural (HMF), which is formed by converting glucose and its transformation products. Further, the regulation of acid can significantly shorten the induction period of hydrochar, even under low-temperature conditions (<180 °C), and increase the growth rate of hydrochar. However, the degree of hydrochar adhesion varies: the lower the temperature, the greater the degree of its adhesion. Molecular structural analysis demonstrates that hydrochar mainly consists of furan structural domains and aromatic clusters, and its surface is rich in oxygen-containing functional groups. The degree at which hydrochar was aromatized was improved by increasing the reaction temperature (160–220 °C); whereas the regulation of acid reduced it and increased the content of oxygen-containing functional groups on the hydrochar surface. Based on these results, it is proposed that hydrochar formation has five stages and three growth periods with or without acid regulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1996

Computational Nanotechnology at NASA Ames Research Center, 1996

Some forms of nanotechnology appear to have enormous potential to improve aerospace and computer systems; computational nanotechnology, the design and simulation of programmable molecular machines, is crucial to progress. NASA Ames Research Center has begun a computational nanotechnology program including in-house work, external research grants, and grants of supercomputer time. Four goals have been established: (1) Simulate a hypothetical programmable molecular machine replicating itself and building other products. (2) Develop molecular manufacturing CAD (computer aided design) software and use it to design molecular manufacturing systems and products of aerospace interest, including computer components. (3) Characterize nanotechnologically accessible materials of aerospace interest. Such materials may have excellent strength and thermal properties. (4) Collaborate with experimentalists. Current in-house activities include: (1) Development of NanoDesign, software to design and simulate a nanotechnology based on functionalized fullerenes. Early work focuses on gears. (2) A design for high density atomically precise memory. (3) Design of nanotechnology systems based on biology. (4) Characterization of diamonoid mechanosynthetic pathways. (5) Studies of the laplacian of the electronic charge density to understand molecular structure and reactivity. (6) Studies of entropic effects during self-assembly. Characterization of properties of matter for clusters up to sizes exhibiting bulk properties. In addition, the NAS (NASA Advanced Supercomputing) supercomputer division sponsored a workshop on computational molecular nanotechnology on March 4-5, 1996 held at NASA Ames Research Center. Finally, collaborations with Bill Goddard at CalTech, Ralph Merkle at Xerox Parc, Don Brenner at NCSU (North Carolina State University), Tom McKendree at Hughes, and Todd Wipke at UCSC are underway.

Globus, Al↗
DOE OSTI2023

Molecular diversity and evolution of far-red light-acclimated photosystem I

The need to acclimate to different environmental conditions is central to the evolution of cyanobacteria. Far-red light (FRL) photoacclimation, or FaRLiP, is an acclimation mechanism that enables certain cyanobacteria to use FRL to drive photosynthesis. During this process, a well-defined gene cluster is upregulated, resulting in changes to the photosystems that allow them to absorb FRL to perform photochemistry. Because FaRLiP is widespread, and because it exemplifies cyanobacterial adaptation mechanisms in nature, it is of interest to understand its molecular evolution. Here, we performed a phylogenetic analysis of the photosystem I subunits encoded in the FaRLiP gene cluster and analyzed the available structural data to predict ancestral characteristics of FRL-absorbing photosystem I. The analysis suggests that FRL-specific photosystem I subunits arose relatively late during the evolution of cyanobacteria when compared with some of the FRL-specific subunits of photosystem II, and that the order Nodosilineales, which include strains like Halomicronema hongdechloris and Synechococcus sp. PCC 7335, could have obtained FaRLiP via horizontal gene transfer. We show that the ancestral form of FRL-absorbing photosystem I contained three chlorophyll f-binding sites in the PsaB2 subunit, and a rotated chlorophyll a molecule in the A0B site of the electron transfer chain. Along with our previous study of photosystem II expressed during FaRLiP, these studies describe the molecular evolution of the photosystem complexes encoded by the FaRLiP gene cluster.

ancestral sequence reconstruction↗
DOE OSTI2022

Functionalization of Electrodes with Tunable [EMIM] x [Cl] x +1 – Ionic Liquid Clusters for Electrochemical Separations

Functionalization of electrodes with clusters of hydrophobic molecules may improve the energy efficiency and selectivity of electrochemical separations by modulating the desolvation process occurring at the interface. Ionic liquids (IL), which have a broad range of compositions and properties, are potential candidates for controlling the transport, desolvation, and adsorption of target ions at electrochemical interfaces. We report a joint experimental and theoretical investigation of the structure, stability, and selective adsorption properties of the IL clusters 1-ethyl-3-methylimidazolium chloride [EMIM] x [Cl] x+1 - (x = 1 – 10) and demonstrate their ability to adsorb and separate ions from solution. The structure and stability of the IL clusters are determined experimentally using high-mass-resolution electrospray ionization mass spectrometry, collision-induced dissociation, and negative ion photoelectron spectroscopy. Global optimization theory and ab initio molecular dynamics simulations provide molecular-level insight into the bonding and structural fluxionality of these species. Ion soft landing is used to selectively functionalize the surface of highly oriented pyrolytic graphite (HOPG) working electrodes with [EMIM] 1 [Cl] 2 - , [EMIM] 3 [Cl] 4 - , and [EMIM] 5 [Cl] 6 - clusters. Kelvin probe microscopy provides insight into the relative stability of the clusters on HOPG and their effect on the work function of IL-functionalized electrodes. Cyclic voltammetry measurements reveal irreversible adsorption of Fe(CN) 6 4-/3- anions during redox cycling, while electrochemical impedance spectroscopy indicates a substantial decrease in the electron transfer resistance of the IL-functionalized electrodes due to adsorption of Fe(CN) 6 4-/3- . Overall, our findings demonstrate that IL clusters with different size and stoichiometry may be used to increase the efficiency of electrochemical separations, opening new horizons in selective electrode functionalization.

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

Materials Data on Ti2H2CSe4Cl16O5 by Materials Project

Ti2C(OCl2)3(HSeOCl)2(SeCl3)2Cl2 crystallizes in the monoclinic P2/c space group. The structure is two-dimensional and consists of two molecular chlorine molecules; four HSeOCl clusters; four SeCl3 clusters; and one Ti2C(OCl2)3 sheet oriented in the (0, 1, 0) direction. In each HSeOCl cluster, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. Se3+ is bonded in a distorted L-shaped geometry to one O2- and one Cl1- atom. The Se–O bond length is 1.83 Å. The Se–Cl bond length is 2.50 Å. O2- is bonded in a distorted water-like geometry to one H1+ and one Se3+ atom. Cl1- is bonded in a distorted single-bond geometry to one Se3+ atom. In each SeCl3 cluster, Se3+ is bonded in a single-bond geometry to one Cl1- atom. The Se–Cl bond length is 2.76 Å. There are three inequivalent Cl1- sites. In the first Cl1- site, Cl1- is bonded in a 1-coordinate geometry to one Cl1- atom. The Cl–Cl bond length is 2.30 Å. In the second Cl1- site, Cl1- is bonded in a 2-coordinate geometry to one Se3+ and one Cl1- atom. The Cl–Cl bond length is 2.32 Å. In the third Cl1- site, Cl1- is bonded in a distorted linear geometry to two Cl1- atoms. In the Ti2C(OCl2)3 sheet, Ti4+ is bonded in a 5-coordinate geometry to one O2- and four Cl1- atoms. The Ti–O bond length is 2.71 Å. There are a spread of Ti–Cl bond distances ranging from 2.15–2.47 Å. C4+ is bonded in a linear geometry to two equivalent O2- atoms. Both C–O bond lengths are 1.18 Å. There are two inequivalent O2- sites. In the first O2- site, O2- is bonded in a single-bond geometry to one Ti4+ and one C4+ atom. In the second O2- site, O2- is bonded in a bent 150 degrees geometry to two equivalent Cl1- atoms. Both O–Cl bond lengths are 1.93 Å. There are three inequivalent Cl1- sites. In the first Cl1- site, Cl1- is bonded in a single-bond geometry to one Ti4+ atom. In the second Cl1- site, Cl1- is bonded in a water-like geometry to one Ti4+ and one O2- atom. In the third Cl1- site, Cl1- is bonded in an L-shaped geometry to two equivalent Ti4+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on MgFe4(S3O23)2 by Materials Project

MgO10(Fe2S3O16)2(O2)2 crystallizes in the triclinic P-1 space group. The structure is one-dimensional and consists of two molecular oxygen molecules; one MgO10 cluster; and two Fe2S3O16 ribbons oriented in the (-1, 1, 0) direction. In the MgO10 cluster, Mg is bonded in an octahedral geometry to six O atoms. There are a spread of Mg–O bond distances ranging from 1.97–2.48 Å. There are five inequivalent O sites. In the first O site, O is bonded in a distorted bent 150 degrees geometry to one Mg and one O atom. The O–O bond length is 1.23 Å. In the second O site, O is bonded in a distorted bent 150 degrees geometry to one Mg and one O atom. The O–O bond length is 1.23 Å. In the third O site, O is bonded in a single-bond geometry to one Mg atom. In the fourth O site, O is bonded in a single-bond geometry to one O atom. In the fifth O site, O is bonded in a single-bond geometry to one O atom. In each Fe2S3O16 ribbon, there are two inequivalent Fe sites. In the first Fe site, Fe is bonded to five O atoms to form FeO5 trigonal bipyramids that share corners with three SO4 tetrahedra. There are a spread of Fe–O bond distances ranging from 1.76–2.01 Å. In the second Fe site, Fe is bonded in a 6-coordinate geometry to six O atoms. There are a spread of Fe–O bond distances ranging from 1.93–2.24 Å. There are three inequivalent S sites. In the first S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.46–1.53 Å. In the second S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.45–1.53 Å. In the third S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.46–1.54 Å. There are sixteen inequivalent O sites. In the first O site, O is bonded in a bent 150 degrees geometry to one Fe and one S atom. In the second O site, O is bonded in a single-bond geometry to one S atom. In the third O site, O is bonded in a distorted bent 120 degrees geometry to one Fe and one S atom. In the fourth O site, O is bonded in a single-bond geometry to one S atom. In the fifth O site, O is bonded in a distorted bent 150 degrees geometry to one Fe and one S atom. In the sixth O site, O is bonded in a single-bond geometry to one S atom. In the seventh O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the eighth O site, O is bonded in a single-bond geometry to one S atom. In the ninth O site, O is bonded in a single-bond geometry to one S atom. In the tenth O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the eleventh O site, O is bonded in a single-bond geometry to one S atom. In the twelfth O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the thirteenth O site, O is bonded in a bent 120 degrees geometry to two Fe atoms. In the fourteenth O site, O is bonded in a single-bond geometry to one Fe atom. In the fifteenth O site, O is bonded in a distorted L-shaped geometry to one Fe and one O atom. The O–O bond length is 1.26 Å. In the sixteenth O site, O is bonded in a 2-coordinate geometry to one Fe and one O atom.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

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

Is the Trotterized UCCSD Ansatz Chemically Well-Defined?

The variational quantum eigensolver (VQE) has emerged as one of the most enticing near-term quantum algorithms that can be used to simulate many-body systems such as molecular electronic structures. Serving as an attractive ansatz in the VQE algorithm, unitary coupled cluster (UCC) theory has seen a renewed interest in recent literature. Yet, unlike the original classical UCC theory, implementation on a quantum computer requires a finite-order Suzuki-Trotter decomposition to separate the exponentials of the large sum of Pauli operators. While previous literature has recognized the non-uniqueness of different orderings of the operators in the Trotterized form of UCC methods, the question of whether or not different orderings matter at the chemical scale has not been addressed. In this letter, we explore the effect of operator ordering on the Trotterized UCCSD ansatz, as well as the much more compact k-UpCCGSD ansatz recently proposed by Lee et al. We observe a significant, system-dependent variation in the energies of Trotterizations with different operator orderings. The energy variations occur on a chemical scale, sometimes on the order of hundreds of kcal/mol. This letter establishes the need to define not only the operators present in the ansatz, but also the order in which they appear. This is necessary for adhering to the quantum chemical notion of a ``model chemistry'', in addition to the general importance of scientific reproducibility. As a final note, we suggest a useful strategy to select out of the combinatorial number of possibilities, a single well-defined and effective ordering of the operators.

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

Observations of grain-boundary phase transformations in an elemental metal

The theory of grain boundary (the interface between crystallites, GB) structure has a long history 1 and the concept of GBs undergoing phase transformations was proposed 50 years ago 2,3 . The underlying assumption was that multiple stable and metastable states exist for different GB orientations 4,5,6 . The terminology ‘complexion’ was recently proposed to distinguish between interfacial states that differ in any equilibrium thermodynamic property 7 . Different types of complexion and transitions between complexions have been characterized, mostly in binary or multicomponent systems 8,9,10,11,12,13,14,15,16,17,18,19 . Simulations have provided insight into the phase behaviour of interfaces and shown that GB transitions can occur in many material systems 20,21,22,23,24 . However, the direct experimental observation and transformation kinetics of GBs in an elemental metal have remained elusive. In this work we demonstrate atomic-scale GB phase coexistence and transformations at symmetric and asymmetric $[11\bar{1}]$ tilt GBs in elemental copper. Atomic-resolution imaging reveals the coexistence of two different structures at Σ19b GBs (where Σ19 is the density of coincident sites and b is a GB variant), in agreement with evolutionary GB structure search and clustering analysis 21,25,26 . We also use finite-temperature molecular dynamics simulations to explore the coexistence and transformation kinetics of these GB phases. Our results demonstrate how GB phases can be kinetically trapped, enabling atomic-scale room-temperature observations. Our work paves the way for atomic-scale in situ studies of metallic GB phase transformations, which were previously detected only indirectly 9,15,27,28,29 , through their influence on abnormal grain growth, non-Arrhenius-type diffusion or liquid metal embrittlement.

36 MATERIALS SCIENCE↗