Source author record

L. Silva

L. Silva appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

61works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

61 published item(s)

preprint2022arXiv

EOS-ESTM: A flexible climate model for habitable exoplanets

Rocky planets with temperate conditions provide the best chance for discovering habitable worlds and life outside the Solar System. In the last decades, new instrumental facilities and large observational campaigns have been driven by the quest for habitable worlds. Climate models aimed at studying the habitability of rocky planets are essential tools to pay off these technological and observational endeavours. In this context, we present EOS-ESTM, a fast and flexible model aimed at exploring the impact on habitability of multiple climate factors, including those unconstrained by observations. EOS-ESTM is built on ESTM, a seasonal-latitudinal energy balance model featuring an advanced treatment of the meridional and vertical transport. The novel features of EOS-ESTM include: (1) parameterizations for simulating the climate impact of oceans, land, ice, and clouds as a function of temperature and stellar zenith distance; (2) a procedure (EOS) for calculating the radiative transfer in atmospheres with terrestrial and non-terrestrial compositions illuminated by solar- and non-solar-type stars. By feeding EOS-ESTM with Earth's stellar, orbital and planetary parameters we derive a reference model that satisfies a large number of observational constraints of the Earth's climate system. Validation tests of non-terrestrial conditions yield predictions that are in line with comparable results obtained with a hierarchy of climate models. The application of EOS-ESTM to planetary atmospheres in maximum greenhouse conditions demonstrates the possibility of tracking the snowball transition at the outer edge of the HZ for a variety of planetary parameters, paving the road for multi-parametric studies of the HZ.

preprint2022arXiv

The impact of very massive stars on the chemical evolution of extremely metal-poor galaxies

In recent observations of extremely metal-poor low-mass starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. We investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as Pair-Instability Supernovae (PISN). We run chemical evolution models with yields that account for the contribution by PISN. We use both the recent non-rotating stellar yields from Goswami et al. 2021, and new yields from rotating very massive stars calculated on purpose in this work. We also search for the best initial mass function (IMF) that is able to reproduce the observations. We can reproduce the observations by adopting a bi-modal IMF and by including an initial burst of rotating very massive stars. Only with a burst of very massive stars can we reproduce the almost solar Fe/O ratios at the estimated young ages. We also confirm that rotation is absolutely needed to concomitantly reproduce the observed N/O ratios. These results stress the importance of very massive stars in galactic chemical evolution studies and strongly support a top-heavy initial mass function in the very early evolutionary stages of metal poor starburst galaxies.

preprint2020arXiv

Directed, elliptic and higher order flow harmonics of protons, deuterons and tritons in Au+Au collisions at $\sqrt{s_{NN}} = 2.4$ GeV

Flow coefficients $v_{n}$ of the orders $n = 1 - 6$ are measured with the High-Acceptance DiElectron Spectrometer (HADES) at GSI for protons, deuterons and tritons as a function of centrality, transverse momentum and rapidity in Au+Au collisions at $\sqrt{s_{NN}} = 2.4$ GeV. Combining the information from the flow coefficients of all orders allows to construct for the first time, at collision energies of a few GeV, a multi-differential picture of the angular emission pattern of these particles. It reflects the complicated interplay between the effect of the central fireball pressure on the emission of particles and their subsequent interaction with spectator matter. The high precision information on higher order flow coefficients is a major step forward in constraining the equation-of-state of dense baryonic matter.

preprint2020arXiv

Identical pion intensity interferometry at sqrt(s_{NN})=2.4 GeV

High-statistics $π^-π^-$ and $π^+π^+$ femtoscopy data are presented for Au+Au collisions at $\sqrt{s_\mathrm{NN}}=2.4$ GeV, measured with HADES at SIS18/GSI. The experimental correlation functions allow the determination of the space-time extent of the corresponding emission sources via a comparison to models. The emission source, parametrized as three-dimensional Gaussian distribution, is studied in dependence on pair transverse momentum, azimuthal emission angle with respect to the reaction plane, collision centrality and beam energy. For all centralities and transverse momenta, a geometrical distribution of ellipsoidal shape is found in the plane perpendicular to the beam direction with the larger extension perpendicular to the reaction plane. For large transverse momenta, the corresponding eccentricity approaches the initial eccentricity. The eccentricity is smallest for most central collisions, where the shape is almost circular. The magnitude of the tilt angle of the emission ellipsoid in the reaction plane decreases with increasing centrality and increasing transverse momentum. All source radii increase with centrality, largely exhibiting a linear rise with the number of participants, irrespective of transverse momentum. A substantial charge-sign difference of the source radii is found, appearing most pronounced at low transverse momentum. The extracted source parameters are consistent with the extrapolation of their energy dependence down from higher energies.

preprint2020arXiv

Proton number fluctuations in $\sqrt{s_{NN}}$ = 2.4 GeV Au+Au collisions studied with HADES

We present an analysis of proton number fluctuations in $\sqrt{s_{NN}}$ = 2.4 GeV Au+Au collisions measured with the High-Acceptance DiElectron Spectrometer (HADES) at GSI. With the help of extensive detector simulations done with IQMD transport model events including nuclear clusters, various nuisance effects influencing the observed proton cumulants have been investigated. Acceptance and efficiency corrections have been applied as a function of fine grained rapidity and transverse momentum bins, as well as considering local track density dependencies. Next, the effects of volume changes within particular centrality selections have been considered and beyond-leading-order corrections have been applied to the data. The efficiency and volume corrected proton number moments and cumulants Kn of orders n = 1, . . . , 4 have been obtained as a function of centrality and phase-space bin, as well as the corresponding correlators C_n . We find that the observed correlators show a power-law scaling with the mean number of protons, i.e. $C_n \propto <N>^n$, indicative of mostly long-range multi-particle correlations in momentum space. We also present a comparison of our results with Au+Au collision data obtained at RHIC at similar centralities, but higher $\sqrt{s_{NN}}$.

preprint2020arXiv

Two-Pion production in the second resonance region in $π^- p$ collisions with HADES

