Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constrained parameter spaces”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Leaf 13 C data constrain the uncertainty of the carbon dynamics of temperate forest ecosystems

Stable carbon isotope discrimination occurred in plant biophysical and biogeochemical processes can help understand plant physiology and soil biogeochemistry with respect to carbon cycling. Here, we incorporated the photosynthetic carbon isotope discrimination into a process-based land surface model (iTEM) to test if stable carbon isotope composition (δ 13 C) can impose additional constraint on model parameters. Sequential data assimilation was implemented at six eddy covariance flux tower sites using carbon flux observations including gross primary productivity (GPP) and net ecosystem exchange (NEE) with and without considering foliar δ 13 C (δ 13 C f ) measurement constraints, respectively. Our model-data fusion showed that δ 13 C f can provide useful constraint on photosynthetic (V cmax25 , the maximum rate of carboxylation at 25°C) and stomatal (g1, the slope of stomatal function) parameters as well as the posterior carbon fluxes. When including δ 13 C f measurement, g1 spatially varies among the six sites and is significantly correlated with annual precipitation. We incorporated the statistical relationship between g1 and annual precipitation into iTEM, which is then used to quantify the regional carbon dynamic in temperate forest ecosystems of the Northern Hemisphere. Compared with the simulation only conditioned on carbon flux observations, regional carbon flux estimations performed slightly better against the FLUXCOM products and the uncertainties of modeled carbon fluxes were reduced by 27%. Our study demonstrated that δ 13 C f data constrains carbon flux uncertainties across space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The growth of density perturbations in the last ~10 billion years from tomographic large-scale structure data

In order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from z ~ 2 using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range 0.2 < z < 0.7, where it is lower than predicted by . This difference is mostly driven by current cosmic shear data, although the growth histories reconstructed from different data combinations are consistent with each other, and we find no evidence of systematic deviations in any particular experiment. In spite of the tension with , the data are well-described by the model, albeit with a lower value of S 8 ≡ σ 8 (Ω m /0.3) 0.5 . As part of our analysis, we find constraints on this parameter of S 8 = 0.7781 ± 0.0094 (68% confidence level), reaching almost percent-level errors comparable with CMB measurements, and 3.4σ away from the value found by .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Cosmological constraints on light but massive relics

Many scenarios of physics beyond the standard model predict new light, weakly coupled degrees of freedom, populated in the early universe and remaining as cosmic relics today. Due to their high abundances, these relics can significantly affect the evolution of the universe. For instance, massless relics produce a shift ΔN eff to the cosmic expectation of the effective number of active neutrinos. Massive relics, on the other hand, additionally become part of the cosmological dark matter in the later universe, though their light nature allows them to freely stream out of potential wells. This produces novel signatures in the large-scale structure (LSS) of the universe, suppressing matter fluctuations at small scales. Here, we present the first general search for such light (but massive) relics (LiMRs) with cosmic microwave background (CMB) and LSS data, scanning the 2D parameter space of their masses m X and temperatures $T^{(0)}_{X}$ today. In the conservative minimum-temperature ($T^{(0)}_{X}$ = 0.91 K) scenario, we rule out Weyl (and higher-spin) fermions -- such as the gravitino -- with m X ≥ 2.26 eV at 95% C.L., and set analogous limits of m X ≤ 11.2,1.06,1.56 eV for scalar, vector, and Dirac-fermion relics. This is the first search for LiMRs with joint CMB, weak-lensing, and full-shape galaxy data; we demonstrate that weak-lensing data is critical for breaking parameter degeneracies, while full-shape information presents a significant boost in constraining power relative to analyses with only baryon acoustic oscillation parameters. Under the combined strength of these datasets, our constraints are the tightest and most comprehensive to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

FolpsD: combining EFT and phenomenological approaches for joint power spectrum and bispectrum analyses

