Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical mean field theory”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Observational Tests of Recent MHD Turbulence Perspectives

This grant seeks to analyze the Heliospheric Missions data to test current theories on the angular dependence (with respect to mean magnetic field direction) of magnetohydrodynamic (MHD) turbulence in the solar wind. Solar wind turbulence may be composed of two or more dynamically independent components. Such components include magnetic pressure-balanced structures, velocity shears, quasi-2D turbulence, and slab (Alfven) waves. We use a method, developed during the first two years of this grant, for extracting the individual reduced spectra of up to three separate turbulence components from a single spacecraft time series. The method has been used on ISEE-3 data, Pioneer Venus Orbiter, Ulysses, and Voyager data samples. The correlation of fluctuations as a function of angle between flow direction and magnetic-field direction is the focus of study during the third year.

Ghosh, Sanjoy↗
NASA NTRS1984

Planetary Mach cones - Theory and observation

The asymptotic behavior of planetary bow shocks and the ability of gas dynamic theory to describe it are studied. Spacecraft observations at Venus, earth, and Mars are used to model the shapes and positions of their distant bow waves. The measured planetary Mach cone angles are compared with the mean sonic and MHD fast wave Mach numbers at 0.7, 1.0, and 1.5 AU to assess the downstream influence of the interplanetary magnetic field. Numerical gas dynamic flow solutions extending to 50 obstacle radii behind each planet are obtained and tested against the observed location of the downstream bow shock.

Slavin, J. A.↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Crossflow Instability on a Wedge-Cone at Mach 3.5

As a follow-on activity to the HyBoLT flight experiment, a six degree half angle wedge-cone model at zero angle of attack has been employed to experimentally and computationally study the boundary layer crossflow instability at Mach 3.5 under low disturbance freestream conditions. Computed meanflow and linear stability analysis results are presented along with corresponding experimental Pitot probe data. Using a model-mounted probe survey apparatus, data acquired to date show a well defined stationary crossflow vortex pattern on the flat wedge surface. This effort paves the way for additional detailed, calibrated flow field measurements of the crossflow instability, both stationary and traveling modes, and transition-to-turbulence under quiet flow conditions as a means of validating existing stability theory and providing a foundation for dynamic flight instrumentation development.

Beeler, George B.↗
DOE OSTI2021

CrI 3 revisited with a many-body ab initio theoretical approach

CrI 3 has recently been shown to exhibit low-dimensional, long-range magnetic ordering from few layers to single layers of CrI 3 . The properties of CrI 3 bulk and few layered systems are uniquely defined by a combination of short-range intralayer and long-range interlayer interactions; including strong correlations, exchange and spin orbit coupling. Unfortunately, both the long-range van der Waals interactions, which are driven by dynamic, many-body electronic correlations and the competing strong intralayer correlations, present a formidable challenging for the local or semi local mean-field approximations employed in workhorse electronic structure approaches like density functional theory. In this work, we employ a sophisticated many-body approach that can simultaneously describe long- and short-range correlations. Herein, we establish that the fixed-node diffusion Monte Carlo (FNDMC) method reproduces the experimental interlayer separation distance of bulk CrI 3 for the high-temperature, monoclinic phase, with a reliable prediction of the interlayer binding energy. We subsequently employed, the FNDMC results to benchmark the accuracy of several density functional theory exchange-correlation approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissipative Superradiant Spin Amplifier for Enhanced Quantum Sensing

Quantum metrology protocols exploiting ensembles of N two-level systems and Ramsey-style measurements are ubiquitous. However, in many cases excess readout noise severely degrades the measurement sensitivity; in particular in sensors based on ensembles of solid-state defect spins. We present a dissipative “spin-amplification” protocol that allows us to dramatically improve the sensitivity of such schemes, even in the presence of realistic intrinsic dissipation and noise. Our method is based on exploiting collective (i.e., superradiant) spin decay, an effect that is usually seen as a nuisance because it limits spin-squeezing protocols. We show that our approach can allow a system with a highly imperfect spin readout to approach standard-quantum-limit-like scaling in N within a factor of 2, without needing to change the actual readout mechanism. Our ideas are compatible with several state-of-the-art experimental platforms where an ensemble of solid-state spins (NV centers, SiV centers) is coupled to a common microwave or mechanical mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Principles of two-dimensional terahertz spectroscopy of collective excitations: The case of Josephson plasmons in layered superconductors

Two-dimensional terahertz spectroscopy (2DTS), a terahertz analog of nuclear magnetic resonance, is a new technique poised to address many open questions in complex condensed matter systems. The conventional theoretical framework used ubiquitously for interpreting multidimensional spectra of discrete quantum level systems is, however, insufficient for the continua of collective excitations in strongly correlated materials. Here, in this work, we develop a theory for 2DTS of a model collective excitation, the Josephson plasma resonance in layered superconductors. Starting from a mean-field approach at temperatures well below the superconducting phase transition, we obtain expressions for the multidimensional nonlinear responses that are amenable to intuition derived from the conventional single-mode scenario. We then consider temperatures near the superconducting critical temperature T c , where dynamics beyond mean-field become important and conventional intuition fails. As fluctuations proliferate near T c , the dominant contribution to nonlinear response comes from an optical parametric drive of counterpropagating Josephson plasmons, which gives rise to 2D spectra that are qualitatively different from the mean-field predictions. As such, and in contrast to one-dimensional spectroscopy techniques, such as third harmonic generation, 2DTS can be used to directly probe thermally excited finite-momentum plasmons and their interactions. Our theory can readily be tested in cuprates, and we discuss implications beyond the present context of Josephson plasmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Generative deep-learning reveals collective variables of Fermionic systems

