Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Study of dijet events with large rapidity separation in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV

The cross sections for inclusive and Mueller-Navelet dijet production are measured as a function of the rapidity separation between the jets in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV for jets with transverse momentum p T > 35 GeV and rapidity |y| < 4.7. Various dijet production cross section ratios are also measured. A veto on additional jets with p T > 20 GeV is introduced to improve the sensitivity to the effects of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution. The measurement is compared with the predictions of various Monte Carlo models based on leading-order and next-to-leading-order calculations including the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi leading-logarithm (LL) parton shower as well as the LL BFKL resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Characterization of Two Extraordinary AmLi Neutron Sources

Two legacy AmLi sources, MRC-101 and MRC-103, are routinely used for neutron detector characterization and testing campaigns that Pacific Northwest National Laboratory manages for the Defense Threat Reduction Agency within the U.S. Department of Defense. However, applicability of source measurement data has been substantially limited due to the lack of original manufacturer documentation and cross comparison to other sources. This paper summarizes results of the measurement campaign conducted to better quantify source radionuclide content, neutron energy spectrum and emission rate of the AmLi sources. The source characterization study consisted of three major consecutive measurement exercises. First, gamma spectrometry measurements using a high purity germanium detector were used to estimate the radionuclide content and potential impurities of the AmLi sources. Then neutron spectral measurements were taken to determine the unique neutron energy distribution. Finally, assessments of the total neutron yields of the sources were made with measurements in a neutron well counter. The neutron detection efficiency for these measurements was determined via detailed Monte Carlo modeling based on the neutron energy distribution measured with the rotating neutron spectrometer. The most exciting and unexpected finding of the study was a significant contamination of the source material with beryllium, contributing up to 17% to the total neutron emission rate and correspondingly changing the neutron energy spectrum. Also, the obtained lithium-only portion of the measured energy distribution of the sources was compared against several AmLi spectra previously reported by different researchers.

AWCC, AmLi, beryllium contamination, neutron spect↗
DOE OSTI2024

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Unraveling the Early-Stage Ordering of Krypton Solid Bubbles in Molybdenum: A Multimodal Study

Self-organization of defects such as fission gas bubbles in materials can lead to high inventory capacity for fission gas storage and help mitigate swelling caused by fission gases in nuclear fuel materials under radiation in nuclear reactors. Here, we report the physical mechanism of self-organization of krypton (Kr) gas bubbles in molybdenum (Mo) under ion implantation. The ion fluence and temperature-dependent formation of Kr solid bubble superlattice (SBS) in Mo were investigated by using both synchrotron-based small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM). Early stage self-organization of gas bubbles is observed at a fluence of 2.5 × 10 16 ions/cm 2 at temperatures of 300–400 °C. The bubble lattice constant increases with increasing implantation temperature from 300 to 400 °C. Both experiments and atomic kinetic Monte Carlo modeling indicate that the Kr solid bubbles are weakly ordered in comparison to previously studied helium (He) gas bubble superlattice (GBS) while the lattice constant are relatively smaller for Kr SBS compared to that of He GBS. The irradiation conditions suggest that spinodal decomposition, which is a form of phase separation, probably precedes gas bubble ordering in Mo. Altogether, our work sheds light on the formation mechanism of noble gas superlattice toward the development of radiation-tolerant materials which are important for the design of advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Dissociation Pathways of BCl 3 on the Si(100)-2 × 1 Surface

BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