Pion induced reactions provide unique opportunities for an unambiguous description of baryonic resonances and their coupling channels by means of a partial wave analysis. Using the secondary pion beam at SIS18, the two pion production in the second resonance region has been investigated to unravel the role of the $N(1520) \frac{3}{2}^-$ resonance in the intermediate $ρ$ production. Results on exclusive channels with one pion ($π^{-}p$) and two pions ($π^{+}π^{-}n$, $π^{0}π^{-}p$) in the final state measured in the $π^{-}-p$ reaction at four different pion beam momenta (0.650, 0.685, 0.733, and 0.786 GeV/c) are presented. The excitation function of the different partial waves and $Δπ$, $Nσ$ and $Nρ$ isobar configurations is obtained, using the Bonn-Gatchina partial wave analysis. The $N(1520) \frac{3}{2}^-$ resonance is found to dominate the $Nρ$ final state with the branching ratio $BR=12.2 \pm 1.9\%$.

preprint2019arXiv

Measurement of the cross section for hard exclusive $π^0$ leptoproduction

We report on a measurement of hard exclusive $π^0$ muoproduction on the proton by COMPASS using 160 GeV/$c$ polarised $μ^+$ and $μ^-$ beams of the CERN SPS impinging on a liquid hydrogen target. From the average of the measured $μ^+$ and $μ^-$ cross sections, the virtual-photon proton cross section is determined as a function of the squared four-momentum transfer between initial and final proton in the range $0.08\,(\text{GeV/}c)^2 < |t| < 0.64\,(\text{GeV/}c)^2$. The average kinematics of the measurement are $\langle Q^2 \rangle =2.0\; {(\text{GeV}/c)^2}$, $\langle ν\rangle = 12.8\; {\text{GeV}}$, $\langle x_{Bj} \rangle = 0.093 $ and $\langle -t \rangle = 0.256\; {(\text{GeV}/c)^2} $. Fitting the azimuthal dependence reveals a combined contribution by transversely and longitudinally polarised photons of $(8.1 \ \pm \ 0.9_{\text{stat}}{}_{- \ 1.0}^{+ \ 1.1}\big\rvert_{\text{sys}})\,{\text{nb}}/{(\text{GeV}/c)^{2}}$, as well as transverse-transverse and longitudinal-transverse interference contributions of $(-6.0 \pm 1.3_{\text{stat}}{}_{- \ 0.7}^{+ \ 0.7}\big\rvert_{\text{sys}})\,{\text{nb}}/{(\text{GeV}/c)^{2}}$ and $(1.4 \pm 0.5_{\text{stat}}{}_{- \ 0.2}^{+ \ 0.3}\big\rvert_{\text{sys}})\,{\text{nb}}/{(\text{GeV}/c)^{2}}$, respectively. Our results provide important input for modelling Generalised Parton Distributions. In the context of the phenomenological Goloskokov-Kroll model, the statistically significant transverse-transverse interference contribution constitutes clear experimental evidence for the chiral-odd GPD $\overline{E}_T$.

preprint2017arXiv

AWAKE readiness for the study of the seeded self-modulation of a 400\,GeV proton bunch

AWAKE is a proton-driven plasma wakefield acceleration experiment. % We show that the experimental setup briefly described here is ready for systematic study of the seeded self-modulation of the 400\,GeV proton bunch in the 10\,m-long rubidium plasma with density adjustable from 1 to 10$\times10^{14}$\,cm$^{-3}$. % We show that the short laser pulse used for ionization of the rubidium vapor propagates all the way along the column, suggesting full ionization of the vapor. % We show that ionization occurs along the proton bunch, at the laser time and that the plasma that follows affects the proton bunch. %

preprint2016arXiv

Observation of a new narrow axial-vector meson $a_1(1420)$

The COMPASS collaboration at CERN has measured diffractive dissociation of 190 GeV/$c$ pions into the $π^-π^-π^+$ final state using a stationary hydrogen target. A partial-wave analysis (PWA) was performed in bins of $3π$ mass and four-momentum transfer using the isobar model and the so far largest PWA model consisting of 88 waves. A narrow $J^{PC} = 1^{++}$ signal is observed in the $f_0(980)\,π$ channel. We present a resonance-model study of a subset of the spin-density matrix selecting $3π$ states with $J^{PC} = 2^{++}$ and $4^{++}$ decaying into $ρ(770)\,π$ and with $J^{PC} = 1^{++}$ decaying into $f_0(980)\,π$. We identify a new $a_1$ meson with mass $(1414^{+15}_{-13})$ MeV/$c^2$ and width $(153^{+8}_{-23})$ MeV/$c^2$. Within the final states investigated in our analysis, we observe the new $a_1(1420)$ decaying only into $f_0(980)\,π$, suggesting its exotic nature.

preprint2016arXiv

The $\bf{Λp}$ interaction studied via femtoscopy in p + Nb reactions at $\mathbf{\sqrt{s_{NN}}=3.18} ~\mathrm{\bf{GeV}}$

We report on the first measurement of $pΛ$ and $pp$ correlations via the femtoscopy method in p+Nb reactions at $\mathrm{\sqrt{s_{NN}}=3.18} ~\mathrm{GeV}$, studied with the High Acceptance Di-Electron Spectrometer (HADES). By comparing the experimental correlation function to model calculations, a source size for $pp$ pairs of $r_{0,pp}=2.02 \pm 0.01(\mathrm{stat})^{+0.11}_{-0.12} (\mathrm{sys}) ~\mathrm{fm}$ and a slightly smaller value for $pΛ$ of $r_{0,Λp}=1.62 \pm 0.02(\mathrm{stat})^{+0.19}_{-0.08}(\mathrm{sys}) ~\mathrm{fm}$ is extracted. Using the geometrical extent of the particle emitting region, determined experimentally with $pp$ correlations as reference together with a source function from a transport model, it is possible to study different sets of scattering parameters. The $pΛ$ correlation is proven sensitive to predicted scattering length values from chiral effective field theory. We demonstrate that the femtoscopy technique can be used as valid alternative to the analysis of scattering data to study the hyperon-nucleon interaction.

preprint2015arXiv

