Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classical 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Systematic input scheme for many-boson Hamiltonians with applications to the two-dimensional 𝜙 4 theory

We develop a novel, systematic input scheme for many-boson Hamiltonians in order to solve field theory problems within the light-front Hamiltonian formalism via quantum computing. We present our discussion of this input scheme based on the light-front Hamiltonian of the two-dimensional ϕ 4 theory. In our input scheme, we employ a set of quantum registers, where each register encodes the occupation of a distinct boson mode as binaries. We squeeze the boson operators of each mode and present the Hamiltonian in terms of unique combinations of the squeezed-boson operators. We design the circuit modules for these unique combinations. Based on these circuit modules, we block encode the many-boson Hamiltonian utilizing the idea of quantum walk. For demonstration purposes, we present the spectral calculations of the Hamiltonian utilizing the hybrid quantum-classical symmetry-adapted quantum Krylov subspace diagonalization algorithm based on our input scheme, where the quantum computations are performed with the IBM Qiskit quantum simulator. The results of the hybrid calculations agree with exact results. Here, we can incorporate the input scheme in this work with the input scheme for many-fermion Hamiltonians; they jointly offer new pathways to solving the structure and dynamics of more general field theory problems on future fault-tolerant quantum computers.

Ab initio calculations↗
DOE OSTI2022

Computational Modeling of Battery Materials

This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.

25 ENERGY STORAGE↗
DOE OSTI2022

Using Hyperoptimized Tensor Networks and First-Principles Electronic Structure to Simulate the Experimental Properties of the Giant {Mn 84 } Torus

The single-molecule magnet {Mn 84 } is a challenge to theory because of its high nuclearity. Here, we directly compute two experimentally accessible observables, the field-dependent magnetization up to 75 T and the temperature-dependent heat capacity, using parameter-free theory. In particular, we use first-principles calculations to derive short- and long-range exchange interactions and compute the exact partition function of the resulting classical Potts and Ising spin models for all 84 Mn S = 2 spins to obtain observables. The latter computation is made possible by using hyperoptimized tensor network contractions, a technique developed to simulate quantum supremacy circuits. We also synthesize the magnet and measure its heat capacity and magnetization, observing qualitative agreement between theory and experiment and identifying an unusual bump in the heat capacity and a plateau in the magnetization. Our work also identifies some limitations of current theoretical modeling in large magnets, such as sensitivity to small, long-range exchange couplings.

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

Functional equations and separation of variables for exact g -function

The g-function is a measure of degrees of freedom associated to a boundary of two-dimensional quantum field theories. In integrable theories, it can be computed exactly in a form of the Fredholm determinant, but it is often hard to evaluate numerically. In this paper, we derive functional equations — or equivalently integral equations of the thermodynamic Bethe ansatz (TBA) type — which directly compute the g-function in the simplest integrable theory; the sinh-Gordon theory at the self-dual point. The derivation is based on the classic result by Tracy and Widom on the relation between Fredholm determinants and TBA, which was used also in the context of topological string. We demonstrate the efficiency of our formulation through the numerical computation and compare the results in the UV limit with the Liouville CFT. As a side result, we present multiple integrals of Q-functions which we conjecture to describe a universal part of the g-function, and discuss its implication to integrable spin chains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Accretion of matter and spectra of binary X-ray sources in massive gravity

We study low-mass binary X-ray sources involving stellar mass black holes within massive gravity. Regarding the accretion disk, we adopt the standard model for an optically thick, cool, and geometrically thin disk by Shakura–Sunyaev. For the gravitational field generated by the black hole, we consider the analogue of the Schwarzschild–de Sitter space–time of Einstein’s theory in massive gravity, for which we found an additional term linear in the radial coordinate. Then, we compute the radial velocity, the energy density and the pressure as a function of the radial coordinate, and the X-ray emission’s soft spectral component expected from the disk. We also investigated in detail the impact of this new geometry. Our result indicates that by using observed spectra from confirmed X-ray binaries involving astrophysical black holes, we can put strong constraints on alternative theories of gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Theory of gradient drift instabilities in low-temperature, partially magnetised plasmas

A fluid dispersion theory in partially magnetised plasmas is analysed to examine the conditions under which large-wavelength modes develop in Penning-type configurations, that is, where an electric field is imposed perpendicular to a homogeneous magnetic field. The fluid dispersion relation assuming a slab geometry shows that two types of low-frequency, gradient drift instabilities occur in the direction of the E×B and diamagnetic drifts. One type of instability, observed when the equilibrium electric field and plasma density gradient are in the same direction, is similar to the classic modified Simon–Hoh instability. A second instability is found for conditions in which (i) the diamagnetic drift is in the direction opposite to the E×B drift and (ii) the magnitude of the diamagnetic drift is sufficiently larger than the electron thermal speed. The present fluid dispersion theory suggests that the rotating spokes driven by such fluid instabilities propagate in the same direction as the diamagnetic drift, which can be in the same direction as or opposite to the E×B drift, depending on the plasma conditions. This finding may account for the observation, in some plasma devices, of the rotation of large-scale structures in both the E×B and -E×B directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2022

Representations and strategies for transferable machine learning improve model performance in chemical discovery

