Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model 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,675 records · Page 93

The CARMENES search for exoplanets around M dwarfs: Photospheric parameters of target stars from high-resolution spectroscopy. II. Simultaneous multiwavelength range modeling of activity insensitive lines

In this paper, we present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52–0.96 μ m) and NIR wavelength range (0.96–1.71 μ m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a χ 2 method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and NIR spectra combined are in better agreement than those derived from the same spectra in the NIR alone.

79 ASTRONOMY AND ASTROPHYSICS↗

A STRONGLY LENSED MASSIVE ULTRACOMPACT QUIESCENT GALAXY AT z {approx} 2.4 IN THE COSMOS/UltraVISTA FIELD

We report the discovery of a massive ultracompact quiescent galaxy that has been strongly lensed into multiple images by a foreground galaxy at z 0.960. This system was serendipitously discovered as a set of extremely K{sub s} -bright high-redshift galaxies with red J - K{sub s} colors using new data from the UltraVISTA YJHK{sub s} near-infrared survey. The system was also previously identified as an optically faint lens/source system using the COSMOS Advanced Camera for Surveys (ACS) imaging by Faure et al. Photometric redshifts for the three brightest images of the source galaxy determined from 27-band photometry place the source at z = 2.4 {+-} 0.1. We provide an updated lens model for the system that is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model implies that the magnification of the three brightest images is a factor of 4-5. We use the lens model, combined with the K{sub s} -band image, to constrain the size and Sersic profile of the galaxy. The best-fit model is an ultracompact galaxy (R{sub e} = 0.64{sup +0.08}{sub -0.18} kpc, lensing-corrected), with a Sersic profile that is intermediate between a disk and a bulge profile (n 2.2{sup +2.3}{sub -{sub 0.9}}), albeit with considerable uncertainties on the Sersic profile. We present aperture photometry for the source galaxy images that have been corrected for flux contamination from the central lens. The best-fit stellar population model is a massive galaxy (log(M{sub star}/M{sub Sun }) = 10.8{sup +0.1}{sub -0.1}, lensing-corrected) with an age of 1.0{sup +1.0}{sub -0.4} Gyr, moderate dust extinction (A{sub v} = 0.8{sup +0.5}{sub -0.6}), and a low specific star formation rate (log(SSFR) <-11.0 yr{sup -1}). This is typical of massive ''red-and-dead'' galaxies at this redshift and confirms that this source is the first bona fide strongly lensed massive ultracompact quiescent galaxy to be discovered. We conclude with a discussion of the prospects of finding a larger sample of these galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-F-BRD-08: A Novel Technique to Derive a Clinically-Acceptable Beam Model for Proton Pencil-Beam Scanning in a Commercial Treatment Planning System

Purpose: To evaluate the accuracy and limitations of a commercially-available treatment planning system’s (TPS’s) dose calculation algorithm for proton pencil-beam scanning (PBS) and present a novel technique to efficiently derive a clinically-acceptable beam model. Methods: In-air fluence profiles of PBS spots were modeled in the TPS alternately as single-(SG) and double-Gaussian (DG) functions, based on fits to commissioning data. Uniform-fluence, single-energy-layer square fields of various sizes and energies were calculated with both beam models and delivered to water. Dose was measured at several depths. Motivated by observed discrepancies in measured-versus-calculated dose comparisons, a third model was constructed based on double-Gaussian parameters contrived through a novel technique developed to minimize these differences (DGC). Eleven cuboid-dose-distribution-shaped fields with varying range/modulation and field size were subsequently generated in the TPS, using each of the three beam models described, and delivered to water. Dose was measured at the middle of each spread-out Bragg peak. Results: For energies <160 MeV, the DG model fit square-field measurements to <2% at all depths, while the SG model could disagree by >6%. For energies >160 MeV, both SG and DG models fit square-field measurements to <1% at <4 cm depth, but could exceed 6% deeper. By comparison, disagreement with the DGC model was always <3%. For the cuboid plans, calculation-versus-measured percent dose differences exceeded 7% for the SG model, being larger for smaller fields. The DG model showed <3% disagreement for all field sizes in shorter-range beams, although >5% differences for smaller fields persisted in longer-range beams. In contrast, the DGC model predicted measurements to <2% for all beams. Conclusion: Neither the TPS’s SG nor DG models, employed as intended, are ideally suited for routine clinical use. However, via a novel technique to be presented, its DG model can be tuned judiciously to yield acceptable results.

60 APPLIED LIFE SCIENCES↗

Analysis of rupture area of aftershocks caused by twin earthquakes (Case study: 11 April 2012 earthquakes of Aceh-North Sumatra)

The 11 April 2012 earthquakes off-shore Aceh-North Sumatra are unique events for the history of Indonesian earthquake. It is unique because that they have similar magnitude, 8.5 Mw and 8.1 Mw; close to epicenter distance, similar strike-slip focal mechanism, and occuring in outer rise area. The purposes of this research are: (1) comparing area of earthquakes base on models and that of calculation, (2) fitting the shape and the area of earthquake rupture zones, (3) analyzing the relationship between rupture area and magnitude of the earthquakes. Rupture area of the earthquake fault are determined by using 4 different formulas, i.e. Utsu and Seki (1954), Wells and Coppersmith (1994), Ellsworth (2003), and Christophersen and Smith (2000). The earthquakes aftershock parameters are taken from PGN (PusatGempabumiNasional or National Earthquake Information Center) of BMKG (Indonesia Agency Meteorology Climatology and Geophysics). The aftershock epicenters are plotted by GMT’s software. After that, ellipse and rectangular models of aftershock spreading are made. The results show that: (1) rupture areas were calculated using magnitude relationship which are larger than the the aftershock distributions model, (2) the best fitting model for that earthquake aftershock distribution is rectangular associated with Utsu and Seki (1954) formula, (3) the larger the magnitude of the earthquake, the larger area of the fault.

58 GEOSCIENCES↗

Latent Variable Network Modeling

lvnm Is a library for fitting latent variable network models within the R programming language. The code focuses on the Stochastic Block Model, Weighted Stochastic Block Model, and Dynamic Stochastic Block Model. Features include simulation methods, estimation methods, diagnostic plots and anomaly detection routines for network data.

Dabbs, Beau↗

Empirical Estimation of Shortest Route Length along U.S. Interstate Highways Based on Great Circle Distance

In this study, 98 regression models were specified for easily estimating shortest distances based on great circle distances along the U.S. interstate highways nationwide and for each of the continental 48 states. This allows transportation professionals to quickly generate distance, or even distance matrix, without expending significant efforts on complicated shortest path calculations. For simple usage by all professionals, all models are present in the simple linear regression form. Only one explanatory variable, the great circle distance, is considered to calculate the route distance. For each geographic scope (i.e., the national or one of the states), two different models were considered, with and without the intercept. Based on the adjusted R-squared, it was observed that models without intercepts generally have better fitness. Additionally, all these models generally have good fitness with the linear regression relationship between the great circle distance and route distance. At the state level, significant variations in the slope coefficients between the state-level models were also observed. Furthermore, a preliminary analysis of the effect of highway density on this variation was conducted.

33 ADVANCED PROPULSION SYSTEMS↗

Substitution Models of Protein Evolution with Selection on Enzymatic Activity

Abstract Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme–substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

Ferreiro, David↗

Photosynthetic responses of switchgrass to light and CO 2 under different precipitation treatments

Switchgrass ( Panicum virgatum L .) is a prominent bioenergy crop with robust resilience to environmental stresses. However, our knowledge regarding how precipitation changes affect switchgrass photosynthesis and its responses to light and CO 2 remains limited. To address this knowledge gap, we conducted a field precipitation experiment with five different treatments, including −50%, −33%, 0%, +33%, and +50% of ambient precipitation. To determine the responses of leaf photosynthesis to CO 2 concentration and light, we measured leaf net photosynthesis of switchgrass under different CO 2 concentrations and light levels in 2020 and 2021 for each of the five precipitation treatments. We first evaluated four light and CO 2 response models (i.e., rectangular hyperbola model, nonrectangular hyperbola model, exponential model, and the modified rectangular hyperbola model) using the measurements in the ambient precipitation treatment. Based on the fitting criteria, we selected the nonrectangular hyperbola model as the optimal model and applied it to all precipitation treatments, and estimated model parameters. Overall, the model fit field measurements well for the light and CO 2 response curves. Precipitation change did not influence the maximum net photosynthetic rate ( P max ) but influenced other model parameters including quantum yield ( α ), convexity ( θ ), dark respiration ( Rd ), light compensation point ( LCP ), and saturated light point ( LSP ). Specifically, the mean P max of five precipitation treatments was 17.6 μmol CO 2 m −2 s −1 , and the ambient treatment tended to have a higher P max . The +33% treatment had the highest α , and the ambient treatment had lower θ and LCP , higher Rd , and relatively lower LSP . Furthermore, precipitation significantly influenced all model parameters of CO 2 response. The ambient treatment had the highest P max , largest α , and lowest θ , R d , and CO 2 compensation point LCP . Overall, this study improved our understanding of how switchgrass leaf photosynthesis responds to diverse environmental factors, providing valuable insights for accurately modeling switchgrass ecophysiology and productivity.

09 BIOMASS FUELS↗

Hyperfine interactions of constituent quarks and the mass spectrum of tetraquark states

The simplest form of the Hyperfine Interaction (HFI) of quarks is a color-spin interaction, sometimes called Fermi-Breit (FB) HFI. For conventional hadrons, this HFI for quark-antiquark (qq) pairs is twice as strong as the same HFI for quark-quark (qq or qq) pairs, but this difference is absorbed into the different HFI constants for mesons and baryons, which are free parameters of the fits. For tetraquarks, however, qqHFI is only half as strong as the qq/qqHFI and both kinds of pairs coexist within the same state. Hence it is worth investigating a model with different fit constants for qq and qq/q pairs. We have investigated 3 models of spin-color HFI for light scalar tetraquarks: 2 models with single constants for both kinds of pairs and one model with different fit constants. The light scalar nonet is chosen because of a lack of other experimentally observed multiplets interpretable as pure exotic states. Since even the nonet contains only 3 particles with distinct masses (i. e. 3 fit points), the constituent quark masses had to be taken from conventional hadron fits obtained in our previous work. The two-constant model did not prove to be a viable alternative to the one-constant models. © 2019 Author(s).

60 APPLIED LIFE SCIENCES↗

Oscillatory dynamics of an intravenous glucose tolerance test model with delay interval

We report that type 2 diabetes mellitus (T2DM) has become prevalent pandemic disease in view of the modern life style. Both diabetic population and health expenses grow rapidly according to American Diabetes Association. Detecting the potential onset of T2DM is an essential focal point in the research of diabetes mellitus. The intravenous glucose tolerance test (IVGTT) is an effective protocol to determine the insulin sensitivity, glucose effectiveness, and pancreatic β-cell functionality, through the analysis and parameter estimation of a proper differential equation model. Delay differential equations have been used to study the complex physiological phenomena including the glucose and insulin regulations. In this paper, we propose a novel approach to model the time delay in IVGTT modeling. This novel approach uses two parameters to simulate not only both discrete time delay and distributed time delay in the past interval, but also the time delay distributed in a past sub-interval. Normally, larger time delay, either a discrete or a distributed delay, will destabilize the system. However, we find that time delay over a sub-interval might not. We present analytically some basic model properties, which are desirable biologically and mathematically. We show that this relatively simple model provides good fit to fluctuating patient data sets and reveals some intriguing dynamics. Moreover, our numerical simulation results indicate that our model may remove the defect in well known Minimal Model, which often overestimates the glucose effectiveness index. Delay differential equations (DDEs) have been frequently used to study complex dynamics observed in nature. More recently, they are used to understand intriguing physiological phenomena such as those expressed by glucose and insulin interaction. We propose a simple set of delay differential equations to model an intravenous glucose tolerance test. This model uses two parameters to simulate not only both discrete time delay and distributed time delay in the past interval, but also the time delay distributed in a past sub-interval. We show that this relatively simple model provides good fit to fluctuating patient data sets and reveals some intriguing dynamics. Lastly and most importantly, our model may remove the defect in the well known Minimal Model (MM) which often overestimates the glucose effectiveness (GE) index.

59 BASIC BIOLOGICAL SCIENCES↗

The Halo Mass of Optically Luminous Quasars at z ≈ 1–2 Measured via Gravitational Deflection of the Cosmic Microwave Background

We measure the average deflection of cosmic microwave background photons by quasars at $\langle z\rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range 0.9 ≤ z ≤ 2.2 with absolute i-band magnitudes of M i ≤ –24 (K-corrected to z = 2). A stack of nearly 200,000 targets reveals an 8σ detection of Planck's estimate of the lensing convergence toward the quasars. We fit the signal with a model comprising a Navarro–Frenk–White density profile and a two-halo term accounting for correlated large-scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass ${\mathrm{log}}_{10}({M}_{{\rm{h}}}/{h}^{-1}\,{M}_{\odot })=12.6\pm 0.2$ and linear bias b = 2.7 ± 0.3 at $\langle z\rangle =1.7$, in excellent agreement with clustering studies. Here, we also report a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than M i lesssim –26 reside in halos of typical mass ${M}_{{\rm{h}}}\approx {10}^{13}\,{h}^{-1}\,{M}_{\odot }$, scaling roughly as ${M}_{{\rm{h}}}\propto {L}_{\mathrm{opt}}^{3/4}$ at ${M}_{i}\lesssim -24$ mag, in good agreement with physically motivated quasar demography models. Although we acknowledge that this luminosity dependence is a marginal result, the observed M h –L opt relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates toward high Eddington ratios: the weak trend of M h with L opt observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit.

79 ASTRONOMY AND ASTROPHYSICS↗

An efficient parallel sampling technique for Multivariate Poisson-Lognormal model: Analysis with two crash count datasets

Our study investigates the Multivariate Poisson-lognormal (MVPLN) model that jointly models crash frequency and severity accounting for correlations. The ordinary univariate count models analyze crashes of different severity level separately ignoring the correlations among severity levels. The MVPLN model is capable to incorporate the general correlation structure and takes account of the over dispersion in the data that leads to a superior data fitting. But, the traditional estimation approach for MVPLN model is computationally expensive, which often limits the use of MVPLN model in practice. In this work, a parallel sampling scheme is introduced to improve the original Markov Chain Monte Carlo (MCMC) estimation approach of the MVPLN model, which significantly reduces the model estimation time. Two MVPLN models are developed using the pedestrian vehicle crash data collected in New York City from 2002 to 2006, and the highway-injury data from Washington State (5-year data from 1990 to 1994) The Deviance Information Criteria (DIC) is used to evaluate the model fitting. The estimation results show that the MVPLN models provide a superior fit over univariate Poisson-lognormal (PLN), univariate Poisson, and Negative Binomial models. Moreover, the correlations among the latent effects of different severity levels are found significant in both datasets that justifies the importance of jointly modeling crash frequency and severity accounting for correlations.

97 MATHEMATICS AND COMPUTING↗

Radio frequency models of novae in eruption. I. The free-free process in bipolar morphologies

Observations of novae at radio frequencies provide us with a measure of the total ejected mass, density profile, and kinetic energy of a nova eruption. The radio emission is typically well characterized by the free-free emission process. Most models to date have assumed spherical symmetry for the eruption, although for as long as there have been radio observations of these systems, it has been known that spherical eruptions are too simplistic a geometry. In this paper, we build bipolar models of the nova eruption, assuming the free-free process, and show the effects of varying different parameters on the radio light curves. The parameters considered include the ratio of the minor- to major-axis, the inclination angle, and shell thickness. We also show the uncertainty introduced when fitting spherical-model synthetic light curves to bipolar-model synthetic light curves. We find that the optically thick phase rises with the same power law (S {sub ν}∝t {sup 2}) for both the spherical and bipolar models. In the bipolar case, there is a 'plateau' phase—depending on the thickness of the shell as well as the ratio of the minor- to major-axis—before the final decline, which follows the same power law (S {sub ν}∝t {sup –3}) as in the spherical case. Finally, fitting spherical models to the bipolar-model synthetic light curves requires, in the worst-case scenario, doubling the ejected mass, more than halving the electron temperature, and reducing the shell thickness by nearly a factor of 10. This implies that in some systems we have been over-predicting the ejected masses and under-predicting the electron temperature of the ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

Contemporaneous broadband observations of three high-redshift BL Lac objects

We have collected broadband spectral energy distributions (SEDs) of three BL Lac objects 3FGL J0022.1−1855 (z = 0.689), 3FGL J0630.9−2406 (z ≳ 1.239), and 3FGL J0811.2−7529 (z = 0.774), detected by Fermi with relatively flat gigaelectronvolt spectra. By observing simultaneously in the near-infrared to hard X-ray band, we can well characterize the high end of the synchrotron component of the SED. Thus, fitting the SEDs to synchro-Compton models of the dominant emission from the relativistic jet, we can constrain the underlying particle properties and predict the shape of the gigaelectronvolt Compton component. Standard extragalactic background light (EBL) models explain the high-energy absorption well, with poorer fits for high-ultraviolet models. The fits show clear evidence for EBL absorption in the Fermi spectrum of our highest-redshift source 3FGL J0630.9−2406. While synchrotron self-Compton models adequately describe the SEDs, the situation may be complicated by possible external Compton components. For 3FGL J0811.2−7529, we also discover a nearby serendipitous source in the X-ray data, which is almost certainly another lower synchrotron peak frequency (ν{sub pk}{sup sy}) BL Lac, that may contribute flux in the Fermi band. Since our sources are unusual high-luminosity, moderate ν{sub pk}{sup sy} BL Lacs, we compare these quantities and the Compton dominance, the ratio of peak inverse Compton to peak synchrotron luminosities (L{sub pk}{sup IC}/L{sub pk}{sup sy}), with those of the full Fermi BL Lac population.

79 ASTRONOMY AND ASTROPHYSICS↗

The Thermal and Kinematic Sunyaev–Zeldovich Effect in Galaxy Clusters and Filaments Using Multifrequency Temperature Maps of the Cosmic Microwave Background: A399–A401 Cluster Pair Case Study

We present a multifrequency and multi-instrument methodology to study the physical properties of galaxy clusters and cosmic filaments using cosmic microwave background observations. Our approach enables simultaneous measurement of both the thermal (tSZ) and kinematic Sunyaev–Zeldovich (kSZ) effects, incorporates relativistic corrections, and models astrophysical foregrounds such as thermal dust emission. We do this by jointly fitting a single physical model across multiple maps from multiple instruments at different frequencies, rather than fitting a model to a single Compton-y map. We demonstrate the success of this method by fitting the A399–A401 galaxy cluster pair and filament system using archival data from the Planck satellite and new, targeted deep data from the Atacama Cosmology Telescope, covering 11 different frequencies over 14 maps from 30 GHz to 545 GHz. Our tSZ results are consistent with previous work using Compton-y maps. We measure the line-of-sight peculiar velocities of the cluster–filament system using the kSZ effect and find statistical uncertainties on individual cluster peculiar velocities of ≲600 km s −1 , which are competitive with current state-of-the-art measurements. Additionally, we measure the optical depth of the filament component with a signal-to-noise of 8.5σ and reveal hints of its morphology. This modular approach is well-suited for application to future instruments across a wide range of millimeter and submillimeter wavebands.

Gill, Ajay S. [National Research Council of Canada↗

Superhorizon isocurvature fluctuations relax tensions

Here, we present a new class of models that have potential to alleviate tensions present in the cosmological data today. We postulate about variation in the sound horizon scale on superhorizon scales, i.e., on scales that are larger than that of the present observable low-redshift Universe (≳ 1Gpc) while at the same time smaller than the largest scales probed by the cosmic microwave background (CMB) (≲ 10Gpc). In this scenario, CMB peaks are naturally smoothed as preferred by the Planck data, while at the same time the low-redshift baryon acoustic oscillation calibration is partially decoupled from the CMB. Taking superhorizon variations in baryon fraction as an example and using approximate modeling, we find improvement in the best fit Planck power spectrum model Δχ 2 ~ 6 for 1 extra degree of freedom with the relevant extension parameter 10 3 σ b = 2.10 ± 0.60, implying about 10% variations in baryon fraction across the Universe. At the same time, S 8 drops by about 1 sigma, easing tension with weak lensing surveys. We find that the combination of Planck 2018 data, eBOSS BAO data, and Riess et al. distance ladder Hubble parameter determination produce a satisfactory fit in our model if we allow for a phantom dark energy equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