Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

Probing the limits of cosmological information from the Lyman- α forest 2-point correlation functions

The standard cosmological analysis with the Lyα forest relies on a continuum fitting procedure that suppresses information on large scales and distorts the three-dimensional correlation function on all scales. In this work, we present the first cosmological forecasts without continuum fitting distortion in the Lyα forest, focusing on the recovery of large-scale information. Using idealized synthetic data, we compare the constraining power of the full shape of the Lyα forest auto-correlation and its cross-correlation with quasars using the baseline continuum fitting analysis versus the true continuum. We find that knowledge of the true continuum enables a ∼ 10% reduction in uncertainties on the Alcock-Paczyński (AP) parameter and the matter density, Ω m . We also explore the impact of large-scale information by extending the analysis up to separations of 240 h -1 Mpc along and across the line of sight. The combination of these analysis choices can recover significant large-scale information, yielding up to a ∼ 15% improvement in AP constraints. This improvement is analogous to extending the Lyα forest survey area by ∼ 40%.

Lyman alpha forest↗

Core plasma fueling by fast inward particle transport after hydrogen pellet injection in Wendelstein 7-X

A large database of more than 1000 individual cryogenic hydrogen pellets injected into Wendelstein 7-X for plasma fueling was analyzed to improve the understanding of the three phases of the process: the ablation, deposition and transport of the pellet material. Kilohertz-sampled electron density and temperature measurements revealed a more complex drift behavior than predicted by numerical code simulation. It could be explained by the poloidal plasma E r x B- drift rotation, which plays a significant role in stellarators, but was not previously considered in pellet injection codes like HPI2. The drift results in a fast poloidal rotation of the pellet material around the plasma core, leading to an almost homogeneous deposition over the involved flux surfaces regardless of magnetic high and low field side injection geometry. Additionally, a novel fast inward directed transport mechanism (‘FIT-effect’) was observed. The effect occurs on timescales of tens of milliseconds and cannot be explained by neoclassical transport or diffusion. It might be linked to the turbulence pinch recently found in Wendelstein 7-X. When the FIT-effect occurs, the pellet particles are rapidly transferred from the deposition flux surfaces to the plasma core, causing the plasma density profile to peak, which is beneficial for confinement in Wendelstein 7-X. The large pellet injection database was statistical analyzed with regard to pellet and plasma parameters, which delivered some starting points towards developing an understanding of the physics behind the FIT-effect. The results indicate, that plasma core fueling via pellet injection is largely independent of the injection geometry in stellarators under certain conditions, reducing the technical complexity of the injection system.

Wendelstein 7-X↗

Use of Neutrino Scattering Events with Low Hadronic Recoil to Inform Neutrino Flux and Detector Energy Scale

Charged-current neutrino interactions with low hadronic recoil ("low-nu") have a cross-section that is approximately constant versus neutrino energy. These interactions have been used to measure the shape of neutrino fluxes as a function of neutrino energy at accelerator-based neutrino experiments such as CCFR, NuTeV, MINOS and MINERvA. In this paper, we demonstrate that low-nu events can be used to measure parameters of neutrino flux and detector models and that utilization of event distributions over the upstream detector face can discriminate among parameters that affect the neutrino flux model. From fitting a large sample of low-nu events obtained by exposing MINERvA to the NuMI medium-energy beam, we find that the best-fit flux parameters are within their a priori uncertainties, but the energy scale of muons reconstructed in the MINOS detector is shifted by 3.6% (or 1.8 times the a priori uncertainty on that parameter). These fit results are now used in all MINERvA cross-section measurements, and this technique can be applied by other experiments operating at MINERvA energies, such as DUNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗

Deep energy-pressure regression for a thermodynamically consistent EOS model