Strategies for machine-learning (ML)-accelerated discovery that are general across material composition spaces are essential, but demonstrations of ML have been primarily limited to narrow composition variations. By addressing the scarcity of data in promising regions of chemical space for challenging targets such as open-shell transition-metal complexes, general representations and transferable ML models that leverage known relationships in existing data will accelerate discovery. Over a large set (~1000) of isovalent transition-metal complexes, we quantify evident relationships for different properties (i.e., spin-splitting and ligand dissociation) between rows of the Periodic Table (i.e., 3d/4d metals and 2p/3p ligands). We demonstrate an extension to the graph-based revised autocorrelation (RAC) representation (i.e., eRAC) that incorporates the group number alongside the nuclear charge heuristic that otherwise overestimates dissimilarity of isovalent complexes. To address the common challenge of discovery in a new space where data are limited, we introduce a transfer learning approach in which we seed models trained on a large amount of data from one row of the Periodic Table with a small number of data points from the additional row. We demonstrate the synergistic value of the eRACs alongside this transfer learning strategy to consistently improve model performance. Analysis of these models highlights how the approach succeeds by reordering the distances between complexes to be more consistent with the Periodic Table, a property we expect to be broadly useful for other material domains.

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

Path integral approach to quantum anomalies in interacting models

The prediction and subsequent discovery of topological semimetal phases of matter in solid state systems has instigated a surge of activity investigating the exotic properties of these unusual materials. Among these are transport signatures which can be attributed to the chiral anomaly; the breaking of classical chiral symmetry in a quantum theory. This remarkable quantum phenomenon, first discovered in the context of particle physics has now found new life in condensed-matter physics, connecting topological quantum matter and band theory with effective field theoretic models. In this paper we investigate the interplay between interactions and the chiral anomaly in field theories inspired by semimetals using Fujikawa's path integral method. Starting from models in one spatial dimension we discuss how the presence of interactions can affect the consequences of the chiral anomaly leading to renormalization of excitations and their transport properties. This is then generalized to the three-dimensional case where we show that the anomalous response of the system, namely, the chiral magnetic and quantum Hall effects, are modified by the presence of interactions. These properties are investigated further through the identification of anomalous modes which exist within interacting Weyl semimetals. These massive excitations are nonperturbative in nature and are a direct consequence of the chiral anomaly. The effects of interactions on mixed axial-gravitational anomalies are then investigated and the conditions required for interaction effects to be observed are discussed.

Anomalies↗
DOE OSTI2023

Revisiting Experimental Signatures of the Ponderomotive Force

The classical theory of single-electron dynamics in focused laser pulses is the foundation of both the relativistic ponderomotive force (RPF), which underlies models of laser-collective-plasma dynamics, and the discovery of novel strong-field radiation dynamics. Despite this bedrock importance, consensus eludes the community as to whether acceleration of single electrons in vacuum has been observed in experimental conditions. We analyze an early experiment on the RPF with respect to several features that were neglected in modeling and that can restore consistency between theory predictions and experimental data. The right or wrong pulse profile function, laser parameters, or initial electron distribution can each make or break the agreement between predictions and data. The laser phase at which the electron’s interaction with the pulse begins has a large effect, explaining why much larger energies are achieved by electrons liberated in the focal region by photoionization from high-Z atoms and by electrons ejected from a plasma mirror. Finally, we compute the difference in a typical electron spectrum arising from fluctuating focal spot size in state-of-the-art ultra-relativistic laser facilities. Our results emphasize the importance of thoroughly characterizing laser parameters in order to achieve quantitatively accurate predictions and the precision required for discovery science.

Hegelich, Bjorn Manuel↗
DOE OSTI2022

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Infrared spectroscopy of RG–Co + (H 2 O) complexes (RG = Ar, Ne, He): The role of rare gas “tag” atoms

RG n –Co + (H 2 O) cation complexes (RG = Ar, Ne, He) are generated in a supersonic expansion by pulsed laser vaporization. Here, complexes are mass-selected using a time-of-flight spectrometer and studied with infrared laser photodissociation spectroscopy, measuring the respective mass channels corresponding to the elimination of the rare gas “tag” atom. Spectral patterns and theory indicate that the structures of the ions with a single rare gas atom have this bound to the cobalt cation opposite the water moiety in a near-C 2v arrangement. The O–H stretch vibrations of the complex are shifted compared to those of water because of the metal cation charge-transfer interaction; these frequencies also vary systematically with the rare gas atom attached. The efficiencies of photodissociation also vary with the rare gas atoms because of their widely different binding energies to the cobalt cation. The spectrum of the argon complex could only be measured when at least three argon atoms were attached. In the case of the helium complex, the low binding energy allows the spectra to be measured for the low-frequency H–O–H scissors bending mode and for the O–D stretches of the deuterated analog. The partially resolved rotational structure for the antisymmetric O–H and O–D stretches reveals the temperature of these complexes (6 K) and establishes the electronic ground state. The helium complex has the same 3 B 1 ground state as the tag-free complex studied previously by Metz and co-workers (2013), but the A rotational constant is contaminated by vibrational averaging from the bending motion of the helium.

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

Boson Fermion Nucleus Explanation for Fission Fragment Anisotropy: Monograph #13

Since 1952, anisotropy has been observed in the paths of fission fragments following the assimilation of an incident neutron or photon in the nucleus. Fission fragment flux favors certain angles relative to the path of the incident particle, constituting anisotropy. The amount of anisotropy is dependent on the energy of the incident particle. Fission fragment anisotropy can be explained by the disruption of nuclear structure combined with the requirement of a specific rotational orientation in theta and phi of that structure relative to the incident neutron or photon flux. A specific and limited area in the nucleus can assimilate the incident particle and induce upset within the structure, followed by scission and fission of the nucleus. The orientation and nature of the structure upset dictates fission fragment trajectory and is responsible for fission fragment anisotropy. This concept can be explained and supported by the theory of boson fermion nucleus (BFN) structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