Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control Samples”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

3D printed carbon fiber reinforced carbon as an energy efficient alternative to graphite for EFAS tooling

As Electric Field Assisted Sintering (EFAS) gains more industrial acceptance and use, it becomes more important to develop more efficient means to implement this technology. To this aim, 3D printed continuous carbon fiber reinforced carbon (CCC) was manufactured and fabricated into tooling for EFAS systems as an alternative to traditional graphite tooling. The impact of fiber orientation on the thermal and electrical properties of the CCC was characterized. Sample material was sintered in Tokai G535 graphite tooling, under common processing conditions and compared with CCC tooling. There was nearly 50 % energy savings compared to graphite while maintaining equivalent sample density and microstructure plus keeping ram temperatures 39 % cooler. This is due to spatial control of generated heat and thermal diffusivity within the molds, by means of fiber orientation anisotropy. Finite element modeling of the tooling design supported the experimental results as well as displays the effect of optimization of this 3D printed CCC material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Controlled Experiments and Optimized Theory of Absorption Spectra of Li Metal and Salts

Investigation of Li metal and ionic compounds through experimental and theoretical spectroscopy has been of tremendous interest due to their prospective applications in Li-metal and Li-ion batteries. Li K-edge soft X-ray absorption spectroscopy (sXAS) provides the most direct spectroscopic characterization; unfortunately, due to the low core-level energy and the highly reactive surface, Li-K sXAS of Li metal has been extremely challenging, as evidenced by many controversial reports. Here, through controlled and ultra-high energy resolution experiments of two kinds of in situ prepared samples, we report the intrinsic Li-K sXAS of Li-metal that displays a prominent leading peak that has not been revealed before. Furthermore, theoretical simulations show that, due to the low number of valence electrons in Li, the Li-K sXAS is strongly affected by the response of the valence electrons to the core hole. We successfully reproduce the Li-K sXAS by state-of-the-art calculations with considerations of a number of relevant parameters such as temperature, energy resolution, and, especially, contributions from transitions which are forbidden in the single-particle treatment. Such a comparative experimental and theoretical investigation is further extended to a series of Li ionic compounds, which highlight the importance of considering the total and single-particle energies for obtaining an accurate alignment of the spectra. Furthermore, our work provides the first reliable Li-K sXAS of the Li metal surface with advanced theoretical calculations. The experimental and theoretical results provide a critical benchmark for studying Li chemistry in both metallic and ionic states.

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

Chemical Activation of Frontal Ring-Opening Metathesis Polymerization

In this study, we present an innovative approach to frontal ring-opening metathesis polymerization (FROMP) through chemical activation. By locally concentrating Grubbs-type initiators, we accelerate exothermic polymerization at the activation site. Sufficient thermal energy initiates a polymerization front that travels through the remainder of the sample even with reduced initiator concentrations in the bulk resin. Initiation location and timing are controlled by adjusting the initiator reactivity, the properties of the transport solvent, and the air–resin interface. Our technique enables multipoint front activation, either through synchronized injection or by timing the relative reactivity of various initiators. This accessible method of chemical activation enables on-demand production of polymer and composite materials using only the thermal energy generated by the reaction, removing the need for external power sources and enabling deployment in remote or infrastructure-limited environments.

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

Flat band carrier confinement in magic-angle twisted bilayer graphene

Abstract Magic-angle twisted bilayer graphene has emerged as a powerful platform for studying strongly correlated electron physics, owing to its almost dispersionless low-energy bands and the ability to tune the band filling by electrostatic gating. Techniques to control the twist angle between graphene layers have led to rapid experimental progress but improving sample quality is essential for separating the delicate correlated electron physics from disorder effects. Owing to the 2D nature of the system and the relatively low carrier density, the samples are highly susceptible to small doping inhomogeneity which can drastically modify the local potential landscape. This potential disorder is distinct from the twist angle variation which has been studied elsewhere. Here, by using low temperature scanning tunneling spectroscopy and planar tunneling junction measurements, we demonstrate that flat bands in twisted bilayer graphene can amplify small doping inhomogeneity that surprisingly leads to carrier confinement, which in graphene could previously only be realized in the presence of a strong magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rapid sensing of hidden objects and defects using a single-pixel diffractive terahertz sensor

