Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deuteron”

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 199 records · Page 11

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗

A Short-Lived Isotope Counting System (SLICS) to measure reaction cross sections producing light, beta-decaying nuclei on OMEGA

A new activation detector optimized for short-lived isotopes was designed, installed, and tested on OMEGA. The Short-Lived Isotope Counting System (SLICS) can be positioned as close as 35 cm from the OMEGA target chamber center in a ten-inch-manipulator. Stable target samples can be mounted on the front end of the detector assembly to undergo triton-, deuteron-, or neutron-induced reactions that produce short-lived reaction products that beta decay. The emitted betas are detected in a phoswich configuration to eliminate the gamma background. The detector is optimized for isotopes with half-lives from 150 ms to 20 s, filling a gap in OMEGA’s current activation detectors for isotopes with half-lives >5 min. In a pilot study, SLICS achieved an efficiency of 24.5% ± 2.0% for counting 8Li, with a total of 505.8 ± 4.4 qualifying events. Using the existing triton laser-ion acceleration platform on OMEGA, this detector can be used to measure rarely studied tritium-induced reactions on light nuclei with an estimated uncertainty of 12%. As a result, the design of the detector and its performance during a first experiment employing deuterium as a surrogate will be discussed.

Beta spectrometry↗

Evaluation of a deuterated detector for gamma-ray beam intensity measurements

Laser-Compton scattering sources of quasi-monochromatic γ-ray beams have been in operation for over 30 years, with new facilities being under construction or proposed. To sustain a high-quality experimental physics program, continuous monitoring of these beams is essential by measuring the spectral, temporal, and spatial characteristics. In this work, we introduce a new method for determining the γ-ray beam intensity by using photodisintegration of deuterons in an active target detector approach. Measurements were performed using deuterated liquid scintillator detectors at several incident gamma energies between 4 and 13 MeV at the High Intensity γ-ray Source, Durham, USA. The results demonstrate that an instrument based on a deuterated active detector can be used to monitor the intensity of the γ-ray beam with high efficiency and operational simplicity.

Compton scattering↗

Insertion products in the reaction of carbonyl oxide Criegee intermediates with acids: Chloro(hydroperoxy)methane formation from reaction of CH 2 OO with HCl and DCl

The reactions of carbonyl oxides (“Criegee intermediates”) with organic and inorganic acids have been postulated to proceed predominantly by an insertion mechanism, and the hydroperoxy-substituted ester products from reactions with organic acids have been spectroscopically observed. Here this work characterizes the products from one of the simplest reactions of carbonyl oxides with inorganic acids, that of formaldehyde oxide, CH 2 OO, with hydrogen chloride, which occurs via a 1,2-insertion in the H–Cl bond. Reactions of both HCl and DCl isotopologues yield product signal at the mass of the insertion product chloro(hydroperoxy)methane and a dissociative ionisation peak at the mass of the protonated (or deuteronated) Criegee intermediate. The isotopic composition of the insertion product has been measured for reaction mixtures where both HCl isotopologues are present, and the H/D ratio of the product is consistently higher (by a factor of 1.6 ± 0.3) than that of the reactants. This isotope selectivity in the products has smaller uncertainty than the ratio of measured rate coefficients and suggests a normal (k H > k D ) kinetic isotope effect in the reaction. Theoretical kinetics calculations predict a small normal kinetic isotope effect for the overall reaction (k H / k D = 1.35 at 20 Torr N 2 and k H / k D = 1.2 for 1 atm N 2 ) but predict a substantial inverse kinetic isotope effect (k D > k H ) for the stabilisation fraction, in disagreement with the experimental observation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Observation of d(t,n) α Neutrons Following d(d,p)t Reactions in a Deuterium Gas Cell: An Attempt to Repeat Ruhlig’s 1938 Observation of Secondary Reactions

