Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampled value”

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 127 records · Page 7

A New Method for Deriving Composition of S-type Asteroids from Noisy and Incomplete Near-infrared Spectra

The surface composition of S-type asteroids can be determined using band parameters extracted from their near-infrared (NIR) spectra (0.7–2.50 μm) along with spectral calibrations derived from laboratory samples. In the past, these empirical equations have been obtained by combining NIR spectra of meteorite samples with information about their composition and mineral abundance. For these equations to give accurate results, the characteristics of the laboratory spectra they are derived from should be similar to those of asteroid spectral data (i.e., similar signal-to-noise ratio (S/N) and wavelength range). Here we present new spectral calibrations that can be used to determine the mineral composition of ordinary chondrite-like S-type asteroids. Contrary to previous work, the S/N of the ordinary chondrite spectra used in this study has been decreased to recreate the S/N typically observed among asteroid spectra, allowing us to obtain more realistic results. In addition, the new equations have been derived for five wavelength ranges encompassed between 0.7 and 2.50 μm, making it possible to determine the composition of asteroids with incomplete data. The new spectral calibrations were tested using band parameters measured from the NIR spectrum of asteroid (25143) Itokawa, and comparing the results with laboratory measurements of the returned samples. We found that the spectrally derived olivine and pyroxene chemistry, which are given by the molar contents of fayalite (Fa) and ferrosilite (Fs), are in excellent agreement with the mean values measured from the samples (Fa{sub 28.6±1.1} and Fs{sub 23.1±2.2}), with a maximum difference of 0.6 mol% for Fa and 1.4 mol% for Fs.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating the ecological fallacy through sampling distributions constructed from finite populations

Correlation coefficients and linear regression values computed from group averages can differ from correlation coefficients and linear regression values computed using individual scores. This observation known as the ecological fallacy often assumes that all the individual scores are available from a population. In many situations, one must use a sample from the larger population. In such cases, the computed correlation coefficient and linear regression values will depend on the sample that is chosen and the underlying sampling distribution. The sampling distribution of correlation coefficients and linear regression values for group averages will be identical to the sampling distribution for individuals for normally distributed variables for random samples drawn from infinitely large continuous distributions. However, data that is acquired in practice is often acquired when sampling without replacement from a finite population. Our objective is to demonstrate through Monte Carlo simulations that the sampling distributions for correlation and linear regression will also be similar for individuals and group averages when sampling without replacement from normally distributed variables. These simulations suggest that when a random sample from a population is selected, the correlation coefficients and linear regression values computed from individual scores will not be more accurate in estimating the entire population values compared to samples when group averages are used as long as the sample size is the same.

97 MATHEMATICS AND COMPUTING↗

Whole genome sequencing of Mycobacterium bovis directly from clinical tissue samples without culture

Advancement in next generation sequencing offers the possibility of routine use of whole genome sequencing (WGS) for Mycobacterium bovis (M. bovis) genomes in clinical reference laboratories. To date, the M. bovis genome could only be sequenced if the mycobacteria were cultured from tissue. This requirement for culture has been due to the overwhelmingly large amount of host DNA present when DNA is prepared directly from a granuloma. To overcome this formidable hurdle, we evaluated the usefulness of an RNA-based targeted enrichment method to sequence M. bovis DNA directly from tissue samples without culture. Initial spiking experiments for method development were established by spiking DNA extracted from tissue samples with serially diluted M. bovis BCG DNA at the following concentration range: 0.1 ng/μl to 0.1 pg/μl (10 –1 to 10 –4 ). Library preparation, hybridization and enrichment was performed using SureSelect custom capture library RNA baits and the SureSelect XT HS2 target enrichment system for Illumina paired-end sequencing. The method validation was then assessed using direct WGS of M. bovis DNA extracted from tissue samples from naturally (n = 6) and experimentally (n = 6) infected animals with variable Ct values. Direct WGS of spiked DNA samples achieved 99.1% mean genome coverage (mean depth of coverage: 108×) and 98.8% mean genome coverage (mean depth of coverage: 26.4×) for tissue samples spiked with BCG DNA at 10 –1 (mean Ct value: 20.3) and 10 –2 (mean Ct value: 23.4), respectively. The M. bovis genome from the experimentally and naturally infected tissue samples was successfully sequenced with a mean genome coverage of 99.56% and depth of genome coverage ranging from 9.2× to 72.1×. The spoligoyping and M. bovis group assignment derived from sequencing DNA directly from the infected tissue samples matched that of the cultured isolates from the same sample. Our results show that direct sequencing of M. bovis DNA from tissue samples has the potential to provide accurate sequencing of M. bovis genomes significantly faster than WGS from cultures in research and diagnostic settings.

