Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference cells”

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 883 records · Page 49

ASTM E 1559 method for measuring material outgassing/deposition kinetics has applications to aerospace, electronics, and semiconductor industries

The American Society for Testing and Materials has published a new standard test method for characterizing time and temperature-dependence of material outgassing kinetics and the deposition kinetics of outgassed species on surfaces at various temperatures. This new ASTM standard, E 1559(1), uses the quartz crystal microbalance (QCM) collection measurement approach. The test method was originally developed under a program sponsored by the United States Air Force Materials Laboratory (AFML) to create a standard test method for obtaining outgassing and deposition kinetics data for spacecraft materials. Standardization by ASTM recognizes that the method has applications beyond aerospace. In particular, the method will provide data of use to the electronics, semiconductor, and high vacuum industries. In ASTM E 1559 the material sample is held in vacuum in a temperature-controlled effusion cell, while its outgassing flux impinges on several QCM's which view the orifice of the effusion cell. Sample isothermal total mass loss (TML) is measured as a function of time from the mass collected on one of the QCM's which is cooled by liquid nitrogen, and the view factor from this QCM to the cell. The amount of outgassed volatile condensable material (VCM) on surfaces at higher temperatures is measured as a function of time during the isothermal outgassing test by controlling the temperatures of the remaining QCM's to selected values. The VCM on surfaces at temperatures in between those of the collector QCM's is determined at the end of the isothermal test by heating the QCM's at a controlled rate and measuring the mass loss from the end of the QCM's as a function of time and temperature. This reevaporation of the deposit collected on the QCM's is referred to as QCM thermogravimetric analysis. Isothermal outgassing and deposition rates can be determined by differentiating the isothermal TML and VCM data, respectively, while the evaporation rates of the species can be obtained as a function of temperature by differentiating the QCM thermogravimetric analysis data.

Garrett, J. W.↗

Hypertonicity-induced transmitter release at Drosophila neuromuscular junctions is partly mediated by integrins and cAMP/protein kinase A

The frequency of quantal transmitter release increases upon application of hypertonic solutions. This effect bypasses the Ca(2+) triggering step, but requires the presence of key molecules involved in vesicle fusion, and hence could be a useful tool for dissecting the molecular process of vesicle fusion. We have examined the hypertonicity response at neuromuscular junctions of Drosophila embryos in Ca(2+)-free saline. Relative to wild-type, the response induced by puff application of hypertonic solution was enhanced in a mutant, dunce, in which the cAMP level is elevated, or in wild-type embryos treated with forskolin, an activator of adenylyl cyclase, while protein kinase A (PKA) inhibitors decreased it. The response was also smaller in a mutant, DC0, which lacks the major subunit of PKA. Thus the cAMP/PKA cascade is involved in the hypertonicity response. Peptides containing the sequence Arg-Gly-Asp (RGD), which inhibit binding of integrins to natural ligands, reduced the response, whereas a peptide containing the non-binding sequence Arg-Gly-Glu (RGE) did not. A reduced response persisted in a mutant, myospheroid, which expresses no integrins, and the response in DC0 was unaffected by RGD peptides. These data indicate that there are at lease two components in the hypertonicity response: one that is integrin mediated and involves the cAMP/PKA cascade, and another that is not integrin mediated and does not involve the cAMP/PKA cascade.

NASA Discipline Cell Biology↗

Integrated Cryogenic Experiment (ICE) microsphere investigation

The main objective is to determine the performance of microsphere insulation in a 0-g environment and compare its performance to reference insulations such as multilayer insulation. The Lockheed Helium Extended-Life Dewar (HELD) is used to provide superfluid-helium cold sink for the experiment. The use of HELD allows the low-g dynamic properties of Passive Orbital Disconnect Struts (PODS) to be characterized and provides a flight demonstration of the PODS system. The thermal performance of microspheres in 1 and 0 g was predicted, a flight experiment was designed to determine microsphere thermal performance, and the interface was also designed between the experimental package and the shuttle through HELD and the Hitchhiker-M carrier. A single test cell was designed and fabricated. The cell was filled with uncoated glass microspheres and tested with a liquid-nitrogen cold sink. The data were found to agree with predictions of microsphere performance in 1 g.