In order to benchmark methods used to calculate reaction-in-flight fusion reactions in inertial confinement fusion and address issues related to the first claimed observation of d(t,n)α reactions in 1938, secondary d(t,n)α reactions have been observed following d(d,p)t reactions in deuterium gas. A pulsed 200-nA, 2.2-MeV deuterium beam from the Triangle Universities Nuclear Laboratory FN tandem accelerator was injected into a cylindrical multiatmosphere deuterium gas target. The incident beam traversed along the target cylinder’s 3-cm symmetry axis after its passage through a Havar entrance foil. Two different Havar foil thicknesses were used to obtain 1.5- and 0.6-MeV deuteron beams entering the deuterium cell. The cylinder’s radius was 2 cm to allow for d(d,p)t tritons emitted perpendicular to the beam to range out in the deuterium gas. The neutron emission from the cell was observed via its time of flight to a liquid scintillator placed at various angles to the beam direction, at a distance of 243 cm. Pulse-shape-discrimination techniques were used to separate neutron and gamma-ray signals seen in the liquid scintillator. The observed probability of ~2 × 10 –4 for inducing secondary d(t,n)α fusion in the gas cell per d(d,p)t reaction is consistent with theoretical expectations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact restoration of Galilei invariance in density functional calculations with quantum Monte Carlo

We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gluon EMC effects in nuclear matter

In this work, we investigate the gluonic structure of nuclei within a mean-field model of nuclear structure based upon the modification of the structure of a bound nucleon, with the nucleon described by the Nambu–Jona-Lasinio model. This approach has been shown to reproduce the European Muon Collaboration (EMC) effect, involving the ratio of the spin-independent structure functions of a heavier nucleus to that of the deuteron. It also predicts a significant nuclear modification for the spin structure functions, known as the polarized EMC effect. Here we report sizeable nuclear modifications of the gluon distributions (a ‘gluon EMC effect’) for the ratios of both the unpolarized and polarized gluon distributions in nuclear matter to those of a free nucleon.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding PDF uncertainty in W boson mass measurements*

Abstract We study the dependence of the transverse mass distribution of charged leptons and the missing energy on parton distributions (PDFs) adapted to W boson mass measurements at the CDF and ATLAS experiments. We compare the shape variations of the distribution induced by different PDFs and find that the spread of predictions from different PDF sets can be significantly larger than the PDF uncertainty predicted by a specific PDF set. We suggest analyzing the experimental data using up-to-date PDFs to gain a better understanding of the PDF uncertainties in W boson mass measurements. We also perform a series of Lagrange multiplier scans to identify the constraints on the transverse mass distribution imposed by individual data sets in the CT18 global analysis. In the case of the CDF measurement, the distribution is mostly sensitive to d -quark PDFs in the intermediate x region, which are largely constrained by DIS and Drell-Yan data on deuteron targets and Tevatron lepton charge asymmetry data.

79 ASTRONOMY AND ASTROPHYSICS↗

Further investigation on the fusion of 6 Li with 209 Bi target at near-barrier energies

The complete and incomplete fusion cross sections for 6 Li+ 209 Bi were measured using the in-beam γ-ray method around the Coulomb barrier. The cross sections of (deuteron captured) incomplete fusion (ICF) products were re-quantified experimentally for this reaction system. Here, the results reveal that the ICF cross section is equivalent to that of complete fusion (CF) above the Coulomb barrier and dominant near or below the Coulomb barrier. A theoretical calculation based on the continuum discretized coupled channel (CDCC) method was performed for the aforementioned CF and ICF cross sections; the result is consistent with the experiments. The universal fusion function (UFF) was also compared with the measured CF cross section for different barrier parameters, demonstrating that the CF suppression factor is significantly influenced by the choice of potential, which can reflect both dynamic and static effects of breakup on the fusion process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simultaneous measurements of unstable and stable Alfvén eigenmodes in JET

