Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational experiment”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Accuracy of XAS theory for unraveling structural changes of adsorbates: CO on Ni(100)

Studying surface reactions using ultrafast optical pump and X-ray probe experiments relies on accurate calculations of X-ray spectra of adsorbates for correct identification of the spectral signatures and their dynamical evolution. We show that experimental XAS can be well reproduced for the different binding sites in a static prototype system CO/Ni(100) at a standard DFT (GGA) level of theory, using a plane-wave basis and pseudopotentials. This validates its utility in analyzing ultrafast X-ray probe experiments. The accuracy of computed relative core level binding energies is about 0.2 eV, representing a lower limit for which spectral features can be resolved with this method. We also show that the commonly used Z +1 approximation gives very good core binding energy shifts overall. However, we find a discrepancy for CO adsorbed in the hollow site, which we assign to the significantly stronger hybridization in hollow bonding than in on-top.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Lessons learned from first-principles calculations of transition metal oxides

Transition metal oxides (TMOs) are an important class of materials with diverse applications, ranging from memristors to photoelectrochemical cells. First-principles calculations are critical for understanding these complex materials at an atomic level and establishing relationships between atomic and electronic structures, particularly for probing quantities difficult or inaccessible to experiment. Here, we discuss computational strategies used to understand TMOs by focusing on two examples, a photoanode material, BiVO4, and an oxide for low-power electronics, La1−xSrxCoO3. We highlight key aspects required for the modeling of TMOs, namely, the descriptions of how oxygen vacancies, extrinsic doping, the magnetic state, and polaron formation impact their electronic and atomic structures and, consequently, many of the observed properties.

Vo, Hien (ORCID:0000000158315599)↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Photoenzymatic Csp 3 –Csp 3 bond formation via enzyme-templated radical–radical coupling

Cross-couplings are essential reactions in modern chemical synthesis, enabling the rapid construction of complex molecules from simple precursors. Transition metal catalysts are prized for these transformations because their reactivity and selectivity can be tuned via judicious selection of the metal and ligand. Although enzymes offer analogous opportunities for tuning via protein engineering, their application to cross-coupling remains limited, as nature relies on alternative paradigms for building molecular complexity. Here, we report the cross-coupling of alkyl halides and benzylic carboxylic acids using an engineered flavin-dependent lactate monooxygenase—a photoenzyme. The enzyme achieves this feat by exploiting the redox versatility of the flavin cofactor. Stoichiometric experiments, ultrafast spectroscopy, and computational studies support a mechanism in which photoexcited flavin quinone initiates the reaction via oxidative decarboxylation to generate a benzylic radical. The resulting flavin semiquinone can reduce the alkyl halide to form a second organic radical within the protein active site, which rapidly engages in C(sp 3 )–C(sp 3 ) bond formation. A variant was engineered to control the stereochemical outcome of this radical–radical coupling event, highlighting the ability of the protein to alter the energetic barrier for a mechanistic step that is traditionally understood to be near barrierless. This work demonstrates that the scope for nonnative reaction mechanisms in biocatalysis far exceeds previously established bounds and has potential to solve a variety of reactivity challenges in cross-coupling chemistry.

biocatalysis↗
DOE OSTI2021

Separate-Effects Tests for Studying Temperature-Gradient-Driven Cracking in UO 2 Pellets

We report a variety of normal operation and accident scenarios can generate thermal stresses large enough to cause cracking in light-water reactor (LWR) fuel pellets. Cracking of fuel pellets can lead to reduced heat removal, higher centerline temperatures, and localized stress in cladding, all of which impact fuel performance. It is important to experimentally characterize the thermal and mechanical behavior in the pellet before and after cracking to improve cracking models in fuel performance codes. However, in-reactor observation and measurement of cracking is very challenging due to the harsh environment and logistics. Recently, an experimental pellet cracking test stand was developed for separate effects testing of normal operations and accident temperature conditions, using thermal imaging to capture the pellet surface temperature for evaluation of thermal stresses and optical imaging to capture the evolution of cracking in real time. Experiments were performed using depleted uranium dioxide (UO 2 ) pellets, which are useful for collecting data valuable for development and validation of cracking models. A combination of induction and resistance heating was used to create an average temperature gradient of 236°C/cm and 193°C/cm before and after cracking respectively. Characterization of the pellets were done before as well as after cracking. The cracking patterns are somewhat different than those expected in a typical reactor because of the differences in thermal conditions and pellet microstructure. However, if the actual conditions of these experiments are reproduced in computational models, these out-of-pile tests on UO 2 pellets provide relevant data for modeling purposes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The scientific potential and technological challenges of the High-Luminosity Large Hadron Collider program