Spradley, I.↗

Assessment of Low-Level Pu Isotope Ratio Measurements Using Multicollector Inductively Coupled Plasma Mass Spectrometry (MC-ICP-MS/MS) Equipped with a Pre-Mass Filter

We present an initial investigation into the performance of a multicollector inductively coupled plasmamass spectrometer equipped with a pre-mass filter (Neoma MC-ICP-MS/MS) for making plutonium (Pu) isotope ratio measurements on solutions containing low level (i.e., pg mL –1 ) Pu concentrations. This assessment was achieved by comparison of the 240 Pu/ 239 Pu, 241 Pu/ 239 Pu, and 242 Pu/ 239 Pu ratios attained over a one month period on the MC-ICP-MS/MS with the long-term (∼1 year) performance observed on the predecessor MC-ICP-MS (Neptune Plus) instrument each equipped with an equipped with an APEXΩ desolvating nebulizer for repeated measurements of certified reference materials from New Brunswick Program Office (NBL PO) CRM 136a and CRM 137. The MC-ICP-MS/MS performance of repeated measurement of CRM 136a (n = 20) resulted in mean values of 240 Pu/ 239 Pu = 0.1448 ± 0.0006, 241 Pu/ 239 Pu = 0.00371 ± 0.00006, and 242 Pu/ 239 Pu = 0.00682 ± 0.00006 (k = 2). The CRM 137 (n = 20), analyzed during the same analytical sessions, produced mean values for 240 Pu/ 239 Pu = 0.2414 ± 0.0006, 241 Pu/ 239 Pu = 0.00464 ± 0.00007, and 242 Pu/ 239 Pu = 0.0157 ± 0.0001 (k = 2). These results closely align with the certificate values for CRM 136a and CRM 137 and are within the k = 2 envelopes defined by the long-term performance of the traditional MC-ICP-MS approach (Neptune Plus). Examination of the performance of the various Pu isotope ratios as a function of total Pu content revealed accurate results (<3% relative difference, or RD) above ∼50 fg total Pu. The results presented here demonstrate the capability of the MC-ICP-MS/MS making accurate and precise low level Pu isotopic measurements. While the intent of this work was not to investigate the functionality of the collision cell, the pre-mass filter was employed. Future studies are warranted to investigate the entire capability of the MC-ICP-MS/MS collision cell and pre-cell mass filter optimization for performing low level Pu isotope measurements, even in mixed matrix samples.

CRM↗

The importance of cycle-by-cycle data in performing rapid battery technology development and validation

Lithium-ion battery (LiB) technology is playing a crucial role in transforming the predominantly fossil fuel-based transportation and stationary storage sectors to achieve a low-carbon economy. Rapid innovation in the LiB materials to electrode to cell design is happening to satisfy the performance, life, and safety metrics required by those myriads of applications. Lately, advanced analytics, such as machine-learning or artificial intelligence (ML/AI) techniques, are being used more frequently to aid in expedited LiB technology development, performance validation, and life prediction. The success of these techniques often relies on a large volume of well-defined and high-quality battery test data. On the other hand, most battery developers and research and development (R&D) communities are still following a classical approach to develop batteries, which is running calendar- and/or cycle-aging tests, performing reference performance tests (RPTs), and conducting post-mortem analyses periodically without paying attention to the wealth of data often not collected during the calendar or cycle life aging tests. This sparse data collection approach is time- and resource-intensive, requiring data capture and evaluation of months to years of RPT data to diagnose accurate battery state of performance, health, and safety. Even so, the underlying aging modes and mechanisms can be missed. If collected properly, battery test data during cycling or calendaring can be efficiently combined with ML/AI techniques to create powerful tools in the rapid diagnosis of battery state of performance, health, and safety along with insights into underlying aging modes and mechanisms. In this report, we discuss the importance of effective cycle-by-cycle (CBC) data collection with example case studies. Within a reasonable timeframe, RPT data are often inadequate in capturing many of the crucial battery aging dynamics, which often predominantly show up in CBC test data. Finally, we also show examples of ML/AI techniques that use CBC data in rapid diagnosis and projection of LiB state of health (SOH) to motivate the scientific community in collecting and using CBC data to facilitate expeditious technology development and validation.

