Search NASA⌕ Search

SEARCH · Search NASA

Results for “crystal cell parameters”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Possible Radiation-Induced Transition from Monazite-(Ce) to Xenotime-(Y)

This study examines two pegmatitic monazite samples (2a and 4b, these numbers are related to a previous study) to determine their crystal chemistry and effects of internal radiation damage using synchrotron high-resolution powder X-ray diffraction and electron-probe micro-analysis. Both the huttonite and cheralite substitutions are discussed. Rietveld structure refinement of sample 2a shows three different phases [2a = monazite-(Ce), 2b = monazite-(Ce), and 2c = xenotime-(Y)] with distinct structural parameters. The changes among the unit-cell parameters between the two monazite-(Ce) phases is more pronounced in the a followed by the b and c unit-cell parameters. Sample 4a is a single-phase monazite-(Sm) that contains 0.164 apfu Th. Phase 2c with space group I41/amd arises from redistribution of La, Ce, Pr, Nd, Sm, Gd, Dy, Si, and Y atoms from those in monazite (space group P21/n). A possible cause for the phase transition from monazite-(Ce) to xenotime-(Y) is α-radiation events over a long geological time. However, other chemical processes cannot be ruled out as a cause for the transition.

structural variations↗
NASA NTRS1996

Crystallography of the NiHfSi Phase in a NiAl (0.5 Hf) Single-Crystal Alloy

Small additions of Hf to conventionally processed NiAl single crystals result in the precipitation of a high density of cuboidal G-phase along with a newly identified silicide phase. Both of these phases form in the presence of Si which is not an intentional alloying addition but is a contaminant resulting from contact with the ceramic shell molds during directional solidification of the single-crystal ingots. The morphology, crystal structure and Orientation Relationship (OR) of the silicide phase in a NiAl (0.5 at.%Hf) single-crystal alloy have been determined using transmission electron microscopy, electron microdiffraction and energy dispersive X-ray spectroscopy. Qualitative elemental analysis and indexing of the electron microdiffraction patterns from the new phase indicate that it is an orthorhombic NiHfSi phase with unit cell parameters, a = 0.639 nm, b = 0.389 nm and c = 0.72 nm, and space group Pnma. The NiHfSi phase forms as thin rectangular plates on NiAl/111/ planes with an OR that is given by NiHfSi(100))(parallel) NiAl(111) and NiHfSi zone axes(010) (parallel) NiAl zone axes (101). Twelve variants of the NiHfSi phase were observed in the alloy and the number of variants and rectangular morphology of NiHfSi plates are consistent with symmetry requirements. Quenching experiments indicate that nucleation of the NiHfSi phase in NiAI(Hf) alloys is aided by the formation of NiAl group of zone axes (111) vacancy loops that form on the NiAl /111/ planes.

Garg, A.↗
NASA NTRS2013

Crystal-Chemical Analysis of Soil at Rocknest, Gale Crater

The CheMin instrument on the Mars Science Laboratory rover Curiosity performed X-ray diffraction analysis on Martian soil [1] at Rocknest in Gale Crater. In particular, crystalline phases from scoop 5 were identified and analyzed with the Rietveld method [2]. Refined unit-cell parameters are reported in Table 1. Comparing these unit-cell parameters with those in the literature provides an estimate of the chemical composition of the crystalline phases. For instance, Fig. 1 shows the Mg-content of Fa-Fo olivine as a function of the b unit-cell parameter using literature data. Our refined b parameter is indicated by the black triangle.

Morrison, S. M.↗
DOE OSTI2023

Conversion efficiency improvement of ELO GaAs solar cell, deposited on water soluble sacrificial buffer

Here, we demonstrate the improvement of power conversion efficiency of an epitaxially-lifted-off single junction GaAs solar cell deposited on a water-soluble sacrificial buffer architecture by introducing an additional germanium (Ge) interlayer between GaAs and the fluoride buffer. The epitaxial lift-off (ELO) technique has been extensively used to separate III-V device layers from their single crystal GaAs substrates. However, conventional ELO requires the use of concentrated hydrofluoric acid (HF) for extended times to etch out the sacrificial layer, subsequently degrading the surface roughness of the parent wafer. As a result, the wafer has to undergo expensive and intensive chemical mechanical polishing (CMP) processes, costing about 25% of a pristine 6-inch GaAs substrate. In our previous work, we demonstrated a method to eliminating the need for CMP post-processing by using water-assisted ELO (H 2 O-ELO). A water-soluble, 3-layer buffer architecture was developed using alkaline earth compounds. However, devices suffered low performance due to a high defect density in the GaAs active layers. A Ge interlayer was introduced to provide a more favorable surface energy for GaAs growth, which leads to an improvement of GaAs crystal quality. The Ge deposition conditions were optimized to achieve a high-quality Ge layer on triple-layer fluoride buffer. Single junction GaAs devices fabricated on Ge/(Ca,Sr)F 2 /BaF 2 /(Ca,Sr)F 2 showed improvement of solar cell performance parameters featuring increases in V oc by 23.7%, and fill factor (F.F.) by 4.9%, which results in an overall improvement of power conversion efficiency from 10.3% to 12.69%.