Are z>2 Herschel galaxies proto-spheroids?

We present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.

preprint2015arXiv

AWAKE, The Advanced Proton Driven Plasma Wakefield Acceleration Experiment at CERN

The Advanced Proton Driven Plasma Wakefield Acceleration Experiment (AWAKE) aims at studying plasma wakefield generation and electron acceleration driven by proton bunches. It is a proof-of-principle R&D experiment at CERN and the world's first proton driven plasma wakefield acceleration experiment. The AWAKE experiment will be installed in the former CNGS facility and uses the 400 GeV/c proton beam bunches from the SPS. The first experiments will focus on the self-modulation instability of the long (rms ~12 cm) proton bunch in the plasma. These experiments are planned for the end of 2016. Later, in 2017/2018, low energy (~15 MeV) electrons will be externally injected to sample the wakefields and be accelerated beyond 1 GeV. The main goals of the experiment will be summarized. A summary of the AWAKE design and construction status will be presented.

preprint2015arXiv

Interplay among transversity induced asymmetries in hadron leptoproduction

In the fragmentation of a transversely polarized quark several left-right asymmetries are possible for the hadrons in the jet. When only one unpolarized hadron is selected, it exhibits an azimuthal modulation known as Collins effect. When a pair of oppositely charged hadrons is observed, three asymmetries can be considered, a di-hadron asymmetry and two single hadron asymmetries. In lepton deep inelastic scattering on transversely polarized nucleons all these asymmetries are coupled with the transversity distribution. From the high statistics COMPASS data on oppositely charged hadron-pair production we have investigated for the first time the dependence of these three asymmetries on the difference of the azimuthal angles of the two hadrons. The similarity of transversity induced single and di-hadron asymmetries is discussed. A new analysis of the data allows to establish quantitative relationships among them, providing for the first time strong experimental indication that the underlying fragmentation mechanisms are all driven by a common physical process.

preprint2015arXiv

Longitudinal double spin asymmetries in single hadron quasi-real photoproduction at high $p_T$

We measured the longitudinal double spin asymmetries $A_{LL}$ for single hadron muo-production off protons and deuterons at photon virtuality $Q^2$ < 1(GeV/$\it c$)$^2$ for transverse hadron momenta $p_T$ in the range 0.7 GeV/$\it c$ to 4 GeV/$\it c$ . They were determined using COMPASS data taken with a polarised muon beam of 160 GeV/$\it c$ or 200 GeV/$\it c$ impinging on polarised $\mathrm{{}^6LiD}$ or $\mathrm{NH_3}$ targets. The experimental asymmetries are compared to next-to-leading order pQCD calculations, and are sensitive to the gluon polarisation $ΔG$ inside the nucleon in the range of the nucleon momentum fraction carried by gluons $0.05 < x_g < 0.2$.

preprint2015arXiv

Measurement of the charged-pion polarisability

The COMPASS collaboration at CERN has investigated pion Compton scattering, $π^-γ\rightarrow π^-γ$, at centre-of-mass energy below 3.5 pion masses. The process is embedded in the reaction $π^-\mathrm{Ni}\rightarrowπ^-γ\;\mathrm{Ni}$, which is initiated by 190\,GeV pions impinging on a nickel target. The exchange of quasi-real photons is selected by isolating the sharp Coulomb peak observed at smallest momentum transfers, $Q^2<0.0015$\,(GeV/$c$)$^2$. From a sample of 63\,000 events the pion electric polarisability is determined to be $α_π =\ (\,2.0\ \pm\ 0.6_{\mbox{\scriptsize stat}}\ \pm\ 0.7_{\mbox{\scriptsize syst}}\,) \times 10^{-4}\,\mbox{fm}^3$ under the assumption $α_π=-β_π$, which relates the electric and magnetic dipole polarisabilities. It is the most precise measurement of this fundamental low-energy parameter of strong interaction, that has been addressed since long by various methods with conflicting outcomes. While this result is in tension with previous dedicated measurements, it is found in agreement with the expectation from chiral perturbation theory. An additional measurement replacing pions by muons, for which the cross-section behavior is unambigiously known, was performed for an independent estimate of the systematic uncertainty.

preprint2015arXiv

Modeling the surface temperature of Earth-like planets

We introduce a novel Earth-like planet surface temperature model (ESTM) for habitability studies based on the spatial-temporal distribution of planetary surface temperatures. The ESTM adopts a surface Energy Balance Model complemented by: radiative-convective atmospheric column calculations, a set of physically-based parameterizations of meridional transport, and descriptions of surface and cloud properties more refined than in standard EBMs. The parameterization is valid for rotating terrestrial planets with shallow atmospheres and moderate values of axis obliquity (epsilon >= 45^o). Comparison with a 3D model of atmospheric dynamics from the literature shows that the equator-to-pole temperature differences predicted by the two models agree within ~5K when the rotation rate, insolation, surface pressure and planet radius are varied in the intervals 0.5 <= Omega/Omega_o <= 2, 0.75 <= S/S_o <= 1.25, 0.3 <= p/(1 bar) <= 10, and 0.5 <= R/R_o <= 2, respectively. The ESTM has an extremely low computational cost and can be used when the planetary parameters are scarcely known (as for most exoplanets) and/or whenever many runs for different parameter configurations are needed. Model simulations of a test-case exoplanet (Kepler-62e) indicate that an uncertainty in surface pressure within the range expected for terrestrial planets may impact the mean temperature by ~60 K. Within the limits of validity of the ESTM, the impact of surface pressure is larger than that predicted by uncertainties in rotation rate, axis obliquity, and ocean fractions. We discuss the possibility of performing a statistical ranking of planetary habitability taking advantage of the flexibility of the ESTM.

preprint2015arXiv

Path to AWAKE: Evolution of the concept

This report describes the conceptual steps in reaching the design of the AWAKE experiment currently under construction at CERN. We start with an introduction to plasma wakefield acceleration and the motivation for using proton drivers. We then describe the self-modulation instability --- a key to an early realization of the concept. This is then followed by the historical development of the experimental design, where the critical issues that arose and their solutions are described. We conclude with the design of the experiment as it is being realized at CERN and some words on the future outlook. A summary of the AWAKE design and construction status as presented in this conference is given in [1].

