Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid crystals”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A universal wet-chemistry synthesis of solid-state halide electrolytes for all-solid-state lithium-metal batteries

Solid-state halide electrolytes have gained revived research interests owing to their high ionic conductivity and high-voltage stability. However, synthesizing halide electrolytes from a liquid phase is extremely challenging because of the vulnerability of metal halides to hydrolysis. In this work, ammonium-assisted wet chemistry is reported to synthesize various solid-state halide electrolytes with an exceptional ionic conductivity (>1 microsiemens per centimeter). Microstrain-induced localized microstructure change is found to be beneficial to lithium ion transport in halide electrolytes. Furthermore, the interfacial incompatibility between halide electrolytes and lithium metal is alleviated by transforming the mixed electronic/ionic conductive interface into a lithium ion–conductive interface. Such all-solid-state lithium-metal batteries (ASSLMBs) demonstrate a high initial coulombic efficiency of 98.1% based on lithium cobalt oxide and a high discharge capacity of 166.9 microampere hours per gram based on single-crystal LiNi 0.6 Mn 0.2 Co 0.2 O 2 . This work provides universal approaches in both material synthesis and interface design for developing halide-based ASSLMBs.

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

Beating absorption in solid-state high harmonics

Since the new millennium coherent extreme ultra-violet and soft x-ray radiation has revolutionized the understanding of dynamical physical, chemical and biological systems at the electron’s natural timescale. Unfortunately, coherent laser-based upconversion of infrared photons to vacuum-ultraviolet and soft x-ray high-order harmonics in gaseous, liquid and solid targets is notoriously inefficient. In dense nonlinear media, the limiting factor is strong re-absorption of the generated high-energy photons. Here we overcome this limitation by generating high-order harmonics from a periodic array of thin one-dimensional crystalline silicon ridge waveguides. Adding vacuum gaps between the ridges avoids the high absorption loss of the bulk and results in a ~ 100-fold increase of the extraction depth. As the grating period is varied, each high harmonic shows a different and marked modulation, indicating their waveguiding in the vacuum slots with reduced absorption. Looking ahead, our results enable bright on-chip coherent short-wavelength sources and may extend the usable spectral range of traditional nonlinear crystals to their absorption windows. Potential applications include on-chip chemically-sensitive spectro-nanoscopy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Existence of a flat band in highly correlated metal BaCo 2 P 2

We report on the observation of a flat band situated at the Fermi level E F along with the structural, electrical transport, and magnetic properties of BaCo 2 P 2 that crystallizes in the ThCr 2 Si 2 -type body-center tetragonal structure. This compound has the largest inter-layer pnictide (Pn) distance d Pn–Pn as well as the largest c/a ratio among all the known ACo 2 Pn 2 (A = alkaline earth metal) compounds, where a and c are the tetragonal lattice parameters. Hence, the magnetic and electronic properties of this compound are expected to have a quasi-two-dimensional character. Despite the evidence of the presence of sizable magnetic interactions, magnetic susceptibility χ(T) of BaCo 2 P 2 does not show magnetic ordering down to 1.8 K. The material shows good metallic conduction with a large residual resistivity ratio ρ 300 K /ρ 1.8 K ≈ 70 and a Fermi liquid behavior at low temperature. Kadowaki–Woods ratio R KW of BaCo 2 P 2 suggests the presence of sizable electronic correlations within this system. Additionally, a large many-body enhancement of 2.3 of the experimental density of states D γ (E F ) over the band-structure D band (E F ) is inferred to arise from sizable electron-electron and/or electron-phonon interactions leading to a substantial deviation from the free-electron behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nitrate and Nitrite Incompatibility with Hydroxide Ions in Concentrated NaOH Solutions: Implications for Hydroxide and Gibbsite Reactivity in Alkaline Nuclear Waste

