Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice dynamics”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dipolar Bose-Hubbard model

We study a simple model of interacting bosons on a d-dimensional cubic lattice whose dynamics conserves both total boson number and total boson dipole moment. This model provides a simple framework in which several remarkable consequences of dipole conservation can be explored. As a function of chemical potential and hopping strength, the model can be tuned between gapped Mott insulating phases and various types of gapless condensates. The condensed phase realized at large hopping strengths, which we dub a Bose-Einstein insulator, is particularly interesting: despite having a Bose condensate, it is insulating, and despite being an insulator, it is compressible.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗
DOE OSTI2023

Precisely computing phonons via irreducible derivatives

Computing phonons from first principles is typically considered a solved problem, yet inadequacies in existing techniques continue to yield deficient results in systems with sensitive phonons. Here, in this study, we circumvent this issue using the lone irreducible derivative (LID) and bundled irreducible derivative (BID) approaches to computing phonons via finite displacements, where the former optimizes precision via energy derivatives and the latter provides the most efficient algorithm using force derivatives. A condition number optimized basis for BID is derived which guarantees the minimum amplification of error. Additionally, a hybrid LID-BID approach is formulated, in which select irreducible derivatives computed using LID replace BID results. We illustrate our approach on two prototypical systems with sensitive phonons: the shape memory alloy AuZn and metallic lithium. Comparing our resulting phonons in the aforementioned crystals to calculations in the literature reveals nontrivial inaccuracies. Our approaches can be fully automated, making them well suited for both niche systems of interest and high-throughput approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. As a result, it also paves the way for many-body computations of the phonon spectrum in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Occupational disorder as the origin of flattening of the acoustic phonon branches in the clathrate Ba 8 ⁢Ga 16 ⁢Ge 30

In the search for high-performance thermoelectrics, materials such as clathrates have drawn attention due to having both glasslike low phonon thermal conductivity and crystal-like high electrical conductivity. Ba 8 ⁢Ga 16 ⁢Ge 30 (BGG) has a loosely bound guest Ba atom trapped inside rigid Ga-Ge cage structures. Avoided crossings between acoustic phonons and the flat guest atom branches have been proposed to be the source of the low lattice thermal conductivity of BGG. Ga-Ge site disorder with Ga and Ge exchanging places in different unit cells has also been reported. We used time-of-flight neutron scattering to measure the complete phonon spectrum in a large single crystal of BGG and compared these results with predictions of density functional theory to elucidate the effect of the disorder on heat-carrying phonons. Experimental results agreed much better with the calculation assuming the disorder than with the calculation assuming the ordered configuration. Although the atomic masses of Ga and Ge are nearly identical, we found that disorder strongly reduces phonon group velocities, which significantly reduces thermal conductivity. Finally, our work points to a path towards optimizing thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural evolution of iodine on approach to the monatomic state

Here, we applied single-crystal x-ray diffraction and Raman spectroscopy in a diamond-anvil cell up to 36 GPa and first-principles theoretical calculations to study the molecular dissociation of solid iodine at high pressure. Unlike previously reported, we find that the familiar Cmce molecular phase transforms to a Cmc2 1 molecular structure at 16 GPa, and then to an incommensurate polymeric Fmmm⁡(00⁢$\gamma$)⁢s⁢00 structure, which can be viewed as a mixture of molecular I 2 and zigzag chains of three iodine atoms, at 20 GPa. The available data are consistent with metallization of iodine at the boundary between the commensurate Cmc⁢2 1 structure and the incommensurate Fmmm⁡(00⁢$\gamma$)s⁢00 structure at 20 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Rearrangement collision theory of phonon-driven exciton dissociation

Understanding the processes governing the dissociation of excitons to free charge carriers in semiconductors and insulators is of central importance for photovoltaic applications. Dyson's S-matrix formalism provides a framework for computing scattering rates between quasiparticle states derived from the same underlying Hamiltonian, often reducing to familiar Fermi's “golden rule” like expressions at first order. By presenting a rigorous formalism for multichannel scattering, we extend this approach to describe scattering between composite quasiparticles and, in particular, the process of exciton dissociation mediated by the electron-phonon interaction. Subsequently, we derive rigorous expressions for the exciton dissociation rate, a key quantity of interest in optoelectronic materials, which enforce correct energy conservation and may be readily used in ab initio calculations. We apply our formalism to a three-dimensional model system to compare temperature-dependent exciton rates obtained for different scattering channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Soft phonon and the central peak at the cubic-to-tetragonal phase transition in SrTiO 3