preprint2015arXiv

The Spin Structure Function $g_1^{\rm p}$ of the Proton and a Test of the Bjorken Sum Rule

New results for the double spin asymmetry $A_1^{\rm p}$ and the proton longitudinal spin structure function $g_1^{\rm p}$ are presented. They were obtained by the COMPASS collaboration using polarised 200 GeV muons scattered off a longitudinally polarised NH$_3$ target. The data were collected in 2011 and complement those recorded in 2007 at 160\,GeV, in particular at lower values of $x$. They improve the statistical precision of $g_1^{\rm p}(x)$ by about a factor of two in the region $x\lesssim 0.02$. A next-to-leading order QCD fit to the $g_1$ world data is performed. It leads to a new determination of the quark spin contribution to the nucleon spin, $ΔΣ$ ranging from 0.26 to 0.36, and to a re-evaluation of the first moment of $g_1^{\rm p}$. The uncertainty of $ΔΣ$ is mostly due to the large uncertainty in the present determinations of the gluon helicity distribution. A new evaluation of the Bjorken sum rule based on the COMPASS results for the non-singlet structure function $g_1^{\rm NS}(x,Q^2)$ yields as ratio of the axial and vector coupling constants $|g_{\rm A}/g_{\rm V}| = 1.22 \pm 0.05~({\rm stat.}) \pm 0.10~({\rm syst.})$, which validates the sum rule to an accuracy of about 9\%.

preprint2014arXiv

A high-statistics measurement of transverse spin effects in dihadron production from muon-proton semi-inclusive deep-inelastic scattering

A measurement of the azimuthal asymmetry in dihadron production in deep-inelastic scattering of muons on transversely polarised proton (NH$_{3}$) targets are presented. They provide independent access to the transversity distribution functions through the measurement of the Collins asymmetry in single hadron production. The data were taken in the year $2010$ with the COMPASS spectrometer using a $160\,\mbox{GeV}/c$ muon beam of the CERN SPS, increasing by a factor of about three the available statistics of the previously published data taken in the year $2007$. The measured sizeable asymmetry is in good agreement with the published data. An approximate equality of the Collins asymmetry and the dihadron asymmetry is observed, suggesting a common physical mechanism in the underlying fragmentation.

preprint2014arXiv

Collins and Sivers asymmetries in muonproduction of pions and kaons off transversely polarised proton

Measurements of the Collins and Sivers asymmetries for charged pions and charged and neutral kaons produced in semi-inclusive deep-inelastic scattering of high energy muons off transversely polarised protons are presented. The results were obtained using all the available COMPASS proton data, which were taken in the years 2007 and 2010. The Collins asymmetries exhibit in the valence region a non-zero signal for pions and there are hints of non-zero signal also for kaons. The Sivers asymmetries are found to be positive for positive pions and kaons and compatible with zero otherwise.

preprint2014arXiv

Combining Physical galaxy models with radio observations to constrain the SFRs of high-z dusty star forming galaxies

We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs.

preprint2014arXiv

GRASIL-3D: an Implemention of Dust Effects in the SEDs of Simulated Galaxies

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

preprint2014arXiv

Measurement of azimuthal hadron asymmetries in semi-inclusive deep inelastic scattering off unpolarised nucleons

Spin-averaged asymmetries in the azimuthal distributions of positive and negative hadrons produced in deep inelastic scattering were measured using the CERN SPS muon beam at $160$ GeV/c and a $^6$LiD target. The amplitudes of the three azimuthal modulations $\cosϕ_h$, $\cos2ϕ_h$ and $\sinϕ_h$ were obtained binning the data separately in each of the relevant kinematic variables $x$, $z$ or $p_T^{\,h}$ and binning in a three-dimensional grid of these three variables. The amplitudes of the $\cos ϕ_h$ and $\cos 2ϕ_h$ modulations show strong kinematic dependencies both for positive and negative hadrons.

preprint2014arXiv

Measurement of radiative widths of $a_2(1320)$ and $π_2(1670)$

The COMPASS Collaboration at CERN has investigated the reaction $π^- γ\rightarrow π^-π^-π^+$ embedded in the Primakoff reaction of $190~\textrm{GeV}$ pions scattering in the Coulomb field of a lead target, $π^- \text{Pb} \rightarrow π^-π^-π^+ \text{Pb}$. Exchange of quasi-real photons is selected by isolating the sharp Coulomb peak observed at momentum transfer below $0.001~(\text{GeV}/c)^2$. Using a partial-wave analysis the amplitudes and relative phases of the $a_2(1320)$ and $π_2(1670)$ mesons have been extracted, and the Coulomb and the diffractive contributions have been disentangled. Measuring absolute production cross sections we have determined the radiative width of the $a_2(1320)$ to be $Γ_0(a_2(1320) \rightarrow πγ) = (358 \pm 6_{\textrm{stat}} \pm 42_{\textrm{syst}})~\textrm{keV}$. As the first measurement, $Γ_0(π_2(1670) \rightarrow πγ) = (181 \pm 11_{\textrm{stat}} \pm 27_{\textrm{syst}})~\textrm{keV} \cdot (\textrm{BR}^{\textrm{PDG}}_{f_2 π}/\textrm{BR}_{f_2 π})$ is obtained for the radiative width of the $π_2(1670)$, where in this analysis the branching ratio $\textrm{BR}^{\textrm{PDG}}_{f_2 π}=0.56$ has been used. We compare these values to previous measurements and theoretical predictions.

preprint2014arXiv

Odd and Even Partial Waves of $ηπ^-$ and $η'π^-$ in $π^-p\toη^{(\prime)}π^-p$ at $191\,\textrm{GeV}/c$