Here, we report the novel experimental observation of both unstable and stable toroidicity-induced Alfvén eigenmodes (TAEs) measured simultaneously in a JET tokamak plasma. The three-ion-heating scheme (D-DNBI-3He) is employed to accelerate deuterons to MeV energies, thereby destabilizing TAEs with toroidal mode numbers n = 3–5, each decreasing in mode amplitude. At the same time, the Alfvén eigenmode active diagnostic resonantly excites a stable n = 6 TAE with total normalized damping rate –γ/ω 0 ≈ 1%–4%. Hybrid kinetic-MHD modeling with codes NOVA-K and MEGA both find eigenmodes with similar frequencies, mode structures, and radial locations as in experiment. NOVA-K demonstrates good agreement with the n = 3, 4, and 6 TAEs, matching the damping rate of the n = 6 mode within uncertainties and identifying radiative damping as the dominant contribution. Improved agreement is found with MEGA for all modes: the unstable n = 3–5 and stable n = 2, 6 modes, with the latter two stabilized by higher intrinsic damping and lower fast ion drive, respectively. While some discrepancies remain to be resolved, this unique validation effort gives us confidence in TAE stability predictions for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗

Design of the Multi-Ion Injector Linac for the JLAB EIC (JLEIC)

An Electron Ion Collider (EIC) is the highest priority for future U.S. accelerator-based nuclear physics facility following the completion of the Facility for Rare Isotope Beams (FRIB). Two laboratories are competing to host the future EIC: Brookhaven National Lab. (BNL) and Jefferson Lab. (JLab). The baseline design of JLab's Electron Ion Collider (JLEIC) ion complex comprises a pulsed superconducting (SC) linac injector capable of accelerating all ions from protons to lead, where proton and light ion beams can be polarized. After reviewing the design requirements for the injector linac, important design choices such as the room-temperature (RT) section design, the transition energy between the RT and SC sections and the stripping energy for heavy ions will be discussed. The design of the different linac sections will be presented as well as the results of end-to-end beam dynamics simulations for polarized deuterons and unpolarized lead ions.

Mustapha, B.↗

Di-nucleons do not form bound states at heavy pion mass

We perform a high-statistics lattice QCD calculation of the low-energy two-nucleon scattering amplitudes. To address discrepancies in the literature, the calculation is performed at a heavy pion mass in the limit that the light quark masses are equal to the physical strange quark mass, 𝑚 𝜋 = 𝑚 𝐾 ≃ 714 MeV. Using a state-of-the-art momentum space method, we rule out the presence of a bound di-nucleon in both the isospin 0 (deuteron) and 1 (di-neutron) channels, in contrast with many previous results that made use of compact hexaquark creation operators. To diagnose the discrepancy, we add such hexaquark interpolating operators to our basis and find that they do not affect the determination of the two-nucleon finite-volume spectrum, and thus they do not couple to deeply bound di-nucleons that are missed by the momentum-space operators. Furthermore, we perform a high-statistics calculation of the HAL QCD potential on the same gauge ensembles and find qualitative agreement with our main results. We conclude that di-nucleons do not form bound states at heavy pion masses and that previous identification of deeply bound di-nucleons must have arisen from a misidentification of the spectrum from off-diagonal elements of a correlation function.

Physics - Physics of elementary particles and fiel↗

Toward scalable quantum computations of atomic nuclei

We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3 He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2 H and 3 He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Furthermore, based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Short-depth circuits for efficient expectation-value estimation

