Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electromagnetic field calculation”

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 217 records · Page 12

Quasielastic electromagnetic scattering cross sections and world data comparisons in the GENIE Monte Carlo event generator

The usage of Monte Carlo neutrino event generators (MC$\textit{ν}$EGs) is a norm within the high-energy $\textit{ν}$ scattering community. The relevance of quasielastic (QE) energy regimes to $\textit{ν}$ oscillation experiments implies that accurate calculations of $\textit{ν}$A cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MC$\textit{ν}$EGs utilize highly phenomenological, parametrized models of QE scattering cross sections. Moreover, a culture of validation of MC$\textit{ν}$EGs against prolific electron ($\textit{e}$) scattering data has been historically lacking. In this work, we implement new eA cross sections obtained from nuclear ab initio approaches in GENIE, the primary MC$\textit{ν}$EG utilized by the FNAL community. In particular, we utilize results from quantum MC methods which solve the many-body nuclear problem in the short-time approximation (STA), allowing consistent retention of two-nucleon dynamics which are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world QE electromagnetic data, finding agreement with available data below ~2 GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study $\textit{A}$ ≤ 12 nuclei, however, its semi-inclusive multibody identity components are exportable to other many-body computational techniques such as auxiliary field diffusion MC which can reach $\textit{A}$ ≤ 40 systems while continuing to realize the factorization contained within the STA’s multinucleon dynamics. Together, these developments promise to make future experiments such as DUNE more accurate in their assessment of MC$\textit{ν}$EG systematics, $\textit{ν}$ properties, and potentially empower the discovery of physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Here, we present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses

The decreasing uncertainties in theoretical predictions and experimental measurements of several hadronic observables related to weak processes, which in many cases are now smaller than O(1%), require theoretical calculations to include subleading corrections that were neglected so far. Precise determinations of leptonic and semi-leptonic decay rates, including QED and strong isospin-breaking effects, can play a central role in solving the current tensions in the first-row unitarity of the CKM matrix. In this work we present the first RBC/UKQCD lattice calculation of the isospin-breaking corrections to the ratio of leptonic decay rates of kaons and pions into muons and neutrinos. The calculation is performed at fixed lattice spacing (a –1 ≃ 1.730 GeV) on a 48 3 × 96 volume with N f = 2 + 1 dynamical quarks close to the physical point and domain wall fermions in the Möbius formulation are employed. Long-distance QED interactions are included according to the QED L prescription and the crucial role of finite-volume electromagnetic corrections in the determination of leptonic decay rates, which produce a large systematic uncertainty, is extensively discussed. Finally, we study the different sources of uncertainty on |V us |/|V ud | and observe that, if finite-volume systematics can be reduced, the error from isospin-breaking corrections is potentially sub-dominant in the final precision of the ratio of the CKM matrix elements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

All-multiplicity one-loop amplitudes in Born-Infeld electrodynamics from generalized unitarity

We initiate a study of non-supersymmetric Born-Infeld electrodynamics in 4d at the quantum level. Explicit all-multiplicity expressions are calculated for the purely rational one-loop amplitudes in the self-dual (+ + … +) and next-to-self-dual (- + … +) helicity sectors. Using a supersymmetric decomposition, d-dimensional unitarity cuts of the integrand factorize into tree-amplitudes in a 4d model of Born-Infeld photons coupled to a massive complex scalar. The two-scalar tree-amplitudes needed to construct the Born-Infeld integrand are computed using two complimentary approaches: (1) as a double-copy of Yang-Mills coupled to a massive adjoint scalar with a dimensionally reduced form of Chiral Perturbation Theory, and (2) by imposing consistency with low-energy theorems under a reduction from 4d to 3d and T-duality. The Born-Infeld integrand is integrated in d = 4 - 2ϵ dimensions at order Ο(ϵ 0 ) using the dimension-shifting formalism. We comment on the implications for electromagnetic duality in quantum Born-Infeld theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

16 O Electroweak Response Functions from First Principles

