Search NASA⌕ Search

SEARCH · Search NASA

Results for “sample average approximation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The polarized-signal density matrix: A practical way to recover molecular frame information from isotropic samples

We present a novel approach to model ultrafast time-dependent nonlinear optical polarization sensitive signals emitted from randomly oriented molecules. By projecting the laboratory-frame analyzer polarization axis into the molecular-frame and linking that axis with the density matrix through a tensor product, we demonstrate an approach to find a specific molecular orientation that yields a good approximation to simulated four-wave mixing signals produced by the same model but with averaging over molecular orientation.

Thurston, Richard L↗
DOE OSTI2024

A Method for Dynamic Kolsky Bar Compression at High Temperatures: Application to Ti-6Al-4V

Abstract An experimental apparatus for measuring the dynamic behavior of materials subjected to strain rates on the order of 10 $$^3$$ 3 s $$^{-1}$$ - 1 and temperatures up to 800°C with a unique triple actuation system is developed in this work. This system is based on the traditional Kolsky (or split-Hopkinson pressure) bar design, with the addition of an external furnace used to heat the specimen to the desired temperature. A synchronized triple pneumatic actuation system is used to control the motion and timing of the the sample, incident, and transmitted bars. The cold contact time (CCT), or the time during which the heated sample is in contact with the room temperature bars before compression, is measured experimentally and carefully controlled to minimize the development of a temperature gradient across the sample and avoid heating of the bars. Experiments are performed in conjunction with ultra high speed imaging and 2D digital image correlation (DIC), as well as high speed thermal imaging. To verify the viability of the proposed system, experiments were conducted on Ti-6Al-4V (wt.%) at temperatures from 25°C up to and 800°C at an average strain rate of approximately 1200 s $$^{-1}$$ - 1 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Grain size dependent high-pressure elastic properties of ultrafine micro/nanocrystalline grossular

We have performed sound velocity and unit cell volume measurements of three synthetic, ultrafine micro/nanocrystalline grossular samples up to 50 GPa using Brillouin spectroscopy and synchrotron X-ray diffraction. The samples are characterized by average grain sizes of 90 nm, 93 nm and 179 nm (hereinafter referred to as samples Gr90, Gr93, and Gr179, respectively). The experimentally determined sound velocities and elastic properties of Gr179 sample are comparable with previous measurements, but slightly higher than those of Gr90 and Gr93 under ambient conditions. However, the differences diminish with increasing pressure, and the velocity crossover eventually takes place at approximately 20–30 GPa. The X-ray diffraction peaks of the ultrafine micro/nanocrystalline grossular samples significantly broaden between 15–40 GPa, especially for Gr179. The velocity or elasticity crossover observed at pressures over 30 GPa might be explained by different grain size reduction and/or inhomogeneous strain within the individual grains for the three grossular samples, which is supported by both the pressure-induced peak broadening observed in the X-ray diffraction experiments and transmission electron microscopy observations. The elastic behavior of ultrafine micro/nanocrystalline silicates, in this case, grossular, is both grain size and pressure dependent.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improved Measurement of C P Violation Parameters in B s 0 → J / ψ K + K − Decays in the Vicinity of the ϕ ( 1020 ) Resonance

The decay-time-dependent C P asymmetry in B s 0 → J / ψ ( → μ + μ − ) K + K − decays is measured using proton-proton collision data, corresponding to an integrated luminosity of 6 fb − 1 , collected with the LHCb detector at a center-of-mass energy of 13 TeV. Using a sample of approximately 349 000 B s 0 signal decays with an invariant K + K − mass in the vicinity of the ϕ ( 1020 ) resonance, the C P -violating phase ϕ s is measured, along with the difference in decay widths of the light and heavy mass eigenstates of the B s 0 − B ¯ s 0 system, Δ Γ s , and the difference of the average B s 0 and B 0 meson decay widths, Γ s − Γ d . The values obtained are ϕ s = − 0.039 ± 0.022 ± 0.006 rad , Δ Γ s = 0.0845 ± 0.0044 ± 0.0024 ps − 1 , and Γ s − Γ d = − 0.005 6 − 0.0015 + 0.0013 ± 0.0014 ps − 1 , where the first uncertainty is statistical and the second systematic. These are the most precise single measurements to date and are consistent with expectations based on the Standard Model and with the previous LHCb analyses of this decay. These results are combined with previous independent LHCb measurements. The phase ϕ s is also measured independently for each polarization state of the K + K − system and shows no evidence for polarization dependence. © 2024 CERN, for the LHCb Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Outlier classification using autoencoders: application for fluctuation driven flows in fusion plasmas