25 ENERGY STORAGE↗

MPAS - Ocean Simulation Quality for Variable-Resolution North American Coastal Meshes

Climate model components utilizing unstructured meshes enable variableresolution, regionally enhanced simulations within global domains. Here we investigate the relationship between mesh quality and simulation statistics using the JIGSAW unstructured meshing library and the Model for Prediction Across ScalesOcean (MPASOcean) with a focus on Gulf Stream dynamics. In the base configuration, the refined region employs 8 km cells that extend 400 km from the coast of North America. This coastal refined region is embedded within a lowresolution global domain, with cell size varying latitudinally between 30 and 60 km. The resolution transition region between the refined region and background mesh is 600 km wide. Three sensitivity tests are conducted: 1) the quality of meshes is intentionally degraded so that horizontal cells are progressively more distorted; 2) the transition region from high to low resolution is steepened; and 3) resolution of the coastal refinement region is varied from 30 km to 8 km. Overall, the ocean simulations are shown to be robust to mesh resolution and quality alterations. Meshes that are substantially degraded still produce realistic currents, with Southern Ocean transports within 0.4% and Gulf Stream transports within 12% of highquality mesh results. The narrowest transition case of 100 km did not produce any spurious effects. Refined regions with high resolution produce eddy kinetic energy and sea surface height variability that are similar to the highresolution reference simulation. These results provide heuristics for the design criteria of variableresolution climate model domains.

Hoch, Kristin E.↗

Physicochemical Heterogeneity in Silicon Anodes from Cycled Lithium-Ion Cells

The severe capacity fade of lithium-ion cells with silicon-dominant anodes has hindered their widescale commercialization. In this work, we link cell capacity fade to the heterogeneous physicochemical evolution of silicon anodes during battery cycling. Through a multi-length scale characterization approach, we demonstrate that silicon particles near the anode surface react differently than those near the copper current collector. In particular, near the anode surface we find an amorphized wispy silicon encased in a highly fluorinated matrix of electrolyte-reduction products. In contrast, closer to the current collector, the silicon retains more of its initial morphology and structure, suggesting the presence of isolated particles. Here, the results show that accessibility of active silicon to lithium-ions varies across the anode matrix. Material and cell designs, which minimize electrode expansion resulting from the in-filling of pores with the solid electrolyte interphase (SEI), are needed to enhance anode homogeneity during the electrochemical cycling.

25 ENERGY STORAGE↗

The August Krogh principle applies to plants

The Krogh principle refers to the use of a large number of animals to study the large number of physiological problems, rather than limiting study to a particular organism for all problems. There may be organisms that are more suited to study of a particular problem than others. This same principle applies to plants. The authors are concerned with the recent trend in plant biology of using Arabidopsis thaliana as the "organism of choice." Arabidopsis is an excellent organism for molecular genetic research, but other plants are superior models for other research areas of plant biology. The authors present examples of the successful use of the Krogh principle in plant cell biology research, emphasizing the particular characteristics of the selected research organisms that make them the appropriate choice.

NASA Discipline Plant Biology↗

Thin Cation-Exchange Layers Enable High-Current-Density Bipolar Membrane Electrolyzers via Improved Water Transport

