Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Using Wannier functions to improve solid band gap predictions in density functional theory

Enforcing a straight-line condition of the total energy upon removal/addition of fractional electrons on eigen states has been successfully applied to atoms and molecules for calculating ionization potentials and electron affinities, but fails for solids due to the extended nature of the eigen orbitals. Here we have extended the straight-line condition to the removal/addition of fractional electrons on Wannier functions constructed within the occupied/unoccupied subspaces. It removes the self-interaction energies of those Wannier functions, and yields accurate band gaps for solids compared to experiments. It does not have any adjustable parameters and the computational cost is at the DFT level. This method can also work for molecules, providing eigen energies in good agreement with experimental ionization potentials and electron affinities. Our approach can be viewed as an alternative approach of the standard LDA+U procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Two mechanisms of nanoparticle generation in picosecond laser ablation in liquids: the origin of the bimodal size distribution

The synthesis of chemically clean and environmentally friendly nanoparticles through pulsed laser ablation in liquids has shown a number of advantages over conventional chemical synthesis methods and has evolved into a thriving research field attracting laboratory and industrial applications. The fundamental understanding of processes leading to the nanoparticle generation, however, still remains elusive. In particular, the origin of bimodal nanoparticle size distributions in femto- and picosecond laser ablation in liquids, where small nanoparticles (several nanometers) with narrow size distribution are commonly observed to coexist with larger (tens to hundreds of nanometers) ones, has not been explained so far. In this paper, joint computational and experimental efforts are applied to understand the mechanisms of nanoparticle formation in picosecond laser ablation in liquids and to explain the bimodal nanoparticle size distributions. The results of a large-scale atomistic simulation reveal the critical role of the dynamic interaction between the ablation plume and the liquid environment, leading to the generation of large nanoparticles through a sequence of hydrodynamic instabilities at the plume-liquid interface and a concurrent nucleation and growth of small nanoparticles in an expanding metal-liquid mixing region. The computational predictions are supported by a series of stroboscopic videography experiments showing the emergence of small satellite bubbles surrounding the main cavitation bubble generated in single pulse experiments. Lastly, carefully timed double pulse irradiation triggers expansion of secondary cavitation bubbles indicating, in accord with the simulation results, the presence of localized sites of laser energy deposition (possibly large nanoparticles) injected into the liquid at the early stage of the bubble formation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact of Detector Simulation in Particle Physics Collider Experiments - Highlights

Detector simulation has become fundamental to the success of modern high-energy physics (HEP) experiments. For example, the Geant4-based simulation applications developed by the ATLAS and CMS experiments played a major role for them to produce physics measurements of unprecedented quality and precision with faster turnaround, from data taking to journal submission, than any previous hadron collider experiment. The material presented here contains highlights of a recent review on the impact of detector simulation in particle physics collider experiments published in Ref. [1]. It includes examples of applications to detector design and optimization, software development and testing of computing infrastructure, and modeling of physics objects and their kinematics. The cost and economic impact of simulation in the CMS experiment is also presented. A discussion on future detector simulation needs, challenges and potential solutions to address them is included at the end.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Accelerating End-to-End Deep Learning for Particle Reconstruction using CMS open data

Machine learning algorithms are gaining ground in high energy physics for applications in particle and event identification, physics analysis, detector reconstruction, simulation and trigger. Currently, most data-analysis tasks at LHC experiments benefit from the use of machine learning. Incorporating these computational tools in the experimental framework presents new challenges. This paper reports on the implementation of the end-to-end deep learning with the CMS software framework and the scaling of the end-to-end deep learning with multiple GPUs. The end-to-end deep learning technique combines deep learning algorithms and low-level detector representation for particle and event identification. We demonstrate the end-to-end implementation on a top quark benchmark and perform studies with various hardware architectures including single and multiple GPUs and Google TPU.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Electronic structures and geometries of the XF{sub 3} (X = Cl, Br, I, At) fluorides

The potential energy surfaces of the group 17 XF{sub 3} (X = Cl, Br, I, At) fluorides have been investigated for the first time with multiconfigurational wave function theory approaches. In agreement with experiment, bent T-shaped C{sub 2v} structures are computed for ClF{sub 3}, BrF{sub 3}, and IF{sub 3}, while we predict that an average D{sub 3h} structure would be experimentally observed for AtF{sub 3}. Electron correlation and scalar relativistic effects strongly reduce the energy difference between the D{sub 3h} geometry and the C{sub 2v} one, along the XF{sub 3} series, and in the X = At case, spin-orbit coupling also slightly reduces this energy difference. AtF{sub 3} is a borderline system where the D{sub 3h} structure becomes a minimum, i.e., the pseudo-Jahn-Teller effect is inhibited since electron correlation and scalar-relativistic effects create small energy barriers leading to the global C{sub 2v} minima, although both types of effects interfere.

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