The evaluation of expectation values Tr[$\rho \mathcal{O}$] for some pure state $\rho$ and Hermitian operator $\mathcal{O}$ is of central importance in a variety of quantum algorithms. Near-optimal techniques have been developed in the past and require a number of measurements N approaching the Heisenberg limit N = O(1/ ε) as a function of target accuracy ε. The use of quantum phase estimation (QPE) requires, however, long circuit depths C = O(1/ ε) making its implementation difficult on near-term noisy devices. The more direct strategy of operator averaging is usually preferred as it can be performed using N = O(1/ε 2 ) measurements and no additional gates aside from those needed for the state preparation. In this work we use a simple but realistic model to describe the bound state of a neutron and a proton (the deuteron) to show that the latter strategy can require an overly large number of measurements in order to achieve a prefixed relative target accuracy ε r . Further, we propose to overcome this problem using a single step of QPE and classical postprocessing. This approach leads to a circuit depth C = O (ε μ ) (with μ ≥ 0) and to a number of measurements N = O (1/ε 2 + ν ) for 0 < ν ≤ 1 and a much smaller prefactor. We provide detailed descriptions of two implementations of our strategy for ν = 1 and ν ≈ 0.5 and derive appropriate conditions that a particular problem instance has to satisfy in order for our method to provide an advantage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Minimally entangled state preparation of localized wave functions on quantum computers

Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires $O(2^{nQ})$ entangling gates on a quantum computer with $n_Q$ qubits per site. It is conceivable that instead initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a noninteracting scalar field theory, initializing to a symmetric exponential wave function requires $n_Q – 1$ entangling gates, the minimal number necessary to create an entangled state of all $n_Q$ qubits. This is compared with the $2^{nQ – 1} + n_Q – 3 + δ_{nQ, 1}$ required for a symmetric Gaussian wave function. In this work, we explore the initialization of one-site $(n_Q = 4)$, two-site $(n_Q = 3)$, and three-site ( n Q = 3 ) noninteracting scalar field theories with symmetric exponential wave functions using IBM's quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations attainable with $n_Q = 3,4$, these tensor-product wave functions are found to have large overlap with a Gaussian wave function and provide a suitable low-noise initialization for improvement and Somma inflation. Further, in performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum $(|\textbf{0}\gg^{\bigotimes nQ})$ and in medium ($n_Q –1$ qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wave function, two- and three-dimensional Cartesian-space wave functions, and nonrelativistic multinucleon systems built on a localized eigenbasis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum simulation of nuclear inelastic scattering

We present a time-dependent quantum algorithm for nuclear inelastic scattering in the time-dependent basis function on qubits approach. This algorithm aims to quantum simulate a subset of the nuclear inelastic scattering problems that are of physical interest in which the internal degrees of freedom of the reaction system are excited by time-dependent external interactions. We expect that our algorithm will enable an exponential speedup in simulating the dynamics of the subset of the inelastic scattering problems, which would also be advantageous for the applications to more complicated scattering problems. For a demonstration problem, we solve for the Coulomb excitation of the deuteron, where the quantum simulations are performed with IBM Qiskit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

12 C(p,p') 12 C Reaction ( E p =19.5–30 MeV) for Active Interrogation of Special Nuclear Material

Passive detection of special nuclear material (SNM) is challenging due to its inherently low rate of spontaneous emission of penetrating radiation, the relative ease of shielding, and the fluctuating and frequently overwhelming background. Active interrogation, the use of external radiation to increase the emission rate of characteristic radiation from SNM, has long been considered to be a promising method to overcome those challenges. Current active-interrogation systems that incorporate radiography tend to use bremsstrahlung beams, which can deliver high radiation doses. Low-energy ion-driven nuclear reactions that produce multiple monoenergetic photons may be used as an alternative. The 12 C ( p , p ' ) 12 C reaction is one such reaction that could produce large yields of highly penetrating 4.4- and 15.1-MeV gamma rays. This reaction does not directly produce neutrons below the approximately 19.7 MeV threshold, and the 15.1-MeV gamma-ray line is well matched to the photofission cross section of 235 U and 238 U . In this article, we report the measurements of thick-target gamma-ray yields at 4.4 and 15.1 MeV from the 12 C ( p , p ' ) 12 C reaction at proton energies of 19.5, 25, and 30 MeV. Measurements are made with two 3 -in. EJ-309 cylindrical liquid scintillation detectors and thermoluminescent dosimeters placed at 0 ° and 90 ° , with an additional 1.5 -in. Na I ( Tl ) cylindrical scintillation detector at 0 ° . We estimate the highest yields of the 4.4- and 15.1-MeV gamma rays of 1.65 × 10 10 and 4.47 × 10 8 sr - 1 μ C - 1 at a proton energy of 30 MeV, respectively. The yields in all experimental configurations are greater than in a comparable deuteron-driven reaction that produces the same gamma-ray energies— 11 B ( d , n γ ) 12 C . However, a significant increase of the neutron radiation dose accompanies the proton energy increase from 19.5 to 30 MeV.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