Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulation”

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

NbZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase Transitions of Eutectic High Entropy Alloy AlCoCrFeNi 2.1 Under Shock Compression

High entropy alloys (HEAs) are a new class of metals that exhibit unique mechanical performance. Among HEAs, additively manufactured eutectic high entropy alloys (AM EHEAs) have recently emerged as candidate materials for use in extreme conditions due to their simultaneous high strength and ductility. However, the deformation and structural evolution of AM EHEAs under conditions of high pressure have not been well characterized, limiting their use in extreme applications. Dynamic compression experiments and molecular dynamics simulations are presented to study the structural evolution of AM EHEA AlCoCrFeNi 2.1 when compressed to pressures up to 400 GPa. In situ X-ray diffraction measurements capture the appearance of face-centered cubic and body-centered cubic phases at different pressure conditions, with pure- and mixed-phase regions. Furthermore, understanding the phase stability and structural evolution of the AM EHEA offers new insights to guide the development of high-performance complex materials for extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Highly-Cyclable Room-Temperature Phosphorene Polymer Electrolyte Composites for Li Metal Batteries

Despite significant interest toward solid-state electrolytes owing to their superior safety in comparison to liquid-based electrolytes, sluggish ion diffusion and high interfacial resistance limit their application in durable and high-power density batteries. Here, a novel quasi-solid Li + ion conductive nanocomposite polymer electrolyte containing black phosphorous (BP) nanosheets is reported. The developed electrolyte is successfully cycled against Li metal (over 550 h cycling) at 1 mA cm -2 at room temperature. The cycling overpotential is dropped by 75% in comparison to BP-free polymer composite electrolyte indicating lower interfacial resistance at the electrode/electrolyte interfaces. Molecular dynamics simulations reveal that the coordination number of Li + ions around (trifluoromethanesulfonyl)imide (TFSI - ) pairs and ethylene-oxide chains decreases at the Li metal/electrolyte interface, which facilitates the Li + transport through the polymer host. Here, density functional theory calculations confirm that the adsorption of the LiTFSI molecules at the BP surface leads to the weakening of N and Li atomic bonding and enhances the dissociation of Li + ions. This work offers a new potential mechanism to tune the bulk and interfacial ionic conductivity of solid-state electrolytes that may lead to a new generation of lithium polymer batteries with high ionic conduction kinetics and stable long-life cycling.

2D materials↗
DOE OSTI2020

Visualizing Energy Transfer at Buried Interfaces in Layered Materials Using Picosecond X-Rays

Understanding the fundamentals of nanoscale heat propagation is crucial for next-generation electronics. For instance, weak van der Waals bonds of layered materials are known to limit their thermal boundary conductance (TBC), presenting a heat dissipation bottleneck. In this paper we present a new non-destructive method to probe heat transport in nanoscale crystalline materials using time-resolved x-ray measurements of photo-induced thermal strain. This technique directly monitors time-dependent temperature changes in the crystal and the subsequent relaxation across buried interfaces by measuring changes in the c-axis lattice spacing after optical excitation. We investigate films of five different layered transition metal dichalcogenides MoX 2 [X= S, Se, Te] and WX 2 [X= S, Se] as well as graphite and a W-doped alloy of MoTe 2 , on c-plane sapphire substrates, finding TBC values in the range 10-30 MWm -2 K -1 at room temperature. In conjunction with molecular dynamics simulations, we show that the high thermal resistances are a consequence of weak interfacial van der Waals bonding and low phonon irradiance. This work paves the way for an improved understanding of thermal bottlenecks in emerging three-dimensional heterogeneously integrated technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Can Polyelectrolyte Mechanical Properties be Directly Modulated by an Electric Field? A Molecular Dynamics Study

