Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 577 records · Page 32

Toward alcohol synthesis from CO hydrogenation on Cu(111)-supported MoS 2 – predictions from DFT+KMC

In the quest for cheap and efficient catalysts for alcohol synthesis from syngas, a material of interest is single-layer MoS2 owing to its low cost, abundancy, and flexible structure. Because of the inertness of its basal plane, however, it is essential to find ways that make it catalytically active. Herein, by means of density functional theory based calculations of reaction pathways and activation energy barriers and accompanying kinetic Monte Carlo simulations, we show that while S vacancy row structures activate the MoS 2 basal plane, further enhancement of chemical activity and selectivity can be achieved by interfacing the MoS 2 layer with a metallic support. When defect-laden MoS2 is grown on Cu(111), there is not only an increase in the active region (surface area of active sites) but also charge transfer from Cu to MoS 2 , resulting in a shift of the Fermi level such that the frontier states (d orbitals of the exposed Mo atoms) appear close to it, making the MoS 2 /Cu(111) system ready for catalytic activity. Finally, our calculated thermodynamics of reaction pathways lead to the conclusion that the Cu(111) substrate promotes both methanol and ethanol as the products, while kinetic Monte Carlo simulations suggest a high selectivity toward the formation of ethanol.

2D materials↗

Gold-tantalum alloy films deposited by high-density-plasma magnetron sputtering

Gold-tantalum alloy films are of interest for biomedical and magnetically-assisted inertial confinement fusion applications. Here, we systematically study the effects of substrate tilt (0°–80°) and negative substrate bias (0–100 V) on properties of ≲3-μm-thick films deposited by high-power impulse magnetron sputtering (HiPIMS) from a Au–Ta alloy target (with 80 at. % of Ta). Results reveal that, for all the substrate bias values studied, an increase in substrate tilt leads to a monotonic decrease in film thickness, density, residual compressive stress, and electrical conductivity. Larger substrate bias favors the formation of a body-centered cubic phase, with films exhibiting lower column tilt and higher density, electrical conductivity, and residual compressive stress. Furthermore, these changes are attributed to metal atom ionization effects, based on the lack of correlation with distributions of landing energies and incident angles of depositing species as calculated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. By varying substrate tilt and bias in HiPIMS deposition, properties of Au–Ta alloy films can be controlled in a very wide range, including residual stress from –2 to +0.5 GPa, density from 12 to 17 g/cm 3 , and the electrical resistivity from 50 to 4500 μΩ cm, enabling optimum deposition conditions to be selected for specific applications.

36 MATERIALS SCIENCE↗

Hysteresis curves reveal the microscopic origin of cooperative CO 2 adsorption in diamine-appended metal–organic frameworks

Diamine-appended metal-organic frameworks (MOFs) of the form Mg 2 (dobpdc)(diamine) 2 adsorb CO 2 in a cooperative fashion, exhibiting an abrupt change in CO 2 occupancy with pressure or temperature. This change is accompanied by hysteresis. While hysteresis is suggestive of a first-order phase transition, we show that hysteretic temperature-occupancy curves associated with this material are qualitatively unlike the curves seen in the presence of a phase transition; they are instead consistent with CO 2 chain polymerization, within one-dimensional channels in the MOF, in the absence of a phase transition. Here, our simulations of a microscopic model reproduce this dynamics, providing a physical understanding of cooperative adsorption in this industrially important class of materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modulation of intensity emerging from zero effort (MIEZE) with extended Fourier time at large scattering angle

Modulation of Intensity Emerging from Zero Effort (MIEZE) is a neutron resonant spin echo technique that allows one to measure time correlation scattering functions in materials by implementing radio-frequency (RF) intensity modulation at the sample and the detector. The technique avoids neutron spin manipulation between the sample and the detector and, thus, could find applications in cases where the sample depolarizes the neutron beam. However, the finite sample size creates a variance in the path length between the locations where scattering and detection happen, which limits the contrast in intensity modulation that one can detect, in particular, toward long correlation times or large scattering angles. We propose a modification to the MIEZE setup that will enable one to extend those detection limits to longer times and larger angles. Furthermore, we use Monte Carlo simulations of a neutron scattering beamline to show that by tilting the RF flippers in the primary spectrometer with respect to the beam direction, one can shape the wave front of the intensity modulation at the sample to compensate for the path variance from the sample and the detector. The simulation results indicate that this change enables one to operate a MIEZE instrument at much increased RF frequencies, thus improving the effective energy resolution of the technique. For the MIEZE instrument simulated, it shows that for an incident beam with the maximum divergence of 0.33°, the maximum Fourier time can be increased by a factor of 3.

