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 145 records · Page 8

Gravothermal collapse and the diversity of galactic rotation curves

The rotation curves of spiral galaxies exhibit a great diversity that challenges our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of σ/m ≈ $\mathscr{O}$(1) cm 2 /g, the predicted dark matter central densities are a good match to the observed densities in galaxies. Here, in this work, we explore a regime with a larger cross section of σ/m ≈ 20−40 cm 2 /g in dwarf galactic halos. We will show that such strong dark matter self-interactions can further amplify the diversity of halo densities inherited from their assembly history. High concentration halos can enter the gravothermal collapse phase within 10 Gyr, resulting in a high density, while low concentration ones remain in the expansion phase and have a low density. We fit the rotation curves of 14 representative low surface brightness galaxies and demonstrate how the large range of observed central densities is naturally accommodated in the strong SIDM regime of σ/m ≈ 20 – 40 cm 2 /g. Galaxies that are outliers in the previous studies, due to their high halo central densities, are no longer outliers in this SIDM regime as their halos would be in the collapse phase. For galaxies with a low density, the SIDM fits are robust to the variation of the cross section. Our findings open up a new window for testing gravothermal collapse, the unique signature of strong dark matter self-interactions, and exploring a broader SIDM model space. As an example, we illustrate how the larger cross sections favored by our fits, together with upper limits from strong lensing observations in clusters, pick out the preferred SIDM model space for a dark matter particle coupled to a light gauge boson in the Born regime.

dark matter↗

SNIa Cosmology Analysis Results from Simulated LSST Images: From Difference Imaging to Constraints on Dark Energy

Abstract The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10 6 transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg 2 ) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg 2 of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are u = 23.66 mag, g = 24.69 mag, r = 24.06 mag, i = 23.45 mag, z = 22.54 mag, and y = 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg –2 in g band, and falling to 300 deg –2 in y band. The photometry has biases <1% for magnitudes 19.5 < m < 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts z < 1, we recover the input cosmology parameters to within statistical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotation curves and dynamical masses of MaNGA barred galaxies

In this paper we analyse a sample of 46 barred galaxies of MaNGA. Our goal is to investigate the stellar kinematics of these galaxies and obtain their rotation curves. Additionally, we aim to derive the total stellar and dynamical masses, as well as the maximum rotation velocity, in order to examine their distributions and scaling relations. Using the Pipe3D dataproducts publicly available we obtained the rotation curves, which were fitted considering two components of an axisymmetric Miyamoto–Nagai gravitational potential. We found a wide range of the maximum rotation velocities (117–340 ${\rm km\, s^{-1}}$), with a mean value of 200 ${\rm km\, s^{-1}}$. In addition we found that the total stellar and dynamical masses are in the range of log(M star /M ⊙ ) = 10.1−11.5, with a mean value of log(M star /M ⊙ ) = 10.8, and log(M dyn /M ⊙ ) = 10.4−12.0, with a mean value of log(M dyn /M ⊙ ) = 11.1, respectively. We found a strong correlation between dynamical mass and maximum velocity, between maximum velocity and magnitude, and between stellar mass and maximum velocity. According to these results, barred galaxies exhibit similar behaviour to that of normal spiral galaxies with respect to these relations, as well as in terms of the distribution of their dynamical mass and maximum rotation velocity. However, we found that the distribution of stellar masses of barred galaxies is statistically different from other samples including non-barred galaxies. Finally, analysing the galaxies that show nuclear activity, we find no difference with the rotation curves of normal galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of a compact fast-neutron spectrometer for nuclear emergency response applications

We have developed a Compact Fast Neutron Spectrometer (CFNS) for passive assay of special nuclear material (SNM) through the observation of fast neutrons. The CFNS consists of eight organic glass scintillators (OGS) coupled to silicon photomultipliers and a waveform digitizer, which are integrated within a human-portable box. The CFNS determines the neutron energy profile by spectrum unfolding using the Maximum-Likelihood Expectation Maximization method. The detector acquisition system was optimized to have a dynamic range of up to 10 MeV neutron energy. Bulk special nuclear material (SNM) measurements from the National Criticality Experiments Research Center were analyzed for SNM validation/examination. Additionally, the results show that the CFNS can be used to distinguish between fission and (α, n) neutron emitters, regardless of intervening material type (Cu and polyethylene) and thickness, by taking the ratio of neutron counts at different regions in the unfolded energy spectrum. Additionally, by fitting an exponential curve to the unfolded energy spectrum of PuO 2 and Pu neutron emitters, the CFNS showed the ability of distinguishing between pure Pu oxide, pure Pu metal and mixed oxide-metal configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High-temperature 2D ferromagnetism in conjugated microporous porphyrin-type polymers

