Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Yield”

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 1,027 records · Page 57

Reconstruction of neutrino events in the Accelerator Neutrino Neutron Interaction Experiment. Part I

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) was designed to reconstruct neutrino events from the Fermilab Booster Neutrino Beam (BNB) with the parallel goals of measuring neutron production in interactions with oxygen and serving as a testbed for new technology. The ANNIE detector consists of a 26-ton water Cherenkov target tank instrumented with conventional photomultiplier tubes (PMTs), a downstream tracking muon spectrometer, and an upstream double wall of plastic scintillator to serve to veto charged particles incoming from neutrino events that occur upstream of the experimental setup. ANNIE has also deployed multiple Large-Area Picosecond PhotoDetectors (LAPPDs) and a test vessel of water-based liquid scintillator (WbLS). This paper describes the event reconstruction performance of the detector before implementation of these novel technologies, which will serve as a baseline against which their impact can be measured. That said, even the techniques used for event reconstruction using only the conventional PMT array and muon spectrometer are significantly different than those used in other water Cherenkov detectors due to the small size of ANNIE (which makes nanosecond-scale timing not as useful as in a large detector) and the availability of reconstruction information from the tracking muon spectrometer. We demonstrate that combining the information from these two elements into a single fit using only pattern recognition yields a muon vertex uncertainty of 60 cm, a directional uncertainty of 13.2 degrees, and energy reconstruction uncertainty of about 10% for BNB muon neutrino Charged Current Zero Pion (CC0π) events.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of Microstructural Composition, Porosity, and Microcracks on the Elastic Moduli of Nuclear Graphites

The elastic responses of nuclear graphite not only depend on the graphitic content itself but are largely dictated by the microstructural constitution of the material. The types of raw materials combined with the manufacturing processes used to produce the graphite yield the microstructural content that primarily includes graphite filler, graphitized pitch binder, and voids/defects that typically occupy approximately 20% of the volume. Among these microstructural components, porosity and microcracking (considered to be part of voids/defects) heavily influence the overall properties of the material, including the elastic moduli. It is widely accepted that the primary effect of oxidation is to increase porosity, but the related effect on the moduli cannot be explained satisfactorily by simply noting changes to porosity. Models describing the elastic moduli of porous, polycrystalline graphite materials have been developed to interpret experimental determinations of Young's modulus and shear modulus in oxidized graphites. Beyond porosity, the moduli are heavily influenced by microcracks, and their effects can be assessed using physical property models based on orientation distribution coefficients and microcrack-modified, single-crystal moduli to represent the combined effects of crystallite orientation and microcracking. This work demonstrates how most modulus measurements in nuclear graphites can be understood using relatively simple models that describe the effects of porosity and microcracking. We also present directions to be pursued to account for microstructure-related effects that occur as a result of neutron irradiation.

Spicer, James B.↗

Integral Measurement of the 237 Np($n$, γ) 238 Np Cross Section with Fast Neutrons

Accurate knowledge of the 237 Np($n$; γ) 238 Np cross section at fast neutron energies is important for applied nuclear science. The presently available experimental data has large disagreements in the fast neutron region. A measurement of the 237 Np($n$; γ) 238 Np cross section was performed using a well characterized fast neutron source and the results was compared with previous measurements and current nuclear data evaluations. This will provide an integral measurement that can be used as a benchmark for current evaluations. Multiple samples of 237 Np were irradiated in the Godiva-IV critical assembly. Following the irradiation, the samples placed in a γ-ray counting setup and the γ-rays emitted from the decay of 238 Np were measured over a time period of approximately 7 days. Multiple γ-ray decay branches of 238 Np were observed. The observed activity of 238 Np was used to calculate the amount of 238 Np produced during the irradiation via the 237 Np($n$; γ) 238 Np and an integral cross section of 359(10) mb was measured. The 238 Np half-life has been measured with a result of 50.34(7) hours. An integral measurement of the 237 Np($n$; γ) 238 Np cross section has been measured. This result has been compared to the calculated integral cross section using the ENDF/BVIII.0 and JENDL-4.0 evaluated cross sections and found to be in good agreement with ENDF/BVIII.0 value of 367 mb. The JENDL-4.0 evaluation yields a value of 312 mb, which is significantly lower than the measured cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Large Destabilization of (TiVNb)-Based Hydrides via (Al, Mo) Addition: Insights from Experiments and Data-Driven Models

