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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

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

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗
DOE OSTI2022

Higher-group symmetries and weak gravity conjecture mixing

In four-dimensional axion electrodynamics, a Chern-Simons coupling of the form $θF$ Λ $F$ leads to a higher-group global symmetry between background gauge fields. At the same time, such a Chern-Simons coupling leads to a mixing between the Weak Gravity Conjectures for the axion and the gauge field, so that the charged excitations of a Weak Gravity Conjecture-satisfying axion string will also satisfy the Weak Gravity Conjecture for the gauge field. In this paper, we argue that these higher-group symmetries and this phenomenon of Weak Gravity Conjecture mixing are related to one another. We show that this relationship extends to supergravities in 5, 6, 7, 8, 9, and 10 dimensions, so higher-dimensional supergravity is endowed with precisely the structure needed to ensure consistency with emergent higher-group symmetries and with the Weak Gravity Conjecture. We further argue that a similar mixing of Weak Gravity Conjectures can occur in two-term Chern-Simons theories or in theories with kinetic mixing, though the connection with higher-group symmetries here is more tenuous, and accordingly the constraints on effective field theory are not as sharp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The Classical Equations of Motion of Quantized Gauge Theories, Part I: General Relativity

In this and a companion paper, we show that quantum field theories with gauge symmetries permit a broader class of classical dynamics than typically assumed. In this article, we show that the dynamics extracted from the path integral or Hamiltonian formulation of general relativity allows for classical states that do not satisfy the full set of Einstein's equations. This amounts to loosening the Hamiltonian and momentum constraints that are imposed on the initial state. Nevertheless, the quantum theory permits gauge invariant time evolution of these states. The time evolution of these states is such that at the classical level the full set of Einstein's equations would appear to hold, with the physical effects of these states being attributable to an auxiliary, covariantly conserved energy-momentum tensor with no internal degrees of freedom. We derive the generalized Einstein equations for these states and show that a homogeneous and isotropic initial background state contributes to expansion identical to cold dark matter. The inhomogeneous components of this state could source curvature perturbations that grow linearly at linear order. This auxiliary contribution to Einstein's equations could have either sign and thus provide a trivial way to violate the null energy condition, enabling novel gravitational dynamics such as cosmic bounces and wormholes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Stress Tensor flows, birefringence in non-linear electrodynamics and supersymmetry

We identify the unique stress tensor deformation which preserves zero-birefringence conditions in non-linear electrodynamics, which is a 4d 4 d version of the T\overline{T} T T ¯ operator. We study the flows driven by this operator in the three Lagrangian theories without birefringence - Born-Infeld, Plebanski, and reverse Born-Infeld - all of which admit ModMax-like generalizations using a root- T\overline{T} T T ¯ -like flow that we analyse in our paper. We demonstrate one way of making this root- T\overline{T} T T ¯ -like flow manifestly supersymmetric by writing the deforming operator in \mathcal{N} = 1 𝒩 = 1 superspace and exhibit two examples of superspace flows. We present scalar analogues in d = 2 d = 2 with similar properties as these theories of electrodynamics in d = 4 d = 4 . Surprisingly, the Plebanski-type theories are fixed points of the classical T\overline{T} T T ¯ -like flows, while the Born-Infeld-type examples satisfy new flow equations driven by relevant operators constructed from the stress tensor. Finally, we prove that any theory obtained from a classical stress-tensor-squared deformation of a conformal field theory gives rise to a related “subtracted” theory for which the stress-tensor-squared operator is a constant.

Physics↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

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

Geometric origin of the energy-momentum tensor improvement terms

In a flat background, the canonical energy momentum tensor of Lorentz and conformally invariant matter field theories can be improved to a symmetric and traceless tensor that gives the same conserved charges. We argue that the geometric origin of this improvement process is unveiled when the matter theory is coupled to metric-affine gravity. In particular, we show that the Belinfante-Rosenfeld improvement terms correspond to the matter theory’s hypermomentum. The improvement terms in conformally invariant matter theories are also related to the hypermomentum; however, a general proof would require an extended investigation. We demonstrate our results through various examples, such as the free massless scalar, the Maxwell field, Abelian p-forms, the Dirac field, and a nonunitary massless scalar field. Possible applications of our method for theories that break Lorentz or special conformal invariance are briefly discussed.

classical solutions in field theory↗
DOE OSTI2025

Scalable quantum simulations of scattering in scalar field theory on 120 qubits

