Search NASA⌕ Search

SEARCH · Search NASA

Results for “BRILLOUIN ZONES”

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

Interorbital nematicity and the origin of a single electron Fermi pocket in FeSe

The electronic structure of the enigmatic iron-based superconductor FeSe has puzzled researchers since spectroscopic probes failed to observe the expected electron pocket at the Y point in the 1-Fe Brillouin zone. It has been speculated that this pocket, essential for an understanding of the superconducting state, is either absent or incoherent. In this study, we perform a theoretical study of the preferred nematic order originating from nearest-neighbor Coulomb interactions in an electronic model relevant for FeSe. We find that at low temperatures the dominating nematic components are of interorbital d x z – d x y and d y z – d x y character, with spontaneously broken amplitudes for these two components. This interorbital nematic order naturally leads to distinct hybridization gaps at the X and Y points of the 1-Fe Brillouin zone, and may thereby produce highly anisotropic Fermi surfaces with only a single electron pocket at one of these momentum-space locations. The associated superconducting gap structure obtained with the generated low-energy electronic band structure from spin-fluctuation mediated pairing agrees well with that measured experimentally. Finally, from a comparison of the computed spin susceptibility to available neutron scattering data, we discuss the necessity of additional self-energy effects, and explore the role of orbital-dependent quasiparticle weights as a minimal means to include them.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Topology invisible to eigenvalues in obstructed atomic insulators

Here we consider the extent to which symmetry eigenvalues reveal the topological character of bands. Specifically, we compare distinct atomic limit phases (band representations) that share the same irreducible representations (irreps) at all points in the Brillouin zone and, therefore, appear equivalent in a classification based on eigenvalues. We derive examples where such “irrep-equivalent” phases can be distinguished by a quantized Berry phase or generalization thereof. These examples constitute a generalization of the Su-Schrieffer-Heeger chain: neither phase is topological, in the sense that localized Wannier functions exist, yet there is a topological obstruction between them. We refer to two phases as “Berry obstructed atomic limits” if they have the same irreps, but differ by Berry phases. This is a distinct notion from eigenvalue obstructed atomic limits, which differ in their symmetry irreps at some point in the Brillouin zone. We compute exhaustive lists of elementary band representations that are irrep equivalent, in all space groups, with and without time-reversal symmetry and spin-orbit coupling, and use group theory to derive a set of necessary conditions for irrep equivalence. Finally, we conjecture, and in some cases prove, that irrep-equivalent elementary band representations that are not equivalent can be distinguished by a topological invariant.

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

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2019

Effective dynamics of two-dimensional Bloch electrons in external fields derived from symmetry

We develop a comprehensive theory for the effective dynamics of Bloch electrons based on symmetry. We begin with a scheme to systematically derive the irreducible representations (IRs) characterizing the Bloch eigenstates in a crystal. Starting from a tight-binding (TB) approach, we decompose the TB basis functions into localized symmetry-adapted atomic orbitals and crystal-periodic symmetry-adapted plane waves. Each of these two subproblems is independent of the details of a particular crystal structure and it is largely independent of the relevant aspects of the other subproblem, hence permitting for each subproblem an independent universal solution. Taking monolayer MoS 2 and few-layer graphene as examples, we tabulate the symmetrized p and d orbitals as well as the symmetrized plane-wave spinors relevant for these crystal structures. The symmetry-adapted basis functions block-diagonalize the TB Hamiltonian such that each block yields eigenstates transforming according to one of the IRs of the group of the wave vector G k . For many crystal structures, it is possible to define multiple distinct coordinate systems such that for wave vectors k at the border of the Brillouin zone the IRs characterizing the Bloch states depend on the coordinate system, i.e., these IRs of G k are not uniquely determined by the symmetry of a crystal structure. The different coordinate systems are related by a coordinate shift that results in a rearrangement of the IRs of G k characterizing the Bloch states. We illustrate this rearrangement with three coordinate systems for MoS 2 and trilayer graphene. The freedom to choose different distinct coordinate systems can simplify the symmetry analysis of the Bloch states. Given the IRs of the Bloch states in one coordinate system, a rearrangement lemma yields immediately the IRs of the Bloch states in the other coordinate systems. The rearrangement of the IRs in different coordinate systems does not affect observable physics such as selection rules or the effective Hamiltonians for the dynamics of Bloch states in external fields. Using monolayer MoS 2 as an example, we combine the symmetry analysis of its bulk Bloch states with the theory of invariants to construct a generic multiband Hamiltonian for electrons near the K point of the Brillouin zone. Here, the Hamiltonian includes the effect of spin-orbit coupling, strain, and external electric and magnetic fields. Invariance of the Hamiltonian under time reversal yields additional constraints for the allowed terms in the Hamiltonian and it determines the phase (real or imaginary) of the prefactors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Polar catastrophe and the structure of KT a 1 - x N b x O 3 surfaces: Results from elastic and inelastic helium atom scattering