Abstract In this paper, we aim to explore novel machine learning (ML) techniques to facilitate and accelerate the construction of universal equation-Of-State (EOS) models with a high accuracy while ensuring important thermodynamic consistency. When applying ML to fit a universal EOS model, there are two key requirements: (1) a high prediction accuracy to ensure precise estimation of relevant physics properties and (2) physical interpretability to support important physics-related downstream applications. We first identify a set of fundamental challenges from the accuracy perspective, including an extremely wide range of input/output space and highly sparse training data. We demonstrate that while a neural network (NN) model may fit the EOS data well, the black-box nature makes it difficult to provide physically interpretable results, leading to weak accountability of prediction results outside the training range and lack of guarantee to meet important thermodynamic consistency constraints. To this end, we propose a principled deep regression model that can be trained following a meta-learning style to predict the desired quantities with a high accuracy using scarce training data. We further introduce a uniquely designed kernel-based regularizer for accurate uncertainty quantification. An ensemble technique is leveraged to battle model overfitting with improved prediction stability. Auto-differentiation is conducted to verify that necessary thermodynamic consistency conditions are maintained. Our evaluation results show an excellent fit of the EOS table and the predicted values are ready to use for important physics-related tasks.

97 MATHEMATICS AND COMPUTING↗

Integrating AlphaFold and deep learning for atomistic interpretation of cryo-EM maps

Abstract Interpretation of cryo-electron microscopy (cryo-EM) maps requires building and fitting 3D atomic models of biological molecules. AlphaFold-predicted models generate initial 3D coordinates; however, model inaccuracy and conformational heterogeneity often necessitate labor-intensive manual model building and fitting into cryo-EM maps. In this work, we designed a protein model-building workflow, which combines a deep-learning cryo-EM map feature enhancement tool, CryoFEM (Cryo-EM Feature Enhancement Model) and AlphaFold. A benchmark test using 36 cryo-EM maps shows that CryoFEM achieves state-of-the-art performance in optimizing the Fourier Shell Correlations between the maps and the ground truth models. Furthermore, in a subset of 17 datasets where the initial AlphaFold predictions are less accurate, the workflow significantly improves their model accuracy. Our work demonstrates that the integration of modern deep learning image enhancement and AlphaFold may lead to automated model building and fitting for the atomistic interpretation of cryo-EM maps.

59 BASIC BIOLOGICAL SCIENCES↗

A semblance measure for model comparison

Algorithmic and computational advances have made it possible that geophysical survey and earth model design can be aided by many systematic trial inverse-modelling runs with synthetic data. Such may, for example, come up in machine-learning approaches. Automated image appraisal pertaining to such applications will involve common statistical tests for goodness-of-data fit as a primary evaluation method. However, solution non-uniqueness may render multiple images equivalent in terms of their data fit, requiring secondary categorizers. A logical choice for classifying synthetic-imaging results quantifies the goodness of model fit where a known reference model replaces the observational input. The task of model intercomparison in terms of measuring the resemblance to the reference model poses challenges to common distance-based metrics like root mean square error and mean absolute error. First, distance-based metrics can introduce spurious contributions when smooth models with fuzzy target contours are to be compared against a sharp reference. Second, large differences due to parameter-estimation overshoots can dominate distance metrics. Here, we propose a remedy that is referred to as semblance and is based on the idea of logistic functions, where a binary-dependent variable adds non-zero or zero accumulation terms for the, respectively, passing or failing of preset target thresholds. This classifying approach is amenable to an objective where model feature recognition is primary. Numerical comparisons to distance-based metrics provide evidence for the advantages of the semblance in view of this objective. Geophysical imaging in conjunction with machine-learning is seen as a benefitting upcoming application area.

58 GEOSCIENCES↗

Dosage optimization for reducing tumor burden using a phenotype-structured population model with a drug-resistance continuum

Abstract Drug resistance is a significant obstacle to effective cancer treatment. To gain insights into how drug resistance develops, we adopted a concept called fitness landscape and employed a phenotype-structured population model by fitting to a set of experimental data on a drug used for ovarian cancer, olaparib. Our modeling approach allowed us to understand how a drug affects the fitness landscape and track the evolution of a population of cancer cells structured with a spectrum of drug resistance. We also incorporated pharmacokinetic (PK) modeling to identify the optimal dosages of the drug that could lead to long-term tumor reduction. We derived a formula that indicates that maximizing variation in plasma drug concentration over a dosing interval could be important in reducing drug resistance. Our findings suggest that it may be possible to achieve better treatment outcomes with a drug dose lower than the levels recommended by the drug label. Acknowledging the current limitations of our work, we believe that our approach, which combines modeling of both PK and drug resistance evolution, could contribute to a new direction for better designing drug treatment regimens to improve cancer treatment.