Exclusive production of $ηπ^-$ and $η'π^-$ has been studied with a $191\,\textrm{GeV}/c$ $π^-$ beam impinging on a hydrogen target at COMPASS (CERN). Partial-wave analyses reveal different odd/even angular momentum ($L$) characteristics in the inspected invariant mass range up to $3\,\textrm{GeV}/c^2$. A striking similarity between the two systems is observed for the $L=2,4,6$ intensities (scaled by kinematical factors) and the relative phases. The known resonances $a_2(1320)$ and $a_4(2040)$ are in line with this similarity. In contrast, a strong enhancement of $η'π^-$ over $ηπ^-$ is found for the $L=1,3,5$ waves, which carry non-$q\bar q$ quantum numbers. The $L=1$ intensity peaks at $1.7\,\textrm{GeV}/c^2$ in $η'π^-$ and at $1.4\,\textrm{GeV}/c^2$ in $ηπ^-$, the corresponding phase motions with respect to $L=2$ are different.

preprint2014arXiv

Search for exclusive photoproduction of Z$_c^{\pm}$(3900) at COMPASS

A search for the exclusive production of the $Z_c^{\pm}(3900)$ hadron by virtual photons has been performed in the channel $Z_c^{\pm}(3900)\rightarrow J/ψπ^{\pm}$. The data cover the range from 7 GeV to 19 GeV in the centre-of-mass energy of the photon-nucleon system. The full set of the COMPASS data set collected with a muon beam between 2002 and 2011 has been used. An upper limit for the ratio $BR(Z_c^{\pm}(3900)\rightarrow J/ψπ^{\pm} )\times σ_{ γ~N \rightarrow Z_c^{\pm}(3900)~ N} /σ_{γ~N \rightarrow J/ψ~ N}$ of $3.7\times10^{-3}$ has been established at the confidence level of 90%.

preprint2014arXiv

Spin alignment and violation of the OZI rule in exclusive $ω$ and $ϕ$ production in pp collisions

Exclusive production of the isoscalar vector mesons $ω$ and $ϕ$ is measured with a 190 GeV$/c$ proton beam impinging on a liquid hydrogen target. Cross section ratios are determined in three intervals of the Feynman variable $x_{F}$ of the fast proton. A significant violation of the OZI rule is found, confirming earlier findings. Its kinematic dependence on $x_{F}$ and on the invariant mass $M_{p\mathrm{V}}$ of the system formed by fast proton $p_\mathrm{fast}$ and vector meson $V$ is discussed in terms of diffractive production of $p_\mathrm{fast}V$ resonances in competition with central production. The measurement of the spin density matrix element $ρ_{00}$ of the vector mesons in different selected reference frames provides another handle to distinguish the contributions of these two major reaction types. Again, dependences of the alignment on $x_{F}$ and on $M_{p\mathrm{V}}$ are found. Most of the observations can be traced back to the existence of several excited baryon states contributing to $ω$ production which are absent in the case of the $ϕ$ meson. Removing the low-mass $M_{p\mathrm{V}}$ resonant region, the OZI rule is found to be violated by a factor of eight, independently of $x_\mathrm{F}$.

preprint2014arXiv

The COMPASS Setup for Physics with Hadron Beams

The main characteristics of the COMPASS experimental setup for physics with hadron beams are described. This setup was designed to perform exclusive measurements of processes with several charged and/or neutral particles in the final state. Making use of a large part of the apparatus that was previously built for spin structure studies with a muon beam, it also features a new target system as well as new or upgraded detectors. The hadron setup is able to operate at the high incident hadron flux available at CERN. It is characterised by large angular and momentum coverages, large and nearly flat acceptances, and good two and three-particle mass resolutions. In 2008 and 2009 it was successfully used with positive and negative hadron beams and with liquid hydrogen and solid nuclear targets. This article describes the new and upgraded detectors and auxiliary equipment, outlines the reconstruction procedures used, and summarises the general performance of the setup.

preprint2014arXiv

The Main Sequence and the Fundamental Metallicity Relation in MaGICC Galaxies: Evolution and Scatter

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

preprint2013arXiv

Hadron Transverse Momentum Distributions in Muon Deep Inelastic Scattering at 160 GeV/$c$

Multiplicities of charged hadrons produced in deep inelastic muon scattering off a $^6$LiD target have been measured as a function of the DIS variables $x_{Bj}$, $Q^2$, $W^2$ and the final state hadron variables $p_T$ and $z$. The $p_T^2$ distributions are fitted with a single exponential function at low values of $p_T^2$ to determine the dependence of $\langle p_T^2 \rangle$ on $x_{Bj}$, $Q^2$, $W^2$ and $z$. The $z$-dependence of $\langle p_T^2 \rangle$ is shown to be a potential tool to extract the average intrinsic transverse momentum squared of partons, $\langle k_{\perp}^2 \rangle$, as a function of $x_{Bj}$ and $Q^2$ in a leading order QCD parton model.

preprint2013arXiv

Measurement of the Cross Section for High-p_T Hadron Production in Scattering of 160 GeV/c Muons off Nucleons

The differential cross section for production of charged hadrons with high transverse momenta in scattering of 160\,GeV/$c$ muons off nucleons at low photon virtualities has been measured at the COMPASS experiment at CERN. The results, which cover transverse momenta from 1.1\,GeV/$c$ to 3.6\,GeV/$c$, are compared to a perturbative Quantum Chromodynamics (pQCD) calculation, in order to evaluate the applicability of pQCD to this process in the kinematic domain of the experiment. The shape of the calculated differential cross section as a function of transverse momentum is found to be in good agreement with the experimental data, but the absolute scale is underestimated by next-to-leading order (NLO) pQCD. The inclusion of all-order resummation of large logarithmic threshold corrections reduces the discrepancy from a factor of three to four to a factor of two. The dependence of the cross section on the pseudo-rapidity and on virtual photon energy fraction is investigated. Finally the dependence on the charge of the hadrons is discussed.

preprint2013arXiv

Study of Σ(1385) and Ξ(1321) hyperon and antihyperon production in deep inelastic muon scattering

Large samples of Λ, Σ(1385) and Ξ(1321) hyperons produced in deep-inelastic muon scattering off a ^6LiD target were collected with the COMPASS experimental setup at CERN. The relative yields of Σ(1385)^+, Σ(1385)^-, \barΣ(1385)^-, \barΣ(1385)^+, Ξ(1321)^-, and \barΞ(1321)^+ hyperons decaying into Λ(\barΛ)πwere measured. The heavy hyperon to Λand heavy antihyperon to \barΛ yield ratios were found to be in the range 3.8% to 5.6% with a relative uncertainty of about 10%. They were used to tune the parameters relevant for strange particle production of the LEPTO Monte Carlo generator.

preprint2013arXiv

The evolution of the dust and gas content in galaxies

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

preprint2013arXiv

Transverse target spin asymmetries in exclusive $ρ^0$ muoproduction

Exclusive production of $ρ^0$ mesons was studied at the COMPASS experiment by scattering 160 GeV/$c$ muons off transversely polarised protons. Five single-spin and three double-spin azimuthal asymmetries were measured as a function of $Q^2$, $x_{Bj}$, or $p_{T}^{2}$. The $\sin ϕ_S$ asymmetry is found to be $-0.019 \pm 0.008(stat.) \pm 0.003(syst.)$. All other asymmetries are also found to be of small magnitude and consistent with zero within experimental uncertainties. Very recent calculations using a GPD-based model agree well with the present results. The data is interpreted as evidence for the existence of chiral-odd, transverse generalized parton distributions.

preprint2012arXiv

Comparison of star formation rates from Halpha and infrared luminosities as seen by Herschel

We empirically test the relation between the SFR(LIR) derived from the infrared luminosity, LIR, and the SFR(Ha) derived from the Ha emission line luminosity using simple conversion relations. We use a sample of 474 galaxies at z = 0.06 - 0.46 with both Ha detection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 μm). We derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very clear trend between E(B - V) and LIR that allows to estimate extinction values for each galaxy even if the Ha emission line measurement is not reliable. We calculate the LIR by integrating from 8 up to 1000 μm the SED that is best fitting our data. We compare SFR(Ha) with the SFR(LIR). We find a very good agreement between the two SFR estimates, with a slope of m = 1.01 \pm 0.03 in the SFR(LIR) vs SFR(Ha) diagram, a normalization constant of a = -0.08 \pm 0.03 and a dispersion of sigma = 0.28 dex.We study the effect of some intrinsic properties of the galaxies in the SFR(LIR)-SFR(Ha) relation, such as the redshift, the mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR comparison. The mean ratio of the two SFR estimators log[SFR(LIR)/SFR(Ha)] varies by approx. 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(O/H) + 12 < 9.2). This effect is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. Considering different morphological types, we find a very good agreement between the two SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much steeper (m = 1.62 \pm 0.43), but it is still consistent with 1 at the 1.5 sigma level, because of the reduced statistics of this sub-sample.

preprint2012arXiv

Cosmic star formation rate: a theoretical approach

The cosmic star formation rate (CSFR), is an important clue to investigate the history of the assembly and evolution of galaxies. Here, we develop a method to study the CSFR from a purely theoretical point of view. Starting from detailed models of chemical evolution, we obtain the histories of star formation of galaxies of different morphological types. These histories are then used to determine the luminosity functions of the same galaxies by means of a spectro-photometric code. We obtain the CSFR under different hypothesis. First, we study the hypothesis of a pure luminosity evolution scenario, in which all galaxies are supposed to form at the same redshift and then evolve only in luminosity. Then we consider scenarios in which the number density or the slope of the LFs are assumed to vary with redshift. After comparison with available data we conclude that a pure luminosity evolution does not provide a good fit to the data, especially at very high redshift, although many uncertainties are still present in the data. On the other hand, a variation in the number density of ellipticals and spirals as a function of redshift can provide a better fit to the observed CSFR. We also explore cases of variable slope of the LFs with redshift and variations of number density and slope at the same time. We cannot find any of those cases which can improve the fit to the data respect to the solely number density variation. Finally, we compute the evolution of the average cosmic metallicity in galaxies with redshift.

preprint2012arXiv

D* and D Meson Production in Muon Nucleon Interactions at 160 GeV/c

The production of D* and D mesons in inelastic scattering of 160 GeV/c muons off a ^6LiD target has been investigated with the COMPASS spectrometer at CERN for 0.003 (GeV/c)^2 < Q^2 < 10 (GeV/c)^2 and 3x10^-5< x_Bj < 0.1. The study is based on 8100 events where a D^0 or anti D^0 is detected subsequently to a D*+ or D*- decay, and on 34000 events, where only a D^0 or anti D^0 is detected. Kinematic distributions of D*, D and K*_2 are given as a function of their energy E, transverse momentum p_T, energy fraction z, and of the virtual photon variables nu, Q^2 and x_Bj. Semi-inclusive differential D* production cross-sections are compared with theoretical predictions for D* production via photon-gluon fusion into open charm. The total observed production cross-section for D*+/- mesons with laboratory energies between 22 and 86 GeV is 1.9 nb. Significant cross-section asymmetries are observed between D*+ and D*- production for nu<40 GeV and z>0.6.

preprint2012arXiv

Exclusive rho^0 muoproduction on transversely polarised protons and deuterons

The transverse target spin azimuthal asymmetry A_UT in hard exclusive production of rho^0 mesons was measured at COMPASS by scattering 160 GeV/c muons off transversely polarised protons and deuterons. The measured asymmetry is sensitive to the nucleon helicity-flip generalised parton distributions E^q, which are related to the orbital angular momentum of quarks in the nucleon. The Q^2, x_B and p_t^2 dependence of A_UT is presented in a wide kinematic range. Results for deuterons are obtained for the first time. The measured asymmetry is small in the whole kinematic range for both protons and deuterons, which is consistent with the theoretical interpretation that contributions from GPDs E^u and E^d approximately cancel.

preprint2012arXiv

Experimental investigation of transverse spin asymmetries in muon-p SIDIS processes: Collins asymmetries