47 OTHER INSTRUMENTATION↗

Oblique angle deposition of boron carbide films by magnetron sputtering

Many applications of boron carbide (B4C) films entail deposition on non-planar substrates, necessitating a better understanding of oblique angle deposition phenomena. Here, we systematically study the effect of substrate tilt on properties of B 4 C films with thicknesses up to 10 μm deposited by direct current magnetron sputtering. Results show that all films are amorphous and columnar with an average column width of ~100 nm, independent of substrate tilt. Column tilt angles are limited to ~20° even for substrate tilt of 80°. Film density, residual stress, and the refractive index weakly (within ≲20%) depend on substrate tilt. Furthermore, oxygen impurities bond preferentially with carbon atoms in inter-columnar regions. Substrate tilt has a major effect on mechanical properties that decrease by ~50%, suggesting weak interconnection between nano-columns. Implications of these observations for the deposition onto non-planar substrates are discussed.

36 MATERIALS SCIENCE↗

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗

Extension of charged-particle spectrometer capabilities for diagnosing implosions on OMEGA, Z, and the NIF

New designs and a new analysis technique have been developed for an existing compact charged-particle spectrometer on the NIF and OMEGA. The new analysis technique extends the capabilities of this diagnostic to measure arbitrarily shaped ion spectra down to 1 MeV with yields as low as 10 6 . Three different designs are provided optimized for the measurement of DD protons, T 3 He deuterons, and 3 He 3 He protons. The designs are highly customizable, and a generalized framework is provided for optimizing the design for alternative applications. Additionally, the understanding of the detector’s response and uncertainties is greatly expanded upon. A new calibration procedure is also developed to increase the precision of the measurements.

47 OTHER INSTRUMENTATION↗

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model validation for inferred hot-spot conditions in National Ignition Facility experiments

Progress toward ignition requires accurately diagnosing current conditions and assessing proximity metrics for implosion experiments on the National Ignition Facility. Hot-spot conditions are not directly measured, but rather inferred, often using simple 0- and 1D models [P. Patel et al., Phys. Plasmas 27, 050901 (2020)]. In this work, we present a detailed accuracy validation exercise using a set of ~20,000 2D simulations encompassing a variety of performance and degradation levels. We find good agreement between the model-inferred pressure and the simulated burn-weighted pressure at peak neutron production and also present results on the precision of inferred quantities using a Markov-Chain Monte Carlo algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-yield and high-angular-fluence neutron generation from deuterons accelerated by laser-driven collisionless shock

Since a bright collimated neutron source is an essential tool for global security missions and fundamental scientific research. In this paper, we study a compact high-yield and high-angular-fluence neutron source particularly suitable for high-energy neutron applications utilizing the breakup reaction of laser-driven deuterons in a 9 Be converter. The neutron generation scaling from such a reaction is used to guide the choice and optimization of the acceleration process for bulk ions in a low density CD 2 foam. In particular, the collisionless shock acceleration mechanism is exploited with proper choice in the laser and target parameter space to accelerate these ions toward energies above the temperature of the distribution. Particle-in-cell and Monte Carlo simulations are coupled here to investigate this concept and possible adverse effects as well as the contribution from the surface ions accelerated and the optimal converter design. The simulation results indicated that our design can be a practical approach to increase both the neutron yield and angular fluence of laser-driven neutron sources, reaching >10 11 neutron/pulse (or >10 8 neutron/J) and >10 11 neutron/sr (or >10 8 neutron/sr/J) with present-day kJ-class high-power lasers. Such developments will advance fundamental neutron science, high precision radiography, and other global security applications with laser-driven sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A scintillator attenuation spectrometer for intense gamma-rays

