1. science
2. physics

Initial Conditions for Imperfect Dark Matter
arXiv:1507.00291v1 [gr-qc] 30 Jun 2015
Sabir Ramazanov∗
Service de Physique Th´eorique, Universit´e Libre de Bruxelles (ULB),
CP225 Boulevard du Triomphe, B-1050 Bruxelles, Belgium
July 2, 2015
Abstract
We discuss initial conditions for the recently proposed Imperfect Dark Matter (Modified Dust). We show that they are adiabatic under fairly moderate assumptions about
the cosmological evolution of the Universe at the relevant times.
1
Introduction and Summary
Planck data favors adiabatic initial conditions at the onset of the hot Big Bang [1]. This
corresponds to the situation in the early Universe, when the number densities of different
particle species, e.g., dark matter (DM) particles and photons, are universally distributed in
space. That picture, comfortably accomodated in the ΛCDM cosmology supplemented by
the short stage of an inflationary expansion, can be less transparent in more exotic setups. In
particular, adiabaticity of initial conditions is obscured in the case, if DM has a non-particle
origin.
In the present paper, we continue to discuss Imperfect Dark Matter (IDM) scenario [2, 3]
(Modified Dust in Ref. [2]). The action of the model is given by
Z
γ(ϕ)
λ µν
2
4 √
(1)
(g ∂µ ϕ∂ν ϕ − 1) +
(ϕ) ,
S = d x −g
2
2
(see Ref. [4] for generalizations). Here λ and ϕ are two scalars, and γ(ϕ) is some function
of the field ϕ. The Lagrange multiplier λ enforces the constraint g µν ∂µ ϕ∂ν ϕ = 1. This
equation defines the geodesics followed by the collisionless particles, with the field ϕ being
the velocity potential. The analogy with the pressureless perfect fluid, dust, is complete in
∗
e-mail: Sabir.Ramazanov@ulb.ac.be
1
Cyclic universe from Loop Quantum Gravity
Francesco Cianfrani† , Jerzy Kowalski-Glikman†, Giacomo Rosati†
†
Institute for Theoretical Physics, University of Wroclaw,
Pl. Maksa Borna 9, Pl–50-204 Wroclaw, Poland.
(Dated: July 2, 2015)
We discuss how a cyclic model for the flat universe can be constructively derived from Loop Quantum
Gravity. This model has a lower bounce, at small values of the scale factor, which shares many
similarities with that of Loop Quantum Cosmology. We find that quantum gravity corrections can
be also relevant at energy densities much smaller than the Planckian one and that they can induce
an upper bounce at large values of the scale factor.
arXiv:1507.00226v1 [gr-qc] 1 Jul 2015
PACS numbers: 04.60.Pp
The great expectation we have of a quantum theory
of gravity is that it will be able to tame the singularities arising in General Relativity. In particular, it should
provide us with a consistent, non-singular scenario for
the initial stages of our universe. It is hoped that this
goal can be achieved in the framework of a quantum cosmological model, constructed by assuming homogeneity
and isotropy in the full Quantum Gravity theory, and
thus reducing it effectively to a finite dimensional model.
Unfortunately, there is no unique way to perform such
reduction to quantum cosmology, and the final outcome
depends on which path is chosen.
Loop Quantum Gravity (LQG) [1–3] is a framework
for a non-perturbative canonical quantization of GR.
Loop Quantum Cosmology (LQC) [4, 5] attempts to
quantize the cosmological model by applying the techniques of LQG to [a] low-dimensional system obtained
by symmetry-reduction of the phase-space of classical GR
(minisuperspace). Since the dynamics in minisuperspace
is much simpler than in the full theory, this procedure
makes it possible to solve some issues of LQG and to realize a consistent quantum scenario for FRW models and
for its most relevant anisotropic extensions (Bianchi type
I [6, 7], II [8] and IX [9]). The most important result of
these investigations is that the initial singularity is replaced by a bounce occurring when the energy density
reaches the critical value [10, 11].
A different path to quantum cosmology is proposed
in Quantum Reduced Loop Gravity (QRLG) [12–14], in
which one infers the cosmological sector of LQG directly
on a quantum level. This procedure consists of several
steps. First, one imposes quantum gauge-fixing conditions by applying the Gupta-Bleuler method, reducing
the number of variables and simplifying the expression
of geometric operators. Then, one analyzes the quantum Hamiltonian neglecting the non-local terms which
arise because of the gauge fixings. It is at this stage that
one imposes the homogeneity requirement, since the nonlocal contributions to the Hamiltonian can be avoided
only for homogeneous configurations (e.g. Bianchi I
model). Finally, the expectation value of the quantum
Hamiltonian on semiclassical states is computed [15, 16]
and it is taken as the generator of the effective semiclassical dynamics.
The models of quantum universe arising from LQC
and QRLG differ from the classical one by the presence
of two kind of corrections, holonomy and inverse-volume
ones, the former being responsible for the bounce replacing the initial singularity. The magnitude of the corrections is determined by a regulator. In LQC the regulator is treated as a free function of phase-space variables,
usually fixed by the demand that the bounce occurs at
Planckian energy density [17] (the so-called “improved
regularization scheme”). This choice is motivated by the
requirement that quantum gravity effects cannot be relevant at length scales much larger than the Planckian
ones. In QRLG, the regulator equals the third root of
the inverse number of nodes of the fundamental graph
at which the states are based. The most natural choice
is to take the regulator to be a constant, since the semiclassical analysis is performed with a non-graph-changing
Hamiltonian, which does not add or remove nodes. However, in principle one could look for different definitions
of the Hamiltonian and the semiclassical states, which
could lead to non-constant regulators and, as a result, to
different regularizations.
Within the LQC research program, the choice of a
constant regulator, also known as the µ0 regularization
scheme, was adopted in the original proposal (see for instance [19]), but it has been abandoned, in favor of the
improved regularization scheme, due to possible unusual
phenomenological implications [20, 21]. However a systematic analysis of the phenomenological downsides was
still missing. In particular, the main reason for discarding
a constant regulator is that it implies the energy density
at which the bounce occurs to be proportional to the inverse square of the scale factor, so that the quantum gravity effects can be in principle relevant at energy densities
much smaller than Planckian one. This is due to the fact
that in LQG a continuous geometry is described by a discrete graph, having a fixed number of nodes, with SU (2)
quantum numbers at each link. The physical length scale
is determined by both the spin numbers and the fiducial
discrete geometry associated with the graph, whose lattice size equals the third root of the inverse number of
nodes. As the universe grows, spin numbers grow, the
fiducial lattice size stays fixed but the physical lattice
size of the geometry keeps growing.
IPMU-15-0099
July 2, 2015
arXiv:1507.00119v1 [hep-ph] 1 Jul 2015
Dynamics of Peccei-Quinn Breaking Field after Inflation
and
Axion Isocurvature Perturbations
Keisuke Harigaya(a,b) , Masahiro Ibe(a,b) ,
Masahiro Kawasaki(a,b) and Tsutomu T. Yanagida(b)
(a) ICRR,
(b) Kavli
University of Tokyo, Kashiwa, Chiba 277-8582, Japan
IPMU (WPI), UTIAS, University of Tokyo, Kashiwa, Chiba 277-8583, Japan
The Peccei-Quinn mechanism suffers from the problem of the isocurvature perturbations. The isocurvature perturbations are suppressed if the
Peccei-Quinn breaking scale is large during inflation. The oscillation of the
Peccei-Quinn breaking field after inflation, however, leads to the formation
of domain walls due to the parametric resonance effect. In this paper, we
discuss the evolution of the Peccei-Quinn breaking field after inflation in detail, and propose a model where the parametric resonance is ineffective and
hence domain walls are not formed. We also discuss consistency of our model
with supersymmetric theory.
Modified Eddington-inspired-Born-Infeld Gravity with a Trace Term
Che-Yu Chen
1,2,3 a
, Mariam Bouhmadi-L´opez4,5,6,7,b and Pisin Chen
1,2,3,8c
1
arXiv:1507.00028v1 [gr-qc] 30 Jun 2015
Department of Physics, National Taiwan University, Taipei, Taiwan 10617
2
LeCosPA, National Taiwan University, Taipei, Taiwan 10617
3
Graduate Institute of Astrophysics, National Taiwan University, Taipei, Taiwan 10617
4
Departamento de F´ısica, Universidade da Beira Interior, 6200 Covilh˜
a, Portugal
5
Centro de Matem´
atica e Aplica¸c˜
oes da Universidade da Beira Interior (CMA-UBI), 6200 Covilh˜
a, Portugal
6
Department of Theoretical Physics, University of the Basque Country UPV/EHU, P.O. Box 644, 48080 Bilbao, Spain
7
IKERBASQUE, Basque Foundation for Science, 48011, Bilbao, Spain
8
Kavli Institute for Particle Astrophysics and Cosmology,
SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, U.S.A.
(Dated: July 2, 2015)
In this paper, a modified Eddington-inspired-Born-Infeld (EiBI) theory with a pure trace term
gµν R being added to the determinantal action is analysed from a cosmological point of view. It
corresponds to the most general action constructed from a rank two tensor that contains up to first
order terms in curvature. This term can equally be seen as a conformal factor multiplying the metric
gµν . This very interesting type of amendment has not been considered within the Palatini formalism
despite the large amount of works on the Born-Infeld-inspired theory of gravity. This model can
provide smooth bouncing solutions which were not allowed in the EiBI model for the same EiBI
coupling. Most interestingly, for a radiation filled universe there are some regions of the parameter
space that can naturally lead to a de Sitter inflationary stage without the need of any exotic matter
field. Finally, in this model we discover a new type of cosmic “quasi-sudden” singularity, where the
cosmic time derivative of the Hubble rate becomes very large but finite at a finite cosmic time.
PACS numbers: 98.80.-k, 98.80.Jk, 04.50.Kd, 04.20.Dw
Keywords: modified theories of gravity, early time universe, Eddington-inspired-Born-Infeld theories, cosmic
singularities
I.
INTRODUCTION
Undeniably, Einstein’s theory of general relativity
(GR) has been an extremely successful theory for around
a century [1]. However, the theory is expected to break
down at some points at very high energies where quantum
effects are expected to become crucial, such as in the past
expansion of the Universe where GR predicts a big bang
singularity [2]. This is one of the motivations for looking for possible modified theories of gravity, which are
hoped to not only be able to preserve the huge achievements of GR, but also to shed some light on smoothing
the singularities predicted in GR. Such theories could be
seen as effective/phenomenological approaches of a more
fundamental quantum theory of gravity.
In 1934 Born and Infeld proposed a non-linear action
for classical electrodynamics, which is characterised by
its grand success in solving the divergence of the selfenergy of point-like charges [3]. This action for electromagnetism has an elegant determinantal structure which
reads:
Z
q
1
SBI =
(1.1)
d4 x |gµν + κFµν |,
κ
with gµν and Fµν being the metric tensor and the electromagnetic field strength. Here κ is a dimensional constant.
a
b
c
b97202056@ntu.edu.tw
mbl@ubi.pt. On leave of absence from UPV and IKERBASQUE.
pisinchen@phys.ntu.edu.tw
Note that the Born-Infeld action (1.1) recovers Maxwell
action for small amplitudes.
Since the proposal of the Born-Infeld action, modified
theories of gravity with a Born-Infeld-inspired action, initiated by the pioneering work of Deser and Gibbons [4],
have received much attention; c.f. Refs. [5–10]. These
theories of gravity not only maintain the properties of
GR for small curvatures, but provide various interesting
deviations from GR at high curvature regimes. Most of
these works start with a general gravitational action of
the form
Z
q
1
d4 x |gµν + κGµν (R, Rαβ , Rµναβ )|, (1.2)
Sdet =
κ
where Gµν represents a linear Ricci term Rµν , plus higher
order curvature terms containing R, Rµν , and Rµναβ .
Rather than amending the higher curvature terms in Gµν ,
we can also modify the action by multiplying the metric
term in the determinant by a factor 1 + α1 R, where α1 is
a constant with length square dimensions. This approach
does not lead to a violation of the requirements in the low
curvature limit, i.e. we can still recover GR at low energies. Actually, a gravitational action with this pure trace
term has been considered in Ref. [5] within a pure metric
formalism. The theory thus inevitably suffers from the
presence of troublesome fourth order derivatives in the
field equations or ghost instabilities [4].
In order to keep the theory free from aforementioned
problems, alternative theories formulated within the
Palatini formalism and teleparallel representation have
arXiv:1506.09057v1 [nucl-th] 30 Jun 2015
European Physical Journal Web of Conferences will be set by the publisher
DOI: will be set by the publisher
c Owned by the authors, published by EDP Sciences, 2015
Symmetry energy systematics and its high density behavior
Lie-Wen Chen1,a
1
Department of Physics and Astronomy and Shanghai Key Laboratory for Particle Physics and Cosmology,
Shanghai Jiao Tong University, Shanghai 200240, China
Abstract. We explore the systematics of the density dependence of nuclear
matter symmetry energy in the ambit of microscopic calculations with various
energy density functionals, and find that the symmetry energy from subsaturation density to supra-saturation density can be well determined by three
characteristic parameters of the symmetry energy at saturation density ρ0 , i.e.,
the magnitude Esym (ρ0 ), the density slope L and the density curvature Ksym .
This finding opens a new window to constrain the supra-saturation density behavior of the symmetry energy from its (sub-)saturation density behavior. In
particular, we obtain L = 46.7 ± 12.8 MeV and Ksym = −166.9 ± 168.3 MeV as
well as Esym (2ρ0 ) ≈ 40.2±12.8 MeV and L(2ρ0 ) ≈ 8.9±108.7 MeV based on the
present knowledge of Esym (ρ0 ) = 32.5 ± 0.5 MeV, Esym (ρc ) = 26.65 ± 0.2 MeV
and L(ρc ) = 46.0 ± 4.5 MeV at ρc = 0.11 fm−3 extracted from nuclear mass and
the neutron skin thickness of Sn isotopes. Our results indicate that the symmetry energy cannot be stiffer than a linear density dependence. In addition,
we also discuss the quark matter symmetry energy since the deconfined quarks
could be the right degree of freedom in dense matter at high baryon densities.
1 Introduction
The nuclear matter symmetry energy, which essentially characterizes the isospin dependent
part of the equation of state (EOS) of asymmetric nuclear matter, is important for understanding many questions in nuclear physics and astrophysics, including the nuclear effective
interactions in asymmetric nuclear matter, the structure and stability of exotic nuclei, the reaction dynamics induced by rare isotopes, the nature and evolution of neutron stars, and the
mechanism of supernova explosion [1–7]. The symmetry energy also plays an important role
in some interesting issues of new physics beyond the standard model [8–12]. During the last
decade, a lot of experimental, observational and theoretical efforts have been devoted to constraining the density dependence of the symmetry energy [13–18]. While significant progress
has been made in determining the density behavior of the symmetry energy around saturation density ρ0 (∼ 0.16 fm−3 ), its supra-saturation density behavior is still poorly known and
remains the most uncertain property of isospin asymmetric nuclear matter. Theoretically,
many experimental and observational probes have been proposed to extract information on
the supra-saturation density behavior of the symmetry energy [13]. In terrestrial laboratories,
a e-mail: lwchen@sjtu.edu.cn
arXiv:1506.08573v1 [hep-th] 29 Jun 2015
First law of thermodynamics for dynamical
apparent horizons and the entropy of Friedmann
universes
Stefano Viaggiu,
Dipartimento di Matematica, Universit`a di Roma “Tor Vergata”,
Via della Ricerca Scientifica, 1, I-00133 Roma, Italy.
E-mail: viaggiu@axp.mat.uniroma2.it
June 30, 2015
Abstract
Recently, we have generalized the Bekenstein-Hawking entropy formula for black holes embedded in expanding Friedmann universes. In
this letter, we begin the study of this new formula to obtain the first
law of thermodynamics for dynamical apparent horizons. In this regard
we obtain a generalized expression for the internal energy U together
with a distinction between the dynamical temperature TD of apparent
horizons and the related one due to thermodynamics formulas. Remarkable, when the expression for U is applied to the apparent horizon of the universe, we found that this internal energy is a constant
of motion. Our calculations thus show that the total energy of our
spatially flat universe including the gravitational contribution, when
calculated at the apparent horizon, is an universal constant that can
be set to zero from simple dimensional considerations. This strongly
support the holographic principle.
Keywords: Bekenstein-Hawking entropy; black holes; expanding universes;
internal energy.
1
Introduction
Since of the Hawking’s discovery that black holes radiate [1], his entropy formula represented a cornerstone for a sound formulation of quantum gravity.
1
Mon. Not. R. Astron. Soc. 000, 000-000 (0000)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00314v1 [astro-ph.HE] 26 Jun 2015
Non-axisymmetric magnetic modes of neutron stars with
purely poloidal magnetic fields
Hidetaka Asai, Umin Lee, and Shijun Yoshida
Astronomical Institute, Tohoku University, Sendai 980-8578, Japan
Received / Accepted
ABSTRACT
We calculate non-axisymmetric oscillations of neutron stars magnetized by purely
poloidal magnetic fields. We use polytropes of index n = 1 and 1.5 as a background
model, where we ignore the equilibrium deformation due to the magnetic field. Since
separation of variables is not possible for the oscillation of magnetized stars, we employ finite series expansions for the perturbations using spherical harmonic functions.
Solving the oscillation equations as the boundary and eigenvalue problem, we find two
kinds of discrete magnetic modes, that is, stable (oscillatory) magnetic modes and unstable (monotonically growing) magnetic modes. For isentropic models, the frequency
or the growth rate of the magnetic modes is exactly proportional to BS , the strength
of the field at the surface. The oscillation frequency and the growth rate are affected
by the buoyant force in the interior, and the stable stratification tends to stabilize the
unstable magnetic modes.
Key words:
1
– stars: magnetic fields – stars: neutron – stars: oscillations.
INTRODUCTION
Neutron stars are believed to have strong magnetic fields. The field strength BS at the surface is estimated to be ∼ 1012 G for
pulsars and ∼ 1015 G for magnetar candidates, but we do not have any good knowledge of the configuration of magnetic fields
inside the star. Recently, quasi-periodic oscillations (QPOs) in the tail of the giant X/γ-ray flare were detected from the soft γ
ray repeaters (SGRs) (e.g., Israel et al. 2005; Strohmayer & Watts 2005). SGRs belong to what we call magnetars. Giant flares
have so far been observed only from three SGRs, that is, SGR 0526-66, 1900+14, and 1806-20, and just once for each of the
SGRs, indicating that giant flares in magnetars are quite rare events. The QPOs are now regarded as a manifestation of global
oscillations of the underlying neutron stars, and it is expected that they can be used for seismological studies of the magnetars.
However, it is not an easy task to properly carry out modal analysis of magnetar candidates because of their extremely strong
magnetic fields, which can significantly modify the modal property of the stars. Because separation of variables between the
radial and angular coordinates is not possible for the perturbations in magnetized stars, we have to employ series expansions
to represent the perturbations in linear analysis, which possibly makes the analysis difficult, particularly when singular points
inherent to the governing equations appear in the interior of the star (see, e.g., Asai, Lee, & Yoshida 2015).
Many authors have calculated axisymmetric (m = 0) oscillation modes of magnetized stars having a purely poloidal
magnetic field, where m is the azimuthal wave number of modes. Using a toy model, Glampedakis, Samuelsson, and Andersson
(2006) and Levin (2006, 2007) investigated global oscillation modes residing in the fluid core and in the solid crust both
threaded by a magnetic field, and they showed that frequency resonance between modes in the core and in the crust could be
important. Lee (2008) and Asai & Lee (2014) carried out normal mode calculations of axisymmetric toroidal modes and found
discrete modes. van Hoven & Levin (2011, 2012) suggested the existence of discrete modes in the gaps between frequency
continua, using spectral method. Besides the normal mode analyses mentioned above, however, most authors have used MHD
simulations to investigate the modal properties of magnetized stars (e.g., Cerd´
a-Dur´
an et al. 2009; Colaiuda & Kokkotas 2011;
Gabler et al. 2011, 2012; Sotani et al. 2008). Sotani & Kokkotas (2009) calculated axisymmetric polar-Alfv´en oscillations and
found that continuous frequency spectra are not formed. Colaiuda & Kokkotas (2012) calculated mixed polar and axial
oscillation modes of magnetized star having both poloidal and toroidal magnetic fields for axisymmetric perturbations. They
suggested that the oscillation spectra can be significantly modified by the toroidal magnetic field component. Gabler et al.
(2013a) calculated magnetic oscillation modes for various magnetic field configurations (e.g., purely poloidal, purely toroidal,
mixed poloidal and toroidal).
For non-axisymmtric (m 6= 0) modes of magnetized stars, however, there are only a few numerical studies of global
oscillations. For purely toroidal magnetic field configuration, Lander et al. (2010), Passamonti & Lander (2013), and Asai,
Mon. Not. R. Astron. Soc. 000, 1–?? (2014)
Printed 2 July 2015
(MN LATEX style file v2.2)
Hu 1-2: a metal-poor bipolar planetary nebula with fast collimated
outflows
arXiv:1507.00284v1 [astro-ph.SR] 1 Jul 2015
X. Fang1? , M. A. Guerrero1 , L. F. Miranda1 , A. Riera2 , P. F. Vel´azquez3 , and A. C. Raga3
1 Instituto
de Astrof´ısica de Andaluc´ıa (IAA, CSIC), Glorieta de la Astronom´ıa s/n, E-18008 Granada, Spain
de F´ısica i Enginyeria Nuclear, EUETIB, Universitat Polit´ecnica de Catalunya, Comte d’Urgell 187, E-08036 Barcelona, Spain
3 Instituto de Ciencias Nucleares, Universidad Nacional Aut´
onoma de M´exico, Apdo. Postal 70-543, CP: 04510, D. F., Mexico
2 Departament
Accepted . Received
ABSTRACT
We present narrow-band optical and near-IR imaging and optical long-slit spectroscopic
observations of Hu 1-2, a Galactic planetary nebula (PN) with a pair of [N II]-bright, fastmoving (>340 km s−1 ) bipolar knots. Intermediate-dispersion spectra are used to derive physical conditions and abundances across the nebula, and high-dispersion spectra to study the
spatio-kinematical structure. Generally Hu 1-2 has high He/H (≈0.14) and N/O ratios (≈0.9),
typical of Type I PNe. On the other hand, its abundances of O, Ne, S, and Ar are low as
compared with the average abundances of Galactic bulge and disc PNe. The position-velocity
maps can be generally described as an hour-glass shaped nebula with bipolar expansion, although the morphology and kinematics of the innermost regions cannot be satisfactorily explained with a simple, tilted equatorial torus. The spatio-kinematical study confines the inclination angle of its major axis to be within 10◦ of the plane of sky. As in the irradiated
bow-shocks of IC 4634 and NGC 7009, there is a clear stratification in the emission peaks of
[O III], Hα, and [N II] in the northwest (NW) knot of Hu 1-2. Fast collimated outflows in PNe
exhibit higher excitation than other low-ionization structures. This is particularly the case for
the bipolar knots of Hu1-2, with He II emission levels above those of collimated outflows
in other Galactic PNe. The excitation of the knots in Hu 1-2 is consistent with the combined
effects of shocks and UV radiation from the central star. The mechanical energy and luminosity of the knots are similar to those observed in the PNe known to harbor a post-common
envelope (post-CE) close binary central star.
Key words: ISM: abundances – ISM: jets and outflows – ISM: planetary nebulae: individual:
Hu 1-2
1
INTRODUCTION
Planetary nebulae (PNe) represent the last stages in the evolution
of low- and intermediate-mass stars, before they turn into white
dwarfs. A large fraction of PNe presents complex morphologies,
including axisymmetric shells, multipolar lobes, and collimated
structures, in sharp contrast with the spherical envelopes typically
seen around asymptotic giant branch (AGB) stars (e.g., Olofsson
et al. 2010). Somehow, the spherical AGB envelope is transformed
into the complex PN shell. Hubble Space Telescope (HST) observations have led to the suggestion that the shaping of the most structured young PNe (or proto-planetary nebulae, proto-PNe) could be
attributed to the action of high-speed collimated outflows or jets
that operate during the late-AGB and/or early post-AGB phase, as
they interact with the intrinsically spherical AGB circumstellar envelope (Sahai & Trauger 1998). These observations indicate that
?
E-mail: fangx@iaa.es
collimated outflows probably also play a crucial role in the shaping
and dynamical evolution of PNe (e.g., Dennis et al. 2008; HuarteEspinosa et al. 2012).
Collimated outflows or jet-like structures have been identified
in a significant number of PNe. The origin of collimated outflows
is still uncertain, although they are most likely related to the evolution of binary systems, the action of magnetic fields, or both (e.g.,
Soker 2006; De Marco 2009). Particularly, binary interactions that
involve accretion (e.g., Soker 1998; Soker & Livio 1994; ReyesRuiz & L´opez 1999; Nordhaus & Blackman 2006; Blackman &
Lucchini 2014) have been considered to be a plausible engine to
produce the ubiquitous collimated outflows in PNe and proto-PNe
(e.g., Balick & Frank 2002). Recently, Tocknell, De Marco & Wardle (2014) constrained the physical properties of the common envelope (CE) interaction using the observed masses and kinematics
of jets in four post-CE PNe.
Detailed studies of the spectra of jet-like structures in PNe
are scarce, but they generally report strong low-excitation emis-
Astronomy & Astrophysics manuscript no. Galli2015_final_printer
July 2, 2015
c
ESO
2015
Evolution of the T Tauri star population in the Lupus association
P. A. B. Galli1 , C. Bertout2 , R. Teixeira1 , and C. Ducourant3
1
2
3
Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão, 1226 - Cidade Universitária,
05508-900, São Paulo - SP, Brazil.
e-mail: galli@astro.iag.usp.br
Institut d’Astrophysique, 98bis, Bd. Arago, 75014 Paris, France.
Observatoire Aquitain des Sciences de l ’Univers, CNRS-UMR 5804, BP 89, Floirac, France.
arXiv:1507.00281v1 [astro-ph.SR] 1 Jul 2015
ABSTRACT
Aims. In a recent study, we derived individual distances for 109 pre-main sequence stars that define the Lupus kinematic association of
young stars. Here, we use these new distances to derive the masses and ages of Lupus T Tauri stars with the aim of better constraining
the lifetime of their circumstellar disks.
Methods. Using the photometric and spectroscopic information available in the literature, we computed the photospheric luminosity
of 92 T Tauri stars in the Lupus association. Then, we estimated their masses and ages from theoretical evolutionary models. Based
on Monte Carlo simulations and statistical tests, we compare the mass and age distribution of the classical T Tauri stars (CTTS) and
weak-line T Tauri (WTTS) in our sample.
Results. We show that the CTTSs are on average younger than the WTTSs and that the probability that both T Tauri subclasses are
drawn from the same mass and age parental distribution is very low. Our results favor the scenario proposed earlier for the TaurusAuriga association, where the CTTSs evolve into WTTSs when their disks are fully accreted by the star. Based on an empirical disk
model, we find that the average disk lifetime for the T Tauri stars in the Lupus association is τd = 3 × 106 (M∗ /M )0.55 yr.
Conclusions. We find evidence that the average lifetime of the circumstellar disks in the Lupus association is shorter than in the
Taurus-Auriga association and discuss the implications of this result.
Key words. stars: formation, stars: pre-main sequence, stars: circumstellar matter
1. Introduction
T Tauri stars (TTSs) are late-type low-mass (M . 2M ) premain sequence stars that are usually associated with molecular
clouds. They tend to be found in young associations (Joy 1945;
Herbig 1962; Bertout 1989). They are classified into one of two
types: (i) the classical T Tauri stars (CTTSs) that show strong
emission lines indicative of accretion from a circumstellar disk,
and (ii) the weak-line T Tauri stars (WTTSs) that exhibit enhanced magnetic activity and show no signs of accretion. To explain the existence of both subgroups, it has been hypothesized
for a long time (Walter et al. 1988; Bertout 1989) that the CTTSs
become WTTSs when their disks dissipate. Although this is not
a problem for the older population (t > 10 Myr) of WTTSs,
this picture implies that the timescale of disk dissipation should
be a few Myr to justify the CTTSs and WTTSs coexisting (see,
e.g., Fig. 15 of Kenyon & Hartmann 1995) in the HertzsprungRussell diagram (HRD). Indeed, previous studies (Haisch et al.
2001; Takagi et al. 2014) indicate that the disk lifetime in different clusters is typically in the range 3-6 Myr. However, Fedele
et al. (2010) demonstrate in a more detailed study that mass accretion and dust dissipation in protoplanetary disks occur on a
typical timescale of 2.3-3.0 Myr (see Fig. 4 of their paper), which
implies a rapid disk evolution.
Some progress in determining the lifetime of circumstellar
disks of TTSs was achieved by Bertout et al. (2007), who used
the kinematic parallaxes derived from a kinematic study of the
Taurus-Auriga T association (Bertout & Genova 2006) and calculated the stellar luminosities, masses, and ages to readdress
the question of the relationship between CTTSs and WTTSs.
As a result, they obtained, for the first time, the lifetime of
a circumstellar disk in terms of the mass of the parent star,
τd = 4 × 106 (M/M )0.75 yr, and demonstrated that the observed
mass and age distributions of the TTSs in the Taurus-Auriga association can be explained by assuming that the CTTSs evolve
into WTTSs when their disks are fully accreted by the star.
While supporting the current scenario of disk evolution for
the Taurus-Auriga T association, this first result needs to be confirmed and compared to other star-forming regions (SFRs). One
reason for extending this study to other SFRs is that the origin of
the initial mass function remains unclear, but is expected to depend on the intrinsic properties of each SFR, such as clustering,
mass segregation, and binarity (Bonnell et al. 2007). Moreover,
the evolution of circumstellar (protoplanetary) disks depends on
environmental criteria (Rosotti et al. 2014), and a comparison
with other SFRs could provide important clues for constraining
the timescale of planet formation in different environments. For
example, Mariñas et al. (2013) conclude that the dependence of
disk fraction on stellar mass in Taurus and the NGC 2264 open
cluster is different (see Fig. 10 of their paper), and they argue
that the lower disk fractions in NGC 2264 might be due to the
much higher stellar density than in Taurus.
In the following, we perform a study similar to the one by
Bertout & Genova (2006) for the TTSs in the Lupus SFR. The
Lupus complex of molecular clouds is located in an environment
surrounded by the more massive stars from the nearby ScorpiusCentaurus OB association (de Zeeuw et al. 1999). This makes
the Lupus SFR an interesting target for comparative studies with
Taurus-Auriga, because of the different environmental and dynamical effects (e.g., solar winds and UV radiation) caused by
Article number, page 1 of 11
c
ESO
2015
Astronomy & Astrophysics manuscript no. Tam_etalFinal
July 2, 2015
Coronal heating in multiple magnetic threads
K. V. Tam1, 2 , A. W. Hood1 , P. K. Browning3 , and P. J. Cargill1, 4
1
2
3
4
arXiv:1507.00259v1 [astro-ph.SR] 1 Jul 2015
5
School of Mathematics and Statistics, University of St. Andrews, St. Andrews, Fife, KY16 9SS, U.K.
Space Science Institute, Macau University of Science and Technology, Avenida Wai Long, Taipa, Macau
Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13
9PL, UK
Space and Atmospheric Physics, The Blackett Laboratory, Imperial College, London SW7 2BW, UK
e-mail: awh@st-andrews.ac.uk
July 2, 2015
ABSTRACT
Context. Heating the solar corona to several million degrees requires the conversion of magnetic energy into thermal
energy. In this paper, we investigate whether an unstable magnetic thread within a coronal loop can destabilise a
neighbouring magnetic thread.
Aims. By running a series of simulations, we aim to understand under what conditions the destabilisation of a single
magnetic thread can also trigger a release of energy in a nearby thread.
Methods. The 3D magnetohydrodynamics code, Lare3d, is used to simulate the temporal evolution of coronal magnetic
fields during a kink instability and the subsequent relaxation process. We assume that a coronal magnetic loop consists
of non-potential magnetic threads that are initially in an equilibrium state.
Results. The non-linear kink instability in one magnetic thread forms a helical current sheet and initiates magnetic
reconnection. The current sheet fragments, and magnetic energy is released throughout that thread. We find that,
under certain conditions, this event can destabilise a nearby thread, which is a necessary requirement for starting an
avalanche of energy release in magnetic threads.
Conclusions. It is possible to initiate an energy release in a nearby, non-potential magnetic thread, because the energy
released from one unstable magnetic thread can trigger energy release in nearby threads, provided that the nearby
structures are close to marginal stability.
Key words. Sun: corona - Sun: magnetic fields - magnetohydrodynamics (MHD) - magnetic reconnection - coronal
heating - avalanche
1. Introduction
Energy release in the solar corona is now known to occur
over a wide range of spatial and temporal scales from large
eruptive flares, through compact flares, to microflares and
the (hypothesised) nanoflares that may be responsible for
heating the quiescent corona. At microflare energies and
above (> 1027 ergs), these events have been observed and
they form a power law distribution E −δ where E is the inferred energy, and δ is in the range 1.4 − 1.8 (Crosby et al.
1993; Hannah et al. 2008). Evidence for nanoflares is only
now becoming clear with modern instrumentation: energies close to 1024 ergs are inferred (Warren et al. 2011;
Testa et al. 2013; Cargill 2014). A task for theoretical models is to account for this range of event sizes, with the added
caveat that larger events should be more ‘difficult’ to trigger
than smaller ones.
The largest solar flares are commonly discussed in terms
of an associated eruption and ejection of plasma from the
Sun (e.g. Priest & Forbes (2000) and Nishida et al. (2009)).
However, the applicability of this scenario to small (compact) flares, microflares and smaller events is questionable,
and other mechanisms need to be identified. One such
fundamental process is the kink instability of a twisted
magnetic field that forms a loop-like structure in the
corona. This is an ‘ideal’ MHD instability and it has fast
(Alfvénic) growth during its initial phase. There is observational evidence of the kink instability (Liu et al. 2007;
Liu & Alexander 2009; Nandy et al. 2008; Srivastava et al.
2010) along with simulations of its development in a single
magnetic loop (Linton et al. 1996; Galsgaard & Nordlund
1997; Baty et al. 1998; Gerrard et al. 2004). A very attractive feature of its non-linear evolution is the rapid formation and dissipation of multiple regions of intense current
(Browning & Van der Linden 2003; Browning et al. 2008;
Hood et al. 2009; Bareford et al. 2011). Such non-linear development of magnetic reconnection can result in a turbulent state (Browning & Lazarian 2013) that is also a good
facilitator of particle acceleration (Gordovskyy et al. 2013).
The field subsequently relaxes to a lower energy state, a
process that can be described by helicity-conserving, Taylor
relaxation (Taylor 1974; Browning & Van der Linden 2003;
Browning et al. 2008; Bareford et al. 2013). The energy released in the instability is of the order of 1028 ergs, which
is characteristic of a microflare or a nanoflare storm. It has
been shown that random footprint driving within a single
loop leads to a temporal sequence of heating events of various sizes, which can take a power law form (Bareford et al.
2010, 2011).
It is likely that the magnetic footpoints of a large coronal loops are not located in a single magnetic source but
Article number, page 1 of 16
Infrared-excess Source DSO/G2 Near the Galactic
Center: Theory vs. Observations
M. Zajaˇcek and A. Eckart
University of Cologne, First Institute of Physics, Cologne, Germany
Max-Planck Institute for Radioastronomy, Bonn, Germany
arXiv:1507.00237v1 [astro-ph.GA] 1 Jul 2015
F. Peissker and G. D. Karssen
University of Cologne, First Institute of Physics, Cologne, Germany
V. Karas
Astronomical Institute of the Academy of Sciences of the Czech Republic
Abstract. We monitored the Dusty S-cluster object (DSO/G2) during its closest
approach to the Galactic Center supermassive black hole in 2014 with ESO
VLT/SINFONI. We report on our findings, i.e. ionized-hydrogen emission from
the DSO that remains spatially compact before and after the peribothron passage.
The detection of DSO/G2 object as a compact single-peak emission line source is
in contradiction with the original hypothesis of a core-less cloud that is necessarily
tidally stretched, hence producing double-peak emission line profile around the
pericentre passage. This strengthens the evidence that the DSO/G2 source is a
dust-enshrouded young star. The accretion of material from the circumstellar disc
onto the stellar surface can contribute significantly to the emission of Brγ line as
well as to the observed large line width of the order of 10 angstroms.
Introduction
The fast-moving infrared excess source Dusty S-cluster object (DSO; also named G2) within
the cluster of high-velocity S-stars was discovered by Gillessen et al. (2012) and subsequently
independently monitored and analyzed (Eckart et al., 2013; Gillessen et al., 2013a,b; Witzel
et al., 2014; Pfuhl et al., 2015; Valencia-S. et al., 2015) to detect a potential effect on the
activity of the supermassive black hole (hereafter SMBH) as well as to determine the character
of the object.
The main debate concerning the nature of the object is whether it is a core-less gas and
dust cloud or a compact stellar source. While a core-less cloud of a few Earth masses (Gillessen
et al., 2012) is necessarily tidally stretched and strongly interacts with the surrounding medium,
a star keeps its compact nature and continues to orbit along the original trajectory. The orbital
elements of the DSO are currently constrained well, see Valencia-S. et al. (2015) for a recent fit.
The source has passed the pericentre in 2014 at the distance of rp = a(1 − e) = 0.033 pc × (1 −
0.976) ≈ 2000 rs , where rs is the Schwarzschild radius (rs ≡ 2GM• /c2 ≈ 2.95 × 105 M• /M cm;
we consider M• = 4 × 106 M for the SMBH mass associated with the compact radio source
Sgr A*). Any stellar object associated with the DSO is not tidally disrupted, since the critical
tidal radius for the SMBH is rt = R? (3M• /M? )1/3 ≈ 13(R? /R )(M? /M )−1/3 rs , where R? is
the stellar radius. However, the tidal radius for an extended, unbound cloud with the radius of
Rcl ≈ 30 mas ≈ 5 × 104 R and the mass of three Earth masses (∼ 10−5 M ) (Gillessen et al.,
2012) is ∼ 3 × 107 rs , hence it should undergo the full tidal disruption at the pericentre. Thus,
there is a clear distinction between the star and the cloud concerning the orbital stability. The
predictions of evolution for different scenarios and geometries are analysed in Zajaˇcek, Karas &
Eckart (2014) and references therein.
In this contribution we further elaborate on the model of a dust-enshrouded star (MurrayClay & Loeb, 2012; Scoville & Burkert, 2013; Ballone et al., 2013; Zajaˇcek, Karas & Eckart,
c ESO 2015
Astronomy & Astrophysics manuscript no. Deharveng˙astroph
July 2, 2015
Bipolar H ii regions - Morphology and star formation in their vicinity
I - G319.88+00.79 and G010.32−00.15
L. Deharveng1 , A. Zavagno1 , M.R. Samal1 , L.D. Anderson2 , G. LeLeu1 , D. Brevot1 , A. Duarte-Cabral3 , S. Molinari4 ,
M. Pestalozzi4 , J.B. Foster5 , J.M. Rathborne6 , and J.M. Jackson7
1
2
3
4
arXiv:1507.00215v1 [astro-ph.GA] 1 Jul 2015
5
6
7
Aix Marseille Universit´e, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France
Department of Physics and Astronomy, West Virginia University, Morgantown, WV 26506, USA
School of Physics and Astronomy, University of Exeter, Stocker Road, Exeter EX4 4QL, UK
INAF-Istituto Fisica Spazio Interplanetario, Via Fosso del Cavaliere 100, I-00133 Roma, Italy
Yale Center for Astronomy and Astrophysics, Yale University, New Haven, CT 06520
CSIRO Astronomy and Space Science, PO Box 76 Epping NSW 1710 Australia
Institute for Astrophysical Research, Boston University, Boston, MA 02215, USA
Received 18/03/2014; Accepted
ABSTRACT
Aims. Our goal is to identify bipolar H ii regions and to understand their morphology, their evolution, and the role they play in the
formation of new generations of stars.
Methods. We use the Spitzer-GLIMPSE, -MIPSGAL, and Herschel-Hi-GAL surveys to identify bipolar H ii regions, looking for
(ionized) lobes extending perpendicular to dense filamentary structures. We search for their exciting star(s) and estimate their distances
using near-IR data from the 2MASS or UKIDSS surveys. Dense molecular clumps are detected using Herschel-SPIRE data, and we
estimate their temperature, column density, mass, and density. MALT90 observations allow us to ascertain their association with the
central H ii region (association based on similar velocities). We identify Class 0/I young stellar objects (YSOs) using their Spitzer and
Herschel-PACS emissions. These methods will be applied to the entire sample of candidate bipolar H ii regions to be presented in a
forthcoming paper.
Results. This paper focuses on two bipolar H ii regions, one that is especially interesting in terms of its morphology, G319.88+00.79,
and one in terms of its star formation, G010.32−00.15. Their exciting clusters are identified and their photometric distances estimated
to be 2.6 kpc and 1.75 kpc, respectively; thus G010.32−00.15 (known as W31 north) lies much closer than previously assumed. We
suggest that these regions formed in dense and flat structures that contain filaments. They have a central ionized region and ionized
lobes extending perpendicular to the parental cloud. The remains of the parental cloud appear as dense (more than 104 cm−3 ) and
cold (14–17 K) condensations. The dust in the photodissociation regions (in regions adjacent to the ionized gas) is warm (19–25 K).
Dense massive clumps are present around the central ionized region. G010.32-00.14 is especially remarkable because five clumps of
several hundred solar masses surround the central H ii region; their peak column density is a few 1023 cm−2 , and the mean density in
their central regions reaches several 105 cm−3 . Four of them contain at least one massive YSO (including an ultracompact H ii region
and a high-luminosity Class I YSO); these clumps also contain extended green objects (EGOs) and Class II methanol masers. This
morphology suggests that the formation of a second generation of massive stars has been triggered by the central bipolar H ii region.
It occurs in the compressed material of the parental cloud.
Key words. ISM: H ii regions – ISM: dust – ISM: individual objects: G010.32−00.15, G319.88+00.79 – Stars: formation
1. Introduction
Hi-GAL, the Herschel infrared GALactic Plane survey (Molinari
et al. 2010a), gives an overall view of the distribution of the
different phases of the Galactic interstellar medium (ISM).
According to Molinari et al. (2010b), “The outstanding feature
emerging from the first images is the impressive and ubiquitous ISM filamentary nature.” Now that all of the Galactic plane
has been surveyed by Hi-GAL we still have the same view of a
very filamentary structure, especially at Herschel-SPIRE wavelengths. Several papers present a preliminary analysis of filaments (for example, Arzoumanian et al. 2011, 2013; Hill et al.
2011; Peretto et al. 2013) and suggest that the filaments are the
main birth site of prestellar cores (for example, Molinari et al.
2010b; Andr´e et al. 2010).
Send offprint requests to: lise.deharveng@oamp.fr
The Spitzer-GLIMPSE survey has shown a “bubbling
Galactic plane” (Churchwell et al. 2006, 2007; Simpson et al.
2012). These bubbles observed at 8.0 µm enclose classical H ii
regions (Deharveng et al. 2010; Anderson et al. 2011). The exact morphology of most of these bubbles is still unknown; are
they 3D spherical structures or 2D rings? The vicinity of 43
bubbles has been observed by Beaumont & Williams (2010) in
the 12 CO (3-2) transition. These observations suggest that these
bubbles are not spherical structures, but rings, and that H ii regions form in flat molecular clouds (with a thickness of a few
parsecs). If massive stars form in filamentary 1D or 2D structures, thus dense filaments or sheets, so in a medium where density gradients are present, many bipolar nebulae should form and
be observed. The Green Bank Telescope H ii Region Discovery
Survey (Bania et al. 2010; resolution 8200 ) shows the opposite
situation, the rarity of bipolar nebulae, hence the conclusion of
Anderson et al. (2011), based on the comparison of the radio
1
Mon. Not. R. Astron. Soc. 000, 000–000 (0000)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00214v1 [astro-ph.GA] 1 Jul 2015
Strong Candidate for AGN Feedback: VLT/X-shooter
Observations of BALQSO SDSS J0831+0354
Carter
Chamberlain,1 Nahum Arav,1 Chris Benn2
1
Department of Physics, Virginia Tech, Blacksburg, VA 24061: carterch@vt.edu
Newton Group, Apartado 321, 38700 Santa Cruz de La Palma, Spain
2 Isaac
2 July 2015
ABSTRACT
We measure the location and energetics of a S iv BALQSO outflow. This ouflow has
a velocity of 10,800 km s−1 and a kinetic luminosity of 1045.7 erg s−1 , which is 5.2%
of the Eddington luminosity of the quasar. From collisional excitation models of the
observed S iv/S iv* absorption troughs, we measure a hydrogen number density of
nH =104.3 cm−3 , which allows us to determine that the outflow is located 110 pc from
the quasar. Since S iv is formed in the same ionization phase as C iv, our results can
be generalized to the ubiquitous C iv BALs. Our accumulated distance measurements
suggest that observed BAL outflows are located much farther away from the central
source than is generally assumed (0.01-0.1 pc).
Key words: galaxies: quasars – galaxies: individual (SDSS J083126.15+035408.0) –
line: formation – quasars: absorption lines
1
INTRODUCTION
Broad absorption line (BAL) outflows are detected as
absorption troughs that are blueshifted in the restframe spectrum of 20–40% of quasars (Hewett & Foltz
2003; Ganguly & Brotherton 2008; Knigge et al. 2008;
Dai, Shankar & Sivakoff 2008). From their detection rate,
we deduce that these outflows cover on average ∼20–40%
of the solid angle around the quasar. Such large opening
angles allow for efficient interaction with the surrounding
medium. As shown by simulations, the mass, momentum and especially the energy carried by these outflows
can play an important role in the evolution of galaxies
and their environments (e.g. Scannapieco & Oh 2004;
Levine & Gnedin 2005; Hopkins et al. 2006; Cattaneo et al.
2009; Ciotti, Ostriker & Proga 2009, 2010; Ostriker et al.
2010; Gilkis & Soker 2012; Choi et al. 2014). Theoretical
studies show that such interactions can provide an explanation for a variety of observations: the self-regulation
of the growth of the supermassive black hole and of the
galactic bulge, curtailing the size of massive galaxies,
and the chemical enrichment of the intergalactic medium.
These processes are part of the so-called AGN feedback
(e.g. Silk & Rees 1998; Di Matteo, Springel & Hernquist
2005;
Germain, Barai & Martel
2009;
Hopkins, Murray & Thompson
2009;
Elvis
2006;
Zubovas & Nayakshin 2014, and references therein).
The importance of quasar outflows to AGN feedback
depends on the mass-flow rate (M˙ ) and kinetic luminosity
(E˙ k ) of the outflowing material. An E˙ k value of at least 0.5%
Table 1. Physical properties of energetic quasar outflows.
˙
v
R
M
Log E˙ k E˙ k /LdEdd
Object
(km/s) (pc) (M⊙ /yr) (erg/s)
(%)
SDSS J0831+03541 -10800 110a
135
45.7
5.2
SDSS J1106+19392
-8250
320a
390
46.0
12
HE 0238-19043
-5000 3400b
140
45.0
0.7
SDSS J0838+29554
-5000 3300c
300
45.4
2.3
SDSS J0318-06005
-4200 6000c
120
44.8
0.13
R from: a high-ionization S iv*/S iv; b high-ionization O iv*/O iv;
c low-ionization Si ii*/Si ii
d See Section 5.1 for L
Edd determination.
References. (1) this work; (2) Borguet et al. 2013; (3)
Arav et al. 2013; (4) Moe et al. 2009; (5) Dunn et al. 2010.
(Hopkins & Elvis 2010) or 5% (Scannapieco & Oh 2004) of
the Eddington luminosity (LEdd ) is deemed sufficient to produce the aforementioned feedback effects. A crucial parameter needed to determine E˙ k is the distance R to the outflow from the central source. Lacking spatial image information, we deduce R from the value of the ionization parameter (UH , see Section 4) once the hydrogen number density (nH ) of the gas is known. Our group has determined
nH (leading to R, M˙ and E˙ k ) for several quasar outflows
(e.g. Moe et al. 2009; Dunn et al. 2010; Bautista et al. 2010;
Aoki et al. 2011; Borguet et al. 2012b) by utilizing absorption lines from excited states of singly-ionized species (e.g.
Fe ii* and Si ii*).
Absorption troughs from singly-ionized species classify
an outflow as a LoBAL. The lower detection rate of LoBALQSO in spectroscopic surveys (3–7%) compared with
20–40% for C iv BALQSO (Dai, Shankar & Sivakoff 2012)
arXiv:1507.00204v1 [astro-ph.IM] 1 Jul 2015
The Microchannel X-ray Telescope on Board the
SVOM Satellite
D. Götz∗a , C. Adamib , S. Basab , V. Beckmannc , V. Burwitzd , R. Chipauxe , B. Cordiera ,
P. Evans f , O. Godetg , R. Goosmannh , N. Meidingerd , A. Meurisa , C. Motchh , K.
Nandrad , P. O’Brien f , J. Osborne f , E. Perinatii , A. Raud , R. Willingale f , K. Mercier j ,
F. Gonzalez j
a CEA
Saclay – DSM/Irfu/Service d’Astrophysique, F-91191, Gif-sur-Yvette, France
Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326,
13388, Marseille, France
c François Arago Centre, APC, Université Paris Diderot, CNRS/IN2P3, CEA/Irfu, Observatoire
de Paris, Sorbonne Paris Cité, 13 Rue Watt, 75013 Paris, France
d MPE, Giessenbachstrasse, D-85748 Garching, Germany
e CEA Saclay – DSM/Irfu/SEDI, F-91191, Gif-sur-Yvette, France
f Dept. of Physics & Astronomy – University of Leicester, Leicester LE1 7RH, GB
g Institut de Recherche en Astrophysique and Planétologie (IRAP), Université de Toulouse, UPS, 9
Avenue du colonel Roche, F-31028 Toulouse Cedex 4, France
h Observatoire astronomique de Strasbourg, Université de Strasbourg, CNRS, UMR 7550, 11 rue
de l’Université, F-67000 Strasbourg, France
i Institut für Astronomie und Astrophysik Tübingen, Eberhard Karls Univ. Tübingen, Germany
j CNES, 18 avenue Edouard Belin, F-31400 Toulouse, France
E-mail: diego.gotz@cea.fr
b Aix
We present the Micro-channel X-ray Telescope (MXT), a new narrow-field (about 1◦ ) telescope
that will be flying on the Sino-French SVOM mission dedicated to Gamma-Ray Burst science,
scheduled for launch in 2021. MXT is based on square micro pore optics (MPOs), coupled with
a low noise CCD. The optics are based on a “Lobster Eye” design, while the CCD is a focal plane
detector similar to the type developed for the seven eROSITA telescopes. MXT is a compact and
light (<35 kg) telescope with a 1 m focal length, and it will provide an effective area of about 45
cm2 on axis at 1 keV. The MXT PSF is expected to be better than 4.2 arc min (FWHM) ensuring a
localization accuracy of the afterglows of the SVOM GRBs to better than 1 arc min (90% c.l. with
no systematics) provided MXT data are collected within 5 minutes after the trigger. The MXT
sensitivity will be adequate to detect the afterglows for almost all the SVOM GRBs as well as to
perform observations of non-GRB astrophysical objects. These performances are fully adapted to
the SVOM science goals, and prove that small and light telescopes can be used for future small
X-ray missions.
Swift: 10 Years of Discovery,
2-5 December 2014
La Sapienza University, Rome, Italy
∗ Speaker.
c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence.
http://pos.sissa.it/
D. Götz
The Microchannel X-ray Telescope
1. Introduction
SVOM (Space based astronomical Variable Object Monitor) is a Sino-French mission dedicated to Gamma-Ray Bursts (GRBs) studies developed in cooperation by the French Space Agency
(CNES), the Chinese National Space Administration (CNSA) and the Chinese Academy of Sciences (CAS). SVOM will be launched in 2021, and will carry two wide field instruments – the
coded mask telescope ECLAIRs (4–150 keV, see [1] and Schanne et al. these proceedings), providing the GRB triggers, and a non-imaging gamma-ray spectrometer GRM (50 keV–5 MeV), and
two narrow field instruments – a visible telescope, VT, and the MXT X-ray telescope, that are
pointed at the ECLAIRs error box after an autonomous platform slew to observe the GRB afterglows. SVOM alerts will be transmitted to ground through a VHF network and will be delivered to
the scientific community within 30 seconds from detection in most of the cases. SVOM will also
include a set of ground based instruments: the GWACs, an ensemble of wide angle cameras that
will continuously monitor the central part of the field of view of ECLAIRs, in order to catch the
GRB prompt emission in the visible band, and two dedicated robotic ground follow-up telescopes
(GFTs). The latter will be provided by the two partners and installed in China and Mexico. For
more information about the SVOM mission, see [2] and Cordier et al. these proceedings.
2. The Telescope and the Optics
The Microchannel X-ray Telescope (MXT) is developed under CNES responsibility in close
collaboration with CEA/Irfu, the University of Leicester, and the Max-Planck-Institut für Extraterrestrische Physik (MPE). It will be composed of the main subsystems indicated in Fig. 1, namely
an optic based on square micro-channels, a camera, a carbon fibre structure, and a radiator.
The optical scheme of the MXT is based on a “Lobster Eye” grazing incidence geometry as
defined by [3], and optimised for narrow-field use. Reflections occur within the square pores of
square micro-pore glass plates (MPOs), which have a thickness of 1-2 mm and a side length of 40
mm; these plates, produced by Photonis, are spherically slumped to a radius twice the focal length
of the MXT.
The MPOs will be coated with Ir and mounted on a round frame, see Fig. 2, in order to obtain
a complete X-ray optic with a 1 m focal length and coverage of the 64×64 arc min field of view
defined by the detector. The front of the optic will be coated with 70 nm of Al for thermal reasons
and in order to reduce the optical load on the detector.
The “Lobster Eye” optic was originally developed in order to cover very large areas of the sky
with homogeneous sensitivity. Thanks to its light weight with respect to classical X-ray optics, it
is very well adapted to mini class satellites like SVOM. Its working principle is sketched in Fig. 3.
X-rays entering in the MPOs can either be reflected twice and focused in the central PSF spot, or
reflected just once and focused in the PSF arms. 25% of the incident X-ray flux is focussed in the
central spot, 2×25% in the arms, and the rest in a diffuse patch, see Fig. 3. Thanks to the “Lobster
Eye” geometry the vignetting is very low, of the order of 15% at the edge of the FOV. Current
simulations indicate that the PSF FWHM of such a system is 4.2 arc min at 1 keV.
2
Mon. Not. R. Astron. Soc. 000, 1–12 (0000)
Printed 2 July 2015
(MN LATEX style file v2.2)
The nature of voids: II. Tracing underdensities with biased galaxies
S.
Nadathur1? & S. Hotchkiss2
1
Department of Physics, University of Helsinki and Helsinki Institute of Physics, P.O. Box 64, FIN-00014, University of Helsinki, Finland
of Physics and Astronomy, University of Sussex, Falmer, Brighton, BN1 9QH, UK
2 Department
arXiv:1507.00197v1 [astro-ph.CO] 1 Jul 2015
2 July 2015
ABSTRACT
Cosmic voids may be very useful in testing fundamental aspects of cosmology. Yet observationally voids can only be seen as regions with a deficit of bright galaxies. To study how
biased galaxies trace matter underdensities and how the properties of voids depend on those
of the tracer galaxy population, we use a ΛCDM N-body simulation populated with mock
galaxies based on the halo occupation distribution (HOD) model. We identify voids in these
mocks using the ZOBOV void finder and measure their abundances, sizes, tracer densities, and
dark matter content. To separate the effects of bias from those of sampling density, we do
the same for voids traced by randomly down-sampled subsets of the dark matter particles in
the simulation. We find that galaxy bias reduces the total number of voids by ∼ 50% and can
dramatically change their size distribution. The matter content of voids in biased and unbiased
tracers also differs. Using simulations to accurately estimate the cosmological constraints that
can be obtained from voids therefore requires the use of realistic mock galaxy catalogues.
We discuss aspects of the dark matter content of voids that can be deduced from properties
of the tracer distribution, such as the void size and the minimum tracer number density. In
particular we consider the compensation of the total mass deficit in voids and find that the
distinction between over- and under-compensated voids is not a function of void size alone,
as has previously been suggested. However, we find a simple linear relationship between the
average density of tracers in the void and the total mass compensation on much larger scales.
The existence of this linear relationship holds independent of the bias and sampling density of
the tracers. This provides a universal tool to classify void environments and will be important
for the use of voids in observational cosmology.
Key words: cosmology: observations – large-scale structure of Universe – methods: numerical – methods: data analysis
1
INTRODUCTION
Several recent studies have advocated the use of large underdensities in the matter distribution, known as cosmic voids, which can
be important tools for cosmology in a number of different contexts. These include using voids to test modified gravity scenarios
(Li et al. 2012; Clampitt et al. 2013; Cai et al. 2014; Zivick et al.
2014), to probe dark energy (Ryden 1995; Lee & Park 2009; Bos
et al. 2012; Sutter et al. 2015; Pisani et al. 2015), or to detect the effects of warm dark matter (Tikhonov et al. 2009; Yang et al. 2014).
These tests make use of observables such as the abundances, size
distributions, shapes and radial density profiles of voids. This information will be available for very large numbers of voids in future
large-scale survey data such as Euclid (Laureijs et al. 2011).
Weak lensing by voids (Krause et al. 2013; Melchior et al.
2014), and their possible gravitational effects on the CMB via
the integrated Sachs-Wolfe (ISW) effect have also been studied
?
seshadri.nadathur@helsinki.fi
c 0000 RAS
(Granett et al. 2008; Cai et al. 2014; Hotchkiss et al. 2015; Planck
Collaboration et al. 2015). In principle for rare voids such studies probe the extreme tails of the distribution of fluctuations of the
gravitational potential so may provide unique insights. Indeed the
case of possible ISW imprints of voids is particularly intriguing,
since theoretical considerations and simulation results show that
any such detection would be in strong tension with the standard
ΛCDM cosmology (Nadathur et al. 2012; Flender et al. 2013; Cai
et al. 2014; Hotchkiss et al. 2015; Aiola et al. 2015).
Several catalogues of voids found in galaxy redshift surveys
have already been published (Pan et al. 2012; Sutter et al. 2012; Nadathur & Hotchkiss 2014), and more will be available from future
data releases. These catalogues can be used for the various studies described above. However, since the dark matter distribution is
not directly observable, the comparison of observations with theoretical expectation necessary for precision cosmology requires a
detailed understanding of how biased galaxies actually trace matter
underdensities.
A common theoretical model of voids is that based on the ex-
Using the transit of Venus to probe the upper planetary
atmosphere
Fabio Reale1,2 , Angelo F. Gambino1, Giuseppina Micela2 , Antonio Maggio2 , Thomas Widemann3 ,
arXiv:1507.00195v1 [astro-ph.EP] 1 Jul 2015
Giuseppe Piccioni4
1
Dipartimento di Fisica e Chimica, Universit`a di Palermo, Piazza del Parlamento 1, 90134
Palermo, Italy
2
INAF/Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
3
Universit`e de Versailles-Saint-Quentin - ESR/DYPAC EA 2449, Observatoire de Paris - LESIA -
UMR CNRS 8109 5, place Jules-Janssen, 92190 Meudon, France
4
INAF-IAPS (Istituto di Astrofisica e Planetologia Spaziali), via del Fosso del Cavaliere 100,
00133 Roma, Italy
The atmosphere of a transiting planet shields the stellar radiation providing us with a powerful method to estimate its size and density. In particular, because of their high ionization
energy, atoms with high atomic number (Z) absorb short-wavelength radiation in the upper atmosphere, undetectable with observations in visible light. One implication is that the
planet should appear larger during a primary transit observed in high energy bands than in
the optical band. The last Venus transit in 2012 offered a unique opportunity to study this
effect. The transit has been monitored by solar space observations from Hinode and Solar
Dynamics Observatory (SDO). We measure the radius of Venus during the transit in three
˚ UV (1600A,
˚ 1700A),
˚ Extreme Uldifferent bands with subpixel accuracy: optical (4500A),
1
than the optical radius, while the altitude in the UV band is in between (40 − 50 km).
The altitude that we measured in the EUV and X-ray bands corresponds to the height where
the optical thickness along the tangential line of sight reaches the value τ = 1. At this altitude
the solar radiation is mostly absorbed by photoionization of neutral atoms, and hence this is also
the altitude where a peak of the electron density is expected. In particular, the F1 peak in density
profiles of planetary atmospheres is associated to the absorption of EUV radiation 8 , while UV
photons and soft X-rays are absorbed in slightly deeper layers. In fact, we have verified that
the different Venus sizes indicated by UV and EUV data can be explained in term of tangential
column densities at the retrieved altitudes (165-170 km for EUV, and 125-135 km for UV) at
terminator between the day and night hemispheres (solar zenith angle, SZA, 90◦ ). The absorption
path length, l, at a given altitude h can be expressed as li ∼
p
2(Rsf + h)Hi , where Hi is the
atmospheric scale height of each molecular species i, and the optical thickness can be evaluated as
τ (λ) =
P
ni σi (λ)li , where ni and σi (λ) are the molecular density and absorption cross section,
respectively. At 170 km CO2 and CO have approximately the same mixing ratio 21 , while at 130
km CO2 dominates. A simple calculation assuming a spherical geometry shows that the neutral
atmosphere is ∼ 20 times more opaque in the EUV than in the UV at 170 km, and ∼ 80 times at
130 km in the regions probed.
Fig. 2 shows a comparison between our ionospheric altitudes measured at different wavelengths with the prediction from a detailed model that includes also a dependence on the SZA 8 .
There is a very good general agreement, and an accurate one at UV wavelength. However, we
5
DOWNSIZING FROM THE POINT OF VIEW OF MERGING MODEL
(preliminary discussion)
V.M. Kontorovich1,2
1
Institute of Radio Astronomy NASU
4 Chervonopraporna Str., Kharkov, 61002, Ukraine,
2
Kharkov National V.N. Karazin University
4, Svobody square, Kharkov 61022, Ukraine,
E-mail:vkont@rian.kharkov.ua
In four-particle processes оf scattering with transfer of mass, unlike mergers in which mass can only increase,
mass of the most massive galaxies may be reduced. Elementary model describing such process is considered. In this
way, it is supposed to explain observed phenomenon of "dawnsizing" when increasing of characteristic mass the
heaviest galaxies over cosmological time replaces by its reduction.
PACS: 98.65.Fz; 98.80.Bp
1. Introduction
Quite unusual in terms of the paradigm of mergers, but
has long discussed the fact that the maximal galaxy
masses (Shechter parameter ), which grow up with
decreasing of the red shift at large distances, begin to
decrease as you get closer to the present time (see.
Figure 1), seems be in contrary to the model of mergers.
We will show that it is not.
manifests itself as unlimited growth of the MF as we
approach the time of the "explosion" t=tcr. Asymptotics
of solution K(M, t) for the modifide MF (mMF, App B)
u
F  M f near the singularity (outside the physical
area of the power behavior) has the form (for u=2)
characteristical maximal masses).
K ( M , t) 

,
1
1

M M max (t )
1
M max (t ) 
(1)
c (tcr  t )
Figure 1. Observational data of the Hubble ultra deep
field [1] relating to the Shechter parameter М* (in our
consideration it is the relevant characteristical maximal
masses).
In the model of galaxy mergers, built on the basis of
Smoluchowski kinetic equation (KE), are taken into
account only processes of the (paired) mergers, that is
the processes involving three "particles" (App.A, Figure
2). The resulting solutions (App.B) allow to find the
slope of the mass function f ( M , t ) in a wide range of
redshifts [2,3], satisfactorily explain the observational
data of the Hubble ultra deep field [1] (the evolution of
the slopes until limiting redshifts). However, arising at
this "explosive" evolution leads to an unlimited growth
of the maximal mass as we approach the time of the
"explosion" [2]. "Explosive" singularity in the solution
That artefact, associated with the use of instant  shaped source in the Smoluchowski KE , possible to
avoid by the obvious physical regularization [2,3], the
meaning of which is in taking into account a finite rise
time of gravitational instability, which leads to
separation galaxies from the general expansion of the
universe. Mathematically, this was taken into account
by blur of  -function on the right side of KE and its
replacement by a  -shaped stepping stone with a finite
small width  . The values of MF remain finite in the
region of maximal masses too. But still, the very
maximal mass in regularized solutions also increases
indefinitely when approaching the moment of the
explosion [2-4].
As in other similar tasks, taking into account
three-particle processes (in our case - mergers of
galaxies), leading to the explosive evolution, the finite
results occur when four-particle processes come into
play in the vicinity of singularity and, in our case,
describe the scattering with transfer of mass. In this
case, unlike mergers (Figure 2), in which the mass can
only increase, a significant role is played the scattering
processes in which the mass of the most massive
galaxies may be decrease (Figure 3, 4). Below, we
consider a simple model scheme describing the
disaggregation. In this way, it is supposed to explain the
observed "downsizing" phenomenon (Figure 1), where
increasing the heaviest characteristic mass over time
replaced of its decreasing.
c
ESO
2015
Astronomy & Astrophysics manuscript no. Laur2015
July 2, 2015
Period change of massive binaries from combined photometric
and spectroscopic data in Cygnus OB2
Jaan Laur1, 2 , Elmo Tempel1 , Taavi Tuvikene1 , Tõnis Eenmäe1 , and Indrek Kolka1
1
2
Tartu Observatory, Observatooriumi 1, 61602 Tõravere, Estonia
e-mail: jaan.laur@to.ee
Institute of Physics, University of Tartu, Ravila 14c, 51010 Tartu, Estonia
arXiv:1507.00180v1 [astro-ph.SR] 1 Jul 2015
Received May 12, 2015; accepted June 29, 2015
ABSTRACT
Context. Mass loss is an important property in evolution models of massive stars. As up to 90% of the massive stars have a visual
or spectroscopic companion and many of them exhibit mass exchange, mass-loss rates can be acquired through the period study of
massive binaries.
Aims. Using our own photometric observations as well as archival data, we look for variations in orbital periods of seven massive
eclipsing binary systems in the Cygnus OB2 association and estimate their mass-loss rates and stellar parameters.
Methods. We use a Bayesian parameter estimation method to simultaneously fit the period and period change to all available data and
a stellar modelling tool to model the binary parameters from photometric and radial-velocity data.
Results. Four out of the seven selected binaries show non-zero period change values at two-sigma confidence level. We also report
for the first time the eclipsing nature of a star MT059.
Key words. binaries: eclipsing – stars: mass-loss – open clusters: individual: Cygnus OB2 – stars: massive
1. Introduction
˙ in a binary system has been comThe change of a period (P)
monly evaluated by measuring the times of minimum of primary eclipses and then comparing those times to the calculated
ephemerides for the eclipses. The period change value is then acquired by fitting a polynomial model through the observed minus
calculated (O − C) data points. The value of this period change
can indicate mass loss through stellar winds or mass transfer between the binary components (Singh & Chaubey 1986).
Mass loss in binary systems comes mostly from the mass
transfer via Roche lobe overflow when the primary star’s photosphere crosses the inner Lagrange point and starts to exchange
mass with its companion. While the mass is transferred, the centre of mass of the interacting binary system shifts towards the
mass-receiving star. Due to the conservation of angular momentum, this mass transfer causes a change in the system’s orbital
period: the period increases while the mass is transferred to the
more massive star and vice versa. The period also increases when
mass is lost from the system through stellar winds.
Sana et al. (2012) studied a sample of massive O stars and
found that over 70% of them exchange mass with a companion. The speckle interferometric survey of Galactic massive stars
(Mason et al. 2009) expanded upon a decade-old survey and
noted that 75% of the massive stars in clusters and 59% in
the field have either a spectroscopic or a visual companion. As
low inclination or long periods make the binary systems hard
to detect, these values are only the lower limit and, when corrected to the observational bias, the upper limit for binarity fraction in bound systems can reach as high as ∼ 90 ± 10%
(Kiminki & Kobulnicky 2012).
Massive stars have a relatively short lifetime that is influenced by their mass-loss rates. This mass loss determines the
evolution of their temperature and luminosity, and also the supernova explosion in the final evolution phases. Massive stars
provide feedback in the form of ultraviolet (UV) radiation, stellar
winds, and supernovae which impact the evolution of disc galaxies by mixing and driving turbulence of interstellar gas and triggering star formation. They also enrich the interstellar medium
with heavy elements and drive the chemical evolution of galaxies
and the metallicity evolution of the universe. Therefore, knowing
the evolution of massive stars gives us insight into the evolution
of galaxies. The uncertainties in determining the mass-loss rates
for massive stars are still large even in the local Universe (Smith
2014).
Constraining orbital parameters for massive binary systems
is important to test the predictions of theoretical models. In order to determine all the parameters of a single binary system,
one needs to obtain photometric and spectroscopic data at various orbital phases, requiring great observational effort. As it is
time-consuming to obtain those observations and as the eclipsing nature is dependent on the inclination of the system, there
are not many massive binaries with fully constrained parameters
known. In this work we present the stellar parameters as well as
the period and the period change values of seven selected systems in Cygnus OB2.
Stellar association Cygnus OB2 is one of the most massive and richest associations known in our Galaxy, containing at
least 60 O-type stars (Negueruela et al. 2008). Determined from
the stellar evolution models, the age of the cluster is found to
be 1–7 Myr (Wright et al. 2015). It is dynamically young, substructured and unbound association that is located in the molecular cloud complex Cygnus X and its total estimated mass is
(1.6 ± 0.3) × 104 M⊙ (Wright et al. 2014).
Article number, page 1 of 12
Astronomy & Astrophysics manuscript no. Mathis2015-p-acc
July 2, 2015
c
ESO
2015
Letter to the Editor
The variation of tidal dissipation in the convective envelope
of low-mass stars along their evolution
S. Mathis1
Laboratoire AIM Paris-Saclay, CEA/DSM - CNRS - Université Paris Diderot, IRFU/SAp Centre de Saclay, F-91191 Gif-sur-Yvette
Cedex, France e-mail: stephane.mathis@cea.fr
Received ... / accepted ...
arXiv:1507.00165v1 [astro-ph.SR] 1 Jul 2015
ABSTRACT
Context. Since 1995, more than 1500 exoplanets have been discovered around a large diversity of host stars (from M- to A-type stars).
Tidal dissipation in stellar convective envelopes is a key actor that shapes the orbital architecture of short-period systems.
Aims. Our objective is to understand and evaluate how tidal dissipation in the convective envelope of low-mass stars (from M to F
types) depends on their mass, evolutionary stage and rotation.
Methods. Using a simplified two-layer assumption, we compute analytically the frequency-averaged tidal dissipation in their convective envelope. This dissipation is due to the conversion into heat of the kinetic energy of tidal non wave-like/equilibrium flow and
inertial waves because of the viscous friction applied by turbulent convection. Using grids of stellar models allows us to study the
variation of the dissipation as a function of stellar mass and age on the Pre-Main-Sequence and on the Main-Sequence for stars with
masses spanning from 0.4 to 1.4M .
Results. During their Pre-Main-Sequence, all low-mass stars have an increase of the frequency-averaged tidal dissipation for a fixed
angular velocity in their convective envelope until they reach a critical aspect and mass ratios (respectively α = Rc /Rs and β = Mc /Ms
where Rs , Ms , Rc and Mc are the star’s radius and mass and its radiative core’s top radius and mass). Next, the dissipation evolves
on the Main Sequence to an asymptotic value that becomes maximum for 0.6M K-type stars and that decreases by several orders
of magnitude with increasing stellar mass. Finally, the rotational evolution of low-mass stars strengthens the importance of tidal
dissipation during the Pre-Main-Sequence for star-planet and multiple star systems.
Conclusions. As shown by observations, tidal dissipation in stars varies over several orders of magnitude as a function of stellar
mass, age and rotation. We demonstrate that: i) it reaches a maximum value on the Pre-Main-Sequence for all stellar masses; ii) on
the Main-Sequence and at fixed angular velocity, it becomes maximum for 0.6M K-type stars and it decreases with increasing mass.
Key words. hydrodynamics – waves – celestial mechanics – planet-star interactions – stars: evolution – stars: rotation
1. Introduction and context
Since 1995, more than 1500 exoplanets have been discovered
around a large diversity of host stars (e.g. Perryman 2011).
As of today, their mass range spreads from M red dwarfs to
intermediate-mass A-type stars. In this context, tidal dissipation
in the host star has a strong impact on the orbital configuration
of short-period systems. Indeed, the dissipation in rotating turbulent convective envelopes of low-mass stars is believed to play
a key-role for tidal migration, orbits’ circularization and spins’
alignment (e.g. Albrecht et al. 2012; Lai 2012; Ogilvie 2014,
and references therein for hot-Jupiter systems). In these regions,
this is the turbulent friction acting on tidal flows, which dissipates their kinetic energy into heat, that drives the tidal evolution
of star-planet systems. In stellar convection zones, tidal flows
are constituted by a large-scale non-wave like/equilibrium flow
driven by the hydrostatic adjustment of the stellar structure because of the presence of the planetary/stellar companion (Zahn
1966; Remus et al. 2012) and the dynamical tide constituted by
inertial waves, their restoring force being the Coriolis acceleration (e.g. Ogilvie & Lin 2007). In this framework, both the structure and rotation of stars strongly varies along their evolution
(e.g. Siess et al. 2000; Gallet & Bouvier 2013, 2015). Moreover,
as reported by Ogilvie (2014), observations of star-planet and binary star systems show that tidal dissipation varies over several
orders of magnitude. Therefore, one of the key questions that
must be addressed is: how does tidal dissipation in the convective envelope of low-mass stars varies as a function of stellar
mass, evolutionary stage and rotation?
In this work, we study the variations of the frequencyaveraged tidal dissipation in stellar convective envelopes as a
function of the mass, the age and the rotation of stars. In sec.
2, we introduce the assumptions and the formalism that allows
us to analytically evaluate this quantity as a function of the structure and rotation of stars (Ogilvie 2013). In sec. 3, we compute
it as a function of stellar mass and evolutionary stage at fixed angular velocity using grids of stellar models for stars from 0.4 to
1.4M . Next, we discuss the impact of the rotational evolution of
stars. In sec. 4, we present our conclusions and the perspectives
of this work.
2. Tidal dissipation modelling
To analytically evaluate the frequency-averaged tidal dissipation
in the convective envelope of low-mass stars, we here choose to
adopt a simplified two-layer model as in Ogilvie (2013) (see also
Penev et al. 2014, and fig. 1). In this modeling, both the radiative
core and the convective envelope are assumed to be homogeneous with respective constant densities ρc and ρe for the sake of
simplicity. This model features a star A of mass Ms and mean raArticle number, page 1 of 5
Astronomy & Astrophysics manuscript no. pap
July 2, 2015
c
ESO
2015
The naked nuclei of LINERs
Barbara Balmaverde1, 2 and Alessandro Capetti3
1
2
3
Dipartimento di Fisica e Astronomia, Università di Firenze, via G. Sansone 1, 50019 Sesto Fiorentino (Firenze), Italy
INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, I-50125 Firenze, Italy
INAF - Osservatorio Astrofisico di Torino, Via Osservatorio 20, I-10025 Pino Torinese, Italy
arXiv:1507.00164v1 [astro-ph.GA] 1 Jul 2015
ABSTRACT
We analyze HST spectra and Chandra observations of a sample of 21 LINERs, at least 18 of which genuine AGN. We find a correlation
between the X-rays and emission lines luminosities, extending over three orders of magnitude and with a dispersion of 0.36 dex; no
differences emerge between LINERs with and without broad lines, or between radio-loud and radio-quiet sources. The presence
of such a strong correlation is remarkable considering that for half of the sample the X-ray luminosity can not be corrected for
local absorption. This connection is readily understood since the X-ray light is associated with the same source producing the
ionizing photons at the origin of the line emission. This implies that we have a direct view of the LINERs nuclei in the X-rays:
the circumnuclear, high column density structure (the torus) is absent in these sources. Such a conclusion is also supported by midinfrared data. We suggest that this is due to the general paucity of gas and dust in their nuclear regions that causes also their low rate
of accretion and low bolometric luminosity.
Key words. galaxies: active – galaxies: nuclei
1. Introduction
The presence of low ionization nuclear emission line regions
(LINERs), initially identified by Heckman (1980), is very common in galaxies but their nature is still matter of debate. This
is because several processes, beside an active nucleus, produce
a low ionization spectrum, such as shocks (as initially proposed
by Heckman 1980) or evolved stellar populations (Trinchieri &
di Serego Alighieri 1991 and Binette et al. 1994). In high spatial resolution X-ray images, about 60% of LINERs show a nuclear compact source (González-Martín et al. 2009). However
the emission from the SMBH can account for only ∼ 60% of
the ionizing photons and therefore the AGN emission alone is
not enough to power the optical emission lines (Flohic et al.
2006, Eracleous et al. 2010). These results are confirmed by the
CALIFA survey: the integral field spectroscopy images show the
presence of extended regions characterized by a low ionization
spectrum; however, the radial emission-line surface brightness
profiles are inconsistent with ionization by a central point-source
and this is most likely due to the effects of post-AGB stars (Singh
et al. 2013). It is then unclear whether LINERs (or which portion of this heterogeneous class) are powered by accretion onto
a supermassive black hole.
Even considering the LINERs that are genuine active galactic nuclei (AGN), we then still need to understand how they compare to the general AGN population. For example, whether the
orientation based unified model (e.g. Antonucci 1993) applies
also to LINERs, i.e., if the detection of a broad line region (BLR)
is related to the effects of a circumnuclear obscuring structure
(Balmaverde & Capetti 2014).
A widely used tool to investigate the nature of AGN is the ratio between the X-rays and the line luminosities. Maiolino et al.
(1998) found that the LX /L[OIII] ratio differs by a factor ∼100 in
Seyferts type 1 and type 2, the former class showing the largest
Send offprint requests to: balmaverd@arcetri.inaf.it
values. This is interpreted as evidence of nuclear absorption in
the X-rays emission in Seyferts 2, since the [O III] emission is
a robust (and isotropic) estimator of the AGN bolometric luminosity. The study by Risaliti et al. (1999) indicate that in the
Seyferts 2 the local column density, NH,z , exceeds 1023 cm−2 in
all but two of 35 sources considered. Only when absorption is
taken into account, Seyferts 2 show multiwavelength properties
analogous to those of Seyferts 1. For example, Krabbe et al.
(2001) found a tight correlation between the 10 µm (similarly to
the emission lines, the MIR light is radiated almost isotropically)
and the X-ray luminosities, followed by both Seyferts 1 and 2.
In most LINERs the quality of the X-ray spectra is insufficient
to measure NH,z directly and they span a very large wide range
of LX /L[OIII] values (González-Martín et al. 2009), even broader
than in Seyfert galaxies. We still need to establish how LINERs
fit into this scheme.
These issues can be explored in greater depth by studying
the nuclear emission in LINERs at the highest available spatial
resolution. Indeed, a key problem with LINERs is that they are
faint objects and thence subject to strong contamination from the
host galaxy. In this Paper we perform an analysis on a sample
of LINERs based on HST and Chandra data. Our aim is to isolate their genuine nuclear contribution in the X-rays and to perform the best possible decomposition between the galactic and
(if present) the AGN emission lines.
2. Sample selection and data analysis
We consider the 60 galaxies robustly classified as LINERs by
Ho et al. (1997) from their spectroscopic survey of bright nearby
galaxies (Filippenko & Sargent 1985; Ho et al. 1995). We select the objects for which both HST/STIS spectroscopic data and
Chandra observations are available. We find 21 objects, having
discarded 5 sources (namely NGC 4261, NGC 4314, NGC 4594,
NGC 4138 and NGC 4374) due to the poor quality of the HST
Article number, page 1 of 5
Mon. Not. R. Astron. Soc. 000, 1–?? (2002)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00149v1 [astro-ph.HE] 1 Jul 2015
K2 and MAXI observations of Sco X-1 - Evidence for disc
precession?
Pasi
Hakala1⋆ , Gavin Ramsay,2, Thomas Barclay3,4, Phil Charles5
1
Finnish Centre for Astronomy with ESO (FINCA), V¨
ais¨
al¨
antie 20, University Of Turku, FIN-21500 Piikki¨
o, Finland.
Observatory, College Hill, Armagh, BT61 9DG, UK
3 NASA Ames Research Center, M/S 244-40, Moffett Field, CA 94035, USA
4 Bay Area Environmental Research Institute, Inc., 560 Third St. West, Sonoma, CA 95476, USA
5 School of Physics and Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK
2 Armagh
ABSTRACT
Sco X-1 is the archetypal low mass X-ray binary (LMXB) and the brightest persistent extra-solar X-ray source in the sky. It was included in the K2 Campaign 2
field and was observed continuously for 71 days with 1 minute time resolution. In this
paper we report these results and underline the potential of K2 for similar observations of other accreting compact binaries. We reconfirm that Sco X-1 shows a bimodal
distribution of optical ”high” and ”low” states and rapid transitions between them on
timescales less than 3 hours (or 0.15 orbits). We also find evidence that this behaviour
has a typical systemic timescale of 4.8 days, which we interpret as a possible disc
precession period in the system. Finally, we confirm the complex optical vs. X-ray
correlation/anticorrelation behaviour for ”high” and ”low” optical states respectively.
Key words: accretion, accretion discs – X-ray binaries - X-rays: individual: Sco X-1
1
INTRODUCTION
Sco X-1 is the brightest persistent extra-solar X-ray source
¯ and is a low mass X-ray binary
in the sky (Morrison 1967)
(LMXB), where a Roche-lobe-filling secondary star is losing
matter that is accreted by a neutron star via an accretion
disc (Charles & Coe 2006). There are numerous studies on
this 0.787 d orbital period system over the last fifty years
covering almost the entire electromagnetic spectrum. It is
also predicted to be a strong source of gravitational waves
(e.g. Aasi et al. 2014).
As the prototype LMXB, Sco X-1 has been studied for
decades at all wavelengths, but especially at optical and Xray wavelengths. Ilovaisky et al. (1980) found that the optical and X-ray flux was well correlated, especially when in
a bright state, whilst a longer series of optical observations
showed that Sco X-1 changes from a high to low state on
a timescale of a few hours (Hiltner & Mook, 1967). McNamara et al. (2003) presented a year-long optical study of Sco
X-1 and concluded that the optical observations can be accounted for by variations in the mass accretion rate. Sco X-1
is unique amongst LMXBs, because to our knowledge, it is
the only one showing clearly bimodal optical states (Hiltner
& Mook, 1967).
⋆
E-mail:pahakala@utu.fi
c 2002 RAS
The optical light curves of LMXBs are complex and
changes are seen on short time-scales (see for instance Homer
et al. 2001 and Hakala et al. 2009). The interpretation of
these light curves is made more difficult because of data
gaps, the limited duration and cadence of the observations.
Much time and resources have been spent in this endeavour
(e.g. Shih, Charles & Cornelisse 2011 to name but one).
The possibility of uninterrupted photometric observations of
sources by Kepler has therefore led to dramatic discoveries
and results in the field of exo-planets, asteroseismology and
accretion physics. Whilst the original Kepler field contained
several dozen accreting cataclysmic variables (e.g. Howell et
al 2013), it did not contain any known LMXB.
Since the loss of two of its reaction wheels, Kepler has
been re-purposed into K2, and is now making a series of observations of fields along the ecliptic plane, where each field
is observed for ∼70 days (Howell et al. 2014). The fact that
these fields go through the Galactic plane allows the study of
types of objects which were not present in the Kepler field.
In the Campaign 2 field, Sco X-1 was included as one of
the target sources. This paper presents the K2 observations
of Sco X-1 and also simultaneous X-ray data using MAXI.
The unprecedented set of optical observations allows us to
characterise the optical behaviour of Sco X-1 in a way not
previously possible.
On the Parallel and Perpendicular Propagating Motions Visible in
Polar Plumes: An Incubator For (Fast) Solar Wind Acceleration?
arXiv:1507.00143v1 [astro-ph.SR] 1 Jul 2015
Jiajia Liu1,2 , Scott W. McIntosh2 , Ineke De Moortel3 , Yuming Wang1
ABSTRACT
We combine observations of the Coronal Multi-channel Polarimeter (CoMP) and
the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory
(SDO) to study the characteristic properties of (propagating) Alfv´enic motions and
quasi-periodic intensity disturbances in polar plumes. This unique combination of instruments highlights the physical richness of the processes taking place at the base of
the (fast) solar wind. The (parallel) intensity perturbations with intensity enhancements around 1% have an apparent speed of 120 km/s (in both the 171 ˚
A and 193 ˚
A
passbands) and a periodicity of 15 minutes, while the (perpendicular) Alfv´enic wave
motions have a velocity amplitude of 0.5 km/s, a phase speed of 830 km/s, and a shorter
period of 5 minutes on the same structures. These observations illustrate a scenario
where the excited Alfv´enic motions are propagating along an inhomogeneously loaded
magnetic field structure such that the combination could be a potential progenitor of
the magnetohydrodynamic turbulence required to accelerate the fast solar wind.
1.
Introduction
Polar plumes are one of the most striking features in polar coronal holes (e.g., Newkirk and Harvey
1968). Indeed, their relatively long lifetime and high contrast compared to the background coronal hole made them a favorite target of the SOHO era (e.g., Wilhelm 2000;
Banerjee et al. 2009). They are thought of as sources of dense plasma in the fast solar wind
(e.g., Gabriel et al. 2003) that result from the relentless magneto-convective forcing of the
upper solar atmospheric plasma (e.g., Wang 1998).
(Quasi-)periodic upward-propagating intensity perturbations have been observed in various regions of the Sun (e.g., Lites et al. 1999; Yamauchi et al. 2003; Banerjee et al. 2009;
1
Earth and Space Science School, University of Science and Technology of China, NO. 96, JinZhai Road,
Hefei, China
2
High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, CO
80307, USA
3
School of Mathematics and Statistics, University of St Andrews, St Andrews, Fife, KY16 9SS, UK
accepted by ApJ; July 2, 2015
Preprint typeset using LATEX style emulateapj v. 12/16/11
FOURIER ANALYSIS OF BLAZAR VARIABILITY: KLEIN-NISHINA EFFECTS AND THE JET SCATTERING
ENVIRONMENT
Justin D. Finke
U.S. Naval Research Laboratory, Code 7653, 4555 Overlook Ave. SW, Washington, DC, 20375-5352
arXiv:1507.00132v1 [astro-ph.HE] 1 Jul 2015
Peter A. Becker
School of Physics, Astronomy, and Computational Sciences, MS 5C3, George Mason University, 4400 University Drive, Fairfax, VA
22030
accepted by ApJ; July 2, 2015
ABSTRACT
The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier
frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing
the PSDs and time lags produced via a combination of stochastic particle injection and emission via
the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism
used the Thomson cross section and simple δ-function approximations to model the synchrotron and
Compton emissivities. Here we expand upon this work, using the full Compton cross section and
detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed
using the δ-function approximations and those computed using the accurate expressions, provided
the observed photons are produced primarily by electrons with energies exceeding the lower limit of
the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the
cooling timescales of the electrons primarily responsible for the observed emission, and the associated
time lags are related to the difference in electron cooling timescales between the two energy channels,
as expected. If the electron cooling timescales can be determined from the observed time lags and/or
the observed EC PSDs, then one could in principle use the method developed here to determine the
energy of the external seed photon source for EC, which is an important unsolved problem in blazar
physics.
Subject headings: BL Lacertae objects: general — quasars: general — radiation mechanisms: nonthermal — galaxies: active — galaxies: jets
1. INTRODUCTION
Blazars, active galactic nuclei (AGN) with jets aligned
with our line of sight, relativistically beam nonthermal radiation towards the Earth. This emission extends across the electromagnetic spectrum, from radio
to γ rays. Although the observed radio emission probably comes from the superposition of many self-absorbed
synchrotron components spread out over several parsecs (Konigl 1981), the location along the jet where the
higher-frequency emission originates is a matter of some
controversy. This is particularly true for the most powerful class of blazars, the flat spectrum radio quasars (FSRQs) with strong broad emission lines. These objects
are thought to emit γ-rays primarily through Compton
scattering of an external radiation field (known as external Compton or EC). However, the primary external
radiation field for EC is not known, and could be provided by the accretion disk (Dermer & Schlickeiser 1993,
2002), the broad line region (BLR; Sikora et al. 1994),
or the dust torus (Kataoka et al. 1999; Bla˙zejowski et al.
2000). This uncertainty in the seed photon source reflects
the uncertainty in the location along the jet of the primary γ-ray emitting region, which one would expect to
be within ∼ 0.1 pc from the black hole (BH) if the scattering of disk radiation dominates, or between ∼ 0.1 and
justin.finke@nrl.navy.mil
pbecker@gmu.edu
∼ 1 pc if the scattering of BLR emission dominates, or
& 1 pc if the scattering of dust torus emission dominates
(Ghisellini & Tavecchio 2009; Sikora et al. 2009).
The association of γ-ray flares with the ejection of superluminal components (e.g., Marscher et al. 2012) indicates that the flares are coincident with the 43 GHz
core, which is likely located at & a few pc from the black
hole, outside the BLR (although see Nalewajko et al.
2014). The detection of & 100 GeV γ rays from FSRQs also indicates the γ rays must originate from inside the BLR to avoid γγ absorption with BLR photons
(e.g., Aleksi´c et al. 2011). However, rapid γ-ray variability observed in FSRQs such as 3C 454.3 (Tavecchio et al.
2010), PKS 1510−089 (Brown 2013; Saito et al. 2013)
and 4C 21.35 (also known as PKS 1222+21; Aleksi´c et al.
2011) limits the size of the emitting region, and, assuming the emitting region takes up the entire cross section of
a conical jet, implies that it should be within ∼ 0.1 pc of
the black hole, and hence inside the BLR. Dotson et al.
(2012) point out a way to distinguish between these scenarios based on the different cooling timescales at two
energies if the seed photons are from the BLR (e.g. with
dimensionless energy for Lyα ǫ0 = E0 /me c2 = 2 × 10−5 )
or from the dust torus (e.g. ∼ 1000 K dust producing
photons with dimensionless energy ǫ0 = 5 × 10−7) due to
the variation of the Compton cross-section at high energies.
Blazars are bright and highly variable at all wave-
Constraints on the Ultra-High Energy Neutrino Flux from Gamma-Ray Bursts
from a Prototype Station of the Askaryan Radio Array
arXiv:1507.00100v1 [astro-ph.HE] 1 Jul 2015
P. Allisona , J. Auffenbergi , R. Bardb , J. J. Beattya , D. Z. Bessonc,d , C. Borae , C.-C. Chenf , P. Chenf , A. Connollya , J.
P. Daviesg , M. A. DuVernoisi , B. Foxh , P. W. Gorhamh , K. Hansonk , B. Hillh , K. D. Hoffmanb , E. Honga , L.-C. Huf ,
A. Ishiharal , A. Karlei , J. Kelleyi , I. Kravchenkoe , H. Landsmanj , A. Laundriei , C.-J. Lif , T. Liuf , M.-Y. Lui , R.
k
´
Maunub , K. Masel , T. Meuresk , C. Mikih , J. Namf , R. J. Nicholg , G. Nirj , A. OMurchadha
, C. G. Pfendnera , K.
n
b
h
i
m
e
b
Ratzlaff , M. Richman , B. Rotter , P. Sandstrom , D. Seckel , A. Shultz , M. Song , J. Stockhamc , M. Stockhamc ,
M. Sullivand , J. Touartb , H.-Y. Tuf , G. S. Varnerh , S. Yoshidal , R. Youngn , D.Guettao
a Dept.
of Physics and CCAPP, The Ohio State University, 191 W. Woodruff Ave., Columbus, OH 43210, USA
b Dept. of Physics, University of Maryland, College Park, MD 20742, USA
c Dept. of Physics and Astronomy, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, KS 66045, USA
d National Research Nuclear University - Moscow Engineering Physics Institute, 31 Kashirskaya Shosse, Moscow 115409, Russia
e Dept. of Physics and Astronomy, University of Nebraska-Lincoln, 855 N 16th Street, Lincoln, NE 68588, USA
f Dept. of Physics, Grad. Inst. of Astrophys.,& Leung Center for Cosmology and Particle Astrophysics, National Taiwan University, No. 1, Sec. 4,
Roosevelt Road, Taipei 10617, Taiwan (R.O.C.)
g Dept. of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
h Dept. of Physics and Astronomy, University of Hawaii-Manoa, 2505 Correa Rd., Honolulu, HI 96822, USA
i Dept. of Physics and Wisconsin IceCube Particle Astrophysics Center, University of Wisconsin-Madison, 222 W. Washington Ave, Madison, WI
53706, USA
j Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot, 76100, Israel
k Service de physique des particules e
´ l´ementaires, Universit´e Libre de Bruxelles, CP230, boulevard du Triomphe, 1050 Bruxelles, Belgium
l Dept. of Physics, Chiba University, 1-33, Yayoi-cho, Inage-ku, Chiba-shi, Chiba 263-8522, Japan
m Dept. of Physics and Astronomy, University of Delaware, 104 The Green, Newark, DE 19716, USA
n Instrumentation Design Laboratory, University of Kansas, 1251 Wescoe Drive, Lawrence, KS 66045, USA
o ORT Braude, Karmiel 21982, OAR-INAF, Italy
Abstract
We searched for ultra-high energy (UHE) neutrinos from Gamma-Ray Bursts (GRBs) with the Askaryan Radio Array
(ARA) Testbed station’s 2011-2012 data set. Among 589 GRBs monitored by the Gamma Ray Coordinate Network
(GCN) catalog from Jan. 2011 to Dec. 2012 over the entire sky, 57 GRBs were selected for analysis. These GRBs were
chosen because they occurred during a period of low anthropogenic background and high stability of the station and
fell within our geometric acceptance. We searched for UHE neutrinos from 57 GRBs and observed 0 events, which
is consistent with 0.11 expected background events. With this result, we set the limits on the UHE GRB neutrino
fluence and quasi-diffuse flux from 1016 to 1019 eV. This is the first limit on the UHE GRB neutrino quasi-diffuse flux
at energies above 1016 eV.
Keywords: Gamma-Ray Bursts, UHE neutrinos, radio Cherenkov
1. Introduction
GRBs are the most powerful explosions in the Universe, and include the highest redshift objects observed.
The widely accepted phenomenological interpretation
of these cosmological sources is the so-called fireball
(FB) model [1, 2]. In this model the energy carried
Email address: connolly@physics.osu.edu (A. Connolly)
Preprint submitted to Astroparticle Physics
by the hadrons in a relativistic expanding wind (fireball) may be dissipated through internal shocks between
different parts of plasma. These shocks reconvert a
substantial part of the kinetic energy to internal energy. This internal energy is then radiated as prompt
gamma rays by inverse-Compton and synchrotron radiation of shock-accelerated electrons. When the fireball
has swept enough material it collides with its surrounding medium giving rise to reverse and forward shocks.
The latter is believed to be responsible for so-called afJuly 2, 2015
arXiv:1507.00096v1 [astro-ph.CO] 1 Jul 2015
The Hubble Flow of Plateau Inflation
Dries Coone,a,b Diederik Roesta and Vincent Venninc
a Van
Swinderen Institute for Particle Physics and Gravity, University of Groningen, Nijenborgh 4, 9747 AG Groningen, (The Netherlands)
b Theoretische Natuurkunde, Vrije Universiteit Brussel and The International Solvay Institutes Pleinlaan 2, B-1050 Brussels, (Belgium)
c Institute of Cosmology & Gravitation, University of Portsmouth, Dennis Sciama Building,
Burnaby Road, Portsmouth, PO1 3FX, (United Kingdom)
E-mail: a.a.coone@rug.nl, d.roest@rug.nl, vincent.vennin@port.ac.uk
Abstract. In the absence of CMB precision measurements, a Taylor expansion has often been
invoked to parametrize the Hubble flow function during inflation. The standard “horizon
flow” procedure implicitly relies on this assumption. However, the recent Planck results
indicate a strong preference for plateau inflation, which suggests the use of Pad´e approximants
instead. We propose a novel method that provides analytic solutions of the flow equations
for a given parametrization of the Hubble function. This method is illustrated in the Taylor
and Pad´e cases, for low order expansions. We then present the results of a full numerical
treatment scanning larger order expansions, and compare these parametrizations in terms
of convergence, prior dependence, predictivity and compatibility with the data. Finally, we
highlight the implications for potential reconstruction methods.
Keywords: physics of the early universe, inflation
c ESO 2015
Astronomy & Astrophysics manuscript no. draft
July 2, 2015
Environmental dependence of polycyclic aromatic hydrocarbon
emission at z∼0.8
Investigation by observing the RX J0152.7-1357 with AKARI
Kazumi Murata1 , Yusei Koyama1,2 , Masayuki Tanaka2 , Hideo Matsuhara1 , and Tadayuki Kodama2
arXiv:1507.00094v1 [astro-ph.GA] 1 Jul 2015
1
2
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, Sagamihara, Kanagawa 229-8510, Japan
e-mail: murata@ir.isas.jaxa.jp
National Astronomical Observatory of Japan, Osawa 2-21-1, Mitaka, Tokyo 181-8588, Japan
Received April 1, 2015; accepted June 29, 2015
ABSTRACT
We study the environmental dependence of the strength of polycyclic aromatic hydrocarbon (PAH) emission by AKARI observations
of RX J0152.7-1357, a galaxy cluster at z=0.84. PAH emission reflects the physical conditions of galaxies and dominates 8 µm
luminosity (L8), which can directly be measured with the L15 band of AKARI. L8 to infrared luminosity (LIR) ratio is used as a
tracer of the PAH strength. Both photometric and spectroscopic redshifts are applied to identify the cluster members. The L15-banddetected galaxies tend to reside in the outskirt of the cluster and have optically green colour, R − z′ ∼ 1.2. We find no clear difference
of the L8/LIR behaviour of galaxies in field and cluster environment. The L8/LIR of cluster galaxies decreases with specific-starformation rate divided by that of main-sequence galaxies, and with LIR, consistent with the results for field galaxies. The relation
between L8/LIR and LIR is between those at z = 0 and z = 2 in the literature. Our data also shows that starburst galaxies, which have
lower L8/LIR than main-sequence, are located only in the outskirt of the cluster. All these findings extend previous studies, indicating
that environment affects only the fraction of galaxy types and does not affect the L8/LIR behaviour of star-forming galaxies.
Key words. – infrared: galaxies – galaxies: starburst – galaxies: clusters: individual (RX J0152.7-1357)
1. Introduction
Galaxy formation and evolution are highly affected by its environment. In the local Universe, galaxies in clusters tend to be
red, early type (Dressler 1980) and show lower star-forming activities. This star formation - density relation can also be seen
in the distant universe (Quadri et al. 2012), while some studies
suggest that the relation is reversed at z ≥ 1 (Elbaz et al. 2007;
Cooper et al. 2008). Although star formation - density relation
at higher redshift is still under debate, it is clear that the environment must play an important role in galaxy formation.
Recent studies suggest that the environment affects only
the fraction of galaxy types, and star-forming galaxies in different environment have similar properties (Patel et al. 2009a;
Muzzin et al. 2012; Koyama et al. 2013, 2014). Patel et al.
(2009a) suggest that lower star formation rate (SFR) of galaxies
in higher density may reflect the lower fraction of star forming
galaxies, rather than a decrease of SFRs. Koyama et al. (2013)
discusses that SFR - stellar mass (M∗ ) relation of Hα selected
galaxies is independent of the environment since z ∼ 2. On the
other hand, they also show that the dust attenuation and the median SFR of Hα emitters are higher in higher density environment at z = 0.4, which can be seen only with infrared data. Some
studies reported that dusty star-forming galaxies reside in groups
(Koyama et al. 2008, 2011; Tran et al. 2009). These studies imply that dust has important roles in understanding the environmental effects on galaxy formation.
Polycyclic aromatic hydrocarbon (PAH) is dust showing
prominent features at mid-infrared (MIR). They are excited by
UV light from young stars and emit the energy at 3.3, 6.2, 7.7,
8.6 and 11.3 µm (Draine & Li 2007). Among the PAH bands,
7.7 and 8.6 µm features dominate 8 µm luminosity (L8) even in
broad band filters, so that L8 is used as a PAH luminosity tracer.
Due to its energy sources, the PAH emission has been thought
to correlate with the infrared luminosity (LIR), which is correlated with SFR. Recent works suggest that high specific SFR
(sSFR) galaxies tend to show weak PAH emission compared
with main-sequence galaxies (Elbaz et al. 2011; Nordon et al.
2012; Murata et al. 2014). The main causes of the PAH weakness is thought to be destruction of PAHs or a lack of UV photon that excites the PAHs. These studies indicate that the PAH
emission traces the physical conditions of the interstellar matter rather than the SFR. On the other hand, the relation between infrared luminosity and PAH emission evolve with redshift (Takagi et al. 2010; Nordon et al. 2012). It suggests that the
physical conditions of galaxies with given infrared luminosity is
different at different redshift. Hence, if the cluster environment
affects the galaxy properties, the relation between infrared luminosity and PAH emission may be different.
However, the behaviour of the PAH emission in cluster environments has not yet been studied well. This is mainly due to
sparse filter sampling at 8-24 µm in the Spitzer space telescope
(Werner et al. 2004), where 7.7 µm PAH emission of galaxies
at z < 2 is redshifted. In contrast, Japanese AKARI satellite
(Murakami et al. 2007) has continuous wavelength coverage at
2-24 µm with nine photometric bands in the infrared camera
(IRC; Onaka et al.2007). It enables us to measure the PAH emission at z < 2.
The galaxy cluster RX 0152.7-1357 (here after RXJ0152)
at z∼0.84 is one of the suitable targets for investigating an
1
arXiv:1507.00079v1 [astro-ph.GA] 1 Jul 2015
The GALAH Survey and Galactic Archaeology in the next decade
Sarah L. Martell1
1 School
of Physics, University of New South Wales, Sydney, NSW 2052,
Australia; s.martell@ unsw.edu.au
Abstract. The field of Galactic Archaeology aims to understand the origins and evolution of the stellar populations in the Milky Way, as a way to understand galaxy formation and evolution in general. The GALAH (Galactic Archaeology with HERMES)
Survey is an ambitious Australian-led project to explore the Galactic history of star
formation, chemical evolution, minor mergers and stellar migration. GALAH is using
the HERMES spectrograph, a novel, highly multiplexed, four-channel high-resolution
optical spectrograph, to collect high-quality R ∼ 28, 000 spectra for one million stars in
the Milky Way. From these data we will determine stellar parameters, radial velocities
and abundances for up to 29 elements per star, and carry out a thorough chemical tagging study of the nearby Galaxy. There are clear complementarities between GALAH
and other ongoing and planned Galactic Archaeology surveys, and also with ancillary
stellar data collected by of major cosmological surveys. Combined, these data sets will
provide a revolutionary view of the structure and history of the Milky Way.
1.
Galactic Archaeology
Galactic Archaeology is a highly detailed approach to near-field cosmology, using the
structure and history of the Milky Way to build a systematic understanding of general
processes in galaxy evolution. One of the main goals of Galactic Archaeology is to
identify the remnants of ancient building blocks of the Galaxy that are now dispersed.
Chemical tagging, which follows the reasoning that stars that form together at the same
time and place should have matching abundance patterns since they are forming from
the same gas cloud, was first suggested by Freeman & Bland-Hawthorn (2002) as a
way to search the present-day Galaxy for stars that originated together in these building
blocks.
Large-scale Galactic Archaeology surveys like GALAH (De Silva et al. 2015),
Gaia-ESO (Gilmore et al. 2012) and APOGEE (Apache Point Observatory Galaxy
Evolution Experiment, Holtzman et al. 2015) will enable powerful chemical tagging
by acquiring large, homogeneous, high-quality sets of spectroscopic data, and using
those to derive high-precision stellar parameters and abundances. At present, the precision of the parameters and abundances is still approaching its intended level, which
limits the chemical tagging that can be carried out with the data. This is in part a result
of the diversity of stars in each survey: stars with very different effective temperatures,
surface gravities and compositions have very different atmospheres and emergent spectra, and require different and specialized approaches. Such a large data volume requires
automated reduction and analysis, and in such a complex and degenerate problem, the
1
Accepted by the ApJL
arXiv:1507.00078v1 [astro-ph.SR] 1 Jul 2015
EVIDENCE OF THE SOLAR EUV HOT CHANNEL AS A MAGNETIC
FLUX ROPE FROM REMOTE-SENSING AND IN-SITU OBSERVATIONS
H.Q. SONG1 , Y. CHEN1 , J. ZHANG2 , X. CHENG3 , B. Wang1 , Q. HU4 , G. LI4 , AND Y.M.
WANG5
1 Shandong Provincial Key Laboratory of Optical Astronomy and Solar-Terrestrial Environment,
and Institute of Space Sciences, Shandong University, Weihai, Shandong 264209, China
hqsong@sdu.edu.cn
2 School of Physics, Astronomy and Computational Sciences, George Mason University, Fairfax,
VA 22030, USA
3 School of Astronomy and Space Science, Nanjing University, Nanjing, Jiangsu 210093, China
4 Department of Space Science and CSPAR, University of Alabama in Huntsville, Huntsville, AL
35899, USA
5 Key Laboratory of Geospace Environment, University of Science and Technology of China,
Chinese Academy of Sciences (CAS), Hefei, Anhui 230026, China
ABSTRACT
Hot channels (HCs), high temperature erupting structures in the lower corona of the
Sun, have been proposed as a proxy of magnetic flux ropes (MFRs) since their initial
discovery. However, it is difficult to make definitive proof given the fact that there is
no direct measurement of magnetic field in the corona. An alternative way is to use
the magnetic field measurement in the solar wind from in-situ instruments. On 2012
July 12, an HC was observed prior to and during a coronal mass ejection (CME) by
the AIA high-temperature images. The HC is invisible in the EUVI low-temperature
images, which only show the cooler leading front (LF). However, both the LF and an
ejecta can be observed in the coronagraphic images. These are consistent with the high
temperature and high density of the HC and support that the ejecta is the erupted HC.
In the meanwhile, the associated CME shock was identified ahead of the ejecta and
the sheath through the COR2 images, and the corresponding ICME was detected by
ACE, showing the shock, sheath and magnetic cloud (MC) sequentially, which agrees
with the coronagraphic observations. Further, the MC contained a low-ionization-state
center and a high-ionization-state shell, consistent with the pre-existing HC observation
and its growth through magnetic reconnection. All of these observations support that
Accepted to ApJS 2015 June 27
arXiv:1507.00057v1 [astro-ph.SR] 30 Jun 2015
The Brown Dwarf Kinematics Project (BDKP). IV. Radial Velocities of 85
Late-M and L dwarfs with MagE
Adam J. Burgasser1 , Sarah E. Logsdon2 , Jonathan Gagn´e3 , John J. Bochanski4 , Jaqueline K.
Faherty5,6,7 , Andrew A. West8 , Eric E. Mamajek9 , Sarah J. Schmidt10 , and Kelle L. Cruz6,11,12
ABSTRACT
Radial velocity measurements are presented for 85 late M- and L-type very low
mass stars and brown dwarfs obtained with the Magellan Echellette (MagE) spectrograph. Targets primarily have distances within 20 pc of the Sun, with more distant
sources selected for their unusual spectral energy distributions. We achieved precisions
of 2–3,km s−1 , and combined these with astrometric and spectrophotometric data to
calculate U V W velocities. Most are members of the thin disk of the Galaxy, and velocity dispersions indicate a mean age of 5.2±0.2 Gyr for sources within 20 pc. We
find signficantly different kinematic ages between late-M dwarfs (4.0±0.2 Gyr) and L
dwarfs (6.5±0.4 Gyr) in our sample that are contrary to predictions from prior simulations. This difference appears to be driven by a dispersed population of unusually
1
Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA, 92093, USA;
aburgasser@ucsd.edu
2
Department of Physics and Astronomy, UCLA, 430 Portola Plaza, Box 951547, Los Angeles, CA 90095-1547,
USA
3
Institute for Research on Exoplanets (iREx), Universit´e de Montr´eal, D´epartement de Physique, C.P. 6128 Succ.
Centre-ville, Montr´eal, QC H3C 3J7, Canada
4
Rider University, 2083 Lawrenceville Road, Lawrenceville, NJ 08648, USA
5
Department of Terrestrial Magnetism, Carnegie Institution of Washington, Washington, DC 20015, USA
6
Department of Astrophysics, American Museum of Natural History, Central Park West at 79th Street, New York,
NY 10034, USA
7
Hubble Postdoctoral Fellow
8
Department of Astronomy, Boston University, 725 Commonwealth Ave Boston, MA 02215 USA
9
Department of Physics & Astronomy, University of Rochester, Rochester, NY 14627, USA
10
Department of Astronomy, Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA
11
Department of Physics and Astronomy, Hunter College, 695 Park Avenue, New York, NY10065, USA
12
Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016,
USA
Mon. Not. R. Astron. Soc. 000, 1–?? (2013)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00055v1 [astro-ph.SR] 30 Jun 2015
Evidence for a [WR] or WEL–type binary nucleus in the
bipolar planetary nebula Vy 1–2
S.
Akras1,2,3 ⋆ , P. Boumis3, J. Meaburn4, J. Alikakos3, J. A. L´opez2, D. R. Gon¸calves1
1
Observat´
orio do Valongo, Universidade Federal do Rio de Janeiro, Ladeira Pedro Antonio 43, 20080-090, Rio de Janeiro, Brazil
de Astronomia, Universidad National Autonoma de Mexico, Ensenada 22800, Baja California, Mexico.
3 IAASARS, National Observatory of Athens, I. Metaxa & V. Pavlou, Penteli, GR–15236, Athens, Greece.
4 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, M13 9PL, UK
2 Instituto
Received **insert**; Accepted **insert**
ABSTRACT
We present high–dispersion spectroscopic data of the compact planetary nebula
Vy 1–2, where high expansion velocities up to 100 km s−1 are found in the Hα, [N ii]
and [O iii] emission lines. HST images reveal a bipolar structure. Vy 1–2 displays a
bright ring–like structure with a size of 2.4′′ ×3.2′′ and two faint bipolar lobes in the
west–east direction. A faint pair of knots is also found, located almost symmetrically on
opposite sides of the nebula at PA=305◦. Furthermore, deep low–dispersion spectra are
also presented and several emission lines are detected for the first time in this nebula,
such as the doublet [Cl iii] 5517, 5537 ˚
A, [K iv] 6101 ˚
A, C ii 6461 ˚
A, the doublet C iv
5801, 5812 ˚
A. By comparison with the solar abundances, we find enhanced N, depleted
C and solar O. The central star must have experienced the hot bottom burning (CN–
cycle) during the 2nd dredge–up phase, implying a progenitor star of M>3 M⊙ . The
very low C/O and N/O abundance ratios suggest a likely post–common envelope close
binary system.
A simple spherically symmetric geometry with either a blackbody or a H–deficient
stellar atmosphere model is not able to reproduce the ionisation structure of Vy 1–2.
The effective temperature and luminosity of its central star indicate a young nebula
located at a distance of ∼9.7 kpc with an age of ∼3500 years. The detection of stellar
emission lines, C ii 6461 ˚
A, the doublet C ivλλ 5801, 5812 and O iii 5592 ˚
A, emitted
from a H–deficient star, indicates the presence of a late–type Wolf–Rayet or a WEL
type central star.
Key words: ISM: jets and outflows – ISM: abundances – planetary nebula: individual: Vy 1–2 – stars: Wolf–Rayet
1
INTRODUCTION
Planetary nebulae (PNe) are the consequent result of the interaction of stellar winds ejected from low–to–intermediate
mass stars (1-8 M⊙ ). A slow and dense wind (10−4 M⊙ ),
emerges from an Asymptotic Giant Branch (AGB) star, is
swept up into a dense spherically symmetric shell by a second, faster and tenuous wind expelled during the post-AGB
phase (interacting stellar wind model; Kwok et al. 1978).
This model can reasonably explain the formation and evolution of spherical PNe. However, most post–AGB stars and
PNe show highly asymmetric circumstellar envelopes and
complex morphologies (Manchado et al. 2006, Meixner et
⋆
e-mail:akras@astro.ufrj.br
c 2013 RAS
al. 1999, Parker et al. 2006, Boumis et al. 2003, 2006, Miszalski et al. 2008, Lagadec et al. 2011, Sahai 2011, Sabin et
al. 2014).
The deviation of PNe from spherically symmetric structures has its origin in an enhanced equatorial–to–polar density gradient in the AGB circumstellar envelopes resulting
in the formation of aspherical and more complex nebulae.
This equatorial density enhancement is interpreted as the
result of common envelope evolution of close binary systems
(Soker & Livio 1994, Nordhaus & Blackman 2006, De Marco
2009). Although, most of the post–AGB and PNe show complex morphologies, only a small number of them are known
to contain a binary system (Miszalski et al. 2009, De Marco
2009). Recent work by Jones et al. (2014) confirmed the
binary nature of Hen 2–11’s nucleus, providing further evi-
Mon. Not. R. Astron. Soc. 000, 1–17 ()
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00051v1 [astro-ph.SR] 30 Jun 2015
A Consistent Spectral Model of WR 136 and its
Associated Bubble NGC 6888
J.
Reyes-P´erez1? , C. Morisset1 , M. Pe˜
na1 and A. Mesa-Delgado2 .
1
Instituto de Astronom´ıa, Universidad Nacional Aut´
onoma de M´
exico
de Astrof´ısica, Facultad de F´ısica, Pontificia Universidad Cat´
olica de Chile
2 Instituto
ABSTRACT
We analyse whether a stellar atmosphere model computed with the code CMFGEN
provides an optimal description of the stellar observations of WR 136 and simultaneously reproduces the nebular observations of NGC 6888, such as the ionization degree,
which is modelled with the pyCloudy code. All the observational material available (far
and near UV and optical spectra) were used to constrain such models. We found that
even when the stellar luminosity and the mass-loss rate were well constrained, the
stellar temperature T∗ at τ = 20, can be in a range between 70 000 and 110 000 K.
When using the nebula as an additional restriction we found that the stellar models
with T∗ ∼ 70 000 K represent the best solution for both, the star and the nebula.
Results from the photoionization model show that if we consider a chemically homogeneous nebula, the observed N+ /O+ ratios found in different nebular zones can be
reproduced, therefore it is not necessary to assume a chemical inhomogeneous nebula.
Our work shows the importance of calculating coherent models including stellar and
nebular constraints. This allowed us to determine, in a consistent way, all the physical parameters of both the star and its associated nebula. The chemical abundances
derived are 12 + log(N/H) = 9.95, 12 + log(C/H) = 7.84 and 12 + log(O/H) = 8.76
for the star and 12 + log(N/H) = 8.40, 12 + log(C/H) = 8.86 and 12 + log(O/H) =
8.20. Thus the star and the nebula are largely N- and C- enriched and O-depleted.
Key words: circumstellar matter – stars: Wolf-Rayet stars.
1
INTRODUCTION
NGC 6888 is an emission nebula associated with the WolfRayet star WR 136 with spectral type WN6(h), and it is a
typical example of a ring nebula resulting from the continuous interaction between a strong stellar wind and the interstellar medium. Projected on the sky, this nebula seems to
be a filamentary and clumpy gaseous ellipse emitting mainly
in Hα and optical [N II] lines, with an axis-size of 12 × 18 arcmin (Chu, Treffers, & Kwitter 1983), and it is surrounded
by a homogeneous sphere which emits mainly in [O III] 5007
˚
A (Gruendl et al. 2000; Moore, Hester & Scowen 2000). The
combination of its filamentary morphology and its nitrogen enhanced abundance suggests that the ejections, coming
from a previous evolutionary stage of the star WR 136, have
contributed in an important way to the formation of the ring
nebula that we see today (Esteban & V´ılchez 1992, hereafter
EV92; Fern´
andez-Mart´ın et al. 2012, hereafter FM12; MesaDelgado et al. 2014, hereafter M-D14).
?
E-mail: jperez@astro.unam.mx (JRP)
Both the nebula NGC 6888 and its central star WR 136
have been already studied in several wavelength ranges from
X-rays to radio emission. For example, the stellar parameters, such as the stellar temperature, mass-loss rate and terminal velocity, have been calculated by Crowther & Smith
(1996) and Hamann, Grafener & Liermann (2006, hereafter
HGL06) by using stellar atmosphere codes and comparing
the results with observed ultraviolet, optical and infrared
stellar spectra; while the wind velocity law and clumpiness were studied by Ignace, Quigley & Cassinelli (2003).
Skinner et al. (2010) showed that WR 136 emits in X-ray
energies. In particular, these authors found soft emission
(kT ∼ 0.56 keV) coming from outside the wind and hard
emission (kT ∼ 2.64 keV) coming from the gas inside the
wind.
The complete hot bubble in the interior of NGC 6888
was observed with the Space Telescope Suzaku (Zhekov &
Park 2011), confirming the existence of nitrogen enriched
X-ray emitting plasma which consists of two thermal components, one emiting soft X-rays and the other one, hard
X-rays. Based on this, the authors proposed a new X-ray
**Volume Title**
ASP Conference Series, Vol. **Volume Number**
**Author**
c **Copyright Year** Astronomical Society of the Pacific
arXiv:1507.00034v1 [astro-ph.IM] 30 Jun 2015
The MUSE Data Reduction Pipeline:
Status after Preliminary Acceptance Europe
Peter M. Weilbacher1 , Ole Streicher1 , Tanya Urrutia1 , Arlette
P´econtal-Rousset2 , Aur´elien Jarno2 , Roland Bacon2
1 Leibniz-Institut
f¨ur Astrophysik Potsdam (AIP), An der Sternwarte 16,
D-14482 Potsdam, Germany
2 Universit´
e
de Lyon, Lyon, F-69003, France; Universit´e Lyon 1, Observatoire
de Lyon, 9 avenue Charles Andr´e, Saint-Genis Laval, F-69230, France; CNRS,
UMR 5574, Centre de Recherche Astrophysique de Lyon; Ecole Normale
Sup´erieure de Lyon, Lyon, F-69007, France
Abstract. MUSE, a giant integral field spectrograph, is about to become the newest
facility instrument at the VLT. It will see first light in February 2014. Here, we summarize the properties of the instrument as built and outline functionality of the data reduction system, that transforms the raw data that gets recorded separately in 24 IFUs by 4k
CCDs, into a fully calibrated, scientifically usable data cube. We then describe recent
work regarding geometrical calibration of the instrument and testing of the processing
pipeline, before concluding with results of the Preliminary Acceptance in Europe and
an outlook to the on-sky commissioning.
1.
The instrument
MUSE (Bacon et al. 2012) is a wide-field integral field spectrograph built for the ESO
Very Large Telescope. It passed its Preliminary Acceptance in Europe (PAE) on Sept. 10th,
2013, and is currently being shipped to Chile for installation at the telescope. MUSE
consists of 24 image slicer integral field units (IFUs) in as many subsections of the field
of view. Spatial and spectral sampling of 0.′′ 2 per spaxel and 1.25 Å per pixel over a
field of ∼ 1′ × 1′ coupled with a spectral resolution between 1600 and 3500 ensure
the usefulness of MUSE for a large variety of scientific observations. For the study of
the early universe the long wavelength range of 480 to 930 nm allows the detection of
Lyman-α emission in the redshift range 3 . z . 6.6. A high throughput of &35% with
a peak of ∼55% at 7000 Å in conjunction with high stability will enable users to carry
out the deepest spectroscopic observations of the universe. In the future, MUSE will
be coupled with the adaptive optics module GALACSI (Str¨obele et al. 2012), and reach
even fainter limits with ground-layer correction. A narrow-field mode will use laser
tomography for high-order correction and deliver significant strehl ratios in the red part
of the spectrum over a field of 7.′′ 5 × 7.′′ 5.
1
c
ESO
2015
Astronomy & Astrophysics manuscript no. PaperIV_arXiv
July 2, 2015
Gaia FGK benchmark stars: abundances of α and iron-peak
elements ⋆ ⋆⋆
arXiv:1507.00027v1 [astro-ph.SR] 30 Jun 2015
P. Jofré1 , U. Heiter2 , C. Soubiran3 , S. Blanco-Cuaresma3, 4 , T. Masseron1 , T. Nordlander2 , L. Chemin3 , C. C. Worley1 ,
S. Van Eck5 , A. Hourihane1 , G. Gilmore1 , V. Adibekyan6 , M. Bergemann1, 7 , T. Cantat-Gaudin8 , E. Delgado-Mena6 , J.
I. González Hernández15, 16 , G. Guiglion9 , C. Lardo10 , P. de Laverny9 , K. Lind2 , L. Magrini11 , S. Mikolaitis9, 12 , D.
Montes15 , E. Pancino13, 14 , A. Recio-Blanco9 , R. Sordo8 , S. Sousa6 , H. M. Tabernero15 , and A. Vallenari8
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, United Kingdom
Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden
Univ. Bordeaux, CNRS, Laboratoire d’Astrophysique de Bordeaux (UMR 5804), F-33270, Floirac, France
Observatoire de Genéve, Université de Genéve, CH-1290 Versoix, Switzerland
Institut d’Astronomie et d’Astrophysique, U. Libre de Bruxelles, CP 226, Boulevard du Triomphe, 1050 Bruxelles, Belgium
Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal
Max-Planck Institute for Astronomy, 69117, Heidelberg, Germany
INAF, Osservatorio di Padova, Università di Padova, Vicolo Osservatorio 5, Padova 35122, Italy
Laboratoire Lagrange (UMR7293), Univ. Nice Sophia Antipolis, CNRS, Observatoire de la Côte d’Azur, 06304 Nice, France
Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK
INAF/Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5 50125 Firenze, Italy
Institute of Theoretical Physics and Astronomy, Vilnius University, A. Goštauto 12, LT-01108 Vilnius, Lithuania
INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 40127 Bologna, Italy
ASI Science Data Center, via del Politecnico s/n, 00133 Roma, Italy
Dpto. Astrofísica, Facultad de CC. Físicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
Instituto de Astrofí sica de Canarias, E-38205 La Laguna, Tenerife, Spain
ABSTRACT
Context. In the current era of large spectroscopic surveys of the Milky Way, reference stars for calibrating astrophysical parameters
and chemical abundances are of paramount importance.
Aims. We determine elemental abundances of Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Co and Ni for our predefined set of Gaia FGK benchmark
stars.
Methods. By analysing high-resolution and high-signal to noise spectra taken from several archive datasets, we combined results
of eight different methods to determine abundances on a line-by-line basis. We perform a detailed homogeneous analysis of the
systematic uncertainties, such as differential versus absolute abundance analysis, as well as we assess errors due to NLTE and the
stellar parameters in our final abundances.
Results. Our results are provided by listing final abundances and the different sources of uncertainties, as well as line-by-line and
method-by-method abundances.
Conclusions. The Gaia FGK benchmark stars atmospheric parameters are already being widely used for calibration of several
pipelines applied to different surveys. With the added reference abundances of 10 elements this set is very suitable to calibrate the
chemical abundances obtained by these pipelines.
1. Introduction
Much of our understanding on the structure and evolution of the
Milky Way comes nowadays from the analysis of large stellar
spectroscopic surveys. After the revolution of the way to pursue Galactic science with the low-resolution spectra from SDSS
data (see Ivezi´c et al. 2012, for a review), new surveys are ongoing. These have much higher resolution than SDSS allowing
not only to determine the stellar parameters of the stars more
Send offprint requests to:
P. Jofré, e-mail: pjofre@ast.cam.ac.uk
⋆
Based on NARVAL and HARPS data obtained within the Gaia
DPAC (Data Processing and Analysis Consortium) and coordinated by
the GBOG (Ground-Based Observations for Gaia) working group and
on data retrieved from the ESO-ADP database.
⋆⋆
Tables XX are only available in electronic form at the CDS
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/
precisely but also the chemical abundances of several individual elements. Examples of such projects are the Gaia-ESO Survey (Gilmore et al. 2012; Randich et al. 2013, hereafter GES),
RAVE (Steinmetz et al. 2006), APOGEE (Allende Prieto et al.
2008), GALAH (De Silva et al. 2015) and the future billion of
stars from the Radial Velocity Spectrograph (RVS) from Gaia.
Furthermore, several groups have collected over the years large
samples of stars, creating “independent" surveys for the same
purpose of unraveling the structure and chemical enrichment history of our Galaxy (e.g. Fuhrmann 2011; Adibekyan et al. 2012;
Ramírez et al. 2013; Bensby et al. 2014, and references therein).
To parametrise these data properly in an automatic way, and
furthermore, to link the data between the different surveys in a
consistent way, good standard calibrators are needed. To this
aim, we have defined a sample, the Gaia FGK benchmark stars
(hereafter GBS), which includes 34 FGK stars of a wide range of
metallicities and gravities. These stars should be representative
Article number, page 1 of 65
Draft version July 2, 2015
Preprint typeset using LATEX style emulateapj v. 5/14/03
THE SUBARU CORONAGRAPHIC EXTREME ADAPTIVE OPTICS SYSTEM:
ENABLING HIGH-CONTRAST IMAGING ON SOLAR-SYSTEM SCALES
N. Jovanovic1 , F. Martinache2 , O. Guyon1,3,4 , C. Clergeon1 , G. Singh1,5 , T. Kudo1 ,V. Garrel6 ,
K. Newman4,7 , D. Doughty1,4 , J. Lozi1 , J. Males3,8 , Y. Minowa1 , Y. Hayano1 , N. Takato1 , J. Morino9 ,
J. Kuhn10 , E. Serabyn10 , B. Norris11 , P. Tuthill11 , G. Schworer5,11 , P. Stewart11 , L. Close3 , E. Huby5,12 ,
G. Perrin5 , S. Lacour5 , L. Gauchet5 , S. Vievard5 , N. Murakami13 , F. Oshiyama13 , N, Baba13 , T. Matsuo14 ,
J. Nishikawa9 , M. Tamura9,15 , O. Lai1,6 , F. Marchis16 , G. Duchene17,18 , T. Kotani9 , and J. Woillez19
arXiv:1507.00017v1 [astro-ph.IM] 30 Jun 2015
Draft version July 2, 2015
ABSTRACT
The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument is a multipurpose highcontrast imaging platform designed for the discovery and detailed characterization of exoplanetary
systems and serves as a testbed for high-contrast imaging technologies for ELTs. It is a multi-band
instrument which makes use of light from 600 to 2500 nm allowing for coronagraphic direct exoplanet
imaging of the inner 3 λ/D from the stellar host. Wavefront sensing and control are key to the
operation of SCExAO. A partial correction of low-order modes is provided by Subaru’s facility adaptive optics system with the final correction, including high-order modes, implemented downstream
by a combination of a visible pyramid wavefront sensor and a 2000-element deformable mirror. The
well corrected NIR (y-K bands) wavefronts can then be injected into any of the available coronagraphs, including but not limited to the phase induced amplitude apodization and the vector vortex
coronagraphs, both of which offer an inner working angle as low as 1 λ/D. Non-common path, loworder aberrations are sensed with a coronagraphic low-order wavefront sensor in the infrared (IR). Low
noise, high frame rate, NIR detectors allow for active speckle nulling and coherent differential imaging,
while the HAWAII 2RG detector in the HiCIAO imager and/or the CHARIS integral field spectrograph (from mid 2016) can take deeper exposures and/or perform angular, spectral and polarimetric
differential imaging. Science in the visible is provided by two interferometric modules: VAMPIRES
and FIRST, which enable sub-diffraction limited imaging in the visible region with polarimetric and
spectroscopic capabilities respectively. We describe the instrument in detail and present preliminary
results both on-sky and in the laboratory.
Subject headings: Astronomical Instrumentation, Extrasolar Planets
1. INTRODUCTION
1
National Astronomical Observatory of Japan, Subaru Telescope, 650 North A’Ohoku Place, Hilo, HI, 96720, U.S.A.
2 Observatoire de la Cote d’Azur, Boulevard de l’Observatoire,
Nice, 06304, France
3 Steward Observatory, University of Arizona, Tucson, AZ,
85721, U.S.A.
4 College of Optical Sciences, University of Arizona, Tucson, AZ
85721, USA
5 LESIA, Observatoire de Paris, Meudon, 5 Place Jules Janssen,
92195, France.
6 Gemini Observatory, c/o AURA, Casilla 603, La Serena, Chile
7 NASA Ames Research Center, Moffett Field, CA 94035, USA
8 NASA Sagan Fellow
9 National Astronomical Observatory of Japan, 2-21-1 Osawa,
Mitaka, Japan
10 Jet Propulsion Laboratory, California Institute of Technology,
4800 Oak Grove Dr, Pasadena, CA 91109, USA
11 Sydney Institute for Astronomy (SIfA), Institute for Photonics and Optical Science (IPOS), School of Physics,
University of Sydney, NSW 2006, Australia
12 Dpartement d?Astrophysique, Gophysique et Ocanographie,
Universit de Lige, 17 Alle du Six Aot, 4000 Lige, Belgium
13 Division of Applied Physics, Faculty of Engineering, Hokkaido
University, Kita-13, Nishi-8, Kita-ku, Sapporo, Hokkaido 060-8628,
Japan
14 Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-ku, Kyoto 606-8502 Japan
15 Department of Astronomy, University of Tokyo, 7-3-1 Hongo,
Bunkyo, Tokyo 113-0033, Japan
16 Carl Sagan Center at the SETI Institute, Mountain View, CA
94043, USA
17 Astronomy Department, University of California, Berkeley,
CA 94720-3411, USA
18 University Grenoble Alpes & CNRS, Institut de Planetologie
The field of high-contrast imaging is advancing at
a great rate with several extreme adaptive optics systems having come online in 2014, including the Gemini
Planet Imager (GPI) (Macintosh 2014), the SpectroPolarimetric High-contrast Exoplanet REsearch instrument (SPHERE) (Beuzit et al. 2008), and the focus of
this work, the Subaru Coronagraphic Extreme Adaptive
Optics (SCExAO) system which join the already running
P1640 (Dekany et al. 2013). These systems all share a
similar underlying architecture: they employ a high order
wavefront sensor (WFS) and a deformable mirror (DM)
to correct for atmospheric perturbations enabling high
Strehl ratios in the near-infrared (NIR) (> 90%), while a
coronagraph is used to suppress on-axis starlight downstream. The primary motivation for such instrumentation is the direct detection of planetary mass companions
at contrasts of 10−5 –10−6 with respect to the host star,
at small angular separations (down to 1–5 λ/D) from the
host star.
The era of exoplanetary detection has resulted in
∼ 1500 planets so far confirmed (Han et al. 2014).
The majority of these were detected via the transit
et d’Astrophysique de Grenoble (IPAG), Grenoble F-3800, France
19 European Southern Observatory (ESO), Karl-SchwarzschildStr. 2, Garching 85748, Germany
Electronic address: jovanovic.nem@gmail.com
Accepted for publication in The Astronomical Journal
arXiv:1507.00016v1 [astro-ph.SR] 30 Jun 2015
Star Formation Across the W3 Complex
Carlos G. Rom´an-Z´
un
˜ iga1, Jason E. Ybarra1 , Guillermo D. Megias2 , Mauricio Tapia1 ,
Elizabeth A. Lada3 and Jo˜ao F. Alves4
ABSTRACT
We present a multi-wavelength analysis of the history of star formation in the
W3 complex. Using deep, near-infrared ground-based images, combined with images obtained with Spitzer and Chandra observatories, we identified and classified young embedded sources. We identified the principal clusters in the complex,
and determined their structure and extension. We constructed extinction-limited
samples for five principal clusters, and constructed K-band luminosity functions
(KLF) that we compare with those of artificial clusters with varying ages. This
analysis provided mean ages and possible age spreads for the clusters. We found
that IC 1795, the centermost cluster of the complex, still hosts a large fraction
of young sources with circumstellar disks. This indicates that star formation was
active in IC 1795 as recently as 2 Myr ago, simultaneous to the star forming activity in the flanking embedded clusters, W3-Main and W3(OH). A comparison
with carbon monoxide emission maps indicates strong velocity gradients in the
gas clumps hosting W3-Main and W3(OH) and show small receding clumps of gas
at IC 1795, suggestive of rapid gas removal (faster than the T Tauri timescale) in
the cluster forming regions. We discuss one possible scenario for the progression
of cluster formation in the W3 complex. We propose that early processes of gas
collapse in the main structure of the complex could have defined the progression of cluster formation across the complex with relatively small age differences
from one group to another. However, triggering effects could act as catalysts
for enhanced efficiency of formation at a local level, in agreement with previous
studies.
1
Instituto de Astronom´ıa, Universidad Nacional Aut´onoma de M´exico, Unidad Acad´emica en Ensenada,
Km 103 Carr. Tijuana–Ensenada, Ensenada 22860, Mexico
2
Facultad de F´ısica. Universidad de Sevilla. Dpto. F´ısica At´
omica, Molecular y Nuclear, Sevilla 41080,
Spain
3
Astronomy Department, University of Florida, 211 Bryant Space Sciences Center, FL 32611, USA
4
Institute of Astronomy, University of Vienna, T¨
urkenschanzstr. 17, 1180 Vienna, Austria
Mon. Not. R. Astron. Soc. 000, 1–?? (2015)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00012v1 [astro-ph.GA] 30 Jun 2015
SImulator of GAlaxy Millimeter/submillimeter Emission
(S´
IGAME?): CO emission from massive z = 2 main sequence
galaxies
Karen P. Olsen1 †, Thomas R. Greve,2 Christian Brinch,3,4 , Jesper Sommer-Larsen,1,5,6
Jesper
Rasmussen,1,7 Sune Toft,1 Andrew Zirm1
1
Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark
Physics and Astronomy, University College London
3 Centre for Star and Planet Formation (Starplan) and Niels Bohr Institute, University of Copenhagen,
Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark
4 DeIC, Technical University of Denmark, Building 309, DK-2800 Kgs. Lyngby, Denmark
5 Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany
6 Marie Kruses Skole, Stavnsholtvej 29-31, DK-3520 Farum, Denmark
7 Department of Physics, Technical University of Denmark, Building 309, DK-2800 Kgs. Lyngby, Denmark
2 Dept.
ABSTRACT
We present S´
IGAME (SImulator of GAlaxy Molecular Emission), a new numerical code
designed to simulate the 12 CO rotational line emission spectrum of galaxies. Using subgrid physics recipes to post-process the outputs of smoothed particle hydrodynamics
(SPH) simulations, a molecular gas phase is condensed out of the initial hot and partly
ionised SPH gas and distributed in Giant Molecular Cloud (GMCs). The GMCs are
subjected to far-UV radiation fields and cosmic ray ionisation rates which scale with
the local star formation rate volume density, thereby ensuring that the thermal state
of the gas is directly coupled to the in situ star formation conditions. Level populations
as well as line radiative transport of the CO rotational lines are solved for with the 3-D
radiative transfer code LIME. We have applied S´
IGAME to cosmological SPH simulations
of three disk galaxies at z = 2 with stellar masses in the range ∼ 0.5 − 2 × 1011 M
and star formation rates ∼ 40 − 140 M yr−1 , for which we predict a low-excitation
gas with CO intensity peaks at the CO J = 3 − 2 transition and total CO(3 − 2)
luminosities within the range of observations of corresponding star-forming galaxies
at z ∼ 1 − 2.5. Global αCO factors range from 1.4 to 1.6 M pc−2 (K km s−1 )−1 , i.e.
about a third of the Galactic value. On resolved scales, the model galaxies display an
increase in CO(J − (J − 1))/CO(1 − 0) brightness temperature line ratios at J > 3,
but a decrease in αCO towards the central regions, in agreement with observations
of nearby galaxies. Adopting a steeper GMC radial density profile or a more shallow
mass spectrum leads to increased αCO factors, though still below the Galactic value.
The inclusion of high pressure (Pext /kB > 104 K cm−3 ) environments descreases line
ratios at high-J.
Key words: galaxies: high-redshift – galaxies: ISM – ISM: clouds – ISM: lines and
bands – methods: numerical – radiative transfer
1
INTRODUCTION
Inferring the physical properties of the molecular gas in
galaxies is an important prerequisite for understanding their
star formation. The rotational lines of the carbon monoxide (CO) molecule has proven to be a valuable tracer of
?
S´
IGAME means ‘follow me’ in Spanish.
† E-mail: karen@dark-cosmology.dk
molecular gas properties, its rotational levels probing a wide
range of densities and easily excited by collisions with H2
molecules in the interstellar medium (ISM). The CO spectral line energy distribution (SLED) has therefore been used
extensively to probe the molecular gas amount of galaxies
at high and low redshift (see review by Carilli et al. 2013).
With facilities such as the Atacama Large Millimeter Array (ALMA) and the Jansky Very Large Array (JVLA) it is
now feasible to observe CO out to z >
∼ 6, thus paving the way
D RAFT VERSION J ULY 2, 2015
Preprint typeset using LATEX style emulateapj v. 5/2/11
SOFIA OBSERVATIONS OF SN 2010JL: ANOTHER NON-DETECTION OF THE 9.7 µM SILICATE DUST FEATURE
B RIAN J. W ILLIAMS , 1 O RI D. F OX , 2
arXiv:1507.00008v1 [astro-ph.HE] 30 Jun 2015
Draft version July 2, 2015
ABSTRACT
We present photometric observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA)
at 11.1 µm of the Type IIn supernova (SN IIn) 2010jl. The SN is undetected by SOFIA, but the upper limits
obtained, combined with new and archival detections from Spitzer at 3.6 & 4.5 µm allow us to characterize
the composition of the dust present. Dust in other Type IIn SNe has been shown in previous works to reside
in a circumstellar shell of material ejected by the progenitor system in the few millenia prior to explosion.
Our model fits show that the dust in the system shows no evidence for the strong, ubiquitous 9.7 µm feature
from silicate dust, suggesting the presence of carbonaceous grains. The observations are best fit with 0.010.05 M⊙ of carbonaceous dust radiating at a temperature of ∼ 550 − 620 K. The dust composition may reveal
clues concerning the nature of the progenitor system, which remains ambiguous for this subclass. Most of the
single star progenitor systems proposed for SNe IIn, such as luminous blue variables, red supergiants, yellow
hypergiants, and B[e] stars, all clearly show silicate dust in their pre-SN outflows. However, this post-SN result
is consistent with the small sample of SNe IIn with mid-IR observations, none of which show signs of emission
from silicate dust in their IR spectra.
Subject headings: dust, extinction — supernovae: individual(SN 2010jl) — supernovae: general
1. INTRODUCTION
Type IIn supernovae (SNe IIn) are named for the narrow hydrogen emission lines produced in a dense, slowlymoving, pre-existing circumstellar medium (CSM). Typically
attributed to the progenitor winds, the CSM has densities
and velocities that suggest extremely high mass-loss rates, in
some cases up to ∼ 0.1M⊙ yr−1 (Fox et al. 2011, hereafter
F11; Kiewe et al. 2012). The progenitors of SNe IIn are unknown. The implied mass-loss rates are often compared to the
episodic dense winds observed in some massive stars, such as
Luminous Blue Variables (LBVs; e.g., Smith et al. 2011, and
references therein).
SNe IIn account for more than half of all known SNe with
late-time (> 100 days) infrared (IR) emission (F11), implying the presence of warm dust. F11 showed, via a warm mission (3.6 and 4.5 µm only) Spitzer survey of 10 SNe IIn, that
the observed emission is consistent with pre-existing CSM
dust, heated by the optical emission generated by the forward shock interaction with the dense CSM. This is in contrast to most core-collapse (CC) SNe with IR emission, where
observations are consistent with a small amount of newlyformed ejecta dust (e.g., Kotak et al. 2005; Sugerman et al.
2006; Williams et al. 2008). Since the dust is produced in
the outflow of the progenitor system, characterizing this dust
may offer a clue about the pre-SN system. However, as F11
showed, Spitzer observations at 3.6 and 4.5 µm are sufficient
only to establish the presence of a warm dust component, and
not to determine the composition of the dust.
It is well-known that massive stars make dust in their stellar
winds, but what type of dust depends on the star’s mass and
evolutionary state. In general, stars make either silicate or carbon dust in their outflows, where the determining factor is the
C/O ratio in the stellar atmosphere. The presence or absence
1 CRESST and X-ray Astrophysics Laboratory, NASA/GSFC, 8800
Greenbelt Road, Greenbelt, MD, Code 662, brian.j.williams@nasa.gov
2 Department of Astronomy, University of California, Berkeley, CA
94720-3411
of a broad, strong emission feature at 9.7 µm that only arises
from silicate dust provides a method by which to distinguish
the two grain populations, provided images or spectra exist in
the 8-13 µm range.
1.1. Potential SNe IIn Progenitors and their Dust Signatures
The LBV scenario is a promising candidate for the progenitor systems of SNe IIn. Gal-Yam & Leonard (2009)
identify a hypergiant progenitor undergoing LBV-like mass
loss in pre-explosion images of SN 2005gl. However, it
has several complications. Most LBVs achieve observed
˙ < 10−4 M⊙ yr−1 in their S Doradus state
rates of only M
(Humphreys & Davidson 1994). Standard radiation and linedriven wind models similarly have difficulty reproducing
such substantial rates (e.g., Smith & Owocki 2006). Classical stellar evolution models suggest these massive stars
should evolve into the Wolf-Rayet (WR) phase before exploding (Heger et al. 2003). More recently, proposed solutions
include swept up Wolf-Rayet winds (Dwarkadas 2011), binary accretion (Kashi 2010), and gravity-wave driven winds
(Quataert & Shiode 2012), but these all have limitations of
their own.
While often variable and unpredictable in brightness, LBVs
are typically very luminous in the IR due to the significant
amounts of dust formed in their outflows. Though hardly a
homogeneous bunch, most LBVs are consistent with some
sort of silicate dust in their spectra. Some, such as HR Car
(Umana et al. 2009), show strong features from amorphous
silicate grains, while others, like the Large Magellanic Cloud
(LMC) LBV R71, exhibit features of crystalline silicates
(Morris et al. 2008). Additionally, a few LBVs like HDE
316285 and HD 168625 (Morris et al. 2008; Umana et al.
2009), show emission features from polycyclic aromatic hydrocarbons (PAHs), though whether these PAHs are formed
in the stellar outflow or are pre-existing in the interstellar
medium (ISM) is unknown. In η Carinae, Smith (2010) report that the IR spectra may be best fit with corundum dust
(Al2 O3 ), which also shows a strong ∼ 10 µm feature. We are
Draft version July 2, 2015
Preprint typeset using LATEX style emulateapj v. 5/2/11
VISCOUS BOUNDARY LAYERS OF RADIATION-DOMINATED, RELATIVISTIC JETS. II. THE
FREE-STREAMING JET MODEL
Eric R. Coughlin1 and Mitchell C. Begelman1
arXiv:1507.00007v1 [astro-ph.HE] 30 Jun 2015
JILA, University of Colorado and National Institute of Standards and Technology, 440 UCB, Boulder, CO 80309
Draft version July 2, 2015
ABSTRACT
We analyze the interaction of a radiation-dominated jet and its surroundings using the equations
of radiation hydrodynamics in the viscous limit. In a previous paper we considered the two-stream
scenario, which treats the jet and its surroundings as distinct media interacting through radiation
viscous forces. Here we present an alternative boundary layer model, known as the free-streaming jet
model – where a narrow stream of fluid is injected into a static medium – and present solutions where
the flow is ultrarelativistic and the boundary layer is dominated by radiation. It is shown that these
jets entrain material from their surroundings and that their cores have a lower density of scatterers
and a harder spectrum of photons, leading to observational consequences for lines of sight that look
“down the barrel of the jet.” These jetted outflow models may be applicable to the jets produced
during long gamma-ray bursts and super-Eddington phases of tidal disruption events.
Subject headings: galaxies: jets – gamma-ray bursts: general – radiation: dynamics – relativistic
processes
1. INTRODUCTION
Particles, magnetic fields, and radiation all contribute
to the propulsion of relativistic jets, though the relative
contribution of each is still an open matter of debate. In
certain situations, however, the mechanism responsible
for launching the jet operates simultaneously with the
release of a large amount of energy in the form of radiation, making it likely that photons dominate the bulk
energetics. This scenario occurs, for example, during the
super-Eddington phase of jetted tidal disruption events
(TDEs), such as Swift J1644+57 (Burrows et al. 2011;
Zauderer et al. 2011) and Swift J2058+05 (Cenko et al.
2012). Radiation-dominated jets should also be present
in the collapsar picture of long gamma-ray bursts (GRBs;
Rees & Meszaros 1992; Woosley 1993; Mészáros & Rees
1993; MacFadyen & Woosley 1999; Piran 2004), where
the energy released in the form of gamma-rays is ultimately derived from accretion onto a black hole, the
associated accretion luminosity exceeding the Eddington limit of the hole by more than ten orders of magnitude. In both of these cases, the propagation of
the radiation-dominated jet is modulated by the presence of a radiation pressure-supported environment;
for super-Eddington TDEs, this environment is in the
form of a highly inflated, quasi-spherical torus of fallback debris (Coughlin & Begelman 2014a), while a “cocoon” of shocked jet material (Morsony et al. 2007;
López-Cámara et al. 2013) and the overlying stellar envelope itself (Matzner 2003; Woosley & Heger 2006) serve
as the confining medium for GRBs.
In a companion paper (Coughlin & Begelman 2015,
hereafter Paper I), we presented a model that describes
the viscous interaction of a radiation-dominated, relativistic jet with its surrounding medium. In that analysis
we treated the jet and its surroundings as two separate
eric.coughlin@colorado.edu, mitch@jila.colorado.edu
1 Department of Astrophysical and Planetary Sciences, University of Colorado, UCB 391, Boulder, CO 80309
fluids, interacting with one another via small anisotropies
in the comoving radiation field that are explicitly accounted for in the equations of radiation hydrodynamics in the viscous limit. This two-stream approximation, in agreement with the non-relativistic analysis of
Arav & Begelman (1992), demonstrates the manner in
which the shear between the two fluids carves out a region of low density material within the boundary layer
between them. We also deduced the dependence of the
boundary layer thickness on the asymptotic properties of
the jet and the ambient medium.
These models also show, however, that the contact discontinuity separating the jet and its surroundings, necessary for maintaining their respective identities, results in
the likely unphysical vanishing of the density of scatterers
along that surface of separation. The contact discontinuity also prevents the jet from entraining ambient material; since the jet in the two-stream model we considered
had an infinite amount of momentum, any entrainment
or lack thereof is formally inconsequential to the evolution of the system. However, realistic jets – those with
finite extent – will almost certainly engulf more material as they expand into their surroundings; because the
total amount of momentum in the system, which is realistically finite, must be conserved, that entrainment will
then result in an overall deceleration of the outflow that
cannot be captured in the two-stream model.
In view of these unphysical properties of the twostream treatment – the vanishing of the mass density
of scatterers along the contact discontinuity and the lack
of entrainment – we present here an alternative boundary layer scenario to describe the interaction of a relativistic, radiation-dominated jet with its surroundings.
This “free-streaming jet” model, which has a well-known
counterpart in the non-relativistic, incompressible limit
(see Chapter 10, Section 12 of Kundu & Cohen 2008), assumes that the jet is injected through a narrow opening
into a static, homogenous medium, and that far enough
from that opening the entire system can be modeled as
Astronomy & Astrophysics manuscript no. orion
July 2, 2015
c
ESO
2015
A MUSE map of the central Orion Nebula (M 42) ?,??
Peter M. Weilbacher1 , Ana Monreal-Ibero2 , Wolfram Kollatschny3 , Adam Ginsburg4 , Anna F. McLeod4 , Sebastian
Kamann3 , Christer Sandin1 , Ralf Palsa4 , Lutz Wisotzki1 , Roland Bacon7 , Fernando Selman8 , Jarle Brinchmann5 ,
Joseph Caruana1 , Andreas Kelz1 , Thomas Martinsson5, 9, 10 , Arlette Pécontal-Rousset7 , Johan Richard7 , and Martin
Wendt1, 6
1
2
arXiv:1507.00006v1 [astro-ph.GA] 30 Jun 2015
3
4
5
6
7
8
9
10
Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany
e-mail: pweilbacher@aip.de
GEPI, Observatoire de Paris, CNRS, Université Paris-Diderot, Place Jules Janssen, 92190 Meudon, France
Institut für Astrophysik, Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
ESO, European Southern Observatory, Karl-Schwarzschild Str. 2, 85748 Garching bei München, Germany
Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands
Institut für Physik und Astronomie, Universität Potsdam, D-14476 Golm, Germany
CRAL, Observatoire de Lyon, CNRS, Université Lyon 1, 9 Avenue Charles André, F-69561 Saint Genis Laval Cedex, France
ESO, European Southern Observatory, 3107 Alonso de Córdova, Santiago, Chile
Instituto de Astrofísica de Canarias (IAC), E-38205 La Laguna, Tenerife, Spain
Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain
Received May 13, 2015; accepted XXX YY, 2015
ABSTRACT
We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE
instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a
0 × 4.9
0 (corresponding to ∼ 0.76 × 0.63 pc2 ) and a contiguous wavelength coverage of 4595 . . . 9366 Å, spatially
spatial size of ∼ 5.9
sampled at 000. 2. We provide two versions with a sampling of 1.25 Å and 0.85 Å in dispersion direction. Together with variance cubes
these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information
content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed
magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a twodimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios.
We estimate physical properties of the Orion Nebula, using the emission line ratios [N ii] and [S iii] (for the electron temperature T e )
and [S ii] and [Cl iii] (for the electron density Ne ), and show two-dimensional images of the velocity measured from several bright
emission lines.
Key words. HII regions – ISM: individual objects: M 42 – open clusters and associations: individual: Trapezium cluster
1. Introduction
An H ii region is a diffuse nebula whose gas is heated and ionized by the ultraviolet radiation of early-type massive stars (see
Shields 1990; Osterbrock & Ferland 2005). H ii regions are typically found in the arms of spiral galaxies and/or irregular galaxies and present spectra with strong emission lines visible even at
cosmological distances. Galactic H ii regions in particular, can be
seen as small-scale versions of the extreme events of star formation occurring in starburst galaxies (e. g. Weilbacher et al. 2003;
Alonso-Herrero et al. 2009; García-Marín et al. 2009; Cairós
et al. 2015). As such, they are laboratories offering an invaluable opportunity to study the interplay between recent and/or ongoing star formation – in particular massive stars – and their surrounding interstellar medium, including gas and dust, at a high
level of detail.
One of the best-studied Galactic H ii regions (and the closest)
is the Orion Nebula (M 42), visible to the naked eye. It is often
one of the first objects targeted with a new instrument; first, to
?
Data products are available at http://muse-vlt.eu/science.
Based on observations made with ESO telescopes at the La Silla
Paranal Observatory under program ID 60.A-9100(A).
??
see if something new can be discovered, and second, to use the
plethora of existing observations for comparison, to validate a
new system. A review about the nebula and its stellar content
can be found in O’Dell (2001). Spectroscopic studies of the ionized gas, confined to one or several slit positions, have partially
characterised the Orion’s emission spectrum (e. g. Baldwin et al.
1991; Pogge et al. 1992; Osterbrock et al. 1992; Mesa-Delgado
et al. 2008; O’Dell & Harris 2010).
However, H ii regions are rarely as simple as the textbooklike Strömgren (1939) spheres and Orion is no exception. Indeed, M 42 is thought to be only a thin blister of ionized gas at
the near side of a giant molecular cloud (Zuckerman 1973; Israel
1978; van der Werf et al. 2013). To make most of the opportunity
to observe an H ii region at the level of detail offered by Orion,
spatially resolved maps with high-quality spectral information in
terms of depth, spatial and spectral resolution are desirable.
Probably, the most efficient way to gather this information
nowadays is the use of Integral Field Spectroscopy. Sánchez
et al. (2007) released a first dataset based on this technique to
the community (using the PPak-IFU of PMAS, Kelz et al. 2006),
mapping most of the Huygens region – the central part of the
Article number, page 1 of 16
c
ESO
2015
Astronomy & Astrophysics manuscript no. 30052015
July 2, 2015
Dust attenuation in z ∼ 1 galaxies from Herschel and 3D-HST
Hα measurements
A. Puglisi1 , G. Rodighiero1 , A. Franceschini1 , M. Talia2, 3 , A. Cimatti2 , I. Baronchelli1 , E. Daddi8 , A.
Renzini4 , K. Schawinski5 , C. Mancini1, 4 , J. Silverman10 , C. Gruppioni3 , D. Lutz6 , S. Berta7 , and S. J. Oliver9
1
arXiv:1507.00005v1 [astro-ph.GA] 30 Jun 2015
2
Dipartimento di Fisica e Astronomia, Università di Padova, vicolo dell’Osservatorio 2, I-35122, Padova, Italy
e-mail: annagrazia.puglisi@studenti.unipd.it
Dipartimento di Fisica e Astronomia, Università di Bologna, Via Ranzani 1, I-40127, Bologna, Italy
3
INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, I-40127, Bologna, Italy
4
INAF-Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio, 5, I-35122, Padova, Italy
Institute for Astronomy, Departement of Physics, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich,
Switzerland
5
6
Max-Planck-Institut für extraterrestrische Physik, Giessenbachstrasse, D-85748 Garching, Germany
7
Max-Planck-Institut für Extraterrestrische Physik (MPE), Postfach 1312, 85741, Garching, Germany
Laboratoire AIM-Paris-Saclay, CEA/DSM-CNRS-Université Paris Diderot, Irfu/Service d’Astrophysique, CEA
Saclay, Orme des Merisiers, F-91191 Gif sur Yvette, France
Astronomy Centre, Dept. of Physics & Astronomy, University of Sussex, Brighton BN1 9QH, UK
Kavli Institute for the Physics and Mathematics of the Universe, Todai Institute for Advanced Study, the University
of Tokyo, Kashiwa, Japan 277-8583 (Kavli IPMU, WPI)
8
9
10
3th May 2015
ABSTRACT
We combined the spectroscopic information from the 3D-HST survey with the PEP/Herschel data to characterize the
Hα dust attenuation properties of a sample of 79 normal star-forming galaxies at 0.7 ≤ z ≤ 1.5 in the GOODS-S
field. The sample was selected in the far-IR, at λ=100 and/or 160 µm, and only includes galaxies with a secure Hα
detection (S/N>3). From the low resolution 3D-HST spectra we measured z and F(Hα) for the whole sample, rescaling
the observed flux by a constant factor of 1.2 to remove the contamination by [NII]. The stellar masses (M? ), infrared
(LIR ) and UV luminosities (LU V ) were derived from the SEDs by fitting multi-band data from GALEX near-UV to
SPIRE 500 µm. We derived the continuum extinction (Estar (B-V)) from both the IRX=LIR /LU V ratio and the UVslope, β, and found an excellent agreement among them. Galaxies in the sample have 2.6 × 109 ≤ M? ≤ 3.5 × 1011 M ,
intense infrared luminosity (LIR > 1.2 × 1010 L ), high level of dust obscuration (0.1 ≤Estar (B-V)≤ 1.1) and strong Hα
emission (with typical fluxes Fobs (Hα) ≥ 4.1×10−17 erg/s/cm2 , uncorrected for dust extinction). The nebular extinction
was estimated by comparing the observed SFRHα and the SFRU V . We obtained f =Estar (B-V)/Eneb (B-V)=0.93±0.06,
i.e. higher than the canonical value measured in the local Universe. This result could be partially due to the adopted
selection criteria, picking up the most obscured but also the Hα brightest sources. The derived dust correction produces
a good agreement between Hα and IR+UV SFRs in our sample for objects with SFR& 20 M /yr and M? & 5 × 1010
M , while sources with lower SFR and M? seem to require a smaller f -factor (i.e. higher Hα extinction correction). Our
results then imply that the nebular extinction for our sample is comparable to the optical-UV continuum extinction
and suggest that the f -factor is a function of both M? and SFR, according with previous studies.
Key words. galaxies: high-redshift – galaxies: star-formation – dust, extinction
1. Introduction
The rate at which a galaxy converts gas into stars, the Star
Formation Rate (SFR) is a fundamental quantity to characterize the evolutionary stage of a galaxy itself. The SFR can
be measured with several methods, depending on the wavelength of the available data for a given sample (e.g. the UV
luminosity, the far-IR emission, the recombination lines or
the radio-continuum, Kennicutt 1998; Madau & Dickinson
2014) although each of the various SFR tracers suffers from
some uncertainties. Concerning the UV and optical indica-
tors, the main source of uncertainty is due to the extinction
by dust, that strongly absorbs the flux emitted by stars at
UV and optical wavelengths and re-emits it in the far-IR.
Hence, an accurate quantification of the impact of dust on
the galaxy integrated emission is crucial to precisely evaluate the SFR.
The dust distribution inside a galaxy can be described by
a two component model (e.g. Charlot & Fall 2000), including a diffuse, optically thin component that attenuates the
light from the whole stellar populations (the InterStellar
Article number, page 1 of 18
Mon. Not. R. Astron. Soc. 000, 1–14 (2015)
Printed 2 July 2015
(MN LATEX style file v2.2)
arXiv:1507.00002v1 [astro-ph.SR] 30 Jun 2015
Not a galaxy: IRAS 04186+5143, a new young stellar
cluster in the outer Galaxy
J. L. Yun1,2? , Davide Elia3 †, A. A. Djupvik4 ‡, J. M. Torrelles5 §, S. Molinari3 ¶
1 Instituto
de Astrof´ısica e Ciˆ
encias do Espa¸co – Universidade de Lisboa, Observat´
orio Astron´
omico de Lisboa, Tapada da Ajuda,
1349-018 Lisboa, Portugal
2 Departament d’Astronomia i Meteorologia, Institut de Ci`
encias del Cosmos, Universitat de Barcelona (IEEC-UB), Mart´ı i Franqu`
es 1,
E-08028 Barcelona, Spain
3 INAF - Istituto di Astrofisica e Planetologia Spaziali, via Fosso del Cavaliere 100, 00133, Rome, Italy
4 Nordic Optical Telescope, Nordic Optical Telescope, Rambla Jos´
e Ana Fern´
andez P´
erez, 7, ES-38711 Bre˜
na Baja, Spain
5 Institut de Ci`
encies de l’Espai (CSIC-IEEC) and Institut de Ci`
encies del Cosmos (UB-IEEC), Mart´ı i Franqu`
es 1, E-08028 Barcelona,
Spain
Accepted 2015 June 24. Received 2015 June 23; in original form 2015 April 29
ABSTRACT
We report the discovery of a new young stellar cluster in the outer Galaxy located at
the position of an IRAS PSC source that has been previously mis-identified as an external galaxy. The cluster is seen in our near-infrared imaging towards IRAS 04186+5143
and in archive Spitzer images confirming the young stellar nature of the sources detected. There is also evidence of sub-clustering seen in the spatial distributions of
young stars and of gas and dust.
Near- and mid-infrared photometry indicates that the stars exhibit colours compatible with reddening by interstellar and circumstellar dust and are likely to be lowand intermediate-mass YSOs with a large proportion of Class I YSOs.
Ammonia and CO lines were detected, with the CO emission well centred near
the position of the richest part of the cluster. The velocity of the CO and NH3 lines
indicates that the gas is Galactic and located at a distance of about 5.5 kpc, in the
outer Galaxy.
Herschel data of this region characterise the dust environment of this molecular
cloud core where the young cluster is embedded. We derive masses, luminosities and
temperatures of the molecular clumps where the young stars reside and discuss their
evolutionary stages.
Key words: stars: formation – infrared: stars – submillimetre: ISM – ISM: clouds –
ISM: Individual objects: IRAS 04186+5143 – ISM: dust, extinction.
1
INTRODUCTION
It has been well established that star formation occurs across
the Galactic disc and at different Galactocentric distances.
Both in the inner and in the outer Galaxy, young stellar
clusters still partly embedded in the dense gas and dust in
molecular clouds have been found (e.g. Tapia et al. 1991;
Strom, Strom & Merrill 1993; McCaughrean & Stauffer
1994; Horner, Lada & Lada 1997; Luhman et al. 1998; San-
?
E-mail: yun@oal.ul.pt
† E-mail: davide.elia@ifsi-roma.inaf.it
‡ E-mail: amanda@not.iac.es
§ E-mail: torrelles@ieec.cat
(The ICC (UB) is a CSIC-Associated Unit through the ICE)
¶ E-mail: sergio.molinari@ifsi-roma.inaf.it
tos et al. 2000). They represent current active star formation sites. The star formation activity seen throughout the
Galactic disc is possible due to the relatively large amounts
of dust that shield young forming stars from the heating of
the external interstellar radiation field. The dust produces
high values of extinction resulting in lines of sight across the
Galactic disc that are highly opaque in the optical wavelengths.
Among other effects, the large values of dust extinction
along the Galactic disc make difficult the task of achieving a complete census of the Milky Way neighbour satellite
galaxies. Even with our best instruments, we may not have
found and catalogued correctly all the stellar systems components of the Milky Way and its neighbour galaxies. Conversely, Galactic stellar systems can be wrongly classified as
extragalactic neighbours. As an example of this fact, Martin