Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “body wave”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Unexpected Rise in Nuclear Collectivity from Short-Range Physics

We discover a surprising relation between the collective motion of nucleons within atomic nuclei, traditionally understood to be driven by long-range correlations, and short-range nucleon-nucleon interactions. Specifically, we find that quadrupole collectivity in low-lying states of 6 Li and 12 C, calculated with state-of-the-art ab initio techniques, is significantly influenced by two opposing 𝑆-wave contact couplings that subtly alter the surface oscillations of one largely deformed nuclear shape, without changing that shape’s overall contribution within the nucleus. The results offer new insights into the nature of emergent nuclear collectivity and its link to the underlying nucleon-nucleon interaction at short distances.

Ab initio calculations↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Massive twistor worldline in electromagnetic fields

We study the (ambi-)twistor model for spinning particles interacting via electromagnetic field, as a toy model for studying classical dynamics of gravitating bodies including effects of both spins to all orders. We compute the momentum kick and spin kick up to one-loop order and show precisely how they are encoded in the classical eikonal. The all-orders-in-spin effects are encoded as a dynamical implementation of the Newman-Janis shift, and we find that the expansion in both spins can be resummed to simple expressions in special kinematic configurations, at least up to one-loop order. We confirm that the classical eikonal can be understood as the generator of canonical transformations that map the in-states of a scattering process to the out-states. We also remark that cut contributions for converting worldline propagators from time-symmetric to retarded amount to the iterated action of the leading eikonal at one-loop order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Covariant spectator theory of n p scattering: Deuteron form factors

The deuteron form factors are calculated using two model wave functions obtained from the 2007 CST high precision fits to $np$ scattering data. Included in the calculation are a new class of isoscalar $np$ interaction currents which are automatically generated by the nuclear force model used in these fits. If the nuclear model WJC2 is used, a precision fit ($\chi^2$/datum $\eqsim1$) to the Sick Global Analysis (GA) of all $ed$ elastic scattering data can be obtained by adjusting the unknown off-shell nucleon form factors $F_3(Q^2)$ (discussed before) and $F_4(Q^2)$ (introduced in this paper), and predicting the high $Q^2$ behavior of the neutron charge form factor $G_{En}(Q^2)$ well beyond the region where it has been measured directly. Relativistic corrections, isoscalar interaction currents, and off-shell effects are defined, discussed, and their size displayed. A rationale for extending $ed$ elastic scattering measurements to higher $Q^2$ is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

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

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Filtering states with total spin on a quantum computer

Starting from a general wave-function described on a set of spins or qubits, we propose several quantum algorithms to extract the components of this state on eigenstates of total spin S 2 and its z component S z . The method plays the role of total spin projection and gives access to the amplitudes of the initial state on a total spin basis. Here, the different algorithms have various degrees of sophistication depending on the requested tasks. They can either solely project onto the subspace with good total spin or completely uplift the degeneracy in this subspace. In the former case, when the projection on the total spin j is made, we show that the number of operations for the projection can be reduced from a quadratic to a linear dependence in j. After each measurement, the state collapses to one of the spin eigenstates that could be used for postprocessing. For this reason, we call the method total quantum spin filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Kolmogorov-Arnold wavefunctions

Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.

1-dimensional systems↗
DOE OSTI2023

Pair Density Wave Order from Electron Repulsion

A pair density wave (PDW) is a superconductor whose order parameter is a periodic function of space, without an accompanying spatially uniform component. Since PDWs are not the outcome of a weak-coupling instability of a Fermi liquid, a generic pairing mechanism for PDW order has remained elusive. Here, we describe and solve models having robust PDW phases. To access the intermediate coupling limit, we invoke large-$N$ limits of Fermi liquids with repulsive BCS interactions that admit saddle point solutions. We show that the requirements for long-range PDW order are that the repulsive BCS couplings must be nonmonotonic in space and that their strength must exceed a threshold value. We obtain a phase diagram with both finite temperature transitions to PDW order and a $T = 0$ quantum critical point, where non-Fermi liquid behavior occurs.

2-dimensional systems↗
DOE OSTI2021

Temperature and Density on the Forsterite Liquid‐Vapor Phase Boundary

Abstract The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid‐vapor phase boundary. The entire shock‐and‐release thermodynamic path must be calculated accurately in order to predict the post‐impact structures of planetary bodies. Forsterite (Mg 2 SiO 4 ) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid‐vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under‐predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.

58 GEOSCIENCES↗
DOE OSTI2020

