Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

PySAGES: flexible, advanced sampling methods accelerated with GPUs

Abstract Molecular simulations are an important tool for research in physics, chemistry, and biology. The capabilities of simulations can be greatly expanded by providing access to advanced sampling methods and techniques that permit calculation of the relevant underlying free energy landscapes. In this sense, software that can be seamlessly adapted to a broad range of complex systems is essential. Building on past efforts to provide open-source community-supported software for advanced sampling, we introduce PySAGES, a Python implementation of the Software Suite for Advanced General Ensemble Simulations (SSAGES) that provides full GPU support for massively parallel applications of enhanced sampling methods such as adaptive biasing forces, harmonic bias, or forward flux sampling in the context of molecular dynamics simulations. By providing an intuitive interface that facilitates the management of a system’s configuration, the inclusion of new collective variables, and the implementation of sophisticated free energy-based sampling methods, the PySAGES library serves as a general platform for the development and implementation of emerging simulation techniques. The capabilities, core features, and computational performance of this tool are demonstrated with clear and concise examples pertaining to different classes of molecular systems. We anticipate that PySAGES will provide the scientific community with a robust and easily accessible platform to accelerate simulations, improve sampling, and enable facile estimation of free energies for a wide range of materials and processes.

Chemistry↗
DOE OSTI2011

Two-color interferometry and shadowgraphy characterization of an orthogonal double-pulse laser ablation

Orthogonal double pulse laser ablation on aluminum target at atmospheric pressure was studied by time-resolved optical emission spectroscopy, shadowgraphy and two-color interferometry. Studies were realized as a function of the distance from the first pulse plasma induced in air to the sample surface. The evolution of the electron and neutral particles densities were determined by two-color interferometry. Line intensification mechanisms are discussed based on the interferometry results and the plumes interaction dynamics. It was found that the sample heating by the preablation pulse could contribute to the observed emission enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Prior data assisted compressed sensing: A novel MR imaging strategy for real time tracking of lung tumors

Purpose: Hybrid radiotherapy-MRI devices promise real time tracking of moving tumors to focus the radiation portals to the tumor during irradiation. This approach will benefit from the increased temporal resolution of MRI's data acquisition and reconstruction. In this work, the authors propose a novel spatial-temporal compressed sensing (CS) imaging strategy for the real time MRI–-prior data assisted compressed sensing (PDACS), which aims to improve the image quality of the conventional CS without significantly increasing reconstruction times. Methods: Conventional 2D CS requires a random sampling of partial k-space data, as well as an iterative reconstruction that simultaneously enforces the image's sparsity in a transform domain as well as maintains the fidelity to the acquired k-space. PDACS method requires the additional acquisition of the prior data, and for reconstruction, it additionally enforces fidelity to the prior k-space domain similar to viewsharing. In this work, the authors evaluated the proposed PDACS method by comparing its results to those obtained from the 2D CS and viewsharing methods when performed individually. All three methods are used to reconstruct images from lung cancer patients whose tumors move and who are likely to benefit from lung tumor tracking. The patients are scanned, using a 3T MRI, under free breathing using the fully sampled k-space with 2D dynamic bSSFP sequence in a sagittal plane containing lung tumor. These images form a reference set for the evaluation of the partial k-space methods. To create partial k-space, the fully sampled k-space is retrospectively undersampled to obtain a range of acquisition acceleration factors, and reconstructed with 2D-CS, PDACS, and viewshare methods. For evaluation, metrics assessing global image artifacts as well as tumor contour shape fidelity are determined from the reconstructed images. These analyses are performed both for the original 3T images and those at a simulated 0.5T equivalent noise level. Results: In the 3.0T images, the PDACS strategy is shown to give superior results compared to viewshare and conventional 2D CS using all metrics. The 2D-CS tends to perform better than viewshare at the low acceleration factors, while the opposite is true at the high acceleration factors. At simulated 0.5T images, PDACS method performs only marginally better than the viewsharing method, both of which are superior compared to 2D CS. The PDACS image reconstruction time (0.3 s/image) is similar to that of the conventional 2D CS. Conclusions: The PDACS method can potentially improve the real time tracking of moving tumors by significantly increasing MRI's data acquisition speeds. In 3T images, the PDACS method does provide a benefit over the other two methods in terms of both the overall image quality and the ability to accurately and automatically contour the tumor shape. MRI's data acquisition may be accelerated using the simpler viewsharing strategy at the lower, 0.5T magnetic field, as the marginal benefit of the PDACS method may not justify its additional reconstruction times.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

An algorithm to correct saturated mass spectrometry ion abundances for enhanced quantitation and mass accuracy in omic studies

The mass accuracy and peak intensity of ions detected by mass spectrometry (MS) measurements are essential to facilitate compound identification and quantitation. However, high concentration species can easily cause problems if their ion intensities reach beyond the limits of the detection system, leading to distorted and non-ideal detector response (e.g. saturation), and largely precluding the calculation of accurate m/z and intensity values. Here we present an open source computational method to correct peaks above a defined intensity (saturated) threshold determined by the MS instrumentation such as the analog-to-digital converters or time-to-digital converters used in conjunction with time-of-flight MS. In this method, the isotopic envelope for each observed ion above the saturation threshold is compared to its expected theoretical isotopic distribution. The most intense isotopic peak for which saturation does not occur is then utilized to re-calculate the precursor m/z and correct the intensity, resulting in both higher mass accuracy and greater dynamic range. The benefits of this approach were evaluated with proteomic and lipidomic datasets of varying complexities. After correcting the high concentration species, reduced mass errors and enhanced dynamic range were observed for both simple and complex omic samples. Specifically, the mass error dropped by more than 50% in most cases with highly saturated species and dynamic range increased by 1-2 orders of magnitude for peptides in a blood serum sample.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

An instrument for in situ time-resolved X-ray imaging and diffraction of laser powder bed fusion additive manufacturing processes

In situ X-ray-based measurements of the laser powder bed fusion (LPBF) additive manufacturing process produce unique data for model validation and improved process understanding. Synchrotron X-ray imaging and diffraction provide high resolution, bulk sensitive information with sufficient sampling rates to probe melt pool dynamics as well as phase and microstructure evolution. Here, we describe a laboratory-scale LPBF test bed designed to accommodate diffraction and imaging experiments at a synchrotron X-ray source during LPBF operation. We also present experimental results using Ti-6Al-4V, a widely used aerospace alloy, as a model system. Both imaging and diffraction experiments were carried out at the Stanford Synchrotron Radiation Lightsource. Melt pool dynamics were imaged at frame rates up to 4 kHz with a ~1.1 μm effective pixel size and revealed the formation of keyhole pores along the melt track due to vapor recoil forces. Diffraction experiments at sampling rates of 1 kHz captured phase evolution and lattice contraction during the rapid cooling present in LPBF within a ~50 × 100 μm area. In conclusion, we also discuss the utility of these measurements for model validation and process improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of cation distribution on the magnetic and hyperfine behaviour of nanocrystalline Co doped Ni–Zn ferrite (Ni{sub 0.4}Zn{sub 0.4}Co{sub 0.2}Fe{sub 2}O{sub 4})

Highlights: • Cation distribution in tetrahedral and octahedral sites of spinel Ni{sub 0.4}Zn{sub 0.4}Co{sub 0.2}Fe{sub 2}O{sub 4}. • Structural analysis of observed X-ray diffraction pattern using Rietveld method. • Study of hyperfine behaviour using Mössbauer spectroscopy. • Static and dynamic magnetic measurements. • Correlation of cation distributions obtained from Rietveld analysis with the results of magnetic and Mössbauer effect measurements. - Abstract: Nanoparticles of Ni{sub 0.4}Zn{sub 0.4}Co{sub 0.2}Fe{sub 2}O{sub 4} are prepared by a simple co-precipitation method. The as dried sample is heat treated at 400, 500, 600, 700 and 800 °C to obtain different sizes of nanoparticles. The crystallographic phase of the samples is confirmed analyzing observed X-ray diffraction (XRD) by Rietveld method. Hyperfine parameters of the samples are derived from room temperature (RT) Mössbauer spectra of the samples. Magnetic properties of the samples are investigated by static and dynamic hysteresis loops. Different magneto-crystalline parameters are calculated from the variation of magnetization with temperature (M–T curve) under zero field cooled (ZFC) and field cooled (FC) conditions of the as dried sample. The cation distribution estimated from Rietveld analysis are correlated with the results of magnetic and Mössbauer effect measurements. The observed high value of saturation magnetization (72.7 emu/g at RT) of the sample annealed at 800 °C would be interesting for applications in different electromagnetic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

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

In Situ Chemical Monitoring and Imaging of Contents within Microfluidic Devices Having a Porous Membrane Wall Using Liquid Microjunction Surface Sampling Probe Mass Spectrometry

The ability to observe dynamic chemical processes (e.g., signaling, transport, etc.) in vivo or in situ using nondestructive chemical imaging opens a new door to understanding the complex dynamics of developing biological systems. With the advent of “biology-on-a-chip” devices has come the ability to monitor dynamic chemical processes in a controlled environment, using these engineered habitats to capture key features of natural systems while allowing visual observation of system development. Having the capability to spatially and temporally map the chemical signals within these devices may yield new insights into the forces that drive biosystem development. Here, a porous membrane sealed microfluidic device was designed to allow normal microfluidic operation while enabling continuous, location specific sampling and chemical characterization by liquid microjunction surface sampling probe mass spectrometry (LMJ-SSP MS). LMJ-SSP was used to extract fluids with nL-to-μL/min flow rates directly from selected areas of the microfluidic device without negatively impacting the device function. These extracts were subsequently characterized using MS. This technique was used to acquire MS images of the entirety of several multi-input microfluidic devices having different degrees of fluid mixing. LMJ-SSP MS imaging visualized the spatial distribution of chemical components within the microfluidic channels and could visualize chemical reactions occurring in the device. These microfluidic devices with a porous membrane wall are wholly compatible with the construction of biology-on-a-chip devices. This ultimately would enable correlation of biosystem physical structure with an evolving chemical environment

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Pairing a Global Optimization Algorithm with EXAFS to Characterize Lanthanide Structure in Solution

Ensemble-average sampling of structures from ab initio molecular dynamics (AIMD) simulations can be used to predict theoretical extended X-ray absorption fine structure (EXAFS) signals that closely match experimental spectra. However, AIMD simulations are time-consuming and resource-intensive, particularly for solvated lanthanide ions, which often form multiple nonrigid geometries with high coordination numbers. Here, to accelerate the characterization of lanthanide structures in solution, we employed the Northwest Potential Energy Surface Search Engine (NWPEsSe), an adaptive-learning global optimization algorithm, to efficiently screen first-shell structures. As case studies, we examine two systems: Eu(NO 3 ) 3 dissolved in acetonitrile with a terpyridine ligand (terpyNO 2 ), and Nd(NO 3 ) 3 dissolved in acetonitrile. The theoretical spectra for structures identified by NWPEsSe were compared to both experimental and AIMD-derived EXAFS spectra. The NWPEsSe algorithm successfully identified the proper solvation structure for both Eu(NO 3 ) 3 (terpyNO 2 ) and Nd(NO 3 )(acetonitrile) 3 , with the calculated EXAFS signals closely matching the experimental spectra for the Eu-ligand complex and showing good similarity for the Nd salt; the better agreement with the ligand-containing structure is attributed to a less dynamic coordination environment due to the rigid ligand. The key advantage of the global optimization algorithm lies in its ability to sample the coordination environment across the potential energy surface and reduce the time required to identify structures from generally a month to within a week. Additionally, this approach is versatile and can be adapted to characterize main-group metal complexes.

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

Structured illumination with thermal imaging (SI-TI): A dynamically reconfigurable metrology for parallelized thermal transport characterization

The recent push for the “materials by design” paradigm requires synergistic integration of scalable computation, synthesis, and characterization. Among these, techniques for efficient measurement of thermal transport can be a bottleneck limiting the experimental database size, especially for diverse materials with a range of roughness, porosity, and anisotropy. Traditional contact thermal measurements have challenges with throughput and the lack of spatially resolvable property mapping, while non-contact pump-probe laser methods generally need mirror smooth sample surfaces and also require serial raster scanning to achieve property mapping. Here, we present structured illumination with thermal imaging (SI-TI), a new thermal characterization tool based on parallelized all-optical heating and thermometry. Experiments on representative dense and porous bulk materials as well as a 3D printed thermoelectric thick film (~50 μm) demonstrate that SI-TI (1) enables paralleled measurement of multiple regions and samples without raster scanning; (2) can dynamically adjust the heating pattern purely in software, to optimize the measurement sensitivity in different directions for anisotropic materials; and (3) can tolerate rough (~3 μm) and scratched sample surfaces. Here, this work highlights a new avenue in adaptivity and throughput for thermal characterization of diverse materials.

42 ENGINEERING↗
DOE OSTI2022

Improving Seasonal Forecast Using Probabilistic Deep Learning

The path toward realizing the potential of seasonal forecasting and its socioeconomic benefits relies on improving general circulation model (GCM) based dynamical forecast systems. To improve dynamical seasonal forecasts, it is crucial to set up forecast benchmarks, and clarify forecast limitations posed by model initialization errors, formulation deficiencies, and internal climate variability. With huge costs in generating large forecast ensembles, and limited observations for forecast verification, the seasonal forecast benchmarking and diagnosing task proves challenging. Here, we develop a probabilistic deep learning-based statistical forecast methodology, drawing on a wealth of climate simulations to enhance seasonal forecast capability and forecast diagnosis. By explicitly modeling the internal climate variability and GCM formulation differences, the proposed Conditional Generative Forecasting (CGF) methodology enables bypassing crucial barriers in dynamical forecast, and offers a top-down viewpoint to examine how complicated GCMs encode the seasonal predictability information. We apply the CGF methodology for global seasonal forecast of precipitation and 2 m air temperature, based on a unique data set consisting 52,201 years of climate simulation. Results show that the CGF methodology can faithfully represent the seasonal predictability information encoded in GCMs. We successfully apply this learned relationship in real-world seasonal forecast, achieving competitive performance compared to dynamical forecasts. Using this CGF as benchmark, we reveal the impact of insufficient forecast spread sampling that limits the skill of the considered dynamical forecast system. Finally, we introduce different strategies for composing ensembles using the CGF methodology, highlighting the potential for leveraging the strengths of multiple GCMs to achieve advantgeous seasonal forecast.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Project Dinos I: A joint lensing–dynamics constraint on the deviation from the power law in the mass profile of massive ellipticals

The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Adaptive Dynamic Event Tree in RAVEN code

RAVEN is a software tool that is focused on performing statistical analysis of stochastic dynamic systems. RAVEN has been designed in a high modular and pluggable way in order to enable easy integration of different programming languages (i.e., C++, Python) and coupling with other applications (system codes). Among the several capabilities currently present in RAVEN, there are five different sampling strategies: Monte Carlo, Latin Hyper Cube, Grid, Adaptive and Dynamic Event Tree (DET) sampling methodologies. The scope of this paper is to present a new sampling approach, currently under definition and implementation: an evolution of the DET me

Dynamic Event Tree↗
DOE OSTI2015

X-ray scaling relations of “core” and “coreless” E and S0 galaxies