The continuous displacive phase transition in SrTiO 3 near 𝑇 𝑐 ≈105 K features a central elastic peak in neutron-scattering investigations at temperatures above 𝑇 𝑐 , i.e., before the corresponding soft phonon mode is overdamped upon cooling. The origin of this central peak is still not understood. Here, in this work, we report an inelastic x-ray scattering investigation of the cubic-to-tetragonal phase transition in SrTiO 3 . We compare quantitatively measurements of the soft phonon mode on two differently grown samples and discuss the findings regarding results from thermodynamic and transport probes such as specific heat and thermal conductivity. Furthermore, we use inelastic x-ray scattering to perform elastic scans with both high momentum and milli-electronvolt energy-resolution and, thus, are able to separate elastic intensities of the central peak from low-energy quasielastic phonon scattering. Our results indicate that the evolution of the soft mode is similar in both samples though the intensities of the central peak differ by a factor of four. Measurements revealing anisotropic correlation lengths on cooling towards 𝑇 𝑐 , indicate that local properties of the crystals to which collective lattice excitations are insensitive are likely at the origin of the central elastic line in SrTiO 3 .

anharmonic lattice dynamics↗
DOE OSTI2021

Role of Static Displacements in Stabilizing Body Centered Cubic High Entropy Alloys

The configurational entropy of high entropy alloys (HEAs) plays little role in the stabilization of one particular crystal structure over another. We show that disorder-induced atomic displacements help stabilize body centered cubic (bcc) structure HEAs with average valences <4.7. These disorder-induced atomic displacements mimic the temperature-induced vibrations that stabilize the bcc structure of group IV elemental metals at high temperatures. The static displacements are significantly larger than for face centered cubic HEAs, approaching values associated with the Lindemann criterion for melting. Chemical disorder in high entropy alloys have a previously unidentified, nonentropic energy contribution that stabilizes a particular crystalline ground state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Autoregressive Neural Network for Simulating Open Quantum Systems via a Probabilistic Formulation

The theory of open quantum systems lays the foundation for a substantial part of modern research in quantum science and engineering. Rooted in the dimensionality of their extended Hilbert spaces, the high computational complexity of simulating open quantum systems calls for the development of strategies to approximate their dynamics. In this Letter, we present an approach for tackling open quantum system dynamics. Using an exact probabilistic formulation of quantum physics based on positive operator-valued measure, we compactly represent quantum states with autoregressive neural networks; such networks bring significant algorithmic flexibility due to efficient exact sampling and tractable density. We further introduce the concept of string states to partially restore the symmetry of the autoregressive neural network and improve the description of local correlations. Efficient algorithms have been developed to simulate the dynamics of the Liouvillian superoperator using a forward-backward trapezoid method and find the steady state via a variational formulation. Our approach is benchmarked on prototypical one-dimensional and two-dimensional systems, finding results which closely track the exact solution and achieve higher accuracy than alternative approaches based on using Markov chain Monte Carlo method to sample restricted Boltzmann machines. Our Letter provides general methods for understanding quantum dynamics in various contexts, as well as techniques for solving high-dimensional probabilistic differential equations in classical setups.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Heteroanionic Ruddlesden-Popper ferroelectrics from anion order and octahedral tilts

We describe a strategy to design ferroelectric heteroanionic materials based on the coupling of anion order and octahedral tilts in n = 1 Ruddlesden-Popper structures, which allows noncentrosymmetric and polar compounds to arise from centrosymmetric anion-ordered structures. Here we investigate the relative phase stabilities of the polymorphs of Sr 2 ScO 3 F and Ca 2 ScO 3 F derived from this coupling using electronic structure calculations. We find that large degrees of octahedral tilting can stabilize different anion orders derived from assembly of [ScO 5 F] 8– octahedra. To further understand the link between octahedral tilting and anion order, we quantitatively separate the contributions of electrostatics and covalent interactions to the stability of tilts. We find that the tilts are driven primarily by covalent interactions and that local out-of-plane polar displacements, induced from the anion order, further stabilize the octahedral tilts through the pseudo-Jahn-Teller effect.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inelastic neutron spectra of polyacrylonitrile-based carbon fibers

We present the vibrational spectra of polyacrylonitrile-based carbon fibers collected using inelastic neutron scattering. We ascertain the behavior of a broad range of vibrational spectra that are optically silent, and demonstrate a direct connection between these modes and thermomechanical properties of the fibers. We show directionally dependent coupling of hydrogen in the carbon fiber matrix, which is directly linked to mechanical properties. Further, we show hydrogen preferentially couples to the midband of the vibrational spectrum, and that there are higher overall mode populations in the traditional Raman D–G intervalley region, suggesting involvement of these modes in tensile strength reduction and transport properties.

36 MATERIALS SCIENCE↗