Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurable radiation levels around individuals externally contaminated by nuclear fallout

During the early response to large-scale radioactive contamination events, people who are potentially affected need to be screened for radioactive contamination and public health staff need to triage individuals who may need immediate decontamination. This is typically done by screening individuals for external contamination using ionising radiation detection equipment. In this study, spatially and temporally dependent isotopic compositions from a simulated nuclear detonation and Monte Carlo methods were used to relate contamination activity levels to the measurable radiation levels at select distances away from an individual with whole-body contamination. Radionuclide-specific air kerma rate coefficients and Geiger–Mueller instrument response coefficients at five select distances from contaminated individuals are presented for 662 radionuclides. Temporally and spatially dependent incident-specific coefficients are presented for a hypothetical surface detonation of a 235 U-fueled device.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Modeling Plutonium Decorporation in a Female Nuclear Worker Treated with Ca-DTPA after Inhalation Intake

The present work models plutonium (Pu) biokinetics in a female former nuclear worker. Her bioassay measurements are available at the US Transuranium and Uranium Registries. The worker was internally exposed to a plutonium-americium mixture via acute inhalation at a nuclear weapons facility. She was medically treated with injections of 1 g Ca-DTPA on days 0, 5, and 14 after the intake. Between days 0 and 20, fecal and urine samples were collected and analyzed for 239 Pu and 241 Am. Subsequently, she was followed up for bioassay monitoring over 14 y, with additional post-treatment urine samples collected and analyzed for 239 Pu. The uniqueness of this dataset is due to the availability of: (1) both early and long-term bioassay data from a female with plutonium intake; (2) data on chelation therapy for a female; and (3) fecal measurement results. Chelation therapy with Ca- and/or Zn-salts of DTPA is known to aid in reducing the internal radiation dose by enhancing the excretion of plutonium and americium from the body. Such enhancement affects plutonium biokinetics in the human body, posing a challenge to the internal dose assessment. The current radiation dose assessment practice is to exclude the data affected by Ca-DTPA from the analysis. The present analysis is the first to explicitly model the chelation-affected bioassay data in a female by using a newly developed chelation model. Thus, the bioassay data collected during and after the Ca-DTPA administrations were used for biokinetic modeling and dose assessment. The Markov Chain Monte Carlo method was used to investigate model parameter uncertainty, based on the bioassay data and assumed prior probability distributions. A χ 2 /nData (number of data points) ≈ 1 was observed in this study, which indicates self-consistency of the data with the model. Results of this study show that the worker’s 239 Pu intake was 12 Bq, with a committed effective dose to the whole-body of 1.2 mSv and a committed equivalent dose to the bone surfaces, liver, and lungs of 37.8, 9.1, and 0.8 mSv, respectively. This study also discusses the worker’s dose reduction due to chelation treatment.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Effect of the orientation of Rydberg atoms on their collisional ionization cross section

Collisional ionization between two Rydberg atoms in relative motion is examined. A classical trajectory Monte Carlo method is used to determine the cross sections associated with Penning ionization. The dependence of the ionization cross section on the magnitude and the direction of orbital angular momentum of the electrons and the direction of the Laplace-Runge-Lenz vector of the electrons is studied. For a given magnitude of angular momentum, there can exist a difference of a factor of up to ~2.5 in the ionization cross section between the orientations with the highest and the lowest ionization cross section. As a result, the case of exchange ionization is examined and its dependence on the magnitude of angular momentum is studied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Thermodynamics and magnetism in the two-dimensional to three-dimensional crossover of the Hubbard model

