Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Feynman”

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 55 records · Page 3

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence for a QCD accelerator in relativistic heavy-ion collisions

Here, we report measurements of forward jets produced in Cu + Au collisions at $\sqrt{s{NN}}$=200 GeV at the Relativistic Heavy Ion Collider. The jet-energy distributions extend to energies much larger than expected by Feynman scaling. This constitutes the first clear evidence for Feynman-scaling violations in heavy-ion collisions. Such high-energy particle production has been in models via QCD string interactions, but so far is untested by experiment. One such model calls this a hadronic accelerator. Studies with a particular heavy-ion event generator ( HIJING ) show that photons and mesons exhibit such very high-energy production in a heavy-ion collision, so a QCD accelerator appropriately captures the physics associated with such QCD string interactions. All models other than HIJING used for hadronic interactions in the study of extensive air showers from cosmic rays either do not include these QCD string interactions or have smaller effects from the QCD accelerator.

43 PARTICLE ACCELERATORS↗

Intrinsic charm and the 𝐷 + − 𝐷 − asymmetry produced in proton-proton collisions

We investigate the contribution of the charm-anticharm (𝑐⁢$\overline{𝑐}$) asymmetry of the proton eigenstate obtained from QCD lattice gauge to the asymmetry of 𝐷 + , 𝐷 − and 𝐷 0 , $\overline{𝐷}$ 0 mesons produced in 𝑝⁢𝑝 collisions at large Feynman variables 𝑥. It is shown that an important tool for establishing the intrinsic charm (IC) content of the proton is the charm hadron-antihadron asymmetry formed in 𝑝⁢𝑝 collisions. Predictions for the asymmetry as a function of 𝑥 for different IC probabilities are presented. We show that the interference of the intrinsic |𝑢⁢𝑢⁢𝑑⁢𝑐⁢$\overline{𝑐}$⟩ Fock state with the standard contribution from the perturbative QCD evolution leads to a large 𝐷 + ⁢𝐷 − asymmetry at large Feynman 𝑥.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Applications of the Landau bootstrap

We advocate a strategy of bootstrapping Feynman integrals from just knowledge of their singular behavior. This approach is complementary to other bootstrap programs, which exploit nonperturbative constraints such as unitarity, or amplitude-level constraints such as gauge invariance. We begin by studying where a Feynman integral can become singular, and the behavior it exhibits near these singularities. We then characterize the space of functions that we expect the integral to evaluate to, in order to formulate an appropriate ansatz. Finally, we derive constraints on where each singularity can appear in this ansatz, and use information about the expansion of the integral around singular points in order to determine the value of all remaining free coefficients. Throughout, we highlight how constraints that have previously only been derived for integrals with generic masses can be extended to integrals involving particles of equal or vanishing mass. We illustrate the effectiveness of this approach by bootstrapping a number of examples, including the four-point double box with a massive internal loop. Published by the American Physical Society 2025

Hannesdottir, Holmfridur S. (ORCID:000000025440208↗

Demonstration of the rodeo algorithm on a quantum computer

The rodeo algorithm is an efficient algorithm for eigenstate preparation and eigenvalue estimation for any observable on a quantum computer. This makes it a promising tool for studying the spectrum and structure of atomic nuclei as well as other fields of quantum many-body physics. The only requirement is that the initial state has sufficient overlap probability with the desired eigenstate. While it is exponentially faster than well-known algorithms such as phase estimation and adiabatic evolution for eigenstate preparation, it has yet to be implemented on an actual quantum device. In this work, we apply the rodeo algorithm to determine the energy levels of a random one-qubit Hamiltonian, resulting in a relative error of 0.08% using mid-circuit measurements on the IBM Q device Casablanca. This surpasses the accuracy of directly-prepared eigenvector expectation values using the same quantum device. We take advantage of the high-accuracy energy determination and use the Hellmann-Feynman theorem to compute eigenvector expectation values for a different random one-qubit observable. For the Hellmann-Feynman calculations, we find a relative error of 0.7%. Here, we conclude by discussing possible future applications of the rodeo algorithm for multi-qubit Hamiltonians.

algorithm↗

Preliminary RAM-RODD results for the MUSiC subcritical configurations

The Measurement of Uranium Subcritical and Critical (MUSiC) was performed at the DOE’s National Criticality Experiments Research Center (NCERC) located in the Nevada National Security Site (NNSS). The measurement utilized the Rocky Flats shells to perform benchmark measurements of similar highly enriched uranium (HEU) systems that span a wide range of reactivities. The Rocky Flats (RF) shells are 93.16% U-235 enriched metal hemishells that can be stacked concentrically. Ten configurations were measured with effective multiplication factors spanning between deeply subcritical (~ 0.64) through delayed critical. Details of the measured configurations are listed in Table 1. This unique set of measurements with its large span of reactivies is being used to determine the range over which neutron noise techniques such as Feynman variance-to-mean, Rossi-alpha, and pulsed neutron source techniques can be accurately employed for a bare HEU system. The results of the measurements will be published as a benchmark in The International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The results will support the growing amount of subcritical benchmark data, such as the SCRaP measurements, that is available to the community. The measurements were performed using several different detector systems for the purposes of cross-validation and obtaining detector independent results. Four detector systems were deployed, three by the NCERC team and one from the University of Michigan. The detector systems included a Neutron Multiplicity Array Detector (NoMAD) system (similar to the MC-15), four small volume 0.635 cm (Ø) × 7.59 cm 3 He detectors (ideal for measuring prompt neutron decay constants due to their fast recovery speed), the Rossi Alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAM-RODD), and the Organic Scintillator Array (OSCAR). RAM-RODD is an array of eight 5.08 cm (Ø) × 5.08 cm EJ-309 organic scintillator detectors. OSCAR is a University of Michigan system and is an array of twelve 5.08 cm (Ø) × 5.08 cm stilbene detectors. Details of measurements performed with OSCAR will be discussed in a separate talk. The focus of this work is preliminary results obtained by RAM-RODD for the 8 subcritical configurations. Additional details on the measurements and on the setup and deployment of RAM-RODD will be discussed. Preliminary neutron noise analysis results including Rossi-alpha and Feynman-alpha will also be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Energy Systems Process Modeling, Analysis, and Experiment Design [Slides]

