Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct Simulation Monte Carlo”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dosimetric and biological impact of activity extravasation of radiopharmaceuticals in PET imaging

The increasing use of nuclear medicine and PET imaging has intensified scrutiny of radiotracer extravasation. To our knowledge, this topic is understudied but holds great potential for enhancing our understanding of extravasation in clinical PET imaging. This work aims to (1) quantify the absorbed doses from radiotracer extravasation in PET imaging, both locally at the site of extravasation and with the extravasation location as a source of exposure to bodily organs and (2) assess the biological ramifications within the injection site at the cellular level. A radiation dosimetry simulation was performed using a whole-body 4D Extended Cardiac-Torso (XCAT) phantom embedded in the GATE Monte Carlo platform. A 10-mCi dose of 18 F-FDG was chosen to simulate a typical clinical PET scan scenario, with 10% of the activity extravasated in the antecubital fossa of the right arm of the phantom. The extravasation volume was modeled as a 5.5 mL rectangle in the hypodermal layer of skin. Absorbed dose contributions were calculated for the first two half-lives, assuming biological clearance thereafter. Dose calculations were performed as absorbed doses at the organ and skin levels. Energy deposition was simulated both at the local extravasation site and in multiple organs of interest and converted to absorbed doses based on their respective masses. Each simulation was repeated ten times to estimate Monte Carlo uncertainties. Biological impacts on cells within the extravasated volume were evaluated by randomizing cells and exposing them to a uniform radiation source of 18 F and 68 Ga. Particle types, their energies, and direction cosines were recorded in phase space files using a separate Geant4 simulation to characterize their entry into the nucleus of the cellular volume. Subsequently, the phase space files were imported into the TOPAS-nBio simulation to assess the extent of DNA damage, including double-strand breaks (DSBs) and single-strand breaks (SSBs). Organ-level dosimetric estimations are presented for 18 F and 68 Ga radionuclides in various organs of interest. With 10% extravasation, the hypodermal layer of the skin received the highest absorbed dose of 1.32 ± 0.01 Gy for 18 F and 0.99 ± 0.01 Gy for 68 Ga. The epidermal and dermal layers received absorbed doses of 0.07 ± 0.01 Gy and 0.13 ± 0.01 Gy for 18 F, and 0.14 ± 0.01 Gy and 0.29 ± 0.01 Gy for 68 Ga, respectively. In the extravasated volume, 18 F caused an average absorbed dose per nucleus of 0.17 ± 0.01 Gy, estimated to result in 10.58 ± 0.50 DSBs and 268.11 ± 12.43 SSBs per nucleus. For 68 Ga, the absorbed dose per nucleus was 0.11 ± 0.01 Gy, leading to an estimated 6.49 ± 0.34 DSBs and 161.24 ± 8.12 SSBs per nucleus. Absorbed doses in other organs were on the order of micro-gray (µGy). The likelihood of epidermal erythema resulting from extravasation during PET imaging is low, as the simulated absorbed doses to the epidermis remain below the thresholds that trigger such effects. Moreover, the organ-level absorbed doses were found to be clinically insignificant across various simulated organs. The minimal DNA damage at the extravasation site suggests that long-term harm, such as radiation-induced carcinogenesis, is highly unlikely.

DNA strand breaks↗
DOE OSTI2021

Nematic order driven by superconducting correlations

The interplay of nematicity and superconductivity has been observed in a wide variety of quantum materials. To explore this interplay, we consider a two-dimensional (2D) array of nematogens, local droplets with Z 3 nematicity, coupled to a network of Josephson junction wires. Using finite temperature classical Monte Carlo simulations, we elucidate the phase diagram of this model and show that the development of superconducting correlations and the directional delocalization of Cooper pairs can promote nematogen ordering, resulting in long-range nematic order. We obtain the transport properties of our model within an effective resistor network picture. We discuss these ideas in the context of the 2D electron gas at the (111) KTaO 3 interface and the doped topological insulators Nb x Bi 2 Se 3 and Cu x Bi 2 Se 3 . Overall, our work makes contact with Phil Anderson’s numerous contributions to broken symmetries driven by the saving of kinetic energy, including double exchange ferromagnetism and the interlayer tunneling theory of high Tc csuperconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Defect evolution in Ni and solid-solution alloys of NiFe and NiFeCoCr under ion irradiation at 16 and 300 K