Here, we present an overview of the High-Luminosity (HL-LHC) program at the Large Hadron Collider (LHC), its scientific potential and technological challenges for both the accelerator and detectors. The HL-LHC program is expected to start circa 2027 and aims to increase the integrated luminosity delivered by the LHC by an order of magnitude at the collision energy of 14 TeV. This requires upgrades to the injector system, accelerator complex and luminosity levelling. The two experiments, ATLAS and CMS, require substantial upgrades to most of their systems in order to cope with the increased interaction rate, and much higher radiation levels than at the current LHC. We present selected examples based on novel ideas and technologies for applications at a hadron collider. Both experiments will replace their tracking systems. We describe the ATLAS pixel detector upgrade featuring novel tilted modules, and the CMS Outer Tracker upgrade with a new module design enabling use of tracks in the level-1 trigger system. CMS will also install state-of-the-art highly segmented calorimeter endcaps. Finally, we describe new picosecond precision timing detectors of both experiments. In addition, we discuss how the upgrades will enhance the physics performance of the experiments, and solve the computing challenges posed by the expected large data sets. The physics program of the HL-LHC is focused on precision measurements probing the limits of the Standard Model (SM) of particle physics and discovering new physics. We present a selection of studies that have been carried out to motivate the HL-LHC program. A central topic of exploration will be the characterization of the Higgs boson. The large HL-LHC data samples will extend the sensitivity of searches for new particles or new interactions whose existence has been hypothesized in order to explain shortcomings of the SM. Finally, we comment on the nature of large scientific collaborations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

GPU coprocessors as a service for deep learning inference in high energy physics

In the next decade, the demands for computing in large scientific experiments are expected to grow tremendously. During the same time period, CPU performance increases will be limited. At the CERN Large Hadron Collider (LHC), these two issues will confront one another as the collider is upgraded for high luminosity running. Alternative processors such as graphics processing units (GPUs) can resolve this confrontation provided that algorithms can be sufficiently accelerated. In many cases, algorithmic speedups are found to be largest through the adoption of deep learning algorithms. We present a comprehensive exploration of the use of GPU-based hardware acceleration for deep learning inference within the data reconstruction workflow of high energy physics. We present several realistic examples and discuss a strategy for the seamless integration of coprocessors so that the LHC can maintain, if not exceed, its current performance throughout its running.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Transient computational and experimental thermal analysis of graphite foam monoblock

A densified graphite foam is being explored for its applicability as plasma-facing material in fusion devices. Three different graphite foam monoblocks are constructed and tested in the Garching Large Divertor Sample Testing Facility. The monoblock samples consist of graphite foam press-fit to a single tube, graphite foam cubes brazed to a single, and graphite foam cubes press-fit to a single tube. The tube is composed of CuCrZr with a steel twisted tape. The monoblocks are exposed to the heat fluxes of 5, 6, and 8 MW/m 2 for 30 s to determine the maximum surface and body temperatures measured with thermocouple for each monoblock design. The press-fit monoblocks are exposed to 8 MW/m 2 for 15 s for 100 cycles to determine the effect of thermal cycling on the contact between the graphite foam and the tube. STAR-CCM+ is used to predict how much the contact between the foam and tubes varies as a result of thermal cycling. In addition, the 6 MW/m 2 loading is modeled in STAR-CCM+ to compare the transient cooldown curves of the computational results to the recorded temperatures at the surface and two different thermocouple locations. This comparison is used to validate the temperature-dependent thermal conductivity and specific heat capacity used in the models of the graphite foam. The computational modeling of the experiment has been used to hypothesize ways to use and improve upon the graphite foam as a suitable material for fusion applications.