With suitable water dissociation (WD) catalysts, bipolar membranes (BPMs) can efficiently dissociate water into H + and OH – at the junction between anion- and cation-exchange layers (AEL and CEL, respectively). First, however, water must be transported through the AEL or CEL and thus against the outward flow of hydrated H + and OH – . This is a challenge intrinsic to the BPM architecture and limits operation to current densities typically less than ~0.5 A·cm –2 . Here we explore how water transport affects durability and performance in reference alkaline and acidic membrane electrolyzers, and we use the insight gained to design BPMs with improved water transport. We demonstrate a thin-CEL BPM (2-μm Nafion CEL|~200 nm TiO 2 |~200 nm NiO + ionomer|50 μm Sustainion AEL) which maintains a pH difference of ~14 units between the anode and cathode for current densities of up to 3.4 A·cm –2 with a total water electrolysis voltage of ~4 V and an estimated WD overpotential of ~1.5 V. Finally, such high-current-density operation is crucial for key emerging BPM applications, including in water and carbon-dioxide electrolyzers and in (regenerative) fuel cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Creating Benchmark Data for Artificial Intelligence and Machine Learning Space Biology Research

To identify an appropriate AI/ML approach for a specific problem, the best practice is to measure algorithm performance through the benchmarking process. A scientific benchmark consists of an AI-ready dataset and a reference implementation on a specific scientific question. The NASA Science Mission Directorate (SMD) has started the “Benchmark Initiative for AI/ML to create scientific benchmark datasets in three applications: 1) scientific benchmarking, which finds the best algorithm for a specific problem; 2) application benchmarking, which measures algorithm performance against a set of parameters; and 3) system benchmarking, which evaluates performance of hardware and software architecture. Currently, there are no standardized datasets available to benchmark AI/ML algorithms in the domain of space biology. In this work, we constructed two AI/ML-ready biological datasets from experiments in space-flown mice: cellular imaging and RNA-seq. First, radiation-exposed immune cells harbor DNA damage foci that can be fluorescently marked to visualize the amount of damage following exposure to ionizing radiation. However, such large datasets are difficult to analyze visually, due to imaging inconsistencies and human bias, and classical image processing approaches can fail on imaging artifacts. AI/ML are therefore exciting alternative, providing the speed of machines and the accuracy of humans. We have made this dataset available at https://registry.opendata.aws/bps_microscopy/. Second, high-throughput nucleic acid sequencing (DNA-seq, RNA-seq) has become widespread in biomedical research due to the growing availability and affordability of these assays. However, most sequencing datasets suffer from high dimensionality and low sample count. In this work, we used a generative adversarial network to synthesize a standardized, AI-ready, publicly available benchmark dataset for space biology RNA-seq data with sufficient space-flown and ground control mouse liver samples from NASA GeneLab. This dataset is available at https://registry.opendata.aws/bps_rnaseq/. These datasets are now fully open the Space Biology community to test their favorite AI/ML approaches.

James Casaletto↗

Theoretical Multinary Perovskite Oxides Dataset

Perovskite oxides (ternary chemical formula ABO3) are a diverse class of materials with applications including heterogeneous catalysis, solid-oxide fuel cells, thermochemical conversion, and oxygen transport membranes. However, their multicomponent chemical space is underexplored due to the immense number of possible compositions. To expand the number of computed compounds we report a dataset of 66,516 theoretical multinary oxides, 59,708 of which are perovskites. First, the AA'BB'O6 compositions were generated in the a-b+a- Glazer tilting mode using the computationally-inexpensive Structure Prediction and Diagnostic Software (SPuDS) program. Next, we optimized these structures with DFT using parameters compatible with the Materials Project (MP) database. Our dataset contains these optimized structures and their formation (ΔHf) and decomposition enthalpies (ΔHd) computed relative to MP tabulated elemental references and competing phases, respectively. This dataset can be mined, used to train machine learning models, and rapidly and systematically expanded by optimizing more SPuDS-generated AA'BB'O6 perovskite structures using MP-compatible DFT calculations

Bare, Zachary↗

Apparatus for Testing Flat Specimens of Thermal Insulation

