Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

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 415 records · Page 23

Multi-objective Bayesian active learning for MeV-ultrafast electron diffraction

Ultrafast electron diffraction using MeV energy beams(MeV-UED) has enabled unprecedented scientific opportunities in the study of ultrafast structural dynamics in a variety of gas, liquid and solid state systems. Broad scientific applications usually pose different requirements for electron probe properties. Due to the complex, nonlinear and correlated nature of accelerator systems, electron beam property optimization is a time-taking process and often relies on extensive hand-tuning by experienced human operators. Algorithm based efficient online tuning strategies are highly desired. Here, we demonstrate multi-objective Bayesian active learning for speeding up online beam tuning at the SLAC MeV-UED facility. The multi-objective Bayesian optimization algorithm was used for efficiently searching the parameter space and mapping out the Pareto Fronts which give the trade-offs between key beam properties. Such scheme enables an unprecedented overview of the global behavior of the experimental system and takes a significantly smaller number of measurements compared with traditional methods such as a grid scan. This methodology can be applied in other experimental scenarios that require simultaneously optimizing multiple objectives by explorations in high dimensional, nonlinear and correlated systems.

43 PARTICLE ACCELERATORS↗

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗

Building workflows for an interactive human-in-the-loop automated experiment (hAE) in STEM-EELS

Exploring the structural, chemical, and physical properties of matter on the nano- and atomic scales has become possible with the recent advances in aberration-corrected electron energy-loss spectroscopy (EELS) in scanning transmission electron microscopy (STEM). However, the current paradigm of STEM-EELS relies on the classical rectangular grid sampling, in which all surface regions are assumed to be of equal a priori interest. However, this is typically not the case for real-world scenarios, where phenomena of interest are concentrated in a small number of spatial locations, such as interfaces, structural and topological defects, and multi-phase inclusions. One of the foundational problems is the discovery of nanometer- or atomic-scale structures having specific signatures in EELS spectra. Herein, we systematically explore the hyperparameters controlling deep kernel learning (DKL) discovery workflows for STEM-EELS and identify the role of the local structural descriptors and acquisition functions in experiment progression. In agreement with the actual experiment, we observe that for certain parameter combinations the experiment path can be trapped in the local minima. We demonstrate the approaches for monitoring the automated experiment in the real and feature space of the system and knowledge acquisition of the DKL model. Based on these, we construct intervention strategies defining the human-in-the-loop automated experiment (hAE). This approach can be further extended to other techniques including 4D STEM and other forms of spectroscopic imaging. The hAE library is available on Github at https://github.com/utkarshp1161/hAE/tree/main/hAE.