In this paper, the structure and dynamics of cleaved (001) surfaces of potassium tantalates doped with niobium, KTa 1-x Nb x O 3 (KTN), with x ranging from 0% to 30%, were measured by helium atom scattering (HAS). Through HAS time-of-flight (TOF) experiments, a dispersionless branch (Einstein phonon branch) with energy of 13-14meV was observed across the surface Brillouin zone in all samples. When this observation is combined with the results from earlier experimental and theoretical studies on these materials, a consistent picture of the stable surface structure emerges: After cleaving the single-crystal sample, the surface should be composed of equal areas of KO and TaO 2 /NbO 2 terraces. The data, however, suggest that K + and O 2- ions migrate from the bulk to the surface, forming a charged KO lattice that is neutralized primarily by additional K + ions bridging pairs of surface oxygens. This structural and dynamic modification at the (001) surface of KTN appears due to its formally charged KO(-1) and TaO 2 /NbO 2 (+1) layers and avoids a “polar catastrophe.” Finally, this behavior is contrasted with the (001) surface behavior of the fluoride perovskite KMnF 3 with its electrically neutral KF and MnF 2 layers.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Slater Insulator in Iridate Perovskites with Strong Spin-Orbit Coupling

The perovskite (Pv) SrIrO 3 is an exotic narrow-band metal owing to a confluence of the strengths of the spin-orbit coupling (SOC) and the electron-electron correlations. It has been proposed that topological and magnetic insulating phases can be achieved by tuning the SOC, Hubbard interactions, and/or lattice symmetry. Here, we report that the substitution of nonmagnetic, isovalent Sn 4+ for Ir 4+ in the SrIr 1-x Sn x O 3 perovskites synthesized under high pressure leads to a metal-insulator transition to an antiferromagnetic phase at TN ≥ 225 K. The continuous change of the cell volume as detected by X-ray diffraction and the l-shape transition of the specific heat on cooling through TN demonstrate that the metal-insulator transition is of second-order. Neutron powder diffraction results indicate that the Sn substitution enlarges an octahedral-site distortion that reduces the SOC relative to the spin-spin exchange interaction and results in the type-G AF spin ordering below TN. Measurement of high-temperature magnetic susceptibility shows the evolution of magnetic coupling in the paramagnetic phase typical of weak itinerant-electron magnetism in the Sn-substituted samples. A reduced structural symmetry in the magnetically ordered phase leads to an electron gap opening at the Brillouin zone boundary below TN in the same way as proposed by Slater.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2014

Elastic waves at periodically-structured surfaces and interfaces of solids

