Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic & molecular 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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Vibrational and Rovibrational Spectroscopy Applied to Astrochemistry

The detection of molecules in astrophysical environments almost always requires remote sensing. While radioastronomical observation and associated rotational spectroscopy are powerful astronomical tools, infrared spectral analysis provides a unique means of examining the observable universe, especially for molecules where permanent dipole moments are small or even non-existent. The molecular vibrations of small molecules are now able to be modeled via quantum chemistry and electronic structure theory conjoined to vibrational analysis to within spectroscopic accuracy in many cases. This chapter will showcase this success and build upon it to show how such advances are now being leveraged to describe molecular vibrations for molecules containing dozens of atoms, electronically excited states, "hot bands," exoplanetary atmospheric opacity data, and even emission of polycyclic aromatic hydrocarbons. All of these are required to prune the interstellar spectral garden of its "weeds" in search of "flowers" that will provide the necessary fingerprints for astronomers to be able to probe the heavens for its past, present, and future secrets.

Ryan C. Fortenberry↗
DOE OSTI2022

Thermal Neutron Scattering Cross Sections for Graphitic Amorphous Carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm 3 . From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical features and the optical graphitic peak at 0.25 eV for the ideal graphite DOS disappear for graphite-like amorphous carbon, which shows good agreement with the expected structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Zintl Concept Applied to Intergrowth Structures: Electron‐Hole Matching, Stacking Preferences, and Chemical Pressures in Pd 5 InAs

Abstract Enumerating the potential stacking sequences of layers is a fundamental way to account for the structure diversity of solid state compounds. In many cases, these stacking variations represent polymorphs with only small energetic differences. Here, we examine a compound for which the preferred stacking pattern instead reveals key aspects about its chemical bonding: Pd 5 InAs. Its structure is based on the intergrowth of slabs of the AuCu 3 and PtHg 2 (or alternatively, fluorite) structure types. Two basic stacking arrangements are available to this compound, represented by the Pd 5 TlAs and HoCoGa 5 structure types. DFT total energy calculations reveal that the former outcompetes the latter by a staggering 0.65 eV/formula unit. Through a combination of DFT‐reversed approximation Molecular Orbital (DFT‐raMO) and DFT‐Chemical Pressure (DFT‐CP) analysis we trace this preference to two factors. First, with DFT‐raMO analysis, we derive a Zintl‐like bonding scheme of Pd 5 InAs. This scheme, along with the inspection of selected crystal orbitals, is then connected to preferred stacking through the coordination environments of the Pd atoms at the interface between the Pd−In and Pd−As layers. In the hypothetical HoCoGa 5 ‐type and observed Pd 5 InAs‐type structures, different Pd coordination environments arise at the interfaces. The hypothetical structure features square planar PdIn 2 As 2 units, in each of which the same 4d‐orbital serves in the Pd sublattice's role as both Lewis acid (for interactions with the As) and Lewis base (for interactions with the In). In the observed structure, tetrahedral PdIn 2 As 2 units occur instead, so that these contradictory roles are distributed to separate 4d‐orbitals, leading to more effective bonding. DFT‐CP analysis illustrates that this driving force for the Pd 5 TlAs‐type arrangement is supplemented by a favorable alignment of the packing tensions in the parent structures. Altogether, the resulting picture demonstrates how the reaction of simple intermetallic structures to form intergrowths can be guided by recognizable chemical interactions.

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

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

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

Materials Data on VH17C4S4(O5F3)4 by Materials Project