We present a theoretical model for the power spectrum and bispectrum of galaxy clustering that exploits the complementarity between small-scale power spectrum information and large-scale bispectrum measurements. We extend the FOLPS code by combining its one-loop EFT galaxy power spectrum with a tree-level galaxy bispectrum projected onto the tripolar spherical harmonics (Sugiyama) basis. To access additional small-scale information, we also consider a line-of-sight damping factor in both statistics, mirroring approaches commonly used in studies of redshift-space distortions. We test the model using DESI DR2 galaxy mocks. Even without damping, the joint analysis of the EFT power spectrum and bispectrum significantly improves constraints and reduces parameter degeneracies relative to power spectrum analyses alone. For LRG-like samples, including the damping further extends the range beyond $k\sim 0.3 \,h \text{Mpc}^{-1}$ in the power spectrum and $k \sim 0.24 \,h \text{Mpc}^{-1}$ in the bispectrum without introducing statistically significant parameter biases. This leads to up to $\sim 30\%$ tighter constraints on $A_s$ and $ω_{cdm}$. For low signal-to-noise tracers such as QSOs, however, the damping parameters are weakly constrained and can absorb noise fluctuations, leading to shifts in inferred parameters. Similar limitations may arise in models where cosmological information is encoded in power-spectrum shape features degenerate with the damping, such as scenarios with massive neutrinos. In contrast, for $w_0w_a$CDM we obtain $15\%$ and $21\%$ tighter constraints on $w_0$ and $w_a$, respectively, yielding a deviation from constant dark energy at slightly more than the $1σ$ level using full-shape information alone. The code is publicly available at https://github.com/cosmodesi/FolpsD

Bansal, P. [Michigan U., MCTP; Michigan U.] (ORCID↗
DOE OSTI2023

Search for Ultraheavy Dark Matter from Observations of Dwarf Spheroidal Galaxies with VERITAS

Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV ≲Mχ ≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; Mχ ≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS γ-ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a γ-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗
DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks. Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol –1 , and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Dark energy survey year 3 results: Cosmology with peaks using an emulator approach

ABSTRACT We constrain the matter density Ωm and the amplitude of density fluctuations σ8 within the ΛCDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We take into account the uncertainty from the remaining, largely unconstrained ΛCDM parameters (Ωb, ns, and h). We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv ~\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68 per cent confidence limits), with a precision of 1.8 per cent, an improvement of 38 per cent compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in S8 is recorded. We find a mild tension of $1.5 \, \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower S8. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment AIA and S8, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Planck constraints and gravitational wave forecasts for primordial black hole dark matter seeded by multifield inflation

Here, we perform a Markov Chain Monte Carlo analysis of a simple yet generic multifield inflation model characterized by two scalar fields coupled to each other and nonminimally coupled to gravity, fit to Planck 2018 cosmic microwave background (CMB) data. In particular, model parameters are constrained by data on the amplitude of the primordial power spectrum of scalar curvature perturbations on CMB scales A s , the spectral index n s , and the ratio of power in tensor to scalar modes r, with a prior that the primordial power spectrum should also lead to primordial black hole (PBH) production sufficient to account for the observed dark matter abundance. We find that n s in particular controls the constraints on our model. Whereas previous studies of PBH formation from an ultraslow-roll phase of inflation have highlighted the need for at least one model parameter to be highly fine-tuned, we identify a degeneracy direction in parameter space such that shifts by ~10% of one parameter can be compensated by comparable shifts in other parameters while preserving a close fit between model predictions and observations. Furthermore, we find this allowed parameter region produces observable gravitational wave signals in the frequency ranges to which upcoming experiments are projected to be sensitive, including Advanced LIGO and Virgo, the Einstein Telescope, Cosmic Explorer, DECIGO, and LISA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Juno mission as a probe of long-range new physics

Orbits of celestial objects, especially the geocentric and heliocentric ones, have been well explored to constrain new long-range forces beyond the Standard Model (SM), often referred to as fifth forces. In this paper, for the first time, we apply the motion of a spacecraft around Jupiter to probe fifth forces that don’t violate the equivalence principle. The spacecraft is the Juno orbiter, and ten of its early orbits already allow a precise determination of the Jovian gravitational field. We use the shift in the precession angle as a proxy to test non-gravitational interactions between Juno and Jupiter. Requiring that the contribution from the fifth force does not exceed the uncertainty of the precession shift inferred from data, we find that a new parameter space with the mass of the fifth-force mediator around 10 −14 eV is excluded at 95% C.L.