Planck 2018 results XI. Polarized dust foregrounds

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization in the quest for the curl-like B-mode polarization from primordial gravitational waves and the low-multipole E-mode polarization associated with the reionization of the Universe. We used the new Planck PR3 maps to characterize Galactic dust emission at high latitudes as a foreground to the CMB polarization and use end-to-end simulations to compute uncertainties and assess the statistical significance of our measurements. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for a set of six nested high-Galactic-latitude sky regions covering from 24 to 71% of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra, which for the largest sky region are α EE = -2.42 ± 0.02 and α BB = -2.54 ± 0.02, respectively. The spectra show that the TE correlation and E/B power asymmetry discovered by Planck extend to low multipoles that were not included in earlier Planck polarization papers due to residual data systematics. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we have determined the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component-separation procedure that is required for measuring the low-ℓ reionization CMB E-mode signal and detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified black-body emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean dust spectral index for dust polarization is β d P = 1.53±0.02. The difference between indices for polarization and total intensity is β d P -β d I = 0.05±0.03. By fitting multi-frequency cross-spectra between Planck data at 100, 143, 217, and 353 GHz, we examine the correlation of the dust polarization maps across frequency. We find no evidence for a loss of correlation and provide lower limits to the correlation ratio that are tighter than values we derive from the correlation of the 217- and 353 GHz maps alone. If the Planck limit on decorrelation for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then frequency decorrelation of dust polarization might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio r ≃ 0.01 at the recombination peak. However, the Planck sensitivity precludes identifying how difficult the component-separation problem will be for more ambitious experiments targeting lower limits on r.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Using spin alignment of inelastically excited nuclei in fast beams to assign spins: The spectroscopy of 13 O as a test case

Excited states in O 13 were investigated using inelastic scattering of an E / A = 69.5 MeV O 13 beam off of a Be 9 target. The excited states were identified in the invariant-mass spectra of the decay products. Both single-proton and sequential two-proton decays of the excited states were examined. Furthermore, for a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave Born-approximation calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced O 13 level scheme is compared to ab initio no-core shell model with continuum predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1 s 1 / 2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Sound speed measurements in silver shock compressed to 300 GPa: Solid-state transition, melting, and liquid-state response

To gain insight into the thermodynamic response of shock compressed Ag states corresponding to the face-centered-cubic (fcc) to body-centered-cubic (bcc) transformation and melting observed using in-situ x-ray diffraction (XRD) measurements, we present longitudinal sound speed results and their analysis in Ag at peak stresses ranging from 60 to 300 GPa. The measured sound speeds increased linearly with density compression to 171 GPa, showing that the sound speeds (and longitudinal moduli) in the fcc and bcc phases are very similar. Between 171 and 218 GPa, the sound speed dropped significantly, consistent with the melting reported using XRD measurements in shock compressed Ag. From 218 to 300 GPa, the increasing sound speeds and Hugoniot states provide a determination of the liquid phase Ag response. In particular, determination of the Grüneisen parameter (Γ) showed that the density – Grüneisen parameter product (ρΓ) for liquid Ag is constant, but differs significantly from that for solid Ag at ambient conditions. Furthermore, the Mie- Grüneisen equation of state can be used to describe the Hugoniot and off-Hugoniot response of liquid Ag.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural study of hcp and liquid iron under shock compression up to 275 GPa

Here we combine nanosecond laser shock compression with in situ picosecond x-ray diffraction to provide structural data on iron up to 275 GPa. We constrain the extent of hcp-liquid coexistence, the onset of total melt, and the structure within the liquid phase. Our results indicate that iron, under shock compression, melts completely by 258(8) GPa. A coordination number analysis indicates that iron is a simple liquid at these pressure-temperature conditions. We also perform texture analysis between the ambient body-centered-cubic (bcc) α, and the hexagonal-closed-packed (hcp) high-pressure ε-phase. We rule out the Rong-Dunlop orientation relationship (OR) between the α and ε - phase s. However, we cannot distinguish between three other closely related ORs: Burger's, Mao-Bassett-Takahashi, and Potter's OR. The solid-liquid coexistence region is constrained from a melt onset pressure of 225(3) GPa from previously published sound speed measurements and full melt [246.5(1.8)–258(8) GPa] from x-ray diffraction measurements, with an associated maximum latent heat of melting of 623 J/g. This value is lower than recently reported theoretical estimates and suggests that the contribution to the earth's geodynamo energy budget from heat release due to freezing of the inner core is smaller than previously thought. Melt pressures for these nanosecond shock experiments are consistent with gas gun shock experiments that last for microseconds, indicating that the melt transition occurs rapidly.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Is the Matrix Completion of Reduced Density Matrices Unique?

Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.

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

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