Understanding the statistics of fluctuation driven flows in the boundary layer of magnetically confined plasmas is desired to accurately model the lifetime of the vacuum vessel components. Mirror Langmuir probes (MLPs) are a novel diagnostic that uniquely allow us to sample the plasma parameters on a time scale shorter than the characteristic time scale of their fluctuations. Sudden large-amplitude fluctuations in the plasma degrade the precision and accuracy of the plasma parameters reported by MLPs for cases in which the probe bias range is of insufficient amplitude. While some data samples can readily be classified as valid and invalid, we find that such a classification may be ambiguous for up to 40% of data sampled for the plasma parameters and bias voltages considered in this study. In this contribution, we employ an autoencoder (AE) to learn a low-dimensional representation of valid data samples. By definition, the coordinates in this space are the features that mostly characterize valid data. Ambiguous data samples are classified in this space using standard classifiers for vectorial data. In this way, we avoid defining complicated threshold rules to identify outliers, which require strong assumptions and introduce biases in the analysis. By removing the outliers that are identified in the latent low-dimensional space of the AE, we find that the average conductive and convective radial heat fluxes are between approximately 5% and 15% lower as when removing outliers identified by threshold values. For contributions to the radial heat flux due to triple correlations, the difference is up to 40%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Patterns and controls on island‐wide aboveground biomass accumulation in second‐growth forests of Puerto Rico

Abstract Understanding the heterogeneity of biomass accumulation in second‐growth tropical forests following land use abandonment is important for informing ecosystem carbon models and forest restoration efforts. There is an urgent need for a broad sample of second‐growth forests to enhance our knowledge of carbon accumulation in human‐dominated landscapes, especially for older forests. Puerto Rico has predominantly second‐growth forests, ranging in age from approximately 25 to more than 80 years. We used an island‐wide sample of airborne lidar from the NASA Goddard Lidar, Hyperspectral, and Thermal (G‐LiHT) Airborne Imager collected on March 2017, forest inventory data, and data on forest age, precipitation, soils, and land use to estimate aboveground biomass stocks in moist and wet, second‐growth tropical forests. Biomass accumulation rates in Puerto Rico were lower, on average, than in other Neotropical forests. Median biomass across >16,700 ha of older second‐growth forests was 105 Mg ha −1 , and sampled biomass rarely surpassed 250 Mg ha −1 . Differences in biomass by age were large and persistent across different substrates and land uses, with a plateau in the pattern of island‐wide biomass accumulation after about 33 years. A spatial regression model showed that multiple factors were related to biomass accumulation, including time since abandonment, geologic substrate, past land use as coffee or pasture, precipitation, topographic wetness index, and slope. Our findings have important consequences for the total carbon storage and expected climate mitigation benefits of large‐scale reforestation efforts, and highlight the value of airborne lidar for quantifying biomass variability in complex tropical landscapes. Abstract in Spanish is available with online material.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Laboratory-Scale Melter Systems for Simulated and Real AP-107 Tank Waste - 20483

