Engineering Papers⌕ Search

Engineering topics

Christensen, A.

Publications and source records attributed to Christensen, A..

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗
DOE OSTI2025

Joint neutrino oscillation analysis from the T2K and NOvA experiments

The landmark discovery that neutrinos have mass and can change type (or flavour) as they propagate—a process called neutrino oscillation—has opened up a rich array of theoretical and experimental questions being actively pursued today. Neutrino oscillation remains the most powerful experimental tool for addressing many of these questions, including whether neutrinos violate charge-parity (CP) symmetry, which has possible connections to the unexplained preponderance of matter over antimatter in the Universe. Oscillation measurements also probe the mass-squared differences between the different neutrino mass states (Δ$m^2$), whether there are two light states and a heavier one (normal ordering) or vice versa (inverted ordering), and the structure of neutrino mass and flavour mixing. Here we carry out the first joint analysis of datasets from NOvA and T2K, the two currently operating long-baseline neutrino oscillation experiments (hundreds of kilometres of neutrino travel distance), taking advantage of our complementary experimental designs and setting new constraints on several neutrino sector parameters. This analysis provides new precision on the $Δm_{32}^2$ mass difference, finding $-2.43_{−0.03}^{+0.04}$ x $10^{-13}$ $\textrm{eV}^2$ in the normal ordering and $-2.43^{+0.03}_{-0.04}$ x $10^{-13}$ $\textrm{eV}^2$ in the inverted ordering, as well as a 3σ interval on δ CP of [−1.38π, 0.30π] in the normal ordering and [−0.92π, −0.04π] in the inverted ordering. The data show no strong preference for either mass ordering, but notably, if inverted ordering were assumed true within the three-flavour mixing model, then our results would provide evidence of CP symmetry violation in the lepton sector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Monte Carlo method for constructing confidence intervals with unconstrained and constrained nuisance parameters in the NOvA experiment

Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. Wilks' theorem provides a simple way to construct confidence intervals on model parameters, but it only applies under certain conditions. These conditions, such as nested hypotheses and unbounded parameters, are often violated in neutrino oscillation measurements and other experimental scenarios. Monte Carlo methods can address these issues, albeit at increased computational cost. In the presence of nuisance parameters, however, the best way to implement a Monte Carlo method is ambiguous. Furthermore, this paper documents the method selected by the NOvA experiment, the profile construction. It presents the toy studies that informed the choice of method, details of its implementation, and tests performed to validate it. It also includes some practical considerations which may be of use to others choosing to use the profile construction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Search for 𝐶⁢𝑃-Violating Neutrino Nonstandard Interactions with the NOvA Experiment

This Letter reports a search for charge-parity (𝐶⁢𝑃) symmetry violating nonstandard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from 𝜈 𝜇 ⁡($\bar{𝜈}$ 𝜇 ) → 𝜈 𝜇 ⁡($\bar{𝜈}$ 𝜇 ) and 𝜈 𝜇 ⁡($\bar{𝜈}$ 𝜇 ) → 𝜈 𝑒 ⁡( $\bar{𝜈}$ 𝑒 ) oscillation channels are used to measure the effect of the NSI parameters 𝜖 𝑒⁢𝜇 and 𝜖 𝑒⁢𝜏 . With 90% CL the magnitudes of the NSI couplings are constrained to be |𝜖 𝑒⁢𝜇 | ≲0.3 and |𝜖 𝑒⁢𝜏 | ≲ 0.4. A degeneracy at |𝜖 𝑒⁢𝜏 | ≈ 1.8 is reported, and we observe that the presence of NSI limits sensitivity to the standard 𝐶⁢𝑃 phase 𝛿 𝐶⁢𝑃 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Measurement of the double-differential cross section of muon-neutrino charged-current interactions with low hadronic energy in the NOvA Near Detector

