Search NASA⌕ Search

SEARCH · Search NASA

Results for “Brillouin”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electronic and magnetic structures that hinder the superconducting state in the collapsed phase of $\mathrm{SrCr_2}$ $\mathrm{As_2}$

The lack of superconductivity in the family of materials $AB_2$ $\mathrm{As_2}$ ($A$ = $\mathrm{Sr, Ba}$ and $B$ = $\mathrm{Cr, Mn}$) under doping or pressure is an intriguing issue, considering the resemblance of these materials to their cousin materials and superconductors $\mathrm{(Ba, Sr) Fe_2}$ $\mathrm{As_2}$. In this context, the suppression of magnetism together with the presence of electron and hole pockets in the Brillouin zone seem to be fundamental ingredients to explain superconductivity. In this paper, we report a tetragonal to collapsed-tetragonal phase transition in $\mathrm{SrCr_2}$ $\mathrm{As_2}$ under high pressure with the appearance of hole and electron pockets in the Brillouin zone. In this collapsed-tetragonal phase, a residual local magnetic moment in the Cr-ion site and a finite global magnetic energy are derived. Importantly, this scenario would suggest why superconductivity is prevented in this compound. Our observations were obtained from the analysis of synchrotron x-ray diffraction measurements of polycrystalline samples of $\mathrm{SrCr_2}$ $\mathrm{As_2}$, and using first-principles simulations, under pressures $P$ in the range $\mathrm{1.4}$ $\mathrm{GPa}$ < $P$ < $\mathrm{20 GPa}$ at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evidence of Dirac Quantum Spin Liquid in YbZn 2⁢ GaO 5

The emergence of a quantum spin liquid (QSL), a state of matter that can result when electron spins are highly correlated but do not become ordered, has been the subject of a considerable body of research in condensed matter physics [1,2]. Spin liquid states have been proposed as hosts for high-temperature superconductivity [3] and can host topological properties with potential applications in quantum information science [4]. The excitations of most quantum spin liquids are not conventional spin waves but rather quasiparticles known as spinons, whose existence is well established experimentally only in one-dimensional systems; the unambiguous experimental realization of QSL behavior in higher dimensions remains challenging. Here, in this work, we investigate the novel compound $YbZn_{2}GaO_{5}$, which hosts an ideal triangular lattice of effective spin-1/2 moments with no detectable inherent chemical disorder. Thermodynamic and inelastic neutron scattering measurements performed on high-quality single crystal samples of $YbZn_{2}GaO_{5}$ exclude the possibility of long-range magnetic ordering down to 0.06 K, demonstrate a quadratic power law for the specific heat and reveal a continuum of magnetic excitations in parts of the Brillouin zone. Both low-temperature thermodynamics and inelastic neutron scattering spectra suggest that $YbZn_{2}GaO_{5}$ is a U(1) Dirac QSL with spinon excitations concentrated at certain points in the Brillouin zone. We advanced these results by performing additional specific heat measurements under finite fields, further confirming the theoretical expectations for a Dirac QSL on the triangular lattice of $YbZn_{2}GaO_{5}$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electronic structure in a rare-earth based nodal-line semimetal candidate PrSbTe

Nodal-line semimetals feature topologically protected band crossings between the bulk valence and conduction bands that extend along a finite dimension in the form of a line or a loop. While ZrSiS and similar materials have attracted extensive research as hosts for the nodal-line semimetallic phase, an alternative avenue has emerged in the form of isostructural rare-earth (RE)-based RESbTe materials. Such systems possess intriguing potentialities for harboring elements of magnetic ordering and electronic correlations owing to the presence of 4f electrons intrinsic to the RE elements. In this study, we have carried out angle-resolved photoemission spectroscopy (ARPES) and thermodynamic measurements in conjunction with first-principles computations on PrSbTe to elucidate its electronic structure and topological characteristics. Magnetic and thermal characterizations indicate the presence of well-localized 4f states with the absence of any discernible phase transition down to 2 K. The ARPES results reveal the presence of gapless Dirac crossings that correspond to a nodal-line along the X-R direction in the three-dimensional Brillouin zone. Furthermore, Dirac crossing that makes up the nodal line, which forms a diamond-shaped nodal plane centered at the center of the Brillouin zone, is also identified within the experimental resolution. In conclusion, this study on the electronic structure of PrSbTe contributes to the understanding of the pivotal role played by spin-orbit coupling in the context of the RESbTe family of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ultra-Long Distance BOTDA Sensor System Employing Hybrid Amplification and Advanced Noise Reduction Techniques