Electrolyte solutions in alkaline nuclear waste contain aluminate, hydroxide, nitrate and nitrite with sodium as the predominant counterion. The salts of these ions are highly soluble, so the liquids are highly concentrated. This study found that there is a substantial incompatibility between the hydroxide and nitrate and/or nitrate ions. This was determined by the observations that adding just one molal of NaNO2 or NaNO3 to saturated NaOH solution precipitation of 12 moles of NaOH·H2O salt, whereas the common ion effect would have expected only about 1 mole to precipitate. Further analysis indicates that the presence of nitrate and nitrite drastically increases the reactivity of sodium hydroxide ions in solution, which likely influences the reactivity of other hydroxide-mediated reactions. This enhanced reactivity is likely because it disrupts large Na+-OH- ion networks because nitrate and nitrite do not fit in those ion networks similarly to how some ions cannot substitute into a foreign crystal lattice. In contrast, the aluminate ion did not have the same large incompatibility with hydroxide.

Nitrite, Nitrate, Aluminate, NaOH*H2O, Ion-aggrega↗
DOE OSTI2025

Extended quantum anomalous Hall states in graphene/hBN moiré superlattices

Electrons in topological flat bands can form new topological states driven by correlation effects. The pentalayer rhombohedral graphene/hexagonal boron nitride (hBN) moiré superlattice was shown to host fractional quantum anomalous Hall effect (FQAHE) at approximately 400 mK, triggering discussions around the underlying mechanism and role of moiré effects. In particular, new electron crystal states with non-trivial topology have been proposed. Here, in this paper, we report electrical transport measurements in rhombohedral pentalayer and tetralayer graphene/hBN moiré superlattices at electronic temperatures down to below 40 mK. We observed two more fractional quantum anomalous Hall (FQAH) states and smaller R xx values in pentalayer devices than those previously reported. In the new tetralayer device, we observed FQAHE at moiré filling factors v = 3/5 and 2/3. With a small current at the base temperature, we observed a new extended quantum anomalous Hall (EQAH) state and magnetic hysteresis, where R xy = h/e 2 and vanishing R xx spans a wide range of v from 0.5 to 1.3. At increased temperature or current, EQAH states disappear and partially transition into the FQAH liquid. Furthermore, we observed displacement field-induced quantum phase transitions from the EQAH states to the Fermi liquid, FQAH liquid and the likely composite Fermi liquid. Our observations established a new topological phase of electrons with quantized Hall resistance at zero magnetic field and enriched the emergent quantum phenomena in materials with topological flat bands.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Additively-manufactured monocrystalline YBCO superconductor

Abstract Single-crystal microstructures enable high-performance YBa 2 Cu 3 O 7-x superconductors which are however limited to simple shapes due to their brittleness. Additive manufacturing can fabricate YBa 2 Cu 3 O 7-x superconductor with complex shapes, albeit with a polycrystalline microstructure. Here, we demonstrate a route to grow single-crystals from 3D-ink-printed, polycrystalline, sintered superconducting YBa 2 Cu 3 O 7-x (YBCO or Y123) + Y 2 BaCuO 5 (Y211), manufacturing objects with complex architectures displaying both high critical current density (J c =2.1 × 10 4 A . cm –2 , 77 K) and high critical temperature (T c = 88-89.5 K). An ink containing precursor powders (Y 2 O 3 , BaCO 3 , and CuO) is 3D-extruded into complex geometries and then reaction-sintered to obtain polycrystalline Y123 + Y211. A seed is then utilized to transform these 3D-printed parts from polycrystal to monocrystal via the melt growth method. The geometric details of 3D-printed parts survive the process without slumping, sagging or collapse, despite the long-term presence of liquid above the peritectic temperature. Origami structures can be created by sheet folding after 3D-printing. This additive approach enables the facile fabrication of superconducting devices with complex shapes and architectures, such as advanced undulator magnets to generate synchrotron radiation and microwave cavities for dark-matter axion search. This work highlights the potential of additive manufacturing for producing monocrystalline cuprate superconductors and opens the door to additive manufacturing of other monocrystalline functional ceramic or semiconductor materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Changes in rock matrix compressibility during deep CO 2 storage

