Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cell behavior”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Behavior of Rapidly Sheared Bubble Suspensions

An experiment to be carried out aboard the International Space Station is described. A suspension consisting of millimeter-sized bubbles in water containing some dissolved salt, which prevents bubbles from coalescing, will be sheared in a Couette cylindrical cell. Rotation of the outer cylinder will produce centrifugal force which will tend to accumulate the bubbles near the inner wall. The shearing will enhance collisions among bubbles creating thereby bubble phase pressure that will resist the tendency of the bubbles to accumulate near the inner wall. The bubble volume fraction and velocity profiles will be measured and compared with the theoretical predictions. Ground-based research on measurement of bubble phase properties and flow in vertical channel are described.

Sangani, A. S.↗
NASA NTRS2017

Behavior of a Salt Water Solution Near the Critical Point Under Microgravity Conditions

The Supercritical Water Mixture (SCWM) experiment conducted on the International Space Station (ISS) is designed to study a salt-water solution in the vicinity of the critical transition with a test cell filled with Na2SO4 (aq) 0.5-w. This salt is a Type II salt which undergoes a dramatic decrease in solubility past waters critical point. The paper discusses the liquid-vapor-salt distribution in the test cell for subcritical, trans-critical, and supercritical regimes. In 0-g, the vapor-liquid interface manifests itself by the presence of a large single vapor bubble that is flattened between the two windows. In contrast with pure water, the vapor bubble remains centered in the test cell and is stable for temperatures up to Tc - 50 mK. The processes leading to critical transition during heat-up, as marked by the break-up of the vapor bubble, depends on the applied heating rate. For temperature ramps 1 mKmin, the bubble interface begins to breakup with the formation of a large number of smaller bubbles. For larger temperature ramps (10mKmin), temperature gradients arise in the cell pushing the bubble toward a window surface. As the critical threshold is crossed small variations in cell temperature ( 1 mK) cause an observable spatial variation in density. However, this spatial variation is not discernible as the temperature is increased further into the supercritical regime. The salt precipitation in the vicinity of the critical crossing appears to be heterogeneous in nature with precipitates observable on the cell windows. Direct imaging as well as light scattering indicate a hexagonal structure of the precipitated salt crystals in the near critical regime. The salt deposition is not fully reversible as the temperature is reduced due to window corrosion effects which traps some of the precipitated salt. Critical temperatures as high as 1.9 K above the value for pure water (647.25 K) have been measured. Finally, results on the critical transition phenomena in the presence of an imposed temperature gradient, demonstrating the remarkable stability of the central bubble, are also presented.

near-critical↗
DOE OSTI2022

Wellbore cement fracture permeability as a function of confining stress and pore pressure

Cemented annulus fractures are a major leakage path in a wellbore system, and their permeability plays an important role in the behavior of fluid flow through a leaky wellbore. The permeability of these fractures is affected by changing conditions including the external stresses acting on the fracture and the fluid pressure within the fracture. Laboratory gas flow experiments were conducted in a triaxial cell to evaluate the permeability of a wellbore cement fracture under a wide range of confining stress and pore pressure conditions. For the first time, an effective stress law that considers the simultaneous effect of confining stress and pore pressure was defined for the wellbore cement fracture permeability. Here the results showed that the effective stress coefficient (λ) for permeability increased linearly with the Terzaghi effective stress ( -p) with an average of λ = 1 in the range of applied pressures. The relationship between the effective stress and fracture permeability was examined using two physical-based models widely used for rock fractures. The results from the experimental work were incorporated into numerical simulations to estimate the impact of effective stress on the interpreted hydraulic aperture and leakage behavior through a fractured annular cement. Accounting for effective stress-dependent permeability through the wellbore length significantly increased the leakage rate at the wellhead compared with the assumption of a constant cemented annulus permeability.

58 GEOSCIENCES↗
DOE OSTI2020

Observation of an Unexpected n -Type Semiconducting Behavior in the New Ternary Zintl Phase Eu 3 InAs 3