An apparatus has been developed to implement an improved method of testing flat-plate specimens of thermal-insulation materials for cryogenic application. The method includes testing under realistic use conditions that could include vacuum and mechanical loading at a pressure up to 70 psi (=0.48 MPa). The apparatus can accommodate a rigid or flexible specimen having thickness up to 1.25 in. (=3.2 cm) and diameters between 6 and 10 in. (about 15.2 and 25.4 cm, respectively). Typical test conditions include boundary temperatures between 77 K and 373 K and vacuum/interstitial gas filling at a pressure between 10(exp -6) torr (=1.3 x 10(exp -4) Pa) and 760 torr (atmospheric pressure =0.1 MPa). The interstitial gas could be N2, He, CO2, or any other suitable gas to which the insulation is expected to be exposed in use. Relative to prior apparatuses and testing methods, this apparatus and the testing method that it implements offer advantages of relative simplicity and ease of use. The basic principle of operation of the apparatus is that of boil-off calorimetry, using liquid nitrogen or any other suitable liquid that boils at a desired temperature below ambient temperature. Comparative rates of flow of heat through the thicknesses of the specimens (heat-leak rates) and apparent-thermal-conductivity values are obtained from tests of specimens. Absolute values of heat-leak rates and apparent thermal conductivities are computed from a combination of (1) the aforementioned comparative values and (2) calibration factors obtained by testing reference specimens of materials that have known thermal-insulation properties. The apparatus includes a full complement of temperature sensors, a vacuum pump and chamber, a monitoring and control system, and tools and fixtures that enable rapid and reliable installation and removal of specimens. A specimen is installed at the bottom of the vacuum chamber, and a cold-mass assembly that includes a tank is lowered into position above and around the specimen (see figure). A spring-based compensating fixture helps to ensure adequate thermal contact with possibly irregular specimen surfaces. For a high-compression test, the springs can be replaced with spacers. A flat circular load cell at the bottom of the chamber measures the compressive load on the specimen. Once the desired compressive-load, temperature, and vacuum/gas-filling conditions are established, testing begins. During a test, all measurements are recorded by use of a portable data-acquisition system and a computer. The total heat-leak rate is measured and calculated as the boil-off flow rate multiplied by the latent heat of vaporization. The parasitic heat leak (to the side of the specimen and to the top and side of the cold-mass tank) is reduced to a small fraction of the total heat leak by use of a combination of multilayer-insulation (MLI) shield rings, reflective film, a fiberglass/epoxy centering ring, and a bulk fill of aerogel beads. This combination eliminates the need for a cryogenic guard chamber used in a typical prior apparatus to reduce the parasitic heat leak.

Fesmire, James E.↗

Examples of X-Ray Characterization Techniques in Energy Storage Research