The need for magnetic 2D materials that are stable to the enviroment and have high Curie temperatures is very important for various electronic and spintronic applications. We have found that two-dimensional porphyrin-type aza-conjugated microporous polymer crystals are such a material (Fe-aza-CMPs). Fe-aza-CMPs are stable to CO, CO 2 , and O 2 atmospheres and show unusual adsorption, electronic, and magnetic properties. Indeed, they are semiconductors with small energy band gaps ranging from 0.27 eV to 0.626 eV. CO, CO 2 , and O 2 molecules can be attached in three different ways where single, double, or triple molecules are bound to iron atoms in Fe-aza-CMPs. For different attachment configurations we find that for CO and CO 2 a uniform distribution of the molecules is most energetically favorable while for O 2 molecules aggregation is most energetically preferable. The magnetic moments decrease from 4 to 2 to 0 for singly, doubly, triply occupied configurations for all gasses respectively. The most interesting magnetic properties are found for O 2 molecules attached to the Fe-aza-CMP. For a single attachment configuration we find that an antiferromagnetic state is favorable. When two O 2 molecules are attached, the calculations show the highest exchange integral with a value of J = 1071 μeV. This value has been verified by two independent methods where in the first method J is calculated by the energy difference between ferromagnetic and anitferromagnetic configurations. The second method is based on the frozen magnon approach where the magnon dispersion curve has been fitted by the Ising model. Finally, for the second method J has been estimated at J = 1100 μeV in excellent agreement with the first method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of a Heat Balance Experiment in the TREAT Reactor with SPARTA and OpenMC

This article presents the interpretation of a heat balance experiment in TREAT using the SPARTA methodology and an OpenMC model. The study examines the effects of control rods and temperature variations on the effective multiplication factor ($k_{\rm eff}$) and correction factors. Results indicate that both core temperature and rod positions have a linear effect on $k_{\rm eff}$ and correction factors in the upper-west channel. A critical hypothesis of SPARTA was validated: reactivities and correction factors can be obtained by linearly combining individual effects. The previously considered negligible effect of temperature on correction factors was found to be significant, necessitating updates to the SPARTA methodology. By fitting an exponential curve to mean temperature data, SPARTA's accuracy and convergence have been improved. Comparisons between SPARTA and OpenMC showed consistent reactivity removal due to temperature increase. While SPARTA effectively interprets complex transients, simpler methods may suffice when correction factors have minimal influence on reactivity calculations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. Here, in this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ~1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations. We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Climate change and the future productivity and distribution of crab in the Bering Sea

Crab populations in the eastern Bering Sea support some of the most valuable fisheries in the United States, but their future productivity and distribution are uncertain. Here, we explore observed changes in the productivity and distribution for snow crab, Tanner crab, and Bristol Bay red king crab. We link historical indices of environmental variation and predator biomass with observed time series of centroids of abundance and extent of crab stock distribution; we also fit stock–recruit curves including environmental indices for each stock. We then project these relationships under forcing from global climate models to forecast potential productivity and distribution scenarios. Our results suggest that the productivity of snow crab is negatively related to the Arctic Oscillation (AO) and positively related to ice cover; Tanner crab’s productivity and distribution are negatively associated with cod biomass and sea surface temperature. Aspects of red king crab distribution and productivity appear to be related to bottom temperature, ice cover, the AO, and/or cod biomass. Projecting these relationships forward with available forecasts suggests that Tanner crab may become more productive and shift further offshore, red king crab distribution may contract and move north, and productivity may decrease for snow crab as the population contracts northward.

54 ENVIRONMENTAL SCIENCES↗

Discovery of a Candidate Binary Supermassive Black Hole in a Periodic Quasar from Circumbinary Accretion Variability

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ~10 –3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67–002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302–102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear inductive response of pinned superconducting vortices in artificial pinning sites

Using the time dependent Ginzburg-Landau equations, we simulate the inductive responses of a variety of thin-film systems containing patterned antidots with different sizes and shapes. The results for all shapes show that the kinetic inductance diverges as the applied current approaches a critical current that is below the BCS depairing current. In conclusion, exploiting the similarity of the observed current-voltage behavior to that of Josephson junctions, we obtain an empirical equation that well fits the inductivity curve as a function of the applied current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Intermediate scattering potential strength in electron-irradiated YBa 2 Cu 3 O 7-δ from London penetration depth measurements