VH12(S2O9)2(CF3)4H5O2 crystallizes in the triclinic P-1 space group. The structure is one-dimensional and consists of four fluoroform molecules; one molecular hydrogen;dihydrate molecule; and one VH12(S2O9)2 ribbon oriented in the (1, 0, 0) direction. In the VH12(S2O9)2 ribbon, V5+ is bonded in an octahedral geometry to six O2- atoms. There are a spread of V–O bond distances ranging from 2.02–2.04 Å. There are six inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a distorted single-bond geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.68 Å) H–O bond length. In the second H1+ site, H1+ is bonded in a single-bond geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.70 Å) H–O bond length. In the third H1+ site, H1+ is bonded in a distorted single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. In the fourth H1+ site, H1+ is bonded in a distorted single-bond geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.69 Å) H–O bond length. In the fifth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the sixth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. There are two inequivalent S+3.50+ sites. In the first S+3.50+ site, S+3.50+ is bonded in a trigonal non-coplanar geometry to three O2- atoms. There is two shorter (1.46 Å) and one longer (1.48 Å) S–O bond length. In the second S+3.50+ site, S+3.50+ is bonded in a trigonal non-coplanar geometry to three O2- atoms. There are a spread of S–O bond distances ranging from 1.45–1.48 Å. There are nine inequivalent O2- sites. In the first O2- site, O2- is bonded in a 1-coordinate geometry to one H1+ and one S+3.50+ atom. In the second O2- site, O2- is bonded in a distorted single-bond geometry to one S+3.50+ atom. In the third O2- site, O2- is bonded in a distorted bent 120 degrees geometry to one H1+ and one S+3.50+ atom. In the fourth O2- site, O2- is bonded in a distorted single-bond geometry to one S+3.50+ atom. In the fifth O2- site, O2- is bonded in a distorted single-bond geometry to one H1+ and one S+3.50+ atom. In the sixth O2- site, O2- is bonded in a 1-coordinate geometry to one S+3.50+ atom. In the seventh O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms. In the eighth O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms. In the ninth O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Americium Oxalate: An Experimental and Computational Investigation of Metal–Ligand Bonding

A novel actinide-containing coordination polymer, [Am(C 2 O 4 )(H 2 O) 3 Cl] (Am-1), has been synthesized and structurally characterized. The crystallographic analysis reveals that the structure is two-dimensional and comprised of pseudo-dimeric Am 3+ nodes that are bridged by oxalate ligands to form sheets. Each metal center is nine-coordinate, forming a distorted capped square antiprism geometry with a C 1 symmetry, and features bound oxalate, aqua, and chloro ligands. The Am 3+ –ligand bonds were probed computationally using the quantum theory of atoms in molecules nd natural localized molecular orbital approaches to investigate the underlying mechanisms and hybrid atomic orbital contributions therein. The analyses indicate that the bonds within Am-1 are predominantly ionic and the 5f shell of the Am 3+ metal centers does not add a significant covalent contribution to the bonds. Further, our bonding assessment is supported by measurements on the optical properties of Am-1 using diffuse reflectance and photoluminescence spectroscopies. The position of the principal absorption band at 507 nm ( 5 L 6' ← 7 F 0' ) is notable because it is consistent with previously reported americium oxalate complexes in solution, indicating similarities in the electronic structure and ionic bonding. Compound Am-1 is an active phosphor, featuring strong bright-blue oxalate-based luminescence with no evidence of metal-centered emission.

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

Materials Data on VH17C4S4(O5F3)4 by Materials Project