Multifunctionality in polymers facilitates their application in emerging technologies. Electrical fields are a preferred stimulus because of the speed and ease of application to bulk polymers. While a wide range of electrically triggered actuators are developed, and electrically controlled adhesion between gels is demonstrated, modification of bulk mechanical properties via electrical stimuli remains elusive. Polymers with covalently incorporated ionic charge (polyelectrolytes) should be well suited to achieving this goal since the mechanical properties depend on electrostatic interactions and these charges are intrinsically susceptible to electric fields. Molecular dynamics simulations are utilized here to investigate whether electric fields can modulate the mechanical properties of polyelectrolytes and to understand the governing mechanisms. Mechanical property modulation by electric field is found to be sensitive to the charge distribution—charges must be tightly attached to the polymer backbone, and responsivity is greater if a single backbone contains both positive and negative charges. The dominant mechanisms are reorientation and stretching of the polymer chains, which also elongate the ionic clusters to maintain strong electrostatic interactions throughout deformation. Furthermore, these insights are critical for future experimental realization of polymers with electric field regulated mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Engineering the Interlayer Spacing by Pre-Intercalation for High Performance Supercapacitor MXene Electrodes in Room Temperature Ionic Liquid

MXenes exhibit excellent capacitance at high scan rates in sulfuric acid aqueous electrolytes, but the narrow potential window of aqueous electrolytes limits the energy density. Organic electrolytes and room-temperature ionic liquids (RTILs) can provide higher potential windows, leading to higher energy density. The large cation size of RTIL hinders its intercalation in-between the layers of MXene limiting the specific capacitance in comparison to aqueous electrolytes. In this work, different chain lengths alkylammonium (AA) cations are intercalated into Ti 3 C 2 T x , producing variation of MXene interlayer spacings (d-spacing). AA-cation-intercalated Ti 3 C 2 T x (AA-Ti 3 C 2 ), exhibits higher specific capacitances, and cycling stabilities than pristine Ti 3 C 2 T x in 1 m 1-ethly-3-methylimidazolium bis-(trifluoromethylsulfonyl)-imide (EMIMTFSI) in acetonitrile and neat EMIMTFSI RTIL electrolytes. Pre-intercalated MXene with an interlayer spacing of ≈2.2 nm, can deliver a large specific capacitance of 257 F g –1 (1428 mF cm –2 and 492 F cm –3 ) in neat EMIMTFSI electrolyte leading to high energy density. Quasi elastic neutron scattering and electrochemical impedance spectroscopy are used to study the dynamics of confined RTIL in pre-intercalated MXene. Furthermore, molecular dynamics simulations suggest significant differences in the structures of RTIL ions and AA cations inside the Ti 3 C 2 T x interlayer, providing insights into the differences in the observed electrochemical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Riddle of Dark LLZO: Cobalt Diffusion in Garnet Separators of Solid‐State Lithium Batteries

Abstract Solid‐state batteries (SSBs) with a Li 7 La 3 Zr 2 O 12 (LLZO) garnet electrolyte are attracting much attention as robust and safe alternative to conventional lithium‐ion batteries. Technical challenges in the practical implementation of garnet SSBs are related to the need for high‐temperature sintering, which often leads to undesirable chemical reactions with the cathode material. While these reactions are well understood for composite cathodes, very little is known about similar processes between cathode and separator during battery fabrication. This work focuses on understanding the processes between the composite LiCoO 2 ‐LLZO cathode and the LLZO separator and how they affect the battery performance. The extensive diffusion of Co‐ions within LLZO, which leads to the often‐observed LLZO darkening, is shown to have a significant impact on ionic conductivity, electronic conductivity, and dendrite stability of the separator. Experimental data coupled with large‐scale molecular dynamics simulations uncover the diffusion mechanism for Co‐ions and identify secondary phases that form during these interactions. In addition to extensive Co‐ion diffusion within the grains, a non‐uniform segregation of Co‐ions at grain boundaries is found leading to the formation of three distinct Co‐containing phases. This work offers a general approach to studying the fundamental ion diffusion processes that occur during the fabrication of oxide SSBs.