Mid-infrared signatures of hydroxyl containing water clusters: Infrared laser Stark spectroscopy of OH–H{sub 2}O and OH(D{sub 2}O){sub n} (n = 1-3)

Small water clusters containing a single hydroxyl radical are synthesized in liquid helium droplets. The OH–H{sub 2}O and OH(D{sub 2}O){sub n} clusters (n = 1-3) are probed with infrared laser spectroscopy in the vicinity of the hydroxyl radical OH stretch vibration. Experimental band origins are qualitatively consistent with ab initio calculations of the global minimum structures; however, frequency shifts from isolated OH are significantly over-predicted by both B3LYP and MP2 methods. An effective Hamiltonian that accounts for partial quenching of electronic angular momentum is used to analyze Stark spectra of the OH–H{sub 2}O and OH–D{sub 2}O binary complexes, revealing a 3.70(5) D permanent electric dipole moment. Computations of the dipole moment are in good agreement with experiment when large-amplitude vibrational averaging is taken into account. Polarization spectroscopy is employed to characterize two vibrational bands assigned to OH(D{sub 2}O){sub 2}, revealing two nearly isoenergetic cyclic isomers that differ in the orientation of the non-hydrogen-bonded deuterium atoms relative to the plane of the three oxygen atoms. The dipole moments for these clusters are determined to be approximately 2.5 and 1.8 D for “up-up” and “up-down” structures, respectively. Hydroxyl stretching bands of larger clusters containing three or more D{sub 2}O molecules are observed shifted approximately 300 cm{sup −1} to the red of the isolated OH radical. Pressure dependence studies and ab initio calculations imply the presence of multiple cyclic isomers of OH(D{sub 2}O){sub 3}.

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

Transferable density functional tight binding for carbon, hydrogen, nitrogen, and oxygen: Application to shock compression

A new parameterization for density functional tight binding (DFTB) theory, lanl31, has been developed for molecules containing carbon, hydrogen, nitrogen, and oxygen. Optimal values for the Hubbard Us, on-site energies, and the radial dependences of the bond integrals and repulsive potentials were determined by numerical optimization using simulated annealing to a modest database of ab initio-calculated atomization energies and interatomic forces. The transferability of the optimized DFTB parameterization has been assessed using the CHNO subset of the QM-9 database [R. Ramakrishnan et al., Sci. Data, 1, 140022 (2014)]. These analyses showed that the errors in the atomization energies and interatomic forces predicted by our model are small and in the vicinity of the di erences between density functional theory calculations with di erent basis sets and exchange-correlation functionals. Good correlations between the molecular dipole moments and HOMO-LUMO gaps predicted by lanl31 and the QM-9 data set are also found. Furthermore, the errors in the atomization energies and forces derived from lanl31 are signi cantly smaller than those obtained from the ReaxFF-lg reactive force eld for organic materials [L. Liu et al., J. Phys. Chem. A, 115, 11016 (2011)]. The lanl31 DFTB parameterization for C, H, N, and O has been applied the to the molecular dynamics simulation of the principal Hugoniot of liquid nitromethane, liquid benzene, liquid nitrogen, pentaerythritol tetranitrate, trinitrotoluene, and cyclotetramethylene tetranitramine. The computed and measured Hugoniot loci are in excellent agreement with experiment and we discuss the sensitivity of the loci to the underestimated shock heating that is a characteristic of classical molecular dynamics simulations.

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

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

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

Deficiencies in compression and yield in x-ray-driven implosions

This paper analyzes x-ray–driven implosions that are designed to be less sensitive to 2-D and 3-D effects in hohlraum and capsule physics. Key performance metrics including the burn-averaged ion temperature, hot-spot areal density, and fusion yield are found to agree with simulations where the design adiabat (internal pressure) is multiplied by a factor of 1.4. Furthermore, these results motivate the development of a simple model for interpreting experimental data, which is then used to quantify how improvements in compression could help achieve ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum–classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Chemistry↗
DOE OSTI2024

Monte Carlo Perturbation Analysis of Fuel Temperature Variations in the MCNP Model of the Annular Core Research Reactor