Geologic carbon storage projects aim to permanently trap large volumes of CO 2 in reservoir rock sealed with low permeability layers. As high-pressure supercritical or liquid CO 2 is injected, hydromechanical and chemical processes caused by the reaction between the rock and acidic mixture of brine and CO 2 are initiated. The compressibility of reservoir rock needs to be properly characterized in order to assess the deformation and stability of the host formations, and there are a number of factors to be considered, including the matrix structure, solid, pores, and fluid. Here this study triggers from a fundamental question whether CO 2 treatment affects the compressibility of the rock matrix and its dominant composing solid minerals. Three different reservoir representatives are selected: Berea sandstone for silica-rich rock, and Apulian limestone and Indiana limestone for calcite-rich rock. Quartz and calcite are the main composing minerals of the reservoir rock, and their crystals are separately examined. Experimental methods are introduced for high-pressure CO 2 treatment of water-saturated materials, and measurements of the unjacketed and solid compressibilities are conducted. No change in the solid compressibility of the sandstone and quartz after CO 2 treatment is observed, whereas it increases by 18–21% for the limestones and by 15% for calcite. The latter observation is inconsistent with the ultrasonic wave velocities measurements and is believed to be due to the local dissolution of the calcite crystal surface. Our results show that only the solid matrix of the limestones becomes more compressible after CO 2 treatment. Consequent microimaging and mercury intrusion porosimetry analyses allowed observations of dissolution and precipitation of calcite, and creation of new connected and non-connected pores. Finally, the changes in limestone solid compressibilities and pore structure could significantly affect the rock properties and behavior during and after CO 2 injection and should be accounted for in the reservoir models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

In Situ Monitoring the Nucleation and Growth of Nanoscale CaCO 3 at the Oil–Water Interface

Interfaces can actively control the nucleation kinetics, orientations, and polymorphs of calcium carbonate (CaCO 3 ). Prior studies have revealed that CaCO 3 formation can be affected by the interplay between chemical functional moieties on solid–liquid or air–liquid interfaces as well as CaCO 3 ’s precursors and facets. Yet little is known about the roles of a liquid–liquid interface, specifically an oil–liquid interface, in directing CaCO 3 mineralization which are common in natural and engineered systems. Here, in this study, by using in situ X-ray scattering techniques to locate a meniscus formed between water and a representative oil, isooctane, we successfully monitored CaCO 3 formation at the pliable isooctane–water interface and systematically investigated the pivotal roles of the interface in the formation of CaCO 3 (i.e., particle size, its spatial distribution with respect to the interface, and its mineral phase). Different from bulk solution, ∼5 nm CaCO 3 nanoparticles form at the isooctane–water interface. They stably exist for a long time (36 h), which can result from interface-stabilized dehydrated prenucleation clusters of CaCO 3 . There is a clear tendency for enhanced amounts and faster crystallization of CaCO 3 at locations closer to isooctane, which is attributed to a higher pH and an easier dehydration environment created by the interface and oil. Our study provides insights into CaCO 3 nucleation at an oil–water interface, which can deepen our understanding of pliable interfaces interacting with CaCO 3 and benefit mineral scaling control during energy-related subsurface operation.

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

Assembly of small silica nanoparticles using lipid-tethered DNA ‘bonds’