Abstract Terahertz waves offer advantages for nondestructive detection of hidden objects/defects in materials, as they can penetrate most optically-opaque materials. However, existing terahertz inspection systems face throughput and accuracy restrictions due to their limited imaging speed and resolution. Furthermore, machine-vision-based systems using large-pixel-count imaging encounter bottlenecks due to their data storage, transmission and processing requirements. Here, we report a diffractive sensor that rapidly detects hidden defects/objects within a 3D sample using a single-pixel terahertz detector, eliminating sample scanning or image formation/processing. Leveraging deep-learning-optimized diffractive layers, this diffractive sensor can all-optically probe the 3D structural information of samples by outputting a spectrum, directly indicating the presence/absence of hidden structures or defects. We experimentally validated this framework using a single-pixel terahertz time-domain spectroscopy set-up and 3D-printed diffractive layers, successfully detecting unknown hidden defects inside silicon samples. This technique is valuable for applications including security screening, biomedical sensing and industrial quality control.

Li, Jingxi (ORCID:0000000165958680)↗
DOE OSTI2020

Creation of discrete active site domains via mesoporous silica poly(styrene) composite materials for incompatible acid–base cascade reactions

This work highlights the design and synthesis of bifunctional mesoporous silicate – polymer composite dual acid–base supported cascade catalysts. Compartmentalization of the two incompatible active sites is sought by segregating acid sites on the silica surface, and base sites within polymer chains and/or polymer domains. The ability to isolate and segregate active sites via control of the mesoporous silica pore size and polymer molecular weight is probed with silica samples functionalized by a grafting-to process. Supplemental activator and reducing agent (SARA) atom transfer radical polymerization is used to synthesize random copolymers containing protected primary amines. Thiol–ene ‘click’ chemistry facilitates silica functionalization via a convergent approach, with the ene-functionalized polymer end group and silica-grafted thiols forming SBA/MCM-SH-poly(styrene-co-2-(4-vinylbenzyl)isoindoline-1,3-dione). Polymer deprotection and thiol oxidation produces primary amine/sulfonic acid containing composite catalysts. With the polymer supported Lewis base and silica grafted Brønsted acid, the two-step deacetalization – Knoevenagel condensation cascade is explored to assess the ability of these polymer/silica hybrids to segregate active sites, allowing both acid and base site accessibility. Six composite catalysts are synthesized and tested in individual and cascade reactions with kinetic results demonstrating that lower molecular weight SBA-15-P1 and MCM-41-P1 catalysts outperform (higher turnover frequencies and initial rates) their higher molecular weight analogues, as well as a polymer-free system containing molecular active sites dispersed on the silica surface. Higher molecular weight composite catalysts perform more poorly due to limited chain solubility, mass transfer limitations, and poor catalyst accessibility. In many cases, the polymer chains effectively thread into the mesopores, with higher molecular weight polymers leading to pore blockage and inhibited mass transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Distinguishing photo-induced oxygen attack on alkyl chain versus conjugated backbone for alkylthienyl-benzodithiophene (BDTT)-based push–pull polymers

