Search NASA⌕ Search

SEARCH · Search NASA

Results for “Curve 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 307 records · Page 17

Bringing 2D Eclipse Mapping out of the Shadows with Leave-one-out Cross Validation

Abstract Eclipse mapping is a technique for inferring 2D brightness maps of transiting exoplanets from the shape of an eclipse light curve. With JWST’s unmatched precision, eclipse mapping is now possible for a large number of exoplanets. However, eclipse mapping has only been applied to two planets, and the nuances of fitting eclipse maps are not yet fully understood. Here, we use Leave-one-out Cross Validation (LOO-CV) to investigate eclipse mapping, with application to a JWST NIRISS/SOSS observation of the ultrahot Jupiter WASP-18b. LOO-CV is a technique that provides insight into the out-of-sample predictive power of models on a data-point-by-data-point basis. We show that constraints on planetary brightness patterns behave as expected, with large-scale variations driven by the phase-curve variation in the light curve and smaller-scale structures constrained by the eclipse ingress and egress. For WASP-18b we show that the need for higher model complexity (smaller-scale features) is driven exclusively by the shape of the eclipse ingress and egress. We use LOO-CV to investigate the relationship between planetary brightness map components when mapping under a positive-flux constraint to better understand the need for complex models. Finally, we use LOO-CV to understand the degeneracy between the competing “hot spot” and “plateau” brightness map models of WASP-18b, showing that the plateau model is driven by the ingress shape and the hot spot model is driven by the egress shape, but preference for neither model is due to outliers or unmodeled signals. Based on this analysis, we make recommendations for the use of LOO-CV in future eclipse-mapping studies.

79 ASTRONOMY AND ASTROPHYSICS↗

The active CGCG 077-102 NED02 galaxy within the Abell 2063 galaxy cluster

Context.Within the framework of investigating the link between the central super massive black holes in the cores of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, which is a member of the CGCG 077-102 galaxy pair within the Abell 2063 cluster of galaxies. Aims.Our goal is to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. Methods.We usedChandraandXMM-Newtonarchival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in active galactic nuclei. Finally, we used public ZTF imaging data to investigate the optical variability. Results.There is no evidence of multiple X-ray sources or extended components within CGCG 077-102 NED02. Single X-ray spectral models fit the source well. We detect significant, nonrandom inter-observation 0.5–10 keV X-ray flux variabilities, for observations separated by ∼4 days for short-term variations and by up to ∼700 days for long-term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. The classification of CGCG 077-102 NED02 is also consistent with its X-ray luminosity of over 10 42 erg s −1 . We do not detect short-term variability in the optical ZTF light curves. However, we find a significant long-term stochastic variability in theg-band that can be well described by the damped random walk model with a best-fit characteristic damping timescale ofτ DRW = 30 −12 +28 days. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, with a long enough history within this massive structure to have been affected by the influence of this cluster for a long time. Conclusions.Our observations point toward a moderately massive black hole in the center of CGCG 077-102 NED02 of ∼10 6 M ⊙ . As compared to another similar pair in the literature, CGCG 077-102 NED02 is not heavily obscured, perhaps because of the surrounding intracluster medium ram-pressure stripping.

Astronomy & Astrophysics↗

The first Hubble diagram and cosmological constraints using superluminous supernovae

In this work, we present the first Hubble diagram of superluminous supernovae (SLSNe) out to a redshift of two, together with constraints on the matter density, Ω M , and the dark energy equation-of-state parameter, w(≡p/ρ). We build a sample of 20 cosmologically useful SLSNe I based on light curve and spectroscopy quality cuts. We confirm the robustness of the peak–decline SLSN I standardization relation with a larger data set and improved fitting techniques than previous works. We then solve the SLSN model based on the above standardization via minimization of the χ2 computed from a covariance matrix that includes statistical and systematic uncertainties. For a spatially flat Λ cold dark matter (ΛCDM) cosmological model, we find $\Omega _{\rm M}=0.38^{+0.24}_{-0.19}$, with an rms of 0.27 mag for the residuals of the distance moduli. For a w 0 w a CDM cosmological model, the addition of SLSNe I to a ‘baseline’ measurement consisting of Planck temperature together with Type Ia supernovae, results in a small improvement in the constraints of w 0 and wa of 4 per cent. We present simulations of future surveys with 868 and 492 SLSNe I (depending on the configuration used) and show that such a sample can deliver cosmological constraints in a flat ΛCDM model with the same precision (considering only statistical uncertainties) as current surveys that use Type Ia supernovae, while providing a factor of 2–3 improvement in the precision of the constraints on the time variation of dark energy, w 0 and w a . This paper represents the proof of concept for superluminous supernova cosmology, and demonstrates they can provide an independent test of cosmology in the high-redshift (z > 1) universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Heat capacity and thermodynamic functions of partially dehydrated sodium and zinc zeolite A (LTA)

