Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “derived lattice”

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

Structural Modifications of Single-Crystal Aragonite CaCO 3 Beginning at ~15 GPa: In Situ Vibrational Spectroscopy and X-Ray Diffraction Evidence

The structural chemistry of carbonates under mantle conditions facilitates our understanding of carbon recycling pathways in the earth’s interior. It also has impacts on the dynamics of mantle–slab interactions. Aragonite is a common calcium carbonate mineral in pelagic marine sediments. The structural chemistry of single-crystal aragonite during successive compression and the behavior of a structural H+ have been investigated by micro-vibrational spectroscopy and synchrotron X-ray diffraction techniques in diamond anvil cells. We describe a reduction of the b-axial compressibility beginning at ~15 GPa, and the related discontinuities in the first-order derivatives of the vibrational modes. The structural modifications of aragonite are manifested by mutations occurring in the pressure relations of the wavenumbers of the O-C-O bending modes, and of the bandwidth and band intensities of the measured internal and external modes. These anomalies are indicative of changes occurring in the force constant of the C-O bonds, and possibly a second-order phase transition. Besides, the [CaO 9 ] polyhedra begin to deform, possibly with some Ca-O bonds becoming elongated and the others shortening. An increase in the co-ordination number for the Ca 2+ sites could be expected under higher pressures. Additionally, the weakening of the OH modes may imply H + -loss from the aragonite lattice above 11.5 GPa.

58 GEOSCIENCES↗
DOE OSTI2023

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation

Understanding the microstructure–property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure–property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure–permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. Here high-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

58 GEOSCIENCES↗
DOE OSTI2025

Bottomonium suppression in pNRQCD and open quantum system approach

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Anomalies of non-invertible symmetries in (3+1)d

Anomalies of global symmetries are important tools for understanding the dynamics of quantum systems. We investigate anomalies of non-invertible symmetries in 3+1d using 4+1d bulk topological quantum field theories given by Abelian two-form gauge theories, with a 0-form permutation symmetry. Gauging the 0-form symmetry gives the 4+1d “inflow” symmetry topological field theory for the non-invertible symmetry. We find a two levels of anomalies: (1) the bulk may fail to have an appropriate set of loop excitations which can condense to trivialize the boundary dynamics, and (2) the “Frobenius-Schur indicator” of the non-invertible symmetry (generalizing the Frobenius-Schur indicator of 1+1d fusion categories) may be incompatible with trivial boundary dynamics. As a consequence we derive conditions for non-invertible symmetries in 3+1d to be compatible with symmetric gapped phases, and invertible gapped phases. Along the way, we see that the defects characterizing \mathbb{Z}_{4} ℤ 4 ordinary symmetry host worldvolume theories with time-reversal symmetry \mathsf{T} 𝖳 obeying the algebra \mathsf{T}^{2}=C 𝖳 2 = C or \mathsf{T}^{2}=(-1)^{F}C, 𝖳 2 = ( − 1 ) F C , with C C a unitary charge conjugation symmetry. We classify the anomalies of this symmetry algebra in 2+1d and further use these ideas to construct 2+1d topological orders with non-invertible time-reversal symmetry that permutes anyons. As a concrete realization of our general discussion, we construct new lattice Hamiltonian models in 3+1d with non-invertible symmetry, and constrain their dynamics.

Córdova, Clay↗
DOE OSTI2022

QCD analysis of Δ S = 0 hadronic parity violation

We present a QCD analysis of the effective weak Hamiltonian at hadronic energy scales for strangeness-nonchanging (ΔS = 0) hadronic processes. Performing a leading-order renormalization group analysis in QCD from the W to the O(2 GeV) energy scale, we derive the pertinent effective Hamiltonian for hadronic parity violation, including the effects of both neutral and charged weak currents. We compute the complete renormalization group evolution of all isosectors and the evolution through heavy-flavor thresholds for the first time. We show that the additional four-quark operators that enter below the W mass scale from QCD operator mixing effects form a closed set, and they result in a 12 x 12 anomalous dimension matrix. Computing the resulting effective Hamiltonian and comparing to earlier results, we affirm the importance of operator mixing effects and find, as an example, that the parity-violating pion-nucleon coupling constant, using the factorization Ansatz and an assessment of the pertinent quark charge of the nucleon in lattice QCD at the 2 GeV scale, is in better agreement with recent experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2022

Multiphysics FISPACT-II and TENDL-2019 simulation: neutron-induced damage metrics