The COMPASS Collaboration at CERN has measured the transverse spin azimuthal asymmetry of charged hadrons produced in semi-inclusive deep inelastic scattering using a 160 GeV positive muon beam and a transversely polarised NH_3 target. The Collins asymmetry of the proton was extracted in the Bjorken x range 0.003<x<0.7. These new measurements confirm with higher accuracy previous measurements from the COMPASS and HERMES collaborations, which exhibit a definite effect in the valence quark region. The asymmetries for negative and positive hadrons are similar in magnitude and opposite in sign. They are compatible with model calculations in which the u-quark transversity is opposite in sign and somewhat larger than the d-quark transversity distribution function. The asymmetry is extracted as a function of Bjorken $x$, the relative hadron energy $z$ and the hadron transverse momentum p_T^h. The high statistics and quality of the data also allow for more detailed investigations of the dependence on the kinematic variables. These studies confirm the leading-twist nature of the Collins asymmetry.

preprint2012arXiv

Experimental investigation of transverse spin asymmetries in muon-p SIDIS processes: Sivers asymmetries

The COMPASS Collaboration at CERN has measured the transverse spin azimuthal asymmetry of charged hadrons produced in semi-inclusive deep inelastic scattering using a 160 GeV positive muon beam and a transversely polarised NH_3 target. The Sivers asymmetry of the proton has been extracted in the Bjorken x range 0.003<x<0.7. The new measurements have small statistical and systematic uncertainties of a few percent and confirm with considerably better accuracy the previous COMPASS measurement. The Sivers asymmetry is found to be compatible with zero for negative hadrons and positive for positive hadrons, a clear indication of a spin-orbit coupling of quarks in a transversely polarised proton. As compared to measurements at lower energy, a smaller Sivers asymmetry for positive hadrons is found in the region x > 0.03. The asymmetry is different from zero and positive also in the low x region, where sea-quarks dominate. The kinematic dependence of the asymmetry has also been investigated and results are given for various intervals of hadron and virtual photon fractional energy. In contrast to the case of the Collins asymmetry, the results on the Sivers asymmetry suggest a strong dependence on the four-momentum transfer to the nucleon, in agreement with the most recent calculations.

preprint2012arXiv

First Measurement of Chiral Dynamics in π^- γ-> π^- π^- π^+

The COMPASS collaboration at CERN has investigated the π^- γ-> π^- π^- π^+ reaction at center-of-momentum energy below five pion masses, sqrt(s) < 5 m(π), embedded in the Primakoff reaction of 190 GeV pions impinging on a lead target. Exchange of quasi-real photons is selected by isolating the sharp Coulomb peak observed at smallest momentum transfers, t' < 0.001 (GeV/c)^2. Using partial-wave analysis techniques, the scattering intensity of Coulomb production described in terms of chiral dynamics and its dependence on the 3π-invariant mass m(3π) = sqrt(s) were extracted. The absolute cross section was determined in seven bins of $\sqrt{s}$ with an overall precision of 20%. At leading order, the result is found to be in good agreement with the prediction of chiral perturbation theory over the whole energy range investigated.

preprint2012arXiv

Leading and Next-to-Leading Order Gluon Polarisation in the Nucleon and Longitudinal Double Spin Asymmetries from Open Charm Muoproduction

The gluon polarisation in the nucleon was measured using open charm production by scattering 160 GeV/c polarised muons off longitudinally polarised protons or deuterons. The data were taken by the COMPASS collaboration between 2002 and 2007. A detailed account is given of the analysis method that includes the application of neural networks. Several decay channels of D^0 mesons are investigated. Longitudinal spin asymmetries of the D meson production cross-sections are extracted in bins of D^0 transverse momentum and energy. At leading order QCD accuracy the average gluon polarisation is determined as (Delta g/G)^LO=-0.06 +/- 0.21 (stat.) +/- 0.08 (syst.) at the scale <mu^2> ~13 (GeV/c)^2 and an average gluon momentum fraction <x>~ 0.11. For the first time, the average gluon polarisation is also obtained at next-to-leading order QCD accuracy as (Delta g/G)^NLO = -0.13 +/- 0.15 (stat.) +/- 0.15 (syst.) at the scale <mu^2> ~ 13 (\GeV/c)^2 and <x> ~ 0.20.

preprint2012arXiv

Leading order determination of the gluon polarisation from DIS events with high-p_T hadron pairs

We present a determination of the gluon polarisation Delta g/g in the nucleon, based on the longitudinal double-spin asymmetry of DIS events with a pair of large transverse-momentum hadrons in the final state. The data were obtained by the COMPASS experiment at CERN using a 160 GeV/c polarised muon beam scattering off a polarised ^6LiD target. The gluon polarisation is evaluated by a Neural Network approach for three intervals of the gluon momentum fraction x_g covering the range 0.04 < x_g < 0.27. The values obtained at leading order in QCD do not show any significant dependence on x_g. Their average is Delta g/g = 0.125 +/- 0.060 (stat.) +/- 0.063 (syst.) at x_g=0.09 and a scale of mu^2 = 3 (GeV/c)^2.

preprint2012arXiv

The complex physics of dusty star-forming galaxies at high redshifts as revealed by Herschel and Spitzer

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $Δ$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $Δ$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

preprint2012arXiv

Transverse spin effects in hadron-pair production from semi-inclusive deep inelastic scattering

First measurements of azimuthal asymmetries in hadron-pair production in deep-inelastic scattering of muons on transversely polarised ^6LiD (deuteron) and NH_3 (proton) targets are presented. The data were taken in the years 2002-2004 and 2007 with the COMPASS spectrometer using a muon beam of 160 GeV/c at the CERN SPS. The asymmetries provide access to the transversity distribution functions, without involving the Collins effect as in single hadron production. The sizeable asymmetries measured on the NH_ target indicate non-vanishing u-quark transversity and two-hadron interference fragmentation functions. The small asymmetries measured on the ^6LiD target can be interpreted as indication for a cancellation of u- and d-quark transversities.

preprint2010arXiv

Modelling the spectral energy distribution of galaxies: introducing the artificial neural network

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

preprint2010arXiv

Predictions for Herschel from LambdaCDM: unveiling the cosmic star formation history

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

preprint2010arXiv

The chemical evolution of elliptical galaxies with stellar and QSO dust production