14 SOLAR ENERGY↗
DOE OSTI2020

Structural Assessment of Fluorine, Chlorine, Bromine, Iodine, and Hydroxide Substitutions in Lead Arsenate Apatites (Mimetites)–Pb 5 (AsO 4 ) 3 X

Five lead-arsenate apatites (mimetites)-Pb 5 (AsO 4 ) 3 X—where X denotes fluorine (F), chlorine (Cl), bromine (Br), iodine (I), and hydroxide (OH)—were synthesized via precipitation from aqueous solutions. The crystal structures were determined through Rietveld refinement of powder synchrotron X-ray data. All the compounds crystallized in the hexagonal class symmetry (space group P6 3 /m). The Rietveld refinement indicated that mimetite-Cl, -Br, -I, and -OH had an anion deficiency at position X. Substitution of halogens in a mimetite structure brought about systematic changes in unit-cell parameters, interatomic distances, and metaprism twist angles φ, proportional to the substituted halogen’s ionic radius. Mimetite-OH did not follow the linear correlations determined within the series. Twist angle φ, a useful device for monitoring changes in apatite topology, ranged from 20.34° for mimetite-F to 11.42° for mimetite-I. The geometric method has been proposed for determining the diameter of hexagonal channels hosting halogens in apatites. A comparison of the results with halogenated pyromorphites showed similar systematic trends: the substitutions in mimetites have comparable effect on the interatomic distances as in their phosphorous analogues.

58 GEOSCIENCES↗
DOE OSTI2024

Crystal structure and synchrotron X-ray powder reference pattern for the porous pillared cyanonickelate, Ni(3-amino-4,4′-bipyridine)[Ni(CN) 4 ]

The structure of Ni(3-amino-4,4′-bipyridine)[Ni(CN) 4 ] (or known as Ni-BpyNH 2 ) in powder form was determined using synchrotron X-ray diffraction and refined using the Rietveld refinement technique (R = 8.8%). The orthorhombic (Cmca) cell parameters were determined to be a = 14.7218(3) Å, b = 22.6615(3) Å, c = 12.3833(3) Å, V = 4131.29(9) Å 3 , and Z = 8. Ni-BpyNH 2 forms a 3-D network, with a 2-D Ni(CN) 4 net connecting to each other via the BpyNH 2 ligands. Further, there are two independent Ni sites on the net. The 2-D nets are connected to each other via the bonding of the pyridine “N” atom to Ni2. The Ni2 site is of six-fold coordination to N with relatively long Ni2–N distances (average of 2.118 Å) as compared to the four-fold coordinated Ni1–C distances (average of 1.850 Å). The Ni(CN) 4 net is arranged in a wave-like fashion. The functional group, –NH 2 , is disordered and was found to be in the m-position relative to the N atom of the pyridine ring. Instead of having a unique position, N has ¼ site occupancy in each of the four m-positions. The powder reference diffraction pattern for Ni-BpyNH 2 was prepared and submitted to the Powder Diffraction File (PDF) at the International Centre of Diffraction Data (ICDD).

36 MATERIALS SCIENCE↗
DOE OSTI2022

To grind or not to grind? The influence of mechanical and thermal treatments on the $\mathrm{Cu/Zn}$ disorder in $\mathrm{Cu}$ 2 $\mathrm{ZnSn}$($\mathrm{S}$ x $\mathrm{Se}$ 1-x ) 4 monograins

We report that Kesterite-type based thin film solar cell technologies are mainly based on polycrystalline absorber layers. A promising low cost alternative technology uses Cu 2 ZnSn(S x Se 1-x ) 4 (CZTSSe) monograins (single crystals of 20–100 μm size) which are fixed in a polymer matrix to form a flexible solar cell. The Cu/Zn disorder is discussed as a possible reason for band tailing causing voltage losses limiting the efficiency of CZTSSe-based devices. The experimental determination of the order parameter Q which is a quantitative measure of Cu/Zn disorder, requires a differentiation between the isoelectronic cations Cu + and Zn 2+ . An in-depth analysis of neutron diffraction data allows the determination of type and concentration of intrinsic point defects including a distinction between Cu and Zn. Neutron diffraction requires large sample volumes, thus monograins offer the unique possibility to correlate structural disorder in CZTSSe with device performance parameters. In this study we tackle the influence of grinding the monograins on stoichiometry deviations, the Cu/Zn disorder as well as intrinsic point defects and optoelectronic properties of CZTSSe monograins. Moreover, an easy methodology based on Raman scattering spectroscopy is proposed for the assessment of Cu/Zn disorder in the CZTSSe compounds.