Zeolite A (LTA) is an industrially important zeolite that exhibits sorption-induced framework flexibility, the thermodynamics of which are poorly understood. In this work, we report heat capacity measurements on zinc and sodium zeolite A from 1.8 to 300 K and compare the heat capacity of water in sodium zeolite A with that of water in other zeolites. The heat capacity of zeolitic water varies significantly depending on the hydration level and identity of the host zeolite, and more tightly bound water exhibits strong inflections in its heat capacity curve. This suggests a combination of effects, including differences in water-framework binding strength and hydration-dependent flexibility transitions. We also report fits of the heat capacity data using theoretical functions, and we report values for $C_{P,m}^°$,$Δ_0^TS_m^°$,$Δ_0^TH_m^°$, and $Φ_m^°$ from 0 to 300 K. These results contribute to a systematic thermodynamic understanding of the effects of cation exchange, guest molecule confinement, and sorbate-dependent flexibility transitions in zeolites.

Geochemistry & Geophysics↗

The Optical Luminosity–Time Correlation for More than 100 Gamma-Ray Burst Afterglows

Gamma-ray bursts (GRBs) are fascinating events due to their panchromatic nature. Their afterglow emission is observed from sub-TeV energies to radio wavelengths. We investigate GRBs that present an optical plateau, leveraging on the resemblance with the X-ray plateau shown in many GRB light curves (LCs). We comprehensively analyze all published GRBs with known redshifts and optical plateau observed mostly by the Neil Gehrels Swift Observatory (Swift). We fit 267 optical LCs and show the existence of the plateau in 102 cases, which is the largest compilation so far of optical plateaus. For 56 Swift GRBs with optical and X-ray plateaus, we compare the rest-frame end time at both wavelengths ( , ), and conclude that the plateau is achromatic between and . We also confirm the existence of the two-dimensional relations between and the optical luminosity at the end of the plateau emission, which resembles the same luminosity–time correlation in X-rays (Dainotti et al. 2013 ). Additionally, the existence of this optical correlation has been demonstrated for the largest sample of optical plateaus in the literature to date. The squared scatter in this optical correlation is smallest for the subset of the Gold GRBs with a decrease in the scatter equivalent to 52.4% when compared to the scatter of the entire GRB sample.

79 ASTRONOMY AND ASTROPHYSICS↗

A mathematical approach to using the forgetting curve to evaluate experience and training factors in human reliability analysis

Traditional human reliability analysis (HRA) methods have difficulty dealing with the dynamic nature of factors such as time and rely on static and expert-judgment-based assessments of performance-shaping factors (PSFs) across limited levels. In this study, we introduce a mathematical approach for dynamically evaluating the experience and training PSF. Our proposed method integrates the psychological concept of the “forgetting curve” to evaluate how PSFs are impacted by the number of trainings and the time elapsed since training. To confirm the validity of the model, we provide experimental data fitted by identifying the quantitative relationship between training and human performance. This research enables dynamic and objective assessments, thus reducing reliance on subjective expert judgment and improving the accuracy of HRA.

99 - GENERAL AND MISCELLANEOUS↗

Three-omega thermal-conductivity measurements with curved heater geometries