The NOvA collaboration reports cross-section measurements for $\nu_{\mu}$ charged-current interactions with low hadronic energy (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. The results are presented as a double-differential cross section as a function of the direct observables of the final-state muon kinematics. Results are also presented as a single-differential cross section as a function of the derived square of the four-momentum transfer, $Q^{2}$, and as a function of the derived neutrino energy. The data correspond to an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino-nucleus interactions with an enhanced fraction of quasi-elastic and two-particle-two-hole (2p2h) interactions. This enhancement allows quantitative comparisons with various nuclear models. We find strong disagreement between data and theory-based models in various regions of the muon kinematic phase space, especially in the forward muon direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Expanding neutrino oscillation parameter measurements in NOvA using a Bayesian approach

NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current ν μ → ν μ (disappearance) and ν μ → ν e (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [] using an alternative statistical approach based on Bayesian Markov chain Monte Carlo. We measure oscillation parameters consistent with the previous results. We also extend our inferences to include the first NOvA measurements of the reactor mixing angle θ 13 , where we find 0.071 ≤ sin 2 2 θ 13 ≤ 0.107 , and the Jarlskog invariant, where we observe no significant preference for the C P -conserving value J = 0 over values favoring C P violation. We use these results to examine the effects of constraints from short-baseline measurements of θ 13 using antineutrinos from nuclear reactors when making NOvA measurements of θ 23 . Our long-baseline measurement of θ 13 is shown to be consistent with the reactor measurements, supporting the general applicability and robustness of the Pontecorvo-Maki-Nakagawa-Sakata framework for neutrino oscillations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Observation of the effect of gravity on the motion of antimatter

Einstein’s general theory of relativity from 1915 remains the most successful description of gravitation. From the 1919 solar eclipse to the observation of gravitational waves, the theory has passed many crucial experimental tests. However, the evolving concepts of dark matter and dark energy illustrate that there is much to be learned about the gravitating content of the universe. Singularities in the general theory of relativity and the lack of a quantum theory of gravity suggest that our picture is incomplete. It is thus prudent to explore gravity in exotic physical systems. Antimatter was unknown to Einstein in 1915. Dirac’s theory appeared in 1928; the positron was observed in 1932. There has since been much speculation about gravity and antimatter. The theoretical consensus is that any laboratory mass must be attracted by the Earth, although some authors have considered the cosmological consequences if antimatter should be repelled by matter. In the general theory of relativity, the weak equivalence principle (WEP) requires that all masses react identically to gravity, independent of their internal structure. Here we show that antihydrogen atoms, released from magnetic confinement in the ALPHA-g apparatus, behave in a way consistent with gravitational attraction to the Earth. Repulsive ‘antigravity’ is ruled out in this case. This experiment paves the way for precision studies of the magnitude of the gravitational acceleration between anti-atoms and the Earth to test the WEP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of substrate and growth method on vanadium dioxide thin films by RF magnetron sputtering: Vanadium metal oxidation vs reactive sputtering

Vanadium dioxide (VO 2 ) undergoes a metal–insulator phase transition at ~70 °C and has attracted substantial interest for potential applications in electronics, including those in neuromorphic computing. The vanadium–oxygen system has a rather complicated phase diagram, and controlling the stoichiometry and the phase of thin films of vanadium oxides is a well-known challenge. Here, we explore the novel combination of two methods of VO 2 thin film deposition using off-axis RF magnetron sputtering on (100)- and (111)-oriented yttria-stabilized zirconia (YSZ) substrates: reactive sputtering of vanadium in an oxygen environment and sputtering of vanadium metal followed by oxidation to VO 2 . Interestingly, the reactive sputtering process on both substrate orientations yields the metastable semiconducting VO 2 (B) phase, which is structurally stabilized by the YSZ surface. The metal sputtering and oxidation process on YSZ produces mainly the equilibrium monoclinic (or M1) phase of VO 2 that exhibits a metal–insulator transition. Using this method, we obtained thin films of (010)-textured polycrystalline VO 2 (M1) that show a metal–insulator transition with an on/off ratio larger than 1000.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Deterministic transport codes play a fundamental role in the modelling and simulation of neutron transport. One of the most common deterministic methods is the method of discrete ordinates, also known as the S method. While offering significant advantages over other deterministic methods or stochastic methods like Monte Carlo, the method of discrete ordinates suffers from non-physical artifacts in its local solution due to its discretization of angle. These artifacts, referred to as ray effects because of their ray-like appearance, tend to be worse in problems with small sources in areas with little scattering. Significant effort has gone into developing methods to mitigate ray effects, such as the first-collision source treatment, which separates the angular flux into the uncollided and collided fluxes and solving them using non-traditional techniques such as ray tracing. One such code capable of doing this is Lawrence Livermore National Laboratory's deterministic transport code ARDRA. Current ray tracing methods typically trace to a set of points inside a zone to compute an overall flux. However, this approach has significant drawbacks, such as a low order of convergence and not being conservative. Therefore, a new method has been developed that traces instead to a set of points on each of a zone's surfaces and computing the currents, before using these to obtain the flux. A comparison between these two ray tracing methods showed significant advantages to the new surface method, including inherent conservation, a higher convergence rate, and an increase in calculable information like leakage. This work performed under the auspices of the U.S. Department of Energy by Lawrence Liver- more National Laboratory under Contract DE-AC52-07NA27344. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Laser cooling of antihydrogen atoms

The photon—the quantum excitation of the electromagnetic field—is massless but carries momentum. A photon can therefore exert a force on an object upon collision 1 . Slowing the translational motion of atoms and ions by application of such a force 2,3 , known as laser cooling, was first demonstrated 40 years ago 4,5 . It revolutionized atomic physics over the following decades 6–8 , and it is now a workhorse in many fields, including studies on quantum degenerate gases, quantum information, atomic clocks and tests of fundamental physics. However, this technique has not yet been applied to antimatter. Here we demonstrate laser cooling of antihydrogen 9 , the antimatter atom consisting of an antiproton and a positron. By exciting the 1S–2P transition in antihydrogen with pulsed, narrow-linewidth, Lyman-α laser radiation 10,11 , we Doppler-cool a sample of magnetically trapped antihydrogen. Although we apply laser cooling in only one dimension, the trap couples the longitudinal and transverse motions of the anti-atoms, leading to cooling in all three dimensions. We observe a reduction in the median transverse energy by more than an order of magnitude—with a substantial fraction of the anti-atoms attaining submicroelectronvolt transverse kinetic energies. We also report the observation of the laser-driven 1S–2S transition in samples of laser-cooled antihydrogen atoms. The observed spectral line is approximately four times narrower than that obtained without laser cooling. The demonstration of laser cooling and its immediate application has far-reaching implications for antimatter studies. A more localized, denser and colder sample of antihydrogen will drastically improve spectroscopic 11–13 and gravitational 14 studies of antihydrogen in ongoing experiments. Furthermore, the demonstrated ability to manipulate the motion of antimatter atoms by laser light will potentially provide ground-breaking opportunities for future experiments, such as anti-atomic fountains, anti-atom interferometry and the creation of antimatter molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Electron cyclotron resonance (ECR) magnetometry with a plasma reservoir

The local magnetic field in a Penning–Malmberg trap is found by measuring the temperatures that result when electron plasmas are illuminated by microwave pulses. Multiple heating resonances are observed as the pulse frequencies are swept. The many resonances are due to electron bounce and plasma rotation sidebands. The heating peak corresponding to the cyclotron frequency resonance is identified to determine the magnetic field. A new method for quickly preparing low density electron plasmas for destructive temperature measurements enables a rapid and automated scan of microwave frequencies. This technique can determine the magnetic field to high precision, obtaining an absolute accuracy better than 1 ppm and a relative precision of 26 ppb. One important application is in situ magnetometry for antihydrogen-based tests of charge-parity-time symmetry and of the weak equivalence principle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

The Dipper Satellite: A Medium-Class Explorer Mission to the Threshold of Space

The Dipper satellite will carry out an unprecedented, systematic, and focused in-situ exploration of the Earth's lower ionosphere and thermosphere below 200 km that will produce a pivotal base of knowledge that will significantly advance our understanding of knowledge that will significantly advance our understanding of how our near-space environment works. The satellite will carry comprehensive in-situ probes to measure vector electric and magnetic fields, plasma density and temperature, ion velocities, ion and neutral composition and winds, energetic particles including suprathermal electrons, gravity waves, and lightning bursts. The satellite will include a propulsion system and tapered body that will provide over 10,000 excursions to altitudes below 200 km with over 3000 dips to altitudes below 150 km. With this instrument complement, spacecraft, and orbit, the Dipper mission will gather the necessary combined electrodynamics and neutral dynamics measurements to provide an understanding of the Earth's critical boundary region where the ionized gases of space and the neutral gases of the atmosphere are coupled, and where impinging forces and momentum are deposited from the magnetosphere above and from the troposphere, stratosphere, and mesosphere below. In exploring those physical processes in the lower ionosphere which can only be measured in-situ, the Dipper mission addresses four main science objectives. The Dipper will: 1) reveal how ion-neutral coupling creates a global system of dynamo electric fields and currents; 2) provide first-hand understanding of how magnetospheric currents close in the ionosphere and reveal the effects on the upper atmosphere of magnetospheric energy and momentum deposition; 3) discover the degree of upwards coupling and energy deposition due to thunderstorm electric fields and determine their significance; 4) determine the dynamics and composition of the Earth's lower thermosphere, including its response to gravity, tidal, and planetary waves on a range of spatial scales. A proposal to design, build, operate, and analyze data from instruments on the Dipper spacecraft within the schedule and budget constraints of NASA's MIDEX program was submitted to NASA in 1998. This presentation summarizes the main features of the mission.