This presentation is to be used in Feynman Center outreach discussions. It describes public information about some Los Alamos interests and activities in modeling and analysis in energy systems including carbon capture and direct air capture. In particular this addresses the CCSI2 project’s open-source toolkit and ongoing partnerships, as well as summarizing the published Feynman Center capability snapshot on direct air capture.

97 MATHEMATICS AND COMPUTING↗

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-temperature many-body perturbation theory for electrons: Algebraic recursive definitions, second-quantized derivation, linked-diagram theorem, general-order algorithms, and grand canonical and canonical ensembles

A comprehensive and detailed account is presented for the finite-temperature many-body perturbation theory for electrons that expands in power series all thermodynamic functions on an equal footing. Algebraic recursions in the style of the Rayleigh–Schrödinger perturbation theory are derived for the grand potential, chemical potential, internal energy, and entropy in the grand canonical ensemble and for the Helmholtz energy, internal energy, and entropy in the canonical ensemble, leading to their sum-over-states analytical formulas at any arbitrary order. For the grand canonical ensemble, these sum-over-states formulas are systematically transformed to sum-over-orbitals reduced analytical formulas by the quantum-field-theoretical techniques of normal-ordered second quantization and Feynman diagrams extended to finite temperature. It is found that the perturbation corrections to energies entering the recursions have to be treated as a nondiagonal matrix, whose off-diagonal elements are generally nonzero within a subspace spanned by degenerate Slater determinants. They give rise to a unique set of linked diagrams—renormalization diagrams—whose resolvent lines are displaced upward, which are distinct from the well-known anomalous diagrams of which one or more resolvent lines are erased. A linked-diagram theorem is introduced that proves the size-consistency of the finite-temperature many-body perturbation theory at any order. General-order algorithms implementing the recursions establish the convergence of the perturbation series toward the finite-temperature full-configuration-interaction limit unless the series diverges. As a result, the normal-ordered Hamiltonian at finite temperature sheds light on the relationship between the finite-temperature Hartree–Fock and first-order many-body perturbation theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A multisite decomposition of the tensor network path integrals

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗

Gauge invariance of radiative jet functions in the position-space formulation of SCET

In subleading powers of soft-collinear effective theory (SCET), the Lagrangian contains couplings between soft quarks and hard-collinear quarks. Matrix elements of the hard-collinear parts of these couplings are radiative jet functions. In the position-space formulation of SCET, the Lagrangians are constructed from operators that appear to be gauge invariant. Nevertheless, we find violations of gauge invariance arise in the hard-collinear sector because gauge transformations can shift the momentum of a hard-collinear quark field from the hard-collinear sector to the soft sector, where the hard-collinear fields, by definition, have no support. The violations of gauge invariance are manifested in perturbation theory in the hard-collinear sector through the absence of certain Feynman diagrams that would be present in full QCD. A consequence of the absence of these diagrams is that the radiative jet functions that follow directly from the position-space Lagrangians are not gauge invariant, and we demonstrate this through explicit calculations in lower-order perturbation theory. We obtain gauge-invariant Lagrangians by adding to existing position-space Lagrangians terms that are proportional to the soft-quark equation of motion. These gauge-invariant Lagrangians are valid for nonzero, as well as zero, quark masses. We also remark briefly on the gauge invariance of certain Lagrangians that have been constructed in the label-momentum formulation of SCET. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. In this work, we introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Future of the MUSiC Experiment Data