14 SOLAR ENERGY↗
DOE OSTI2025

On the Structural Origin of Fast Li-Ion Cycling in Tetragonal Bronze-Type Nb 8 W 9 O 47

Bronze and bronze-derived, mixed metal oxides are materials of growing interest for lithium-ion battery anodes due to their high-rate capabilities and the potential for high energy densities via multielectron redox. We report on the synthesis and electrochemical properties of the tetragonal tungsten bronzetype phase Nb 8 W 9 O 47 and our investigation into the structural evolution of this phase upon lithium (de)intercalation, with particular attention to how it relates to the fast-cycling capabilities. Electrochemical cycling shows Nb 8 W 9 O 47 can achieve greater than one Li + per transition metal at rates of C/2 or slower and maintains a capacity equivalent to 0.65 Li + per transition metal at a rate of 20C. Sequential Rietveld analysis of operando X-ray diffraction data reveals anisotropic changes to the unit cell parameters during discharge, accompanied by an off-centering of the transition metals within their octahedral environments. While higher capacities can be accessed below 1.2 V, cells discharged to 1.0 V exhibit a significant expansion of the unit cell volume and reduced Coulombic efficiencies. This report suggests that rigid crystal structures incapable of undergoing polyhedral rotational deformations can instead exploit displacive distortions of the transition metal redox center to enable fast Li-ion cycling with minimal atomic rearrangements.

25 ENERGY STORAGE↗
DOE OSTI2021

Optimizing Vertical Crystallization for Efficient Perovskite Solar Cells by Buried Composite Layers

Planar-heterojunction perovskite solar cells (PSCs) have experienced rapid evolution in recent years because of the low-temperature processing, suitable alignment, and high mobility of the tin oxide buried contact layer. However, improper SnO 2 surface states and poor crystallinity of the top perovskite films are still the main obstacles for the planar PSCs in which performance always lags behind their mesoporous counterparts. Herein, a new buried contact is reported by introducing graphitic carbon nitride (g-C 3 N 4 ) into the commonly used SnO 2 which performs outstanding transmittance, conductivity, and surface states for a high-quality electron-transporting layer. Moreover, the vertical composition and crystallinity of the top perovskite film are manipulated by rich amino groups on the edge of the g-C 3 N 4 nanosheets which induce the prenucleation of the lead-rich species at the buried interface. Benefiting from the high-quality buried contacts and the optimized perovskite layers, the resultant PSCs achieve a champion efficiency of 21.5% with all photovoltaic parameters enhanced in comparison with their control counterparts (<20%).

planar-heterojunction perovskite solar cells↗
DOE OSTI2023

Physical structure of constitutional isomers influences antiproliferation activity of thiosemicarbazone-alkylthiocarbamate copper complexes

A series of hybrid thiosemicarbazone-alkylthiocarbamate copper complexes with similar electronic environments but distinct physical structures have been prepared, characterized, and evaluated for antiproliferation activity. The complexes include the constitutional isomers (1-phenylpropane-1-imine-(O-ethylthiocarbamato)-2-one-(N-methylthiosemicarbazonato))copper(II) (CuL 1 ) and (1-phenylpropane-1-one-(N-methylthiosemicarbazonato)-2-imine-(O-ethylthiocarbamato))copper(II) (CuL 2 ) along with (1-propane-1-imine-(O-ethylthiocarbamato)-2-one-(N-methylthiosemicarbazonato))copper(II) (CuL 3 ). Complexes CuL 1 and CuL 2 differ in the positions of the pendent thiosemicarbazone (TSC) and alkylthiocarbamate (ATC) moieties on the 1-phenylpropane backbone. Complex CuL 3 employs a propane backbone with the TSC in the 2-position as in CuL 1 . The isomer pair CuL 1 and CuL 2 have equivalent electronic environments with indistinguishable Cu II/I potentials (E 1/2 = -0.86 V vs. ferrocenium/ferrocene) and electron paramagnetic resonance (EPR) spectra (g ∥ = 2.26, g ⊥ = 2.08). The electronic structure of CuL 3 has a similar E 1/2 of -0.84 V and identical EPR parameters to CuL 1, 2 . Single crystal X-ray diffraction studies confirm a consistent donor environment with no substantial variation in the Cusingle bondN or Cusingle bondS bond distances and angles between the complexes. The antiproliferation activities of the CuL 1–3 were evaluated against the lung adenocarcinoma cell line (A549) and nonmalignant lung fibroblast cell line (IMR-90) using the MTT assay. CuL 1 had the highest A549 activity ( A549 EC 50 = 0.065 μM) and selectivity ( IMR-90 EC 50 / A549 EC 50 = 20). The constitutional isomer CuL 2 displayed decreased A549 activity (0.18 μM) and selectivity (10.6). The complex CuL 3 displayed activity (0.009 μM) similar to CuL 1 but with a lack of selectivity (1.0). Cellular copper loading determined by ICP-MS was consistent with the activity and selectivity trends. The complexes CuL 1–3 did not induce reactive oxygen species (ROS) generation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