The ternary arsenides Eu 3 InAs 3 and Sr 3 InAs 3 have been obtained by reactions of the elements in In flux at 1373 K. Structure elucidation by single-crystal X-ray diffraction reveals that Eu 3 InAs 3 and Sr 3 InAs 3 adopt the same orthorhombic structure (space group Pnma, Z = 4, Ca 3 AlAs 3 structure type) with unit cell parameters a = 12.9179(9) Å, b = 4.3990(3) Å, c = 13.9337(10) Å and a = 13.0218(11) Å, b = 4.4364(4) Å, c = 14.1339(12) Å, respectively. The structure consists of linear chains of corner-sharing InAs4 tetrahedra, [InAs 2 As 2/2 ] 6– , and Eu 2+ /Sr 2+ cations. Therefore, both Eu 3 InAs 3 and Sr 3 InAs 3 are valence-precise Zintl phases. As expected from the closed-shell electronic configurations, semiconducting behavior is confirmed by resistivity measurements on single crystals for both and by electronic band structure calculations for Sr 3 InAs 3 . The temperature dependence of resistivity and the computational work are in agreement that Eu 3 InAs 3 and Sr 3 InAs 3 are intrinsic semiconductors with narrow band gaps. Thermopower measurement on single-crystalline samples of Eu 3 InAs 3 shows that in the whole measured temperature range, from 300 to 700 K, the values for the Seebeck coefficient are negative. The observation of a negative Seebeck coefficient with very large absolute value (>400–500 μV K –1 at 700 K) is unexpected among the Zintl phases and suggestive that electrons are the majority charge carriers. Such a rare, n-type charge transport in an undoped compound such as Eu 3 InAs 3 , a material that has not been purposely optimized, could indicate native “defect” chemistry, and not extrinsic doping, as a reason for the unusual behavior. Here, a possible explanation involves a mixed-valent Eu 2+ /Eu 3+ state, which might be inferred from the measured effective paramagnetic moment of 7.2 μB per Eu atom, which is lower than the theoretically predicted value for free-ion moment of 7.9 μ B /Eu.

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

Chloride modification of vacuum-assisted blade-coated perovskite solar cells and mini-modules in ambient environment

Scalable deposition of perovskite thin films under ambient conditions remains a key challenge for the commercialization of perovskite photovoltaics, primarily due to the difficulty of controlling crystallization kinetics and suppressing non-radiative recombination during film formation. In this work, we present an effective strategy to modulate crystallization kinetics and non-radiative recombination in vacuum-quenched blade-coated perovskite films processed in air using chloride-based additives. By systematically comparing PbCl 2 and 3-chloropropylamine hydrochloride (Cl-PACl), we elucidate how different chloride incorporation pathways influence crystallization behavior and optoelectronic properties. In situ transmission and in situ photoluminescence measurements reveal that PbCl 2 acts as a transient additive that slows intermediate/perovskite phase formation while suppressing non-radiative recombination, leading to improved film quality. In contrast, Cl-PACl introduces additional trap states that deteriorate optoelectronic performance. As a result, PbCl 2 -modified perovskite solar cells achieve a champion power conversion efficiency (PCE) of 24.2% and remarkable operational stability, retaining over 90% of their initial performance after 800 hours of continuous illumination. Furthermore, we successfully scaled the process, fabricating a perovskite mini-module that delivered a champion efficiency of 18.0%. This study underscores the critical importance of coupling vacuum-assisted blade coating with transient additive engineering for achieving high-performance, stable, and scalable perovskite photovoltaics processed in ambient air.

14 SOLAR ENERGY↗
DOE OSTI2024