Simulations of collisions of fundamental particles on a quantum computer are expected to have an exponential advantage over classical methods and promise to enhance searches for new physics. Furthermore, scattering in scalar field theory has been shown to be bounded-error quantum polynomial time (BQP) complete, making it a representative problem for which quantum computation is efficient. As a step toward large-scale quantum simulations of collision processes, scattering of wave packets in one-dimensional scalar field theory is simulated using 120 qubits of IBM’s Heron superconducting quantum computer ibm_fez. Variational circuits compressing vacuum preparation, wave packet initialization, and time evolution are determined using classical resources. By leveraging physical properties of states in the theory, such as symmetries and locality, the variational quantum algorithm constructs scalable circuits that can be used to simulate arbitrarily large system sizes. A new strategy is introduced to mitigate errors in quantum simulations, which enables the extraction of meaningful results from circuits with up to 4924 two-qubit gates and two-qubit gate depths of 103. The effect of interactions is clearly seen, and is found to be in agreement with classical matrix product state simulations. Finally, the developments that will be necessary to simulate high-energy inelastic collisions on a quantum computer are discussed.

quantum circuits↗
DOE OSTI2021

Strongly coupled $\mathrm{QFT}$ dynamics via $\mathrm{TQFT}$ coupling

We consider a class of quantum field theories and quantum mechanics, which we couple to $\mathbb{Z}$ N topological QFTs, in order to classify non-perturbative effects in the original theory. The $\mathbb{Z}$ N TQFT structure arises naturally from turning on a classical background field for a $\mathbb{Z}$ N 0- or 1-form global symmetry. In SU(N) Yang-Mills theory coupled to $\mathbb{Z}$ N TQFT, the non-perturbative expansion parameter is exp[-S I /N] = exp[-8π 2 /g 2 N] both in the semi-classical weak coupling domain and strong coupling domain, corresponding to a fractional topological charge configurations. To classify the non-perturbative effects in original SU(N) theory, we must use PSU(N) bundle and lift configurations (critical points at infinity) for which there is no obstruction back to SU(N). These provide a refinement of instanton sums: integer topological charge, but crucially fractional action configurations contribute, providing a TQFT protected generalization of resurgent semi-classical expansion to strong coupling. Monopole-instantons (or fractional instantons) on T 3 x $S^1_L$ can be interpreted as tunneling events in the ’t Hooft flux background in the PSU(N) bundle. The construction provides a new perspective to the strong coupling regime of QFTs and resolves a number of old standing issues, especially, fixes the conflicts between the large-N and instanton analysis. We derive the mass gap at θ = 0 and gaplessness at θ = π in $\mathbb{CP}$ 1 model, and mass gap for arbitrary θ in $\mathbb{CP}$ N-1 , N ≥ 3 on $\mathbb{R}$ 2 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The constructive method for massive particles in QED

The constructive method of determining amplitudes from on-shell pole structure has been shown to be promising for calculating amplitudes in a more efficient way. However, challenges have been encountered when a massless internal photon is involved in the gluing of three-point amplitudes with massive external particles. In this paper, we use the original on-shell method, old-fashioned perturbation theory, to shed light on the constructive method, and show that one can derive the Feynman amplitude by correctly identifying the residue even when there is an internal photon involved.

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

Quantum Simulation of Field Theories Without State Preparation

We propose an algorithm for computing real-time observables using a quantum processor while avoiding the need to prepare the full quantum state. This reduction in quantum resources is achieved by classically sampling configurations in imaginary-time using standard lattice field theory. These configurations can then be passed to quantum processor for time-evolution. This method encounters a signal-to-noise problem which we characterize, and we demonstrate the application of standard lattice QCD methods to mitigate it.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

6D large charge and 2D Virasoro blocks

We compute observables in the interacting rank-one 6D 𝒩 =(2,0) superconformal field theory (SCFT) at large 𝑅-charge. We focus on correlators involving Φ 𝑛 , namely symmetric products of the bottom component of the supermultiplet containing the stress tensor. By using the moduli space effective action and methods from the large-charge expansion, we compute the operator product expansion coefficients ⟨Φ 𝑛 ⁢Φ 𝑚 ⁢Φ 𝑛+𝑚 ⟩ in an expansion in 1/𝑛. The coefficients of the expansion are only partially determined from the 6D perspective, but we manage to fix them order-by-order in 1/𝑛 numerically by utilizing the 6⁢D/2⁢D correspondence. This is made possible by the fact that this 6D observable can be extracted in 2D from a specific double-scaling limit of the vacuum Virasoro block, which can be efficiently computed numerically. We also extend the computation to higher-rank SCFTs, and discuss various applications of our results to 6D as well as 2D.

classical solutions in field theory↗
DOE OSTI2024

Effective Field Theory for Extreme Mass Ratio Binaries

We derive an effective field theory describing a pair of gravitationally interacting point particles in an expansion in their mass ratio, also known as the self-force (SF) expansion. The 0SF dynamics are trivially obtained to all orders in Newton’s constant by the geodesic motion of the light body in a Schwarzschild background encoding the gravitational field of the heavy body. The corrections at 1SF and higher are generated by perturbations about this configuration—that is, the geodesic deviation of the light body and the fluctuation graviton—but crucially supplemented by an operator describing the recoil of the heavy body as it interacts with the smaller companion. Using this formalism we compute new results at third post-Minkowskian order for the conservative dynamics of a system of gravitationally interacting massive particles coupled to a set of additional scalar and vector fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