The Annular Core Research Reactor (ACRR) Monte Carlo N-Particle (MCNP) model is used by ACRR reactor operators and experiment designers at Sandia National Laboratories for a variety of computational calculations ranging from reactor kinetics parameter estimates and safety analyses to experimental planning. To understand the dominant source of uncertainty within the MCNP model, perturbations in temperature were applied to individual ACRR MCNP fuel rods. Fuel rod temperatures were randomly sampled from a uniform distribution from operational temperatures to quantify temperature-related uncertainty effects. Stochastic mixing was used to blend the cross sections of the desired temperatures using the MCNP continuous and Thermal Neutron Scattering Treatment [S(α,β)] libraries in ENDF/B-VII.1. Furthermore, this uncertainty analysis produced a 640 row × 640 column correlation and covariance matrix of the neutron energy spectra. Positive covariance was produced around the 1-MeV region and the 0.2-eV region. Correlation was found in the thermal and fast energy regions, but no correlation was observed in the slowing-down energy region because interactions in this region are not dominated by fuel.

ACRR↗
DOE OSTI2015

Impact of pre-plasma on fast electron generation and transport from short pulse, high intensity lasers

Previous experiments and modeling examining the impact of an underdense, pre-formed plasma in laser-plasma interactions have shown that the fast electrons are generated with energies higher than predicted by ponderomotive scaling. In this paper, we report on experiments using the Texas Petawatt high intensity (150 fs, 1.5 × 10 20 W cm -2 ) laser pulse, which were conducted to examine the mechanism for accelerating these high energy electrons. These experiments gauge the impact a controlled low density pre-formed plasma has on electron generation with a shorter time scale than previous experiments, 150–180 fs. Electron temperatures measured via magnetic spectrometer on experiment were found to be independent of preformed plasma. Supplemental computational results using 1D PIC simulations predict that super-ponderomotive electrons are generated inside a potential well in the pre-plasma. Finally, however, while the potential well is established around 150 fs, the electrons require at least an additional 50 fs to be trapped and heated inside it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Self-modulated laser wakefield accelerators as x-ray sources

The development of a directional, small-divergence, and short-duration picosecond x-ray probe beam with an energy greater than 50 keV is desirable for high energy density science experiments. We therefore explore through particle-in-cell (PIC) computer simulations the possibility of using x-rays radiated by betatron-like motion of electrons from a self-modulated laser wakefield accelerator as a possible candidate to meet this need. Two OSIRIS 2D PIC simulations with mobile ions are presented, one with a normalized vector potential a0 = 1.5 and the other with an a0 = 3. We find that in both cases direct laser acceleration (DLA) is an important additional acceleration mechanism in addition to the longitudinal electric field of the plasma wave. Together these mechanisms produce electrons with a continuous energy spectrum with a maximum energy of 300 MeV for a0 = 3 case and 180 MeV in the a0 = 1.5 case. Forward-directed x-ray radiation with a photon energy up to 100 keV was calculated for the a0 = 3 case and up to 12 keV for the a0 = 1.5 case. The x-ray spectrum can be fitted with a sum of two synchrotron spectra with critical photon energies of 13 and 45 keV for the a0 of 3 and critical photon energies of 0.3 and 1.4 keV for a0 of 1.5 in the plane of polarization of the laser. The full width at half maximum divergence angle of the x-rays was 62 × 1.9 mrad for a0 =3and77×3.8mradfora0 =1.5.

42 ENGINEERING↗
DOE OSTI2022

High temperature superconductivity in the candidate phases of solid hydrogen

As the simplest element in nature, unraveling the phase diagram of hydrogen is a primary task for condensed matter physics. As conjectured many decades ago, in the low-temperature and high-pressure part of the phase diagram, solid hydrogen is expected to become metallic with a high superconducting transition temperature. The metallization may occur via band gap closure in the molecular solid or via a transition to the atomic solid. Recently, a few experimental studies pushed the achievable pressures into the 400–500 GPa range. There are strong indications that at some pressure in this range metallization via either of these mechanisms occurs, although there are disagreements between experimental reports. Furthermore, there are multiple good candidate crystal phases that have emerged from recent computational and experimental studies which may be realized in upcoming experiments. Therefore, it is crucial to determine the superconducting properties of these candidate phases. In a recent study, we reported the superconducting properties of the C2/c-24 phase, which we believe to be a strong candidate for metallization via band gap closure. Here, we report the superconducting properties of the Cmca-12, Cmca-4 and I4 1 /amd-2 phases including the anharmonic effects using a Wannier function-based dense k-point and q-point sampling. Here we find that the Cmca-12 phase has a superconducting transition temperature that rises from 86 K at 400 GPa to 212 K at 500 GPa, whereas the Cmca-4 and I4 1 /amd-2 phases show a less pressure-dependent behavior with their T c in the 74–94 K and 307–343 K ranges, respectively. These properties can be used to distinguish between crystal phases in future experiments. Understanding superconductivity in pure hydrogen is also important in the study of high-T c hydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Uncertainty quantification and propagation in CALPHAD modeling