We present calculations of various electroweak response functions for the 16 O nucleus obtained using coupled-cluster theory in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-leading order and next-to-next-to-leading order in chiral effective field theory and perform a Bayesian analysis to assess uncertainties. Our results are in good agreement with available electron-scattering data at |𝐪|≈326 MeV/c. Additionally, we provide several predictions for the weak response functions in the quasielastic peak region at |𝐪| =300 and 400 MeV/c, which are critical for long-baseline neutrino experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

First Field Measurements of the 15 T Nb 3 Sn Dipole Demonstrator MDPCT1

Within the US Magnet Development Program (MDP) a 15 T Nb 3 Sn dipole was developed, to demonstrate a magnet design for a post-LHC pp Collider. Specifically, the magnet design is based on 60 mm aperture 4-layer shell-type coils, graded between the inner and outer layers to maximize the magnet performance. The cable in the two innermost layers has 28 strands 1.0 mm in diameter and the cable in the two outermost layers has 40 strands 0.7 mm in diameter. Both cables use RRP Nb3Sn wires produced by Bruker-OST. Magnet coils are surrounded by vertically-split thick iron laminations, con-nected by aluminum I-clamps, and a thick stainless-steel skin. The magnet was tested at the Vertical Magnet Test Facility (VMTF) at Fermilab. This paper reports the first results of magnetic measure-ments of the 15 T Nb 3 Sn dipole demonstrator including geometrical harmonics, coil magnetization and iron saturation effects. The ex-perimental data are compared with the magnetic calculations.

43 PARTICLE ACCELERATORS↗

Simultaneous Observations of Electromagnetic Ion Cyclotron (EMIC) Waves and Pitch Angle Scattering During a Van Allen Probes Conjunction

On 22 December 2015, the two Van Allen Probes observed two sets of electromagnetic ion cyclotron (EMIC) wave bursts during a close conjunction when both Probe A and Probe B were separated by 0.57 to 0.68 R E . Here, the EMIC waves occurred during an active period in the recovery phase of a coronal mass ejection-driven geomagnetic storm. Both spacecraft observed EMIC wave bursts that had similar spatial structure within a 1–2 min time delay. The EMIC waves occurred outside the plasmasphere, within ΔL ≈ 1–2 of the plasmapause and within a few degrees in magnetic latitude of the equatorial plane. The spatial structure of the EMIC wave bursts may have been related to the proton drift paths outside the plasmasphere and influenced by total magnetic field strength variations associated with solar wind pressure enhancements. The EMIC waves were observed in a narrow L shell region from L ≈ 4.55–5.32 between 10 and 11 magnetic local time (MLT) on the outbound halves of the spacecraft orbits and from L ≈ 4.82–5.51 between 13 and 14 MLT on the inbound halves of the spacecraft orbits. However, Pc1 pulsations were observed on the ground over a broad range of local times. The anisotropy of the proton pitch angle distributions was enhanced when the EMIC waves were observed. Although the overall radiation belt response during this storm was dominated by acceleration and transport processes, the EMIC waves produced local pitch angle scattering of 13–15 keV protons and 2.1–2.6 MeV electrons, consistent with calculations of the expected resonant energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE↗

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗

Recent progress in low energy neutrino scattering physics and its implications for the standard and beyond the standard model physics

Neutrinos continue to provide a testing ground for the structure of the standard model of particle physics as well as hints towards the physics beyond the standard model. Neutrinos of energies spanning over several orders of magnitude, originating in many terrestrial and astrophysical processes, have been detected via various decay and interaction mechanisms. At MeV scales, there has been one elusive process, until a few years ago, known as coherent elastic neutrino-nucleus scattering (CEvNS) that was theoretically predicted over five decades ago but was never observed experimentally. The recent experimental observation of the CEvNS process by the COHERENT collaboration at a stopped pion neutrino source has inspired physicists across many subfields. This has vital implications for nuclear physics, high-energy physics, astrophysics, and beyond. CEvNS, being a low-energy process, provides a natural window to study light, weakly-coupled, new physics in the neutrino sector. In this review, we intend to provide the current status of low energy neutrino scattering physics and its implications for the standard and beyond the standard model physics. We discuss the general formalism of calculating the tree-level CEvNS cross section and present estimated theoretical uncertainties on the CEvNS cross section stemming from different sources. We also discuss the inelastic scattering of tens of MeV neutrinos that have implications for supernova detection in future neutrino experiments. We discuss how the CEvNS experiments can be used as a testing ground for the Standard Model (SM) weak physics as well as in searching for the Beyond the Standard Model (BSM) physics signals. Any deviation from the SM predicted event rate either with a change in the total event rate or with a change in the shape of the recoil spectrum, could indicate new contributions to the interaction cross-section. The SM implications include the study of weak nuclear form factor and weak mixing angle. The BSM studies include non-standard interactions, neutrino electromagnetic properties, and sterile neutrino searches. In conclusion, stopped pion facilities are also a copious source of neutral and changed mesons that allow study of several dark sector physics scenarios such as vector portal models, leptophobic dark matter as well as axion-like particle searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The photon content of the proton in the CT18 global analysis

