Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2022

BioGeoChemistry at Interfaces (LLNL SFA OBER FY22 Program Management and Performance Report)

The focus of the BioGeoChemistry at Interfaces SFA is to identify and quantify the biogeochemical processes and the underlying mechanisms that control actinide mobility in an effort to reliably predict and control the cycling and migration of actinides in the environment. The research approach includes: (1) Field Studies that capture actinide behavior on the timescale of decades (Research Thrust 1), and (2) Fundamental Laboratory Studies that isolate specific biogeochemical processes observed in the field (Research Thrust 2). These Research Thrusts are underpinned by the unique capabilities and staff expertise at Lawrence Livermore National Laboratory, allowing the BioGeoChemistry at Interfaces SFA to advance our understanding of actinide migration behavior in the environment, and serve as an resource for environmental radiochemistry research internationally. In the past year, our greater focus on transient redox gradients across stratified waters, sediment-water interfaces, and mineral-water interfaces extended our research beyond actinides to address processes controlling cycling of redox-sensitive metals more broadly. Nevertheless, Research Thrusts 1 and 2 are guided by the following broad central hypotheses that were developed during our last program review held at the end of FY18: Thrust 1 Hypothesis: Biogeochemical processes occurring on the timescale of years to decades lead to greater actinide recalcitrance in sediments and limits their migration in surface and groundwater. Thrust 2 Hypothesis: Long-term biogeochemical processes include mineral and surface alteration, which leads to stabilization of actinide surface associations or incorporation into mineral precipitates. Our strategic goal is to use the knowledge gained from our Science Plan to advance our understanding of the behavior of actinides and other radionuclides (e.g. Cs) and provide DOE with the scientific basis for remediation and long-term stewardship of DOE’s legacy sites. More broadly, we will enhance our understanding of transport phenomena in environmental systems sciences with a particular emphasis on environmentally relevant (long-term) timescales. While we retain our focus on actinide biogeochemistry, our increased focus on overall biogeochemical processes occurring at unique Test Bed locations associated with this SFA provides fundamental information on abiotic and biotic redox processes that control the cycling of redox sensitive metals under dynamic and transient conditions.

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

The Effects of Morphology and Hydration on Anion Transport in Self-Assembled Nanoporous Membranes

Ordered nanoporous polymer membranes offer opportunities for systematically probing the mechanisms of ion transport under confinement and for realizing useful materials for electrochemical devices. Here, we examine the impact of morphology and ion hydration on the transport of hydroxide and bromide anions in nanostructured polymer membranes with 1 nm scale pores. We use aqueous lyotropic self-assembly of an amphiphilic monomer, with a polymerizable surfactant to create direct hexagonal (H I ) and gyroid mesophases. UV-induced cross-linking leads to the formation of nanoporous polymers with water continuous channels. The membranes are mechanically robust and chemically durable, resisting degradation during extended exposure to 1 M NaOH solutions. We use a combination of electrochemical impedance spectroscopy, pulsed-field gradient NMR spectroscopy, and molecular simulations to elucidate anion and water transport. The as-prepared hexagonal systems display higher conductivity and lower activation energies for both anions relative to the gyroid system. When compared at equivalent hydration, however, gyroid and hexagonal membranes show similar activation energies, with nearly identical conductivities at ambient temperatures. Both ionic conductivity and water diffusivity increase with increasing hydration. Here, the water uptake as a function of relative humidity for the hexagonal and gyroid mesophases ultimately dictates the water diffusion and magnitude of the ionic conductivity, with the hexagonal system showing overall higher capacity for hydration and thus faster ion transport. The durability of these materials under aggressive alkaline conditions and their relatively high hydroxide ion conductivity suggest that these nanostructured polymers could be of interest as membranes for alkaline fuel cells.

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

Ionic Liquids for Direct Air Capture of CO 2 using Electric‐Field‐Mediated Moisture Gradient Process (Final Technical Report)

The final report provides executive summary, a list of publications, and information on training graduate students and postdoctoral researchers. We carried out computational and experimental research on understanding molecular-level mechanism of how CO 2 is absorbed in a solution containing ethylene glycol as the solvent and KOH as the salt in the presence of ionic liquids and under the influence of electric field. In doing so, we developed an automated high-throughput method which allowed us to measure the solubility of CO 2 in a large number of ionic liquids, considerably speeding up the CO 2 solubility measurement. We also demonstrated how varying the concentration of ionic liquids in ethylene glycol can result in a maximum in ionic conductivity. Reaction of CO 2 and subsequent release results in a 50% reduction when the process is operated at an ionic liquid-ethylene glycol concentration yielding maximum ionic conductivity amongst all the ionic liquid-ethylene glycol combinations studied as a part of this research. We utilized machine learning models to identify unique ionic liquid-solvent combinations with ionic conductivity much higher than that measured for ionic liquid-ethylene glycol combinations. We demonstrated that the rate of CO 2 reaction with KOH in ethylene glycol can be optimized with the type of ionic liquid and its concentration. Overall, the research led to publication of 10 peer-reviewed research articles and several presentations at national conferences. We are also in the process of developing additional manuscripts based on the research carried out as a part of this project. Two graduate students and two postdoctoral researchers were supported on the funding.

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

