Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “developable planes”

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 271 records · Page 15
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 OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Geochemical and geomechanical alteration of siliciclastic reservoir rock by supercritical CO 2 -saturated brine formed during geological carbon sequestration

Geologic carbon sequestration (GCS) is a method for storing CO 2 and mitigating greenhouse gas emissions. During GCS, carbon dioxide dissolves into pore water, resulting in a low-pH brine that can react with reservoir rock minerals. This work evaluates the effects of geochemical reactions on geomechanical integrity of representative siliciclastic reservoir samples obtained from the Mt. Simon formation. Rock samples were aged 4 or 8 weeks in CO 2 -saturated brine under reservoir conditions, and in N 2 -saturated brine as a control. Post-aging, CT scans revealed more extensive micro-fracture development along horizontal bedding planes at grain edges in CO 2 versus N 2 -aged samples. Digital analysis of CO 2 -aged samples showed porosity increase from 8.1% to 15.8%. Scanning electron microscopy unearthed the loss of clay cementation, greater exposure of quartz and K-feldspar grains, and apparent surface roughening (confirmed by laser profilometry) in CO 2 -aged samples, but not in N 2 -aged samples. Fracture toughness as evaluated by scratch testing was reduced by 32.1% after 4 weeks in scCO 2 -saturated brine and 69.5% after 8 weeks. The main reason for weakening appears to be detachment of clays from quartz and feldspar grain surfaces, resulting in weakening of the rock matrix. Rock weakening may alter the geomechanical stability of storage formations.

58 GEOSCIENCES↗
DOE OSTI2016

The interface character distribution of cold-rolled and annealed duplex stainless steel

The interface character distributions (ICDs) of cold-rolled and annealed duplex stainless steel specimens, previously solid-solution-heated (SHT) at 1573 K and 1323 K, are investigated using electron backscatter diffraction (EBSD) and five-parameter analysis (FPA). For the δ-ferrite phase, high concentrations of low angle grain boundaries (LAGBs) are developed, and the boundary planes are predominantly oriented on (111). High angle grain boundaries (HAGBs) with misorientations ranging from 50{sup o} to 60° are mostly pure tilt boundaries, and the boundary planes are primarily located on (112). For the austenite phase, very high density of coherent twin boundaries are introduced, implying that the boundary planes are exactly oriented on (111). The phase boundary character distribution (PBPD) appears to be connected with the K-S and N-W orientation relationships (ORs) terminating on (110) {sub F}‖(111) {sub A} and (110){sub F}‖(hkl){sub A}, respectively. - Highlights: •Five-parameter analysis was used to determine interface character distribution. •ICD results of statistical meaning were achieved. •Initial microstructure has a significant effect on ICD. •Low and high angle GBs in δ phase prefer in (111) and (112) planes, respectively. •K-S phase boundaries terminate on (110){sub F}‖(111){sub A} while N-W ones on (110){sub F}‖(hkl){sub A}.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cold air quench control of local crystallization environment in fully air-processed carbon-based perovskite solar cells

Carbon-based perovskite solar cells (C-PSCs) present a low-cost route to efficient photovoltaic technology. Gas quenching is an essential process commonly used for solvent extraction in scalable module fabrication using solution processes, but film quality remains highly sensitive to the localized environmental processing temperature, especially during intermediate phase transitions before perovskite film high temperature annealing. In this study, we introduce a cold air quench strategy to precisely control the local crystallization temperature that simultaneously promotes uniform crystallization and facilitates strain relaxation in fully air-processed C-PSCs. It is found that implementing dry air at 10 °C as the quenching gas yields a champion power conversion efficiency (PCE) of 20.52 %, with preferential (110) orientation and reduced, homogenized in-plane strain. This newly developed technique of a cost-effective low temperature air quenching process also enhances carrier lifetime, reduces interface recombination, and improves charge extraction. Furthermore, this work presents cold air quenching as a scalable, economic method to control perovskite crystallization and strain in C-PSC fabrication, advancing the industrial viability of high-performance perovskite modules.