Nuclear interactions can be the source of atomic displacement, irradiation-induced defects and transmutation in structural materials. Such quantities are derived from, or can be correlated to, nuclear kinematic simulations of primary atomic energy distributions spectra and the quantification of the numbers of secondary defects produced per primary as a function of the available recoils, residual and emitted, energies. Recoil kinematics of neutral, residual, gas, proton, alpha particle emissions are now more rigorously treated based on recent, complete and enhanced nuclear data parsed in state-of-the-art processing tools. Defect production metrics are the starting point in the complex problem of correlating and simulating the behaviour of materials under irradiation, as experimental information is rare or scattered. Detailed, segregated primary knock-on-atom metrics are now becoming available as the starting point of further simulation processes of isolated and clustered defects in material lattices. This allows more materials, lattices, neutron incident energy ranges, and irradiation conditions to be explored with sufficient data to adequately cover both standard and novel applications and materials: the broader reactor applications landscape. The damage metrics of of materials are systematically explored under typical but different reactor's type environment. The inventory code FISPACT-II combined with the enhanced nuclear data forms of the TENDL-2019 libraries allow one to not only calculate dpa from mostly scattering events but to also properly predict gas production, nuclear heating and transmutation under the same conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Master integrals for two-loop QCD corrections to quark quasi PDFs

We compute the master integrals for the two-loop QCD corrections to quark quasi parton distribution functions (PDFs) in the large momentum effective theory (LaMET). With a proper canonical basis, we derive the analytical results for the three families of master integrals using the method of differential equations. The final expressions for the master integrals are given in terms of Goncharov polylogarithms. These results allow us to extract the two-loop short-distant matching coefficients between quark quasi and lightcone PDFs in LaMET, and are valuable to improve the determination of the nucleon PDFs from first principles in future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Topological flat bands in frustrated kagome lattice CoSn

Electronic flat bands in momentum space, arising from strong localization of electrons in real space, are an ideal stage to realize strongly-correlated phenomena. Theoretically, the flat bands can naturally arise in certain geometrically frustrated lattices, often with nontrivial topology if combined with spin-orbit coupling. Here, we report the observation of topological flat bands in frustrated kagome metal CoSn, using angle-resolved photoemission spectroscopy and band structure calculations. Throughout the entire Brillouin zone, the bandwidth of the flat band is suppressed by an order of magnitude compared to the Dirac bands originating from the same orbitals. The frustration-driven nature of the flat band is directly confirmed by the chiral d-orbital texture of the corresponding real-space Wannier functions. Spin-orbit coupling opens a large gap of 80 meV at the quadratic touching point between the Dirac and flat bands, endowing a nonzero Z2 invariant to the flat band. These findings demonstrate that kagome-derived flat bands are a promising platform for novel emergent phases of matter at the confluence of strong correlation and topology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Light–ferroelectric interaction in two-dimensional lead iodide perovskites

The unique physical properties of two-dimensional (2D) metal halide perovskites (MHPs) such as nonlinear optics, anisotropic charge transport, and ferroelectricity have made these materials promising candidates for multifunctional applications. Recently, fluorine derivatives such as 4,4-difluoropiperidinium lead iodide perovskite or (4,4-DFPD, C 5 H 10 F 2 N) 2 PbI 4 have shown strong ferroelectricity as compared to other 2D MHPs. Although it was previously addressed that the ferroelectricity in MHPs can be affected by illumination, the underlying physical mechanisms of light–ferroelectricity interaction in 2D MHPs are still lacking. Here, we explore the electromechanical responses in 4,4-(DFPD) 2 PbI 4 thin films using advanced scanning probe microscopy techniques revealing ferroelectric domain structures. In this work, hysteretic ferroelectric loops measured by contact-Kelvin probe force microscopy are dependent on domain structures under dark conditions, while ferroelectricity weakens under illumination. The X-ray diffraction patterns exhibit significant changes in preferential orientation of individual lattice planes under illumination. Particularly, the reduced intensity of the (1 1 1) lattice plane under illumination leads to transitioning from a ferroelectric to a paraelectric phase. The instability of positive ions, especially molecular organic cations, is observed under illumination by time-of-flight secondary ion mass spectrometry. The combination of crystallographic orientation and chemical changes under illumination clearly contributes to the origin of light–ferroelectricity interaction in 2D (4,4-DFPD, C 5 H 10 F 2 N) 2 PbI 4 .

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

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Interpolation of compound semiconductor alloy parameters from those of their constituents

Several methods have been proposed for interpolation of the value of physical parameters of quaternary alloys from those of their constituent ternary and binary sub-alloys. These expressions agree when non-linear bowing terms are not required; they differ in how the bowing terms of the bounding ternaries should be utilized. Common interpolation expressions for quaternaries can be generalized into two groups: (1) those that use a linear interpolation of the nearest ternary parameter values and (2) those that interpolate over binary values with a bowing term derived from the bounding ternaries. The second group of methods is equivalent to a polynomial expansion over the alloy’s interpolation space. For compound semiconductor alloys, the geometry of the composition space is the direct sum of the group-III and group-V mixture sub-spaces. The mixture sub-spaces are best described using barycentric coordinates on a regular simplex. A general polynomial expansion of the value of an alloy parameter using barycentric coordinates for the group-III and group-V simplex spaces is described along with an algorithm to generate interpolation expressions for alloys with arbitrary numbers of elements, including quinary and senary alloys. It is shown that a polynomial expansion produces values in closer agreement with the direct gap of quaternaries lattice-matched to common substrates than do approaches using an interpolation of the ternary values, despite a prominent recommendation to the contrary. Finally, a quaternary correction term is described that improves the predicted direct bandgap energies of GaInAsSb for compositions near those lattice matched to InP, InAs, and GaSb.

