Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures

We present molecular dynamics (MD) calculations of shear viscosity for asymmetric mixed plasma for thermodynamic conditions relevant to astrophysical and inertial confinement fusion plasmas. Specifically, we consider mixtures of deuterium and argon at temperatures of 100–500 eV and a number density of 10 25 ions/cc. The motion of 30 000–120 000 ions is simulated in which the ions interact via the Yukawa (screened Coulomb) potential. The electric field of the electrons is included in this effective interaction; the electrons are not simulated explicitly. Shear viscosity is calculated using the Green-Kubo approach with an integral of the shear stress autocorrelation function, a quantity calculated in the equilibrium MD simulations. We systematically study different mixtures through a series of simulations with increasing fraction of the minority high- Z element (Ar) in the D-Ar plasma mixture. In the more weakly coupled plasmas, at 500 eV and low Ar fractions, results from MD compare very well with Chapman-Enskog kinetic results. In the more strongly coupled plasmas, the kinetic theory does not agree well with the MD results. Here, we develop a simple model that interpolates between classical kinetic theories at weak coupling and the Murillo Yukawa viscosity model at higher coupling. Finally, this hybrid kinetics-MD viscosity model agrees well with the MD results over the conditions simulated, ranging from moderately weakly coupled to moderately strongly coupled asymmetric plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Combining phase-field crystal methods with a Cahn-Hilliard model for binary alloys

Diffusion-induced phase transitions typically change the lattice symmetry of the host material. In battery electrodes, for example, Li ions (diffusing species) are inserted between layers in a crystalline electrode material (host). This diffusion induces lattice distortions and defect formations in the electrode. The structural changes to the lattice symmetry affect the host material’s properties. We propose a 2D theoretical framework that couples a Cahn-Hilliard (CH) model, which describes the composition field of a diffusing species, with a phase-field crystal (PFC) model, which describes the host-material lattice symmetry. We couple the two continuum models via coordinate transformation coefficients. We introduce the transformation coefficients in the PFC method to describe affine lattice deformations. These transformation coefficients are modeled as functions of the composition field. Using this coupled approach, we explore the effects of coarse-grained lattice symmetry and distortions on a diffusion-induced phase transition process. In this paper, we demonstrate the working of the CH-PFC model through three representative examples: First, we describe base cases with hexagonal and square symmetries for two composition fields. Next, we illustrate how the CH-PFC method interpolates lattice symmetry across a diffuse phase boundary. Finally, we compute a Cahn-Hilliard type of diffusion and model the accompanying changes to lattice symmetry during a phase transition process.

25 ENERGY STORAGE↗
DOE OSTI2018

Compact interaction potential for van der Waals nanorods

Here, we studied the van der Waals interactions of two finite, solid, cylindrical rods at arbitrary angle and position with respect to each other. An analytic interpolative formula for the interaction potential energy is constructed based on various asymptotic cases. In conclusion, the potential can be readily used for numerical and analytic descriptions of multiwall carbon nanotubes, metallic nanorods, rod-shaped colloids, or any other similar objects with significant van der Waals interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Bicollinear antiferromagnetic order, monoclinic distortion, and reversed resistivity anisotropy in FeTe as a result of spin-lattice coupling

The bicollinear antiferromagnetic order experimentally observed in FeTe is shown to be stabilized by the coupling g ~ 12 between monoclinic lattice distortions and the spin-nematic order parameter with B 2g symmetry, within a three-orbital spin-fermion model studied with Monte Carlo techniques. A finite but small value of g ~ 12 is required, with a concomitant lattice distortion compatible with experiments, and a tetragonal-monoclinic transition strongly first order. Remarkably, the bicollinear state found here displays a planar resistivity with the reversed puzzling anisotropy discovered in transport experiments. Orthorhombic distortions are also incorporated, and phase diagrams interpolating between pnictides and chalcogenides are presented. Here, we conclude that the spin-lattice coupling we introduce is sufficient to explain the challenging properties of FeTe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Theory for Glassy Behavior of Supercooled Liquid Mixtures