The Measurements of Uranium Subcritical and Critical (MUSiC) experiment was a highly enriched uranium (HEU) experiment performed at the National Criticality Experiments Research Center (NCERC) executed between December 2020 and April 2021. The experiment intended to measure criticality and reactor kinetics parameters in a bare HEU system. The experiment concurrently measured radiation signatures from the system while utilizing different neutron source types. This was an attempt to benchmark both detectors and analysis techniques against one another for identical measurements. The experiment consisted of the Rocky Flats HEU hemi-shells constructed into ten configurations spanning between deeply subcritical (about 14 kgs) to supercritical (about 60 kgs). Two of the ten configurations were supercritical and the other eight were subcritical. The two supercritical configurations, often referred to as the critical configurations, were documented into an International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmark evaluation and submitted to the technical review group (TRG). The subcritical configurations of the MUSiC experiment were examined using three different detection systems. The systems include: the NoMAD He-3 neutron detector which consists of 15 He-3 tubes surrounded by a polyethylene matrix, a liquid scintillator system called the Rossi-α Measurement Rapid Organic Discriminating Detector (RAM-RODD), and a trans-stilbene organic scintillator array (OSCAR) provided by the University of Michigan. The NoMAD and RAM-RODD data are planned to be evaluated in two separate ICSBEP evaluations utilizing neutron noise methods such as Feynman Variance-to-Mean, Rossi-α and the pulsed neutron source method. Each of the subcritical configurations were examined using three different source types including: a Cf-252 source in the center, with only the intrinsic neutron source in the HEU, and with an external D-T neutron generator.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of the MUSIC 3 He Multiplicity Data

A measurement campaign called the Measurement of Uranium Subcritical and Critical (MUSiC) was performed on a range of configurations of highly-enriched uranium (HEU) from December 2020 through April of 2021. While part of the focus was to measure reactor kinetics parameters on delayed supercritical systems, an additional focus was performing neutron noise measurements on subcritical configurations from deeply subcritical to nearly delayed critical. Multiple detector systems were used to perform these measurements, such as a 3 He multiplicity detector called the Neutron Multiplicity Array Detector (NoMAD) and a liquid scintillator system called the Rossi-α Measurement Rapid Organic Discriminating Detector (RAM-RODD). Also included were a scintillator system from the University of Michigan and a set of small 3He tubes that have previously been used to measure Rossi-α values on near-critical systems. The focus of this paper will be a comparison of prospective analysis methods for the NoMAD measurements. Previous subcritical measurements at the National Criticality Experiments Research Center (NCERC) submitted to the International Criticality Safety Benchmark Evaluation Project (ICSBEP) used the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This relies on the time correlations of neutron detections to infer the spontaneous fission rate and neutron multiplication of a system through binning the time tagged detections and analyzing resulting histograms of the numbers of counts. However, there are other neutron noise methods that rely on similar processes, such as the Hansen-Dowdy formalism which uses a slightly different methodology to extract multiplication from the neutron multiplicity counting moments. Comparisons can be made between these experimental results and those obtained through simulations to validate or identify deficiencies in analysis, detection methods, or the underlying nuclear data. Different time gating strategies and their effects on count rate uncertainties are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High Multiplication Neutron Noise Measurements Using the 7uPCX Assembly

The Seven Percent Critical Experiment, or 7uPCX, is a system that mimics the physics of light water nuclear reactor systems by using uranium dioxide fuel pins at an enrichment of approximately 7%. 7uPCX is often used for benchmarking and nuclear data purposes. As part of a collaboration between Lawrence Livermore National Laboratory, Los Alamos National Laboratory, Sandia National Laboratories, and the Institut de Radioprotection et de Sûreté Nucléaire, measurements were performed in support of a high-multiplication, subcritical benchmark candidate for a thermal system. The measurements were completed by making subcritical configurations using the fuel loading pattern of a well-documented configuration of the Seven Percent Critical Experiment at Sandia National Laboratories, which exists as a benchmark in the International Criticality Safety Benchmark Evaluation Project handbook. These measurements, which aimed to capture multiplications ranging from approximately 10 to 1,000, also serve as an intercomparison between both the fielded detector systems and analysis methodologies with the goal of better characterizing the detectors and their ability to capture the state of the criticality these types of systems. Los Alamos National Laboratory’s measurements for this collaboration were made with five linked helium-3 based neutron multiplicity detectors placed on the periphery of the reactor tank beyond the infinite reflector thickness of water, and four organic scintillators placed in dry-wells in-reactor near the edge of the upper fuel grid plate. This work captures the Los Alamos National Laboratory measurements, the Rossi-α and Feynman-Y results, provides an intercomparison between the results of the 3 He neutron detectors and the organic scintillators, compares to simulation where possible, and expands on future work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