Single-phase concentrated solid-solution alloys (SP-CSAs) have shown unique chemical complexity at the levels of electrons and atoms, and their defect evolution is expected to be different from conventional dilute alloys. Single crystals of Ni, NiFe and NiFeCoCr are chosen as model systems to understand the chemical complexity on defect formation and damage accumulation in SP-CSAs under ion irradiation. The high-quality crystals were irradiated at 16 and 300 K to different ion fluences, to form irradiated region with little to heavy damages. Additionally, the ion-induced damage was determined using Rutherford backscattering spectrometry technique along a channeling direction (RBS/C) and the level of lattice damage in irradiated Ni and SP-CSAs was quantified from Monte Carlo (MC) simulations. The results are interpreted using the Multi Step Damage Accumulation model to reveal material damage accumulation kinetics. Key findings of the study are that in case of room temperature irradiations the damage level measured for complex alloys at the highest irradiation fluence of 2 × 10 15 cm -2 (~3 dpa) is significantly higher than that obtained for pure nickel samples and suggest two-step damage accumulation process with a defect transformation taking place at a fluence of about 1.5 × 10 15 cm -2 . Furthermore, structural and damage kinetic differences clearly imply that, with increasing degree of chemical complexity and high solid-solution strengthening effects from Ni to NiFe and to NiFeCoCr, the enhanced lattice stiffness resists to randomization of atomic configurations and inhibits the growth of extended defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

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

Designs and Modeling of Architect Multi-material Scintillator System

We present a new conceptual radiation detector, the Architected Multimaterial Scintillator System, that utilizes a scintillator composed of multiple materials arranged in architected structures to enable new capabilities. By structuring differently-dyed materials, the wavelength of the scintillation light encodes additional information in radiation measurements. These structures can be realized through additive manufacture (3D-printing). Two classes of this concept are described and modelled using Monte Carlo simulations to evaluate their performance. The first detector design uses dye microstructures to encode particle tracking information, allowing for directional neutron detection and gamma/neutron discrimination. The second design uses a dye gradient to indicate the position of radiation along the gradient. The simulation results indicate this new concept in radiation detection can achieve strong performance in a variety of capabilities including particle identification, directionality and spectroscopy measurements, and particle position reconstruction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

On the Relationship Between Shallow Cumulus Cloud Field Properties and Surface Solar Irradiance

Abstract Shallow cumulus clouds exhibit highly three‐dimensional (3‐D) spatial structure leading to complex variability in the surface solar irradiance (SSI) beneath. This variability is captured by the typically bimodal shape of the SSI probability density function (PDF). Using large eddy simulation to generate well‐resolved cloud fields and Monte Carlo 3‐D radiative transfer to reproduce realistic associated SSI PDFs, we seek direct relationships between the cloud field properties and the SSI PDF shape. Applying both random forest and artificial neural network algorithms, we find variations in the two modes of the SSI PDF are well predicted by just a handful of cloud field properties. The two algorithms utilize cloud properties similarly, with indistinguishable performance despite their different architectures. These results offer a marked improvement in realism relative to one‐dimensional radiative transfer while bypassing computationally expensive 3‐D radiative transfer, with immediate application to renewable energy assessments, and potential for several other geophysical applications.

58 GEOSCIENCES↗
DOE OSTI2021

High-resolution planar electron beam induced current in bulk diodes using high-energy electrons

Understanding the impact of high-energy electron radiation on device characteristics remains critical for the expanding use of semiconductor electronics in space-borne applications and other radiation harsh environments. Here in this paper, we report on in situ measurements of high-energy electron radiation effects on the hole diffusion length in low threading dislocation density homoepitaxial bulk n-GaN Schottky diodes using electron beam induced current (EBIC) in high-voltage scanning electron microscopy mode. Despite the large interaction volume in this system, quantitative EBIC imaging is possible due to the sustained collimation of the incident electron beam. This approach enables direct measurement of electron radiation effects without having to thin the specimen. Using a combination of experimental EBIC measurements and Monte Carlo simulations of electron trajectories, we determine a hole diffusion length of 264 ± 11 nm for n-GaN. Irradiation with 200 kV electron beam with an accumulated dose of 24 × 10 16 electrons cm -2 led to an approximate 35% decrease in the minority carrier diffusion length.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Role of hot electrons in shock ignition constrained by experiment at the National Ignition Facility