Pfaff, R. F., Jr.↗
NASA NTRS1995

Atmospheric response in aurora experiment: Observations of E and F region neutral winds in a region of postmidnight diffuse aurora

The goal of the Atmospheric Response in Aurora (ARIA) experiment carried out at Poker Flat, Alaska, on March 3, 1992, was to determine the response of the neutral atmosphere to the long-lived, large-scale forcing that is characteristic of the diffuse aurora in the post midnight sector. A combination of chemical release rocket wind measurements, instrumented rocket composition measurements, and ground-based optical measurements were used to characterize the response of the neutral atmosphere. The rocket measurements were made at the end of a 90-min period of strong Joule heating. We focus on the neutral wind measurements made with the rocket. The forcing was determined by running the assimilated mapping of ionospheric electrodynamics (AMIE) analysis procedure developed at the National Center for Atmospheric Research. The winds expected at the latitude and longitude of the experiment were calculated using the spectral thermospheric general circulation model developed at the Danish Meteorological Institute. Comparisons of the observations and the model suggest that the neutral winds responded strongly in two height ranges. An eastward wind perturbation of approximately 100 m/s developed between 140 and 200 km altitude with a peak near 160 km. A southwestward wind with peak magnitude of approximately 150 m/s developed near 115 km altitude. The large amplitude winds at the lower altitude are particularly surprising. They appear to be associated with the upward propagating semidiurnal tide. However, the amplitude is much larger than predicted by any of the tidal models, and the shear found just below the peak in the winds was nominally unstable with a Richardson number of approximately 0.08.

Larsen, M. F.↗