Self-interaction error overbinds water clusters but cancels in structural energy differences

Here, we gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi–Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew–Burke–Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction–corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/ H 2 O from about 42 meV/ H 2 O for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction–corrected SCAN. These findings suggest that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

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

HfZr_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 hafnium-zirconium (Hf-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. 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 Hf 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.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A 2 BCQ 4 ) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q 2 ) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl 2 CdGeSe 4 , Tl 2 CdSnSe 4 , Tl 2 HgSiSe 4 , Tl 2 HgSnS 4 , AuBSe 2 , AuBTe 2 , AuAlTe 2 , AuGaTe 2 , AuInTe 2 , AgAlSe 2 , and AgAlTe 2 are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe 2 and AgSbTe 2 crystals with dopant compositions of AgSb 0.94 Cd 0.06 Te 2 and AgSbTe 1.85 Se 0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

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

HfNb_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 hafnium-niobium (Hf-Nb). 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. 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 Hf and Nb. 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfTa_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 hafnium-tantalum (Hf-Ta). 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. 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 Hf and Ta. 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 sites. 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,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2024

VZr_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 vanadium-zirconium (V-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 V 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.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of Ether-Functionalized Pyrrolidinium Ionic Liquids on Properties and Li + Cation Solvation in Solvate Ionic Liquids

Ionic liquids are tunable solvents composed entirely of ions that have properties desirable as electrolytes for lithium batteries such as non-flammability and a large electrochemical stability window. Solvate ionic liquids are a subclass of ionic liquids that consist of a glyme-based solvent and lithium salt in an equimolar ratio, where Li + cation-glyme solvation interactions result in ionic liquid-like properties. LiG4TFSi is a well-studied solvate ionic liquid consisting of equimolar amounts of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and tetraglyme (G4). In this work, pyrrolidinium ionic liquids with ether-functionalized side chains were synthesized containing either one ether (EO1) moiety or three ether (EO3) moieties and mixed with LiG4TFSI to form a new class of electrolyte mixtures. Their physical and transport properties, as well as ion solvation structures, were characterized by electrochemical, thermal, rheological, and spectroscopic measurements. The conductivity of the electrolyte mixture composed of EO1:LiTFSI:G4 in a 1:1:1 molar ratio is 2.54 mS/cm at 30 °C, compared to 1.53 mS/cm for LiG4TFSI, an increase of 67%. A significant decrease in the conductivity to 0.279 mS/cm is observed for the EO3:LiTFSI:G4 mixture in a 1:1:0.4 molar ratio. Pulsed-field gradient nuclear magnetic resonance (PFG-NMR) measurements revealed that the EO1 cation diffuses significantly faster than the EO3 cation in their respective mixtures. Liquid-state 13 C NMR experiments indicate that Li + cations preferentially coordinate with tetraglyme. Li + cations do not coordinate with the EO1 cation and only coordinate with EO3 ether side chains at lower concentrations of tetraglyme. We hypothesize that the oligoether EO3 cation competes with G4 and TFSI - for lithium cation solvation in G4 deficient compositions, leading to a largely adverse effect on the mass transport properties of the electrolyte.

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

Factors Influencing Preferential Anion Interactions during Solvation of Multivalent Cations in Ethereal Solvents

Most multivalent secondary batteries have employed electrolytes composed of cyclic ether solvents such as tetrahydrofuran or linear glycol ether solvents (glymes) such as 1,2-dimethoxyethane (G1). A robust understanding of multivalent cation solvation tendencies in these classes of solvents provides insight into corresponding structure–property relationships which, in turn, promotes the design and discovery of improved electrolytes. In this work, our goal is to systematically address how electrolyte constituent properties, namely, ether solvent structure and dication size, direct the solvation interactions of divalent electrolytes and their resultant properties. This study utilizes pulsed-field gradient (PFG) nuclear magnetic resonance (NMR) spectroscopy in conjunction with Raman spectroscopy and ionic conductivity measurements to elucidate the preferential interactions between multivalent cations, anions, and solvent molecules along with their correlated ion dynamics. These investigations incorporate two representative divalent cations (Ca 2+ and Zn 2+ ) as well as two ethereal solvent representatives from both the cyclic ether and glyme structural classes. The results reveal that anions coordinate more readily with divalent cations in cyclic ethers than in glymes. Furthermore, the coordination of the anions with Ca 2+ , i.e., contact-ion pair (CIP) formation is more pronounced than with Zn 2+ in a glyme solvent of limited chain length (G1), providing insight into cation size effects that are important for translating solvation behavior across various multivalent electrolytes. Importantly, we find that specific anion coordination is more strongly controlled by solvent structure than by salt concentration in the practical range of 0.1–0.5 M. However, simply reducing these inner-sphere inter-ionic interactions by changing solvent structure does not necessarily de-correlate ionic motion. Instead, concentration-dependent changes in molar ionic conductivity suggest that second-shell interactions, i.e., solvent separated ion pairs (SSIPs), are prevalent in these electrolytes and that the solution dielectric constant, which is increased by the presence of dipolar ion pairs, is critical for controlling these interactions. These findings thus provide a basis for understanding the physical chemistry of multivalent battery electrolytes.

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

Electro-chemo-mechanically Driven Ni Exsolution from (Pr,Ce,Ni)O 2−δ : Controlled Nucleation Density and Enhanced Electrode Kinetics

In situ exsolution of metal nanoparticles is a promising strategy to prepare electrocatalysts with enhanced activity and resistance to agglomeration for efficient chemical transformations and energy conversion. Achieving a high nucleation density of nanoparticles under mild conditions and understanding how to tailor the process is important for performance of these electrodes in electrochemical cells. In this work, we demonstrate facile exsolution of Ni nanoparticles using fluorite-structured (Pr,Ce)O 2−δ as the support oxide, driven by electrochemical potential and aided by the metastability of Ni in the solid solution (elastic driving force). We prepare single-phase oriented thin films of (Pr,Ce,Ni)O 2−δ (NPCO) on (Zr,Y)O 2−δ (YSZ) substrates by pulsed laser deposition. With the aid of a high-throughput electrochemical cell that provides a lateral gradient in Nernst voltage, we apply in situ near-ambient pressure synchrotron X-ray photoelectron spectroscopy and ex situ atomic force microscopy to investigate the impact of electrochemical potential on Ni nucleation density. We find that metallic Ni can be successfully exsolved at 550 °C upon cathodic biasing in 20 mTorr O 2 , and its nucleation density increases with increasing electrochemical driving force/decreasing oxygen chemical potential. We further evaluate the electrochemical performance under highly reducing (fuel electrode) conditions by electrochemical impedance spectroscopy. With the exsolved Ni nanoparticles, the surface exchange coefficient of the NPCO is found to be ∼4× higher than for PCO without exsolution. This work confirms mixed conducting fluorites as beneficial host lattices for facile transition-metal exsolution and suggests the possibility for constructing an all ceria-based electrochemical cell with PCO serving as both the cathode and the anode.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mercury Speciation via Diffusive Gradients Thin-Films Technology

Objectives: Develop and test various diffusive gradient in thin films (DGT) samplers for mercury - SRNL is developing reactive DGTS (rDGTs). Test chemistry to differentiate methylmercury from total and/or inorganic mercury in environmental samples. Fabricate and test rDGT samplers for deployment. Deploy the rDGTs in variable settings. Mercury is a persistent-bioaccumulative-toxic environmental pollutant. Mercury exists in different species such as elemental, inorganic, and methyl. Mercury speciation determine mercury behavior and toxicity. For example, methylmercury strongly accumulates in biota such as fish resulting in potential human health impacts. DGTs are innovative samplers to measure water concentration by diffusion and capture. Types of Hg DGTs used: Total Hg, Inorganic Hg, and Methyl Hg. Deployment Sites: Upper Three Runs, Lower Three Runs, Steele Creek, East Fork Poplar Creek (SRNL Sites), Horizons Creek, NOA Creek, Bruners Creek (ORNL Sites). DGT concentration estimates closely matched alternative measurements for soluble mercury at each Oak Ridge site - DGTs did not respond to particulate mercury. DGTs provide a representative measure for biota uptake since they are left out longer which allows for an average exposure concentration. DMA-80 provided an efficient and quick analysis. Variability in site location made a difference in Hg species. DGT preparations and sampling wasn't as quick and simple as microcolumn technology. Copper reagent degraded agarose gel - more work and cleaner data are needed to assess speciation in rDGTs. Future Work: Use cross-linked polyacrylamide for the collection layer for a stronger plastic. Make a copper diffusion layer without agarose. Test copper separation in lab using realistic stream conditions such as high organic carbon.

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

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exchange-correlation functional challenges in modeling quaternary chalcogenides

The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.

14 SOLAR ENERGY↗
DOE OSTI2025

TaTi_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-titanium (Ta-Ti). 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. 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 Hf and Nb. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