Gehrig, Monica↗
DOE OSTI2026

Efficient CP Rounding Using Alternating Least Squares with QR Decomposition

The CANDECOMP/PARAFAC (CP) decomposition is widely used for analyzing multidimensional data, and the alternating least squares (CP-ALS) algorithm is a common method for its computation. CP rounding is the problem of computing a lower-rank CP decomposition of an input already in a higher-rank CP format. While the normal equations (NE) approach in CP-ALS is efficient for the CP rounding problem and frequently used, it becomes unstable in the presence of ill-conditioned subproblems. This paper presents a new QR-based CP-ALS method for CP rounding that preserves both numerical stability and computational efficiency. Here, our experiments show that the proposed method offers significant speedup over a previous QR-based approach and the Tensor Toolbox's NE-based implementation, particularly for higher-order tensors. Furthermore, our approach demonstrates a marked reduction in error for ill-conditioned problems, with error reductions several orders of magnitude smaller compared to the NE-based method, while achieving faster convergence and more accurate solutions. By using a more numerically stable approach, we can solve more problems in reduced working precision, which enables further reduction in time to solution.

CANDECOMP/PARAFAC↗
DOE OSTI2022

fpikachu96/KRED-data

Exploiting nature’s catalysts for non-natural transformations that are inaccessible to chemocatalysis is highly desirable but challenging. On the one hand, the widespread nicotinamide-dependent oxidoreductases have not been utilized for single-electron-transfer-induced bimolecular cross-couplings; on the other, the addition of catalytic asymmetric radical conjugate to terminal alkenes remains a challenge owing to strong racemic background reaction and unselective termination of prochiral radical species. Here we report a chemomimetic biocatalysitic approach for construction of alpha-carbonyl stereocentres via an unnatural intermolecular conjugate addition of N-(acyloxy)phthalimides-derived radicals with acceptor-substituted terminal alkenes, by combination of visible-light excitation and nicotinamide-dependent ketoreductases (KREDs). Based on protein crystal structure, we engineered KREDs via a semi-rational mutagenesis strategy to improve reaction outcomes with a small and high-quality variants library. Mechanistic investigations combining wet experiments, crystallographic studies and computational simulations demonstrate that the repurposed biocatalyst can suppress racemic background reaction and unselected side reactions, yielding enantioselectivity that is challenging to achieve by chemocatalysis.

Huang, Xiaoqiang↗
DOE OSTI2021

Comparison of NMC estimates with trans-stilbene, EJ-309, and He-3 detection systems

Neutron multiplicity counting (NMC) is a technique for the assay of fissile material. In this work, three detection systems are utilized for active interrogation assay of shells of the Rocky Flats shells (highly enriched uranium, 93% 235U) stacked from 13.25-54.92 kg assemblies. The singles ($R_1$) and doubles ($R_2$) rates are calculated with each system to estimate two sample parameters: $M_L$- the leakage multiplication and F - the sample fission rate. The estimated mass, m, is found by dividing F by a constant activity per unit mass. Since we are interrogating HEU, the α ratio of (α; n) to fission neutrons is taken to be zero. The system of equations to calculate these quantities was originally derived. The Neutron Multiplicity 3 He Array Detector (NOMAD) consists of 15 3 He tubes inside polyethylene and is the traditional, capture-based detection system for NMC. The polyethylene moderates incoming neutrons for thermal capture in individual tubes. The low gamma background, discrete capture signals, and high efficiency of the NOMAD are beneficial for NMC. However, the time to slow down neutrons to thermal energies leads to a system die away time on the scale of microseconds. Comparatively, the Rossi-alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAMRODD) and the Organic Scintillator Array (OSCAR) are scatter-based detection systems. RAMRODD consists of 8, 5.08 by 5.08 cm EJ-309 liquid scintillators in 4 pairs, evenly-spaced with 90 degrees separation about the center of each assembly. OSCAR consists of a single array of 12, 5.08 by 5.08 cm trans-stilbene crystals aligned with the center of the assembly. Scatter-based systems detect fast neutrons without any moderation, leading to system die away times on the scale of tens of nanoseconds. This allows much shorter timing gate widths compared to thermal systems, thus increasing counting statistics of true correlated fission events. However, scatter-based systems are susceptible to neutron cross-talk when an incident neutron scatters off one detector and interacts in an adjacent detector, causing two seemingly correlated detection signals. The equations developed adjust for cross-talk to conduct NMC with scatter-based systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Integrated photonic Galton board and its application for photon counting

The Galton board is a desktop probability machine traditionally used to visualize the principles of statistical physics from classical particles. Here, we demonstrate a photonic Galton board that enables on-chip observation of statistics from single-photon interference. The photonic Galton board consists of a directional coupler matrix terminated by an array of superconducting nanowire detectors to provide spatiotemporal resolution. This design also allows for photon-number-resolving capability, making it suitable for high-speed photon counting. Our results demonstrate the compatibility between single-photon detector array and photonic integrated circuits, paving the way for implementing on-chip large-scale quantum optics experiments and photonic quantum computing.

Qin, Hezheng (ORCID:0009000139348800)↗
DOE OSTI2024

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μ ECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Modeling of the amplified optical stochastic cooling experiment at IOTA

Optical Stochastic Cooling (OSC), a beam cooling technique based on Stochastic Cooling, is in the early stages of experimental development. It uses radiation produced by the beam in an undulator magnet (the pickup) to correct the momentum deviation of particles downstream in another undulator (the kicker). The OSC mechanism was recently demonstrated at Fermilab’s IOTA ring using a passive configuration. However, the cooling rate of OSC can be dramatically increased by first amplifying the undulator radiation before applying the corrective kick. In collaboration with the IOTA experiment, we developed a computational model of the OSC mechanism. This paper presents beam-dynamics simulations of the amplified-OSC configuration. We implement a model of intrabeam scattering and study the effects on beam equilibrium and diffusion rate as a function of bunch charge. Finally, we investigate the phase-space dynamics with various coupling configurations between the transverse and longitudinal planes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Constructing a Testing Application for GWMS

Many of Fermilab’s High Energy Physics experiments require High Throughput Computing to carry out simulations, data reconstruction, and data analysis. Glidein Workflow Management System (GWMS) is a tool that distributes High Throughput Computing. The purpose of GWMS is to provide convenient access to Grid resources and sites. GWMS is coupled to HTCondor. HTCondor is the workload management system that is used for scheduling and job control. Users submit jobs to a local HTCondor queue and GWMS will make sure the job will run on one of the many remote resources.

Neely-Brown, LeRayah↗
DOE OSTI2020

Measuring and Simulating T1 and T2 for Qubits

Quantum computers perform computations exploiting quantum mechanics to a possible advantage, allowing us to prepare and manipulate states that do not have a classical equivalent. In particular, phenomena like superposition and entanglement may enable quantum computers to outperform their classical counterparts in certain applications. Implementing these useful quantum algorithms is contingent upon building accurate quantum hardware that is not affected by noise. Environmental noise decreases coherence time of qubits, meaning that qubits do not stay in a desired state long enough to carry out a complex computation. To that end, harnessing the full power of quantum computers necessitate characterization of noise sources and how they impact a given quantum system. Often times, T1and T2 are used to quantify noise. In this project, we provide an approach as to how T1 and T2 values are calculated and simulated for quantum systems. In addition, we compare simulated values of T1 and T2 with those of real quantum computer’s measurements. IBMQ Experience, an open source software allowing users to simulate and use real quantum hardware, is used. QuTip, a Python-based toolbox offering quantum simulation tools for open quantum systems, is also used.

Youssef, Rahaf↗