The Hanford Waste Treatment and Immobilization Plant (WTP) will process and stabilize waste that is stored in underground tanks on the Hanford Site. At the WTP, the tank waste will be combined with glass-forming and modifying additives to form an aqueous and solid mixture called melter feed, which will then be vitrified in joule-heated melters. A small-scale melter system, the continuous laboratory-scale melter (CLSM), is designed to operate with a continuous feeding process while periodically pouring glass product and collecting off-gas. Two CLSMs were constructed: one in a non-radioactive environment for processing waste simulants and another in a fume hood capable of handling radioactive material for processing actual Hanford tank waste. The waste sample received for vitrification in the CLSM was from Hanford tank 241-AP-107. Approximately 12 liters of melter-feed slurry was vitrified into 7 kilograms of glass product over the course of 10 hours of charging to the CLSM. The partitioning of technetium-99 in the product streams was used to calculate the technetium-99 recovery and help determine when the vitrification process had reached a steady state. The average single-pass retention of technetium-99 in the glass product was then calculated for the processing of AP-107 waste in the CLSM system. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Variation in carbon and nitrogen concentrations among peatland categories at the global scale

Peatlands account for 15 to 30% of the world’s soil carbon (C) stock and are important controls over global nitrogen (N) cycles. However, C and N concentrations are known to vary among peatlands contributing to the uncertainty of global C inventories, but there are few global studies that relate peatland classification to peat chemistry. We analyzed 436 peat cores sampled in 24 countries across six continents and measured C, N, and organic matter (OM) content at three depths down to 70 cm. Sites were distinguished between northern (387) and tropical (49) peatlands and assigned to one of six distinct broadly recognized peatland categories that vary primarily along a pH gradient. Peat C and N concentrations, OM content, and C:N ratios differed significantly among peatland categories, but few differences in chemistry with depth were found within each category. Across all peatlands C and N concentrations in the 10–20 cm layer, were 440 ± 85.1 g kg -1 and 13.9 ± 7.4 g kg -1 , with an average C:N ratio of 30.1 ± 20.8. Among peatland categories, median C concentrations were highest in bogs, poor fens and tropical swamps (446–532 g kg -1 ) and lowest in intermediate and extremely rich fens (375–414 g kg -1 ). The C:OM ratio in peat was similar across most peatland categories, except in deeper samples from ombrotrophic tropical peat swamps that were higher than other peatlands categories. Peat N concentrations and C:N ratios varied approximately two-fold among peatland categories and N concentrations tended to be higher (and C:N lower) in intermediate fens compared with other peatland types. This study reports on a unique data set and demonstrates that differences in peat C and OM concentrations among broadly classified peatland categories are predictable, which can aid future studies that use land cover assessments to refine global peatland C and N stocks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Sampling Microbial Dynamics in the Salish Sea Estuary: Evaluating Methods to Capture Cyanobacteria and Cyanophage

Introduction: Picocyanobacteria from the genera Prochlorococcus and Synechococcus thrive across the globe in aquatic environments, have relatively small genomes, and have growth dynamics regulated by both viral interactions and abiotic conditions, making them excellent model organisms for exploring host-pathogencoevolution. Methods: We developed and refined methods to sample and sequence cyanobacteria, cyanophages, and measured features of their abiotic environment. Results: The protocol described herein can successfully discriminate large-cell eukaryotic organisms, but size fractionation of picocyanobacteria appears to be affected by the presence of free DNA, multicellular structures, and abundant tycheposons. Our preferred final protocol from this exploratory effort included a combination of in-line and single vacuum flask filtrations, which reduced filtration processing time by over threefold in some cases compared to other tested methods, such as a fully in-line sequence or in-site filtrations. We successfully extracted an average of approximately 400–1200 ng for all filter fractions, with some variations between kits. Discussion: The protocol described herein can successfully discriminate large-cell eukaryotic organisms, but size fractionation of picocyanobacteria appears to be affected by the presence of free DNA, multicellular structures, and abundant tycheposons.

Salish Sea↗
DOE OSTI2025

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗
DOE OSTI2025

Electrically Conductive Amine Functionalized Reduced Graphite Oxide Foam for CO 2 Removal from the Air