Design is about making decisions bounded by a quantifiable degree of certainty. In the context of alloy design, Integrated Computational Materials Engineering (ICME) provides the framework whereby performance requirements are ultimately transformed into alloy/processing specifications through the combination of (complex) computational models connecting process-structure-property-performance relationships and experiments. Most ICME approaches consider the models used as deterministic and thus do not provide the means to make alloy design decisions with proper confidence measures. At the root of ICME lie CALPHAD models that describe the thermodynamics and phase stability of phases under specific thermodynamic boundary conditions. To date, the vast majority of efforts within the CALPHAD community have been deterministic in that thermodynamic models and the resulting thermodynamic properties and phase diagram features do not explicitly account for the uncertainties inherent in the model formulation or in the experimental/computational data used. In this contribution, we provide an overview of the state of the field. We review major efforts thus far and we then provide a (brief) tutorial on basic concepts of uncertainty quantification and propagation (UQ/UP) in CALPHAD. We discuss the major features of frequentist and Bayesian interpretations of uncertainty and proceed with a discussion of recent case studies in which UQ has been used to parameterize models for the thermodynamic properties of phases. We follow our discussion by presenting frameworks and demonstrating the propagation of uncertainty in thermodynamic properties and phase diagram predictions and briefly discuss how we can use Bayesian frameworks for rigorous model selection as well as for model fusion. Furthermore, we close our contribution by providing context for what has been done and what remains to be accomplished in order to fully embrace the management of uncertainty in CALPHAD modeling, a foundational element of ICME.

42 ENGINEERING↗
DOE OSTI2017

Grid site availability evaluation and monitoring at CMS

The Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) uses distributed grid computing to store, process, and analyse the vast quantity of scientific data recorded every year. The computing resources are grouped into sites and organized in a tiered structure. Each site provides computing and storage to the CMS computing grid. Over a hundred sites worldwide contribute with resources from hundred to well over ten thousand computing cores and storage from tens of TBytes to tens of PBytes. In such a large computing setup scheduled and unscheduled outages occur continually and are not allowed to significantly impact data handling, processing, and analysis. Unscheduled capacity and performance reductions need to be detected promptly and corrected. CMS developed a sophisticated site evaluation and monitoring system for Run 1 of the LHC based on tools of the Worldwide LHC Computing Grid. For Run 2 of the LHC the site evaluation and monitoring system is being overhauled to enable faster detection/reaction to failures and a more dynamic handling of computing resources. Furthermore, enhancements to better distinguish site from central service issues and to make evaluations more transparent and informative to site support staff are planned.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Reflection thermal diffuse x-ray scattering for quantitative determination of phonon dispersion relations

Synchrotron reflection x-ray thermal diffuse scattering (TDS) measurements, rather than previously reported transmission TDS, are carried out at room temperature and analyzed using a formalism based upon second-order interatomic force constants and long-range Coulomb interactions to obtain quantitative determinations of MgO phonon dispersion relations (h) over bar omega(j) (q), phonon densities of states g((h) over bar omega), and isochoric temperature-dependent vibrational heat capacities c v (T). We use MgO as a model system for investigating reflection TDS due to its harmonic behavior as well as its mechanical and dynamic stability. Resulting phonon dispersion relations and densities of states are found to be in good agreement with independent reports from inelastic neutron and x-ray scattering experiments. Temperature-dependent isochoric heat capacities c v (T), computed within the harmonic approximation from (h) over bar omega(j) (q) values, increase with temperature from 0.4 x 10 -4 eV/atom K at 100 K to 1.4 x 10 -4 eV/atom K at 200 K and 1.9 x 10 -4 eV/atom K at 300 K, in excellent agreement with isobaric heat capacity values c p (T) between 4 and 300 K. We anticipate that the experimental approach developed here will be valuable for determining vibrational properties of heteroepitaxial thin films since the use of grazing-incidence (θ ≲ θ c where θ c is the density-dependent critical angle) allows selective tuning of x-ray penetration depths to ≲ 10 nm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