Blade: A package for block-triangular form improved Feynman integrals decomposition
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.
In conventional reactor noise measurement based on neutron detection, a neutron detector must be placed as close to a core as possible, to observe a largest possible correlation amplitude. For most of reactors other than critical assemblies, however, the placement is hardly possible. We focused on an alternative reactor noise measurement using prompt gamma-rays having a long flight range, to overcome the above constraint. A Feynman-α analysis using a NaI(Tl) scintillation detector was previously performed in UTR-KINKI, however, white noise originating from radioactivation of the scintillator significantly degraded the correlation amplitude. In this study, another Feynman-α analysis using a BGO gamma-ray detector, which never suffers from any radioactivation, is performed. A series of Feynman-α analyses for a critical state and several subcritical states of UTR-KINKI reactor has been carried out using two BGO detectors and a BF{sub 3} neutron detector to determine prompt-neutron decay constant. In the present analysis, the dead-time effect of these detectors and the delayed-neutron effect are considered. The prompt-neutron decay constant determined from Feynman-α analysis based on gamma-ray detection agrees with that from the analysis based on neutron detection. The Feynman-α analysis based on detecting gamma rays appears to be able to determine prompt-neutron decay constant at least as efficiently as the conventional analysis detecting neutrons.
We present a detailed description of the recent idea for a direct decomposition of Feynman integrals onto a basis of master integrals by projections, as well as a direct derivation of the differential equations satisfied by the master integrals, employing multivariate intersection numbers. We discuss a recursive algorithm for the computation of multivariate intersection numbers, and provide three different approaches for a direct decomposition of Feynman integrals, which we dub the straight decomposition, the bottom-up decomposition, and the top-down decomposition. These algorithms exploit the unitarity structure of Feynman integrals by computing intersection numbers supported on cuts, in various orders, thus showing the synthesis of the intersection-theory concepts with unitarity-based methods and integrand decomposition. We perform explicit computations to exemplify all of these approaches applied to Feynman integrals, paving a way towards potential applications to generic multi-loop integrals.
We propose a natural, parameter-free, discrete-variable formulation of Feynman path integrals. We show that for discrete-variable quantum systems, Feynman path integrals take the form of walks on the graph whose weighted adjacency matrix is the Hamiltonian. By working out expressions for the partition function and transition amplitudes of discretized versions of continuous-variable quantum systems, and then taking the continuum limit, we explicitly recover Feynman's continuous-variable path integrals. We also discuss the implications of our result.
We complete the proof of “Feynman rules” for constructing M-point conformal blocks with external and internal scalars in any topology for arbitrary M in any spacetime dimension by combining the rules for the blocks (based on their Witten diagram interpretation) with the rules for the construction of conformal cross ratios (based on the OPE and “flow diagrams”). The full set of Feynman rules leads to blocks as power series of the hypergeometric type in the conformal cross ratios. We then provide a proof by recursion of the Feynman rules which relies heavily on the first Barnes lemma and the decomposition of the topology of interest in comb structures. Finally, we provide a nine-point example to illustrate the rules.
The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.
The effective delayed neutron fraction, β eff , is a kinematic parameter describing the contribution of delayed neutrons to the effective multiplication factor, $\kappa$ eff , and is of great interest to the nuclear criticality safety community and beyond. In fact, $\kappa$ eff presently must be inferred (as opposed to directly measured) since variables like β eff must be obtained from simulation, reference tables, or a dynamic or pulsed measurement where the assembly of fissile material is perturbed and not in steady state. Several recent works have proposed neutron noise techniques to estimate β eff from static measurements (such as benchmark experiments). This work experimentally investigates an approach based on looking at late-time (up to thousands of seconds) correlated neutrons in Feynman histograms. This work utilizes measurements of weapons grade plutonium and highly enriched uranium to construct late-time Feynman histograms and investigate a proposed methodology of estimating β eff from a fit of the distribution. Data is shown from both He-3 detectors and an organic scintillator array. This work is preliminary in nature, demonstrating some experimental feasibility and proposing improvements for future investigations.
We present a loop-by-loop method for computing the differential equations of Feynman integrals using the recently developed dual form formalism. We give explicit prescriptions for the loop-by-loop fibration of multi-loop dual forms. Then, we test our formalism on a simple, but non-trivial, example: the two-loop three-mass elliptic sunrise family of integrals. We obtain an ε-form differential equation within the correct function space in a sequence of relatively simple algebraic steps. In particular, none of these steps relies on the analysis of q-series. Then, we discuss interesting properties satisfied by our dual basis as well as its simple relation to the known ε-form basis of Feynman integrands. The underlying K3-geometry of the three-loop four-mass sunrise integral is also discussed. Finally, we speculate on how to construct a “good” loop-by-loop basis at three-loop.
Richard Feynman was fresh out of Princeton University’s doctoral program when recruited to assist in the creation of the atomic bomb at Los Alamos. In 1943, Lab director J. Robert Oppenheimer wrote that the 24-year-old was, “by all odds the most brilliant young physicist here, and everyone knows this.” Feynman attempted to live and examine life in a state of play and, as such, it’s only fitting to take a look back at the Nobel Prize winning scientist’s years at the wartime Lab as April Fools’ Day approaches.
The Richard P. Feynman Center for Innovation helps businesses and organizations work with Los Alamos National Laboratory to turn research into practical products and solutions. We connect industry partners with new technologies, research expertise, and licensing opportunities. From protecting inventions to building collaborations, the Feynman Center makes it easier to bring scientific advances to the marketplace and create lasting impact.
In this article we review, for a mathematical audience, the computation of (tree-level) scattering amplitudes in Yang-Mills theory in detail. In particular we demonstrate explicitly how the same formulas for six-particle NMHV helicity amplitudes are obtained from summing Feynman diagrams and from computing the canonical form of the n=6, k=1, m=4 amplituhedron.
We generalize the relation between discontinuities of scattering amplitudes and cut diagrams to cover sequential discontinuities (discontinuities of discontinuities) in arbitrary momentum channels. The new relations are derived using time-ordered perturbation theory, and hold at phase-space points where all cut momentum channels are simultaneously accessible. As part of this analysis, we explain how to compute sequential discontinuities as monodromies and explore the use of the monodromy group in characterizing the analytic properties of Feynman integrals. We carry out a number of cross-checks of our new formulas in polylogarithmic examples, in some cases to all loop orders.
We perform a Feynman diagram calculation of the two-loop scattering amplitude for gravitationally interacting massive particles in the classical limit. Conveniently, we are able to sidestep the most taxing diagrams by exploiting the test-particle limit in which the system is fully characterized by a particle propagating in a Schwarzschild spacetime. We assume a general choice of graviton field basis and gauge fixing that contains as a subset the well-known deDonder gauge and its various cousins. As a highly nontrivial consistency check, all gauge parameters evaporate from the final answer. Moreover, our result exactly matches that of Bern et al., here verified up to sixth post-Newtonian order while also reproducing the same unique velocity resummation at third post-Minkowksian order.
The Hellmann–Feynman (HF) theorem provides a way to compute forces directly from the electron density, enabling efficient force calculations for large systems through machine learning (ML) models for the electron density. The main issue holding back the general acceptance of the HF approach for atom-centered basis sets is the well-known Pulay force which, if naively discarded, typically constitutes an error upward of 10 eV/Å in forces. In this work, we demonstrate that if a suitably augmented Gaussian basis set is used for density functional calculations, the Pulay force can be suppressed, and HF forces can be computed as accurately as analytical forces with state-of-the-art basis sets, allowing geometry optimization and molecular dynamics to be reliably performed with HF forces. Our results pave a clear path forward for the accurate and efficient simulation of large systems using ML densities and the HF theorem.
A recent experiment by D. Bhetuwal et al. [Phys. Rev. Lett. 126, 082301 (2021)] used the (e,e'p) reaction on 12 C to search for the effects of color transparency (the absence of final-state interactions). Color transparency was said to be ruled out. Observing the effects of color transparency depends on the ability of a putative point-like-configuration (PLC), formed in a high-momentum transfer coherent reaction, to escape the nucleus without expanding its size. Here, we study the expansion aspect of color transparency using superconformal baryon-meson symmetry and light-front holographic QCD. A new formalism is obtained and used to analyze the recent experiment. The resulting conclusion is that effects of expansion would not be sufficiently significant in causing final-state interactions to occur. Therefore, we conclude that a PLC was not formed. In conclusion, this means that the Feynman mechanism involving virtual photon absorption on a single high-momentum quark is responsible for the high-momentum electromagnetic form factor of the proton.
We reformulate the analysis of singularities of Feynman integrals in a way that can be practically applied to perturbative computations in the standard model in dimensional regularization. After highlighting issues in the textbook treatment of Landau singularities, we develop an algorithm for classifying and computing them using techniques from computational algebraic geometry. We introduce an algebraic variety called the principal Landau determinant, which captures the singularities even in the presence of massless particles or UV/IR divergences. We illustrate this for 114 example diagrams, including a cutting-edge 2-loop 5-point nonplanar QCD process with multiple mass scales. Published by the American Physical Society 2024
We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.