Lithium-ion batteries have revolutionized the portable electronics and transportation sectors. Their performance is often critically dependent on the crystal structures of the anode and cathode electrode materials, which must enable the transport and reversible storage of lithium ions into and out of the lattice. Because lithium is a low-Z element, characterization of materials for lithium-ion batteries can be particularly challenging. Regardless, X-ray techniques enable analysis of material structures to better understand how battery materials perform and degrade, particularly when combined with other materials characterization and electrochemical characterization techniques. While X-ray techniques are most often used in battery research for phase identification of crystal structures, X-ray characterization techniques are also used for a wide variety of other purposes. I will discuss several examples from my research with various collaborators on several projects that highlight the impact that X-ray characterization techniques can have on battery research. The first example will focus on low-temperature microwave-assisted solvothermal synthesis of vanadium-doped LiFePO4 cathode materials for lithium-ion batteries. (1,2) Through a combination of electrochemical and materials characterization, we determined that low temperature synthesis resulted in metastable phases that enabled incorporation of higher dopant levels than resulting from high-temperature synthesis of thermodynamically stable phases. Rietveld refinement of X-ray diffraction data enabled understanding of how lattice parameters changed with doping levels and synthesis temperature. X-ray absorption near edge spectroscopy enabled understanding of the vanadium and iron oxidation states to confirm how vacancies in the structure caused by doping were charge compensated. This was important to understand because the literature suggests doping can improve LiFePO4 electrical conductivity, which improves battery charge and discharge rates. The second example will focus on understanding residual strain in lithium metal anodes. Lithium-ion batteries typically use graphite anodes, but the charge-storage capacity can be theoretically improved ~10x by using lithium metal as the anode material instead. However, lithium anodes suffer from growth of high-aspect-ratio features, such as dendrites, that can pierce nanoporous polymer separators and lead to short circuits and fires. External pressure is commonly applied to cells to enable better morphological control. We hypothesized that applied pressure may promote strain and possibly work hardening during electrochemical cycling, which motivated us to look for evidence of residual strain in lithium metal cycled under applied pressure using X-ray diffraction and sin2(..psi..) analysis. We found that lithium electrodeposited under high pressure exhibited in-plane compressive strain and that that lithium electrodeposited under low pressure did not. (3) The residual strain that accompanies electrodeposition under high pressure may lead to work hardening, which may explain how a soft metal like lithium can puncture separators and why higher pressure does not always decrease short circuits. (4-6) References: 1) Harrison, K. L.; Manthiram, A. Microwave-Assisted Solvothermal Synthesis and Characterization of Metastable LiFe1- x (VO) x PO4 Cathodes. Inorganic chemistry 2011, 50(8), 3613-3620. 2) Harrison, K. L.; Bridges, C. A.; Paranthaman, M. P.; Segre, C. U.; Katsoudas, J.; Maroni, V. A.; Idrobo, J. C.; Goodenough, J. B.; Manthiram, A. Temperature Dependence of Aliovalent-Vanadium Doping in LiFePO4 Cathodes. Chemistry of Materials 2013, 25(5), 768-781. 3) Rodriguez, M. A.; Harrison, K. L.; Goriparti, S.; Griego, J. J.; Boyce, B. L.; Perdue, B. R. Use of a Be-Dome Holder for Texture and Strain Characterization of Li Metal Thin Films via Sin2 (..psi..) Methodology. Powder Diffraction 2020, 35(2), 89-97. 4) Jungjohann, K. L.; Gannon, R. N.; Goriparti, S.; Randolph, S. J.; Merrill, L. C.; Johnson, D. C.; Zavadil, K. R.; Harris, S. J.; Harrison, K. L. Cryogenic Laser Ablation Reveals Short-Circuit Mechanism in Lithium Metal Batteries. ACS Energy Letters 2021, 6(6), 2138-2144. 5) Harrison, K. L.; Merrill, L. C.; Long, D. M.; Randolph, S. J.; Goriparti, S.; Christian, J.; Warren, B.; Roberts, S. A.; Harris, S. J.; Perry, D. L. Cryogenic Electron Microscopy Reveals That Applied Pressure Promotes Short Circuits in Li Batteries. Iscience 2021, 24(12). 6) Harrison, K. L.; Goriparti, S.; Merrill, L. C.; Long, D. M.; Warren, B.; Roberts, S. A.; Perdue, B. R.; Casias, Z.; Cuillier, P.; Boyce, B. L. Effects of Applied Interfacial Pressure on Li-Metal Cycling Performance and Morphology in 4 M LiFSI in DME. ACS Applied Materials & Interfaces 2021, 13(27), 31668-31679.

batteries↗

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↗

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↗

Single-crystal elasticity of (Al,Fe)-bearing bridgmanite up to 82 GPa