Rapid regeneration of CO 2 adsorbents is critical to improving the productivity of direct air capture (DAC) systems. In this study, we codesigned a material to have appropriate electrical conductivity and CO 2 adsorption properties to enable efficient CO 2 capture from air. Specifically, we present a poly(ethylenimine) (PEI)-impregnated thermally annealed graphite oxide (TAGO900) foam adsorbent tailored for vacuum-assisted electrically driven thermal swing adsorption (V-ETSA). This structured adsorbent leverages the high electrical conductivity of the reduced graphite oxide framework to enable fast and direct heating of the adsorbent material by electrical resistance heating (Joule heating). An optimal sample, 40 wt % PEI (molecular weight 25k) impregnated TAGO900, shows the best balance between adsorption capacity (1.54 mmol g –1 ) and adsorption rates (0.016 mmol g –1 min –1 ) using fixed bed breakthrough experiments at 25 °C and 70% RH using 50 sccm 400 ppm of CO 2 /N 2 flow. Compared to conventional temperature vacuum swing adsorption (TVSA), the V-ETSA approach achieves substantially faster CO 2 desorption, achieving average desorption rates (including cooling time) of 0.09 mmol g –1 min –1 ─approximately 2.5 times faster than TVSA under similar operating conditions. The maximum desorption rate reaches 0.23 mmol/g/min during the desorption stage. These results underscore the importance of the direct heating strategy, such as Joule heating, for fast and highly productive vacuum swing adsorption in DAC systems.

amines↗
DOE OSTI2026

Photoinduced correlations in stochastic dynamics of a solid-state ionic conductor

Photoexcitation by ultrashort laser pulses plays a crucial role in controlling reaction pathways, creating nonequilibrium material properties, and probing complex molecular dynamics. The photoresponse following a laser pulse is generally nonidentical between exposures due to spatiotemporal fluctuations or the stochastic nature of dynamical pathways. However, most ultrafast pump-probe experiments struggle to distinguish intrinsic sample fluctuations from extrinsic apparatus noise, often missing deviations from the averaged response. Leveraging the stability and high photon flux of time-resolved X-ray micro-diffraction at a synchrotron, we characterized stochastic photoinduced dynamics in a solid-state ionic conductor. By analyzing temporal evolutions of the lattice parameter of a single grain, we found that shot-to-shot fluctuations are not independent. Instead, correlations exist between nonequilibrium lattice trajectories following adjacent shots, with a characteristic correlation length of approximately 1500 shots, corresponding to an energy barrier of 0.4 ± 0.1 eV, close to the activation energy of lithium-ion diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analysis of the Surface Morphology and Chemical Composition of Zr-Nb3Sn Alloys with Zr Different Concentrations

The inclusion of zirconium (Zr) in niobium-tin (Nb3Sn) significantly enhances the performance of Nb3Sn radiofrequency cavities in high magnetic fields. This research project is dedicated to characterizing the surface properties and chemical composition of Zr-doped Nb3Sn. Two different concentrations of Zr were used for doping: approximately 0.5% and 24%. Various spectroscopy techniques were employed to analyze how the surface morphology and chemical composition of Nb3Sn change with increasing Zr content. The findings of this study indicate that the grain size decreases as the Zr concentration increases. In samples with 24% Zr, nanometer-sized particles, likely oxides, were observed. X-ray photoelectron spectroscopy (XPS) data revealed that the thickness of niobium oxides (NbOx) decreases with increasing Zr, while the thickness of tin oxides (SnOx) increases. Additionally, grain size distribution analysis showed that the average grain size is around 5300 nm , with a grain area density of about 165 grains/μm .

Sue, Micah↗
DOE OSTI2021

LaNi x Fe 1– x O 3–δ as a Robust Redox Catalyst for CO 2 Splitting and Methane Partial Oxidation