Life Sciences & Biomedicine - Other Topics↗

Black hole spin in X-ray binaries: giving uncertainties an $\mathcal{f}$

The two established techniques for measuring black hole spin in X-ray binaries often yield conflicting results, which must be resolved before either method may be deemed robust. In practice, black hole spin measurements based on fitting the accretion disc continuum effectively do not marginalize over the colour-correction factor f col . This factor parametrizes spectral hardening of the disc continuum by the disc atmosphere, whose true properties are poorly constrained. We incorporate reasonable systematic uncertainties in f col into the eight (non-maximal) black hole spin measurements vetted by the disc continuum fitting community. In most cases, an f col uncertainty of ±0.2–0.3 dominates the black hole spin error budget. We go on to demonstrate that plausible departures in f col values from those adopted by the disc continuum fitting practitioners can bring the discrepant black hole spins into agreement with those from iron line modelling. Systematic uncertainties in f col , such as the effects of strong magnetization, should be better understood before dismissing their potentially dominant impact on the black hole spin error budget.

4U 1543–47↗

Understanding the extreme luminosity of DES14X2fna

We present DES14X2fna, a high-luminosity, fast-declining Type IIb supernova (SN IIb) at redshift z = 0.0453, detected by the Dark Energy Survey (DES). DES14X2fna is an unusual member of its class, with a light curve showing a broad, luminous peak reaching M r ≃ –19.3 mag 20 d after explosion. This object does not show a linear decline tail in the light curve until ≃60 d after explosion, after which it declines very rapidly (4.30 ± 0.10 mag 100 d –1 in the r band). By fitting semi-analytic models to the photometry of DES14X2fna, we find that its light curve cannot be explained by a standard 56 Ni decay model as this is unable to fit the peak and fast tail decline observed. Inclusion of either interaction with surrounding circumstellar material or a rapidly-rotating neutron star (magnetar) significantly increases the quality of the model fit. Furthermore, we also investigate the possibility for an object similar to DES14X2fna to act as a contaminant in photometric samples of SNe Ia for cosmology, finding that a similar simulated object is misclassified by a recurrent neural network (RNN)-based photometric classifier as an SN Ia in ~1.1–2.4 per cent of cases in DES, depending on the probability threshold used for a positive classification.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling the stellar halo with RR-Lyrae stars

ABSTRACT A seven-parameter distribution function (DF) is fitted to $20\, 000$ RR-Lyrae stars for which only astrometric data are available. The observational data are predicted by the DF in conjunction with the gravitational potential of a self-consistent model Galaxy defined by DFs for the dark halo, the bulge, and a four-component disc. Tests of the technique developed to deal with missing line-of-sight velocities show that adding such velocities tightens constraints on the DF only slightly. The recovered model of the RR-Lyrae population confirms that the population is flattened and has a strongly radially biased velocity distribution. At large radii, its density profile tends to ρ ∼ r−4.5 but no power law provides a good fit inside the solar sphere. The model is shown to provide an excellent fit to the data for stars brighter than r = 16.5 but at certain longitudes it predicts too few faint stars at Galactocentric radii $\sim 20\, \mathrm{kpc}$, possibly signalling that the halo is not axisymmetric. The DF is used to predict the velocity distribution of BHB stars for which space velocities are available. The z components are predicted successfully but too much anisotropy in the vRvϕ plane is expected.

Li, Chengdong↗

Optical variability of quasars with 20-yr photometric light curves

ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the growth rate of structure from the small-scale clustering of the luminous red galaxy sample

ABSTRACT We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\, h^{-1}{\rm Mpc}$ with a model based on the aemulus cosmological emulator to measure the growth rate of cosmic structure, parametrized by fσ8. We obtain a measurement of fσ8(z = 0.737) = 0.408 ± 0.038, which is 1.4σ lower than the value expected from 2018 Planck data for a flat ΛCDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1\text{--}60\, h^{-1}{\rm Mpc}$ modelled by the aemulus cosmological emulator and find a 4.5σ tension in the amplitude of the halo velocity field with the Planck + ΛCDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.

79 ASTRONOMY AND ASTROPHYSICS↗

The DEHVILS survey overview and initial data release: high-quality near-infrared Type Ia supernova light curves at low redshift