We have re-examined the two X-ray scaling relations of early-type galaxies (ETGs), L{sub X,GAS}–L{sub K} and L{sub X,GAS}–T{sub GAS}, using 61 ATLAS{sup 3D} E and S0 galaxies observed with Chandra (including ROSAT results for a few X-ray bright galaxies with extended hot gas). With this sample, which doubles the number of ETGs available for study we confirm the strong, steep correlations reported by Boroson et al. Moreover, the larger sample allows us to investigate the effect of structural and dynamical properties of ETGs in these relations. Using the sub-sample of 11 “genuine” E galaxies with central surface brightness cores, slow stellar rotations and old stellar populations, we find that the scatter of the correlations is strongly reduced, yielding an extremely tight relation of the form L{sub X,GAS} ∼ T{sub GAS}{sup 4.5±0.3}. The rms deviation is only 0.13 dex. For the gas-rich galaxies in this sample (L{sub X,GAS} > 10{sup 40} erg s{sup −1}), this relation is consistent with recent simulations for velocity dispersion supported E galaxies. However, the tight L{sub X,GAS}–T{sub GAS} relation of genuine E galaxies extends down into the L{sub X} ∼ 10{sup 38} erg s{sup −1} range, where simulations predict the gas to be in outflow/wind state, with resulting L{sub X} much lower than observed in our sample. The observed correlation may suggest the presence of small bound hot halos even in this low luminosity range. At the high luminosity end, the L{sub X,GAS}–T{sub GAS} correlation of core elliptical galaxies is similar to that found in samples of cD galaxies and groups, but shifted down toward relatively lower L{sub X,GAS} for a given T{sub GAS}. In particular cDs, central dominant galaxies sitting at the bottom of the potential well imposed by the group dark matter, have an order of magnitude higher L{sub X,GAS} than our sample core galaxies for the same L{sub K} and T{sub GAS}. We suggest that enhanced cooling in cDs, which have higher hot gas densities and lower entropies, could lower T{sub GAS} to the range observed in giant Es; this conclusion is supported by the presence of extended cold gas in several cDs. Instead, in the sub-sample of coreless ETGs—these galaxies also tend to show stellar rotation, a flattened galaxy figure and rejuvenation of the stellar population—L{sub X,GAS} and T{sub GAS} are not correlated. The L{sub X}−T{sub GAS} distribution of coreless ETGs is a scatter diagram clustered at L{sub X,GAS} < 10{sup 40} erg s{sup −1}, similar to that reported for the hot interstellar medium (ISM) of spiral galaxies, suggesting that both the energy input from star formation and the effect of galactic rotation and flattening may disrupt the hot ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Grain size dependence of dynamic mechanical behavior of AZ31B magnesium alloy sheet under compressive shock loading

The effects of grain size on the dynamic deformation behavior of rolled AZ31B alloy at high strain rates were investigated. Rolled AZ31B alloy samples with grain sizes of 6, 18 and 37 μm, were subjected to shock loading tests using Split Hopkinson Pressure Bar at room temperature and at a strain rate of 1100 s{sup −} {sup 1}. It was found that a double-peak basal texture formed in the shock loaded samples. The strength and ductility of the alloy under the high strain-rate compressive loading increased with decreasing grain size. However, twinning fraction and strain hardening rate were found to decrease with decreasing grain size. In addition, orientation imaging microscopy showed a higher contribution of double and contraction twins in the deformation process of the coarse-grained samples. Using transmission electron microscopy, pyramidal dislocations were detected in the shock loaded sample, proving the activation of pyramidal slip system under dynamic impact loading. - Highlights: • A double-peak basal texture developed in all shock loaded samples. • Both strength and ductility increased with decreasing grain size. • Twinning fraction and strain hardening rate decreased with decreasing grain size. • ‘g.b’ analysis confirmed the presence of dislocations in shock loaded alloy.

36 MATERIALS SCIENCE↗