NbZr_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 niobium-zirconium (Nb-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 Nb 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 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. 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, 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 OSTI2024

Static and shock compression studies of eutectic high-entropy alloy AlCoCrFeNi 2.1 to ultrahigh pressures

The high-entropy alloy with composition AlCoCrFeNi 2.1 , additively manufactured with the laser powder-bed fusion technique, has a far-from-equilibrium BCC/FCC eutectic nanolamellar structure. We studied the high-pressure response of this alloy under both static compression and high-strain rate shock compression. The response to static compression using a diamond anvil cell was studied at pressures up to 302 GPa with synchrotron x-ray diffraction at the advanced photon source. The high-pressure FCC-only phase of the EHEA previously observed by Pope et al. [AIP Adv. 13, 035124 (2023)] is found to be stable up to the highest pressure achieved in this study with a volume compression of V/V 0 = 0.587 at ambient temperature. The shock experiments were performed by using GEKKO XII lasers at the Institute of Laser Engineering, Osaka University. The principal Hugoniot equation-of-state of the EHEA was measured up to a pressure of 515 GPa and a compression of V/V 0 = 0.613. Additionally, the thermal equation of state of the EHEA was measured up to 6.2 GPa and 1623 K using a large-volume Paris–Edinburgh cell to obtain the temperature dependence of bulk modulus and thermal expansion coefficients. The melting temperature for EHEA AlCoCrFeNi 2.1 at a pressure of 5.6 GPa was measured to be 1648 ± 25 K. These results can be used to refine stochastic (or special) quasi-random structure (SQS) models for high-pressure high-temperature behavior of high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog↗
DOE OSTI2022

Oxygen and Proton Transport in Flooded Graphene Pores with N-Dopants and Defects

Reactant transport is an important consideration in the design of ideal electrode structures. For the oxygen reduction reaction catalyzed by Pt/C in proton exchange membrane fuel cell cathodes, O 2 and H + delivery to Pt surfaces and H 2 O transport away are required. Some Pt nanoparticles may only be accessible via micropores that are too small for ionomer molecules to enter, possibly requiring flooding for H + transport. Here, to test if these “buried” Pt particles can play a role in activity through this proposed transport mechanism, we have performed atomic-scale simulations based on reactive force field molecular dynamics with an emphasis on confinement below 20 Å. Diffusion coefficients as a function of the molar concentration and local environment have been evaluated in water domains confined in two-dimensional graphene nanochannels of various channel heights representing a morphological model for micropores in proton exchange fuel cell cathodes. Our study shows that local atomic-scale structures can strongly modify H + , O 2 , and H 2 O transport rates in flooded micropores less than 20 Å in size. We find that there is a critical crossover in diffusion behavior around the 20 Å spacing with larger pores having bulk-like diffusion properties and confinement below 20 Å monotonically decreases diffusion rates. As pore size decreases, we observe locally dispersed water regions that ultimately strand reactants from long-distance transport. These findings suggest that flooded micropores may in fact be viable transport pathways for relevant reactants and products if the pore walls on opposite sides remain separated by ≥10 Å separation. Furthermore, the confinement effect is so strong that N-doping and C-vacancy defects in the C pore wall have only a minimal impact on diffusion rates and their effects are counterintuitively more apparent at larger spacings. These findings provide valuable insight regarding cathode performance and the role “stranded” catalyst particles may play in fuel cell cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2024

A Green, Fire–Retarding Ether Solvent for Sustainable High–Voltage Li–Ion Batteries at Standard Salt Concentration

Lithium-ion batteries (LIBs) are increasingly encouraged to enhance their environmental friendliness and safety while maintaining optimal energy density and cost-effectiveness. Although various electrolytes using greener and safer glyme solvents have been reported, the low charge voltage (usually lower than 4.0 V vs Li/Li + ) restricts the energy density of LIBs. Herein, tetraglyme, a less-toxic, non-volatile, and non-flammable ether solvent, is exploited to build safer and greener LIBs. It is demonstrated that ether electrolytes, at a standard salt concentration (1 m), can be reversibly cycled to 4.5 V vs Li/Li + . Anchored with Boron-rich cathode-electrolyte interphase (CEI) and mitigated current collector corrosion, the LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) cathode delivers competitive cyclability versus commercial carbonate electrolytes when charged to 4.5 V. Synchrotron spectroscopic and imaging analyses show that the tetraglyme electrolyte can sufficiently suppress the overcharge behavior associated with the high-voltage electrolyte decomposition, which is advantageous over previously reported glyme electrolytes. The new electrolyte also enables minimal transition metal dissolution and deposition. NMC811||hard carbon full cell delivers excellent cycling stability at C/3 with a high average Coulombic efficiency of 99.77%. In conclusion, this work reports an oxidation-resilient tetraglyme electrolyte with record-high 4.5 V stability and enlightens further applications of glyme solvents for sustainable LIBs by designing Boron-rich interphases.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A single oscillating proto-hypothalamic neuron gates taxis behavior in the primitive chordate Ciona