ABSTRACT While the sample of optical Type Ia supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good ‘standard candles,’ is still ≲ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from 0.002 < z < 0.09 as a part of the Dark Energy, H0, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT’s WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for Y, J, and H-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands (Y,J, and H) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, snoopy, and B ayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128 to 0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM Y-band magnitudes.

Astronomy & Astrophysics↗

The Dark Energy Survey Year 3 and eBOSS: constraining galaxy intrinsic alignments across luminosity and colour space

ABSTRACT We present direct constraints on galaxy intrinsic alignments (IAs) using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift z ∼ 0.2–1.0, luminous red galaxies from eBOSS at z ∼ 0.8, and also an SDSS-III BOSS CMASS sample at z ∼ 0.5. We measure two-point IA correlations, which we fit using a model that includes lensing, magnification, and photometric redshift error. Fitting on scales 6 Mpc h−1 < rp < 70 Mpc h−1, we make a detection of IAs in each sample, at 5σ–22σ (assuming a simple one-parameter model for IAs). Using these red samples, we measure the IA–luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best-fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim 2\!-\!18~{{\ \rm per\ cent}}$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of emission-line galaxies from eBOSS at z ∼ 0.8. We constrain the non-linear alignment amplitude to be $A_1=0.07^{+0.32}_{-0.42}$ (|A1| < 0.78 at 95 per cent CL).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Empirically Driven multiwavelength K-c orrections at low redshift

K -corrections – a necessary ingredient for converting between flux in observed bands to flux in rest-frame bands – are critical for comparing galaxies at differing redshifts. These corrections often rely on fits to empirical or theoretical spectral energy distribution (SED) templates of galaxies. However, templates can only produce reliable K-corrections in regimes where SED models are robust. For instance, the templates utilized in some popular software packages are not well-constrained in some bands (e.g. WISE W4 in KCORRECT ), which results in ill-behaved K -corrections. We address this shortcoming by developing an empirically driven approach to K -corrections that limits the dependence on SED templates. We perform a polynomial fit for the K -correction as a function of a galaxy’s rest-frame colour determined in a pair of well-constrained bands (e.g. 0 (g − r)) and redshift, exploiting the fact that galaxy SEDs can be approximated as a one-parameter family at low redshift. For bands well-constrained by SED templates, our empirically driven K -corrections yield results comparable to the SED fitting methods used by KCORRECT and the GSWLC-M2 catalogue (the updated medium-deep GALEX–SDSS–WISE Legacy Catalogue). However, our method dramatically outperforms Kcorrect derived K -corrections for WISE W4 . Our method is also robust to incorrect template assumptions outside of the optical bands and enforces that the K -correction must be zero at z = 0. Our K -corrected photometry and code are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

The DESI one-per cent survey: exploring the halo occupation distribution of luminous red galaxies and quasi-stellar objects with AbacusSummit

ABSTRACT We present the first comprehensive halo occupation distribution (HOD) analysis of the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey luminous red galaxy (LRG) and Quasi Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 h−1 Mpc in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of $f_\mathrm{sat} = 11\pm 1~{y{\ \mathrm{per\,cent}}}$, a mean halo mass of $\log _{10}\overline{M}_h/M_\odot =13.40^{+0.02}_{-0.02}$, and a linear bias of $b_\mathrm{lin} = 1.93_{-0.04}^{+0.06}$. For LRGs in 0.6 < z < 0.8, we find $f_\mathrm{sat}=14\pm 1~{{\ \mathrm{per\,cent}}}$, $\log _{10}\overline{M}_h/M_\odot =13.24^{+0.02}_{-0.02}$, and $b_\mathrm{lin}=2.08_{-0.03}^{+0.03}$. For QSOs, we infer $f_\mathrm{sat}=3^{+8}_{-2}\mathrm{per\,cent}$, $\log _{10}\overline{M}_h/M_\odot = 12.65^{+0.09}_{-0.04}$, and $b_\mathrm{lin} = 2.63_{-0.26}^{+0.37}$ in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high fidelity galaxy mocks, forming the basis of systematic tests for DESI Y1 cosmological analyses. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealling significant and interesting trends in mean halo mass, linear bias, and satellite fraction.

79 ASTRONOMY AND ASTROPHYSICS↗

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