Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HISTORY”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Effects of large-scale environment on the assembly history of central galaxies

We examine whether large-scale environment affects the mass assembly history of central galaxies. To facilitate this, we constructed dark matter halo merger trees from a cosmological N-body simulation and calculated the formation and evolution of galaxies using a semi-analytic method. We confirm earlier results that smaller halos show a notable difference in formation time with a mild dependence on large-scale environment. However, using a semi-analytic model, we found that on average the growth rate of the stellar mass of central galaxies is largely insensitive to large-scale environment. Although our results show that the star formation rate (SFR) and the stellar mass of central galaxies in smaller halos are slightly affected by the assembly bias of halos, those galaxies are faint and the difference in the SFR is minute, therefore it is challenging to detect it in real galaxies given the current observational accuracy. Future galaxy surveys, such as the BigBOSS experiment and the Large Synoptic Survey Telescope, which are expected to provide observational data for fainter objects, will provide a chance to test our model predictions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Delayed star formation in isolated dwarf galaxies: Hubble space telescope star formation history of the Aquarius dwarf irregular

We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope Advanced Camera for Surveys. The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ≈10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ≈10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ≈ 2). The star formation rate increased dramatically ≈6-8 Gyr ago (z ≈ 0.7-1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M {sub H} {sub I}/M {sub *}, dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degrees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A uniform history for galaxy evolution

Recent observations indicate a remarkable similarity in the properties of evolving galaxies at fixed mass and redshift, prompting us to consider the possibility that most galaxies may evolve with a common history encompassing star formation, quasar accretion, and eventual quiescence. We quantify this by defining a synchronization timescale for galaxies as a function of mass and redshift that characterizes the extent to which different galaxies of a common mass are evolving in the same manner at various cosmic epochs. We measure this synchronization timescale using nine different star-forming galaxy observations from the literature and Sloan Digital Sky Survey quasar observations spanning 0 < z ≲ 6. Surprisingly, this synchronization timescale is a constant, approximately 1.5 Gyr for all combinations of mass and time. We also find that the ratio between the stellar mass of galaxies turning off star formation and black hole mass of turnoff quasars is approximately 30:1, much lower than the 500:1 for quiescent galaxies at low redshift. As a result, we propose a model in which the star-forming 'main sequence', analogous quasar behavior, and other observations form a galactic evolution 'main sequence', in which star formation occurs earliest, followed by supermassive black hole accretion, and feedback between the two are dominated by deterministic rather than stochastic processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A TALE OF A RICH CLUSTER AT z ∼ 0.8 AS SEEN BY THE STAR FORMATION HISTORIES OF ITS EARLY-TYPE GALAXIES

We present a detailed stellar population analysis for a sample of 24 early-type galaxies (ETGs) belonging to the rich cluster RX J0152.7-1357 at z = 0.83. We have derived the age, metallicity, abundance pattern, and star formation history (SFH) for each galaxy individually to further characterize this intermediate-z reference cluster. We then study how these stellar population parameters depend on the local environment. This provides a better understanding on the formation timescales and subsequent evolution of the substructures in this cluster. We have also explored the evolutionary link between z ∼ 0.8 ETGs and those in the local universe by comparing the trends that the stellar population parameters followed with galaxy velocity dispersion at each epoch. We find that the ETGs in Coma are consistent with being the (passively evolving) descendants of the ETG population in RX J10152.7-1357. Furthermore, our results favor a downsizing picture, where the subclumps centers were formed first. These central parts contain the most massive galaxies, which formed the bulk of their stars in a short, burst-like event at high z. On the contrary, the cluster outskirts are populated with less-massive, smaller galaxies that show a wider variety of SFHs. In general, they present extended star formation episodes over cosmic time, which seems to be related to their posterior incorporation into the cluster around 4 Gyr after the initial event of formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE STAR-FORMING HISTORIES OF THE NUCLEUS, BULGE, AND INNER DISK OF NGC 5102: CLUES TO THE EVOLUTION OF A NEARBY LENTICULAR GALAXY {sup ,} {sup ,}

Long slit spectra recorded with the Gemini Multi-Object Spectrograph on Gemini South are used to examine the star-forming history (SFH) of the lenticular galaxy NGC 5102. Structural and supplemental photometric information are obtained from archival Spitzer [3.6] images. Absorption features at blue and visible wavelengths are traced out along the minor axis to galactocentric radii ∼60 arcsec (∼0.9 kpc), sampling the nucleus, bulge, and disk components. Comparisons with model spectra point to luminosity-weighted metallicities that are consistent with the colors of resolved red giant branch stars in the disk. The nucleus has a luminosity-weighted age at visible wavelengths of ∼1{sub −0.1}{sup +0.2} Gyr, and the integrated light is dominated by stars that formed over a time period of only a few hundred Myr. For comparison, the luminosity-weighted ages of the bulge and disk are ∼2{sub −0.2}{sup +0.5} Gyr and 10{sub −2}{sup +2} Gyr, respectively. The g' – [3.6] colors of the nucleus and bulge are consistent with the spectroscopically based ages. In contrast to the nucleus, models that assume star-forming activity spanning many Gyr provide a better match to the spectra of the bulge and disk than simple stellar population models. Isophotes in the bulge have a disky shape, hinting that the bulge was assembled from material with significant rotational support. The SFHs of the bulge and disk are consistent with the bulge forming from the collapse of a long-lived bar, rather than from the collapse of a transient structure that formed as the result of a tidal interaction. It is thus suggested that the progenitor of NGC 5102 was a barred disk galaxy that morphed into a lenticular galaxy through the buckling of its bar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A NEW CHEMICAL EVOLUTION MODEL FOR DWARF SPHEROIDAL GALAXIES BASED ON OBSERVED LONG STAR FORMATION HISTORIES

We present a new chemical evolution model for dwarf spheroidal galaxies (dSphs) in the local universe. Our main aim is to explain both their observed star formation histories and metallicity distribution functions simultaneously. Applying our new model for the four local dSphs, that is, Fornax, Sculptor, Leo II, and Sextans, we find that our new model reproduces the observed chemical properties of the dSphs consistently. Our results show that the dSphs have evolved with both a low star formation efficiency and a large gas outflow efficiency compared with the Milky Way, as suggested by previous works. Comparing the observed [α/Fe]-[Fe/H] relation of the dSphs with the model predictions, we find that our model favors a longer onset time of Type Ia supernovae (i.e., 0.5 Gyr) than that suggested in previous studies (i.e., 0.1 Gyr). We discuss the origin of this discrepancy in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

P-MaNGA: GRADIENTS IN RECENT STAR FORMATION HISTORIES AS DIAGNOSTICS FOR GALAXY GROWTH AND DEATH

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000 Å break (D{sub n}(4000)), Hδ absorption (EW(Hδ{sub A})), and Hα emission (EW(Hα)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D{sub n}(4000) measured in the central spaxel of each datacube exceeds 1.6. For each spaxel we generate both 2D maps and radial profiles of D{sub n}(4000), EW(Hδ{sub A}), and EW(Hα). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D{sub n}(4000) decreases, while both EW(Hδ{sub A}) and EW(Hα) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outward as it moves to the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The star formation histories of local group dwarf galaxies. III. Characterizing quenching in low-mass galaxies

We explore the quenching of low-mass galaxies (10{sup 4} ≲ M{sub ⋆} ≲ 10{sup 8} M{sub ⊙}) as a function of lookback time using the star formation histories (SFHs) of 38 Local Group dwarf galaxies. The SFHs were derived by analyzing color–magnitude diagrams of resolved stellar populations in archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. We find: (1) lower-mass galaxies quench earlier than higher-mass galaxies; (2) inside of R{sub virial} there is no correlation between a satellite’s current proximity to a massive host and its quenching epoch; and (3) there are hints of systematic differences in the quenching times of M31 and Milky Way (MW) satellites, although the sample size and uncertainties in the SFHs of M31 dwarfs prohibit definitive conclusions. Combined with results from the literature, we qualitatively consider the redshift evolution (z = 0–1) of the quenched galaxy fraction over ∼7 dex in stellar mass (10{sup 4} ≲ M{sub ⋆} ≲ 10{sup 11.5} M{sub ⊙}). The quenched fraction of all galaxies generally increases toward the present, with both the lowest and highest-mass systems exhibiting the largest quenched fractions at all redshifts. In contrast, galaxies between M{sub ⋆} ∼ 10{sup 8}–10{sup 10} M{sub ⊙} have the lowest quenched fractions. We suggest that such intermediate-mass galaxies are the least efficient at quenching. Finally, we compare our quenching times with predictions for infall times for low-mass galaxies associated with the MW. We find that some of the lowest-mass satellites (e.g., CVn II, Leo IV) may have been quenched before infall, while higher-mass satellites (e.g., Leo I, Fornax) typically quench ∼1–4 Gyr after infall.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The panchromatic Hubble Andromeda Treasury. XI. The spatially resolved recent star formation history of M31

We measure the recent star formation history (SFH) across M31 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color–magnitude diagrams in ∼9000 regions that are ∼100 pc × 100 pc in projected size, covering a 0.5 square degree area (∼380 kpc{sup 2}, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ∼500 Myr, most notably in M31's 10 kpc star-forming ring. We find that the 10 kpc ring is at least 400 Myr old, showing ongoing star formation (SF) over the past ∼500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring’s position is constant throughout this time, and is stationary at the level of 1 km s{sup −1}, although there is evidence for broadening of the ring due to the diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring’s position makes a purely collisional ring origin highly unlikely. Besides the well-known 10 kpc ring, we observe two other ring-like features. There is an outer ring structure at 15 kpc with concentrated SF starting ∼80 Myr ago. The inner ring structure at 5 kpc has a much lower star formation rate (SFR) and therefore lower contrast against the underlying stellar disk. It was most clearly defined ∼200 Myr ago, but is much more diffuse today. We find that the global SFR has been fairly constant over the last ∼500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ∼500 Myr, ∼60% of all SF has occurred in the 10 kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is 0.28 ± 0.03 M{sub ⊙} yr{sup −1} with an average deprojected intensity of 7.3×10{sup −4} M{sub ⊙} yr{sup −1} kpc{sup −2}, which yields a total SFR of ∼0.7 M{sub ⊙} yr{sup −1} when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

GAS INFLOW AND OUTFLOW HISTORIES IN DISK GALAXIES AS REVEALED FROM OBSERVATIONS OF DISTANT STAR-FORMING GALAXIES

We investigate gas inflow and outflow histories in Milky Way-like disk galaxies, to get new insights into the baryonic processes in galaxy formation and evolution. For this purpose, we solve the equations for the evolution of the surface mass densities of gas and metals at each radius in a galactic disk, based on the observed structural properties of distant star-forming galaxies, including the redshift evolution of their stellar mass distribution, their scaling relation between the mass of baryonic components, star formation rate (SFR), and chemical abundance, as well as the supposed evolution of their radial metallicity gradients (RMGs). We find that the efficiency of gas inflow for a given SFR decreases with time and that the inflow rate is always nearly proportional to the SFR. For gas outflow, although its efficiency for a given SFR is a decreasing function of time, similar to gas inflow, the outflow rate is not necessarily proportional to the SFR and the relation between the outflow rate and SFR strongly depends on the evolution of the adopted RMG. We also find that the results on the outflow rate can be reproduced in the framework of a momentum-driven (energy-driven) wind mechanism if the RMG is steepening (flattening) with time. Therefore if the well-measured RMGs and their evolution for Milky Way-like galaxies are obtained from future observations, then our results will be useful to constrain the main driving mechanism for their galactic outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PATCHY BLAZAR HEATING: DIVERSIFYING THE THERMAL HISTORY OF THE INTERGALACTIC MEDIUM

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron–positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature–density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large-scale simulations to determine the impact of inhomogeneous heating. We show that because of blazar clustering, blazar heating is inhomogeneous for z ≳ 2. At high redshift, the temperature–density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature–density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large-scale simulations and provides a cost-effective alternative to subgrid models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Hubble Tarantula Treasury Project. II. The star-formation history of the starburst region NGC 2070 in 30 Doradus

We present a study of the recent star formation (SF) of 30 Doradus in the Large Magellanic Cloud (LMC) using the panchromatic imaging survey Hubble Tarantula Treasury Project. In this paper we focus on the stars within 20 pc of the center of 30 Doradus, the starburst region NGC 2070. We recovered the SF history by comparing deep optical and near-infrared color–magnitude diagrams (CMDs) with state-of-the-art synthetic CMDs generated with the latest PAdova and TRieste Stellar Evolution Code (PARSEC) models, which include all stellar phases from pre-main-sequence to post-main-sequence. For the first time in this region we are able to measure the SF using intermediate- and low-mass stars simultaneously. Our results suggest that NGC 2070 experienced prolonged activity. In particular, we find that the SF in the region (1) exceeded the average LMC rate ≈ 20 Myr ago, (2) accelerated dramatically ≈ 7 Myr ago, and (3) reached a peak value 1–3 Myr ago. We did not find significant deviations from a Kroupa initial mass function down to 0.5 M{sub ⊙}. The average internal reddening E(B−V) is found to be between 0.3 and 0.4 mag.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

YOUNG STELLAR OBJECTS IN LYNDS 1641: DISKS, ACCRETION, AND STAR FORMATION HISTORY

We investigate the young stellar objects (YSOs) in the Lynds 1641 (L1641) cloud using multi-wavelength data including Spitzer, WISE, the Two Micron All Sky Survey, and XMM covering {approx}1390 YSOs across a range of evolutionary stages. In addition, we targeted a sub-sample of YSOs for optical spectroscopy with the MMT/Hectospec and the MMT/Hectochelle. We use these data, along with archival photometric data, to derive spectral types, extinction values, masses, ages, and accretion rates. We obtain a disk fraction of {approx}50% in L1641. The disk frequency is almost constant as a function of stellar mass with a slight peak at log (M{sub *}/M{sub Sun }) Almost-Equal-To -0.25. The analysis of multi-epoch spectroscopic data indicates that the accretion variability of YSOs cannot explain the two orders of magnitude of scatter for YSOs with similar masses. Forty-six new transition disk (TD) objects are confirmed in this work, and we find that the fraction of accreting TDs is lower than for optically thick disks (40%-45% versus 77%-79%, respectively). We confirm our previous result that the accreting TDs have a median accretion rate similar to normal optically thick disks. We confirm that two star formation modes (isolated versus clustered) exist in L1641. We find that the diskless YSOs are statistically older than the YSOs with optically thick disks and the TD objects have a median age that is intermediate between those of the other two populations. We tentatively study the star formation history in L1641 based on the age distribution and find that star formation started to be active 2-3 Myr ago.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

STAR FORMATION HISTORY OF THE MILKY WAY HALO TRACED BY THE OOSTERHOFF DICHOTOMY AMONG GLOBULAR CLUSTERS

In our recent investigation of the Oosterhoff dichotomy in the multiple population paradigm, we have suggested that the RR Lyrae variables in the globular clusters (GCs) of Oosterhoff groups I, II, and III are produced mostly by first, second, and third generation stars (G1, G2, and G3), respectively. Here we show, for the first time, that the observed dichotomies in the inner and outer halo GCs can be naturally reproduced when these models are extended to all metallicity regimes, while maintaining reasonable agreements in the horizontal-branch type versus [Fe/H] correlations. In order to achieve this, however, specific star formation histories are required for the inner and outer halos. In the inner halo GCs, the star formation commenced and ceased earlier with a relatively short formation timescale between the subpopulations (∼0.5 Gyr), while in the outer halo, the formation of G1 was delayed by ∼0.8 Gyr with a more extended timescale between G1 and G2 (∼1.4 Gyr). This is consistent with the dual origin of the Milky Way halo. Despite the difference in detail, our models show that the Oosterhoff period groups observed in both outer and inner halo GCs are all manifestations of the “population-shift” effect within the instability strip, for which the origin can be traced back to the two or three discrete episodes of star formation in GCs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Dark before light: testing the cosmic expansion history through the cosmic microwave background

The cosmic expansion history proceeds in broad terms from a radiation dominated epoch to matter domination to an accelerated, dark energy dominated epoch. We investigate whether intermittent periods of acceleration (from a canonical, minimally coupled scalar field) are possible in the early universe — between Big Bang nucleosynthesis (BBN) and recombination and beyond. We establish that the standard picture is remarkably robust: anisotropies in the cosmic microwave background consistent with ΛCDM will exclude any extra period of accelerated expansion between 1 ≤ z∼<10{sup 5} (corresponding to 5 × 10{sup −4}eV ≤ T∼<25eV)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Cosmic histories of star formation and reionization: an analysis with a power-law approximation

With a simple power-law approximation of high-redshift (∼>3.5) star formation history, i.e., ρ-dot {sub *}(z)∝[(1+z)/4.5]{sup −α}, we investigate the reionization of intergalactic medium (IGM) and the consequent Thomson scattering optical depth for cosmic microwave background (CMB) photons. A constraint on the evolution index α is derived from the CMB optical depth measured by the Wilkinson Microwave Anisotropy Probe (WMAP) experiment, which reads α ≈ 2.18 lg N{sub γ}−3.89, where the free parameter N{sub γ} is the number of the escaped ionizing ultraviolet photons per baryon. At the same time, the redshift z{sub f} at which the IGM is fully ionized can also be expressed as a function of α as well as N{sub γ}. By further taking into account the implication of the Gunn-Peterson trough observations to quasars for the full reionization redshift, i.e., 6∼<z{sub f}∼<7, we obtain 0.3∼<α∼<1.3 and 80∼<N{sub γ}∼<230. For a typical number of ∼ 4000 of ionizing photons released per baryon of normal stars, the fraction of these photons escaping from the stars, f{sub esc}, can be constrained to within the range of (2.0–5.8)%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Neutrinos in the holographic dark energy model: constraints from latest measurements of expansion history and growth of structure

The model of holographic dark energy (HDE) with massive neutrinos and/or dark radiation is investigated in detail. The background and perturbation evolutions in the HDE model are calculated. We employ the PPF approach to overcome the gravity instability difficulty (perturbation divergence of dark energy) led by the equation-of-state parameter w evolving across the phantom divide w=−1 in the HDE model with c<1. We thus derive the evolutions of density perturbations of various components and metric fluctuations in the HDE model. The impacts of massive neutrino and dark radiation on the CMB anisotropy power spectrum and the matter power spectrum in the HDE scenario are discussed. Furthermore, we constrain the models of HDE with massive neutrinos and/or dark radiation by using the latest measurements of expansion history and growth of structure, including the Planck CMB temperature data, the baryon acoustic oscillation data, the JLA supernova data, the Hubble constant direct measurement, the cosmic shear data of weak lensing, the Planck CMB lensing data, and the redshift space distortions data. We find that ∑ m{sub ν}<0.186 eV (95% CL) and N{sub eff}=3.75{sup +0.28}{sub −0.32} in the HDE model from the constraints of these data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A general reconstruction of the recent expansion history of the universe

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary (n+1)-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and H(z) data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. We perform the MC analysis generating mock catalogs and computing their best-fit. For each fiducial model, we test different reconstructions using, in each case, more than 10{sup 4} catalogs in a total of about 5× 10{sup 5}. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most 10% of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit H{sub 0} and the sound horizon r{sub s}(z{sub d}) at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift 1.3 (using only JLA) and 2.3 (JLA+BAO+H(z)). We obtain that the standard cosmological model agrees within 3σ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