The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is therefore important to study he chemical evolution of galaxies by taking into account dust evolution. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, respectively. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and intermediate mass stars, supernovae Ia, supernovae II, QSOs and both dust destruction and accretion processes. By means of such a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical downsizing. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the QSO SDSS J114816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and during the early development of the galactic wind.

preprint2010arXiv

The dust content of high-z submillimeter galaxies revealed by Herschel

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

preprint2010arXiv

Unusual PAH Emission in Nearby Early-Type Galaxies: A Signature of an Intermediate Age Stellar Population?

We present the analysis of Spitzer-IRS spectra of four early-type galaxies, NGC 1297, NGC 5044, NGC 6868, and NGC 7079, all classified as LINERs in the optical bands. Their IRS spectra present the full series of H2 rotational emission lines in the range 5--38 microns, atomic lines, and prominent PAH features. We investigate the nature and origin of the PAH emission, characterized by unusually low 6 -- 9/11.3 microns inter-band ratios. After the subtraction of a passive early type galaxy template, we find that the 7 -- 9 microns spectral region requires dust features not normally present in star forming galaxies. Each spectrum is then analyzed with the aim of identifying their components and origin. In contrast to normal star forming galaxies, where cationic PAH emission prevails, our 6--14 microns spectra seem to be dominated by large and neutral PAH emission, responsible for the low 6 -- 9/11.3 microns ratios, plus two broad dust emission features peaking at 8.2 microns and 12 microns. Theses broad components, observed until now mainly in evolved carbon stars and usually attributed to pristine material, contribute approximately 30-50% of the total PAH flux in the 6--14 microns region. We propose that the PAH molecules in our ETGs arise from fresh carbonaceous material which is continuously released by a population of carbon stars, formed in a rejuvenation episode which occurred within the last few Gyr. The analysis of the MIR spectra allows us to infer that, in order to maintain the peculiar size and charge distributions biased to large and neutral PAHs, this material must be shocked, and excited by the weak UV interstellar radiation field of our ETG.

preprint2009arXiv

Collins and Sivers asymmetries for pions and kaons in muon-deuteron DIS

The measurements of the Collins and Sivers asymmetries of identified hadrons produced in deep-inelastic scattering of 160 GeV/c muons on a transversely polarised 6LiD target at COMPASS are presented. The results for charged pions and charged and neutral kaons correspond to all data available, which were collected from 2002 to 2004. For all final state particles both the Collins and Sivers asymmetries turn out to be small, compatible with zero within the statistical errors, in line with the previously published results for not identified charged hadrons, and with the expected cancellation between the u- and d-quark contributions.

preprint2009arXiv

Modelling the dusty universe I: Introducing the artificial neural network and first applications to luminosity and colour distributions

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

preprint2009arXiv

The role of the synchrotron component in the mid infrared spectrum of M 87

We study in detail the mid-infrared Spitzer-IRS spectrum of M 87 in the range 5 to 20 micron. Thanks to the high sensitivity of our Spitzer-IRS spectra we can disentangle the stellar and nuclear components of this active galaxy. To this end we have properly subtracted from the M 87 spectrum, the contribution of the underlying stellar continuum, derived from passive Virgo galaxies in our sample. The residual is a clear power-law, without any additional thermal component, with a zero point consistent with that obtained by high spatial resolution, ground based observations. The residual is independent of the adopted passive template. This indicates that the 10 micron silicate emission shown in spectra of M 87 can be entirely accounted for by the underlying old stellar population, leaving little room for a possible torus contribution. The MIR power-law has a slope alpha ~ 0.77-0.82 (S$_ν\proptoν^{-α}$), consistent with optically thin synchrotron emission.

preprint2008arXiv

AMAZE. I. The evolution of the mass-metallicity relation at z>3

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.

preprint2007arXiv

The evolution of the mass-metallicity relation at z~3

We present preliminary results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z~3 by means of deep near-IR spectroscopy. Gas metallicities and stellar masses are measured for an initial sample of nine star forming galaxies at z~3.3. When compared with previous surveys, the mass-metallicity relation inferred at z~3.3 shows an evolution significantly stronger than observed at lower redshifts. There are also some indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered as the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.3 is difficult to reconcile with the predictions of some hierarchical evolutionary models. We shortly discuss the possible implications of such discrepancies.

preprint2006arXiv

Quasar Luminosity Functions from Joint Evolution of Black Holes and Host Galaxies

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

preprint2006arXiv

The SCUBA HAlf Degree Extragalactic Survey (SHADES) -- II. Submillimetre maps, catalogue and number counts

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

preprint2005arXiv

Can the faint sub-mm galaxies be explained in the Lambda-CDM model?

We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $Λ$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.

preprint2005arXiv

The evolution of actively star forming galaxies in the mid infrared

In this paper we analyze the evolution of actively star forming galaxies in the mid-infrared (MIR). This spectral region, characterized by continuum emission by hot dust and by the presence of strong emission features generally ascribed to polycyclic aromatic hydrocarbon (PAH) molecules, is the most strongly affected by the heating processes associated with star formation and/or active galactic nuclei (AGN). Following the detailed observational characterization of galaxies in the MIR by ISO, we have updated the modelling of this spectral region in our spectro-photometric model GRASIL (Silva et al. 1998). In the diffuse component we have updated the treatment of PAHs according to the model by Li & Draine (2001). As for the dense phase of the ISM associated with the star forming regions, the molecular clouds, we strongly decrease the abundance of PAHs as compared to that in the cirrus, basing on the observational evidences of the lack or weakness of PAH bands close to the newly formed stars, possibly due to the destruction of the molecules in strong UV fields. The robustness of the model is checked by fitting near infrared to radio broad band spectra and the corresponding detailed MIR spectra of a large sample of galaxies (Lu et al. 2003), at once. With this model, we have analyzed the larger sample of actively star forming galaxies by Dale et al. (2000). We show that the observed trends of galaxies in the ISO-IRAS-Radio color-color plots can be interpreted in terms of different evolutionary phases of star formation activity, and the consequent different dominance in the spectral energy distribution (SED) of the diffuse or dense phase of the ISM.