Here, we introduce a model for glassy dynamics in supercooled liquid mixtures. Given the relaxation behavior of individual supercooled liquids, the model predicts the relaxation times of their mixtures as temperature is decreased. The model is based on dynamical facilitation theory for glassy dynamics, which provides a physical basis for relaxation and vitrification of a supercooled liquid. This is in contrast to empirical linear interpolations such as the Gordon-Taylor equation typically used to predict glass transition temperatures of liquid mixtures. To understand the behavior of supercooled liquid mixtures we consider a multicomponent variant of the kinetically constrained East model in which components have a different energy scale and can also diffuse when locally mobile regions, i.e., excitations, are present. Using a variational approach we confirm an effective single component model with a single effective energy scale that best approximates a mixture. When scaled by this single effective energy, we show that experimental relaxation times of many liquid mixtures all collapse onto the "parabolic law" predicted by dynamical facilitation theory. The model can be used to predict transport properties and glass transition temperatures of mixtures of glassy materials, with implications in atmospheric chemistry, biology, and pharmaceuticals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Piezoelectric Electron-Phonon Interaction from Ab Initio Dynamical Quadrupoles: Impact on Charge Transport in Wurtzite GaN

First-principles calculations of $\textit{e}$-ph interactions are becoming a pillar of electronic structure theory. However, the current approach is incomplete. The piezoelectric (PE) $\textit{e}$-ph interaction, a long-range scattering mechanism due to acoustic phonons in noncentrosymmetric polar materials, is not accurately described at present. Current calculations include short-range $\textit{e}$-ph interactions (obtained by interpolation) and the dipolelike Frölich long-range coupling in polar materials, but lack important quadrupole effects for acoustic modes and PE materials. Here we derive and compute the long-range $\textit{e}$-ph interaction due to dynamical quadrupoles, and apply this framework to investigate $\textit{e}$-ph interactions and the carrier mobility in the PE material wurtzite GaN. In this work, we show that the quadrupole contribution is essential to obtain accurate $\textit{e}$-ph matrix elements for acoustic modes and to compute PE scattering. Our work resolves the outstanding problem of correctly computing $\textit{e}$-ph interactions for acoustic modes from first principles, and enables studies of $\textit{e}$-ph coupling and charge transport in PE materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Size-Dependent Nucleation in Crystal Phase Transition from Machine Learning Metadynamics

In this Letter, we present a framework that combines machine learning potential (MLP) and metadynamics to investigate solid-solid phase transition. Here, based on the spectral descriptors and neural networks regression, we develop a scalable MLP model to warrant an accurate interpolation of the energy surface where two phases coexist. Applying it to the simulation of B4–B1 phase transition of GaN under 50 GPa with different model sizes, we observe sequential change of the phase transition mechanism from collective modes to nucleation and growths. When the size is at or below 128 000 atoms, the nucleation and growth appear to follow a preferred direction. At larger sizes, the nuclei occur at multiple sites simultaneously and grow to microstructures by passing the critical size. The observed change of the atomistic mechanism manifests the importance of statistical sampling with large system size in phase transition modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bootstrapping the String Kawai-Lewellen-Tye Kernel

We show that a generalized version of the 4-point string theory KLT double-copy map is the most general solution to the minimal-rank double-copy bootstrap in effective field theory. This follows from significant restrictions of the 4-point map resulting from the 6-point bootstrap analysis. The generalized 4-point double-copy map is defined by a function with only two parameters times a simple function that is symmetric in s,t,u. The two parameters can be interpreted as independent choices of α', one for each of the two sets of amplitudes double-copied with the map. Specifically, each of those two sets of amplitudes must obey either the string monodromy relations or the field theory KK & BCJ relations; there are no other options. We propose a closed form of the new double-copy map that interpolates between the original KLT string double-copy and the open & closed string period integrals. The construction clarifies the "single-valued projection" property of the Riemann zeta-function values for the 4-point string theory double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Monopole-Fermion Scattering and the Solution to the Semiton–Unitarity Puzzle

We study Polchinski’s “fermion-rotor system” as an accurate description of charged Weyl fermions scattering on a magnetic monopole core in the limit of zero gauge coupling. Traditionally it was thought such scattering could lead to fractional particle numbers (“semitons”). By direct calculations we show those semitonic processes are in fact free propagation, facilitated by composite fermion-rotor operators interpolating the “forbidden” states, effectively “recovering” both ingoing and outgoing states in every lowest partial wave. Nonsemitonic Callan-Rubakov processes are unchanged.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hybrid SO(10) Axion Model without Quality Problem