Shock ignition is a scheme for direct drive inertial confinement fusion that offers the potential for high gain with the current generation of laser facility; however, the benefits are thought to be dependent on the use of low adiabat implosions without laser–plasma instabilities reducing drive and generating hot electrons. A National Ignition Facility direct drive solid target experiment was used to calibrate a 3D Monte Carlo hot-electron model for 2D radiation-hydrodynamic simulations of a shock ignition implosion. The α = 2.5 adiabat implosion was calculated to suffer a 35% peak areal density decrease when the hot electron population with temperature T ℎ = 55 keV and energy E ℎ = 13 kJ was added to the simulation. Optimizing the pulse shape can recover ~1/3 of the peak areal density lost due to a change in shock timing. Here, despite the harmful impact of laser–plasma instabilities, the simulations indicate shock ignition as a viable method to improve performance and broaden the design space of near ignition high adiabat implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Development of a predictive capability of short-pulse laser-driven broadband x-ray radiography

High intensity, short-pulse laser interaction with a solid metal target produces broadband hard x-rays potentially for various applications of x-ray radiography. In this work, experimental benchmarking of numerical modelling for short-pulse laser-driven broadband x-ray radiography is presented. Angular dependent x-ray spectra are first calculated with a hybrid particle-in-cell code, Large Scale Plasma (LSP), using fast electron parameters inferred from an analysis of measured bremsstrahlung signals. Subsequently, a calculated x-ray spectrum in the direction of radiography is used in photon transport calculations using a Monte Carlo code, Particle and Heavy Ion Transport code System (PHITS), to simulate a radiographic image including a modelled 3D test object, an x-ray attenuation filter and an image plate detector. Simulated radiographic images are compared with measurements obtained in an experiment using a 50-TW Leopard short-pulse laser at the University of Nevada Reno. Results show that simulations reproduce the experimental images well for three different attenuation filters (plastic, aluminium, and brass), while 1D transmission profiles for the plastic and aluminium filters are quantitatively in good agreement. The modelling approach established in this work could be used as a predictive tool to simulate radiographic images of complex 3D solid objects at any arbitrary angular position or to optimize experimental components such as the source spectrum, x-ray attenuation filters and a detector type depending on a radiographic object without carrying out radiographic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Monte Carlo Modeling and Design of Photon Energy Attenuation Layers for >10× Quantum Yield Enhancement in Si-Based Hard X-ray Detectors

High-energy (>20 keV) X-ray photon detection at high quantum yield, high spatial resolution, and short response time has long been an important area of study in physics. Scintillation is a prevalent method but limited in various ways. Directly detecting high-energy X-ray photons has been a challenge to this day, mainly due to low photon-to-photoelectron conversion efficiencies. Commercially available state-of-the-art Si direct detection products such as the Si charge-coupled device (CCD) are inefficient for >10 keV photons. Here, we present Monte Carlo simulation results and analyses to introduce a highly effective yet simple high-energy X-ray detection concept with significantly enhanced photon-to-electron conversion efficiencies composed of two layers: a top high-Z photon energy attenuation layer (PAL) and a bottom Si detector. We use the principle of photon energy down conversion, where high-energy X-ray photon energies are attenuated down to ≤10 keV via inelastic scattering suitable for efficient photoelectric absorption by Si. Our Monte Carlo simulation results demonstrate that a 10–30× increase in quantum yield can be achieved using PbTe PAL on Si, potentially advancing high-resolution, high-efficiency X-ray detection using PAL-enhanced Si CMOS image sensors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Impact of solvent forces and broken symmetry on the assembly of designed proteins at a liquid-solid interface