Temperature-dependent London penetration depth, λ(T), of a high quality optimally-doped YBa 2 Cu 3 O 7-δ single crystal was measured using tunnel-diode-resonator technique. Controlled artificial disorder was induced by low-temperature 20 K irradiation by 2.5 MeV electrons at two large doses of 3.8 x 10 19 and 5.3 x 10 19 electrons per cm 2 . The irradiation caused significant suppression of the superconductor's critical temperature, T c , from 94.6 K to 90.0 K, and to 78.7 K, respectively. The low-temperature behavior of λ(T) evolves from a T-linear in pristine state to a T 2 behavior after irradiation, expected for a line-nodal d-wave superconductor. However, the original theory that explained such behavior assumed a unitary limit of the scattering potential, whereas usually in normal metals and semiconductors, Born scattering is sufficient to describe the experiment. To estimate the scattering potential strength, we calculated the superfluid density, ρ s (t = T/T c ) = λ 2 (0)/λ 2 (T), varying the amount and strength of non-magnetic scattering using a self-consistent t-matrix theory. Fitting the obtained curves to a power law, ρ s = 1 - Rt n , and to a polynomial, ρ s = 1 - At - Bt 2 , and comparing the coefficients n in one set and A and B in another with the experimental values, we estimate the phase shift to be around 70 and 65°, respectively. We correlate this result with the evolution of the density of states with nonmagnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

jeanloz_method

A simple Jupyter notebook program for fitting isothermal compression curves to pressure-versus-volume data.

Brennan, Matthew↗

Spectral Variability of VHS J1256–1257b from 1 to 5 μm

Multiwavelength time-resolved observations of rotationally modulated variability from brown dwarfs and giant exoplanets are the most effective method for constraining their heterogeneous atmospheric structures. In a companion paper, we reported the discovery of strong near-infrared variability in HST/WFC3/G141 light curves of the very red L-dwarf companion VHS J1256–1257b. In this paper, we present a follow-up 36 hr Spitzer/IRAC Channel 2 light curve together with an in-depth analysis of the Hubble space telescope (HST) and the Spitzer data. The combined data set provides time-resolved light curves of VHS1256b sampling 1.1 to 4.5 μm. The Spitzer light curve is best fit with a single sine wave with a period of 22.04 ± 0.05 hr and a peak-to-peak amplitude of 5.76 ± 0.04%. Combining the period with a previously measured projected rotational velocity (v sin i), we find that VHS1256b is most consistent with equatorial viewing geometry. The HST/G141+Spitzer spectral energy distribution favors a model with a T{sub eff} of 1000 K and low surface gravity with disequilibrium chemistry. The spectral variability of VHS1256b is consistent with predictions from partly cloudy models, suggesting that heterogeneous clouds are the dominant source of the observed modulations. We find evidence at the 3.3σ level for amplitude variations within the 1.67 μm CH{sub 4} band, which is the first such detection for a variable L-dwarf. We compare the HST/G141 time-resolved spectra of three red L-dwarfs with high-amplitude near-infrared rotational modulations and find that although their time-averaged spectra are similar, their spectroscopic variabilities exhibit notable differences. This diversity reinforces the advantage of time-resolved spectroscopic observations for understanding the atmospheres of brown dwarfs and directly imaged exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗

A Differential Measurement of Circumstellar Extinction for AA Tau’s 2011 Dimming Event

AA Tau is a classical T Tauri star with a highly inclined, warped circumstellar disk. For decades, AA Tau exhibited photometric and spectroscopic variabilities that were successfully modeled as occultations of the primary star by circumstellar material. In 2011, AA Tau entered an extended faint state, presumably due to enhanced levels of circumstellar dust. We use two sets of contemporaneous optical–near-IR (NIR) spectra of AA Tau, obtained in 2008 and 2014 December, to directly measure the wavelength-dependent extinction associated with the dust enhancement driving AA Tau’s 2011 optical fade. Taken alone, AA Tau’s apparent optical–NIR increased extinction curve cannot be fit well with standard extinction laws. At optical wavelengths, AA Tau’s dimming event is consistent with predictions of common extinction models for an increase of A {sub V} = 2, but no such model reproduces AA Tau’s color–color excess at NIR wavelengths. We show that veiling emission accounts for the apparent anomalous NIR extinction: after including this veiling flux, AA Tau’s dimming event is consistent with a standard A {sub V} = 2 extinction law across the full optical–NIR range. We also report an increase in AA Tau’s mid-IR flux since its 2011 fade and suggest that an increase in the height of AA Tau’s inner disk is the most likely explanation for both the IR brightening and the additional extinction along the line of sight. In addition to informing our understanding of AA Tau, this analysis demonstrates that caution should be exercised when inferring extinction (and stellar parameters) from the NIR color–color excess of young stars with measurable NIR veiling fluxes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Lyα Transit Left Undetected: the Environment and Atmospheric Behavior of K2-25b