Single-stranded DNA molecules modified with cholesterol functional groups are physically tethered to silica nanoparticles (diameter 25 nm) that are encapsulated in a lipid bilayer. Such tethering increases the azimuthal mobility of the DNA molecules across the nanoparticle surface and enables nonspecific bonding, eliminating the need for specialized surface chemistries (such as silane or thiol ligands). To induce assembly, double-stranded DNA ‘bridge’ molecules are then added with complementary nucleotides to the DNA ‘anchor’ molecules that are physically tethered to the lipids on the surface of the particles. Assembly is observed to occur at room temperature and without the need for temperature annealing. Using automated liquid handling tools, assemblies are created in high throughput and rapidly characterized using SAXS. It is determined that the relative concentration of DNA-to-silica and the ionic strength of the solution are important parameters that affect the resulting assembly. Analysis of SAXS data is performed using coarse-grained particle dynamics simulations. The results support the spontaneous formation of semi-crystalline particle assemblies by particle condensation, where the interparticle distance is tuned by the sequence of the DNA ‘bridge’ used to link the particles. Crystallinity analysis performed on the resulting simulations, optimized to match SAXS observations, suggest that particle clusters display increased crystallinity in the center of the clusters, but their maximum size remains relatively small (sub-micron) before settling occurs, which limits the extent of crystallization.