The current study reports LaNi 0.5 Fe 0.5 O 3–δ as a robust redox catalyst for CO 2 splitting and methane partial oxidation at relatively low temperatures (~700 °C) in the context of a hybrid redox process. Specifically, perovskite-structured LaNi x Fe 1–x O 3–δ (LNFs) with nine different compositions (x = 0.05–0.5) were prepared and investigated. Among the samples evaluated, LaNi 0.4 Fe 0.6 O 3–δ and LaNi 0.5 Fe 0.5 O 3–δ showed superior redox performance, with ~90% CO 2 and methane conversions and >90% syngas selectivity. The standalone LNFs also demonstrated performance comparable to that of LNF promoted by mixed conductive Ce 0.85 Gd 0.1 Cu 0.05 O 2–δ (CGCO). Long-term testing of LaNi 0.5 Fe 0.5 O 3–δ indicated that the redox catalyst gradually loses its activity over repeated redox cycles, amounting to approximately 0.02% activity loss each cycle, averaged over 500 cycles. This gradual deactivation was found to be reversible by deep oxidation with air. Further characterizations indicated that the loss of activity resulted from a slow accumulation of iron carbide (Fe 3 C and Fe 5 C 2 ) phases, which cannot be effectively removed during the CO 2 splitting step. Reoxidation with air removed the carbide phases, increased the availability of Fe for the redox reactions via solid-state reactions with La 2 O 3 , and decreased the average crystallite size of La 2 O 3 . As a result, reactivating the redox catalyst periodically, e.g., once every 40 cycles, was shown to be highly effective, as confirmed by operating the redox catalyst over 900 cumulative cycles while maintaining satisfactory redox performance.

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

Realistic operation of two residential cordwood-fired outdoor hydronic heater appliances—Part 2: Particle number and size

The use of wood as a fuel for home heating is a concern from an environmental health and safety perspective as biomass combustion appliances emit high concentrations of particulate matter. Wood burning significantly contributes to wintertime particulate matter concentrations in many states in the northern United States. Of particular concern are outdoor wood-fired hydronic heaters. These devices are concerning as they tend to have very large combustion chambers and typical use patterns can result in long periods of low output, which result in an increased particulate matter emission rate relative to high heat output operating conditions. In this study, the performance of two hydronic heaters operating under different combustion conditions, including four different heat output categories approximately corresponding to categories I–IV denoted in Environmental Protection Agency Method 28 Outdoor Wood-fired Hydronic Heaters, and during start-up and reloading events were investigated. Measurements of flue gas particulate number concentration and size for particles with aerodynamic diameters between 0.006 and 10 µm were made using a dilution sampling system. The measured particle number concentration in the flue gas was between 0.71 and 420 million particles per cubic centimeter and was dependent on fuel loading and heat output. For each hydronic heater tested, the highest average particle concentration was found at the beginning of each test during the cold-start condition. Additionally, the majority of the particles had aerodynamic diameters less than 0.100 µm (particles of this size made up between 64% and 97% of all particles) and less than 1% of all particles had aerodynamic diameters greater than 1 µm for all phases. For particles in the accumulation mode, between 0.100 and 1 µm, the mean particle diameter was dependent on fuel loading and heat output.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A differentiable model of the evolution of dark matter halo concentration

ABSTRACT We introduce a new model of the evolution of the concentration of dark matter haloes, c(t). For individual haloes, our model approximates c(t) as a power law with a time-dependent index, such that at early times, concentration has a nearly constant value of c ≈ 3–4, and as cosmic time progresses, c(t) smoothly increases. Using large samples of halo merger trees taken from the Bolshoi–Planck and MultiDark Planck 2 cosmological simulations, we demonstrate that our three-parameter model can approximate the evolution of the concentration of individual haloes with a typical accuracy of 0.1 dex for $t\gtrsim 2\, {\rm Gyr}$ for all Bolshoi–Planck and MultiDark Planck 2 haloes of present-day peak mass $M_{0}\gtrsim 10^{11.5}\, {\rm M}_{\odot }$. We additionally present a new model of the evolution of the concentration of halo populations, which we show faithfully reproduces both average concentration growth and the diversity of smooth trajectories of c(t), including capturing correlations with halo mass and halo assembly history. Our publicly available source code, diffprof, can be used to generate Monte Carlo realizations of the concentration histories of cosmologically representative halo populations. diffprof is differentiable due to its implementation in the jax autodiff library, which facilitates the incorporation of our model into existing analytical halo model frameworks.

Stevanovich, Dash↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