Olesberg, Jonathon T. [Sandia National Laboratorie↗
DOE OSTI2024

Geometric origin of the intrinsic transverse spin transport in a canted-antiferromagnet/heavy-metal heterostructure

We theoretically study the conditions under which an intrinsic spin Nernst effect–a transverse spin current induced by an applied temperature gradient–can occur in a canted-antiferromagnet insulator, such as LaFeO 3 and other materials of the same family. The spin Nernst effect may provide a microscopic mechanism for an experimentally observed anomalous thermovoltage in LaFeO 3 /Pt heterostructures, where spin is transferred across the insulator/metal interface when a temperature gradient is applied to LaFeO 3 parallel to the interface. We find that LaFeO 3 exhibits an intrinsic spin Nernst effect when inversion symmetry is broken on the axes parallel to both the applied temperature gradient and the direction of spin transport, which can result in a spin injection across the insulator/metal interface. Furthermore, our paper provides a general derivation of a symmetry-breaking-induced spin Nernst effect, which may open a path to engineering a finite spin Nernst effect in systems where it would otherwise not arise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

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

Numerical Evaluation of Effective Thermal Conductivity of PCM with Metal Foam Incorporating Buoyancy Effects for Thermal Energy Storage

The thermal energy storage (TES) system has the capability to efficiently preserve thermal energy directly derived from the energy source, minimizing any conversion losses. Especially latent heat storage offers distinct advantages, including a substantial increase in energy storage density and minimization of temperature fluctuations within the plants. However, the phase change material (PCM) employed in latent heat storage has low thermal conductivity. Consequently, various studies are being conducted to enhance heat transfer. One approach to enhance heat transfer involves utilizing metal foam to maximize the heat transfer area. However, modeling metal foam with its intricate structure is a challenging task in numerical analysis. For this reason, ongoing research focuses on simplifying the modeling of metal foam. Nevertheless, fully encompassing all the characteristics of actual metal foam proves to be a challenging task for the simplified analytical model. The objective of this paper is to interpret the simple lattice metal foam analysis model from the perspective of behavior induced by buoyancy. When comparing the analysis results of solid PCM and liquid PCM with the same thermal conductivity under changes in porosity and gravity direction, we conducted an analysis to discern the trends in effective thermal conductivity that are overestimated due to convection. In the analysis, a constant heat flux of 10 kW and a constant surface boundary condition of 350 K were applied, and a sensitivity study regarding the mesh was conducted. The results indicate that, from the perspective of gravity in the simple lattice model, the solid analysis yields an effective thermal conductivity 29-47% higher compared to the liquid analysis. Additionally, as porosity increases, there is an observed increase of 24-33% in effective thermal conductivity.

25 ENERGY STORAGE↗
DOE OSTI2022

Lattice dynamics of high-pressure hydrides studied by inelastic neutron scattering

Due to the small mass and anomalously large neutron scattering cross-section of proton (about 80 barns compared to a few barns for other nuclei), inelastic neutron scattering is considered as one of the most effective tools in studying optical vibrations of hydrogen atoms in metal hydrides. The current review is focused on the binary hydrides of 3d- and 4d-metals of groups VI–VIII, which were produced at high hydrogen pressures of several gigapascals in relatively large quantities of hundreds of mg, quenched to low temperature and studied by INS ex situ at ambient pressure with high statistical accuracy. One of the unusual effects revealed by INS is a strong increase in the strength of the metal-hydrogen interactions with decreasing atomic number of the d-metal accompanied by an increase in the Me-H distance. Based on the available experimental results, the spectra g(E) of the phonon density of states and temperature dependencies CV(T) of the heat capacity at constant volume at T up to 1000 K have been derived in this paper and presented both in the figures and in digital form. This provides the reference data for the theoretical investigations of the crystal structures and compositions of new practically important hydrides giving the opportunity to validate calculation methods by comparing the calculated g(E) and CV(T) with the accurate experimental dependencies for the binary hydrides. Recent INS studies showed [R.A. Klein et al., J. Alloy. Compd. 894 (2022) 162381] that the fingerprints of anomalously short H-H separations of 1.6 Å violating the “2 Å rule” can be easily and unambiguously identified in the complex INS spectra of quaternary hydrides (La,Ce)NiInH 1+x . This makes neutron spectroscopy an attractive means for obtaining valuable data in the search for novel hydrides with a record high hydrogen capacity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