Recently, two photon PDF sets based on implementations of the LUX ansatz into the CT18 global analysis were released. In CT18lux, the photon PDF is calculated directly using the LUX master formula for all scales, \mu μ . In an alternative realization, CT18qed, the photon PDF is initialized at the starting scale, \mu_0 μ 0 , using the LUX formulation and evolved to higher scales \mu(>\mu_0) μ ( > μ 0 ) with a combined QED+QCD kernel at \mathcal{O}(\alpha),~\mathcal{O}(\alpha\alpha_s) 𝒪 ( α ) , 𝒪 ( α α s ) and \mathcal{O}(\alpha^2) 𝒪 ( α 2 ) . In the small- x x region, the photon PDF uncertainty is mainly induced by the quark and gluon PDFs, through the perturbative DIS structure functions. In comparison, the large- x x photon uncertainty comes from various low-energy, nonperturbative contributions, including variations of the inelastic structure functions in the resonance and continuum regions, higher-twist and target-mass corrections, and elastic electromagnetic form factors of the proton. We take the production of doubly-charged Higgs pairs, (H^{++}H^{--}) ( H + + H − − ) , as an example of scenarios beyond the Standard Model to illustrate the phenomenological implications of these photon PDFs at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Photoelectron–residual-ion interaction in angle-resolved streaked shake-up ionization of helium

Streaked photoelectron emission spectra access the correlated dynamics of photoelectrons and residual target electrons with attosecond temporal resolution. Here, we calculated ab initio single-ionization spectra for photoemission from helium atoms by colinearly polarized ultrashort XUV and assisting few-femtosecond IR pulses. Distinguishing direct and shake-up ionization resulting in ground-state and excited (𝑛 = 2, 3) He + residual ions, respectively, we examined the effects of the correlated photoemission dynamics on the photoelectron phase accumulation as a function of the observable photoelectron detection direction and kinetic energy and XUV–IR pulse delay. We tracked the dynamical evolution of the residual ion in relative streaked photo- emission delays and found dominant contributions for shake-up emission from the residual-ion–photoelectron interaction. These are in very good and fair agreement, respectively, for 𝑛 = 2 and 𝑛 = 3 shake-up photoemission along the pulse-polarization directions, with previous experimental and theoretical investigations [Ossiander et al ., Nat. Phys. 13, 280 (2017)] and reveal a strong photoemission-direction dependence for shake-up ionization due to the coupling between the photoelectron and evolving residual-ion charge distribution in the IR-laser field.

Stark effect↗

Electromagnetic duality and D3-brane scattering amplitudes beyond leading order