Chiang, Huat Thart [Univ. of Washington, Seattle, ↗
DOE OSTI2023

Beyond single tetrahedron physics of the breathing pyrochlore compound Ba 3 ⁢Yb 2 ⁢Zn 5 ⁢O 11

Recently, a new class of quantum magnets, the so-called breathing pyrochlore spin systems, have attracted much attention due to their potential to host exotic emergent phenomena. Here, we present magnetometry, heat capacity, thermal conductivity, muon-spin relaxation, and polarized inelastic neutron scattering measurements performed on high-quality single crystal samples of the breathing pyrochlore compound Ba 3 ⁢Yb 2 ⁢Zn 5 ⁢O 11 . We interpret these results using a simplified toy model and provide insight into the low-energy physics of this system beyond the single tetrahedron physics proposed previously.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase equilibria in iron-rich Sm–Fe–Ti and Sm–(Fe,Co)–Ti alloys at 1100–1200 °C

Iron-rich corners of ternary Sm–Fe–Ti phase diagram at 1100 °C and quasi-ternary Sm–Fe 0.8 Co 0.2 –Ti phase diagrams at 1100 and 1200 °C are constructed based on experimental investigation of equilibrated alloys with electron probe microanalysis, X-ray diffraction and thermomagnetic analysis. In addition, the upper boundaries of the temperature ranges of Sm(Fe,[Co,]Ti) 12 and Sm 3 (Fe,[Co,]Ti) 29 phases are determined with differential thermal analysis to update earlier rough estimates. The existence of a high-temperature phase of the Th 2 Ni 17 type, originally reported by Ivanova et al. [J. Alloys Compd. 224 (1995) 29], is confirmed. In the Sm–Fe 0.8 Co 0.2 –Ti system, the composition and equilibria of this hexagonal phase are established for 1200 °C; it is Sm-depleted (≈9.8 at.% Sm) compared to the 2:17 stoichiometry and it coexists with the rhombohedral 2:17 phase. The magnetic anisotropy of the cobalt-substituted Th 2 Ni 17 -type phase is planar, with the easy magnetization direction parallel to [100]. Equilibrium between a Sm-rich liquid phase and the 1:12 phase, which is important for the development of new high-performance permanent magnets, is absent up to 1000 °C, but does exist at 1100 °C (for the 1:12 phase with at least 8.7–8.9 at.% Ti) and at 1200 °C (for the 1:12 phase with as little as 7.4 at.% Ti). The development of magnets may be complicated, however, by an observed tendency of the high-temperature liquid to solidify into ferromagnetic phases including the Th 2 Ni 17 -type phase. The Curie temperatures of the α-(Fe,Ti), Sm(Fe,Ti) 12 , Sm 3 (Fe,Ti) 29 and rhombohedral Sm 2 (Fe,Ti) 17 phases are not only increased by the partial Co substitution for Fe, but their dependence on the Ti concentration is changed by this Co substitution from positive (or, for the 1:12 phase, zero) to negative values.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast relaxation of lattice distortion in two-dimensional perovskites

Direct visualization of ultrafast coupling between charge carriers and lattice degrees of freedom in photoexcited semiconductors has remained a long-standing challenge and is critical for understanding the light-induced physical behaviour of materials under extreme non-equilibrium conditions. Here we obtain a direct visualization of the structural dynamics in monocrystalline 2D perovskites. We achieve this by monitoring the evolution of wavevector-resolved ultrafast electron diffraction intensity following above-bandgap high-density photoexcitation. Our analysis reveals a light-induced ultrafast reduction in antiferro-distortion resulting from a strong interaction between the electron–hole plasma and perovskite lattice, which induces an in-plane octahedra rotation towards a more symmetric phase. Correlated ultrafast spectroscopy performed at the same carrier density as ultrafast electron diffraction reveals that the creation of a dense electron–hole plasma triggers the relaxation of lattice distortion at shorter timescales by modulating the crystal cohesive energy. Then we show that the interaction between carrier gas and lattice can be altered by tailoring the rigidity of the 2D perovskite by choosing an appropriate organic spacer layer.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2022

Thermoplasticity of metallic glasses: Processing and applications

Owing to their extremely low shear viscosity and fragility, bulk metallic glasses (BMGs) are envisioned as new net-shaping material replacing conventionally used crystalline metals for various applications. The first part of this review describes the general technique and kinetics of thermoplastic forming. The second part elaborates on the thermal processing capability of metallic glasses from atomic- to meter-scale. The micro-/nano-structure obtained after various quenching rates followed by temperature, time and heating rate dependent thermoplastic forming is described in the third section. The deformation behavior and flow kinetics of BMGs in terms of composition, mold features, and applied conditions are elaborated in the fourth section. The variation of volume and enthalpy at the supercooled liquid region, strength retention of metallic glasses compared to other conventional metals and alloys, and kinetics of liquid fragility are given in the fifth section. Finite element modeling and molecular dynamics simulations of high-temperature deformation in BMGs are presented in the sixth part. Finally, thermally-formed BMGs used for different applications, including energy, biomedical and micro-optics, are presented in the final part. Altogether, this review provides an overview of shaping capabilities and modifications in the macro-scale properties and short-to-medium range order of BMGs upon thermoplastic forming.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reverse order-disorder transition of Janus particles confined in two dimensions

In this study Janus particles with different patch sizes, confined to two dimensions, generate a series of patterns of interest to the field of nanoscience. Here we observe reverse melting, where for some densities the system melts under cooling. For a broad range of hydrophobic patch sizes ( 60° < θ 0 < 90°), a reentrant transition from solid to liquid and then to an ordered phase emerges as temperature (T) decreases due to the formation of rhombus chains at low T. This reentrant phase has pseudo long-range orientational order but short-range translational order, similar to a hexatic phase. Our work provides guidelines to study the melting and assembly of Janus particles in two dimensions, as well as mechanisms to generate phases with specific symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Orientation-dependent stabilization of MgCr 2 O 4 spinel thin films

AB 2 O 4 normal spinels with a magnetic B site can host a variety of magnetic and orbital frustrations leading to spin-liquid phases and field-induced phase transitions. Here, we report the epitaxial growth of (111)-oriented MgCr 2 O 4 thin films. By characterizing the structural and electronic properties of films grown along the (001) and (111) directions, the influence of growth orientation has been studied. Despite distinctly different growth modes observed during deposition, the comprehensive characterization reveals no measurable disorder in the cation distribution nor multivalency issue for Cr ions in either orientation. Contrary to a naive expectation, the (111) stabilized films exhibit a smoother surface and a higher degree of crystallinity than (001)-oriented films. The preference in growth orientation is explained within the framework of heteroepitaxial stabilization in connection to a significantly lower (111) surface energy. Furthermore, these findings open broad opportunities in the fabrication of two-dimensional kagome-triangular heterostructures with emergent magnetic behavior inaccessible in bulk crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