A new type of compact high-resolution high-sensitivity gamma-ray spectrometer for short-pulse intense gamma-rays (250 keV to 50 MeV) has been developed by combining the principles of scintillators and attenuation spectrometers. The first prototype of this scintillator attenuation spectrometer (SAS) was tested successfully in Trident laser experiments at LANL. Later versions have been used extensively in the Texas Petawatt laser experiments in Austin, TX, and more recently in OMEGA-EP laser experiments at LLE, Rochester, NY. The SAS is particularly useful for high-repetition-rate laser applications. Furthermore, we give a concise description of the design principles, capabilities, and sample preliminary results of the SAS.

47 OTHER INSTRUMENTATION↗

Small angle scattering of diblock copolymers profiled by machine learning

We outline a machine learning strategy for quantitively determining the conformation of AB-type diblock copolymers with excluded volume effects using small angle scattering. Complemented by computer simulations, a correlation matrix connecting conformations of different copolymers according to their scattering features is established on the mathematical framework of a Gaussian process, a multivariate extension of the familiar univariate Gaussian distribution. We show that the relevant conformational characteristics of copolymers can be probabilistically inferred from their coherent scattering cross sections without any restriction imposed by model assumptions. This work not only facilitates the quantitative structural analysis of copolymer solutions but also provides the reliable benchmarking for the related theoretical development of scattering functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PIONEER, a high-resolution single-crystal polarized neutron diffractometer

PIONEER is a high Q-resolution, single-crystal, polarized neutron diffractometer at the Second Target Station (STS), Oak Ridge National Laboratory. It will provide the unprecedented capability of measuring tiny crystals (0.001 mm 3 , i.e., x-ray diffraction size), ultra-thin films (10 nm thickness), and weak structural and magnetic transitions. PIONEER benefits from the increased peak brightness of STS cold-neutron sources and uses advanced Montel mirrors that are able to deliver a focused beam with a high brilliance transfer, a homogeneous profile, and a low background. Monte Carlo simulations suggest that the optimized instrument has a high theoretical peak brilliance of 2.9 × 10 12 n cm -2 sr -1 Å -1 s -1 at 2.5 Å at the sample position, within a 5 × 5 mm 2 region and a ±0.3° divergence range. The moderator-to-sample distance is 60 m, providing a nominal wavelength band of 4.3 Å with a wavelength resolution better than 0.2% in the wavelength range of 1.0–6.0 Å. PIONEER is capable of characterizing large-scale periodic structures up to 200 Å. With a sample-to-detector distance of 0.8 m, PIONEER accommodates various sample environments, including low/high temperature, high pressure, and high magnetic/electric field. Additionally, a large cylindrical detector array (4.0 sr) with a radial collimator is planned to suppress the background scattering from sample environments. Bottom detector banks provide an additional 0.4 sr coverage or can be removed if needed to accommodate special sample environments. We present virtual experimental results to demonstrate the scientific performance of PIONEER in measuring tiny samples.

47 OTHER INSTRUMENTATION↗

Selective amorphization of SiGe in Si/SiGe nanostructures via high energy Si+ implant

Here, the selective amorphization of SiGe in Si/SiGe nanostructures via a 1 MeV Si + implant was investigated, resulting in single-crystal Si nanowires (NWs) and quantum dots (QDs) encapsulated in amorphous SiGe fins and pillars, respectively. The Si NWs and QDs are formed during high-temperature dry oxidation of single-crystal Si/SiGe heterostructure fins and pillars, during which Ge diffuses along the nanostructure sidewalls and encapsulates the Si layers. The fins and pillars were then subjected to a 3 × 10 15 ions/cm 2 1 MeV Si + implant, resulting in the amorphization of SiGe, while leaving the encapsulated Si crystalline for larger, 65-nm wide NWs and QDs. Interestingly, the 26-nm diameter Si QDs amorphize, while the 28-nm wide NWs remain crystalline during the same high energy ion implant. This result suggests that the Si/SiGe pillars have a lower threshold for Si-induced amorphization compared to their Si/SiGe fin counterparts. However, Monte Carlo simulations of ion implantation into the Si/SiGe nanostructures reveal similar predicted levels of displacements per cm 3 . Molecular dynamics simulations suggest that the total stress magnitude in Si QDs encapsulated in crystalline SiGe is higher than the total stress magnitude in Si NWs, which may lead to greater crystalline instability in the QDs during ion implant. The potential lower amorphization threshold of QDs compared to NWs is of special importance to applications that require robust QD devices in a variety of radiation environments.