VH12(S2O9)2(CF3)4H5O2 crystallizes in the triclinic P-1 space group. The structure is one-dimensional and consists of four fluoroform molecules; one molecular hydrogen;dihydrate molecule; and one VH12(S2O9)2 ribbon oriented in the (1, 0, 0) direction. In the VH12(S2O9)2 ribbon, V5+ is bonded in an octahedral geometry to six O2- atoms. There are two shorter (2.03 Å) and four longer (2.04 Å) V–O bond lengths. There are six inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a distorted linear geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.61 Å) H–O bond length. In the second H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the third H1+ site, H1+ is bonded in a distorted single-bond geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.67 Å) H–O bond length. In the fourth H1+ site, H1+ is bonded in a distorted single-bond geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.69 Å) H–O bond length. In the fifth H1+ site, H1+ is bonded in a distorted linear geometry to two O2- atoms. There is one shorter (1.00 Å) and one longer (1.62 Å) H–O bond length. In the sixth H1+ site, H1+ is bonded in a distorted linear geometry to two O2- atoms. There is one shorter (1.01 Å) and one longer (1.58 Å) H–O bond length. There are two inequivalent S+3.50+ sites. In the first S+3.50+ site, S+3.50+ is bonded in a trigonal non-coplanar geometry to three O2- atoms. There is two shorter (1.46 Å) and one longer (1.47 Å) S–O bond length. In the second S+3.50+ site, S+3.50+ is bonded in a trigonal non-coplanar geometry to three O2- atoms. There is two shorter (1.46 Å) and one longer (1.47 Å) S–O bond length. There are nine inequivalent O2- sites. In the first O2- site, O2- is bonded in a 2-coordinate geometry to one H1+ and one S+3.50+ atom. In the second O2- site, O2- is bonded in a bent 120 degrees geometry to one H1+ and one S+3.50+ atom. In the third O2- site, O2- is bonded in a distorted bent 120 degrees geometry to one H1+ and one S+3.50+ atom. In the fourth O2- site, O2- is bonded in a 1-coordinate geometry to one H1+ and one S+3.50+ atom. In the fifth O2- site, O2- is bonded in a distorted single-bond geometry to one S+3.50+ atom. In the sixth O2- site, O2- is bonded in a distorted bent 120 degrees geometry to one H1+ and one S+3.50+ atom. In the seventh O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms. In the eighth O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms. In the ninth O2- site, O2- is bonded in a distorted trigonal planar geometry to one V5+ and two H1+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnesium ions mitigate metastable states in the regulatory landscape of mRNA elements