K2-25b is a Neptune-sized exoplanet (3.45R {sub ⊕}) that orbits its M4.5 host with a period of 3.48 days. Due to its membership in the Hyades Cluster, the system has a known age (727 ± 75 Myr). K2-25b’s youth and its similarities with Gl 436b suggested that K2-25b could be undergoing strong atmospheric escape. We observed two transits of K2-25b at Lyα using HST/STIS in order to search for escaping neutral hydrogen. We were unable to detect an exospheric signature, but placed an upper limit of (R {sub p}/R {sub ⋆})∣{sub Ly α} < 0.56 at 95% confidence by fitting the light curve of the Lyα red wing, or < 1.20 in the blue wing. We reconstructed the intrinsic Lyα profile of K2-25 to determine its Lyα flux, and analyzed XMM-Newton observations to determined its X-ray flux. Based on the total X-ray and extreme ultraviolet irradiation of the planet (8763 ± 1049 erg s{sup −1} cm{sup −2}), we estimated the maximum energy-limited mass-loss rate of K2-25b to be 10.6{sub −6.13}{sup +15.2}×10{sup 10} g s{sup −1} (0.56M {sub ⊕} per 1 Gyr), five times larger than the similarly estimated mass-loss rate of Gl 436b (2.2 × 10{sup 10} g s{sup −1}). The photoionization time is about 3 hr, significantly shorter than Gl 436b’s 14 hr. A nondetection of a Lyα transit could suggest K2-25b is not significantly losing its atmosphere, or factors of the system are resulting in the mass loss being unobservable (e.g., atmosphere composition or the system’s large high-energy flux). Further observations could provide more stringent constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Identifying RR Lyrae Variable Stars in Six Years of the Dark Energy Survey

We present a search for RR Lyrae stars using the full six-year data set from the Dark Energy Survey covering ∼5000 deg2 of the southern sky. Using a multistage multivariate classification and light-curve template-fitting scheme, we identify RR Lyrae candidates with a median of 35 observations per candidate. We detect 6971 RR Lyrae candidates out to ∼335 kpc, and we estimate that our sample is >70% complete at ∼150 kpc. We find excellent agreement with other wide-area RR Lyrae catalogs and RR Lyrae studies targeting the Magellanic Clouds and other Milky Way satellite galaxies. We fit the smooth stellar halo density profile using a broken-power-law model with fixed halo flattening (q = 0.7), and we find strong evidence for a break at with an inner slope of and an outer slope of . We use our catalog to perform a search for Milky Way satellite galaxies with large sizes and low luminosities. Using a set of simulated satellite galaxies, we find that our RR Lyrae-based search is more sensitive than those using resolved stellar populations in the regime of large (r h ≳ 500 pc), low-surface-brightness dwarf galaxies. A blind search for large, diffuse satellites yields three candidate substructures. The first can be confidently associated with the dwarf galaxy Eridanus II. The second has a distance and proper motion similar to the ultrafaint dwarf galaxy Tucana II but is separated by ∼5 deg. The third is close in projection to the globular cluster NGC 1851 but is ∼10 kpc more distant and appears to differ in proper motion.

79 ASTRONOMY AND ASTROPHYSICS↗

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

Structure evolution of Bi 4 O 7 under high pressure

X-ray diffraction and Raman spectroscopy measurements were used to study the structure evolution of Bi 4 O 7 during compression. Two isostructural phase transitions were observed under high pressure, with the first transition to phase II occurring at approximately 3.0 (2) GPa and the second to phase III occurring at approximately 13.0 (5) GPa. The pressure–volume ( P – V ) curve was found to be discontinuous at those two pressure transition boundaries. By fitting the P – V curves in three different stages, the bulk moduli were estimated to be approximately 101 (2), 68 (2) and 172 (3) GPa. Such large differences in the bulk moduli confirm two isostructural phase changes in the three stages. An irreversible amorphization process was found to start at approximately 20 (1) GPa, and the sample gradually becomes fully amorphous during compression at 30 (1) GPa.

Chemistry↗