new light particles↗
DOE OSTI2025

Renormalized classical theory of quantum magnets

Here, we derive a renormalized classical spin (RCS) theory for 𝑆 >1/2 quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP 1 phase space of dipolar SU(2) coherent states. When the spin Hamiltonian $\hat{ℋ}$(𝑆) is linear in the spin operators $\hat{𝑺}$ 𝑗 for each lattice site 𝑗, the RCS Hamiltonian $\tilde{ℋ}$ cl coincides with the usual classical model ℋ cl = lim 𝑆→∞⁡ $\hat{ℋ}$(𝑆). In the presence of nonlinear terms, however, the RCS theory is more accurate than ℋ cl . For the many materials modeled by spin Hamiltonians with (nonlinear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron-scattering data.

magnetic anisotropy↗
DOE OSTI2025

Clarity through the neutrino fog: constraining new forces in dark matter detectors

The PandaX-4T and XENONnT experiments present indications of Coherent Elastic Neutrino Nucleus Scattering (CEνNS) from 8 B solar neutrinos at 2.6σ and 2.7σ, respectively. This constitutes the first observation of the neutrino “floor” or “fog”, an irreducible background that future dark matter searches in terrestrial detectors will have to contend with. Here, we first discuss the contributions from neutrino–electron scattering and from the Migdal effect in the region of interest of these experiments, and we argue that they are non-negligible. Second, we make use of the recent PandaX-4T and XENONnT data to derive novel constraints on light scalar and vector mediators coupling to neutrinos and quarks. We demonstrate that these experiments already provide world-leading laboratory constraints on new light mediators in some regions of parameter space.

neutrino interactions↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

ACCESS: Confirmation of No Potassium in the Atmosphere of WASP-31b

We present a new optical (400–950 nm) transmission spectrum of the hot Jupiter WASP-31b (M = 0.48 M {sub J}; R = 1.54 R {sub J}; P = 3.41 days), obtained by combining four transit observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet, as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na i and K i, with particular focus on K i, for which there have been two previously published disagreeing results. Observations with Hubble Space Telescope (HST)/STIS detected K i, whereas ground-based low- and high-resolution observations did not. We use equilibrium and nonequilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (α=5.3{sub −3.1}{sup +2.9}), high-altitude clouds at a log cloud top pressure of −3.6{sub −2.1}{sup +2.7} bars, and possible muted H{sub 2}O features. We find that our observations support other ground-based claims of no K i. Clouds are likely why signals like H{sub 2}O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (α = 8.3 ± 1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Rotation‐Controlled Diurnal Evolution of Uranus' Asymmetric Bow Shock at Equinox

Uranus possesses the most extreme magnetic and rotational geometry in the solar system, resulting in a uniquely dynamic and asymmetric interaction between its magnetosphere and the solar wind. Here we investigate the diurnal evolution of the Uranian bow shock (BS) at equinox using global multifluid magnetohydrodynamic simulations constrained by Voyager 2 observations. The BS structure and variability are quantified using the stand‐off distance, terminator distance, flaring parameter, and BS parameter, enabling a systematic assessment of its global asymmetry over one planetary rotation. Our results show that the BS exhibits pronounced diurnal asymmetry that is primarily controlled by planetary rotation. Even under steady upstream solar wind conditions, the BS undergoes periodic expansion and contraction, reflecting rotation‐driven reconfiguration of the planetary magnetospheric topology. This diurnal modulation represents a remarkably large asymmetry in planetary BS geometry within the solar system. These results establish planetary rotation as the intrinsic primary driver of Uranus' BS dynamics at equinox and provide essential constraints for the design and formulation of future space missions to the ice giants. These findings also offer a unique framework for understanding shock physics in the abundant population of ice‐giant exoplanets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