Residing in the 5' untranslated region of the mRNA, the 2'-deoxyguanosine (2'-dG) riboswitch mRNA element adopts an alternative structure upon binding of the 2'-dG molecule, which terminates transcription. RNA conformations are generally strongly affected by positively charged metal ions (especially Mg 2+ ). We have quantitatively explored the combined effect of ligand (2'-dG) and Mg 2+ binding on the energy landscape of the aptamer domain of the 2'-dG riboswitch with both explicit solvent all-atom molecular dynamics simulations (99 μsec aggregate sampling for the study) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments. We show that both ligand and Mg 2+ are required for the stabilization of the aptamer domain; however, the two factors act with different modalities. The addition of Mg 2+ remodels the energy landscape and reduces its frustration by the formation of additional contacts. In contrast, the binding of 2'-dG eliminates the metastable states by nucleating a compact core for the aptamer domain. Mg 2+ ions and ligand binding are required to stabilize the least stable helix, P1 (which needs to unfold to activate the transcription platform), and the riboswitch core formed by the backbone of the P2 and P3 helices. Mg 2+ and ligand also facilitate a more compact structure in the three-way junction region.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermodynamics and kinetics of interconversion between platinum nanoparticles and cations in zeolites

In metal-containing zeolites, sintering and redispersion processes exercise control over the identity of metal site structures, but the thermodynamic and kinetic factors that influence these molecular-level processes are not completely understood. Furthermore, we assess the ability of first principles informed free energy models (for supported and unsupported nanoparticles) and kinetic models integrating Ostwald ripening with atom trapping to describe the interconversion between Pt cations and nanoparticles encapsulated in chabazite (CHA) zeolites. Density functional theory-derived thermodynamic phase diagrams show that the interconversion between cations, favored in oxidizing environments, and particles, favored in reducing environments, is fully reversible within a wide range of their respective conditions (temperatures and pressures) for CHA and several other zeolite topologies. Kinetic Monte Carlo simulations of Pt redispersion are consistent with experimentally observed redispersion kinetics of encapsulated Pt nanoparticles in CHA zeolites, and model results suggest the zeolite host imparts additional stability for Pt nanoparticles. We envision our thermodynamic and kinetic models for Pt-CHA are also capable of describing nanoparticle and cation interconversions for other zeolite frameworks under similar conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In situ spectroscopic insights into the redox and acid-base properties of ceria catalysts

Cerium oxide (ceria) plays an important and fascinating role in heterogeneous catalysis as illustrated by its versatile use as a catalyst, a catalyst support, or a promotor in various oxidation and reduction reactions. Central to these reactions is the rich defect chemistry, facile redox capability, and unusual acid-base properties of ceria. Understanding the unique redox and acid-base properties of ceria is essential to build the structure-catalysis relationship so that improved catalytic functions can be achieved for ceria-based materials. Among the characterization toolbox, spectroscopic approach indisputably stands out for its unparalleled power in offering chemical insights into the surface properties of ceria at atomic and molecular level. In this review, we summarize advances in revealing the redox and acid-base properties of ceria via a variety of spectroscopic methods including optical, X-ray, neutron, electronic and nuclear spectroscopy. Both direct spectroscopy characterization and its coupling with probe molecules are analyzed to illustrate how the nature, strength and density of different surface sites are influenced by the pretreatment, the morphology and size of ceria nanoparticles. Further directions in taking advantage of in situ/operando spectroscopy for better understanding the catalysis of ceria-based materials are detailed in the summary and outlook section.

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

Influence of Polymorphs and Local Defect Structures on NMR Parameters of Graphite Fluorides

In this study, the role of local molecular structure on calculated 13 C and 19 F NMR chemical shifts for graphite fluoride materials was explored by using gauge-including projector augmented wave (GIPAW) computational methods for different periodic crystal polymorphs and density functional theory (DFT) gauge-including atomic orbital (GIAO) computational methods for individual graphite fluoride platelets, i.e., fluorinated graphene (FG). The impact of stacking sequences, d -spacing, and ring conformations on fully fluorinated graphite fluoride structures was investigated. A range of different defects including Stone–Wales, F and C vacancies, void formation, and F inversion were also evaluated using FG structures. These calculations show that distinct chemical shift signatures exist for many of these polymorphs and defects, therefore providing a basis for spectral assignment and development of models describing the mean local CF structure in disordered graphite fluoride materials.

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

Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABABR1a

A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABA B R1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABA B R1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = -1.91 ± 0.66 kcal mol -1 ), indicating a potentially enhanced GABA B R1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABA B R1a-specific therapeutics.

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

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cation Dynamics in Hybrid Halide Perovskites

Hybrid halide perovskite semiconductors exhibit complex, dynamical disorder while also harboring properties ideal for optoelectronic applications that include photovoltaics. However, these materials are structurally and compositionally distinct from traditional compound semiconductors composed of tetrahedrally coordinated elements with an average valence electron count of silicon. The additional dynamic degrees of freedom of hybrid halide perovskites underlie many of their potentially transformative physical properties. Neutron scattering and spectroscopy studies of the atomic dynamics of these materials have yielded significant insights into their functional properties. Specifically, inelastic neutron scattering has been used to elucidate the phonon band structure, and quasi-elastic neutron scattering has revealed the nature of the uncorrelated dynamics pertaining to molecular reorientations. Understanding the dynamics of these complex semiconductors has elucidated the temperature-dependent phase stability and origins of defect-tolerant electronic transport from the highly polarizable dielectric response. Furthermore, the dynamic degrees of freedom of the hybrid perovskites provide additional opportunities for application engineering and innovation.

Materials Science↗
DOE OSTI2022

First Principles Determination of the Potential-of-Zero-Charge in an Alumina-Coated Aluminum/Water Interface Model for Corrosion Applications

The surfaces of most metals immersed in aqueous electrolytes have a several-nanometer-thick oxide/hydroxide surface layer. This gives rise to the existence of both metal|oxide and oxide|liquid electrotlyte interfaces, and makes it challenging to correlate atomic length-scale structures with electrochemical properties such the potential-of-zero-charge (PZC). The PZC has been shown to be correlated the onset potential for pitting corrosion. Here, we conduct large-scale Density Functional Theory and ab initio molecular dynamics to calculate the PZC of a Al(111)| γ -Al 2 O 3 (110)| water double-interface model within the context of aluminum corrosion. By partitioning the multiple interfaces involved into binary components with additive contributions to the overall work function and voltage, we predict the PZC to be -1.53 V vs SHE for this model. Furthermore, we calculate the orbital energy levels of defects like oxygen vacancies in the oxide, which are critical parameters in theories associated with pitting corrosion. We predict that the Fermi level at the PZC lies above the impurity defect levels of the oxygen vacancies, which are therefore uncharged at the PZC. From the PZC estimate, we predict the voltage needed to create oxygen vacancies with net positive charges within a flatband approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