The three-omega method, a powerful technique to measure the thermal conductivity of nanometer-thick films and the interfaces between them, has historically employed straight conductive wires to act as both heaters and thermometers. When investigating stochastically prepared samples such as two-dimensional materials and nanomembranes, residue and excess material can make it difficult to fit the required millimeter-long straight wire on the sample surface. There are currently no available criteria for how diverting three-omega heater wires around obstacles affects the validity of the thermal measurement. In this Letter, we quantify the effect of the wire curvature by performing three-omega experiments with a wide range of frequencies using both curved and straight heater geometries on SiO 2 /Si samples. When the heating wire is curved, we find that the measured Si substrate thermal conductivity changes by only 0.2%. Similarly, we find that wire curvature has no significant effect on the determination of the thermal resistance of an ~65 nm SiO 2 layer, even for the sharpest corners considered here, for which the largest measured ratio of the thermal penetration depth of the applied thermal wave to radius of curvature of the heating wire is 4.3. This result provides useful design criteria for three-omega experiments by setting a lower bound for the maximum ratio of the thermal penetration depth to wire radius of curvature.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Physical Interpretation of Early Battery Life Prediction Models

Early battery life prediction models are most useful for R&D if they help us understand the early changes in battery electrochemical response that correspond with long-term degradation and failure. Linear regression models such as Fused lasso and Partial Least Squares can fit coefficients directly to high-dimensional electrochemical data like capacity-voltage and ΔV–state-of-charge, i.e., Q(V) and ΔV(SOC) curves, learning coefficients that can be physically interpreted. We leverage the ISU-ILCC battery aging data set to learn high-dimensional coefficients for early battery life prediction from traditional slow-rate capacity check data, demonstrating learning on Q(V), d Q· d V −1 , and ΔV(SOC) curves. A thorough study on the dependence of coefficient values on train/test size and data preprocessing methods is made, demonstrating the reliability of high-dimensional regression approaches unless very small amounts of data are used for model training. For this data set, coefficients from Q(V) and d Q· d V −1 models highlight changes in electrode stoichiometry due to lithium loss, while ΔV(SOC) coefficients highlight changes in positive electrode diffusivity due to particle cracking as well as electrode stoichiometry shifts. By directly interpreting the coefficients of a regression model, we make physical insights into battery degradation mechanisms without requiring the assumptions of traditional battery data analysis methods.

25 ENERGY STORAGE↗

Definition of a new (Doniach-Sunjic-Shirley) peak shape for fitting asymmetric signals applied to reduced graphene oxide/graphene oxide XPS spectra

The existence of asymmetry in X-ray photoelectron spectroscopy (XPS) photoemission lines is widely accepted, but line shapes designed to accommodate asymmetry are generally lacking in theoretical justification. Here in this work, we present a new line shape for describing asymmetry in XPS signals that is based on two facts. First, the most widely known line shape for fitting asymmetric XPS signals that has a theoretical basis, referred to as the Doniach-Sunjic (DS) line shape, suffers from a mathematical inconvenience, which is that for asymmetric shapes the area beneath the curve (above the x-axis) is infinite. Second, it is common practice in XPS to remove the inelastically scattered background response of a peak in question with the Shirley algorithm. The new line shape described herein attempts to retain the theoretical virtues of the DS line shape, while allowing the use of a Shirley background, with the consequence that the resulting line shape has a finite area. To illustrate the use of this Doniach-Sunjic-Shirley (DSS) line shape, a set of spectra obtained from varying amounts of graphene oxide (GO) and reduced GO on a patterned, heterogeneous surface are fit and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗

Machine learning for neutron reflectometry data analysis of two-layer thin films

Neutron reflectometry (NR) is a powerful tool for probing thin films at length scales down to nanometers. We investigated the use of a neural network to predict a two-layer thin film structure to model a given measured reflectivity curve. Application of this neural network to predict a thin film structure revealed that it was accurate and could provide an excellent starting point for traditional fitting methods. Employing prediction-guided fitting has considerable potential for more rapidly producing a result compared to the labor-intensive but commonly-used approach of trial and error searches prior to refinement. A deeper look at the stability of the predictive power of the neural network against statistical fluctuations of measured reflectivity profiles showed that the predictions are stable. We conclude that the approach presented here can provide valuable assistance to users of NR and should be further extended for use in studies of more complex n-layer thin film systems. This result also opens up the possibility of developing adaptive measurement systems in the future.

36 MATERIALS SCIENCE↗

Photometry of Outer Solar System Objects from the Dark Energy Survey. I. Photometric Methods, Light-curve Distributions, and Trans-Neptunian Binaries