The era of protein design has enabled the creation of hybrid protein-inorganic interfaces, leading to both surface-directed self-assembly of de novo protein architectures and protein-directed formation of inorganic materials. However, the resulting patterns of protein assembly are often unexpected, implying that essential interactions are not accounted for in current design platforms. Here, we use high-speed atomic force microscopy (AFM) analyzed through machine learning to follow the assembly of protein nanorods in aqueous electrolytes on two types of mica exhibiting disparate symmetry elements, which are imprinted on the overlying hydration structure. Using Monte Carlo simulations, we reproduce the observed phases and show that an observed smectic phase, previously thought to be unstable for non-interacting rods in two dimensions, emerges when crystal symmetry introduces a directional bias. The findings demonstrate the importance of incorporating solvent forces as modulated by the hydration structure inherent to interfacial systems when designing protein assemblies at liquid-crystal interfaces. Coupling physics-based simulations that can account for these factors to de novo protein design algorithms can lead to improved design platforms for bio-inspired, hybrid materials.

bioinspired materials↗
DOE OSTI2022

Towards the ionizing radiation induced bond dissociation mechanism in oxygen, water, guanine and DNA fragmentation: a density functional theory simulation

The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy and radio-biology, and critical to space projects. In this study, we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneously direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage and introduces a systematic methodology for fine-tuning the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Controlling the Self-Assembly of DNA Origami Octahedra via Manipulation of Inter-Vertex Interactions

Recent studies have demonstrated novel strategies for the organization of nanomaterials into three-dimensional (3D) ordered arrays with prescribed lattice symmetries using DNA-based self-assembly strategies. In one approach, the nanomaterial is sequestered into DNA origami frames or “material voxels” and then coordinated into ordered arrays based on the voxel geometry and the corresponding directional interactions based on its valency. While the lattice symmetry is defined by the valency of the bonds, a larger-scale morphological development is affected by assembly processes and differences in energies of anisotropic bonds. To facilely model this assembly process, we investigate the self-assembly behavior of hard particles with six interacting vertices via theory and Monte Carlo simulations and exploration of corresponding experimental systems. We demonstrate that assemblies with different 3D crystalline morphologies, but the same lattice symmetry can be formed depending on the relative strength of vertex-to-vertex interactions in orthogonal directions. We observed three distinct assembly morphologies for such systems: cube-like, sheet-like, and cylinder-like. A simple analytical theory inspired by well-established ideas in the areas of protein crystallization, based on calculating the second virial coefficient of patchy hard spheres, captures the simulation results and thus represents a straightforward means of modeling this self-assembly process. To complement the theory and simulations, experimental studies were performed to investigate the assembly of octahedral DNA origami frames with varying binding energies at their vertices. In conclusion, x-ray scattering confirms the robustness of the formed nanoscale lattices for different binding energies, while both optical and electron microscopy imaging validated the theoretical predictions on the dependence of the distinct morphologies of assembled state on the interaction strengths in the three orthogonal directions.

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

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Monte Carlo calculations of the extraction of scintillation light from cryogenic $N$-type GaAs

We report the high scintillation luminosity of n-type GaAs at 10 °K is surprising because (1) with a refractive index of 3.5, escape is inhibited by total internal reflection and (2) narrow-beam experiments at 90 °K report infrared absorption coefficients of several per cm. This paper presents Monte Carlo calculations showing that the high luminosity at 10 °K can be explained if (1) the narrow-beam absorption is almost all optical scattering and (2) the absolute absorption coefficient is below 0.1 per cm. Sixteen surface reflector configurations are simulated for a range of internal scattering and absolute absorption coefficients, and these can guide the design of cryogenic scintillating GaAs targets for the direct detection of dark matter. The discussion section presents a possible infrared scattering mechanism based on the metallic nature of n-type GaAs. A supplement file describes (1) the Monte Carlo program steps in detail and (2) how narrow-beam and integrating sphere experiments can measure the cryogenic optical scattering and absolute absorption coefficients.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Improving Detector Systematic Uncertainties Through Data-Driven Machine Learning

Detector simulation in liquid argon time projection chambers (LArTPCs) is a constant challenge. In particular the modeling of electrons response on wires is highly nontrivial. However, new machine learning techniques exist which can be leveraged to ameliorate these concerns. We present a novel methodology to attempt to learn from cosmic muon data in the ICARUS detector how reconstructed wire signals are influenced by features of the hits such as location, particle direction, angle relative to the wire plane, etc. A model can then be generated which can apply the learned mapping to Monte Carlo events to create a more data-like simulation sample. By creating such a sample we expect to reduce the systematic uncertainties at ICARUS due to our detector modeling.

Hausner, Harry [Fermilab] (ORCID:0000000188932280)↗