\section{Introduction}
Cosmology is entering in a precision era where the interplay with
particle physics is becoming more and more important.  A noteworthy
example is the effective number of light degrees of freedom
appea\-ring in the different extensions of the Standard Model (SM).
Any extra radiation component in the Universe is usually parametrized,
independently of its statistics, in terms of an effective number of
neutrino species, $N_\textrm{eff}=N_\textrm{eff}^\textrm{SM}+\Delta
N_{\textrm{eff}}$ \cite{Steigman:1977kc}, where
$N_\textrm{eff}^\textrm{SM}$ stands for the number of active neutrinos
in the SM\footnote{In the standard cosmological model with three
neutrino flavors and zero chemical potential we have
$N_\textrm{eff}^\textrm{SM}=3$ at BBN, and
$N_\textrm{eff}^\textrm{SM}=3.046$ at CMB. The small excess in the
last case with respect to the LEP result  \cite{Beringer:1900zz} is
due to the entropy transfer between neutrino species and the thermal
bath during  electron-positron annihilation
\cite{Lesgourgues:2006nd}.}.

The strongest constraints on the effective number of neutrino species
come from Big Bang Nucleosynthesis (BBN). A non-standard value of
$N_\textrm{eff}$ increases the expansion rate, which  results on an
enhancement of the primordial helium abundance. Assuming zero lepton
asymmetry, the number of effective degrees of freedom at BBN turns out
to be $N_\textrm{eff}=3.71^{+0.47}_{-0.45}$ ($68 \%$ C.L.)
\cite{Steigman:2012ve}.  Note that, although the existence of extra
species is somehow favored, the obtained value is still compatible
with the SM prediction within the $95\%$ C.L.

Some constraints on $N_\textrm{eff}$ can be also obtained from the
analysis of the Cosmic Microwave Background (CMB), although the
current limits are significantly weaker than those of BBN. The
combined analysis of WMAP7 results,  Hubble constant measurements and
baryon acoustic oscillations \cite{Komatsu:2010fb} provides a value
$N_\textrm{eff}=4.34^{+0.86}_{-0.88}$ ($68 \%$ C.L.). Similar and
complementary results for smaller CMB scales have been also reported
by the Atacama Cosmo\-logy Telescope \cite{Dunkley:2010ge} and the
South Pole Telescope \cite{Keisler:2011aw}.  It is interesting to
notice the dependence of the e\-ffective number of neutrino species on
the priors considered in the different Bayesian analysis existing in
the literature. While in some references the SM value, $N^{\rm
SM}_{\textrm{eff}}=3$, is ruled out at 95\% C.L.
\cite{GonzalezGarcia:2010un,Hou:2011ec},  in others, such as
\cite{Reid:2009nq}, it is not. Besides, if the helium abundance
obtained from CMB measures is taken into account, together with the
most precise primordial deuterium abundance \cite{Pettini:2012ph}, the
BBN result becomes perfectly consistent with the SM one at the
$2\sigma$ level, $N_\textrm{eff}=3.22\pm 0.55$ \cite{Steigman:2012ve}.
The number of extra degrees of freedom is therefore an open question
to be solved by the Planck satellite, which is expected to determine $
N_\textrm{eff}$ with an accuracy of $\sim 0.3$ at $2\sigma$
\cite{Ichikawa:2008pz}, breaking thereby the degeneracies with nonzero
neutrino masses and dynamical dark energy \cite{Hamann:2010pw}.

In order to account for the apparent radiation excess one can consider
several possibilities. It could be, for instance, the indication of an
extra sterile neutrino \cite{Hamann:2010bk,Hannestad:2012ky}, of relic
gravitational waves \cite{Giovannini:2002qw}, or  arise from other
exotic possibilities such as a decaying particle
\cite{Ichikawa:2007jv,Zhang:2007zzh,Fischler:2010xz}, the
interaction between dark energy and dark matter \cite{Bjaelde:2012wi} or the reheating of the
 neutrino thermal bath \cite{Boehm:2012gr}.
In this Le\-tter we will consider a different possibility within the
minimalistic framework of Higgs-Dilaton cosmology
\cite{Shaposhnikov:2008xb,GarciaBellido:2011de,Shaposhnikov:2008xi}.
This  constitutes an extension of the Higgs inflation idea
\cite{Bezrukov:2007ep}, where the Standard Model Higgs doublet $H$ is
non-minimally coupled to gravity.  The novel ingredient  of
Higgs-Dilaton cosmology is the invariance of the action under scale
transformations. This extra symmetry leads to the absence of any
dimensional parameters or scales\footnote{In particular it forbids the
appearance of a cosmological constant term in the action. In
Higgs-Dilaton cosmology,  the late dark energy dominated period of the
Universe is reco\-vered, at the level of the equations of motion, by
repla\-cing General Relativity with Unimodular Gravity. However, both
the inflationary and prehea\-ting stages considered in this Letter take
place in field space regions where the dark energy contribution is
completely ne\-gligible. We will thus omit  this point here. The
reader is referred to Ref. \cite{GarciaBellido:2011de} for details
about the phenomenological consequences of  Unimodular Gravity in the
Higgs-Dilaton scenario.}. The simplest phenomenologically viable
theory of this kind requires the existence of a new scalar singlet
under the SM gauge group  \cite{Shaposhnikov:2008xb}, the dilaton
$\chi$, non-minimally coupled to gravity. It corresponds to the
Goldstone boson associated with the spontaneous symmetry breaking of
scale invariance and it is therefore massless. This property makes it
a potential candidate for contributing to the effective number of
relativistic degrees of freedom at BBN and recombination. Indeed, this
cosmological test seems to be the only available probe for determining
the existence of the dilaton particle. The coupling between the dilaton
and all the SM fields (apart from the Higgs) is forbidden by quantum
numbers, which, together with the Goldstone boson nature of this
particle, excludes the possibility of a direct detection in an
accelerator experiment \cite{Shaposhnikov:2008xb}.

In this Letter we study the (p)reheating stage in Higgs-Dilaton
cosmology, paying special attention to the dilaton production. This
Letter is organized as follows. In Sections \ref{sec:model} and
\ref{sec:gaugesector} we review the Higgs-Dilaton inflationary model
and show that, given the hierarchical structure of the non-minimal
couplings to gravity, the production of SM particles takes place, up
to some small corrections, as in the simplest Higgs inflationary
scenario \cite{Bezrukov:2008ut,GarciaBellido:2008ab,Bezrukov:2011sz}.
The difference between the two models is described in Section
\ref{sec:dilatonproc}, where we compute the dilaton production,
compare it with the total energy density of SM particles at the end of
the preheating stage and determine its contribution to the effective
number of relativistic degrees of freedom. The  conclusions are
presented in Section \ref{sec:conclusions}.

\section{Higgs-Dilaton inflation}\label{sec:model}

We start by reviewing the Higgs-Dilaton model
\cite{Shaposhnikov:2008xb,GarciaBellido:2011de}. In the unitary gauge
$H^T=(0,h/\sqrt{2})$, it is described by the following Lagrangian
density
\be{general-theory}
\frac{\mathscr L}{\sqrt{-g}}=
\frac{1}{2}(\xi_\chi \chi^2 +\xi_h h^2)R-
\frac{1}{2}\left(\partial \chi\right)^2-U(\chi,h) \ ,
\end{equation}
where we have omitted the part of the SM Lagrangian not involving the
Higgs potential, $\mathscr L_{\text{SM}[\lambda\rightarrow0]}$. The
values of the non-minimal couplings  to gravity can be determined from
CMB observations and turn out to be highly hierarchical ($\xi_\chi\sim
10^{-3}, \xi_h\sim 10^3-10^5$)
\cite{GarciaBellido:2011de,Bezrukov:2009db}. The scale-invariant
potential $U(\chi,h)$ is given by
\be{jord-pot}
U(\chi,h)=\frac{\lambda}{4}
\left(h^2-\frac{\alpha}{\lambda}\chi^2 \right)^2+\beta \chi^4\,,
\end{equation}
with $\lambda$ the self-coupling of the Higgs field. The parameters
$\alpha$ and $\beta$ must be properly tuned in order to reproduce the
correct hierarchy between the  electroweak, Planck and cosmological
constant scales.  In particular, we must require  $\beta\ll\alpha \ll
1$. The smallness of all the couplings involving the dilaton field
gives rise to an approximate shift symmetry
$\chi\rightarrow\chi+\textrm{const.}$, which, as described in Ref.
\cite{Bezrukov:2010jz}, has important consequences for the analysis
of quantum e\-ffects. For the typical energy scales involved in the
(p)reheating stage we can safely set $\alpha=\beta=0$ in all the
following developments.