Ciona larvae display a number of behaviors, including negative phototaxis. In negative phototaxis, the larvae first perform short spontaneous rhythmic casting swims. As larvae cast in a light field, their photoreceptors are directionally shaded by an associated pigment cell, providing a phototactic cue. This then evokes an extended negative taxis swim. We report here that the larval forebrain of Ciona has a previously uncharacterized single slow-oscillating inhibitory neuron (neuron cor-assBVIN78) that projects to the midbrain, where it targets key interneurons of the phototaxis circuit known as the photoreceptor relay neurons. The anatomical location, gene expression and oscillation of cor-assBVIN78 suggest homology to oscillating neurons of the vertebrate hypothalamus. Ablation of cor-assBVIN78 results in larvae showing extended phototaxis-like swims, but which occur in the absence of phototactic cues. These results indicate that cor-assBVIN78 has a gating activity on phototaxis by projecting temporally-oscillating inhibition to the photoreceptor relay neurons. However, in intact larvae the frequency of cor-assBVIN78 oscillation does not match that of the rhythmic spontaneous swims, indicating that the troughs in oscillations do not themselves initiate swims, but rather that cor-assBVIN78 may modulate the phototaxis circuit by filtering out low level inputs while restricting them temporally to the troughs in inhibition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Tunable Orbital Ferromagnetism at Noninteger Filling of a Moiré Superlattice

Here, the flat bands resulting from moiré superlattices exhibit fascinating correlated electron phenomena such as correlated insulators, ( Nature 2018, 556 (7699), 80–84), ( Nature Physics 2019, 15 (3), 237) superconductivity, ( Nature 2018, 556 (7699), 43–50), ( Nature 2019, 572 (7768), 215–219) and orbital magnetism. ( Science 2019, 365 (6453), 605–608), ( Nature 2020, 579 (7797), 56–61), ( Science 2020, 367 (6480), 900–903) Such magnetism has been observed only at particular integer multiples of n 0 , the density corresponding to one electron per moiré superlattice unit cell. Here, we report the experimental observation of ferromagnetism at noninteger filling (NIF) of a flat Chern band in a ABC-TLG/hBN moiré superlattice. This state exhibits prominent ferromagnetic hysteresis behavior with large anomalous Hall resistivity in a broad region of densities centered in the valence miniband at n = -2.3n 0 . We observe that, not only the magnitude of the anomalous Hall signal, but also the sign of the hysteretic ferromagnetic response can be modulated by tuning the carrier density and displacement field. Rotating the sample in a fixed magnetic field demonstrates that the ferromagnetism is highly anisotropic and likely purely orbital in character.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Structural reversibility of Cu doped NU-1000 MOFs under hydrogenation conditions

The metal-organic framework (MOF), NU-1000, and its metalated counterparts have found proof-of-concept application in heterogeneous catalysis and hydrogen storage amongst others. A vapor-phase technique, akin to atomic layer deposition (ALD), is used to selectively deposit divalent Cu ions on oxo, hydroxo-bridged hexa-zirconium(IV) nodes capped with terminal –OH and -OH 2 ligands. Subsequent reaction with steam yields node-anchored, CuII-oxo,hydroxo clusters. We find that cluster installation via AIM (= ALD In MOFs) is accompanied by an expansion of MOF mesopore (channel) diameter . We investigated the behavior of the cluster-modified material, termed Cu-AIM-NU-1000, to heat treatment up to 325 °C, at atmospheric pressure with a low flow of H2 into the reaction cell. The response under these conditions revealed two important results: (1) Above 200 °C, the initially installed few-metal-ion clusters reduce to neutral Cu atoms. The neutral atoms migrate from the nodes and aggregate into Cu nanoparticles. While the size of particles formed in the MOF interior is constrained by the width of mesopores (ca. 3 nm), those formed on the exterior surface of the MOF can grow as large as ca. 8 nm. (2) Reduction and release of Cu atoms from the MOFs nodes is accompanied NU-1000 undergoes dynamic structural transformation as it reverts back to its original dimension following the release. These results show while the MOF framework itself remains intact at 325 °C in an H 2 atmosphere, the small, AIM-installed Cu II -oxo,hydroxo clusters are stable with respect to reduction and conversion to metallic nanoparticles only up to ~200 °C.

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

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Long term predictions of coal fired power plant data using evolved recurrent neural networks

Abstract This work presents an investigation into the ability of recurrent neural networks (RNNs) to provide long term predictions of time series data generated by coal fired power plants. While there are numerous studies which have used artificial neural networks (ANNs) to predict coal plant parameters, to the authors’ knowledge these have almost entirely been restricted to predicting values at the next time step, and not farther into the future. Using a novel neuro-evolution strategy called Evolutionary eXploration of Augmenting Memory Models (EXAMM), we evolved RNNs with advanced memory cells to predict per-minute plant parameters and per-hour boiler parameters up to 8 hours into the future. These data sets were challenging prediction tasks as they involve spiking behavior in the parameters being predicted. While the evolved RNNs were able to successfully predict the spikes in the hourly data they did not perform very well in accurately predicting their severity. The per-minute data proved even more challenging as medium range predictions miscalculated the beginning and ending of spikes, and longer range predictions reverted to long term trends and ignored the spikes entirely. We hope this initial study will motivate further study into this highly challenging prediction problem. The use of fuel properties data generated by a new Coal Tracker Optimization (CTO) program was also investigated and this work shows that their use improved predictive ability of the evolved RNNs.

Automation & Control Systems↗
NASA NTRS1991

Examination of the discharge mechanism of Li/CF(x) cells: Comparison of the electrochemical reduction mechanisms of PTFE and (CF(x))n by lithium

Swelling of (CF(x))(sub n) electrodes in commercial Li/(CF(x))(sub n) cells presents a limiting factor in cell design optimization. Examination of cathodes from such cells, after discharge, reveals a relation between cell operating temperatures and cathode swelling. Attempts to explain the swelling using the prevailing model for the cathode reaction have failed. A more suitable reaction mechanism is proposed based on the observed behavior of (CF(x))(sub n) electrodes on discharge and a comparison of the reaction products of (CF(x))(sub n) and polytetra fluoroethylene (PTFE) with lithium amalgams. The proposed mechanism is in agreement with the experimental data found in the literature.

Baxam, Carl C.↗
NASA NTRS2002

MAC/GMC 4.0 User's Manual: Keywords Manual

This document is the second volume in the three volume set of User's Manuals for the Micromechanics Analysis Code with Generalized Method of Cells Version 4.0 (MAC/GMC 4.0). Volume 1 is the Theory Manual, this document is the Keywords Manual, and Volume 3 is the Example Problem Manual. MAC/GMC 4.0 is a composite material and laminate analysis software program developed at the NASA Glenn Research Center. It is based on the generalized method of cells (GMC) micromechanics theory, which provides access to the local stress and strain fields in the composite material. This access grants GMC the ability to accommodate arbitrary local models for inelastic material behavior and various types of damage and failure analysis. MAC/GMC 4.0 has been built around GMC to provide the theory with a user-friendly framework, along with a library of local inelastic, damage, and failure models. Further, applications of simulated thermo-mechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated in MAC/GMC 4.0. Finally, classical lamination theory has been implemented within MAC/GMC 4.0 wherein GMC is used to model the composite material response of each ply. Consequently, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. This volume describes the basic information required to use the MAC/GMC 4.0 software, including a 'Getting Started' section, and an in-depth description of each of the 22 keywords used in the input file to control the execution of the code.

Bednarcyk, Brett A.↗