14 SOLAR ENERGY↗
DOE OSTI2018

Small-Angle Twist Grain Boundaries as Sinks for Point Defects

It is known that grain boundaries (GBs) provide sinks for defects induced into a solid by irradiation. At the same time radiation can change the atomic structure and chemistry of GBs, which in turn impacts the ability of GBs to continue absorbing defects. Although a number of studies have been reported for tilt GBs acting as defect sinks, the questions of how twist GBs evolve to absorb non-equilibrium concentrations of defects and whether GBs saturate as defect sinks for typical irradiation conditions have remained largely unanswered. Here, we used a combination of molecular dynamics and grand canonical Monte Carlo simulations to determine how twist GBs accommodate point defects. We used SiC and {001} and {111} twist GBs as model systems. We found that diffusion of defects along GBs in this material is slow and for most experimentally relevant conditions point defects will accumulate at twist GBs, driving structural and chemical evolution of these interfaces. During irradiation, screw dislocations within GB planes absorb interstitials by developing mixed dislocation segments that climb. Formation of mixed dislocations occurs either by nucleation of interstitial loops or by faulting/unfaulting of stacking faults. Furthermore both types of twist GBs can accommodate a high density of interstitials without losing the crystalline structure, irrespectively of the interstitial flux.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Equilibration of the planar modes of ultracold two-dimensional ion crystals in a Penning trap

In this work, planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane motion of the ions into two branches: high-frequency cyclotron modes dominated by kinetic energy and low-frequency E × B modes dominated by potential energy associated with thermal position displacements. Using an eigenmode analysis we extract the equilibration rate between these two branches as a function of the ratio of the frequencies that characterize the two branches and observe this equilibration rate to be exponentially suppressed as the ratio increases. Under experimental conditions relevant for current work at NIST, the predicted equilibration time is orders of magnitude longer than any relevant experimental timescales. We also study the coupling rate dependence on the initial energy of the planar modes and the number of ions. In addition, we show how increasing the rotating wall strength improves crystal stability. These details of in-plane mode dynamics help set the stage for developing strategies to efficiently cool the in-plane modes and improve the performance of single plane ion crystals for quantum information processing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

La 139 NMR investigation of the charge and spin order in a La 1.885 Sr 0.115 CuO 4 single crystal

139 La NMR is suited for investigations into magnetic properties of La 2 CuO 4 -based cuprates in the vicinity of their magnetic instabilities, owing to the modest hyperfine interactions between 139 La nuclear spins and Cu electron spins. We report comprehensive 139 La NMR measurements on a single-crystal sample of high-T c superconductor La 1.885 Sr 0.115 CuO 4 in a broad temperature range across the charge and spin order transitions (T charge ≃ 80 K, T neutron spin ≃ T c = 30 K). From the high-precision measurements of the linewidth for the nuclear spin I z = + 1 / 2 to -1/2 central transition, we show that paramagnetic line broadening sets in precisely at T charge due to enhanced spin correlations within the CuO 2 planes. Additional paramagnetic line broadening ensues below ~35 K, signaling that Cu spins in some segments of CuO 2 planes are on the verge of three-dimensional magnetic order. A static hyperfine magnetic field arising from ordered Cu moments along the ab plane, however, begins to develop only below T$μSR\atop{spin}$ = 15 – 20 K, where earlier muon spin rotation measurements detected Larmor precession for a small volume fraction (~20 % ) of the sample. Based on the measurement of 139 La nuclear-spin-lattice relaxation rate 1/T 1 , we also show that charge order triggers enhancement of low-frequency Cu spin fluctuations inhomogeneously; a growing fraction of 139 La sites is affected by enhanced low-frequency spin fluctuations toward the eventual magnetic order, whereas a diminishing fraction continues to exhibit a behavior analogous to the optimally superconducting phase even below T charge . In conclusion, these 139 La NMR results corroborate our recent 63 Cu NMR observation that a very broad, anomalous winglike signal gradually emerges below T charge , whereas the normally behaving, narrower main peak is gradually wiped out [T. Imai et al., Phys. Rev. B 96, 224508 (2017)]. Furthermore, we show that the enhancement of low-energy spin excitations in the low-temperature regime below Tneutron spin (≃ Tc) depends strongly on the magnitude and orientation of the applied magnetic field.