Campbell, Quinn T. [Sandia National Laboratories (↗
DOE OSTI2023

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Overcoming the Entropy Penalty of Direct Air Capture for Efficient Gigatonne Removal of Carbon Dioxide

Atmospheric carbon poses an existential threat to civilization via global climate change. Hundreds of gigatonnes of carbon dioxide must be removed from earth’s atmosphere in the next three decades, necessitating a low-cost, energy-efficient process to extract low concentrations of carbon dioxide for conversion to a stable material permanently stored for thousands of years. In this work, the challenge of removing gigatonnes of CO 2 is described via the scale of effort and the thermodynamics of collecting and reducing this diffuse chemical, the accumulation of which imparts a substantial entropy penalty on any atmospheric carbon capture process. The methods of CO 2 reduction combined with upstream direct air capture (DAC) including absorption, membrane separation, and adsorption are compared with biomass torrefaction and permanent burial (BTB). A Monte Carlo model assesses the mass, energy, and economics of the full process of biomass torrefaction from biomass collection and transport to stable carbon burial to determine that 95% of scenarios could remove carbon for less than $200 per CO 2 -tonne-equivalent. Torrefied carbon is further discussed for its long-term stability and availability at the scale required to substantially mitigate the threat of climate change.

biomass↗
DOE OSTI2024

Amorphous Oxyhalide Matters for Achieving Lithium Superionic Conduction

The recently surged halide-based solid electrolytes (SEs) are great candidates for high-performance all-solid-state batteries (ASSBs), due to their decent ionic conductivity, wide electrochemical stability window, and good compatibility with high-voltage oxide cathodes. In contrast to the crystalline phases in halide SEs, amorphous components are rarely understood but play an important role in Li-ion conduction. Here, we reveal that the presence of amorphous component is common in halide-based SEs that are prepared via mechanochemical method. The fast Li-ion migration is found to be associated with the local chemistry of the amorphous proportion. Taking Zr-based halide SEs as an example, the amorphization process can be regulated by incorporating O, resulting in the formation of corner-sharing Zr–O/Cl polyhedrons. This structural configuration has been confirmed through X-ray absorption spectroscopy, pair distribution function analyses, and Reverse Monte Carlo modeling. The unique structure significantly reduces the energy barriers for Li-ion transport. As a result, an enhanced ionic conductivity of (1.35 ± 0.07) × 10 –3 S cm –1 at 25 °C can be achieved for amorphous Li 3 ZrCl 4 O 1.5 . In addition to the improved ionic conductivity, amorphization of Zr-based halide SEs via incorporation of O leads to good mechanical deformability and promising electrochemical performance. In conclusion, these findings provide deep insights into the rational design of desirable halide SEs for high-performance ASSBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Chemical Framework to Design Linear-like Relaxors toward Capacitive Energy Storage

ABO 3 -type perovskite relaxor ferroelectrics (RFEs) have emerged as the preferred option for dielectric capacitive energy storage. However, the compositional design of RFEs with high energy density and efficiency poses significant challenges owing to the vast compositional space and the absence of general rules. Here, we present an atomic-level chemical framework that captures inherent characteristics in terms of radius and ferroelectric activity of ions. By categorizing A/B-site ions as host framework, rattling, ferroelectrically active, and blocking ions and assembling these four types of ions with specific criteria, linear-like relaxors with weak locally correlated and highly extendable unit-cell polarization vectors can be constructed. As example, we demonstrate two new compositions of Bi 0.5 K 0.5 TiO 3- based and BaTiO 3- based relaxors, showing extremely high recoverable energy densities of 17.3 and 12.1 J cm –3 , respectively, both with a high efficiency of about 90%. Further, the role of different types of ions in forming heterogeneous polar structures is identified through element-specific local structure analysis using neutron total scattering combined with reverse Monte Carlo modeling. Further, our work not only opens up new avenues toward rational compositional design of high energy storage performance lead-free RFEs but also sheds light on atomic-level manipulation of functional properties in compositionally complex ferroelectrics.

25 ENERGY STORAGE↗
DOE OSTI2023

Local inversion-symmetry breaking in a bismuthate high-Tc superconductor

Abstract The doped perovskite BaBiO 3 exhibits a maximum superconducting transition temperature ( T c ) of 34 K and was the first high- T c oxide to be discovered, yet pivotal questions regarding the nature of both the metallic and superconducting states remain unresolved. Although it is generally thought that superconductivity in the bismuthates is of the conventional s -wave type, the pairing mechanism is still debated, with strong electron-phonon coupling and bismuth valence or bond disproportionation possibly playing a role. Here we use diffuse x-ray scattering and Monte Carlo modeling to study the local structure of Ba 1-x K x BiO 3 across its insulator-metal boundary. We find no evidence for either long- or short-range disproportionation, which resolves a major conundrum, as disproportionation and the related polaronic effects are likely not relevant for the metallic and superconducting states. Instead, we uncover nanoscale structural correlations that break inversion symmetry, with far-reaching implications for the electronic physics. This unexpected finding furthermore establishes that the bismuthates belong to the broader classes of materials with hidden spin-orbit coupling and a tendency towards inversion-breaking displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neutron spectroscopy of plutonium using a handheld detection system

The ability to distinguish multiple forms of plutonium from one another, such as oxide and metal, is paramount in areas of nuclear nonproliferation and international safeguards. In its metal form, plutonium can be readily used in a nuclear weapon, while oxide forms are associated with nuclear reactor fuel. Oxide-based plutonium forms emit neutrons with an energy spectrum that is significantly different from the fission neutrons that are emitted from plutonium metal. Organic scintillation detectors output pulses that are proportional to the neutron energy deposited, and therefore present a means of distinguishing these plutonium forms based on their energy spectra. In this work, metal and oxide forms of plutonium were measured using a handheld detection system based on an organic glass scintillator. Monte Carlo modeling of these experiments was performed to provide insight into the origin of the features in the observed light output spectra. Through analysis of multiple regions of these spectra, in a matter of minutes we were able to unambiguously discriminate oxide and metal plutonium forms from one another and from a plutonium-beryllium neutron source, which was considered for comparison because these sources are commonly used in industrial applications. The ability to discriminate weapons-usable material from nuclear reactor fuel has applications in nuclear treaty verification and safeguards.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Structural regulation-induced Li-electron disentanglement for stabilized oxygen redox of Li-excess disordered rock-salt cathode materials

Since the discovery of its electrochemical activity, Li-excess disordered rock-salt (DRX) cathode material has received worldwide attention as it sets up a new way to exploit oxygen redox beyond the conventional layered structure with late-3d transition metals. However, the intricate structure-function relationship in the disordered lattice of the DRX material fogs the researcher's lens on the underlying redox mechanisms. Here, in this study, we employ a synergistic approach combining neutron total scattering with reverse Monte Carlo modeling and density functional theory calculations to unravel the landscape of oxygen redox reactions in DRX. Redox activities are evaluated in diverse oxygen clusters (OLi x TM 6-x ) and the spatial distribution of these clusters in the model DRX structure (Li 1.16 Ti 0.37 Ni 0.37 Nb 0.1 O 2 and Li 1.2 Ti 0.35 Ni 0.35 Nb 0.1 O 1.8 F 0.2 ) is explicitly determined. The results unveil that by regulating the short-range ordering between cations, fluorine atoms can effectively decouple the location of Li extraction and electron depletion. Such disentanglement between the Li reservoir and electron reservoir in the DRX lattice could play a pivotal role in protecting the oxidized oxygen and preserving the lattice framework during cycling. Through a tentatively designed non-fluorinated DRX oxide realizing similar Li-electron decoupling, an obvious enhancement of the cycling capability can be achieved without compromising the capacity release.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Combinatorial deposition of Au–Bi alloys via high-rate magnetron sputtering

Gold–bismuth (Au–Bi) alloy films are promising candidate materials for inertial confinement fusion (ICF) hohlraums due to their high laser-to-x-ray conversion efficiency, particularly compared with Au, Ta–Au, and Bi hohlraums. However, the fabrication of uniform and dense Au–Bi alloy films remains a challenge. Here, we use a combination of Monte-Carlo modeling and experiments to demonstrate that the microstructure and properties of Au–Bi alloy films can be greatly improved when direct-current magnetron sputtering with high deposition rates of ⩾5 μm h -1 is used. Resultant films are ~90% of their maximum theoretical densities and have low O content of <1 at.%. Films with Bi content above 40 at. % exhibit high electrical resistivity >100μΩ cm, making them suitable for both magnetized and non-magnetized ICF schemes.

Materials science↗
DOE OSTI2024

Correlations of baryon and charge stopping in heavy ion collisions*

Baryon numbers are theorized to be carried by valence quarks in the standard QCD picture of the baryon structure. Another theory proposed an alternative baryon number carrier, a non-perturbative Y-shaped configuration of the gluon field, called the baryon junction in the 1970s. However, neither of these theories has been verified experimentally. Recently, searching for the baryon junction by investigating the correlation of net-charge and net-baryon yields at midrapidity in heavy-ion collisions has been suggested. Here, this paper presents studies of such correlations in collisions of various heavy ions from oxygen to uranium with the UrQMD Monte Carlo model. The UrQMD model implements valence quark transport as the primary means of charge and baryon stopping at midrapidity. Detailed studies are also conducted for isobaric $_{40}^{96}{\rm{Zr}}$ + $_{40}^{96}{\rm{Zr}}$ and $_{44}^{96}{\rm{Ru}}$ + $_{44}^{96}{\rm{Ru}}$ collisions. We found a universal trend of charge stopping with respect to baryon stopping and discovered that the charge stopping is always greater than the baryon stopping. This study provides a model baseline in valence quark transport for what is expected in net-charge and net-baryon yields at the midrapidity of relativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Beam modulation and bump-on-tail effects on Alfvén eigenmode stability in DIII-D

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Furthermore, detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic 'kicks' by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