We study here the (p)reheating of the universe after Higgs-Dilaton
inflation. As emphasized in Ref. \cite{GarciaBellido:2011de}, particle
production is more easily a\-nalyzed in the so-called Einstein-frame,
where the Higgs and dilaton fields are minimally coupled to gravity.
Performing a conformal redefinition of the metric,  $\tilde
g_{\mu\nu}=\Omega^2g_{\mu\nu}$, with conformal factor
$\Omega^2=M_P^{-2}(\xi_\chi \chi^2 +\xi_h h^2)$, we obtain
\be{einst-theory}
\frac{\mathscr L}{\sqrt{-\tilde g}}=
\frac{M_P^2}{2}\tilde R-\frac{1}{2}\tilde K(\chi,h)-\tilde U(\chi,h) \ .
\end{equation}
Here $\tilde K(\chi,h)$ is a non-canonical kinetic term in the basis
$(\phi^1, \phi^2)=(\chi,h)$
\be{kinetic}
\tilde K(\chi,h) = \frac{\kappa_{ij}}{\Omega^2}
\tilde g^{\mu\nu}\partial_\mu\phi^i\partial_\nu\phi^j\,,
\end{equation}
with
\be{kappa} \kappa_{ij}=\left(\delta_{ij}+
\frac{3}{2}M_P^2\frac{\partial_i\Omega^2\partial_j\Omega^2}
{\Omega^2} \right) \,,
\end{equation}
and $\tilde U(\chi,h)\equiv U(\chi,h)/\Omega^4$ is the Einstein-frame
potential.  In order to diagonalize the kinetic term we can make use
of the conserved Noether's current a\-ssociated to scale invariance. It
can be easily shown, via the homogeneous Friedmann and Klein-Gordon equations in
the slow-roll approximation, that the field combination
$(1+6\xi_\chi)\chi^2+(1+6\xi_h)h^2$ is time-independent in the absence
of any explicit symmetry breaking term \cite{GarciaBellido:2011de}.
This conservation suggests a field redefinition to polar variables in
the $(h,\chi)$ plane
\be{var-rho}
r= \frac{M_P}{2}\log\left[\frac{(1+6\xi_\chi)\chi^2+
(1+6\xi_h)h^2}{M_P^2} \right] \,,
\end{equation}
\be{var-theta}
\tan\theta=\sqrt{\frac{1+6\xi_h}{1+6\xi_\chi}}\frac{h}{\chi} \,.
\end{equation}
In terms of the new coordinates, the kinetic term \eqref{kinetic} becomes dia\-gonal,
although non-canonical,
\ba{polar-kin}
\tilde K &=& \left( \frac{1+6\xi_h}{\xi_h}\right)\frac{1}{\sin^2\theta+\varsigma\cos^2\theta}
(\partial r)^2 \nonumber \\&+&\frac{M_P^2 \ \varsigma}{\xi_\chi}
\frac{\tan^2\theta+\mu}{\cos^2\theta(\tan^2\theta+\varsigma)^2}(\partial \theta)^2 \, ,
\end{eqnarray}
where we have defined
\be{var-defs1}
\mu=\frac{\xi_\chi}{\xi_h} \ \
\ \text{and} \ \ \  \varsigma=\frac{(1+6\xi_h)\xi_\chi}{(1+6\xi_\chi)\xi_h} \ .
\end{equation}
The dilatonic field $r$ is massless, as corresponds to the Goldstone
boson associated with the spontaneously broken scale symmetry. The
inflationary potential depends only on the angular variable $\theta$
and it is symmetric around $\theta=0$
\be{polar-pot}
\tilde U(\theta)= \frac{\lambda M_P^4}{4\xi_h^2}
\left(\frac{\sin^2\theta}
{\sin^2\theta+\varsigma\cos^2\theta}\right)^2 \ .
\end{equation}
It can be easily seen that during the (p)reheating stage the values of
the oscillating field $\theta$ are much larger than the $\mu$
parameter for a large number of oscillations,  $\tan^2\theta\gg\mu$.
This allows us  to neglect the $\mu$ term in Eq. \eqref{polar-kin} and
perform an extra field redefinition
\be{var-r-theta}
\rho= \gamma^{-1}r\,,\ \ \ \
\vert\phi\vert= \phi_0-\frac{M_P}{a}\tanh^{-1}\left[\sqrt{1-\varsigma}\cos\theta\ \right] \,,
\end{equation}
with
\be{var-defs2}
\gamma= \sqrt{\frac{\xi_\chi}{1+6\xi_\chi}}
\ \ \ \text{and} \ \ \ a=\sqrt{\frac{\xi_\chi(1-\varsigma)}{\varsigma}} \ .
\end{equation}
The variable $\phi$ is periodic and defined in the compact interval
$\phi \in\left[-\phi_0,\phi_0\right]$, where $\phi_0=
M_P/a\tanh^{-1}\left[\sqrt{1-\varsigma}\ \right] $ corresponds to the
value of the  field $\phi$ at the beginning of inflation. As happened
in Higgs inflation  \cite{Bezrukov:2008ut,GarciaBellido:2008ab}, the
absolute value in Eq.
\eqref{var-r-theta} is required for $\phi$ to maintain the symmetry of the initial $\theta$ field around the minimum of the potential.
In terms of these variables the Lagrangian \eqref{einst-theory} takes a
very simple form
\be{angul-theory}
\frac{\mathscr L}{\sqrt{-\tilde g}}=\frac{M_P^2}{2}\tilde R -
\frac{e^{2b\left(\phi\right)}}{2} (\partial \rho)^2-
\frac{1}{2}(\partial \phi)^2-\tilde U(\phi)\ ,
\end{equation}
that has been widely studied in the literature \cite{GarciaBellido:1995fz,GarciaBellido:1995qq,DiMarco:2002eb,Lachapelle:2008sy}. In our case, the coefficient in front of the dilaton kinetic term is given by
\be{exp2b}
e^{2b\left(\phi\right)}\equiv\varsigma\cosh^2\left[a \kappa\left(\phi_0-\vert\phi\vert \right)\right]\,,
\end{equation}
with $\kappa$ the inverse of the reduced Planck mass $M_P$.  The inflationary potential \eqref{polar-pot} becomes
\be{defin-pot}
\tilde U(\phi)=\frac{\lambda M_P^4}{4\xi_h^2 (1-\varsigma)^2}
\left(1- \varsigma\cosh^2[a\kappa\left(\phi_0-\vert\phi\vert\right)]\right)^2\,.
\end{equation}
\begin{figure}
\centering
\includegraphics[scale=0.7]{HDcomparisonH.eps}
\caption{Comparison between the Higgs-Dilaton inflatio\-nary potential
(blue continuous line) obtained from Eq.\eqref{defin-pot} and the
corresponding one for the  Higgs inflation mode (red dotted line). In
spite of the slight differences in the upper inflationary region, they
nicely agree in the lower part, where  the (p)reheating stage takes
place. Here $U_0=\lambda M_P^4/(4\xi_h^2)$ and $\kappa=M_P^{-1}$.}
 \label{comparisonHDH}
\end{figure}
\section{SM particle production}\label{sec:gaugesector}
 As shown in Fig.~\ref{comparisonHDH}, the shape of the Higgs-Dilaton
potential \eqref{defin-pot}  clearly resembles that of the simplest
Higgs inflationary scenario. In spite of the slight differences, both
of them present an exponentially flat region for large field values
and nicely agree for small ones. Indeed, the relation between them
becomes explicit if we approximate Eq. (\ref{defin-pot}) for small
$\phi$. The potential around the minimum behaves, in a good
approximation, as the standard chaotic potential\footnote{It is
interesting to notice at this point that this statement is only valid
for values of the field such that $\tan^2\theta \gg \mu$. For very
small values of the angular variable $\phi$ we recover the standard
$\lambda \phi^4$ Higgs potential. As it happens in Higgs inflation
\cite{Bezrukov:2008ut,GarciaBellido:2008ab}, this region turns out to
be extremely small, being completely irrelevant for the study of the
(p)reheating stage.}
 \begin{equation}\label{HDpot}
