Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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,351 records · Page 75

A pectin methyltransferase modulates polysaccharide dynamics and interactions in Arabidopsis primary cell walls: Evidence from solid-state NMR

Plant cell walls contain cellulose embedded in matrix polysaccharides. Understanding carbohydrate structures and interactions is critical to the production of biofuel and biomaterials using these natural resources. Here we present a solid-state NMR study of cellulose and pectin in 13 C-labeled cell walls of Arabidopsis wild-type and mutant plants. Using 1D 13 C and 2D 13 C– 13 C correlation experiments, we detected a highly branched arabinan structure in qua2 and tsd2 samples, two allelic mutants for a pectin methyltransferase. Both mutants show close physical association between cellulose and the backbones of pectic homogalacturonan and rhamnogalacturonan-I. Relaxation and dipolar order parameters revealed enhanced microsecond dynamics due to polymer disorder in the mutants, but restricted motional amplitudes due to tighter pectin-cellulose associations. These molecular data shed light on polymer structure and packing in these two pectin mutants, helping to elucidate how pectin could influence cell wall architecture at the nanoscale, cell wall mechanics, and plant growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials

Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.

14 SOLAR ENERGY↗

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗

Demonstration of an AI-driven workflow for autonomous high-resolution scanning microscopy

Abstract Modern scanning microscopes can image materials with up to sub-atomic spatial and sub-picosecond time resolutions, but these capabilities come with large volumes of data, which can be difficult to store and analyze. We report the Fast Autonomous Scanning Toolkit (FAST) that addresses this challenge by combining a neural network, route optimization, and efficient hardware controls to enable a self-driving experiment that actively identifies and measures a sparse but representative data subset in lieu of the full dataset. FAST requires no prior information about the sample, is computationally efficient, and uses generic hardware controls with minimal experiment-specific wrapping. We test FAST in simulations and a dark-field X-ray microscopy experiment of a WSe 2 film. Our studies show that a FAST scan of <25% is sufficient to accurately image and analyze the sample. FAST is easy to adapt for any scanning microscope; its broad adoption will empower general multi-level studies of materials evolution with respect to time, temperature, or other parameters.

47 OTHER INSTRUMENTATION↗

Constraining primordial non-Gaussianity from DESI DR1 quasars and Planck PR4 CMB lensing

We present the first measurement of local-type primordial non-Gaussianity from the cross-correlation between 1.2 million spectroscopically confirmed quasars from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI) and the Planck PR4 CMB lensing reconstructions. The analysis is performed in three tomographic redshift bins covering 0.8 < z < 3.5, covering a sky fraction of ∼20%. We adopt a catalog-based pseudo-C ℓ estimator and apply linear imaging weights validated on noiseless mocks. Compared to previous analyses using photometric quasar samples, our results benefit from the high purity of the DESI spectroscopic sample, the reduced noise of PR4 lensing, and the absence of excess large-scale power in the spectroscopic quasar auto-correlation. Fitting simultaneously for the non-Gaussianity parameter f NL and the linear bias amplitude in each redshift bin, we obtain f NL = 2 +28 -34 for a response parameter p = 1.6, and f NL = 6 +20 -24 for p = 1.0. These results improve the constraints on f NL by ∼35% compared to the previous analysis based on the Legacy Imaging Survey DR9. Additionally, we derive an optimal weighting scheme to maximize the constraining power. In this case, and assuming p = 1.6, we obtain f NL = 19 +25 -31 . Our results demonstrate the statistical power of DESI quasars for probing inflationary physics, and highlight the promise of future DESI data releases.

cosmological parameters from CMBR↗

Microwave-assisted synthesis of iron sulfide motifs for electrochemical applications

Abstract The syntheses of FeS 2 and Fe 3 S 4 nanomaterials were optimized using a novel facile, surfactant-free, and microwave-assisted, one-pot synthesis method, run under ambient and reasonably mild reaction conditions. Synthetic parameters, such as metal precursor salt identity, reaction time, reaction temperature, metal:sulfur molar ratios, and solvent combinations, were all systematically investigated and optimized. A series of FeS 2 (pyrite) samples was initially fabricated using thioacetamide (TAA) as the sulfur precursor to generate a distinctive, uniform octahedra-based morphology. Switching the sulfur precursor from TAA to L-cysteine resulted in a corresponding transformation in not only chemical composition from FeS 2 to an iron thiospinel structure, Fe 3 S 4 (otherwise known as greigite), but also an associated morphological evolution from octahedra to nanosheet aggregates. The study of these materials has enabled crucial insights into the formation mechanisms of these materials under a relatively non-conventional microwave-assisted setting. Furthermore, in separate experiments, multi-walled carbon nanotubes (MWNTs) and graphene were added in with underlying metal sulfide species to create conductive Fe–S/MWNT composites and Fe–S/graphene composites, respectively. The method of addition of either MWNTs or graphene was also explored, wherein an ‘ ex-situ ’ synthetic procedure was found to be the least disruptive means of attachment and immobilization onto iron sulfide co-reagents as a means of preserving the latter’s inherent composition and morphology. The redox acidity for the parent material and associated composites demonstrates the utility of our as-developed synthetic methods for creating motifs relevant for electrochemical applications, such as energy storage.