This paper presents a simple treatment of elastic wave scattering at periodically structured surfaces and interfaces of solids, and the existence and nature of surface acoustic waves (SAW) and interfacial (IW) waves at such structures. Our treatment is embodied in phenomenological models in which the periodicity resides in the boundary conditions. These yield zone folding and band gaps at the boundary of, and within the Brillouin zone. Above the transverse bulk wave threshold, there occur leaky or pseudo-SAW and pseudo-IW, which are attenuated via radiation into the bulk wave continuum. These have a pronounced effect on the transmission and reflection of bulk waves. We provide examples of pseudo-SAW and pseudo-IW for which the coupling to the bulk wave continuum vanishes at isloated points in the dispersion relation. These supersonic guided waves correspond to embedded discrete eigenvalues within a radiation continuum. We stress the generality of the phenomena that are exhibited at widely different scales of length and frequency, and their relevance to situations as diverse as the guiding of seismic waves in mine stopes, the metrology of periodic metal interconnect structures in the semiconductor industry, and elastic wave scattering by an array of coplanar cracks in a solid.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence of distinct nonpolar/polar ordering at long/short ranges in relaxor ferroelectrics

Distinct atomic nonpolar/polar ordering is determined at long, intermediate, and short ranges using synchrotron x-ray diffraction, Raman scattering, and pair distribution function data in K 0.5 Na 0.5 NbO 3 -doped BaTiO 3 -based relaxor ferroelectrics, viz., (1 - x) Ba 0.9 Sr 0.1 TiO 3 -x KNN 50. Two different polar phases with distinct symmetries are identified at short/intermediate ranges for x = 0.10, where a monoclinic phase at short ranges gradually transforms into a rhombohedral phase at intermediate ranges. Here, in contrast, a nonpolar phase with average cubic symmetry is found to be stable at long ranges for a wide temperature range (100 ≤ T ≤ 500 K). The polar behavior of short/intermediate-range ordering is quantified in terms of the amplitude of ferroelectric phonon mode $Γ^-_4$ (q = 0, 0, 0) (which corresponds to the zone center of the cubic Brillouin zone). The amplitude of the $Γ^-_4$ phonon mode increases with decreasing temperature, suggesting the enhancement of ferroelectric polarization at short/intermediate ranges on lowering temperatures. Enhanced ferroelectric polarization at low temperature is also visible and evident by ferroelectrostriction and is quantified by the increase in magnitude of spontaneous volume ferroelectrostriction. Consequently, Zero thermal expansion is observed at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Origin of superconductivity in layered centrosymmetric LaNiGa{sub 2}

We have examined the origin of superconductivity in layered centrosymmetric LaNiGa{sub 2} by employing a linear response approach based on the density functional perturbation theory. Our results indicate that this material is a conventional electron-phonon superconductor with intermediate level of coupling strength, with the electron-phonon coupling parameter of 0.70, and the superconducting critical temperature of 1.90 K in excellent accordance with experimental value of 1.97 K. The largest contribution to the electron-phonon coupling comes from the La d and Ga p electrons near the Fermi energy and the B{sub 3g} phonon branch resulting from vibrations of these atoms along the Γ-Z symmetry line in the Brillouin zone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Observation of Dirac-like semi-metallic phase in NdSb

The search of new topological phases of matter is one of the new directions in condensed matter physics. Recent experimental realizations of Dirac semimetal phases pave the way to look for other exotic phases of matter in real materials. In this paper, we present a systematic angle-resolved photoemission spectroscopy (ARPES) study of NdSb, a potential candidate for hosting a Dirac semi-metal phase. Our studies reveal two hole-like Fermi surface pockets present at the zone center ($\Gamma$ ) point as well as two elliptical electron-pockets present in the zone corner (X) point of the Brillouin zone (BZ). Interestingly, Dirac-like linearly dispersive states are observed about the zone corner (X) point in NdSb. Our first-principles calculations agree with the experimentally observed bands at the $\Gamma$ point. Moreover, the Dirac-like state observed in NdSb may be a novel correlated state, not yet predicted in calculations. Finally, our study opens a new direction to look for Dirac semi-metal states in other members of the rare earth monopnictide family.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Electronic and vibrational properties of lithium doped graphene