Pratiush, Utkarsh [Univ. of Tennessee, Knoxville, ↗

Machine Learning of Performance Space Mapping for the DPD Simulation of Drug Delivery to Endothelial Cells

Despite huge effort over the years, the design of functionalised nanocarriers (NCs) for targeted drug delivery to endothelial cells is still to be completely unveiled. Dissipative particle dynamics (DPD) simulation is used to study the adhesion of NCs to endothelial cells under the influence of series of parameters such as shape, size, and ligand density of the NCs. However, preparing a performance space mapping that illustrates the penetration depths of NCs as a function of variations in their properties requires simulations of all possible NCs with the above-mentioned properties, which is not feasible. This challenge was addressed by leveraging a Gaussian process regression (GPR)-informed active learning strategy and an extensive exploration of numerous samples, each representing different properties of NCs. The performance space mapping reveals that NCs with rod and disc shapes exhibit superior penetration capabilities compared to those with a spherical shape. Furthermore, it demonstrates that smaller-sized rod-shaped NCs and larger-sized disc-shaped NCs tend to achieve better penetration. When considering smaller NCs, the influence of ligand density appears to be limited. On the contrary, for larger NCs, an increase in ligand density correlates with greater penetration depth, underscoring its substantial role in shaping their penetration capabilities.

36 MATERIALS SCIENCE↗

Accelerating cosmological inference with Gaussian processes and neural networks – an application to LSST Y1 weak lensing and galaxy clustering

ABSTRACT Studying the impact of systematic effects, optimizing survey strategies, assessing tensions between different probes and exploring synergies of different data sets require a large number of simulated likelihood analyses, each of which cost thousands of CPU hours. In this paper, we present a method to accelerate cosmological inference using emulators based on Gaussian process regression and neural networks. We iteratively acquire training samples in regions of high posterior probability which enables accurate emulation of data vectors even in high dimensional parameter spaces. We showcase the performance of our emulator with a simulated 3×2 point analysis of LSST-Y1 with realistic theoretical and systematics modelling. We show that our emulator leads to high-fidelity posterior contours, with an order of magnitude speed-up. Most importantly, the trained emulator can be re-used for extremely fast impact and optimization studies. We demonstrate this feature by studying baryonic physics effects in LSST-Y1 3×2 point analyses where each one of our MCMC runs takes approximately 5 min. This technique enables future cosmological analyses to map out the science return as a function of analysis choices and survey strategy.

Astronomy & Astrophysics↗

Vacuum-assisted extrusion to reduce internal porosity in large-format additive manufacturing

Large-scale 3D printing of polymer composite structures has gained popularity and seen extensive use over the last decade. Much of the research related to improving the mechanical properties of 3D-printed parts has focused on exploring new materials and optimizing print parameters to improve geometric control and minimize voids between printed beads. However, porosity at the microstructural level (within the printed bead) has been much less studied although it is almost universally observed at levels of 4 %-10 % when using fiber reinforced materials. This study introduces a vacuum-assist approach that minimizes internal porosity by removing ambient air from the interstitial space between pellets in the hopper and acts as a negative pressure vent for gases that evolve during the initial stages of single-screw extrusion. Vacuum-assisted extrusion was able to reduce porosity below 2 % across a wide range of processing parameters, moisture content, fiber reinforcements, and printing platforms. Specifically, when printing on a large-format extruder (Strangpresse Model-30), the vacuum-assisted extrusion reduced internal porosity by 35–75 % compared to conventional non-vacuum extrusion, and only pores with length scale > 2 microns are affected. The success of this approach prompted the design of a patent-pending continuous vacuum hopper relevant for large-scale 3D printing on commercial systems.

36 MATERIALS SCIENCE↗

Correlation between the gas-phase metallicity and ionization parameter in extragalactic H II regions

The variations of the metallicity and ionization parameter in H II regions are usually thought to be the dominant factors that produce the variations we see in the observed emission line spectra. There is an increasing amount of evidence that these two quantities are physically correlated, although the exact form of this correlation is debatable in the literature. Simulated emission line spectra from photoionized clouds provide important clues about the physical conditions of H II regions and are frequently used for deriving metallicities and ionization parameters. Through a systematic investigation on the assumptions and methodology used in applying photoionization models, we find that the derived correlation has a strong dependence on the choice of model parameters. On the one hand, models that give consistent predictions over multiple emission-line ratios yield a positive correlation between the metallicity and ionization parameter for the general population of H II regions or star-forming galaxies. On the other hand, models that are inconsistent with the data locus in high-dimensional line ratio space yield discrepant correlations when different subsets of line ratios are used in the derivation. The correlation between the metallicity and ionization parameter has a secondary dependence on the surface density of the star formation rate (SFR), with the higher SFR regions showing a higher ionization parameter but weaker correlations. The existence of the positive correlation contradicts the analytical wind-driven bubble model for H II regions. We explore assumptions in both dynamical models and photoionization models, and conclude that there is a potential bias associated with the geometry. However, this is still insufficient to explain the correlation. Mechanisms that suppress the dynamical influence of stellar winds in realistic H II regions might be the key to solving this puzzle, though more sophisticated combinations of dynamical models and photoionization models to test are required.

79 ASTRONOMY AND ASTROPHYSICS↗

Emulating galaxy and peculiar velocity clustering on non-linear scales

We explore the potential of cross-correlating galaxies and peculiar velocities on non-linear scales to enhance cosmological constraints. Leveraging the ABACUSSUMMIT simulation suite and the halo occupation distribution (HOD) formalism, we trained emulator models to describe the non-linear clustering of galaxies and velocities in redshift space. Our analysis demonstrates that combining galaxy and peculiar velocity clustering provides tighter constraints on both HOD and cosmological parameters, particularly on σ8 and w0. We further applied our models to realistic mock catalogues, reproducing the expected density and peculiar velocity errors of type-Ia supernovae, Tully-Fisher and fundamental plane measurements for the combined ZTF and DESI measurements. While systematic biases arise in the HOD parameters, the cosmological constraints remain unbiased, yielding a 3.8% precision measurement on fσ 8 compared to 4.7% when using galaxy clustering alone. We demonstrate that while combining tracers with realistic velocity measurements still yields an improvement, the gains are diminished, highlighting the need for further efforts to reduce velocity measurement uncertainties and correct observational systematics on small scales.

79 ASTRONOMY AND ASTROPHYSICS↗

The Simons Observatory: gain, bandpass and polarization-angle calibration requirements for B-mode searches

We quantify the calibration requirements for systematic uncertainties for next-generation ground-based observatories targeting the large-angle B-mode polarization of the Cosmic Microwave Background, with a focus on the Simons Observatory (SO). We explore uncertainties on gain calibration, bandpass center frequencies, and polarization angles, including the frequency variation of the latter across the bandpass. We find that gain calibration and bandpass center frequencies must be known to percent levels or less to avoid biases on the tensor-to-scalar ratio r on the order of Δ r∼10-3, in line with previous findings. Polarization angles must be calibrated to the level of a few tenths of a degree, while their frequency variation between the edges of the band must be known to O(10) degrees. Given the tightness of these calibration requirements, we explore the level to which residual uncertainties on these systematics would affect the final constraints on r if included in the data model and marginalized over. We find that the additional parameter freedom does not degrade the final constraints on r significantly, broadening the error bar by O(10%) at most. We validate these results by reanalyzing the latest publicly available data from the collaboration within an extended parameter space covering both cosmological, foreground and systematic parameters. Finally, our results are discussed in light of the instrument design and calibration studies carried out within SO.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Perspectives and Challenges in Bolide Infrasound Processing and Interpretation: A Focused Review with Case Studies

Infrasound sensing plays a critical role in the detection and analysis of bolides, offering passive, cost-effective global monitoring capabilities. Key objectives include determining the timing, location, and yield of these events. Achieving these goals requires a robust approach to detect, analyze, and interpret rapidly moving elevated sources such as bolides (also re-entry). In light of advancements in infrasonic methodologies, there is a need for a comprehensive overview of the characteristics that distinguish bolides from other infrasound sources and methodologies for bolide infrasound analysis. This paper provides a focused review of key considerations and presents a unified framework to enhance infrasound processing approaches specifically tailored for bolides. Three representative case studies are presented to demonstrate the practical application of infrasound processing methodologies and deriving source parameters while exploring challenges associated with bolide-generated infrasound. These case studies underscore the effectiveness of infrasound in determining source parameters and highlight interpretative challenges, such as variations in signal period measurements across different studies. Future research should place emphasis on improving geolocation and yield accuracy. This can be achieved through rigorous and systematic analyses of large, statistically significant samples of such events, aiming to resolve interpretative inconsistencies and explore the causes for variability in signal periods and back azimuths. The topic described here is also relevant to space exploration involving planetary bodies with atmospheres, such as Venus, Mars, and Titan.

79 ASTRONOMY AND ASTROPHYSICS↗

Properties of first-order hadron-quark phase transition from inverting neutron star observables

By inverting the observational data of several neutron star observables in the three-dimensional parameter space of the constant speed of sound (CSS) model while fixing all hadronic equation-of-state parameters at their currently known most probable values, we constrain the three parameters of the CSS model and their correlations. Using two lower radius limits of $R$ $2.01$ = $11.41$ km and $R$ $2.01$ = $12.2$ km for PSR J0740+6620 obtained from two independent analyses using different approaches by the Neutron Star Interior Composition Explorer (NICER) Collaboration, the speed of sound squared $c$ $^2_{QM}$ in quark matter is found to have a lower limit of 0.35 or 0.43 in units of $c$ $2$ , above its conformal limit of $c$ $^2_{QM}$ < $1/3$. Furthermore, an approximately linear correlation between the first-order hadron-quark transition density $ρ$ $t$ and its strength $Δε$ is found. Moreover, the presence of twin stars is deemed improbable by the present work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-pressure high-temperature synthesis and characterization of B 10 C

The boron-rich boron carbide materials have been traditionally synthesized by adding boron powder to B 4 C material and subjecting it to hot pressing sintering for materials composition containing 8.8–20 at. % carbon in boron (composition range of B 10.4 C to B 4 C). Our study explores a synthesis route for B 10 C starting from high-purity boron and carbon and direct conversion under high pressure and high temperature (HPHT) conditions of 2000 °C and 6–8 GPa. Synthesis was verified via x-ray diffraction analysis, showing the conversion of the high-purity boron and carbon powder mixture into a hexagonal B 10 C structure (R-3m space group) with lattice parameters of a = b = 5.6115 Å and c = 12.197 Å. The concentration of boron was measured through x-ray photoelectron spectroscopy, confirming the B 10 C ratio. The measured nanoindentation mean hardness of B 10 C was 40 GPa. Raman spectroscopy of the HPHT synthesized sample shows characteristic vibrational breathing modes of boron icosahedron and an additional intense band at a vibrational frequency of 380 cm –1 . This Raman band, which appears notably weaker in earlier studies and B 4 C samples, is assigned to the linear chain of B–B–B and attributed to the maximal incorporation of boron within the hexagonal structure.

36 MATERIALS SCIENCE↗

Limitations to the ‘basic’ HOD model and beyond

ABSTRACT We use the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based halo occupation distribution (HOD) approach to modelling the galaxy–halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (basic) HOD model, we find that, on average, the ‘basic’ HOD model underpredicts the real-space correlation function in the TNG300 box by ∼15 per cent on scales of $1 \,\,\lt\,\, r \,\,\lt\,\, 20 \ {\rm Mpc}\, h^{-1}$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ‘secondary’ parameters with which to augment the ‘basic’ HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, β, and the product of the half-mass radius and the velocity dispersion, σ2Rhalfmass, are the three most effective measures of assembly bias that help reconcile the ‘basic’ HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch, and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy–halo connection can be constructed if we combine both mass and local environment information about the halo.

42 ENGINEERING↗

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗

Energy transfer in random-matrix ensembles of Floquet Hamiltonians

Here, we explore the statistical properties of energy transfer in ensembles of doubly driven random-matrix Floquet Hamiltonians based on universal symmetry arguments. The energy-pumping efficiency distribution P⁡($\overline{E}$) is associated with the Hamiltonian parameter ensemble and the eigenvalue statistics of the Floquet operator. For specific Hamiltonian ensembles, P⁡($\overline{E}$) undergoes a transition which cannot be associated with a symmetry breaking of the instantaneous Hamiltonian. The Floquet eigenvalue spacing distribution indicates the considered ensembles constitute generic nonintegrable Hamiltonian families. As a step towards Hamiltonian engineering, we develop a machine-learning classifier to understand the relative parameter importance in resulting high-conversion efficiency. We propose random Floquet Hamiltonians as a general framework to investigate frequency conversion effects in a class of generic dynamical processes beyond adiabatic pumps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

Mapping structural heterogeneity at the nanoscale with scanning nano-structure electron microscopy (SNEM)

Here, in this work, we explore the use of scanning electron diffraction (also known as 4D-STEM) coupled with electron atomic pair distribution function analysis (ePDF) to understand the local order (structure and chemistry) as a function of position in a complex multicomponent system, a hot rolled, Ni-encapsulated, Zr 65 Cu 17.5 Ni 10 Al 7.5 bulk metallic glass (BMG), with a spatial resolution of 3 nm. We show that it is possible to gain insight into the chemistry and chemical clustering/ordering tendency in different regions of the sample, including in the vicinity of nano-scale crystallites that are identified from virtual dark field images and in heavily deformed regions at the edge of the BMG. In addition to simpler analysis, unsupervised machine learning was used to extract partial PDFs from the material, modeled as a quasi-binary alloy, and map them in space. These maps allowed key insights not only into the local average composition, as validated by EELS, but also a unique insight into chemical short-range ordering tendencies in different regions of the sample during formation. The experiments are straightforward and rapid and, unlike spectroscopic measurements, don’t require energy filters on the instrument. We spatially map different quantities of interest (QoI’s), defined as scalars that can be computed directly from positions and widths of ePDF peaks or parameters refined from fits to the patterns. We developed a flexible and rapid data reduction and analysis software framework that allows experimenters to rapidly explore images of the sample on the basis of different QoI’s. The power and flexibility of this approach are explored and described in detail. Because of the fact that we are getting spatially resolved images of the nanoscale structure obtained from ePDFs we call this approach scanning nano-structure electron microscopy (SNEM), and we believe that it will be powerful and useful extension of current 4D-STEM methods.

36 MATERIALS SCIENCE↗