electrochemical properties↗
DOE OSTI2024

Operando Investigation of the Molecular Origins of Dipole Switching in P(VDF‐TrFE‐CFE) Terpolymer for Large Adiabatic Temperature Change

Abstract Relaxor ferroelectric polymers exhibiting a giant electrocaloric effect (ECE) can potentially be used to create next‐generation solid‐state coolers. Under an electric field, poly(vinylidene fluoride‐trifluoroethylene‐chlorofluoroethylene) terpolymer goes through a large dipolar entropy change producing a high adiabatic temperature change (Δ T ECE ). This work resolves the molecular origins of the large entropy change behind the electric field‐induced dipole switching. A Fourier transform infrared spectroscopy equipped with a high voltage source is used to operandoly observe the characteristic molecular vibrational modes. A short‐range trans (T) conformation of the CF 2 ‐CH 2 dyads interrupted by a gauche (G) conformation, e.g., TTTG in the terpolymer chain, undergoes a dynamic transformation that leads to a corresponding Δ T ECE whenever an electric field is applied. The molecular dynamics simulation also proves that the energy barrier that the transformation from TTTGs into a long T sequence overcomes is smaller than that for all other conformations. A mixed solvent system is used to obtain T3G‐enriched terpolymer films exhibiting a 4.02 K Δ T ECE at 60 MV m −1 and these films are employed to manufacture a 2‐layer‐cascaded cooling device that achieves a 6.7 K temperature lift, the highest reported value for a 2‐layer cascaded device made of fluoropolymers.

Zhu, Yuan↗
DOE OSTI2024

Transforming an Ionic Conductor into an Electronic Conductor via Crystallization: In Situ Evolution of Transference Numbers and Structure in (La,Sr)(Ga,Fe)O 3–x Perovskite Thin Films

Mixed-conducting perovskites are workhorse electrochemically active materials, but typical high-temperature processing compromises their catalytic activity and chemo-mechanical integrity. Low-temperature pulsed laser deposition of amorphous films plus mild thermal annealing is an emerging route to form homogeneous mixed conductors with exceptional catalytic activity, but little is known about the evolution of the oxide-ion transport and transference numbers during crystallization. Here the coupled evolution of ionic and electronic transport behavior and structure in room-temperature-grown amorphous (La,Sr)(Ga,Fe)O 3-x films as they crystallize is explored. In situ ac-impedance spectroscopy with and without blocking electrodes, simultaneous capturing-synchrotron-grazing-incidence X-ray diffraction, dc polarization, transmission electron microscopy, and molecular dynamics simulations to evaluate isothermal and non-isothermal crystallization effects and the role of grain boundaries on transference numbers is combined. Ionic conductivity increases by ≈2 orders of magnitude during crystallization, with even larger increases in electronic conductivity. Consequently, as crystallinity increases, LSGF transitions from a predominantly ionic conductor to a predominantly electronic conductor. The roles of evolving lattice structural order, microstructure, and defect chemistry are examined. Grain boundaries appear relatively nonblocking electronically but significantly blocking ionically. The results demonstrate that ionic transference numbers can be tailored over a wide range by tuning crystallinity and microstructure without having to change the cation composition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanism of Antiferroelectricity in Polycrystalline ZrO 2

The size and electric field dependent induction of polarization in antiferroelectric ZrO 2 is the key to several technological applications that are unimaginable a decade ago. However, the lack of a deeper understanding of the mechanism hinders progress. Molecular dynamics simulations of polycrystalline ZrO 2 , based on machine-learned interatomic forces with near ab initio quality, shed light on the fundamental mechanism of the size effect on the transition fields. Stress in the oxygen sublattice is the most important factor. The so constructed interatomic forces allow the calculation of the transition fields as a function of the ZrO 2 film thickness and predict the ferroelectricity at large thickness. The simulation results are validated with electrical and piezo response force microscopy measurements. The results allow a clear interpretation of the properties of the double-hysteresis loops as well as the construction of the free energy landscape of ZrO 2 grains.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nanoscale Editing of Multi and Single Layer Tungsten Disulfide via Gas‐Assisted Focused Electron Beam Induced Etching for Device Prototyping

