Laboratoire Univers et Théories

The study of stellar populations in galaxies is entering a new era with the availability of large and high-quality data bases of both observed galactic spectra and state-of-the-art evolutionary synthesis models. In this paper we investigate the power of spectral synthesis as a means to estimate the physical properties of galaxies. Spectral synthesis is nothing more than the decomposition of an observed spectrum in terms of a superposition of a base of simple stellar populations of various ages and metallicities, producing as output the star formation and chemical histories of a galaxy, its extinction and velocity dispersion. Our implementation of this method uses the recent models of Bruzual & Charlot and observed spectra in the 3650-8000 range. The reliability of this approach is studied by three different means: (1) simulations, (2) comparison with previous work based on a different technique, and (3) analysis of the consistency of results obtained for a sample of galaxies from the Sloan Digital Sky Survey (SDSS).

What are the thermodynamic properties of matter at extreme densities, even exceeding nuclear matter density severely? How can we describe the composition of matter for such conditions, the resulting pressure, and the maximum mass of cold neutron stars? How is this affected by finite temperatures, as they occur in core collapse supernovae and in compact star mergers? This review addresses these points within the framework of constraints from experiments as well as astronomical observations.

We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). A robust GWB detection is conditioned upon resolving the Hellings-Downs angular pattern in the pairwise cross-correlation of the pulsar timing residuals. Additionally, the GWB is expected to yield the same (common) spectrum of temporal correlations across pulsars, which is used as a null hypothesis in the GWB search. Such a common-spectrum process has already been observed in pulsar timing data. We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of 4%. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of 60 and a false alarm probability of about 0.1% (≳3 σ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The latest EPTA data from new generation observing systems show non-negligible evidence for the GWB. At the same time, the inferred spectrum is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of (2.5 ± 0.7) × 10 −15 at a reference frequency of 1 yr −1 . Further investigation of these issues is required for reliable astrophysical interpretations of this signal. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years.

We have re-evaluated empirical expressions for the abundance determination of N, O, Ne, S, Cl, Ar and Fe taking into account the latest atomic data and constructing an appropriate grid of photoionization models with state-of-the art model atmospheres. Using these expressions we have derived heavy element abundances in the ~310 emission-line galaxies from the Data Release 3 of the Sloan Digital Sky Survey (SDSS) with an observed Hβ flux F(Hβ) > 10-14 erg s-1 cm-2 and for which the [O iii] λ4363 emission line was detected at least at a 2σ level, allowing abundance determination by direct methods. The oxygen abundance 12 + log O/H of the SDSS galaxies lies in the range from ~7.1 (/30) to ~8.5 (0.7 ). The SDSS sample is merged with a sample of 109 blue compact dwarf (BCD) galaxies with high quality spectra, which contains extremely low-metallicity objects. We use the merged sample to study the abundance patterns of low-metallicity emission-line galaxies. We find that extremely metal-poor galaxies (12 + log O/H < 7.6, i.e. Z < /12) are rare in the SDSS sample. The α element-to-oxygen abundance ratios do not show any significant trends with oxygen abundance, in agreement with previous studies, except for a slight increase of Ne/O with increasing metallicity, which we interpret as due to a moderate depletion of O onto grains in the most metal-rich galaxies. The Fe/O abundance ratio is smaller than the solar value, by up to 1 dex at the high metallicity end. We also find that Fe/O increases with decreasing Hβ equivalent width EW(Hβ). We interpret this as a sign of strong depletion onto dust grains, and gradual destruction of those grains on a time scale of a few Myr. All the galaxies are found to have log N/O > –1.6, implying that they have a different nature than the subsample of high-redshift damped Lyα systems with log N/O of ~–2.3 and that their ages are larger than 100–300 Myr. We confirm the apparent increase in N/O with decreasing EW(Hβ), already shown in previous studies, and explain it as the signature of gradual nitrogen ejection by massive stars from the most recent starburst.

International audience

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar data set spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar system ephemeris errors, obtaining a robust 95 per cent upper limit on the dimensionless strain amplitude A of the background of A < 3.0 10 -15 at a reference frequency of 1 yr -1 and a spectral index of 13/3, corresponding to a background from inspiralling supermassive black hole binaries, constraining the GW energy density to gw (f)h 2 < 1.1 10 -9 at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of 5 10 -9 Hz. Finally, we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95 per cent upper limits on the string tension, G/c 2 , characterizing a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit G/c 2 < 1.3 10 -7 , identical to that set by the Planck Collaboration, when combining Planck and high-cosmic microwave background data from other experiments. For a stochastic relic background, we set a limit of relic gw (f )h 2 < 1.2 10 -9 , a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

We present the revised ``Meudon'' model of Photon Dominated Region (PDR code), presently available on the web under the Gnu Public Licence at: http://aristote.obspm.fr/MIS. General organisation of the code is described down to a level that should allow most observers to use it as an interpretation tool with minimal help from our part. Two grids of models, one for low excitation diffuse clouds and one for dense highly illuminated clouds, are discussed, and some new results on PDR modelisation highlighted.

A numerous population of weak line galaxies (WLGs) is often left out of statistical studies on emission-line galaxies (ELGs) due to the absence of an adequate classification scheme, since classical diagnostic diagrams, such as [O &lt;scp&gt;iii&lt;/scp&gt;]/Hβ versus [N &lt;scp&gt;ii&lt;/scp&gt;]/Hα (the BPT diagram), require the measurement of at least four emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between star-forming (SF) galaxies and active galactic nuclei (AGN) hosts and between Seyferts and LINERs to diagrams that are more economical in terms of line quality requirements. By doing this, we rescue from the classification limbo a substantial number of sources and modify the global census of ELGs. More specifically, (1) we use the Sloan Digital Sky Survey Data Release 7 to constitute a suitable sample of 280 000 ELGs, one-third of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al., Kauffmann et al. and Stasi&amp;nacute;ska et al. to distinguish SF galaxies and AGN hosts and Kewley et al. to distinguish Seyferts from LINERs. (3) We transpose these classification schemes to alternative diagrams keeping [N &lt;scp&gt;ii&lt;/scp&gt;]/Hα as a horizontal axis, but replacing Hβ by a stronger line (Hα or [O &lt;scp&gt;ii&lt;/scp&gt;]), or substituting the ionization-level sensitive [O &lt;scp&gt;iii&lt;/scp&gt;]/Hβ ratio with the equivalent width of Hα(&lt;it&gt;W&lt;/it&gt;&lt;inf&gt;Hα&lt;/inf&gt;). Optimized equations for the transposed divisory lines are provided. (4) We show that nothing significant is lost in the translation, but that the new diagrams allow one to classify up to 50 per cent more ELGs. (5) Introducing WLGs in the census of galaxies in the local Universe increases the proportion of metal-rich SF galaxies and especially LINERs. In the course of this analysis, we were led to make the following points. (i) The Kewley et al. BPT line for galaxy classification is generally ill-used. (ii) Replacing [O &lt;scp&gt;iii&lt;/scp&gt;]/Hβ by &lt;it&gt;W&lt;/it&gt;&lt;inf&gt;Hα&lt;/inf&gt; in the classification introduces a change in the philosophy of the distinction between LINERs and Seyferts, but not in its results. Because the &lt;it&gt;W&lt;/it&gt;&lt;inf&gt;Hα&lt;/inf&gt; versus [N &lt;scp&gt;ii&lt;/scp&gt;]/Hα diagram can be applied to the largest sample of ELGs without loss of discriminating power between Seyferts and LINERs, we recommend its use in further studies. (iii) The dichotomy between Seyferts and LINERs is washed out by WLGs in the BPT plane, but it subsists in other diagnostic diagrams. This suggests that the right wing in the BPT diagram is indeed populated by at least two classes, tentatively identified with bona fide AGN and ‘retired’ galaxies that have stopped forming stars and are ionized by their old stellar populations.

We report on the high-precision timing of 42 radiomillisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 yr. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TEMPONEST yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semimajor axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler. However, we measure an average uncertainty of 80 per cent (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binaryMSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600-3053 and J1918-0642, implying pulsar and companion masses m(p) = 1.22(-0.35)(+0.5) M-circle dot, m(c) = 0.21(-0.04)(+0.06) M-circle dot and m(p) = 1.25(-0.4)(+0.6)M(circle dot), m(c) = 0.23(-0.05)(+0.07) M-circle dot, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012+5307 and J1909-3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600-3053 and J1909-3744.

The highly stable spin of neutron stars can be exploited for a variety of (astro)physical investigations. In particular, arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such 'pulsar timing arrays' (PTAs) have been set up around the world over the past decades and collectively form the 'International' PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015–2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14 000 deg 2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends, spherical-overdensity and phase-space-based algorithms. We

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Using the code autostructure, extensive calculations of inner-shell atomic data have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni. The results are used to obtain updated opacities from the Opacity Project (OP). A number of other improvements on earlier work have also been included. Rosseland-mean opacities from the OP are compared with those from OPAL. Differences of 5–10 per cent occur. The OP gives the ‘Z-bump’, at log(T) ≃ 5.2, to be shifted to slightly higher temperatures. The opacities from the OP, as functions of temperature and density, are smoother than those from OPAL. The accuracy of the integrations used to obtain mean opacities can depend on the frequency mesh used. Tests involving variation of the numbers of frequency points show that for typical chemical mixtures the OP integrations are numerically correct to within 0.1 per cent. The accuracy of the interpolations used to obtain mean opacities for any required values of temperature and density depends on the temperature–density meshes used. Extensive tests show that, for all cases of practical interest, the OP interpolations give results correct to better than 1 per cent. Prior to a number of recent investigations which have indicated a need for downward revisions in the solar abundances of oxygen and other elements, there was good agreement between properties of the Sun deduced from helioseismology and from stellar evolution models calculated using OPAL opacities. The revisions destroy that agreement. In a recent paper, Bahcall et al. argue that the agreement would be restored if opacities for the regions of the Sun with 2 × 106≲T≲ 5 × 106 K (0.7–0.4 R⊙) were larger than those given by OPAL by about 10 per cent. In the region concerned, the present results from the OP do not differ from those of OPAL by more than 2.5 per cent.

Abstract The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where each detector has a 'xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple 'metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives.

The Planck and Herschel missions are currently measuring the far-infrared to millimeter emission of dust, which combined with existing IR data, will for the first time provide the full spectral energy distribution (SED) of the galactic interstellar medium dust emission, from the mid-IR to the mm range, with an unprecedented sensitivity and down to spatial scales 30 . Such a global SED will allow a systematic study of the dust evolution processes (e.g. grain growth or fragmentation) that directly affect the SED because they redistribute the dust mass among the observed grain sizes. The dust SED is also affected by variations of the radiation field intensity. Here we present a versatile numerical tool, DustEM, that predicts the emission and extinction of dust grains given their size distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED that will allow us to study dust evolution. We present a coherent set of observations for the DHGL SED, which has been obtained by correlating the IR and HI-21 cm data. The dust components in our DHGL model are (i) polycyclic aromatic hydrocarbons; (ii) amorphous carbon and (iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains (a < 10 nm) relative to the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck (and similarly to Herschel) channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is publically available.

This paper considers the techniques to distinguish normal star forming (NSF) galaxies and active galactic nuclei (AGN) hosts using optical spectra. The observational data base is a set of 20 000 galaxies extracted from the Sloan Digital Sky Survey, for which we have determined the emission line intensities after subtracting the stellar continuum obtained from spectral synthesis. Our analysis is based on photoionization models computed using the stellar ionizing radiation predicted by population synthesis codes (essentially Starburst 99) and, for the AGNs, a broken power-law spectrum. We explain why, among the four classical emission line

ABSTRACT We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power-law strain spectrum of the form $h_c = A(f/1\, \mathrm{yr}^{-1})^{\alpha }$, we found strong evidence for a spectrally similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times 10^{-15}$ and spectral index α = −0.5 ± 0.5, where the uncertainties represent 95 per cent credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of α = −2/3, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times 10^{-15}$. None the less, no significant evidence of the Hellings–Downs correlations that would indicate a gravitational-wave origin was found. We also analysed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.

We have re-evaluated empirical expressions for the abundance determination of N, O, Ne, S, Cl, Ar and Fe taking into account the latest atomic data and constructing an appropriate grid of photoionization models with state-of-the art model atmospheres. Using these expressions we have derived heavy element abundances in the $\\sim$ 310 emission-line galaxies from the Data Release 3 of the Sloan Digital Sky Survey (SDSS)with an observed Hbeta flux F(Hbeta)&gt; 1E-14 erg s^{-1} cm^{-2} and for which the [O III] 4363 emission line was detected at least at a 2sigma level, allowing abundance determination by direct methods. The oxygen abundance 12 + log O/H of the SDSS galaxies lies in the range from ~ 7.1 (Zsun/30) to 8.5 (0.7 Zsun). The SDSS sample is merged with a sample of 109 blue compact dwarf (BCD) galaxies with high quality spectra, which contains extremely low-metallicity objects. We use the merged sample to study the abundance patterns of low-metallicity emission-line galaxies. We find that extremely metal-poor galaxies (12 + log O/H &lt; 7.6, i.e. Z &lt; Zsun/12) are rare in the SDSS sample. The alpha element-to-oxygen abundance ratios do not show any significant trends with oxygen abundance, in agreement with previous studies, except for a slight increase of Ne/O with increasing metallicity, which we interpret as due to a moderate depletion of O onto grains in the most metal-rich galaxies. The Fe/O abundance ratio is smaller than the solar value, by up to 1 dex at the high metallicity end. (abridged)

Abstract The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in 2019 April and lasting six months, O3b starting in 2019 November and lasting five months, and O3GK starting in 2020 April and lasting two weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org . The main data set, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages.