Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing 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 919 records · Page 51

Predissociation-based measurements of bond dissociation energies: US 2 , OUS, and USe

The uranium-containing molecules US 2 , OUS, and USe have been investigated using a pulsed laser ablation supersonic beam molecular source with time-of-flight mass spectrometric detection. Spectra have been recorded using the resonant two-photon ionization method over the spectroscopic range from 277 to 238 nm. These species have a myriad of excited electronic states in this spectroscopic region, leading to spectra that are highly congested and appear quasicontinuous. Sharp predissociation thresholds are observed, allowing precise bond dissociation energies to be measured. In the case of the triatomic molecules, it was necessary to use one laser for excitation and a delayed laser for ionization in order to observe a sharp predissociation threshold that allowed a precise bond dissociation energy to be measured. The resulting thermochemical values are D 0 (SU-S) = 4.910 ± 0.003 eV, D 0 (OU-S) = 5.035 ± 0.004 eV, and D 0 (USe) = 4.609 ± 0.009 eV. These results provide the first measurement of D 0 (USe) and reduce the error limits in the previous values of D 0 (SU-S) and D 0 (OU-S) by a factor of more than 70.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

nautilus : boosting Bayesian importance nested sampling with deep learning

ABSTRACT We introduce a novel approach to boost the efficiency of the importance nested sampling (INS) technique for Bayesian posterior and evidence estimation using deep learning. Unlike rejection-based sampling methods such as vanilla nested sampling (NS) or Markov chain Monte Carlo (MCMC) algorithms, importance sampling techniques can use all likelihood evaluations for posterior and evidence estimation. However, for efficient importance sampling, one needs proposal distributions that closely mimic the posterior distributions. We show how to combine INS with deep learning via neural network regression to accomplish this task. We also introduce nautilus, a reference open-source python implementation of this technique for Bayesian posterior and evidence estimation. We compare nautilus against popular NS and MCMC packages, including emcee, dynesty, ultranest, and pocomc, on a variety of challenging synthetic problems and real-world applications in exoplanet detection, galaxy SED fitting and cosmology. In all applications, the sampling efficiency of nautilus is substantially higher than that of all other samplers, often by more than an order of magnitude. Simultaneously, nautilus delivers highly accurate results and needs fewer likelihood evaluations than all other samplers tested. We also show that nautilus has good scaling with the dimensionality of the likelihood and is easily parallelizable to many CPUs.

97 MATHEMATICS AND COMPUTING↗

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗

Optical and magneto-optical properties of ferromagnetic monolayer CrBr 3 : A first-principles G W and G W plus Bethe-Salpeter equation study

The discovery of atomically thin two-dimensional (2D) magnetic semiconductors has triggered enormous research interest recently. In this paper, we use first-principles many-body perturbation theory to study a prototypical 2D ferromagnetic semiconductor, monolayer chromium tribromide (CrBr 3 ). With broken time-reversal symmetry, spin-orbit coupling, and excitonic effects included through the full-spinor GW and GW plus Bethe-Salpeter equation (GW-BSE) methods, we compute the frequency-dependent layer polarizability tensor and dielectric function tensor that govern the optical and magneto-optical (MO) properties. In addition, we provide a detailed theoretical formalism for simulating magnetic circular dichroism, MO Kerr effect, and Faraday effect, demonstrating the approach with monolayer CrBr 3 . Due to reduced dielectric screening in 2D and the localized nature of the Cr d orbitals, we find strong self-energy effects on the quasiparticle band structure of monolayer CrBr 3 that give a 3.8 eV indirect bandgap. Also, excitonic effects dominate the low-energy optical and MO responses in monolayer CrBr 3 where a large exciton binding energy of 2.3 eV is found for the lowest bright exciton state with excitation energy at 1.5 eV. We further find that the MO signals demonstrate strong dependence on the excitation frequency and substrate refractive index. Furthermore, our theoretical framework for modeling optical and MO effects could serve as a powerful theoretical tool for future study of optoelectronic and spintronics devices consisting of van der Waals 2D magnets.

36 MATERIALS SCIENCE↗

Integrated analysis of X-ray diffraction patterns and pair distribution functions for machine-learned phase identification

Abstract To bolster the accuracy of existing methods for automated phase identification from X-ray diffraction (XRD) patterns, we introduce a machine learning approach that uses a dual representation whereby XRD patterns are augmented with simulated pair distribution functions (PDFs). A convolutional neural network is trained directly on XRD patterns calculated using physics-informed data augmentation, which accounts for experimental artifacts such as lattice strain and crystallographic texture. A second network is trained on PDFs generated via Fourier transform of the augmented XRD patterns. At inference, these networks classify unknown samples by aggregating their predictions in a confidence-weighted sum. We show that such an integrated approach to phase identification provides enhanced accuracy by leveraging the benefits of each model’s input representation. Whereas networks trained on XRD patterns provide a reciprocal space representation and can effectively distinguish large diffraction peaks in multi-phase samples, networks trained on PDFs provide a real space representation and perform better when peaks with low intensity become important. These findings underscore the importance of using diverse input representations for machine learning models in materials science and point to new avenues for automating multi-modal characterization.

36 MATERIALS SCIENCE↗

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

A dynamic Bayesian optimized active recommender system for curiosity-driven partially Human-in-the-loop automated experiments

Optimization of experimental materials synthesis and characterization through active learning methods has been growing over the last decade, with examples ranging from measurements of diffraction on combinatorial alloys at synchrotrons, to searches through chemical space with automated synthesis robots for perovskites. In virtually all cases, the target property of interest for optimization is defined a priori with the ability to shift the trajectory of the optimization based on human-identified findings during the experiment is lacking. Thus, to highlight the best of both human operators and AI-driven experiments, here we present the development of a human–AI collaborated experimental workflow, via a Bayesian optimized active recommender system (BOARS), to shape targets on the fly with human real-time feedback. Here, the human guidance overpowers AI at early iteration when prior knowledge (uncertainty) is minimal (higher), while the AI overpowers the human during later iterations to accelerate the process with the human-assessed goal. We showcase examples of this framework applied to pre-acquired piezoresponse force spectroscopy of a ferroelectric thin film, and in real-time on an atomic force microscope, with human assessment to find symmetric hysteresis loops. It is found that such features appear more affected by subsurface defects than the local domain structure. This work shows the utility of human–AI approaches for curiosity driven exploration of systems across experimental domains.

36 MATERIALS SCIENCE↗

Importance of charge self-consistency in first-principles description of strongly correlated systems

First-principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first-principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. The success of such embedding theories, often referred to as second-principles, is commonly measured by the quality of self-energy Σ which is either a function of energy or momentum or both. However, Σ should, in principle, also modify the electronic eigenfunctions and thus change the real space charge distribution. While such practices are not prevalent, some works that use embedding techniques do take into account these effects. In such cases, choice of partitioning, of the parameters defining the correlated Hamiltonian, of double-counting corrections, and the adequacy of low-level Hamiltonian hosting the correlated subspace hinder a systematic and unambiguous understanding of such effects. Further, for a large variety of correlated systems, strong correlations are largely confined to the charge sector. Then an adequate nonlocal low-order theory is important, and the high-order local correlations embedding contributes become redundant. Here we study the impact of charge self-consistency within two example cases, TiSe 2 and CrBr 3 , and show how real space charge re-distribution due to correlation effects taken into account within a first-principles Green’s function-based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

36 MATERIALS SCIENCE↗

Target Selection and Sample Characterization for the DESI LOW-Z Secondary Target Program

Abstract We introduce the DESI LOW- Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine-learning methods, to target low-redshift dwarf galaxies ( z < 0.03) between 19 < r < 21 with high completeness. We employ a convolutional neural network (CNN) to select high-priority targets. The LOW- Z survey has already obtained over 22,000 redshifts of dwarf galaxies ( M * < 10 9 M ⊙ ), comparable to the number of dwarf galaxies discovered in the Sloan Digital Sky Survey DR8 and GAMA. As a spare fiber survey, LOW- Z currently receives fiber allocation for just ∼50% of its targets. However, we estimate that our selection is highly complete: for galaxies at z < 0.03 within our magnitude limits, we achieve better than 95% completeness with ∼1% efficiency using catalog-level photometric cuts. We also demonstrate that our CNN selections z < 0.03 galaxies from the photometric cuts subsample at least 10 times more efficiently while maintaining high completeness. The full 5 yr DESI program will expand the LOW- Z sample, densely mapping the low-redshift Universe, providing an unprecedented sample of dwarf galaxies, and providing critical information about how to pursue effective and efficient low-redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Seeing the whole picture: Methods for getting the most from micro X-ray computed tomography of TRISO nuclear fuel particles

Tristructural isotropic (TRISO) coated fuel particles are a nuclear fuel form under extensive study for use in advanced nuclear reactor concepts. TRISO fuels are subjected to high temperature neutron irradiations and then examined to assess their performance by determining fission product retention and studying morphological changes. Micro X-ray computed tomography is one method of nondestructively studying the effects of TRISO performance. This work addresses the need for image processing to remove X-ray tomographic reconstruction artifacts that prevent the study of TRISO features, as the TRISO particles’ high Z kernel can introduce metal artifacts that degrade the image quality in the surrounding low Z coating layers. These metal artifacts were reduced by imaging the TRISO particles with both high- and low-energy X-rays and applying a mask to the radiographs obtained with low-energy X-rays to digitally remove the dense fuel kernel region. These masked radiographs were then used to produce a tomographic reconstruction which was combined with the tomographic reconstruction of the high-energy data. This enabled the relatively-low-density TRISO buffer layer to be examined in more detail, providing information on irradiation induced dimensional changes of the coatings. This methodology, which helps see the full picture of a TRISO particle, is not limited to nuclear fuels but can be applied to systems that contain highly attenuating material surrounded by less dense materials.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

In‐Silico Device Performance Prediction of Cosensitizer Dye Pairs for Dye‐Sensitized Solar Cells

Abstract Endeavors in the field of dye‐sensitized solar cells (DSCs) have shown great promise when adopting a data‐driven approach to materials discovery, such as successful molecular‐scale predictions of light‐harvesting chromophores. However, predictions of DSC dyes would become much more sophisticated if a molecular‐to‐macroscopic DSC device prediction methodology existed. Thereby, a fully computational pipeline is presented that predicts device‐performance parameters of DSCs which contain varying dye combinations. Optimal pairing of complementary dyes is identified via a data‐driven workflow that affords cosensitized DSCs with maximum power‐conversion efficiencies. Six high‐performing DSC dyes are paired with partner dyes that are screened from a database of 8488 compounds using sequential heuristic filters. Existing models that predict short‐circuit‐current density ( J SC ) and open‐circuit voltage ( V OC ) parameters are adapted to predict singly sensitized and cosensitized DSC performance. The predictions for J sc values of singly sensitized devices match experimental literature values with comparable accuracy to more computationally costly methods. Five out of six dye pairings are predicted to have greater J SC values when cosensitized compared to their corresponding singly sensitized devices, including two pairs that show strong J sc boosts of +13% and +12% when cosensitized. Thus, the prospect of an entirely in‐silico prediction pipeline for DSC performance that can be used to realize the fully automated design of optimized cosensitized DSCs is demonstrated.

14 SOLAR ENERGY↗

A remapping scheme for history-dependent material state variables

Here, we present an interpolation scheme for use with material models involving history-dependent state variables. The interpolation scheme is based on the solution of a constrained optimization problem, to identify the smoothest monotone curve which produces an admissible solution. Unlike classical schemes used for transport equations, the proposed methodology accounts for the interdependence of the history variables and produces admissible interpolations of the material's state variables. After introducing a general solution scheme for the problem, a simplified and more computationally efficient method is derived. Finally, as a validation problem for the proposed methods, we consider the advection of internal state variables belonging to material models with highly nonlinear yield surfaces.

42 ENGINEERING↗

Using Neural Networks to Identify Mixture Components in Hyperspectral Reflectance Data

Neural networks have been employed to identify materials of interest from hyperspectral data (generally imagery) based on their unique spectral signatures. This approach assumes that there is a single material that is standing out from the rest of the spectrum to be identified. However, pixels often contain more than one material, or a material of interest may itself be a mixture of multiple materials. Neural networks are only as good as the data used to train them, and it takes a great deal of work in the laboratory to identify, make, and measure all potential mixtures of interest. Thus, researchers often calculate synthetic spectra using algorithms with varying degrees of fidelity to the physics that govern the interactions between light and multiple materials. In this work, we have (1) adapted a neural network designed to identify mixture components from Raman spectroscopy to work with visible to near‐infrared reflectance data and (2) tested three common mixture algorithms to determine the most accurate and least computationally expensive method to build synthetic training datasets. With our initial test dataset, we have achieved accuracies of > 90% and found that the synthetic training dataset produced using the Hapke mixture model provides the best results.

99 GENERAL AND MISCELLANEOUS↗

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗

Entanglement between two disjoint universes

We use the replica method to compute the entanglement entropy of a universe without gravity entangled in a thermofield-double-like state with a disjoint gravitating universe. Including wormholes between replicas of the latter gives an entropy functional which includes an “island” on the gravitating universe. We solve the back-reaction equations when the cosmological constant is negative to show that this island coincides with a causal shadow region that is created by the entanglement in the gravitating geometry. At high entanglement temperatures, the island contribution to the entropy functional leads to a bound on entanglement entropy, analogous to the Page behavior of evaporating black holes. We demonstrate that the entanglement wedge of the non-gravitating universe grows with the entanglement temperature until, eventually, the gravitating universe can be entirely reconstructed from the non-gravitating one.

2D gravity↗

Sum rules in the standard model effective field theory from helicity amplitudes

The dispersion relation of an elastic 4-point amplitude in the forward direction leads to a sum rule that connects the low energy amplitude to the high energy observables. We perform a classification of these sum rules based on massless helicity amplitudes. With this classification, we are able to systematically write down the sum rules for the dimension-6 operators of the Standard Model Effective Field Theory (SMEFT), some of which are absent in previous literatures. These sum rules offer distinct insights on the relations between the operator coefficients in the EFT and the properties of the full theory that generates them. Their applicability goes beyond tree level, and in some cases can be used as a practical method of computing the one loop contributions to low energy observables. They also provide an interesting perspective for understanding the custodial symmetries of the SM Higgs and fermion sectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Repulsive black holes and higher-derivatives

In two-derivative theories of gravity coupled to matter, charged black holes are self-attractive at large distances, with the force vanishing at zero temperature. However, in the presence of massless scalar fields and four-derivative corrections, zero-temperature black holes no longer need to obey the no-force condition. In this paper, we show how to calculate the long-range force between such black holes. We develop an efficient method for computing the higher-derivative corrections to the scalar charges when the theory has a shift symmetry, and compute the resulting force in a variety of examples. We find that higher-derivative corrected black holes may be self-attractive or self-repulsive, depending on the value of the Wilson coefficients and the VEVs of scalar moduli. Indeed, we find black hole solutions which are both superextremal and self-attractive. Furthermore, we present examples where no choice of higher-derivative coefficients allows for self-repulsive black hole states in all directions in charge space. This suggests that, unlike the Weak Gravity Conjecture, which may be satisfied by the black hole spectrum alone, the Repulsive Force Conjecture requires additional constraints on the spectrum of charged particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