We report the methods of and initial scientific inferences from the extraction of precision photometric information for the >800 trans-Neptunian objects (TNOs) discovered in the images of the Dark Energy Survey (DES). Scene-modeling photometry is used to obtain shot-noise-limited flux measures for each exposure of each TNO, with background sources subtracted. Comparison of double-source fits to the pixel data with single-source fits are used to identify and characterize two binary TNO systems. A Markov Chain Monte Carlo method samples the joint likelihood of the intrinsic colors of each source as well as the amplitude of its flux variation, given the time series of multiband flux measurements and their uncertainties. A catalog of these colors and light-curve amplitudes A is included with this publication. We show how to assign a likelihood to the distribution q(A) of light-curve amplitudes in any subpopulation. Using this method, we find decisive evidence (i.e., evidence ratio <0.01) that cold classical (CC) TNOs with absolute magnitude 6 < H$_{r}$ < 8.2 are more variable than the hot classical (HC) population of the same H$_{r}$, reinforcing theories that the former form in situ and the latter arise from a different physical population. Resonant and scattering TNOs in this H$_{r}$ range have variability consistent with either the HCs or CCs. DES TNOs with H$_{r}$ < 6 are seen to be decisively less variable than higher-H$_{r}$ members of any dynamical group, as expected. More surprising is that detached TNOs are decisively less variable than scattering TNOs, which requires them to have distinct source regions or some subsequent differential processing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

OGLE-2018-BLG-0532Lb: Cold Neptune with Possible Jovian Sibling

We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio q=M{sub planet}/M{sub host}≃1×10{sup −4}. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that χ {sup 2} can be significantly improved by adding either a third lens (3L1S, Δχ {sup 2} = 81) or second source (2L2S, Δχ {sup 2} = 77) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are a host mass M ∼ 0.25 M {sub ⊙}, system distance D {sub L} ∼ 1 kpc, and planet mass m{sub p,1}=8 M{sub ⊕} with projected separation a{sub 1,⊥}=1.4 au. However, there is a relatively bright I = 18.6 (and also relatively blue) star projected within <50 mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental and theoretical investigation of FeCrVAl and related compounds

Abstract We have carried out a combined theoretical and experimental investigation of FeCrVAl, and the effect of Mn and Co doping on its structural, magnetic, and electronic band properties. Our first principles calculations indicate that FeCrVAl, FeCr 0.5 Mn 0.5 VAl, and FeCr 0.5 Co 0.5 VAl exhibit nearly perfect spin polarization, which may be further enhanced by mechanical strain. At the same time, FeCrV 0.5 Mn 0.5 Al and FeCrV 0.5 Co 0.5 Al exhibit a relatively small value of spin polarization, making them less attractive for practical applications. Using arc melting and high vacuum annealing, we synthesized three compounds FeCrVAl, FeCr 0.5 Mn 0.5 VAl, and FeCr 0.5 Co 0.5 VAl, which are predicted to exhibit high spin polarization. The room temperature x-ray diffraction patterns of all samples are fitted with full B2 type disorder with a small amount of FeO 2 secondary phase. All samples show very small saturation magnetizations at room temperature. The thermomagnetic curves M(T) of FeCrVAl and FeCr 0.5 Co 0.5 VAl are similar to that of a paramagnetic material, whereas that of FeCr 0.5 Mn 0.5 VAl indicates ferrimagnetic behavior with the Curie temperature of 135 K. Our findings may be of interest for researchers working on Heusler compounds for spin-based electronic applications.

Physics↗

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗

One Planet or Two Planets? The Ultra-sensitive Extreme-magnification Microlensing Event KMT-2019-BLG-1953

We present the analysis of a very high-magnification (A ∼ 900) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary-lens (2L1S) model with a mass ratio of q ∼ 2 × 10{sup −3} improves the fits by Δχ {sup 2} = 181.8, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with Δχ {sup 2} = 16.0 and 13.5 for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of M{sub host}=0.31{sub −0.17}{sup +0.37} M{sub ⊙} and that the planetary system is located at a distance of D{sub L}=7.04{sub −1.33}{sup +1.10} kpc toward the Galactic center. The mass of the securely detected planet is M{sub p}=0.64{sub −0.35}{sup +0.76} M{sub J}. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by follow-up observations, and thus the event illustrates the need for intensive follow-up observations for very high-magnification events even in the current generation of high-cadence surveys.

79 ASTRONOMY AND ASTROPHYSICS↗