Abstract Thermoelastic properties of mantle candidate minerals are essential to our understanding of geophysical phenomena, geochemistry, and geodynamic evolutions of the silicate Earth. However, the lower-mantle mineralogy remains much debated due to the lack of single-crystal elastic moduli (Cij) and aggregate sound velocities of (Al,Fe)-bearing bridgmanite, the most abundant mineral of the planet, at the lower mantle pressure-temperature (P-T) conditions. Here we report single-crystal Cij of (Al,Fe)-bearing bridgmanite, Mg0.88Fe0.1Al0.14Si0.90O3 (Fe10-Al14-Bgm) with Fe3+/ΣFe = ~0.65, up to ~82 GPa using X-ray diffraction (XRD), Brillouin light scattering (BLS), and impulsive stimulated light scattering (ISLS) measurements in diamond-anvil cells (DACs). Two crystal platelets with orientations of (–0.50, 0.05, –0.86) and (0.65, –0.59, 0.48), that are sensitive to deriving all nine Cij, are used for compressional and shear wave velocity (νP and νS) measurements as a function of azimuthal angles over 200° at each experimental pressure. Our results show that all Cij of singe-crystal Fe10-Al14-Bgm increase monotonically with pressure with small uncertainties of 1–2% (±1σ), except C55 and C23, which have uncertainties of 3–4%. Using the third-order Eulerian finite-strain equations to model the elasticity data yields the aggregate adiabatic bulk and shear moduli and respective pressure derivatives at the reference pressure of 25 GPa: KS = 326 ± 4 GPa, µ = 211 ± 2 GPa, KS′ = 3.32 ± 0.04, and µ′ = 1.66 ± 0.02 GPa. The high-pressure aggregate νS and νP of Fe10-Al14-Bgm are 2.6–3.5% and 3.1–4.7% lower than those of MgSiO3 bridgmanite end-member, respectively. These data are used with literature reports on bridgmanite with different Fe and Al contents to quantitatively evaluate pressure and compositional effects on their elastic properties. Comparing with one-dimensional seismic profiles, our modeled velocity profiles of major lower-mantle mineral assemblages at relevant P-T suggest that the lower mantle could likely consist of about 89 vol% (Al,Fe)-bearing bridgmanite. After considering uncertainties, our best-fit model is still indistinguishable from pyrolitic or chondritic models.

Geochemistry & Geophysics↗

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa

Abstract Background Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials. Results Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions. Conclusion Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

59 BASIC BIOLOGICAL SCIENCES↗

NGEE Arctic Phase 4 Plant Functional Type Framework for Pan-Arctic Vegetation

The NGEE-Arctic research team identified a common set of hierarchical plant functional types (PFTs) for pan-arctic vegetation that we will use across our research activities. Interdisciplinary work within a large team requires agreement regarding levels of functional organization so that knowledge, data, and technologies can be shared and combined effectively. The team has identified plant functional types as a crucial area where such interoperability is needed. PFTs are used to represent plant pools and fluxes within models, summarize observational data, and map vegetation across the landscape. Within each of these applications, varying levels of PFT specificity are needed according to the specific scientific research goal, computational limitations, and data availability. By agreeing on a specific hierarchical framework for grouping variables in our vegetation data, we ensure the resulting research products will be robust, flexible, and scalable. In this document, we lay out the agreed upon PFT framework with definitions and references to existing literature. Table 1 included in the "NGA700_Phase4PFTFramework_about*" file outlines the relationship between NGEE-Arctic Phase 4, Tier 1 PFTs and the PFTs used within prominent arctic literature as well as publications by the NGEE-Arctic team during phases 1-3.This dataset consists of a table detailing a hierarchical PFT framework that spans 4 tiers with the most granular PFTs listed in tier 1 and the most general PFTs in tier 4. The PFTs within each tier has a single column in the dataset where the PFTs are named and a separate column where the characteristics used to define that PFT are listed. Grey fill of the cells is used to indicate where a given PFT starts to “lose” tier 1 details as you look from left to right. Note the excel file has merged cells to indicate grouping of PFTs across the Tiers- it will not translate into a delimited filetype (.csv, .txt, etc) without modification thus the hierarchical PFT framework table is available in three different file formats: 1) NGA700_Phase4PTS.xlsx – maintains the merged cells and grey fill; 2) NGA700_Phase4PTS.csv – merged cells are split, and grey fill is removed; 3) NGA700_Phase4PTS.pdf – image of the table with merged cells and grey fill. Metadata document included as a *.pdf and file-level metadata and data dictionary as *.csv files.

54 ENVIRONMENTAL SCIENCES↗