Focused electron beam induced etching (FEBIE) with XeF 2 (xenon difluoride) precursor is conducted on multi-layer exfoliated WS 2 (tungsten disulfide) and monolayer WS 2 grown by chemical vapor deposition (CVD). The films are characterized by atomic force microscopy (AFM) and Raman and photoluminescence (PL) spectroscopy post-etching. The etch rates/efficiencies are reported as a function of electron beam energy, current, dwell time, and XeF 2 pressure. Bulk film Raman spectra are unchanged post-FEBIE, indicating minimal subsurface damage. Monolayer WS 2 shows a decrease in Raman and PL intensity post-FEBIE, with a dose-to-clear of ≈2 nC µm −2 . The study reveals regimes affected by the various mass transport contributions such as refresh time and the ratio of electrons/XeF 2 . Spontaneous etching was discovered during FEBIE of large patterned areas due to the long frame/refresh times. Density functional theory and ab initio molecular dynamics simulations compares desorption of SF x and WF x molecules from pristine WS 2 basal planes and pore edges, revealing the spontaneous etching is consistent with etching of partially etched monolayers during each frame. Single-line etching width of 21 nm, and patterning flakes into 100 nm wide channels are demonstrated. In conclusion, this work demonstrates the possibility of editing WS 2 flakes into electronic devices of arbitrary dimensions for semiconductor applications.

2D materials↗
DOE OSTI2021

Interfacial Atomistic Mechanisms of Lithium Metal Stripping and Plating in Solid–State Batteries

All-solid-state batteries based on a Li metal anode represent a promising next-generation energy storage system, but are currently limited by low current density and short cycle life. Further research to improve the Li metal anode is impeded by the lack of understanding in its failure mechanisms at lithium–solid interfaces, in particular, the fundamental atomistic processes responsible for interface failure. Here, using large-scale molecular dynamics simulations, the first atomistic modeling study of lithium stripping and plating on a solid electrolyte is performed by explicitly considering key fundamental atomistic processes and interface atomistic structures. In the simulations, the interface failure initiated with the formation of nano-sized pores, and how interface structures, lithium diffusion, adhesion energy, and applied pressure affect interface failure during Li cycling are observed. By systematically varying the parameters of solid-state lithium cells in the simulations, the parameter space of applied pressures and interfacial adhesion energies that inhibit interface failure during cycling are mapped to guide selection of solid-state cells. Furthermore, this study establishes the atomistic modeling for Li stripping and plating, and predicts optimal solid interfaces and new strategies for the future research and development of solid-state Li-metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Quantifying Anionic Diffusion in 2D Halide Perovskite Lateral Heterostructures

Anionic diffusion strongly impacts the stability of halide perovskite materials, but it is still not well understood. Here, a quantitative investigation of in-plane thermally driven anionic inter-diffusion in a series of novel 2D and quasi-2D halide perovskites lateral heterostructures is reported. The calculated diffusion coefficients (D) reveal the inhibition of Br–I inter-diffusion with bulky π-conjugated organic cations compared with short-chain aliphatic organic cations. Furthermore, halide diffusion is found to be faster in quasi-2D (n > 1) than 2D perovskites (n = 1). The increment becomes less apparent as the “n” number increases, akin to the quantum confinement effect observed for band gaps. These trends are rationalized by molecular dynamics simulations of free energy barriers for halide diffusion that reveal mechanisms for suppressing diffusion. We report this work provides important fundamental insights on the anionic migration and diffusion process in halide perovskite materials.

2D materials↗
DOE OSTI2022