59 BASIC BIOLOGICAL SCIENCES↗

Revisiting Néel 60 years on: The magnetic anisotropy of L1 0 FeNi (tetrataenite)

The magnetocrystalline anisotropy energy of atomically ordered L1 0 FeNi (the meteoritic mineral tetrataenite) is studied within a first-principles electronic structure framework. Two compositions are examined: equiatomic Fe 0.5 Ni 0.5 and an Fe-rich composition, Fe 0.56 Ni 0.44 . It is confirmed that, for the single crystals modeled in this work, the leading-order anisotropy coefficient K 1 dominates the higher-order coefficients K 2 and K 3 . To enable comparison with experiment, the effects of both imperfect atomic long-range order and finite temperature are included. While our computational results initially appear to undershoot the measured experimental values for this system, careful scrutiny of the original analysis due to Néel et al. [J. Appl. Phys. 35, 873 (1964)] suggests that our computed value of K 1 is, in fact, consistent with experimental values, and that the noted discrepancy has its origins in the nanoscale polycrystalline, multivariant nature of experimental samples, that yields much larger values of K 2 and K 3 than expected a priori. These results provide fresh insight into the existing discrepancies in the literature regarding the value of tetrataenite’s uniaxial magnetocrystalline anisotropy in both natural and synthetic samples.

36 MATERIALS SCIENCE↗

Quantifying uncertainty in uranium concentration measurements via K-edge densitometry

This study quantifies the uncertainty in uranium concentration predictions of fluoride and chloride-based salts within a steel pipe using K-edge densitometry. Modeling and simulation was conducted with the Monte Carlo N-Particle Transport (MCNP) code. The quality of of this technique’s prediction in a pipe requires proper characterization of the pipe’s thickness, which is dependent on the source size and axial offset from the pipe centerline. The thickness was determined as either the center-line thickness seen by the X-ray source or an average value determined through random sampling. Generally, the predicted concentrations were slightly better at lower offset with the random sampling thickness and using the center-line thickness for the highest offsets. For a line-beam source and varying axial offsets, the relative error of concentration was within 1% of the true value but uncertainty increased by 2 orders of magnitude. Similarly, for no axial offset, the relative error was significantly less than 1% while no trend for uncertainty was found. However, at the largest possible offset for a given source size, the concentrations become erroneous and greater than the allowable 1% relative error. Furthermore, high offsets tended to increase the variance of the transmission spectra by 3 orders of magnitude.

Characterization and Analytical Technique↗

Dark Energy Survey Year 3 results: marginalization over redshift distribution uncertainties using ranking of discrete realizations

ABSTRACT Cosmological information from weak lensing surveys is maximized by sorting source galaxies into tomographic redshift subsamples. Any uncertainties on these redshift distributions must be correctly propagated into the cosmological results. We present hyperrank, a new method for marginalizing over redshift distribution uncertainties, using discrete samples from the space of all possible redshift distributions, improving over simple parametrized models. In hyperrank, the set of proposed redshift distributions is ranked according to a small (between one and four) number of summary values, which are then sampled, along with other nuisance parameters and cosmological parameters in the Monte Carlo chain used for inference. This approach can be regarded as a general method for marginalizing over discrete realizations of data vector variation with nuisance parameters, which can consequently be sampled separately from the main parameters of interest, allowing for increased computational efficiency. We focus on the case of weak lensing cosmic shear analyses and demonstrate our method using simulations made for the Dark Energy Survey (DES). We show that the method can correctly and efficiently marginalize over a wide range of models for the redshift distribution uncertainty. Finally, we compare hyperrank to the common mean-shifting method of marginalizing over redshift uncertainty, validating that this simpler model is sufficient for use in the DES Year 3 cosmology results presented in companion papers.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploration for Blind Geothermal Resources in the State of Hawaii Utilizing Dissolved Noble Gasses in Well Waters