42 ENGINEERING↗
DOE OSTI2016

High-energy tail of the Galactic Center gamma-ray excess

Observations by the Fermi-LAT have uncovered a bright, spherically symmetric excess surrounding the center of the Milky Way Galaxy. The spectrum of the ${\gamma}$-ray excess peaks sharply at an energy ${\sim}2\text{ }\text{ }\mathrm{GeV}$, exhibiting a hard spectrum at lower energies, and falls off quickly above an energy ${\sim}5\text{ }\text{ }\mathrm{GeV}$. The spectrum of the excess above ${\sim}10\text{ }\text{ }\mathrm{GeV}$ is potentially an important discriminator between different physical models for its origin. We focus on observations of the ${\gamma}$-ray excess at energies exceeding 10 GeV, finding: (1) a statistically significant excess remains in the energy range 9.5–47.5 GeV, which is not degenerate with known diffuse emission templates such as the Fermi bubbles, (2) the radial profile of the excess at high energies remains relatively consistent with data near the spectral peak (3) the data above ${\sim}5\text{ }\text{ }\mathrm{GeV}$ prefer a slightly greater ellipticity with a major axis oriented perpendicular to the Galactic plane. Using the recently developed non-Poissonian template fit, we find mild evidence for a point-source origin for the high-energy excess, although given the statistical and systematic uncertainties we show that a smooth origin of the high-energy emission cannot be ruled out. We discuss the implication of these findings for pulsar and dark matter models of the ${\gamma}$-ray excess. Finally we provide a number of updated measurements of the ${\gamma}$-ray excess, utilizing novel diffuse templates and the Pass 8 data set.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

WE-AB-BRA-11: Improved Imaging of Permanent Prostate Brachytherapy Seed Implants by Combining an Endorectal X-Ray Sensor with a CT Scanner

Purpose: To test feasibility of the use of a digital endorectal x-ray sensor for improved image resolution of permanent brachytherapy seed implants compared to conventional CT. Methods: Two phantoms simulating the male pelvic region were used to test the capabilities of a digital endorectal x-ray sensor for imaging permanent brachytherapy seed implants. Phantom 1 was constructed from acrylic plastic with cavities milled in the locations of the prostate and the rectum. The prostate cavity was filled a Styrofoam plug implanted with 10 training seeds. Phantom 2 was constructed from tissue-equivalent gelatins and contained a prostate phantom implanted with 18 strands of training seeds. For both phantoms, an intraoral digital dental x-ray sensor was placed in the rectum within 2 cm of the seed implants. Scout scans were taken of the phantoms over a limited arc angle using a CT scanner (80 kV, 120–200 mA). The dental sensor was removed from the phantoms and normal helical CT and scout (0 degree) scans using typical parameters for pelvic CT (120 kV, auto-mA) were collected. A shift-and add tomosynthesis algorithm was developed to localize seed plane location normal to detector face. Results: The endorectal sensor produced images with improved resolution compared to CT scans. Seed clusters and individual seed geometry were more discernable using the endorectal sensor. Seed 3D locations, including seeds that were not located in every projection image, were discernable using the shift and add algorithm. Conclusion: This work shows that digital endorectal x-ray sensors are a feasible method for improving imaging of permanent brachytherapy seed implants. Future work will consist of optimizing the tomosynthesis technique to produce higher resolution, lower dose images of 1) permanent brachytherapy seed implants for post-implant dosimetry and 2) fine anatomic details for imaging and managing prostatic disease compared to CT images. Funding: LSU Faculty Start-up Funding. Disclosure: XDR Radiography has loaned our research group the digital x-ray detector used in this work. CoI: None.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Peculiarities of the Self-Action of Inclined Wave Beams Incident on a Discrete System of Optical Fibers

Based on a discrete nonlinear Schrödinger equation (DNSE), we studied analytically and numerically the peculiarities of the self-action of one-dimensional quasi-optic wave beams injected into a spatially inhomogeneous medium consisting of a set of equidistant mutually coupled optical fibers. A variational approach allowing the prediction of the global evolution of localized fields with the initially plane phase front was developed. The self-consistent equations are obtained for the main parameters of such beams (the position of the center of mass, the effective width, and linear and quadratic phase-front corrections) in the aberrationless approximation. The case of radiation incident on a periodic system of nonlinear optical fibers at an angle to the axis oriented along them is analyzed in detail. It is shown that for the radiation power exceeding a critical value, the self-focusing of the wave field is observed, which is accompanied by the shift of the intensity maximum followed by the concentration of the main part of radiation only in one of the structural elements of the array under study. In this case, the beams propagate along paths considerably different from linear and the direction of their propagation changes compared to the initial direction. Asymptotic expressions are found that allow us to estimate the self-focusing length and to determine quite accurately the final position of a point with the maximum field amplitude after radiation trapping a channel. The results of the qualitative study of the possible self-channeling regimes for wave beams in a system of weakly coupled optical fibers in the aberrationless approximation are compared with the results of direct numerical simulations within the DNSE framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Measurement of the electric potential and the magnetic field in the shifted analysing plane of the KATRIN experiment

The projected sensitivity of the effective electron neutrino-mass measurement with the KATRIN experiment is below 0.3 eV (90 % CL) after 5 years of data acquisition. The sensitivity is affected by the increased rate of the background electrons from KATRIN’s main spectrometer. A special shifted-analysing-plane (SAP) configuration was developed to reduce this background by a factor of two. The complex layout of electromagnetic fields in the SAP configuration requires a robust method of estimating these fields. We present in this paper a dedicated calibration measurement of the fields using conversion electrons of gaseous 83m Kr, which enables the neutrino-mass measurements in the SAP configuration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
NASA NTRS1973

Texturing in metals as a result of sliding

Sliding friction experiments were conducted with copper, nickel, iron, and cobalt sliding on themselves in air and argon. The resulting wear surfaces were examined with X-ray analysis to determine if surface texturing had occurred as a result of sliding. Results of the investigation indicate that, for the face-centered-cubic metals copper and nickel, a (111) texture develops with the (111) planes tilted 10 deg in the direction of sliding. The body-centered-cubic metal iron exhibited a (110) texture with the (100) direction oriented in the direction of sliding. It also exhibited a 10 deg tilt in the direction of sliding. The environment influenced the results in that the degree of texture observed in argon was less than that seen in air for iron. No texturing was observed for the close-packed-hexagonal metal cobalt. Recrystallization was observed with copper as a result of sliding.

Wheeler, D. R.↗
NASA NTRS1972

Analysis and design of three dimensional supersonic nozzles. Volume 1: Nozzle-exhaust flow field analysis by a reference plane characteristics technique

A second order numerical method employing reference plane characteristics has been developed for the calculation of geometrically complex three dimensional nozzle-exhaust flow fields, heretofore uncalculable by existing methods. The nozzles may have irregular cross sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. The nozzles may have highly nonuniform entrance conditions, the medium considered being an equilibrium hydrogen-air mixture. The program calculates and carries along the underexpansion shock and contact as discrete discontinuity surfaces, for a nonuniform vehicle external flow.

Dash, S.↗