25 ENERGY STORAGE↗

Topological Charge Pumping in Excitonic Insulators

Here, we show that in excitonic insulators with s-wave electron-hole pairing, an applied electric field (either pulsed or static) can induce a p-wave component to the order parameter, and further drive it to rotate in the s + ip plane, realizing a Thouless charge pump. In one dimension, each cycle of rotation pumps exactly two electrons across the sample. Higher dimensional systems can be viewed as a stack of one-dimensional chains in momentum space in which each chain crossing the Fermi surface contributes a channel of charge pumping. Physics beyond the adiabatic limit, including in particular dissipative effects is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarking Optimizers for Qumode State Preparation with Variational Quantum Algorithms

Quantum state preparation involves preparing a target state from an initial system, a process integral to applications such as quantum machine learning and solving systems of linear equations. Recently, there has been a growing interest in qumodes due to advancements in the field and their potential applications. However there is a notable gap in the literature specifically addressing this area. This paper aims to bridge this gap by providing performance benchmarks of various optimizers used in state preparation with Variational Quantum Algorithms. We conducted extensive testing across multiple scenarios, including different target states, both ideal and sampling simulations, and varying numbers of basis gate layers. Our evaluations offer insights into the complexity of learning each type of target state and demonstrate that some optimizers perform better than others in this context. Notably, the Powell optimizer was found to be exceptionally robust against sampling errors, making it a preferred choice in scenarios prone to such inaccuracies. Additionally, the Simultaneous Perturbation Stochastic Approximation optimizer was distinguished for its efficiency and ability to handle increased parameter dimensionality effectively.

Kan, Shuwen [Fordham University]↗

Void-free Copper Electrodeposition in High Aspect Ratio, Full Wafer Thickness Through-Silicon Vias with Endpoint Detection

High density interconnects are required for increased input/output for microelectronics applications, incentivizing the development of Cu electrochemical deposition (ECD) processes for high aspect ratio through-silicon vias (TSVs). This work outlines Cu ECD processes for 62.5 μ m diameter TSVs, etched into a 625 μ m thick silicon substrate, a 10:1 aspect ratio. Cu ECD in high aspect ratio features relies on a delicate balance of electrolyte composition, solution replenishment, and applied voltage. Implementing a CuSO 4 -H 2 SO 4 electrolyte, which contains suppressor and a low chloride concentration, allows for a tunable relationship between applied voltage and localized deposition in the vias. A stepped potential waveform was applied to move the Cu growth front from the bottom of the via to the top. Sample characterization was performed through mechanical cross-sections and X-ray computed tomography (CT) scans. The CT scans revealed small seam voids in the Cu electrodeposit, and process parameters were tuned accordingly to produce void-free Cu features. During the voltage-controlled experiments, measured current data showed a characteristic current minimum, which was identified as an endpoint detection method for Cu deposition in these vias. We believe this is the first report of this novel endpoint detection method for TSV filling.

42 ENGINEERING↗

RxnRover/amlro

AMLRO (Active Machine Learning Reaction Optimizer) is an open-source framework designed to accelerate chemical reaction optimization using active learning with classical machine learning regression models. AMLRO integrates space-filling sampling strategies (e.g., Sobol and Latin Hypercube sampling) with iterative model training, prediction, and experiment selection to efficiently navigate complex reaction spaces. The platform supports multiple regression models, flexible multi-objective definitions, and user-defined parameter bounds, enabling data-efficient optimization from small initial datasets. AMLRO is designed for ease of use by experimentalists and can operate as a standalone decision-support tool or be integrated into closed-loop automated experimentation workflows.