36 MATERIALS SCIENCE↗

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phased plan for the implementation of the time-resolving magnetic recoil spectrometer on the National Ignition Facility (NIF)

The time-resolving magnetic recoil spectrometer (MRSt) is a transformative diagnostic that will be used to measure the time-resolved neutron spectrum from an inertial confinement fusion implosion at the National Ignition Facility (NIF). It uses a CD foil on the outside of the hohlraum to convert fusion neutrons to recoil deuterons. An ion-optical system positioned outside the NIF target chamber energy-disperses and focuses forward-scattered deuterons. A pulse-dilation drift tube (PDDT) subsequently dilates, un-skews, and detects the signal. While the foil and ion-optical system have been designed, the PDDT requires more development before it can be implemented. Therefore, a phased plan is presented that first uses the foil and ion-optical systems with detectors that can be implemented immediately—namely CR-39 and hDISC streak cameras. These detectors will allow the MRSt to be commissioned in an intermediate stage and begin collecting data on a reduced timescale, while the PDDT is developed in parallel. A CR-39 detector will be used in phase 1 for the measurement of the time-integrated neutron spectra with excellent energy-resolution, necessary for the energy calibration of the system. Streak cameras will be used in phase 2 for measurement of the time-resolved spectrum with limited spectral coverage, which is sufficient to diagnose the time-resolved ion temperature. Simulations are presented that predict the performance of the streak camera detector, indicating that it will achieve excellent burn history measurements at current yields, and good time-resolved ion-temperature measurements at yields above 3 × 10 17 . The PDDT will be used for optimal efficiency and resolution in phase 3.

47 OTHER INSTRUMENTATION↗

Quantifying electron temperature distributions from time-integrated x-ray emission spectra

K-shell x-ray emission spectroscopy is a standard tool used to diagnose the plasma conditions created in high-energy-density physics experiments. In the simplest approach, the emissivity-weighted average temperature of the plasma can be extracted by fitting an emission spectrum to a single temperature condition. It is known, however, that a range of plasma conditions can contribute to the measured spectra due to a combination of the evolution of the sample and spatial gradients. In this work, we define a parameterized model of the temperature distribution and use Markov Chain Monte Carlo sampling of the input parameters, yielding uncertainties in the fit parameters to assess the uniqueness of the inferred temperature distribution. Here we present the analysis of time-integrated S and Fe x-ray spectroscopic data from the Orion laser facility and demonstrate that while fitting each spectral region to a single temperature yields two different temperatures, both spectra can be fit simultaneously with a single temperature distribution. We find that fitting both spectral regions together requires a maximum temperature of $1310^{+90}_{-70}$ eV with significant contributions from temperatures down to 200 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a multi-detector, single line-of-sight, time-of-flight system to measure time-resolved neutron energy spectra

In the dynamic environment of burning, thermonuclear deuterium–tritium plasmas, diagnosing the time-resolved neutron energy spectrum is of critical importance. Strategies exist for this diagnosis in magnetic confinement fusion plasmas, which presently have a lifetime of ~10 12 longer than inertial confinement fusion (ICF) plasmas. In this work, we present a novel concept for a simple, precise, and scale-able diagnostic to measure time-resolved neutron spectra in ICF plasmas. The concept leverages general tomographic reconstruction techniques adapted to time-of-flight parameter space, and then employs an updated Monte Carlo algorithm and National Ignition Facility-relevant constraints to reconstruct the time-evolving neutron energy spectrum. Reconstructed spectra of the primary 14.028 MeV n DT peak are in good agreement with the exact synthetic spectra. The technique is also used to reconstruct the time-evolving downscattered spectrum, although the present implementation shows significantly more error.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