High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH 2 ), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. Here, in this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al 0.05 (TiVNb) 0.95–x Mo x (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.

36 MATERIALS SCIENCE↗

Verification and Quantification of Single Event Effects on High Speed SRAM in Terrestrial Environments

As integrated circuits become more sensitive to charged particles and neutrons, anomalous performance due to single event effects (SEE) is a concern and requires experimental verification and quantification. The Center for Applied Radiation Research (CARR) at Prairie View A&M University has developed experiments as a participant in the NASA ER-2 Flight Program, the APEX balloon flight program and the Student Launch Program. Other high altitude and ground level experiments of interest to DoD and commercial applications are being developed. The experiment characterizes the SEE behavior of high speed and high density SRAM's. The system includes a PC-104 computer unit, an optical drive for storage, a test board with the components under test, and a latchup detection and reset unit. The test program will continuously monitor the stored checkerboard data pattern in the SW and record errors. Since both the computer and the optical drive contain integrated circuits, they are also vulnerable to radiation effects. A latchup detection unit with discrete components will monitor the test program and reset the system when necessary. The first results will be obtained from the NASA ER-2 flights, which are now planned to take place in early 1998 from Dryden Research Center in California. The series of flights, at altitudes up to 70,000 feet, and a variety of flight profiles should yield a distribution of conditions for correlating SEES. SEE measurements will be performed from the time of aircraft power-up on the ground throughout the flight regime until systems power-off after landing.

Huff, H.↗

Complex structure of the supernova remnant HB 3

HB 3 is an old, large (84 pc diameter) supernova remnant associated with the W3 H II region/molecular cloud complex. Observations of the imaging proportional counter (IPC) onboard the Einstein X-ray astronomy satellite have been reprocessed to yield a contour map of X-ray brightness and spectra of various regions of this remnant. The measured IPC flux is 2.4 x 10 to the -11th ergs per sq cm per s, giving a 0.2-4 keV luminosity of 1.6 x 10 to the 35th ergs/s for a column densityof 6 x 10 to the 21st per sq cm. The measured X-ray temperatures reveal a decrease from center to limb of the remnant of 1-0.3 keV. HB 3 is in the late adiabatic blast-wave phase of evolution, 30,000 to 50,000 yr old and with an initial blast energy of 3 x 10 to the 50th ergs. The X-ray map is compared with available radio and optical images. In X-rays, HB 3 has two components - a diffuse emission inside the 84 pc radio remnant and a ring of emission at the center of 30 pc in diameter. The diffuse emission is similar to that from other supernova remnants which are moderately obscured (column density, nH approximately 10 to the 22nd per sq cm). Three possibilities for the origin of the ring are explored: (1) a second supernova remnant, (2) a shocked shell in the interstellar medium surrounding HB 3, and (3) reverse-shock heated ejecta. There is no hot neutron star within the remnant.

Leahy, D. A.↗

Photonuclear Physics in SCALE

Radiation shielding modeling involves a wide range of applications, such as power reactors, critical experiments, portable sources, linear accelerators, and fusion systems. Particle energy ranges vary from thermal energies to fast energy spectra. MAVRIC, as the radiation shielding sequence of the SCALE code system [1] developed and maintained by Oak Ridge National Laboratory, is capable of modeling both neutron and photon interactions in matter and is widely used by the US Nuclear Regulatory Commission, Department of Energy, the National Nuclear Security Administration, and radiation shielding practitioners all around the world. Although most applications don’t involve highly energetic photons that can cause photonuclear reactions and create coupled chains of high energy, ionized electrons and bremsstrahlung radiation within a material medium, when present, these reactions can have a significant impact upon the downstream particle fluxes from incident photons above 20 MeV. For example, linear accelerators used for medical isotope production are examples of applications in which high-energetic photons are generated and used for driving or controlling other components of a complex facility. Radiation shielding analyses of such a facility would require accounting for any and all significant sources of radiation. Unfortunately, MAVRIC included no data and methods to model either photonuclear reactions or the bremsstrahlung photons until this work, which enabled accounting for both. Photonuclear reaction cross sections were imported from the Evaluated Nuclear Data File libraries, and corresponding methods were implemented in the SCALE code system; see Section 2 for how the nuclear data are imported and used. Bremsstrahlung photons were modeled by assuming thick-target approximation whereby all generated electrons and positrons due to photon interactions are assumed to be absorbed within the same material in which they were born, yielding a quick particle transport simulation method by using fast lookup tables for generating subsequent gammas; see Section 4 for details. A suite of verification models was used to test implemented methods and associated data. Tests yielded satisfactory results for the new features. Enhancements will be released as part of the next major SCALE release (v7.0), which is anticipated to be deployed in 2025, with beta releases available in the interim.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Application of a Hough Search for Continuous Gravitational Waves on Data from the Fifth LIGO Science Run

We report on an all-sky search for periodic gravitational waves in the frequency range 50-1000 Hertz with the first derivative of frequency in the range −8.9 × 10(exp −10) Hertz per second to zero in two years of data collected during LIGO's fifth science run. Our results employ a Hough transform technique, introducing a chi(sup 2) test and analysis of coincidences between the signal levels in years 1 and 2 of observations that offers a significant improvement in the product of strain sensitivity with compute cycles per data sample compared to previously published searches. Since our search yields no surviving candidates, we present results taking the form of frequency dependent, 95% confidence upper limits on the strain amplitude h(sub 0). The most stringent upper limit from year 1 is 1.0 × 10(exp −24) in the 158.00-158.25 Hertz band. In year 2, the most stringent upper limit is 8.9 × 10(exp −25) in the 146.50-146.75 Hertz band. This improved detection pipeline, which is computationally efficient by at least two orders of magnitude better than our flagship Einstein@Home search, will be important for 'quicklook' searches in the Advanced LIGO and Virgo detector era.

neutron stars↗

Theory of nuclear fission

Atomic nuclei are quantum many-body systems of protons and neutrons held together by strong nuclear forces. Under the proper conditions, nuclei can break into two (sometimes three) fragments which will subsequently decay by emitting particles. This phenomenon is called nuclear fission. Since different fission events may produce different fragmentations, the end-products of all fissions that occurred in a small chemical sample of matter comprise hundreds of different isotopes, including α particles, together with a large number of emitted neutrons, photons, electrons and antineutrinos. The extraordinary complexity of this process, which happens at length scales of the order of a femtometer, mostly takes less than a femtosecond but is not entirely over until all the lingering β decays have completed – which can take years – is a fascinating window into the physics of atomic nuclei. While fission may be more naturally known in the context of its technological applications, it also plays a crucial role in the synthesis of heavy elements in astrophysical environments. In both cases, simulations are needed for the many systems or energies inaccessible to experiments in the laboratory. In this context, the level of accuracy and precision required poses formidable challenges to nuclear theory. Overall, the goal of this article is to provide a comprehensive overview of the theoretical methods employed in the description of nuclear fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗

Alpha-heating analysis of burning plasma and ignition experiments on the National Ignition Facility

A recent experiment conducted on the National Ignition Facility (NIF) described in the study by Abu-Shawareb et al. achieved a fusion yield output of 1.3 MJ from ~ 220 kJ of x-ray energy absorbed by the capsule, demonstrating remarkable progress in the field of laser driven inertial confinement fusion. In the study by A. R. Christopherson [“Effects of charged particle heating on the hydrodynamics of inertially confined plasmas,” Ph.D. thesis (2020)], the plasma conditions needed to claim the onset of ignition and burn propagation were outlined and multiple criterion were provided to assess progress in inertial fusion experiments. In this work, we modify the metrics from A. R. Christopherson to accurately calculate performance metrics for indirect-drive experiments on the NIF. We also show that performance metric trends observed in NIF data are consistent with theory and simulations. This analysis indicates that all the identified criterion for ignition and burn propagation have been exceeded by experiment 210 808.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

“The Goose” Pulsar Wind Nebula of PSR J1016–5857: The Birth of a Plerion

We report the results of X-ray (CXO) and radio (ATCA) observations of the pulsar wind nebula(PWN) powered by the young pulsar PSR J1016–5857, which we dub “the Goose” PWN. In bothbands the images reveal a tail-like PWN morphology which can be attributed to pulsar’s motion. Bycomparing archival and newCXOobservations, we measure the pulsar’s proper motionμ= 28.8±7.3mas yr−1, yielding a projected pulsar velocityv≈440±110 km s−1(atd= 3.2 kpc); its direction isconsistent with the PWN shape. Radio emission from the PWN is polarized, with the magnetic fieldoriented along the pulsar tail. The radio tail connects to a larger radio structure (not seen in X-rays)which we interpret as a relic PWN (also known as a plerion). The spectral analysis of theCXOdatashows that the PWN spectrum softens from Γ = 1.7 to Γ≈2.3−2.5 with increasing distance fromthe pulsar. The softening can be attributed to the rapid synchrotron burn-off, which would explainthe lack of X-ray emission from the older relic PWN. In addition to non-thermal PWN emission, wedetected thermal emission from a hot plasma which we attribute to the host SNR. The radio PWNmorphology and the proper motion of the pulsar suggest that the reverse shock passed through thepulsar’s vicinity and pushed the PWN to one side.

Pulsars↗

A Study of the Non-Thermal X-Ray Emission of Shell-Type Supernova Remnants

The term of the second year of the award is the period from March 15, 2001 to March 14, 2002. As was the specified goal of the second year, we analyzed the spatial and spectral X-ray data for several young supernova remnants. I published a paper about an analysis of the ROSAT, ASCA, and RXTE data for the supernova remnant SN 1006. A copy of this paper is enclosed. As described in the paper, we believe that we accurately modeled the nonthermal X-ray emission from the remnant. The results of this analysis are used to infer properties about the cosmic rays accelerated in the remnant and to argue that the strength of the magnetic field in the remnant is considerably larger than the value of about 10 micro G reported elsewhere. The results were presented at the August 2001 International Cosmic Ray Conference in Hamburg, German),. I began analyzing new Chandra X-ray data for SN 1006. This analysis will yield the first measure of the strength of the magnetic field in the remnant for the first time. Preliminary results support our previous conclusion that the magnetic field strength in the remnant is much larger than 10 micro G. The field strength seems to be about the strength expected based on an equipartition calculation. The result supports recent models that describe the how the shock structure is influenced by the efficient acceleration of cosmic rays. This work will be presented at the April 2002 High Energy Astrophysics Division meeting in Albuquerque and published this summer. A copy of the abstract for the talk is enclosed. I began studying new Chandra X-ray data for the supernova remnant Cas A. The results of this work show that the forward shock is a region where cosmic-ray electrons are accelerated, which is consistent with theoretical expectations. The work was presented at the September 2001 Two Years of Science with Chandra symposium in Washington, DC. A copy of the poster paper is enclosed. Dr. Thomas Pannuti, whose research work is supported by the award, analyzed ROSAT, ASCA, and RXTE data for the supernova remnant G347.3-0.5. The results show for the first time that thermal X-ray emission is produced in the remnant. As expected, the thermal emission is consistent with a model in which the remnant is expanding into a very low density environment. The results also provide an accurate description of the nonthermal emission from the remnant. Dr. Pannuti presented this work at several conferences. A copy of the paper for the proceedings of the August 2001 Neutron Stars in Supernova Remnants symposium in enclosed. The work will be submitted to the Astrophysical Journal in the next few months.

Allen, Glenn E.↗

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Magnetization Process of Atacamite: A Case of Weakly Coupled S = 1/2 Sawtooth Chains

We present a combined experimental and theoretical study of the mineral atacamite Cu 2 Cl(OH) 3 . Density-functional theory yields a Hamiltonian describing anisotropic sawtooth chains with weak 3D connections. Experimentally, we fully characterize the antiferromagnetically ordered state. Here, magnetic order shows a complex evolution with the magnetic field, while, starting at 31.5 T, we observe a plateaulike magnetization at about M sat /2. Based on complementary theoretical approaches, we show that the latter is unrelated to the known magnetization plateau of a sawtooth chain. Instead, we provide evidence that the magnetization process in atacamite is a field-driven canting of a 3D network of weakly coupled sawtooth chains that form giant moments.

1-dimensional spin chains↗