Kulathunga, Dulitha Prasanna [Iowa State Universit↗

Cloud microphysical structure analysis based on high-resolution insitu measurements

This study investigates the microphysical parameters and shapes of droplet size distributions (DSDs) along three aircraft traverses of developing convective clouds during Cloud Aerosol Interactions and Precipitation Enhancement Experiment (CAIPEEX) Phase III 2015 at a sampling frequency of 25 Hz. The droplet number concentration (Nc, cm -3 ), and liquid water content (LWC, g m -3 ) present steep gradients within zones of a few tens of meters near the cloud edges and relatively gentle gradients in the strong updraft zones. Sometimes, the horizontal LWC distribution resembles a trapezoidal shape with steep LWC and Nc gradients near the cloud edges. The LWC maximums (LWC max ) are lower than the adiabatic LWCs, but the high adiabatic fractions in the cloud core indicate low dilution. High LWC/LWC max , largest droplets, and narrow and similarly shaped DSDs are found in the regions of high updrafts. Zones of low LWC/LWC max are found close to the cloud edges, where DSDs are highly diverse, containing both large and small droplets. Finally, we analyze the mixing diagrams. Significant in-phase turbulent fluctuations in LWC and N c were found. The effective radii change slightly across cloud updraft zones but decrease at the zone of low LWC/LWC max ratio close to cloud edges. The spectra of LWC and N c obey Kolmogorov -5/3 turbulence law. The radii of the correlation of LWC and N c in updraft zones are of several tens of meters. Filaments up to 120–175 m in size are also noticed.

54 ENVIRONMENTAL SCIENCES↗

Interferometric focal planes

We propose arrays of integrated interferometers to characterize the mutual intensity on focal planes. While focal coherence measurement does not increase the aperture-limited spatial bandpass, it can increase Shannon information capacity relative to irradiance measurement by increasing the number of degrees of freedom per Nyquist sample. We describe a sampling model for interferometric measurement using arrays of 2-port Mach-Zehnder interferometers and show within this model that interferometric focal planes enable more accurate estimation of prototypical scene parameters.

Brady, David J. (ORCID:0000000156552478)↗

Prospective Study Reveals Host Microbial Determinants of Clinical Response to Fecal Microbiota Transplant Therapy in Type 2 Diabetes Patients

Increasing evidence shows that alterations in gut microbiome (GM) contribute to the development of type 2 diabetes mellitus (T2DM), and fecal microbiota transplantation (FMT) successfully treats various human diseases. However, the benefits of FMT therapy to T2DM patients remain unknown. We enrolled 17 patients with T2DM for nonblinded, one-armed intervention trial of FMT. A total of 20 healthy individuals were recruited as the baseline control. HbA1c% and metabolic parameter change were evaluated in 17 T2DM patients 12 weeks after they received FMT from healthy donors. The GM composition was characterized by 16S rRNA gene amplicon sequencing from fecal samples prior to and 12 weeks after FMT treatment. We found that the GM of T2DM patients was reconstituted by FMT. We observed a statistically significant decrease in HbA1c% (from 7.565 ± 0.148 to 7.190 ± 0.210, p<0.01), blood glucose (from 8.483 ± 0.497 to 7.286 ± 0.454 mmol/L, p<0.01), and uric acid (from 309.4 ± 21.5 to 259.1 ± 15.8 µmol/L, p<0.01) while a significant increase in postprandial C-peptide (from 4.503 ± 0.600 to 5.471 ± 0.728 ng/ml, p<0.01) at 12 weeks after FMT. Closely evaluating the changes in these assays, we found individual variability in response to FMT treatment. Out of 17 T2DM patients, 11 were found to significantly improve T2DM symptoms. The FMT responders have significantly higher levels of the family Rikenellaceae and the genus Anaerotruncus (family Ruminococcaceae) in their pretreated fecal in comparison to nonresponders, which could predict the clinical response with an area under the curve of 0.83. Our findings suggest that certain T2DM patients can potentially benefit from FMT, and the pretreated abundance of Rikenellaceae and Anaerotruncus in the fecal of patients may serve as potential biomarkers for selecting T2DM patients to receive FMT.

59 BASIC BIOLOGICAL SCIENCES↗

The Pantheon+ Analysis: Cosmological Constraints

We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from z = 0.001 to 2.26. This work features an increased sample size from the addition of multiple cross-calibrated photometric systems of SNe covering an increased redshift span, and improved treatments of systematic uncertainties in comparison to the original Pantheon analysis, which together result in a factor of 2 improvement in cosmological constraining power. For a flat ΛCDM model, we find Ω M = 0.334 ± 0.018 from SNe Ia alone. For a flat w 0 CDM model, we measure w 0 = –0.90 ± 0.14 from SNe Ia alone, H 0 = 73.5 ± 1.1 km s –1 Mpc –1 when including the Cepheid host distances and covariance (SH0ES), and w 0 = $-{0.978}_{-0.031}^{+0.024}$ when combining the SN likelihood with Planck constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both w 0 values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a flat w 0 w a CDM universe, and measure w a = $-{0.1}_{-2.0}^{+0.9}$ from Pantheon+ SNe Ia alone, H 0 = 73.3 ± 1.1 km s –1 Mpc –1 when including SH0ES Cepheid distances, and w a = $-{0.65}_{-0.32}^{+0.28}$ when combining Pantheon+ SNe Ia with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one-third of the total uncertainty in the measurement of H 0 and cannot explain the present "Hubble tension" between local measurements and early universe predictions from the cosmological model.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Compaction of crushed salt for safe containment – overview of the KOMPASS project

Abstract. In Germany, rock salt formations are possible host rock candidates for a repository for heat-emitting radioactive waste. The safety concept of a repository in salt bases on a multibarrier system consisting mainly of the geological barrier salt and geotechnical seals ensuring safe containment. Crushed salt will be used for backfilling of cavities and sealing measures in drifts and shafts due to its favourable properties and its easy availability (mined-off material). The creep of the rock salt leads to crushed salt compaction with time. Thereby, the crushed salts' porosity is reduced from the initial porosity of 30 %–40 % to a value comparable to the porosity of undisturbed rock salt (≤1 %). In such low porosity ranges, technical impermeability is assumed. The compaction behaviour of crushed salt is rather complex and involves several coupled THM processes (Kröhn et al., 2017; Hansen et al., 2014). It is influenced by internal properties like humidity and grain size distribution, as well as boundary conditions such as temperature, compaction rate or stress state. However, the current process understanding has some important gaps referring to the material behaviour, experimental database and numerical modelling. It needs to be extended and validated, especially in the low porosity range. The objective of the KOMPASS project was development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. Key results are as follows (KOMPASS Phase 1, 2020): selection of an easily available and permanently producible synthetic crushed salt mixture, acting as a reference material for generic investigations; development and proof of different techniques for producing pre-compacted samples for further investigations; establishment of a tool of microstructure investigation methods to demonstrate the comparability of grain structures of pre-compacted samples with in-situ compacted material for future investigations; execution of various laboratory experiments using pre-compacted samples, e.g. long-term creep tests which deliver reliable information about time- and stress-dependent compaction behaviour; development of a complex experimental investigation strategy to derive necessary model parameters considering individual functional dependencies. Its technical feasibility was successfully verified; benchmarking with various existing numerical models using datasets from three different triaxial long-term tests. The result was not entirely satisfactory; however, the number of influencing factors is small and further validation work has to be done. Overall, the KOMPASS project has made significant progress in the approaches to solving the outstanding question, building the basis for further investigations.

Friedenberg, Larissa↗

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗

The Dark Energy Survey: Cosmology Results with ∼1500 New High-redshift Type Ia Supernovae Using the Full 5 yr Data Set

Abstract We present cosmological constraints from the sample of Type Ia supernovae (SNe Ia) discovered and measured during the full 5 yr of the Dark Energy Survey (DES) SN program. In contrast to most previous cosmological samples, in which SNe are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being an SN Ia, we find 1635 DES SNe in the redshift range 0.10 < z < 1.13 that pass quality selection criteria sufficient to constrain cosmological parameters. This quintuples the number of high-quality z > 0.5 SNe compared to the previous leading compilation of Pantheon+ and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints, we combine the DES SN data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning 0.025 < z < 0.10. Using SN data alone and including systematic uncertainties, we find Ω M = 0.352 ± 0.017 in flat ΛCDM. SN data alone now require acceleration ( q 0 < 0 in ΛCDM) with over 5 σ confidence. We find ( Ω M , w ) = ( 0.264 − 0.096 + 0.074 , − 0.80 − 0.16 + 0.14 ) in flat w CDM. For flat w 0 w a CDM, we find ( Ω M , w 0 , w a ) = ( 0.495 − 0.043 + 0.033 , − 0.36 − 0.30 + 0.36 , − 8.8 − 4.5 + 3.7 ) , consistent with a constant equation of state to within ∼2 σ . Including Planck cosmic microwave background, Sloan Digital Sky Survey baryon acoustic oscillation, and DES 3 × 2pt data gives (Ω M , w ) = (0.321 ± 0.007, −0.941 ± 0.026). In all cases, dark energy is consistent with a cosmological constant to within ∼2 σ . Systematic errors on cosmological parameters are subdominant compared to statistical errors; these results thus pave the way for future photometrically classified SN analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