Synthetic design has enabled increasing power conversion efficiency advances in organic photovoltaics (OPV). One continuing knowledge gap is detailed understanding of single-material chemical (photo)stability. Many considerations are based on prior OPV-relevant donor homopolymer systems rather than the next-generation push–pull architectures. Generally, energetic offsets between the lowest occupied molecular orbital of the donor and molecular oxygen are assumed to dictate kinetics of photo-induced charge transfer to a super oxide radical. Herein we determine the ambient-induced photo-degradation pathways, presented as proposed site-specific arrow-pushing mechanisms, for five donor polymers all containing the same push unit – alkylthienyl-substituted-benzodithiophene (BDTT) – but with chemically distinct pull units. The donor-only polymer films were subject to controlled photobleaching in air as an accelerated degradation approach, coupled with simultaneous monitoring of absorptance. Sample subsets were periodically removed and analyzed with X-ray photoelectron spectroscopy, to evaluate near-surface chemical composition and oxygen additions to hetero reporter atoms on either the BDTT push unit or distinct pull unit (i.e., sulphur and nitrogen). Here, this methodology allows us to distinguish between different mechanisms of bond cleavage and formation. Overall, we find that neither polymer redox properties nor individual push or pull unit stability are sufficient to predict photo-oxidative degradation of these polymers. Rather, there is a greater dependence on the susceptibility of unique structural groups within a polymeric system. Alkyl chain oxygen addition is generally the first attack site and direct sulfur oxidation on the conjugated backbone occurs after saturation of the alkyl chain initiation sites. This work provides a standard that could be used to evaluate relative photo-oxidative (in)stability for new OPV materials quickly – prior to time-consuming device optimization – and demonstrates an effective methodology for correlating optical degradation with chemical structure alterations via spectroscopic signatures to guide synthetic design.

14 SOLAR ENERGY↗
DOE OSTI2025

Generation of strong fields with subcritical density plasmas to study the phase transitions of magnetized warm dense matter

Warm dense matter (WDM) is a regime where Fermi degenerate electrons play an important role in the macroscopic properties of a material. Recent experiments have brought us closer to understanding unmagnetized processes in WDM, but magnetized WDM remains unexplored because kilotesla magnetic fields are required. Although there are examples of field compression generating such fields by imploding pre-magnetized targets, these existing methods give no independent control over the parameters of the magnetized plasma and result in limited laser access for sample creation and diagnosis. In this paper, numerical simulations show that kilotesla magnetic fields can be obtained by shining laser beams onto the inner surface of a cylindrical target, rather than on the outer surface. This approach relies on field compression by a low-density, high-temperature plasma, rather than a high-density, low-temperature plasma, used in the more conventional approach. With this novel configuration, the region of peak magnetic field is mostly free of plasma, hence, other beams can reach a sample placed in the region of the peak field to form WDM and diagnose it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Alloying of Pu-Al with Stainless Steel for Material Disposition

The Savannah River National Laboratory (SRNL) has evaluated several options for the disposition of stainless-steel clad plutonium metal alloy. One of the technologies under consideration is alloying of the material with stainless steel (SS). The resulting SS-Pu alloy would be a non-proliferable waste form consisting of a secondary Pu composition region microencapsulated in the refractory stainless steel. Two 8-kg ingots were made at SS-1.8Zr-0.4Pu alloys (wt %); 8 kg was determined in a previous study to be the maximum mass of SS ingot at the maximum target Pu loading of 350-g that would result in a SS-4.4Pu alloy (wt %). Two smaller 500-g ingots were also produced at SS-1.6Zr-1.4Pu and SS-1.4Pu (wt %). The 500-g ingots evaluated alloying at a higher Pu concentration than in the 8 kg ingots, and they evaluated the necessity of adding Zr metal to incorporate the Pu and control Pu oxidation. Zr addition was found to be unnecessary to incorporate the Pu and control Pu oxidation. Drill turnings were collected from the large and small ingots and metallographic samples were directly cut from the small ingots. Both were analyzed to validate the structure and composition region formation. Chemical analyses of turnings proved that the Pu was dispersed within the SS ingots.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detection of Part Porosity in Additive Manufacturing Using Machining