Brillouin Optical Time Domain Analysis (BOTDA) sensor system play a pivotal role in distributed sensing, which enables precise measurements of strain and temperature across extensive fiber lengths. Nonetheless, challenges emerge as distances grow due to signal attenuation and noise interference resulting in measurement errors. This research offers a comprehensive strategy to extend the sensing range of BOTDA systems beyond 10’s of kilometers while maintaining high spatial resolutions. Such enhanced sensing is realized through the integration of distributed Raman amplification, inline amplification using erbium-doped fiber amplifiers (EDFA), and advanced noise reduction techniques. Leveraging inherent redundancy in measured data as a function of frequency and fiber distance, the non-local means (NLM) filter removes noise while preserving essential physical information. This approach proves particularly advantageous in BOTDA systems, where accurate measurement of Brillouin scattering signals is paramount for long-range sensing, while concurrently safeguarding high spatial resolutions. In summary, this research has shown a holistic exploration of extending BOTDA's distance sensing capabilities up to 150 km with spatial resolutions of 8 meters.

Bhatta, Hari↗
DOE OSTI2020

Noninteracting electrons in a prototypical one-dimensional sinusoidal potential

A prototypical model of a one-dimensional metallic monatomic solid containing noninteracting electrons is studied, where the argument of the cosine potential energy, periodic with the lattice, contains the first reciprocal lattice vector G 1 =2π/a, where a is the lattice constant. The time-independent Schrödinger equation can be written in reduced variables as a Mathieu equation for which numerically exact solutions for the band structure and wave functions are obtained. The band structure has band gaps that increase with increasing amplitude q of the cosine potential. In the extended-zone scheme, the energy gaps decrease with increasing index n of the Brillouin-zone boundary ka=nπ, where k is the crystal momentum of the electron. The wave functions of the band electron are derived for various combinations of k and q as complex combinations of the real Mathieu functions with even and odd parity, and the normalization factor is discussed. The wave functions at the bottoms and tops of the bands are found to be real or imaginary, respectively, corresponding to standing waves at these energies. Irrespective of the wave vector k within the first Brillouin zone, the electron probability density is found to be periodic with the lattice. The Fourier components of the wave functions are derived versus q, which reveal multiple reciprocal-lattice-vector components with variable amplitudes in the wave functions unless q = 0. The magnitudes of the Fourier components are found to decrease exponentially as a power of n for n~3 to 45 for ka = π/2 and q = 2, and a precise fit is obtained to the data. The probability densities and probability currents obtained from the wave functions are also discussed. The probability currents are found to be zero for crystal momenta at the tops and bottoms of the energy bands, because the wave functions for these crystal momenta are standing waves. To conlcude, the band structure is calculated from the central equation and compared to the numerically exact band structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
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

Automating Analysis of Neutron Scattering Time-of-Flight Single Crystal Phonon Data

This article introduces software called Phonon Explorer that implements a data mining workflow for large datasets of the neutron scattering function, S(Q, ω), measured on time-of-flight neutron spectrometers. This systematic approach takes advantage of all useful data contained in the dataset. It includes finding Brillouin zones where specific phonons have the highest scattering intensity, background subtraction, combining statistics in multiple Brillouin zones, and separating closely spaced phonon peaks. Using the software reduces the time needed to determine phonon dispersions, linewidths, and eigenvectors by more than an order of magnitude.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Quantum limit studies in single crystal and pyrolytic graphite

The occupied regions of momentum space in graphite are confined to elongated volumes near the hexagonal Brillouin zone edges. The elongation is along the hexagonal axis (c-axis), so that cross sectional areas perpendicular to that axis are very small. When magnetic fields are applied parallel to the c-axis, only the n = 0 Landau levels remain below the Fermi energy for fields above 7 tesla. The energy bands in graphite can be described by a set of seven parameters (including the Fermi energy as a parameter). In a quantizing magnetic field the Landau levels are described by a secular equation determined by diagonalizing the Hamiltonian. A major feature for the levels along the Brillouin zone edge is that the n = 0 Landau level for electrons is independent of magnetic field, while that for holes is field-dependent.

Woollam, J. A.↗
NASA NTRS1971

Electron dynamics and energy conversion in O-type linear-beam devices.

A general nonlinear interaction theory is used to investigate the effects of transverse fields (i.e., radial circuit fields and radial space-charge fields) in traveling-wave amplifiers for a variety of beam-focusing conditions. Magnetic focusing fields which are periodic or tapered (increased) with distance along the device are considered in addition to uniform magnetic fields. Results are presented for Brillouin flow and near-Brillouin flow, and the minimum magnetic field strength required to effectively constrain the electron beam is determined as a function of the operating parameters for the various focusing systems. Confined flow is also examined for the uniform-field case in order to have a basis of comparison from which the effects of radial motion of the beam electrons can be determined. The results indicate the importance of transverse effects and further yield information on the stability of strongly modulated cylindrical electron beams.-

Detweiler, H. K.↗
NASA NTRS1974

Quantum limit studies in single-crystal and pyrolytic graphite

The occupied regions of momentum space in graphite are confined to very elongated volumes near the hexagonal Brillouin zone edges. The elongation is along the hexagonal axis (c axis) so that cross-sectional areas perpendicular to that axis are very small. Hence, when magnetic fields are applied parallel to the c axis only the n = 0 Landau levels remain below the Fermi energy for fields above 7 T (Woollam, 1971). The energy bands in graphite can be described by a set of seven parameters (McClure, 1971) (including the Fermi energy as a parameter). In a quantizing magnetic field the Landau levels are described by a secular equation determined by diagonalizing the Hamiltonian (McClure, 1960). A major unusual feature for the levels along the Brillouin zone edge is that the n = 0 Landau level for electrons is independent of magnetic field, while that for holes is field dependent.