Invisible axion models that solve the strong C P problem via the Peccei-Quinn (PQ) mechanism typically have a quality problem that arises from quantum gravity effects which violate all global symmetries. These models therefore require extreme fine-tuning of parameters for consistency. We present a new solution to the quality problem in a unified SO ( 10 ) × U ( 1 ) a gauge model, where U ( 1 ) a is an anomaly free axial gauge symmetry. PQ symmetry emerges as an accidental symmetry in this setup, which admits a PQ breaking scale as large as 7 × 10 11 GeV , allowing for the axion to be the cosmological dark matter. We call this a hybrid axion model due to its unique feature that it interpolates between the popular KSVZ and DFSZ axion models. Its predictions for the experimentally measurable axion couplings to the nucleon and electron are distinct from those of the usual models, a feature that can be used to test it. Furthermore, the model has no domain wall problem and it provides a realistic and predictive framework for fermion masses and mixings. Published by the American Physical Society 2025

Babu, K. S. (ORCID:0000000161475155)↗
DOE OSTI2019

Mini volume collapse as evidence for a three-body magnetic polaron in S m 1 - x E u x S

Samarium sulfide (SmS) is a nonmagnetic narrow-gap (0.06 eV) semiconductor that undergoes a transition to a metallic intermediate valence state at 6.5 kbar. Europium sulfide (EuS) is a ferromagnetic semiconductor with a Curie temperature of 16 K and a gap of 1.6 eV. Here we present a study of the lattice constant, magnetic susceptibility, and resistivity of the substitution series Sm 1-x Eu x S for 0 < x < 1. We observe a smooth interpolation of magnetic and transport behavior across the series, consistent with a virtual crystal scenario and Vegard’s law. Surprisingly, however, the lattice constant deviates below Vegard’s law in a manner that suggests parametric control of the Sm-Sm distance by the Eu moment in the manner of a magnetic polaron.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimal function estimation with photonic quantum sensor networks

The problem of optimally measuring an analytic function of unknown local parameters each linearly coupled to a qubit sensor is well understood, with applications ranging from field interpolation to noise characterization. Here we resolve a number of open questions that arise when extending this framework to Mach-Zehnder interferometers and quadrature displacement sensing. In particular, we derive lower bounds on the achievable mean square error in estimating a linear function of either local phase shifts or quadrature displacements. In the case of local phase shifts, these results prove, and somewhat generalize, a conjecture by Proctor []. For quadrature displacements, we extend proofs of lower bounds to the case of arbitrary linear functions. We provide optimal protocols achieving these bounds up to small (multiplicative) constants and describe an algebraic approach to deriving new optimal protocols, possibly subject to additional constraints. Using this approach, we prove necessary conditions for the amount of entanglement needed for any optimal protocol for both local phase and displacement sensing. Published by the American Physical Society 2024

Bringewatt, Jacob↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Data-Driven Compression of Electron-Phonon Interactions

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ( e − ph ) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e − ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e − ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e − ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e − ph interactions. Furthermore, they accelerate state-of-the-art first-principles e − ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e − ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Learning Conservation Laws in Unknown Quantum Dynamics

We present a learning algorithm for discovering conservation laws given as sums of geometrically local observables in quantum dynamics. This includes conserved quantities that arise from local and global symmetries in closed and open quantum many-body systems. The algorithm combines the classical shadow formalism for estimating expectation values of observable and data analysis techniques based on singular value decompositions and robust polynomial interpolation to discover all such conservation laws in unknown quantum dynamics with rigorous performance guarantees. Our method can be directly realized in quantum experiments, which we illustrate with numerical simulations, using closed and open quantum system dynamics in a ℤ 2 gauge theory and in many-body localized spin chains.

Machine learning↗
DOE OSTI2025

Window Observables for Benchmarking Parton Distribution Functions

Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel “window observables” that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken-𝑥, and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different window observables that can be defined within a region of 𝑥 where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.

lattice QCD↗