The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using the determinant quantum Monte Carlo method, a simple realization of anisotropy in the three-dimensional (3D) Hubbard model in which the tunneling between planes, $t_⊥$, is unequal to the intraplane tunneling $\textit{t}$. This model interpolates between the three-dimensional isotropic ($t_⊥ = t$) and two-dimensional (2D;$t_⊥$ = 0) systems. In this work, we show that at fixed interaction strength to tunneling ratio ($\textit{U/t}$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Néel temperature $T_{\text{Néel}}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $\textit{U/t}$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropy's dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Numerically exact full counting statistics of the energy current in the Kondo regime

In this study, we use the inchworm quantum Monte Carlo method to investigate the full counting statistics of particle and energy currents in a strongly correlated quantum dot. Our method is used to extract the heat fluctuations and entropy production of a quantum thermoelectric device, as well as cumulants of the particle and energy currents. The energy-particle current cross correlations reveal information on the preparation of the system and the interplay of thermal and electric currents. Finally, we demonstrate the signature of a crossover from Coulomb blockade to Kondo physics in the energy current fluctuations, and show how the conventional master equation approach to full counting statistics systematically fails to capture this crossover.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Charge density wave order on a $\textit{π}$-flux square lattice

The effect of electron-phonon coupling on Dirac fermions was recently explored numerically on a honeycomb lattice, leading to precise quantitative values for the finite temperature and quantum critical points. In this paper, we use the unbiased determinant quantum Monte Carlo method to study the Holstein model on a half-filled staggered-flux square lattice and compare the results with the honeycomb lattice geometry, presenting results for a range of phonon frequencies 0.1 ≤ $\omega$ ≤ 2.0. In this work, we find that the interactions give rise to charge density wave order, but only above a finite coupling strength $\gamma_{\text{crit}}$. The transition temperature is evaluated and presented in a $T_{c^{–λ}}$ phase diagram. An accompanying mean-field theory calculation also predicts the existence of a quantum phase transition, but at a substantially smaller coupling strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Witnessing entanglement in quantum magnets using neutron scattering

We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF 3 . We apply three entanglement witnesses—one tangle, two tangle, and quantum Fisher information—to its inelastic neutron spectrum and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find that each witness provides direct access to entanglement. Of these, quantum Fisher information is the most robust experimentally and indicates the presence of at least bipartite entanglement up to at least 50 K, corresponding to around 10% of the spinon zone-boundary energy. We apply quantum Fisher information to higher spin-S Heisenberg chains and show theoretically that the witnessable entanglement gets suppressed to lower temperatures as the quantum number increases. Finally, we outline how these results can be applied to higher dimensional quantum materials to witness and quantify entanglement.

1-dimensional systems↗
DOE OSTI2024

Magnetic, charge, and bond order in the two-dimensional Su-Schrieffer-Heeger-Holstein model

Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. As a result, we map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Ab initio many-body study of cobalt adatoms adsorbed on graphene

Many recent calculations have been performed to study a Co atom adsorbed on graphene, with significantly varying results on the nature of the bonding. In this study, we use the auxiliary-field quantum Monte Carlo method and a size-correction embedding scheme to accurately calculate the binding energy of Co on graphene. We find that as a function of the distance $h$ between the Co atom and the sixfold hollow site, there are three distinct ground states corresponding to three electronic configurations of the Co atom. Two of these states provide binding and exhibit a double-well feature with nearly equal binding energy of 0.4 eV at $h = 1.51$ and $h = 1.65$ Å, corresponding to low-spin 2 Co($3d^94s^0$) and high-spin 4 Co($3d^84s^1$), respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Entanglement properties of the antiferromagnetic-singlet transition in the Hubbard model on bilayer square lattices

Here, we calculate the bipartite R enyi entanglement entropy of an L x L x 2 bilayer Hubbard model using a determinantal quantum Monte Carlo method recently proposed by Grover [Phys. Rev. Lett. 111, 130402 (2013)]. Two types of bipartition are studied: (i) One that divides the lattice into two L x L planes, and (ii) One that divides the lattice into two equal-size (L x L=2 x 2) bilayers. Furthermore, we compare our calculations with those for the tight-binding model studied by the correlation matrix method. As expected, the entropy for bipartition (i) scales as L 2 , while the latter scales with L with possible logarithmic corrections. The onset of the antiferromagnet to singlet transition shows up by a saturation of the former to a maximal value and the latter to a small value in the singlet phase. We also comment on the large uncertainties in the numerical results with increasing U, which would have to be overcome before the critical behavior and logarithmic corrections can be quanti ed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Bond patterns and charge-order amplitude in quarter-filled charge-transfer solids

Most quasi-one-dimensional (quasi-1D) quarter-filled organic charge-transfer solids (CTS) with insulating ground states have two thermodynamic transitions: a high-temperature metal-insulator transition followed by a low-temperature magnetic transition. This sequence of transitions can be understood within the 1D Peierls-extended Hubbard (PEH) model. However, in some quasi-1D CTS both transitions occur simultaneously in a direct metal to spin-gapped insulator transition. In this second class of materials the organic stack bond distortion pattern does not follow the pattern of a second dimerization of a dimer lattice. These materials also display charge ordering of a large amplitude below the transition. Using quantum Monte Carlo methods we show that the same PEH model can be used to understand both classes of materials, however, within different parameter regions. We discuss the relevance of our work to experiments on several quarter-filled conductors, focusing in particular on the materials (EDO-TTF) 2 X and (DMEDO-TTF) 2 X.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerically exact full counting statistics of the nonequilibrium Anderson impurity model

The time-dependent full counting statistics of charge transport through an interacting quantum junction is evaluated from its generating function, controllably computed with the inchworm Monte Carlo method. Exact noninteracting results are reproduced; then, we continue to explore the effect of electron-electron interactions on the time-dependent charge cumulants, first-passage time distributions, and n-electron transfer distributions. We observe a crossover in the noise from Coulomb blockade to Kondo-dominated physics as the temperature is decreased. In addition, we uncover long-tailed spin distributions in the Kondo regime and analyze queuing behavior caused by correlations between single-electron transfer events

36 MATERIALS SCIENCE↗
DOE OSTI2018

Mixing A -type and G -type B -site antiferromagnetism in A M n 1 - x F e x O 3 ( A = La , Nd )

Using a classical Heisenberg model with anisotropic nearest-neighbor superexchange and uniaxial anisotropy, I investigate the magnetic order of the transition metals on the B sites in AMn 1-x Fe x O 3 (A=La, Nd; x=[0,1]) with Monte Carlo methods. Magnetic parameters are extracted from preexisting experimental data when available and otherwise estimated with calculations. Superexchange energies are compared to density functional theory and anisotropy energies to ligand field theory, and the potential for doping-induced strain affecting these parameters is considered. I find coexistence regions of different types of magnetic order (A-type and G-type antiferromagnetism), regions with distinct magnetic transitions for the two order parameters, noncollinear ground states for intermediate values of x that decrease the amplitude of the ordered moment in a way consistent with experiment for low dopings, and a reproduction of the experimentally determined, nonmonotonic interpolation of the ordering temperatures. The detailed shape of the phase boundaries in (x,T) space is found to be diagnostic of the magnetic interactions. Finally, these results also suggest a ferrimagnetic ground state for the A 2 MnFeO 6 double perovskite with T N ≈270K for A=La and Nd.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Indirect determination of neutron-capture cross sections for Sm isotopes

Neutron-capture cross sections in the energy range of 0.01–10 MeV for Sm 146 , 148 , 150 , 152 nuclei were calculated using the γ -ray strength and level density functions extracted from the Oslo-type experiments. The uncertainties in the cross-section values were determined using a Monte Carlo method. For the Sm 148 , 150 , 152 isotopes, the calculated cross sections are in a good agreement with the existing experimental data and for the Sm 146 nucleus, an experimental ( n , γ ) cross section is reported for the first time. The results are compared with the ENDF, EAF, and TENDL evaluations. Maxwellian-averaged cross sections were also calculated using the same input γ SF and LD functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