Complex processes of fermionic systems ranging from protein folding to nuclear fission often follow a low-dimensional reaction path parametrized in terms of a few collective variables. In nuclear theory, variables related to the shape of the nuclear density in a mean-field picture are key to describing the large amplitude collective motion of the neutrons and protons. Exploring the adiabatic energy landscape spanned by these degrees of freedom reveals the possible reaction channels while simulating the dynamics in this reduced space yields their respective probabilities. Unfortunately, this theoretical framework breaks down whenever the systems encounters a quantum phase transition with respect to the collective variables. Here, in this study, we introduce a novel generative deep-learning algorithm designed to build reaction paths that ensure that the many-fermion wave function stays differentiable with respect to the collective variables. This approach is applicable to any fermionic system described by a coherent state. We use the case of potential energy curves in the 16 O nucleus within the Hartree-Fock theory to illustrate its main features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Search for in-band transitions in the candidate superdeformed band in Si 28

Background: Superdeformed (SD) bands are suggested by theory around Ca 40 and in lighter alpha-conjugate nuclei such as Mg 24 , Si 28 , and S 32 . Such predictions originate from a number of theoretical models including mean-field models and antisymmetrized molecular dynamics (AMD) calculations. While SD bands have been identified in Ca 40 and its near neighbors, evidence of their existence in the lighter, midshell nuclei is circumstantial at best. Additionally, the key evidence of superdeformation would be the observation of transitions with high B ( E 2 ) transition strengths connecting states in a rotational sequence. This is challenging information to obtain since the bands lie at a high excitation energy and competition from out-of-band decay is dominant. Purpose: The purpose of the present study is to establish a new methodology to circumvent the difficulties in identifying and quantifying in-band transitions through directly populating candidate states in the SD band in Si 28 through inelastic alpha scattering, selecting such states with a spectrometer, and measuring their gamma-ray decay with a large array of high-purity germanium detectors, allowing direct access to electromagnetic transition strengths. Methods: Excited states in Si 28 were populated in the Si 28 ( α , α ' ) reaction using a 130-MeV He 4 beam from the K140 AVF cyclotron at the Research Center for Nuclear Physics. Outgoing alpha particles were analyzed using the Grand Raiden spectrometer positioned at an angle of 9 . 1 ° to favor the population of states with J ≈ 4 . Coincident gamma rays were detected with the CAGRA array of 12 HPGe clover detectors augmented by a set of four large LaBr 3 detectors. Results: Data analysis showed that it was possible to identify additional low-energy transitions in competition with high-energy decays from excited states in Si 28 in the vicinity of 10 MeV. However, while the candidate 4 + SD state at 10.944 MeV was populated, a 1148-keV transition to the candidate 2 + SD state at 9.796 MeV was not observed, and only an upper limit for its transition strength of B ( E 2 ) < 43 W.u. could be established. This contradicts AMD predictions of ≈ 200 W.u. for such a transition. Conclusion: The present study strongly rejects the hypothesis that the candidate set of states identified in Si 28 represents an SD band, which demonstrates the potential of the methodology devised here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hybrid algorithm for the time-dependent Hartree–Fock method using the Yang–Baxter equation on quantum computers *

Abstract The time-dependent Hartree–Fock (TDHF) method is an approach to simulate the mean field dynamics of electrons within the assumption that the electrons move independently in their self-consistent average field and within the space of single Slater determinants. One of the major advantages of performing time dynamics within Hartree–Fock theory is the free fermionic nature of the problem, which makes TDHF classically simulatable in polynomial time. Here, we present a hybrid TDHF implementation for quantum computers. This quantum circuit grows with time; but with our recent work on circuit compression via the Yang–Baxter equation (YBE), the resulting circuit is constant depth. This study provides a new way to simulate TDHF with the aid of a quantum device as well as provides a new direction for the application of YBE symmetry in quantum chemistry simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Solution Structure and Scaling Laws of Cylindrical and Tapered Bottlebrush Polymers

Bottlebrush polymers are a unique class of macromolecular architectures with a plethora of potential industrial and pharmaceutical applications that critically depend on the bottlebrush shape and dimensions. Here, a systematic series of 12 cylindrical and 12 cone-shaped (tapered) bottlebrush polymers with poly(tert-butyl acrylate) (PtBA) or polystyrene (PS) side chains were synthesized using the sequential addition of macromonomers ring-opening metathesis polymerization (SAM-ROMP) grafting-through method. Small-angle neutron scattering (SANS) studies on dilute solutions of the two types of bottlebrush polymers provided noninvasive characterization of their structural dimensions and chain conformations. Simulated SANS traces, generated using coarse-grained molecular dynamics simulations, reproduced the distinctive features observed in the experimental SANS signals and provided necessary validation for data modeling. The combined analysis of experimental and simulated SANS signals yielded key structural and conformational parameters, including the bottlebrush radius, length, and Kuhn length as well as the excluded volume parameter and the correlation length of the polymer side chains. Importantly, the obtained structural parameters followed well-defined scaling laws as a function of the backbone and side chain degrees of polymerization, as predicted by mean field theories. In conclusion, these findings provide clear experimental and computational evidence of the interdependence of structural and conformational properties in an important class of polymer architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