Monitoring of part porosity in additive manufacturing is critical for part qualification and process parameter development. This article describes the monitoring of porosity in additive manufacturing using machining. Grade 91 steel samples with different levels of porosity were additively manufactured using a laser blown-powder process. The samples were machined, and the tool position and cutting power were monitored using the machine controller and an external power transducer, respectively. Results show that the power transducer is sensitive to machining porosities. A novel approach using the standard deviation of cutting power was developed to determine the location of porosities along the length of the machining pass. Test cuts were completed at different machining process parameters. Results showed that the probability of detection was one for pore volume greater than 0.3 mm 3 at different machining parameters. In conclusion, although the detection of porosity using machining is a destructive method for evaluation, it offers a low-cost and fast approach for monitoring porosity in additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-bandwidth image-based predictive laser stabilization via optimized Fourier filters

Controlling the delivery of kHz-class pulsed lasers is of interest in a variety of industrial and scientific applications, from next-generation laser-plasma acceleration to laser-based x-ray emission and high-precision manufacturing. The transverse position of the laser pulse train on the application target is often subject to fluctuations by external drivers (e.g., room cooling and heating systems, motorized optics stages and mounts, vacuum systems, chillers, and/or ground vibrations). For typical situations where the disturbance spectrum exhibits discrete peaks on top of a broad-bandwidth lower-frequency background, traditional PID (proportional-integral-derivative) controllers may struggle, since as a general rule PID controllers can be used to suppress vibrations up to only about 5%–10% of the sampling frequency. Here, a predictive feed-forward algorithm is presented that significantly enhances the stabilization bandwidth in such laser systems (up to the Nyquist limit at half the sampling frequency) by online identification and filtering of one or a few discrete frequencies using optimized Fourier filters. Furthermore, the system architecture demonstrated here uses off-the-shelf CMOS cameras and piezo-electric actuated mirrors connected to a standard PC to process the alignment images and implement the algorithm. To avoid high-end, high-cost components, a machine-learning-based model of the piezo mirror’s dynamics was integrated into the system, which enables high-precision positioning by compensating for hysteresis and other hardware-induced effects. A successful demonstration of the method was performed on a 1 kHz laser pulse train, where externally-induced vibrations of up to 400 Hz were attenuated by a factor of five, far exceeding what could be done with a standard PID scheme.

Natal, Joseph↗
DOE OSTI2026

Mono-mix strategy enables comparative proteomics of a cross-kingdom microbial symbiosis

Cross-kingdom microbial symbioses, such as those between algae and bacteria, are key players in biogeochemical cycles. The molecular changes during initiation and establishment of symbiosis are of great interest, but quantitatively monitoring such changes can be challenging, particularly when the microorganisms differ greatly in size or are intimately associated. Here, we analyze output from label-free, data-dependent acquisition (DDA) LC-MS/MS proteomics experiments investigating the well-studied interaction between the alga Chlamydomonas reinhardtii and the heterotrophic bacterium Mesorhizobium japonicum. We found that detection of bacterial proteins decreased in coculture by 50% proteome-wide due to the abundance of algal proteins. As a result, standard differential expression analysis led to numerous false-positive reports of significantly downregulated proteins, where it was not possible to distinguish meaningful biological responses to symbiosis from artifacts of the reduced protein detection in coculture relative to monoculture. We show that data normalization alone does not eliminate the impact of altered detection on differential expression analysis of the cross-kingdom symbiosis. We assessed two additional strategies to overcome this methodological artifact inherent to DDA proteomics. In the first, we combined algal and bacterial monocultures at a relative abundance that mimicked the coculture, creating a “mono-mix” control to which the coculture could be compared. This approach enabled comparable detection of bacterial proteins in the coculture and the monoculture control. In the second strategy, we enhanced detection of lowly abundant bacterial proteins by using sample fractionation upstream of LC-MS/MS analysis. When these simple approaches were combined, they allowed for meaningful comparisons of nearly 10,000 algal proteins and over 4,000 bacterial proteins in response to symbiosis by DDA. They successfully recovered expected changes in the bacterial proteome in response to algal coculture, including upregulation of sugar-binding proteins and transporters. They also revealed novel proteomic responses to coculture that guide hypotheses about algal-bacterial interactions.