We use on-shell methods to study the non-supersymmetric and supersymmetric low-energy S-matrix on a probe D3-brane, including both the 1-loop contributions of massless states as well as the effects of higher-derivative operators. Our results include: (1) A derivation of the duality invariance of Born-Infeld electrodynamics as the dimensional oxidation of the group of spatial rotations transverse to a probe M2-brane; this is done using a novel implementation of subtracted on-shell recursion. (2) The first explicit loop-level BCJ double-copy in a non-gravitational model, namely the calculation of the 4-point self-dual amplitude of non-supersymmetric Born-Infeld. (3) From previous results for n-point self-dual 1-loop BI amplitudes and the conjectured dimension-shifting relations in Yang-Mills, we obtain an explicit all-multiplicity, at all orders in E, expression for the 1-loop integrand of the MHV sector of N = 4 DBI. (4) For all n > 4, the explicitly integrated duality-violating 1-loop amplitudes (self-dual and next-to-self-dual in pure BI as well as MHV in N = 4 DBI) are shown to be removable at O(ε 0 ) by adding finite local counterterms; we propose that this may be true more generally at 1-loop order. (5) We find that in non-supersymmetric Born-Infeld, not all finite local counterterms needed to restore electromagnetic duality can be constructed using the double-copy with higher-derivative corrections, suggesting a fundamental tension between electromagnetic duality and color-kinematics duality at loop-level. Finally we comment on oxidation of duality symmetries in supergravity and the parallels it has to the M2-brane to D3-brane oxidation demonstrated in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors $G^c_{E,M}(Q^2)$ in the momentum transfer range $0\leq Q^2 \leq 1.4$ $\rm GeV^2$. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment $\mu^c_M=-0.00127(38)_{\rm stat}(5)_{\rm sys}$, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero $G^c_{E}(Q^2)$ indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the $[c(x)-\bar{c}(x)]$ distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Plasmon-induced excitation energy transfer in silver nanoparticle dimers: A real-time TDDFTB investigation

Here, using real-time quantum dynamics calculations, we perform theoretical investigations of light-induced interactions and electronic excitation transfer in a silver nanoparticle dimer. Real-time time-dependent density functional tight-binding (RT-TDDFTB) calculations provide details of the quantum dynamical processes at an electronic/atomistic level with attosecond resolution. The computational efficiency of RT-TDDFTB allows us to examine electronic dynamics up to picosecond time scales. With time scales varying over six orders of magnitude, we provide insight into interactions between the nanoparticle and laser and between nanoparticles. Our results show that the coupling between nanoparticle monomers is dependent on the separation distance between the nanoparticles in the dimer. As the interparticle distance is varied, the dipole–dipole interactions and electronic excitation transfer mechanisms are markedly different. At large distances (from 50 to 20 Å), the energy transfer from NP1 to NP2 becomes more efficient as the interparticle distance decreases. The total dipole moment of the Ag 14 nanoparticle dimer increases linearly at an interparticle distance of 20 Å and reaches its maximum after 1.2 ps. The electronic excitation transfer is also the most efficient at 20 Å. At short distances, back-transfer effects reduce the ability of the dimer and NP1 to accept energy from the incident electric field. We attribute the distance-dependent features of the nanoparticle dimer to the beating between the laser acting on NP1 and the back transfer from NP2 to NP1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Precision Spectroscopy of 20 O Benchmarking Ab Initio Calculations in Light Nuclei

The excited states of unstable 20 O were investigated via γ-ray spectroscopy following the 19 O(d, p) ⁢20 O reaction at 8 A⁢ MeV. By exploiting the Doppler shift attenuation method, the lifetimes of the $2^+_2$ and $3^+_1$ states were firmly established. From the γ-ray branching and E⁢2/M⁢1 mixing ratios for transitions deexciting the $2^+_2$ and $3^+_1$ states, the B⁡(E⁢2) and B⁡(M⁢1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of γ-ray transition probabilities with the valence space in medium similarity renormalization group ab initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab initio approaches using chiral effective field theory forces are challenged by detailed high-precision spectroscopic properties of nuclei. Finally, the reduced transition probabilities were found to be a very constraining test of the performance of the ab initio models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-precision scale setting with the Ω-baryon mass and gradient flow

The gradient-flow scale 𝑤 0 in lattice QCD is determined using the mass of the Ω − baryon to set the physical scale. Nine ensembles using the highly improved staggered quark (HISQ) action with lattice spacings of 0.15 fm down to 0.04 fm are used, seven of which have nearly physical light-quark masses. Electromagnetic corrections to the Ω − mass are defined in order to compute a pure-QCD Ω mass. The final result is 𝑤 0 = 0.17187⁢(68) fm, corresponding to a relative uncertainty of 0.40% and a central value in good agreement with previous calculations in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