An overview of vertebrate mineralization with emphasis on collagen-mineral interaction

The nucleation, growth, and development of mineral crystals through their interaction principally with collagen in normal bone and calcifying tendon have been elaborated by applying a number of different techniques for analysis of the inorganic and organic constituents of these tissues. The methods have included conventional and high voltage electron microscopy, electron diffraction, microscopic tomography and 3D image reconstruction, and atomic force microscopy. This summary presents results of these studies that have now characterized the size, shape, and aspects of the chemical nature of the crystals as well as their orientation, alignment, location, and distribution with respect to collagen. These data have provided the means for understanding more completely the formation and strength of the collagen-mineral composite present in most vertebrate calcifying tissues and, from that information, a basis for the adaptation of such tissues under mechanical constraints. In the context of the latter point, other data are given showing effects on collagen in bone cell cultures subjected to the unloading parameters of spaceflight. Implications of these results may be particularly relevant to explaining loss of bone by humans and other vertebrate animals during missions in space, during situations of extended fracture healing, long-term bedrest, physical immobilization, and related conditions. In a broader sense, the data speak to the response of bone and mineralized vertebrate tissues to changes in gravitational loading and applied mechanical forces in general.

Non-NASA Center↗
NASA NTRS1994

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗
DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

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

Effects of Physical and Chemical Pressure on Charge Density Wave Transitions in LaAg1−xAuxSb2 Single Crystals

The structural characterization and electrical transport measurements at ambient and applied pressures of the compounds of the LaAg1−xAuxSb2 family are presented. Up to two charge density wave (CDW) transitions could be detected upon cooling from room temperature and an equivalence of the effects of chemical and physical pressure on the CDW ordering temperatures was observed with the unit cell volume being a salient structural parameter. As such LaAg1−xAuxSb2 is a rare example of a non-cubic system that exhibits good agreement between the effects of applied, physical, pressure and changes in unit cell volume from steric changes induced by isovalent substitution. Additionally, for LaAg0.54Au0.46Sb2 anomalies in low temperature electrical transport were observed in the pressure range where the lower charge density wave is completely suppressed.

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

Exploring cycling induced crystallographic change in NMC with X-ray diffraction computed tomography

This study presents the application of X-ray diffraction computed tomography for the first time to analyze the crystal dimensions of LiNi 0.33 Mn 0.33 Co 0.33 O 2 electrodes cycled to 4.2 and 4.7 V in full cells with graphite as negative electrodes at 1 μm spatial resolution to determine the change in unit cell dimensions as a result of electrochemical cycling. The nature of the technique permits the spatial localization of the diffraction information in 3D and mapping of heterogeneities from the electrode to the particle level. An overall decrease of 0.4% and 0.6% was observed for the unit cell volume after 100 cycles for the electrodes cycled to 4.2 and 4.7 V. Additionally, focused ion beam-scanning electron microscope cross-sections indicate extensive particle cracking as a function of upper cut-off voltage, further confirming that severe cycling stresses exacerbate degradation. Finally, the technique facilitates the detection of parts of the electrode that have inhomogeneous lattice parameters that deviate from the bulk of the sample, further highlighting the effectiveness of the technique as a diagnostic tool, bridging the gap between crystal structure and electrochemical performance.

25 ENERGY STORAGE↗
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 OSTI2025

Structure of ice VII with Hirshfeld atom refinement

In the refinement of the crystal structures of ice, the best results obtained so far have been with neutron diffraction because the most troublemaking aspects are the hydrogen atoms. In nine out of twenty ice structures, the hydrogen atoms are disordered, which makes proper refinement more difficult. In our previous paper describing the structure of ice VI we proved that, using Hirshfeld atom refinement (HAR) based on synchrotron X-ray data, it is possible to obtain results comparable with those from neutron experiments. In this work, we investigate another structure of high-pressure disordered ice, cubo-ice (ice VII). Single crystals of cubo-ice were grown under pressure in diamond anvil cells. X-ray diffraction measurements were conducted at a synchrotron source facility (APS, University of Chicago, USA) as well as on our regular in-house laboratory diffractometer with Ag radiation. The data collected were further refined with HAR. Comparison of the structural parameters obtained with those derived from neutron diffraction showed very good agreement in terms of bond lengths and fairly good agreement in terms of hydrogen atom ADPs. We were also able to perform unconstrained refinements with various split-atom models.

Chemistry↗