In the frame-work of density functional theory calculation, using planewave pseudopotentials within local density approximation, the electronic and vibrational properties of graphene supercell by adsorption of lithium at three different sites top, hollow and bridge, have been systematically investigated and analyzed. We found that the hollow site is the most favorable site having lowest energy and positive phonon frequency throughout Brillouin zone indicating dynamical stability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1966

Electronic spectrum of crystalline antimony.

Energy levels and interband oscillator strengths of antimony calculated in Brillouin zone by pseudopotential method, predicting polarization effects and spin-orbit splittings

BRILLOUIN ZONE↗
DOE OSTI2015

A framework for solving atomistic phonon-structure scattering problems in the frequency domain using perfectly matched layer boundaries

We present a numerical approach to the solution of elastic phonon-interface and phonon-nanostructure scattering problems based on a frequency-domain decomposition of the atomistic equations of motion and the use of perfectly matched layer (PML) boundaries. Unlike molecular dynamic wavepacket analysis, the current approach provides the ability to simulate scattering from individual phonon modes, including wavevectors in highly dispersive regimes. Like the atomistic Green's function method, the technique reduces scattering problems to a system of linear algebraic equations via a sparse, tightly banded matrix regardless of dimensionality. However, the use of PML boundaries enables rapid absorption of scattered wave energies at the boundaries and provides a simple and inexpensive interpretation of the scattered phonon energy flux calculated from the energy dissipation rate in the PML. The accuracy of the method is demonstrated on connected monoatomic chains, for which an analytic solution is known. The parameters defining the PML are found to affect the performance and guidelines for selecting optimal parameters are given. The method is used to study the energy transmission coefficient for connected diatomic chains over all available wavevectors for both optical and longitudinal phonons; it is found that when there is discontinuity between sublattices, even connected chains of equivalent acoustic impedance have near-zero transmission coefficient for short wavelengths. The phonon scattering cross section of an embedded nanocylinder is calculated in 2D for a wide range of frequencies to demonstrate the extension of the method to high dimensions. The calculations match continuum theory for long-wavelength phonons and large cylinder radii, but otherwise show complex physics associated with discreteness of the lattice. Examples include Mie oscillations which terminate when incident phonon frequencies exceed the maximum available frequency in the embedded nanocylinder, and scattering efficiencies larger than two near the Brillouin zone edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Stabilizing the Metallic State in a Quasi-Two-Dimensional (ET){sub 8}Hg{sub 4}Br{sub 12}(C{sub 6}H{sub 5}Cl){sub 2} Organic Metal under a Pressure of 6 kbar

A quasi-two-dimensional (ET){sub 8}Hg{sub 4}Br{sub 12}(C{sub 6}H{sub 5}Cl){sub 2} organic metal undergoes the metal–dielectric transition at T ~ 90 K upon cooling. An external pressure of P ~ 6 kbar completely recovers the metallic state and makes it possible to observe the Shubnikov–de Hass oscillations with a frequency of F ≈ 215 T and a cyclotron mass of m* ≈ 0.95m{sub 0}. The behavior of the oscillations is consistent with the Fermi surface representation in the form of electron and hole orbits covering the same areas of the first Brillouin zone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

On the influence of tetrahedral covalent-hybridization on electronic band structure of topological insulators from first principles

Based on first-principles calculations, we investigate the influence of tetrahedral covalent-hybridization between main-group and transition-metal atoms on the topological band structures of binary HgTe and ternary half-Heusler compounds, respectively. Results show that, for the binary HgTe, when its zinc-blend structure is artificially changed to rock-salt one, the tetrahedral covalent-hybridization will be removed and correspondingly the topologically insulating band character lost. While for the ternary half-Heusler system, the strength of covalent-hybridization can be tuned by varying both chemical compositions and atomic arrangements, and the competition between tetrahedral and octahedral covalent-hybridization has been discussed in details. As a result, we found that a proper strength of tetrahedral covalent-hybridization is probably in favor to realizing the topologically insulating state with band inversion occurring at the Γ point of the Brillouin zone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