Dupuis, Sunnyjoy [University of California, Berkel↗
DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Sludge Batch 11 Assembly: Tank 35

Savannah River Mission Completion Nuclear Safety and Engineering Integration (SRMC-E) has requested that Savannah River National Laboratory (SRNL) perform Tank 35 characterization analyses in support of Sludge Batch 11 (SB11) assembly. Two Tank 35 samples were delivered to SRNL and composited into a single sample in April 2023. The composite sample was analyzed for the following: density, weight percent solids, chemical composition, radionuclides, supernate corrosion control tests, and x-ray diffraction for burkeite, gibbsite, and boehmite.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Sludge Batch 11 Assembly: Tank 35

Savannah River Mission Completion Nuclear Safety and Engineering Integration (SRMC-E) has requested that Savannah River National Laboratory (SRNL) perform Tank 35 characterization analyses in support of Sludge Batch 11 (SB11) assembly. Two Tank 35 samples were delivered to SRNL and composited into a single sample in April 2023. The composite sample was analyzed for the following: density, weight percent solids, chemical composition, radionuclides, supernate corrosion control tests, and x-ray diffraction for burkeite, gibbsite, and boehmite.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Non-Ionizing Millimeter Waves Non-Thermal Radiation of Saccharomyces cerevisiae—Insights and Interactions

Non-ionizing millimeter-waves (MMW) interact with cells in a variety of ways. Here the inhibited cell division effect was investigated using 85–105 GHz MMW irradiation within the International Commission on Non-Ionizing Radiation Protection (ICNIRP) non-thermal 20 mW/cm 2 safety standards. Irradiation using a power density of about 1.0 mW/cm 2 SAR over 5–6 h on 50 cells/μL samples of Saccharomyces cerevisiae model organism resulted in 62% growth rate reduction compared to the control (sham). The effect was specific for 85–105 GHz range and was energy- and cell density-dependent. Irradiation of wild type and Δrad52 (DNA damage repair gene) deleted cells presented no differences of colony growth profiles indicating non-thermal MMW treatment does not cause permanent genetic alterations. Dose versus response relations studied using a standard horn antenna (~1.0 mW/cm 2 ) and compared to that of a compact waveguide (17.17 mW/cm 2 ) for increased power delivery resulted in complete termination of cell division via non-thermal processes supported by temperature rise measurements. We have shown that non-thermal MMW radiation has potential for future use in treatment of yeast related diseases and other targeted biomedical outcomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

SCILLA Secondary Aerosol Volume Concentration Airborne Data

This data set contains secondary aerosol volume concentration in cubic micrometers per cubic centimeter (um^3/cm^3). The aerosol was generated in an oxidation flow reactor (OFR) and its volume concentration measured with a scanning mobility particle sizer (SMPS). Both components were designed and built at the University of California Riverside. Ambient particles were removed with a Teflon filter upstream of the OFR and then pure ammonium sulfate particles were added to create a stable aerosol surface area on which low-volatility oxidation products condense. High concentrations of hydroxyl radical (OH) formed inside the OFR accelerate the oxidative chemistry that would typically occur over a period of several days in the atmosphere, resulting in the production of secondary aerosol from precursor gases present in the ambient air. Concentrations of added ozone and water vapor were controlled to produce the desired and approximately constant level of photochemical aging. The reactor temperature was controlled, while the pressure was not, and was always slightly lower than that of the sampled outside air. Sampled air was pulled from above the Naval Postgraduate School’s (NPS) Twin Otter aircraft through a rear-facing, ¼” OD PFA Teflon tube. The sample stream was split between the OFR and several gas analyzers (NOx, CO, O3, and H2O). Though not contained in these files, data from an aerosol mass spectrometer intermittently operated downstream of the OFR are also available.

{"secondary aerosol concentration",aerosol_concent↗
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↗