This study is an extension of the Hawaii Play Fairway Analysis (PFA), a statewide geothermal exploration project funded by the United States Department of Energy. Based on results from prior phases of the PFA, this project targeted 66 wells on the islands of Hawaii, Maui, Lanai, Oahu, and Kauai for sampling of dissolved noble gases, trace metals, common ions, and the stable isotopes 2H and 18O. Ultimately, 23 of the 66 well targets were sampled. The noble gas helium is used as an indicator of geothermal heat when excess 3He and/or 4He is present when compared to the atmospheric ratio of those isotopes (R/Ra). R/Ra is minimally affected by dilution and transport, allowing even those wells not perfectly situated over a geothermal system to indicate a geothermal anomaly. R/Ra anomalies are present on every island in this study. There is a strong correlation between R/Ra anomalies and proximity to rift zones and calderas. Across the islands R/Ra ranged from 15-16 on Kilauea's lower east rift zone, which is a mantle plume value, to 0.37 on Lanai, which is a crustal value. The majority of anomalous well samples had R/Ra values consistent with an upper mantle source. Mixing between upper mantle and crustal helium is evident on all islands. Geographically, R/Ra decreases from the high at Kilauea to upper mantle values at Mauna Loa, and remains at upper mantle values for all wells across the other volcanoes, with the exception of two of four sampled wells on Lanai.

15 GEOTHERMAL ENERGY↗

Hydrothermal liquefaction of Cyanidioschyzon merolae and Salicornia bigelovii Torr.: The interaction effect on product distribution and chemistry

Hydrothermal liquefaction (HTL) has proven to be a promising technology for processing wet biomass into valuable biofuel intermediates. We investigated the influence of HTL process parameters on the individual and co-liquefaction of Salicornia bigelovii Torr. (SL) and Cyanidioschyzon merolae (CM) biomass and their respective product properties and characteristics. SL and CM biomasses were subjected to HTL conditions of process temperatures (250–350 °C), 10 wt% solid loading, 30 min reaction time, and co-liquefaction ratios (80–20, 50–50, 20–80, CM-SL, respectively). The individual maximum HTL biocrude yields for CM (34.63 wt%) and SL (7.63 wt%) were obtained at 300 °C. In both cases, biocrude oils had a higher carbon and hydrogen content relative to those of the original biomass. The HHV (higher heating value) values for the biocrude samples ranged from 24.2 to 33.12 MJ/kg, which were significantly higher than those of the original biomasses. A positive synergistic effect was observed during the co-liquefaction of two biomasses, where the experimental yield of biocrude oil was higher than the theoretical value. During the co-liquefaction experiments, a maximum biocrude yield of 32.95 wt% was obtained at 300 °C and 80–20 CM-SL biomass mix ratio. Co-liquefaction led to an increase in the energy content of the co-liquefied biocrude oil and a higher energy recovery rate. Finally, the properties of HTL products and biomasses were determined with TGA, CHN elemental analyzer, GC/MS, bomb calorimeter, HACH Spectrophotometer, and ICP OES.

09 BIOMASS FUELS↗

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A review of imputation strategies for isobaric labeling-based shotgun proteomics

The throughput efficiency and increased depth of coverage provided by isobaric-labeled proteomics measurements have led to increased usage of these techniques. However, the structure of missing data is uniquely different than unlabeled studies. In this review, we compare the efficacy of nine imputation methods on a CPTAC proteomics iTRAQ dataset. Imputation methods were evaluated with regard to accuracy, variability, statistical hypothesis test inference and run time over datasets consisting of varying number of iTRAQ plexes and percentages of missing data. In general, expectation maximization and random forest imputation methods yielded the best performances, and constant-based methods performed poorly consistently across all dataset sizes and percentages of missing values. For datasets with small sample sizes and higher percentages of missing data, results indicate that statistical inference with no imputation may be preferable. Based on the findings in this review, there are core imputation methods that perform higher for isobaric-labeled proteomics data, but great care and consideration as to whether imputation should be used should be given for datasets comprised of a small number of samples, as well as to factors such as computational time and reproducibility of imputation values.

Bramer, Lisa M.↗

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗

Considerations for Defining G-Values for Aluminum-Clad Spent Nuclear Fuel

Sealed-canister dry storage of aluminum-clad spent nuclear fuel (ASNF) generated by research reactors is an alternative to current storage and disposition pathways as directed by the U.S. Department of Energy. The major challenge faced for this storage approach is radiolytic H 2 generation, including from the aluminum (oxy)hydroxide layers on the surface of ASNF. Experimental and modeling activities have been carried out to characterize the radiolytic yield as part of a DOE-sponsored research program to develop the technical basis for ASNF dry storage. The G-value is a commonly way to report results of radiolysis testing and is defined as the radiolytic yield of a species (e.g. molecular hydrogen) per unit radiation energy deposited into the material system. An independent technical review of the ASNF dry storage technical basis performed by Pacific Northwest National Laboratory raised questions about differences in G-value definitions used for experiments on ASNF surrogates consisting of aluminum samples with adherent (oxy)hydroxides compared to G-values reported in prior literature and how the magnitudes compared between different studies. Material systems resembling ASNF pose complications for measuring/defining G-values to predict the evolution of H 2 in a sealed canister, including i) accounting for radiolytic yields potentially arising from multiple sources, i.e., residual free (vapor), physisorbed, and chemisorbed/chemically bound waters; ii) deciding what portions of the multi-material system to include in the absorbed energy (radiation dose) calculation, considering possible energy exchange between materials as well as measurement limitations, and iii) capturing variations in G-value associated with non-linear yield vs. dose curves and/or dependence on the cover gas. This report summarizes previous literature information on radiolytic H 2 generation and associated G-values from mixed-material systems (generally oxides in contact with water or organic compounds) and from (oxy)hydroxides/hydrates to compare with the definitions and values for ASNF surrogate samples containing adherent aluminum (oxy)hydroxides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

First direct 7 Be electron-capture $\mathrm{Q}$-value measurement toward high-precision searches for neutrino physics beyond the Standard Model

Here, we report the first direct measurement of the nuclear electron-capture (EC) decay Q value of 7 Be → 7 Li via high-precision Penning trap mass spectrometry (PTMS). This was performed using the LEBIT Penning trap located at the National Superconducting Cyclotron Laboratory/Facility for Rare Isotope Beams (NSCL/FRIB) using the newly commissioned Batch-Mode Ion-Source (BMIS) to deliver the unstable 7 Be + samples. With a measured value of Q EC = 861.963(23) keV, this result is three times more precise than any previous determination of this quantity. This improved precision and accuracy of the 7 Be EC decay Q value is critical for ongoing experiments that measure the recoiling nucleus in this system as a signature to search for beyond the Standard Model (BSM) neutrino physics using 7 Be-doped superconducting sensors. This experiment has extended LEBIT capabilities, using the first low-energy beam delivered by BMIS at FRIB for PTMS, as well as measuring the lightest-mass isotopes so far with LEBIT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predicting Rare Earth Element Potential in Produced and Geothermal Waters of the United States via Emergent Self-Organizing Maps

This work applies emergent self-organizing map (ESOM) techniques, a form of machine learning, in the multidimensional interpretation and prediction of rare earth element (REE) abundance in produced and geothermal waters in the United States. Visualization of the variables in the ESOM trained using the input data shows that each REE, with the exception of Eu, follows the same distribution patterns and that no single parameter appears to control their distribution. Cross-validation, using a random subsample of the starting data and only using major ions, shows that predictions are generally accurate to within an order of magnitude. Using the same approach, an abridged version of the U.S. Geological Survey Produced Waters Database, Version 2.3 (which includes both data from produced and geothermal waters) was mapped to the ESOM and predicted values were generated for samples that contained enough variables to be effectively mapped. Results show that in general, produced and geothermal waters are predicted to be enriched in REEs by an order of magnitude or more relative to seawater, with maximum predicted enrichments in excess of 1000-fold. Cartographic mapping of the resulting predictions indicates that maximum REE concentrations exceed values in seawater across the majority of geologic basins investigated and that REEs are typically spatially co-associated. The factors causing this co-association were not determined from ESOM analysis, but based on the information currently available, REE content in produced and geothermal waters is not directly controlled by lithology, reservoir temperature, or salinity.

Engle, Mark A. (ORCID:0000000152587374)↗

Origin of the isotopic composition of natural perchlorate: Experimental results for the impact of reaction pathway and initial ClO x reactant

Natural perchlorate (ClO 4 - ) exists in many places on Earth, in lunar regolith, meteorites, and on the surface of Mars. Terrestrial natural ClO 4 - has widely variable Cl and O stable isotopic compositions (δ 37 Cl, δ 18 O, Δ 17 O). The δ 18 O and Δ 17 O values of ClO 4 - from the most hyper-arid locations co-vary. ClO 4 - from less arid areas has relatively little 17 O excess and poor Δ 17 O-δ 18 O correlation. ClO 4 - from the Atacama Desert has unusually low δ 37 Cl (<-10‰) and exhibits a positive correlation between δ 37 Cl and δ 18 O, while the δ 37 Cl of ClO 4 - from all other locations varies between -5 and +7‰ with no δ 37 Cl-δ 18 O covariation. To evaluate the impact of different precursors (ClO x ) and reaction pathways on the isotopic composition of ClO 4 - , we measured the isotopic composition of ClO 4 - produced in the laboratory by UV or O 3 mediated aqueous oxidation of Cl-, OCl-, ClO2-, and ClO2° as well as O 3 mediated oxidation of dry NaCl. ClO x oxidation in aqueous or dry systems enriched in O 3 produced ClO 4 - with Δ 17 O values that generally increased with the number of O atoms required and included evidence that the site-specific 17 O anomaly in O 3 was preferentially transferred to ClO 4 - . Based on the inferred number of O atoms sourced from O 3 , and known Cl and O reaction pathways, it appears that ClO 2 ° and ClO 3 * were required intermediates in the production of ClO 4 - in the O 3 experiments. ClO x aqueous oxidation by UV irradiation produced ClO 4 - with a large range of δ 18 O values and little or no 17 O anomaly. ClO 3 - was produced to a much greater extent than ClO 4 - in all experiments except dry oxidation of NaCl by O 3 . The isotopic composition of ClO 3 - was distinct from that of ClO 4 - produced from the same initial reactants. Combined results of O 3 and UV mediated reactions largely bracketed the range of natural ClO 4 - δ 18 O and Δ 17 O values as well as δ 37 Cl values of non-Atacama natural samples, but no conditions produced the low δ 37 Cl values of Atacama ClO 4 - . Finally, our results indicate that variation in production mechanisms, possibly combined with isotopically variable precursors, could be responsible for much of the observed isotopic variation in natural ClO 4 - and ClO 3 - .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Normative Ranges for Oculomotor and Reaction Time Tests in U.S. Military Service Members and Veterans

AbstractBackground Oculomotor and reaction time tests are frequently used assessments of vestibular symptoms, traumatic brain injury (TBI), or other neurological disorders in both clinical and research contexts. When interpreting these tests it is important to have a reference interval (RI) as a comparison for what constitutes a typical/expected response; however, the current body of research has only limited information regarding normative ranges calculated according to established standards or for a military-specific sample.Purpose The purpose of the present study was to describe RIs for oculomotor and reaction time tests in a cohort of service members and veterans (SMVs) for use as comparators by clinicians and scientists.Research Design Descriptive.Study Sample Participants were prospectively enrolled in the Defense and Veterans Brain Injury Center-Traumatic Brain Injury Center of Excellence 15-year Longitudinal Traumatic Brain Injury Study. Only SMVs without a history of TBI or blast exposure were included in the RI calculations.Data Collection and Analysis The test paradigms included in this analysis were: smooth pursuit, prosaccades, antisaccades, saccades and reaction time, predictive saccades, optokinetic nystagmus, auditory reaction time, and visual reaction time. Nonparametric methods, based on the U.S. Food and Drug Administration's recognized consensus standards, were used to calculate 95% RIs. A comparison between the calculated RIs and those available from previously published research is provided.Results Summary statistics and RIs were calculated for 47 outcome parameters from 13 oculomotor and reaction time tests. Sample sizes and age ranges varied across outcome parameters depending on the availability of reference values for RI calculations. The sample sizes used to calculate RIs ranged from 51 to 69. The age of SMVs included in each RI ranged from 19 to 61 years with mean ages ranging from 37 to 39 years. Similarities/differences between the RIs in the present study and those in previously published research are highly dependent on the outcome parameter; however, in general, the RIs in the present study tended to be somewhat wider.Conclusion The RIs provided in this paper can serve as comparisons for clinicians and scientists who are utilizing these oculomotor and reaction time testing paradigms in similar cohorts of patients or research participants.

Audiology & Speech-Language Pathology↗

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↗

Past, present and future—sample environments for materials research studies in scattering and spectroscopy; a UK perspective

Small Angle X-Ray Scattering (SAXS) and X-Ray Absorption Fine Structure (XAFS) are two techniques that have been employed at synchrotron sources ever since their inception. Over the course of the development of the techniques, the introduction of sample environments for added value experiments has grown dramatically. This article reviews past successes, current developments and an exploration of future possibilities for these two X-ray techniques with an emphasis on the developments in the United Kingdom between 1980-2020.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