Controlled Formation of Conduction Channels in Memristive Devices Observed by X–ray Multimodal Imaging

Neuromorphic computing provides a means for achieving faster and more energy efficient computations than conventional digital computers for artificial intelligence (AI). However, its current accuracy is generally less than the dominant software-based AI. The key to improving accuracy is to reduce the intrinsic randomness of memristive devices, emulating synapses in the brain for neuromorphic computing. Here using a planar device as a model system, the controlled formation of conduction channels is achieved with high oxygen vacancy concentrations through the design of sharp protrusions in the electrode gap, as observed by X-ray multimodal imaging of both oxygen stoichiometry and crystallinity. Classical molecular dynamics simulations confirm that the controlled formation of conduction channels arises from confinement of the electric field, yielding a reproducible spatial distribution of oxygen vacancies across switching cycles. Furthermore, this work demonstrates an effective route to control the otherwise random electroforming process by electrode design, facilitating the development of more accurate memristive devices for neuromorphic computing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Atomic Drill Bit: Precision Controlled Atomic Fabrication of 2D Materials

The ability to deterministically fabricate nanoscale architectures with atomic precision is the central goal of nanotechnology, whereby highly localized changes in the atomic structure can be exploited to control device properties at their fundamental physical limit. Here, an automated, feedback-controlled atomic fabrication method is reported and the formation of 1D–2D heterostructures in MoS 2 is demonstrated through selective transformations along specific crystallographic orientations. The atomic-scale probe of an aberration-corrected scanning transmission electron microscope (STEM) is used, and the shape and symmetry of the scan pathway relative to the sample orientation are controlled. The focused and shaped electron beam is used to reliably create Mo 6 S 6 nanowire (MoS-NW) terminated metallic-semiconductor 1D–2D edge structures within a pristine MoS 2 monolayer with atomic precision. From these results, it is found that a triangular beam path aligned along the zig-zag sulfur terminated (ZZS) direction forms stable MoS-NW edge structures with the highest degree of fidelity without resulting in disordering of the surrounding MoS 2 monolayer. Density functional theory (DFT) calculations and ab initio molecular dynamic simulations (AIMD) are used to calculate the energetic barriers for the most stable atomic edge structures and atomic transformation pathways. These discoveries provide an automated method to improve understanding of atomic-scale transformations while opening a pathway toward more precise atomic-scale engineering of materials.

2D materials↗
DOE OSTI2023

Engineering Relaxor Behavior in (BaTiO 3 ) n /(SrTiO 3 ) n Superlattices

Complex-oxide superlattices provide a pathway to numerous emergent phenomena because of the juxtaposition of disparate properties and the strong interfacial interactions in these unit-cell-precise structures. This is particularly true in superlattices of ferroelectric and dielectric materials, wherein new forms of ferroelectricity, exotic dipolar textures, and distinctive domain structures can be produced. Here, in this study, relaxor-like behavior, typically associated with the chemical inhomogeneity and complexity of solid solutions, is observed in (BaTiO 3 ) n /(SrTiO 3 ) n (n = 4–20 unit cells) superlattices. Dielectric studies and subsequent Vogel–Fulcher analysis show significant frequency dispersion of the dielectric maximum across a range of periodicities, with enhanced dielectric constant and more robust relaxor behavior for smaller period n. Bond-valence molecular-dynamics simulations predict the relaxor-like behavior observed experimentally, and interpretations of the polar patterns via 2D discrete-wavelet transforms in shorter-period superlattices suggest that the relaxor behavior arises from shape variations of the dipolar configurations, in contrast to frozen antipolar stripe domains in longer-period superlattices (n = 16). Moreover, the size and shape of the dipolar configurations are tuned by superlattice periodicity, thus providing a definitive design strategy to use superlattice layering to create relaxor-like behavior which may expand the ability to control desired properties in these complex systems.

36 MATERIALS SCIENCE↗