\tilde U(\phi)\simeq \frac{1}{2} M_{HD}^2 \phi^2 + {\cal O}(\phi^3)\,,
\end{equation}
where the higher order corrections can be safely neglected after a few
oscillations.  The curvature of the potential,
$M_{HD}^2=\left(1+6\xi_\chi\right)M^2$, coincides, up to small
corrections, with that of Higgs inflation $M^2=\frac{\lambda
M_P^2}{3\xi_h^2}$. The same applies to the masses of gauge bosons and
fermions.  We obtain  $ (\tilde m_{{\cal A},f}^2)_{HD}=\tilde
m^2_{{\cal A},f} \left(1+6 \xi_\chi\right)$, with $\tilde m^2_{{\cal
A},f} \simeq \frac{\alpha g^2 M_P}{4\xi_h }|\phi| $ the Einstein-frame
gauge boson and fermion masses in the Higgs inflation model
\cite{Bezrukov:2008ut,GarciaBellido:2008ab}. Here  $g = g_2,
g_2/\cos\theta_w$ and $\sqrt{2}y_f$, for ${\cal A}=W,Z$ bosons and
fermions $f$, respectively.  We see therefore that,  from the point of
view of (p)reheating, all the relevant physical scales in Higgs and
Higgs-Dilaton inflation coincide, up to small corrections proportional
to the small parameter $\xi_\chi$. This allows us to apply the main
results of Refs. \cite{Bezrukov:2008ut,GarciaBellido:2008ab} to the
Higgs-Dilaton case. Let us summarize here those results. In Higgs
inflation, the SM particles are produced through the so-called
Combined Preheating mechanism
~\cite{Bezrukov:2008ut,GarciaBellido:2008ab}. Intermediate $W^{\pm}$
and $Z$ bosons are created from the oscillations of the Higgs at the
bottom of the potential~\eqref{HDpot}, whenever there is a violation
of the adiabaticity condition. While there is no restriction on the
number of created gauge bosons, the direct production of SM fermions
by this mechanism is severely res\-tricted by Fermi-Dirac statistics.
The SM fermions appear as se\-condary products of the weak bosons
created in each zero crossing.  Once produced, the gauge bosons
acquire a large effective mass due to the increasing expectation value
of the Higgs field and decay perturbatively into quarks and leptons
within a semioscillation of the Higgs field. This decay rapidly
depletes their occupation numbers and postpones the development of
parametric re\-sonance. During the first oscillations, the fraction of
energy into SM particles is still very small compared with the energy
in the oscillating Higgs field. A large number of oscillations ($
t\sim 300 \,M_{HD}^{-1}$) will be needed in order to transfer a
significant amount of energy into the  SM bosons and fermions. The
decreasing of the Higgs amplitude due to the expansion of the Universe
eventually reduces the decay rate and parametric resonance becomes the
dominant effect.  At this point, the gauge bosons start to build up
their occupation numbers via bosonic stimulation and reheating occurs
within a few osci\-llations. Soon afterwards, the Universe is filled
with the remnant Higgs condensate and a non-thermal distribution of
fermions and bosons, redshifting as radiation and matter respectively.
From there on until thermalization, the evolution of the system is
highly non-linear and non-perturbative, which makes difficult to make
a clear statement about the subsequent evolution without the use of
numerical lattice
simulations~\cite{Bezrukov:2008ut,GarciaBellido:2008ab} .  However,
thermal equili\-brium is expected to be achieved at a rehea\-ting
temperature $T_0\sim (3-15)\times10^{13}$ GeV, much above the QCD phase
transition scale, $T_\textrm{QCD}\sim 300$ MeV, due to the large SM
couplings  \cite{Bezrukov:2008ut,GarciaBellido:2008ab}.

\section{Dilaton production}\label{sec:dilatonproc}
In addition to the SM fields, the Higgs-Dilaton inflationary scenario
incorporates an extra degree of freedom, the dilaton field $\rho$. The
constancy of the classical background  component  is of course
gua\-ranteed by the scale invariance current conservation, but this
reasoning does not apply to the corresponding quantum excitations. As
suggested in Ref. \cite{Lachapelle:2008sy}, these modes can be excited
in the preheating stage after inflation\footnote{Any dilaton
production previous  to this stage is completely diluted by the
inflationary expansion.} through the non-canonical kinetic term in the
Einstein-frame Lagrangian \eqref{angul-theory}, which mixes quantum
excitations and background solutions. Although the perturbative
dilaton production through this mixing is expected to be very small,
non-perturbative effects might play an important role
\cite{Lachapelle:2008sy}. In this section we estimate the energy
density residing in the dilaton field at the end of the preheating
stage.  Let us start by considering the linearized\footnote{ The
coupled Higgs-Dilaton equations further simplify since the Higgs
fluctuations are not significantly
amplified~\cite{GarciaBellido:2008ab}, and thus can be treated as
decoupled equations.} equations of motion for dilaton perturbations
$\delta \rho_k$ in Fourier space

 \begin{equation}\label{perturbations}
\delta  \ddot\rho_k+\left(3H+2 \dot b\right)\delta\dot\rho_k+\frac{k^2}{a^2}\delta\rho_k=0\,,
 \end{equation}
where we have ignored metric perturbations and taken into account the
constancy of the background field $\rho$  during the Higgs
oscillations at the end of inflation. The function $b=b(\phi)$  plays
the role of an additional oscillatory damping term for the dilaton
perturbations and depends on the absolute value of the inflaton field
$\phi$, cf. Eq. \eqref{exp2b}. This dependence makes it cumbersome
the direct application of the techniques presented in Ref.
\cite{Lachapelle:2008sy} for the study of particle creation in models
with non-canonical kinetic terms. The field redefinitions used there
would imply delta functions coming from the derivatives of the
absolute value, which substantially complicates the analytic and
numerical treatment of the problem. On the other hand, although the
rephrasing of Eq. \eqref{perturbations} as a Hill's equation is
extremely useful for the understanding of the particle creation
mechanisms, it is not necessary for a precise computation in an
expanding background, as the one needed to estimate $N_\textrm{eff}$.
For this reason, we will adopt an alternative approach, dealing only
with non-singular evolution equations in an expanding Universe. Let us
rewrite Eq. \eqref{perturbations} as
 \begin{equation}\label{perturbations2}
\frac{1}{a^3 e^{2b}}\frac{d}{dt}\left(a^3 e^{2b}\frac{d}{dt}\delta\rho_k\right)+\frac{k^2}{a^2}\delta\rho_k=0\,,
 \end{equation}
which, after a redefinition of time, $d\tau=a^{-3} e^{-2b}dt$, can be
recast in the form of a time-dependent harmonic oscillator
 \be{dilatonconf}
\delta\rho''_k+\omega_k^2(\tau)\delta\rho_k=0\,,
 \end{equation}
with frequency $\omega_k^2(\tau)=k^2 a^4 e^{4b}$. Here the prime
denotes derivative with respect to the new time $\tau$.  Choosing an
initial vacuum state with zero particle content\footnote{This corresponds to the initial vacuum  conditions $\delta\rho_k(0)=1/\sqrt{2\omega_k}$ and $\delta\rho'_k(0)=-i\omega_k\delta\rho_k$. } , the number of created
dilatons  is given by
\be{nk}
 n_k+\frac{1}{2}=\frac{1}{2\omega_k}\Big(\vert\delta\rho'_k\vert^2+ \omega_k^2 \vert\delta \rho_k\vert^2\Big)\,,
 \end{equation}
 and its associated energy density,
\be{rhok}
 \hspace{15mm} \rho_\chi=\frac{1}{2\pi^2}\int_0^\infty dk k^2\,\omega_k n_k\,,
 \end{equation}
can be computed numerically by solving Eq. \eqref{dilatonconf},
together with the background evolution equations. The resulting energy
density must be compared with  the energy density in SM particles at
the end of the preheating stage,
$C\equiv\rho_\chi^0/\rho^0_\textrm{SM}$. This quantity can be easily
related to the effective number of light degrees of freedom
$N_\textrm{eff}$. In order to do that, let us note that, once produced, the
dilaton particles are completely decoupled from the SM particles,
being its energy density only diluted by the expansion of the Universe
$\rho_\chi^0 a_{\textrm{0}}^4=\rho_\chi a_{\textrm{f}}^4$.  Here the
subscripts '0' and 'f' stand for the end of the preheating stage and
the BBN epoch respectively. On the other hand, the total entropy of SM
particles after thermalization remains constant,  $s_{\textrm{0}}
a^3_{\textrm{0}}=s_{\textrm{f}} a^3_{\textrm{f}}$. Taking into account
the relation between the entropy density and the number of
relativistic degrees of freedom $g$ at a given temperature,
$s=\frac{\pi^2}{30}g T^3$, the previous expression can be rewritten as
$g_0 T_0^3 a_0^3= g_f T_f^3 a_f^3$.  By combining this expression with the
evolution equation for the dilaton energy density described above and
dividing the result by the energy density stored in a single neutrino
species, $\rho_{\nu}=\frac{\pi^2}{30}g_{\nu}T_{\textrm{f}}^4$, we get
 \be{Neff}
 \Delta N_\textrm{eff}\equiv\left(\frac{\rho_\chi}{\rho_\nu}\right)_\textrm{f}=\frac{g_0}{g_\nu}\left(\frac{g_f}{g_0}\right)^{4/3}C\simeq 2.85\, C\,.
 \end{equation}
In the last equality we have made use of the number of SM degrees of
freedom at the end of inflation ($g_{\textrm{0}}=106.75$) and at BBN
($g_{\textrm{f}}=10.75$). Therefore, we see that, in order to have a
contribution to the effective number of relativistic degrees of
freedom within the reach of the Planck satellite, roughly a $10 \%$ of
the energy density at the end of inflation must be converted into
dilatons. Nevertheless, the transferred fraction turns out to be
signi\-ficantly smaller. Evaluating the numerical solution of
Eq.~\eqref{dilatonconf} at the time at which the energy density in SM
particles roughly equals the initial energy density of the inflaton
field, $\tau(t_0 \sim 300 \,M_{HD}^{-1})$
\cite{GarciaBellido:2008ab}, we get $C\sim 10^{-7}$.  The precise
value of the $C$ parameter weakly depends on the ratio
$\xi_h/\sqrt{\lambda}$, which determines the total energy density
available at the end of inflation \cite{GarciaBellido:2011de}, and it is
quite insensitive to the particular value of the small non-minimal coupling
 $\xi_\chi$.

Although the non-perturbative creation of dilatons due
to the background field $\phi$ turns out to be extremely small, one should consider
 the possi\-bility of producing them as secondary products of the Higgs particles
  created at the preheating and thermalization stages. The Lagrangian
  \eqref{angul-theory} leads  to a number of perturbative processes, such as the decay of Higgs particles
into dilatons ($\phi\to\rho\rho$) or Higgs-Higgs scatterings ($\phi\phi\to\rho\rho$). The
corresponding decay rate and cross-section are given respectively by
\be{gammasigma}
\Gamma \simeq \frac{M_{HD}^3}{192\pi M_P^2} \ \ \ \ \text{and} \ \ \ \
\sigma \simeq \frac{E^2}{576\pi M_P^4}~,
\end{equation}
where $M_\textrm{HD}^2$ is the effective Higgs mass in the region  $M_P/\xi_h < \phi \ll M_P$ (cf. Eq.~\eqref{HDpot}) and $E$ is the Higgs energy in the center-of-mass frame. Assuming this energy to be of the order of the temperature of the thermalized SM plasma, $T_0$, it can be easily seen that the contribution of these processes to the $C$ parameter is of order
\be{cpert}
C \propto \left(\frac{T_0}{M_P}\right)^3\,,
\end{equation}
and therefore much smaller than the non-perturbative contribution found above.

\section{Conclusions}\label{sec:conclusions}
We have considered particle production in a scale-invariant extension
of Higgs inflation known as Higgs-Dilaton inflation. This model
predicts the existence of an extra massless particle -  the dilaton -
which might contribute to the effective number of light degrees of
freedom. After recasting the problem in the appropriate set of
varia\-bles, all the particle masses and energy scales in the model
turn out to coincide, up to small corrections, with those of Higgs
inflation. Gauge bosons and fermions are therefore produced through
the so-called Combined Preheating mechanism.  On the other hand, the
production of dilatons could take place only through the non-canonical
kinetic term for the dilaton field, fed by the absolute value of the
Higgs field. We have shown that the number of non-perturbatively
created particles can be easily eva\-luated thanks to a particular time
redefinition. The dilaton energy density computed this way turns out
to be extremely small, which is translated into an effective number of
relativistic degrees of freedom  very close to the SM one. This result
is not modified by any of the subsequent perturbative processes
involving the Higgs particles produced at the bottom of the potential.
We thus conclude that, in spite of the potential value of BBN and CMB
for testing the existence of the elusive dilaton particle, its
subdominant production at the (p)reheating stage after inflation
makes  it completely undetectable by any particle physics experiment
or cosmological observation. The only remnant of the dilaton field is
a dynamical Dark Energy stage with an equation of state close, but slightly
different, to that of a cosmological constant, which leads to a
power-like expansion of the Universe in the far future
\cite{Shaposhnikov:2008xb}.

\section*{Acknowledgements}
We thank Daniel Zenh\"ausern for collaboration in the early stages of this project and Julien Lesgourgues for useful discussions.
This work was partially supported by the Swiss National Science
Foundation. We also acknowledge  financial support from the Madrid
Regional Government (CAM) under the program  HEPHACOS
S2009/ESP-1473-02, from MICINN under grant  AYA2009-13936-C06-06  and
Consolider-Ingenio 2010 PAU (CSD2007-00060), as well as from the
European  Union Marie Curie Initial Training Network "UNILHC"
PITN-GA-2009-237920.


\section{Introduction}
Ethylene oxide (c-C$_{2}$H$_{4}$O) is the simplest epoxide, a ring-shaped organic compound involving an oxygen atom bonded to two carbon atoms.
Since it is believed to lead to the synthesis of amino acids (Cleaves 2003) and hence potentially influence early pre-biotic pathways (Miller \& Schlesinger 1993), c-C$_{2}$H$_{4}$O may play an important role in the interstellar medium (ISM).

First detected towards the Sgr B2N molecular cloud associated with 10 rotational transitions (Dickens et al. 1997), this molecule has been investigated in a more extended frequency range (662 lines between 262-358 GHz)(Pan et al. 1998) and more recently Bernstein \& Lynch (2009) attributed the unidentified IR band features discovered by Gillett et al. (1973) to ethylene oxide and to another small carbonaceous molecule, cyclopropenylidene. Ikeda et al. (2001) estimated the c-C$_{2}$H$_{4}$O fractional abundances by observing 20 massive star-forming regions and two dark clouds. In fact, ethylene oxide was detected towards 10 of these sources but it was not found in TMC-1 and TMC-1(NH$_{3}$) dark clouds, in contrast to one of its isomers, acetaldehyde (CH$_{3}$CHO) which seems to be ubiquitous. Ikeda et al. (2001) also investigated the relationship between c-C$_{2}$H$_{4}$O abundances and the dust temperature, but without finding any correlation. On the other hand, the evidence of similar trends in the column densities of all C$_{2}$H$_{4}$O isomers and ethanol identified the latter as a possible precursor of ethylene oxide and its isomers. Indeed, a gas-phase pathway for ethylene oxide formation has been proposed (Dickens et al. 1997) with ethanol as one of the reactants:

CH$_{3}$$^{+}$ + C$_{2}$H$_{5}$OH $\rightarrow$ C$_{2}$H$_{5}$O$^{+}$ + CH$_{4}$

C$_{2}$H$_{5}$O$^{+}$ + e$^{-}$ $\rightarrow$ C$_{2}$H$_{4}$O + H

where the formula C$_{2}$H$_{4}$O does not distinguish among the different isomers.
However, Ikeda et al. (2001) performed simulations of c-C$_{2}$H$_{4}$O chemical evolution by including only pure gas-phase reactions at different temperatures and volume densities and they found that the calculated abundances of ethylene oxide decreased when the gas-kinetic temperature was raised. The results of Ikeda et al. (2001) suggested that other reactions such as grain surface reactions may be important for the synthesis of ethylene oxide.

Another gas-phase pathway for the formation of ethylene oxide has been proposed by Turner \& Apponi (2001) which studied the reaction chain to produce c-C$_{2}$H$_{4}$O and one of its isomers, vinyl alcohol (CH$_{2}$CHOH), starting with the interaction between CH$_{3}$$^{+}$ and H$_{2}$CO; unfortunately, none of the intermediates of the latter reaction have been identified yet in the ISM. The abundance ratios among the three isomers, acetaldehyde, ethylene oxide and vinyl alcohol, seem to be very uncertain (Bennet et al. 2005), but generally CH$_{3}$CHO was found to be the most abundant and CH$_{2}$CHOH has been predicted to have a lower abundance than ethylene oxide in the experiments by Bennet et al. (2005). This comparison could be explained in the light of the fact that the formation of ethylene oxide occurs via direct addition of the oxygen across the C-C double bond of ethylene, while vinyl alcohol requires an insertion of an oxygen into the C-H bond. Hudson \& Moore (2003) suggested a solid state route for the formation of CH$_{2}$CHOH involving exposure of H$_{2}$O:C$_{2}$H$_{2}$ ices to protons and photons, but the first proposed surface scheme for ethylene oxide formation dates back to Charnley (2004) who identified acetaldehyde (CH$_{3}$CHO) as a possible precursor.
The contribution of grain-surface reactions to the production of the three isomers with formula  C$_{2}$H$_{4}$O has also been investigated by Bennet et al. (2005): they presented the results of electron irradiation on a mix of carbon dioxide and ethylene ice, under interstellar and cometary conditions, to highlight the effect of non-thermal processing of ices.

Recently, an alternative solid state pathway for the formation of ethylene oxide has been investigated by Ward \& Price (2011) who studied the surface reactivity between C$_{2}$H$_{4}$ and thermal O atoms under astrophysically relevant conditions. They co-deposited the reactants on a graphite surface at various temperatures in the 12-90 K range. In order to identify the products formed, they recorded desorption spectra and they also selectively ionized the three different isomers with a tunable dye laser. Finally, they performed a similar laboratory study with C$_{3}$H$_{6}$ as reactant, but unfortunately although C$_{3}$H$_{6}$O was observed they were unable to identify the isomeric products of the reaction due to the large number of possible molecules with formula C$_{3}$H$_{6}$O. Despite these limitations, the authors assumed that the majority of the C$_{3}$H$_{6}$O formed is propylene oxide by following the same mechanism as for the case of ethylene oxide production. Unlike c-C$_{2}$H$_{4}$O, c-C$_{3}$H$_{6}$O has not yet been detected in the ISM since it is believed to isomerize to propanal or acetone, which have both been observed in the ISM (Hollis et al. 2004 and Combes et al. 1987, Snyder et al. 2002).
To date, the observed abundances of ethylene oxide have not yet been reproduced by chemical models; in fact, to our knowledge, none of the grain-surface pathways for formation for ethylene oxide have been included in the models that simulate the star formation process.

In this paper, we update the UCL\_CHEM gas-grain chemical model with the reaction rates for the production of this species as evaluated by Ward \& Price (2011). By following the same logic as these authors, we extended our study to include the formation of propylene oxide in order to verify if C$_{3}$H$_{6}$O can be produced under ISM conditions even without having any observational evidence of its presence in the interstellar space.
The paper is organized as follows: Section 2 gives details of the model; Section 3 considers the astronomical implications of our study; Section 4 contains a discussion of our results and it also shows a comparison between our theoretical column densities and those taken from the observations. Finally, we give our conclusions in Section 5.

\section{The Chemical Model}
UCL\_CHEM is a time-depth dependent gas-grain chemical model which simulates the abundances of molecules in the ISM and star-forming regions.
The first version of the code was developed by Viti \& Williams (1999) and revisited by the same authors in 2004 (Viti et al. 2004). Our model has been extensively described in the literature (see Occhiogrosso et al. 2011 for more details) so here we only highlight the main points.

The process of star formation is simulated in two steps: during Phase I the gas collapses to a dense core from a diffuse and atomic state; at this stage the gas freezes onto the  grains and surface reactions occur.
During Phase II, the gas and the dust surrounding the central star heat up, causing the sublimation of the icy mantles. The thermal desorption can be set up to occur instantaneously or as a function of time as in Viti et al. (2004). In the latter case, molecules are classified in five different categories as described by the experimental study of Collings et al. (2004): (i) CO-like; (ii) H$_{2}$O-like; (iii) intermediate; (iv) reactive and (v) refractory. In the present work, we consider c-C$_{2}$H$_{4}$O as H$_{2}$O-like (it will co-desorb when the H$_{2}$O ice desorbs) based on the fact that in the Ward \& Price (2011) experiments, its desorption profile shows only one peak (as for the water desorption profile) and its dipole moment is big enough to assume that the species will have a lower desorption energy from the surface (in that specific case, graphite) than from itself. In our model the temperature of the gas is derived as a function of the luminosity (and therefore the age) of the accreting protostar through a simple power law (see Viti et al. 2004 for details).
The model includes two networks: the gas-phase network is based on the UMIST database (Woodall et al. 2000) (even if some coefficients have been modified with those from the most recent database KIDA (Wakelam et al. 2009)) with some updates discussed later (see Section 3), and a grain-surface network with reactions whose rates are experimentally determined where possible.

For this study, we include in our surface reaction network the formation of c-C$_{2}$H$_{4}$O and c-C$_{3}$H$_{6}$O on the grain surface as investigated by Ward \& Price (2011) (see Table 1). In order to analyze their data, these authors developed a simple kinetic model in which molecules can react via Eley-Rideal (ER) or Langmuir-Hinshelwood (LH) mechanisms. In the latter case, both reactants adsorb and thermalize on the surface before a diffusion reaction occurs. We consider the reactions to occur via the LH mechanism that seems of greater astrophysical relevance as stated by Awad et al. (2005) and by the authors themselves (see their paper for the details).

\begin{table}
 \begin{minipage}{126mm}
  \caption{Surface reactions investigated by Ward \& Price (2011)}
  \label{anymode}
  \begin{tabular}{@{}|cc}
  \hline
$^{a}$Reaction & $\alpha$ (cm$^{2}$mol$^{-1}$s$^{-1}$) at 20K\\
   \hline
mO + mC$_{2}$H$_{4}$ $\rightarrow$ mc-C$_{2}$H$_{4}$O & 1.20x10$^{-19}$\\

mO + mC$_{3}$H$_{6}$ $\rightarrow$ mc-C$_{3}$H$_{6}$O & 3.40x10$^{-19}$\\
 \hline
  \end{tabular}
\footnotetext[1]{The m before the molecular formul$\ae$ stands for $mantle$}
 \end{minipage}
\end{table}

Since our code accounts for reactions that occur in three dimensions (i.e. gas-phase reactions), we need to transform the surface rate constants of Ward and Price discussed above. Therefore we first evaluate two quantities: the volume of ice per grain and the number density of grains. In the latter case we assume a dust/mass ratio of 0.01 and a dust density of 2.5 g cm$^{-3}$. Thus, we obtain a grain number density of 3.8x10$^{-5}$ cm$^{-3}$. In order to estimate the volume of ice we consider the grains as spheres and we assume a thick ice layer of 0.3 $\micron$ (Collings et al. 2003). Since the major partner in the ice mixture is water we also account for the fraction of C$_{2}$H$_{4}$ or C$_{3}$H$_{6}$ ice compared to water ice in order to avoid overestimating our rate. We consider that not all the oxygen ice will react with the ethylene or propene ice because of the presence of other ice components in between. We hence estimate a final value of 6x10$^{-15}$ cm$^{3}$s$^{-1}$ at 20 K for the reaction between O and C$_{2}$H$_{4}$ and a rate of 5x10$^{-14}$cm$^{3}$s$^{-1}$ at 20 K for the reaction between O and C$_{3}$H$_{6}$ onto the grains.

In order to reveal the contribution of these rates to the abundance of c-C$_{2}$H$_{4}$O and c-C$_{3}$H$_{6}$O in the ISM, we initially did not include any gas-phase pathway for the formation of these species.

The percentage of depletion of oxygen atoms on the grains is certainly one of the critical parameters in the formation of ethylene and propylene oxide; Caux et al. (1999) found that the fractional abundance of oxygen (relative to hydrogen nuclei) on the grains is of the order of 10$^{-6}$ and this amount depends on the percentage of the oxygen that is depleted in the dense, cold cores. In our model, atomic oxygen freezes on the grains in the form of H$_{2}$O, O and OH. It is known that oxygen containing species in interstellar ices consist mainly of H$_{2}$O and also that the fractional abundances of water ice (relative to the total number of hydrogen nuclei) are of the order of 10$^{-4}$ (Caux et al. 1999); we vary the percentage contribution of each channel by which oxygen containing species accumulate on the grains by varying the efficiency of the different oxygen accretion pathways as shown in Table 2.

\begin{table}
 \begin{minipage}{126mm}
  \caption{List of test chemical models.}
  \label{symbol}
  \begin{tabular}{@{}|cccc}
  \hline
Test & mH$_{2}$O (\%) & mO (\%) & mOH (\%)\\
   \hline
Test1 & 99.5 & 0.5 & 0 \\
Test2 & 90 & 1 & 9 \\
Test3 & 90 & 2 & 8 \\
Test4 & 90 & 5 & 5 \\
Test5 & 90 & 10 & 0 \\
Test6 & 80 & 20 & 0 \\
Test7 & 70 & 15 & 15 \\
   \hline
  \end{tabular}
 \end{minipage}
\end{table}

Since we find 1 \% of oxygen atoms in the ice (mO) is sufficient for the observational values of water ice to be satisfied we use this percentage subsequently for our model.
Figure 1 shows the fractional abundances of our selected species as a function of the time at the end of Phase I for a model with a final density of 10$^{7}$ cm$^{-3}$. During this phase atoms and molecules freeze onto the grains and surface reactions can occur due to the catalytic effects of the icy mantles. The dotted and dashed lines refer to the reactants, atomic oxygen and ethylene respectively, the solid line shows the ethylene oxide formation. Note that the fractional abundance of ethylene at the beginning of Phase I is low since it also mainly forms on the grains.
We mostly focus on ethylene oxide because results show that c-C$_{3}$H$_{6}$O is produced below the threshold of detection; this fact may mean that this molecule is not detectable in the ISM although we stress that we do not have a complete network of reactions (both in gas-phase and on grain surfaces) for c-C$_{3}$H$_{6}$O and its precursors.

\begin{figure}
 \includegraphics[width=60mm, angle=270]{F1.eps}
 \caption{Grain surface formation of selected species during Phase I. The fractional abundances (with respect to the total number of hydrogen nuclei) vary along with the time.}
\end{figure}

We conclude that the reaction between oxygen and ethylene on grains with a reaction rate of 6x10$^{-15}$ cm$^{3}$s$^{-1}$ (evaluated as described above) may be the main formation path for c-C$_{2}$H$_{4}$O.

In order to compare our theoretical results with observations we ran further models where gas-phase reactions were included despite the fact that the routes proposed to date (as mentioned in the Introduction) do not distinguish among the different isomers with formul$\ae$ C$_{2}$H$_{4}$O.
We then include in our model the gas-phase routes of formation for c-C$_{2}$H$_{4}$O listed in Table 3. We insert the C$_{2}$H$_{4}$O species in the UMIST database with the reaction rates taken from the OSU database (http://www.physics.ohio-state.edu/~eric/research.html). Since the literature does not distinguish among the different isomers, we assume the ratios of 2:1.5:1 (as reported by Ikeda et al. 2001) for acetaldehyde, ethylene oxide and vinyl alcohol respectively to evaluate the rate of formation of each product.

\begin{table*}
 \begin{minipage}{126mm}
  \caption{Gas-phase paths of c-C$_{2}$H$_{4}$O formation and destruction. The rates are taken from the OSU database and they have been scaled by considering the ratio 2:1.5:1 as in Ikeda et al. (2001). R indicates a reactant, P a product.}
  \label{symbol}
  \begin{tabular}{@{}|cccccc}
   \hline
R1 & R2 & P1 & P2 & P3 & $k$$\backslash$(cm$^{3}$s$^{-1}$)\\
   \hline
C$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{3}$O$^{+}$ & CH & & 7.3x10$^{-9}$\\
C$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{4}$O$^{+}$ & C & & 7.3x10$^{-9}$\\
H$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{3}$O$^{+}$ & H$_{2}$ & & 2.2x10$^{-10}$\\
H$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{4}$O$^{+}$ & H & & 2.2x10$^{-10}$\\
He$^{+}$ & C$_{2}$H$_{4}$O & HCO$^{+}$ & CH$_{3}$ & He & 6.1x10$^{-9}$\\
He$^{+}$ & C$_{2}$H$_{4}$O & CH$_{3}$$^{+}$ & HCO & He & 6.1x10$^{-9}$\\
He$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{2}$O$^{+}$ & H$_{2}$ & He & 6.1x10$^{-9}$\\
He$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{3}$O$^{+}$ & H & He & 6.1x10$^{-9}$\\
He$^{+}$ & C$_{2}$H$_{5}$OH & C$_{2}$H$_{3}$O$^{+}$ & H$_{2}$ & & 2.9x10$^{-9}$\\
S$_{2}$$^{+}$ & C$_{2}$H$_{5}$OH & H$_{2}$S$_{2}$$^{+}$ & C$_{2}$H$_{4}$O & & 6.4x10$^{-11}$\\
H$_{3}$$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{5}$O$^{+}$ & H$_{2}$ & & 2.7x10$^{-8}$\\
HCO$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{5}$O$^{+}$ & CO & & 1.1x10$^{-8}$\\
H$_{3}$O$^{+}$ & C$_{2}$H$_{4}$O & C$_{2}$H$_{5}$O$^{+}$ & H$_{2}$O & & 1.2x10$^{-8}$\\
O & C$_{2}$H$_{5}$$^{+}$ & C$_{2}$H$_{4}$O$^{+}$ & H & & 7.5x10$^{-11}$\\
CH$_{3}$$^{+}$ & C$_{2}$H$_{4}$O & C$_{3}$H$_{6}$OH$^{+}$ & & & 4.3x10$^{-10}$\\
e$^{-}$ & C$_{2}$H$_{4}$O$^{+}$ & CH$_{3}$ & HCO & & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{4}$O$^{+}$ & C$_{2}$H$_{2}$O & H & H & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{4}$O$^{+}$ & C$_{2}$H$_{2}$O & H$_{2}$ & & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{5}$O$^{+}$ & C$_{2}$H$_{4}$O & H & & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{5}$OH$^{+}$ & C$_{2}$H$_{4}$O & H$_{2}$ & & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{5}$OH$_{2}$$^{+}$ & C$_{2}$H$_{4}$O & H$_{2}$ & H & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{2}$H$_{6}$CO$^{+}$ & C$_{2}$H$_{4}$O & CH$_{2}$ & & 6.5x10$^{-7}$\\
e$^{-}$ & C$_{3}$H$_{6}$OH$^{+}$ & C$_{2}$H$_{4}$O & CH$_{3}$ & & 6.5x10$^{-7}$\\
C$_{2}$H$_{4}$O & photon & CH$_{3}$ & HCO & & 2.2x10$^{-10}$\\
C$_{2}$H$_{4}$O & photon & CH$_{4}$ & CO & & 2.2x10$^{-10}$\\
C$_{2}$H$_{4}$O & photon & C$_{2}$H$_{4}$O$^{+}$ & e$^{-}$ & & 1.5x10$^{-10}$\\
   \hline
  \end{tabular}
 \end{minipage}
\end{table*}

The different paths listed in Table 3 are mostly two-body reactions. Specifically we have charge exchange and ion-neutral reactions, electronic recombination of positive ions with electrons followed by the dissociation of the resulting molecule (dissociative attachment) and finally photo-processes where the dissociation or ionization of neutral species is performed by UV photons from a standard interstellar UV field. In the latter case, the reaction rates depend on the visual extinction of the source so we do not expect it to be efficient at the high visual extinction of hot cores and hot corinos.

Before testing the sensitivity of the c-C$_{2}$H$_{4}$O abundances to the different physical conditions of potential astronomical sources, we investigate another important parameter for the formation of ethylene oxide: the amount of C$_{2}$H$_{4}$ available on grain surfaces. Despite the fact that C$_{2}$H$_{4}$ ice has not been detected to date, several experiments (Hiraoka et al. 2000; Kaiser \& Roessler 1998; Barratta et al. 2003; Bennett et al. 2006) have looked at its formation on grains under physical conditions relevant for the interstellar medium. Using different experimental techniques, C$_{2}$H$_{4}$ was observed as a product following irradiation of CH$_{4}$-containing icy mixtures and C$_{2}$H$_{4}$ was also predicted to be an intermediate in the hydrogenation of C$_{2}$H$_{2}$ (Kaiser \& Roessler 1998; Bennett et al. 2006). In these cases, ethylene formation competes with the formation of the C$_{2}$H$_{6}$ that was found the most abundant molecule in the product mixture. In our model we include the hydrogenation of C$_{2}$H$_{2}$ on grains in Phase I to form C$_{2}$H$_{4}$ as well as C$_{2}$H$_{6}$ (see Table 4).

\begin{table}
 \begin{minipage}{80mm}
  \caption{The main accretion and hydrogenation pathways for C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, C$_{2}$H$_{6}$ on grain surfaces.}
  \label{symbol}
  \begin{tabular}{@{}|cc}
  \hline
Reaction & Branching ratios (\%)\\
   \hline
1 C$_{2}$H$_{4}$ $\rightarrow$ mC$_{2}$H$_{4}$ & 10-1 \\
2 C$_{2}$H$_{6}$ $\rightarrow$ mC$_{2}$H$_{6}$ & 100 \\
3 C$_{2}$H$_{2}$ $\rightarrow$ mC$_{2}$H$_{2}$ & 10 \\
4 C$_{2}$H$_{2}$ $\rightarrow$ mC$_{2}$H$_{4}$ & 40-1 \\
5 C$_{2}$H$_{2}$ $\rightarrow$ mC$_{2}$H$_{6}$ & 50-89 \\
6 C$_{2}$H$_{4}$ $\rightarrow$ mC$_{2}$H$_{6}$ & 90-99 \\
   \hline
  \end{tabular}
 \end{minipage}
\end{table}

Based on the experimental studies performed to date (Barratta et al. 2006), each time C$_{2}$H$_{6}$ is produced, a small experimental signature that has been assigned to icy ethylene can also be observed; it is therefore appropriate to investigate the minimum amount of ethylene on grains necessary to reproduce the relevant c-C$_{2}$H$_{4}$O observations. In a similar manner as described above for the case of oxygen, we vary the branching ratios for the different channels involved in the relevant hydrogenation pathways (see Table 4) in order to determine the minimum concentration of C$_{2}$H$_{4}$ on grains needed for the formation of ethylene oxide. As shown in Table 4, we investigated the variation of the branching of C$_{2}$H$_{2}$ hydrogenation to mC$_{2}$H$_{4}$ (m denotes a molecule in an icy mantle) between 1$\%$ and 40$\%$ (Table 4) and we obtain mC$_{2}$H$_{4}$ fractional abundances (relative to the total hydrogen nuclei) of 9.31x10$^{-11}$ and 2.53x10$^{-9}$ respectively. We find the lower abundance of mC$_{2}$H$_{4}$ still generates viable quantities of c-C$_{2}$H$_{4}$O and so we fixed the branching of reaction 4 to 1$\%$ in all further studies. We then investigated the effect of the role of accretion of gas-phase C$_{2}$H$_{4}$ to form mC$_{2}$H$_{4}$, a process which competes with hydrogenation of C$_{2}$H$_{4}$ to yield mC$_{2}$H$_{6}$ (reaction 1 and 6, Table 4). Here we find that a branching of only 1$\%$ in the direction of mC$_{2}$H$_{4}$ still results in a c-C$_{2}$H$_{4}$O abundance in agreement with observations.

In summary, the above investigations of the hydrogenation pathways showed that a branching of only 1$\%$ towards both reactions 1 and 4 in Table 4 yields c-C$_{2}$H$_{4}$O abundances (3.9x10$^{-12}$ cm$^{-2}$) comparable with observations and these branching ratios were used in all further modelling. Finally, we qualitatively model a sample of high-mass star-forming regions by looking at the influence of the physical parameters of the sources on the ethylene oxide abundances. Models are discussed in Section 3.

\section{Astronomical Implications}

Based on the study by Ikeda et al. (2001) ethylene oxide has been detected towards the Orion compact ridge and several hot cores. We qualitatively model these regions with UCL\_CHEM by varying the final sizes and densities of Phase I and the temperatures in Phase II. In Phase I the non-thermal desorption is included as in Roberts et al. (2007). During Phase II icy mantles sublimate as in Viti et al. (2004).
In order to investigate the sensitivity of the chemistry to the physical parameters of dense cores we run a grid of 16 chemical models by varying: (i) the size of the core, (ii) the final density of the collapsing core, (iii) the temperature, (iv) the adsorption energy (controlled by $\epsilon$) and (v) the percentage of the accreted species onto the grain surfaces.
Table 5 lists the different models and the c-C$_{2}$H$_{4}$O fractional abundances ($f$) (relative to hydrogen nuclei) obtained as outputs (see the last column). We also derive the column density (N) by using the formula below:

\begin{equation}
N(c-C_{2}H_{4}O) = X \times A_{\rm v} \times N(H_{2}),
\end{equation}

where X is the ethylene oxide fractional abundance, $A_{\rm v}$ is the visual extinction and N(H$_{2}$) is equal to 1.6x10$^{21}$ cm$^{-2}$ which is the hydrogen column density in 1 magnitude.

\begin{table*}
 \begin{minipage}{126mm}
  \caption{List of chemical models and their parameters: size, density (n$_{H}$), gas temperature (T) during Phase II, threshold adsorption energy ($\epsilon$), efficiency of the freeze-out (fr) during phase I, the percentage of mantle CO (mCO) given by the freeze-out parameter at the end of Phase I of the chemical model. The last column lists the c-C$_{2}$H$_{4}$O theoretical fractional abundances ($f$) (relative to the total number of hydrogen nuclei) obtained as outputs from each model. We also report its calculated column densities (N).}
  \label{symbol}
  \begin{tabular}{@{}|ccccccccc}
  \hline
 & Size (pc) & n$_{H}$/cm$^{-3}$ & T/K & $\epsilon$ & fr & mCO (\%) & N/cm$^{-2}$ & $f$ \\
   \hline
M1 & 0.02 & 2x10$^{8}$ & 300 & 0.1 & 0.3 & 98 & 1.1x10$^{13}$ & 1.8x10$^{-11}$  \\
M2 & 0.02 & 2x10$^{8}$ & 300 & 0.01 & 0.3 & 98 & 6.2x10$^{12}$ & 1.0x10$^{-11}$ \\
M3 & 0.02 & 2x10$^{8}$ & 300 & 0.1 & 0.0075 & 75 & 1.1x10$^{12}$ & 1.8x10$^{-12}$ \\
M4 & 0.02 & 2x10$^{8}$ & 300 & 0.01 & 0.0075 & 75 & 1.1x10$^{12}$ & 1.7x10$^{-12}$\\
M5 & 0.02 & 2x10$^{7}$ & 200 & 0.1 & 0.3 & 98 & 6.8x10$^{12}$ & 1.1x10$^{-11}$\\
M6 & 0.02 & 2x10$^{7}$ & 200 & 0.01 & 0.3 & 98 & 3.9x10$^{12}$ & 6.4x10$^{-12}$\\
M7 & 0.02 & 2x10$^{7}$ & 200 & 0.1 & 0.0075 & 75 & 6.2x10$^{13}$ & 1.0x10$^{-10}$\\
M8 & 0.02 & 2x10$^{7}$ & 200 & 0.01 & 0.0075 & 75 & 6.1x10$^{13}$ & 9.9x10$^{-11}$\\
M9 & 0.04 & 2x10$^{8}$ & 300 & 0.1 & 0.3 & 98 & 2.2x10$^{13}$ & 1.8x10$^{-11}$\\
M10 & 0.04 & 2x10$^{7}$ & 200 & 0.1 & 0.3 & 98 & 1.5x10$^{13}$ & 1.2x10$^{-11}$\\
M11 & 0.06 & 2x10$^{6}$ & 80 & 0.1 & 0.5 & 98 & 2.7x10$^{14}$ & 1.5x10$^{-10}$\\
M12 & 0.06 & 2x10$^{6}$ & 80 & 0.01 & 0.5 & 98 & 1.5x10$^{14}$ & 8.1x10$^{-11}$ \\
M13 & 0.06 & 2x10$^{6}$ & 80 & 0.1 & 0.25 & 75 & 1.5x10$^{14}$ & 8.1x10$^{-11}$\\
M14 & 0.06 & 2x10$^{6}$ & 80 & 0.01 & 0.25 & 75 & 1.6x10$^{14}$ & 8.7x10$^{-11}$\\
M15 & 0.06 & 2x10$^{6}$ & 80 & 0.1 & 0.1 & 40 & 6.8x10$^{14}$ & 3.7x10$^{-10}$\\
M16 & 0.06 & 2x10$^{6}$ & 80 & 0.01 & 0.1 & 40 & 6.4x10$^{14}$ & 3.5x10$^{-10}$\\
   \hline
  \end{tabular}
 \end{minipage}
\end{table*}

\subsection{Sensitivity to variation in the density, temperature and size core}

Models from M1 to M10 reproduce the typical physical parameters for an hot core model as described in Lerate et al. (2008); models M11-M16 best represent the characteristics for a compact ridge where we set a temperature of 80 K in order to look at the molecular abundances when the water desorption has not yet occurred. We consider a star of 25 M$\odot$ and 0.02-0.06 pc in size, characterized by a visual extinction of about 300-500 mags.
We then vary the density of the core (n$_{H}$) in the 10$^{6}$-10$^{8}$ cm$^{-3}$ range and we also modify the temperature (T) according to the observations towards the different astronomical sources, even if for T $\ge$ 100 K we do not expect any influence of this parameter since, as already explained, the sublimation of ethylene oxide occurs simultaneously with the water desorption (at around 100 K).

If we compare, for instance, M1, M5 and M6, no notable differences in the column densities can be found: with increasing density the values are slightly higher as expected from the faster chemistry in a denser medium.
On the other hand, keeping all the other parameters unchanged, the core size may influence the abundances of c-C$_{2}$H$_{4}$O: the final c-C$_{2}$H$_{4}$O column density calculated from the output of M9 is double the value estimated for M1; the differences in the ethylene oxide abundances can vary up to one order of magnitude (see M1 and M11). This observation is explained by the fact that the larger the core the smaller its visual extinction (A$_{\rm v}$) at a certain density. The magnitude of A$_{\rm v}$ is related to the ability of the UV flux penetrating into the core and the UV flux leads in turn to the photodissociation of species; in other words, in the larger regions we will observe a greater efficiency of molecular photodissociation.

\subsection{Sensitivity to variation in the percentage of the accreted species onto the grain surfaces}

Based on the study by Roberts et al. (2007) we introduce three non-thermal desorption mechanisms during the cold phase of our model (Phase I). Specifically, desorption resulting from: H$_{2}$ formation on grains, direct cosmic ray heating and cosmic ray-induced photodesorption. In particular, the surface formation of molecular hydrogen releases an energy of around 4.5 eV capable of inducing the thermal desorption of molecules from the grains. This is a selective process which involves mainly the most volatile species. The cosmic-ray induced photodesorption is also a selective process that accounts for the number of molecules released to the gas-phase after each cosmic ray impact. These cosmic rays are able to ionize and excite the surrounding gas resulting in the emission of a UV photon, even where the visual extinction is high enough to extinguish any UV radiation propagating in from the outside of the cloud. This internal UV flux impinges on the grain surfaces and dissociates the absorbed molecules (mainly H$_{2}$O). As for the case of H$_{2}$ formation, molecules desorb due to energy released by recombination of OH and O; this type of desorption is non-selective. All the three non-thermal desorption processes are found to be very significant in dark cloud conditions where temperatures are low enough for the thermal desorption to be insignificant. Finally, direct cosmic rays are an efficient heating source able to penetrate deeply into molecular clouds; indeed when these energetic particles impact on the grains, molecules (adsorbed on the surfaces) start to desorb because of the energy transfer from the cosmic rays to the icy mantles.

The efficiency of the non-thermal desorption of molecules depends on the percentage of the accreted species on the grain surfaces, which is in turn linked to the rate of freeze-out. The latter quantity is a free parameter in our models controlled by the product of the sticking coefficient and the average grain surface area. We vary the efficiency of freeze out (fr) in order to obtain different percentages of solid CO within a range of 75-98 \% for the hot core models and 40-98 \% for the compact ridge model. We finally investigate the effects of varying the adsorption energy that is controlled by the $\epsilon$ parameter. The latter is defined as the number of species removed each time an H$_{2}$ molecule is formed. The value of $\epsilon$ is uncertain and may be greater than 1, but, as stated by Roberts et al. (2007), at large values, the desorption will be so efficient as to prevent the build-up of any mantle material until very late times. This justifies the choice of small values for $\epsilon$, 0.01 and 0.1, in order to look at the effect of this free parameter on the desorption rates. As reported in Table 5, the column densities are slightly higher when the H$_{2}$ formation becomes the dominant non-thermal desorption mechanism ($\epsilon$ = 0.1).

\vspace{5mm}
In general, we find that all the changes made to the physical parameters (see Table 5) produce fluctuations of the estimated ethylene oxide column densities within two orders of magnitude; indeed our theoretical column densities are in the 10$^{12}$-10$^{14}$ cm$^{-2}$ range for all models. This observation may seem surprising for the case of the compact ridge models since at 80 K (the temperature adopted for these models) the co-desorption of molecules with water has not yet occurred. This explanation for the constancy of the column density of ethylene oxide may reflect the fact that non-thermal and thermal desorption provide the same contribution to the molecular abundances, but under different physical conditions. Specifically, the chemistry of denser, but smaller regions is mostly influenced by high temperature effects. However, in low density, but more extended sources, the chemistry is driven mainly by non-thermal desorption phenomena.

\section{Results and Discussion}

\begin{figure}
 \includegraphics[width=60mm, angle=270]{F2.eps}
 \caption{Fractional abundances of ethylene oxide in both gas-phase and grain surfaces as a function of time during Phase II.}
\end{figure}

\begin{table*}
 \begin{minipage}{126mm}
  \caption{Comparison between our theoretical column densities (in cm$^{-2}$) with those observed towards several high-mass star-forming regions.}
  \label{symbol}
  \begin{tabular}{@{}|ccccc}
  \hline
 & c-C$_{2}$H$_{4}$O & CH$_{3}$OH & HCOOCH$_{3}$ & CH$_{3}$CN \\
   \hline
Model & 6.2x10$^{13}$ & 2.4x10$^{17}$ & 5.6x10$^{15}$ & 2.2x10$^{16}$\\
Observations & 1.3x10$^{13}$-1.7x10$^{14}$$^{a}$ & 1.4x10$^{16}$-1.4x10$^{17}$$^{a}$ & 1.0x10$^{14}$-1.2x10$^{16}$$^{a}$ & $\ge$ 2.1x10$^{16}$$^{b}$\\
   \hline
  \end{tabular}
\footnotetext[1]{Molecular column density taken from Ikeda et al. (2001)}
\footnotetext[2]{Molecular column density taken from Wilner et al. (1994)}
 \end{minipage}
\end{table*}

Figure 2 reports the evolution of the fractional abundances of ethylene oxide in both gas-phase (solid line) and solid state (dashed line) during Phase II for Model M7 (see Table 5). At the end of the collapse, c-C$_{2}$H$_{4}$O is mostly on the grain surfaces, but its abundance drops rapidly at around 100 K when its co-desorption with water occurs (after 10$^{4}$ yrs). At the same time we observe a small growth in the gas-phase c-C$_{2}$H$_{4}$O abundances thanks to the thermal desorption of solid C$_{2}$H$_{4}$O. The amount of ethylene oxide in the gas-phase remains high at all times due to non-thermal desorption effects.

\begin{figure}
 \includegraphics[width=60mm, angle=270]{F3.eps}
 \caption{Theoretical fractional abundances of three typical tracers of hot cores: CH$_{3}$OH (solid line), HCOOCH$_{3}$ (dashed line) and CH$_{3}$CN (dotted line).}
\end{figure}

We now compare the theoretical column densities with the observations towards several high-mass star-forming regions where c-C$_{2}$H$_{4}$O has been detected. In general, we achieve a good agreement, within an order of magnitude, with all of the observational values for c-C$_{2}$H$_{4}$O reported in Table 3 by Ikeda et al. (2001). Referring to the range of column densities reported by these authors, we consider M7 as the model representing most of the physical characteristics of the regions listed in Ikeda et al. (2001). Results are reported in Table 6 where the calculated and observed column densities for three typical tracers for hot cores: CH$_{3}$OH, HCOOCH$_{3}$ and CH$_{3}$CN are listed. The first two species were detected by Ikeda et al. (2001), the latter molecule was investigated by Wilner et al (1994) which mapped the methyl cyanide lines in the millimeter spectrum of Orion KL region. The trends of the fractional abundances of these species as a function of time are displayed in Figure 3. The c-C$_{2}$H$_{4}$O, CH$_{3}$OH (solid line), HCOOCH$_{3}$ (dashed line) and CH$_{3}$CN (dotted line) theoretical abundances perfectly match their estimated column densities.

This result suggests that the reaction between atomic oxygen and ethylene on the grains is a suitable formation route for ethylene oxide.

\section{Conclusions}
We have investigated the feasibility of the reaction between oxygen ice and ethylene ice as possible route for the formation of ethylene oxide.
For this purpose, we employed a chemical model for a collapsing core by updating the rate for the formation of ethylene oxide by using the recent experimental determination of this quantity of Ward \& Price (2011). We also tested the sensitivity of ethylene oxide fractional abundances to the changes in the main physical parameters of the collapsing core. We finally reproduced the effects of non-thermal and thermal desorption of ethylene oxide during both the collapse and the warm-up phase respectively.
Our work shows a good agreement between our estimated column densities for ethylene oxide (under hot core conditions) and those observed towards several high-mass star-forming regions (Ikeda et al. 2001). In order to verify the reliability of our comparison we also looked at the theoretical and observational abundances of other molecules, such as CH$_{3}$OH, HCOOCH$_{3}$ and CH$_{3}$CN, known as good tracers of hot cores. We conclude that c-C$_{2}$H$_{4}$O can form on the grain surfaces thanks to the reaction between oxygen and ethylene.
We stressed that only 1\% of oxygen atoms depletion has to be in form of oxygen ice in order to produce ethylene oxide on the icy mantle (see Section 2); we also note that there is experimental evidence for a small quantity of ethylene being produced during the hydrogenation of acetylene to form ethane; therefore it is robust to assume that 1$\%$ of acetylene is hydrogenated to solid ethylene as an initial condition of our model. This amount of ethylene is already enough for the observed c-C$_{2}$H$_{4}$O column densities to be reproduced by our theoretical results. Further experimental studies are perhaps needed to extend our study to other cyclic species that have been detected in the ISM, although the complexity of larger molecules may prevent their isomeric speciation as is the case for the recent laboratory study of the formation of C$_{3}$H$_{6}$O.

\section{Acknowledgments}
The research leading to these results has received funding from the [European Community's] Seventh Framework Program [FP7/2007-2013] under grant agreement n$^{\circ}$ 238258.

The authors thank the anonymous referee for useful suggestions that greatly improved the manuscript.