Woollam, J. A.↗
NASA NTRS1984

Recent progress in the measurement of temperature and salinity by optical scattering

Recent progress is described in the use of Brillouin and Raman scattering for the measurement of temperature and salinity in the ocean. The use of Brillouin scattering is described for the measurement of the sound velocity, and the use of Raman scattering is described for the independent measurement of the temperature and salinity. Coupling these techniques permits the assessment of both temperature and salinity. The experimental techniques are described together with the results of recent experiments and an assessment of the errors to be expected.

Collins, D. J.↗
NASA NTRS1987

Mid infrared LHS system packaging using flexible waveguides

As mid IR fiber optic systems are rapidly approaching a reality, so is the feasibility of fiber optic laser heterodyne systems. Laser heterodyne spectroscopy for high resolution monitoring of atmospheric gaseous pollutants is necessarily in the mid IR, the region in which the absorption signature of gaseous species is most prominent. It so happens that the lowest theoretical loss due to Rayleigh-Brillouin scattering also lies in the mid IR. Prospects of highly efficient laser heterodyne systems are thus very good. Such fibers are now beginning to be commercially available, and a test program is being conducted for such fibers with ambient temperature ranging from cryogenic to above room, and stringest mechanical flexibility requirements. Preliminary results are encouraging. A program is being started to explore the possibility of mid IR fiber optic device applications, by taking advantage of this phonon rich region. The potential long interaction length in fibers coupled with predicted extremely low losses point to stimulated Brillouin scattering based devices in the mW range. The generation of backscattered sBs at low laser powers is significant not only as an ultimate power limiting factor for laser transmission in fibers in the mid IR, but also the presence of frequency-shifted multiple order sBs Stokes and antiStokes lines will certainly have severe effect on the laser beats crucial in high resolution heterodyne spectroscopy.

Yu, Chung↗
NASA NTRS1998

Double-Edge Molecular Technique for Doppler Lidar Wind Measurement

The double-edge lidar technique for measuring the wind using molecular backscatter is described. Two high spectral resolution edge filters are located in the wings of the Rayleigh-Brillouin profile. This doubles the signal change per unit Doppler shift, the sensitivity, and gives nearly a factor of two improvement in measurement accuracy. The use of a crossover region is described where the sensitivity of a molecular and aerosol-based measurement are equal. This desensitizes the molecular measurement to the effects of aerosol scattering over a frequency range of +/- 100 m/s. We give methods for correcting for short-term frequency jitter and drift using a laser reference frequency measurement and methods for long-term frequency correction using a servo control system. The effects of Rayleigh-Brillouin scattering on the measurement are shown to be significant and are included in the analysis. Simulations for a conical scanning satellite-based lidar at 355 nm show an accuracy of 2-3 m/s for altitudes of 2 to 15 km for a 1 km vertical resolution, a satellite altitude of 400 km and a 200 km x 200 km spatial resolution. Results of ground based wind measurements are presented.

Flesia, Cristina↗
NASA NTRS2001

Frequency conversion system

A frequency conversion system comprises first and second gain sources providing first and second frequency radiation outputs where the second gain source receives as input the output of the first gain source and, further, the second gain source comprises a Raman or Brillouin gain fiber for wave shifting a portion of the radiation of the first frequency output into second frequency radiation output to provided a combined output of first and second frequencies. Powers are gain enhanced by the addition of a rare earth amplifier or oscillator, or a Raman/Brillouin amplifier or oscillator between the high power source and the NFM device. Further, polarization conversion using Raman or Brillouin wavelength shifting is provided to optimize frequency conversion efficiency in the NFM device.

Sanders, Steven↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photoacoustic 3-D imaging of polycrystalline microstructure improved with transverse acoustic waves

Non-invasive fast imaging of grain microstructure of polycrystalline ceria with sub-micrometric spatial resolution is performed via time-domain Brillouin scattering. The propagation of a nanoacoustic pulse is monitored down to 8 µm deep in a 30 × 30 µm2 area. Grains boundaries are reconstructed in three-dimensions via a twostep processing method, relying on the wavelet synchro-squeezed transform and the alpha shape algorithm. Imaging contrast is improved by taking advantage of stronger sensitivity to anisotropy of transverse acoustic waves, compared with longitudinal waves. Utilization of transverse waves in the image processing reveals additional boundaries that are not discerned using longitudinal waves, confirmed by an electron backscattering diffraction pattern. A buried inclined interface between differently oriented grains is identified by monitoring changes in amplitude (phase) of the portion of the signal associated with transverse (longitudinal) waves. Estimates of the inclination angle of this interface prove the sensitivity of our laser ultrasonic method to image inclined boundaries.

36 MATERIALS SCIENCE↗