arXiv:1502.05296v1 [hep-th] 18 Feb 2015
Lectures on the Cosmological Constant Problem
Antonio Padilla
School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, UK
E-mail: [email protected]
Abstract: These lectures on the cosmological constant problem were prepared for the X Mexican
School on Gravitation and Mathematical Physics. The problem itself is explained in detail, emphasising the importance of radiative instability and the need to repeatedly fine tune as we change our
effective description. Weinberg’s no go theorem is worked through in detail. I review a number of proposals including Linde’s universe multiplication, Coleman’s wormholes, the fat graviton, and SLED, to
name a few. Large distance modifications of gravity are also discussed, with causality considerations
pointing towards a global modification as being the most sensible option. The global nature of the
cosmological constant problem is also emphasized, and as a result, the sequestering scenario is reviewed in some detail, demonstrating the cancellation of the Standard Model vacuum energy through
a global modification of General Relativity.
arXiv:1502.05138v1 [hep-th] 18 Feb 2015
Inflation with Massive Vector Fields
Junyu Liua Yi Wangb Siyi Zhoub
a University
of Science and Technology of China,
Hefei, Anhui 230026, P.R.China
b Department of Physics, The Hong Kong University of Science and Technology,
Clear Water Bay, Kowloon, Hong Kong, P.R.China
E-mail: [email protected], [email protected], [email protected]
Abstract. We investigate the coupling between the inflaton and massive vector fields. All
renormalizable couplings with shift symmetry of the inflaton are considered. The massive
vector can be decomposed into a scalar mode and a divergence-free vector mode. We show
that the former naturally interacts with the inflaton and the latter decouples at tree level.
equil
The model in general predicts fN
L = O(1), while in some regions of the parameter space
large non-Gaussianity can arise.
Explicit Integration of Friedmann’s Equation
with Nonlinear Equations of State
arXiv:1502.05042v1 [gr-qc] 17 Feb 2015
Shouxin Chen∗
Institute of Contemporary Mathematics
School of Mathematics
Henan University
Kaifeng, Henan 475004, PR China
Gary W. Gibbons†
D. A. M. T. P.
University of Cambridge
Cambridge CB3 0WA, U. K.
Yisong Yang‡
Department of Mathematics
Polytechnic School, New York University
Brooklyn, New York 11201, U. S. A
&
NYU-ECNU Institute of Mathematical Sciences
New York University - Shanghai
3663 North Zhongshan Road, Shanghai 200062, PR China
Abstract
This paper is a continuation of our earlier study on the integrability of the
Friedmann equations in the light of the Chebyshev theorem. Our main focus will be
on a series of important, yet not previously touched, problems when the equation
of state for the perfect-fluid universe is nonlinear. These include the generalized
Chaplygin gas, two-term energy density, trinomial Friedmann, Born–Infeld, and
two-fluid models. We show that some of these may be integrated using Chebyshev’s
∗
Email address: [email protected]
Email address: [email protected]
‡
Email address: [email protected]
†
1
Acoustic signal detection through the cross-correlation
method in experiments with different signal to noise ratio
and reverberation conditions
S.Adrián-Martínez, M.Ardid*, M.Bou-Cabo, I.Felis, C.Llorens, J.A.Martínez-Mora,
M.Saldaña
Universitat Politècnica de València, Institut d’Investigació per a la Gestió Integrada de
Zones Costaneres (IGIC). C/Paranimf 1, 46730 Gandia, Spain
*[email protected]
Abstract. The study and application of signal detection techniques based on cross-
correlation method for acoustic transient signals in noisy and reverberant environments
are presented. These techniques are shown to provide high signal to noise ratio, good
signal discernment from very close echoes and accurate detection of signal arrival time.
The proposed methodology has been tested on real data collected in environments and
conditions where its benefits can be shown. This work focuses on the acoustic detection
applied to tasks of positioning in underwater structures and calibration such those as
ANTARES and KM3NeT deep-sea neutrino telescopes, as well as, in particle detection
through acoustic events for the COUPP/PICO detectors. Moreover, a method for
obtaining the real amplitude of the signal in time (voltage) by using cross correlation
has been developed and tested and is described in this work.
Keywords: Acoustic signal detection, cross-correlation method, processing techniques,
positioning, underwater neutrino telescopes, particle detectors.
1
Introduction
Acoustic signal detection has become an object of interest due to its utility and
applicability in fields such as particle detection, underwater communication, medical
issues, etc. The group of Acoustics Applied to Astroparticle Detection from the
Universitat Politècnica de València collaborates with the particle detectors ANTARES
[1], KM3NeT [2] and COUPP/PICO [3]. Acoustic technologies and processing
analyses are developed and studied for positioning, calibration and particle detection
tasks of the detectors.
Acoustic emitters and receivers are used for the positioning systems of underwater
neutrino telescopes ANTARES [4] and KM3NeT [5] in order to monitor the position
of the optical detection modules of these telescopes. The position of optical sensors
need to be monitored with 10 cm accuracy to be able to determine the trajectory of the
muon produced after a neutrino interaction in the vicinity of the telescope from the
Cherenkov light that it produces [6]. An important aspect of the acoustic positioning
New Estimates of Time Delays
in the Gravitationally Lensed System PG1115+080
E. Shimanovskaya,1 , *
arXiv:1502.05392v1 [astro-ph.CO] 18 Feb 2015
1 Moscow
V. Oknyanskiy,1 , **
B. Artamonov1, ***
MV Lomonosov State University, Sternberg Astronomical Institute,
Moscow 119992, Russia
We analyze all publicly available long-term optical observations of the gravitationally lensed quasar PG1115+080 for the purpose of estimating time delays between
its four components. In particular, the light curves of PG1115+080 components
obtained in 2001-2006 at Maidanak observatory (Uzbekistan) [60] are considered. We
find that the linear trend is observed in 2006 in light curves of all four components
with fast variations only in the A1 and C components that can be due to microlensing
and observational errors. Application of the MCCF method [36] to the photometric
+2
data obtained in 2004-2005 gives values of time delays τBC = 22+2
−3 , τAC = 12−1
and τBA = 10+2
−3 days, which are in agreement with the results received earlier by
Schechter [48] and Barkana [2] for 1995-1996 light curves with two different methods
of statistic analysis. However, our estimates of τBA and τBC differ from the values
received by the group of Vakulik based on the same Maidanak data [61]. The ratio
τAC /τBA is equal to ∼ 1.2 that is close to the value, received by Barkana (∼ 1.13)
and predicted by lens models (∼ 1.4), unlike the values received by Schechter (∼ 0.7)
and Vakulik (∼ 2.7).
1. INTRODUCTION
According to General relativity, beams of light deviate from the straight line under the
inﬂuence of the gravitational ﬁeld of a massive object [11]. If a galaxy is located close to the
line of sight between a far quasar and an observer, the light from the quasar passes diﬀerent
*
Electronic address: [email protected]
**
Electronic address: [email protected]
***
Electronic address: [email protected]
Preprint typeset using LATEX style emulateapj v. 04/17/13
HUBBLE SPACE TELESCOPE SPECTROSCOPY OF BROWN DWARFS DISCOVERED WITH THE
WIDE-FIELD INFRARED SURVEY EXPLORER
arXiv:1502.05365v1 [astro-ph.SR] 18 Feb 2015
Adam C. Schneidera , Michael C. Cushinga , J. Davy Kirkpatrickb , Christopher R. Gelinob,c , Gregory N. Maced,e ,
Edward L. Wrightd , Peter R. Eisenhardtf , M. F. Skrutskieg , Roger L. Griffithh , & Kenneth A. Marshi
ABSTRACT
We present a sample of brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE)
for which we have obtained Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) near-infrared
grism spectroscopy. The sample (twenty-two in total) was observed with the G141 grism covering
1.10−1.70 µm, while fifteen were also observed with the G102 grism, which covers 0.90−1.10 µm. The
additional wavelength coverage provided by the G102 grism allows us to 1) search for spectroscopic
features predicted to emerge at low effective temperatures (e.g. ammonia bands) and 2) construct
a smooth spectral sequence across the T/Y boundary. We find no evidence of absorption due to
ammonia in the G102 spectra. Six of these brown dwarfs are new discoveries, three of which are found
to have spectral types of T8 or T9. The remaining three, WISE J082507.35+280548.5 (Y0.5), WISE
J120604.38+840110.6 (Y0), and WISE J235402.77+024015.0 (Y1) are the nineteenth, twentieth, and
twenty-first spectroscopically confirmed Y dwarfs to date. We also present HST grism spectroscopy
and reevaluate the spectral types of five brown dwarfs for which spectral types have been determined
previously using other instruments.
Keywords: stars: low-mass, brown dwarfs: individual (WISEA J032504.52−504403.0, WISEA
J040443.50−642030.0, WISEA J082507.37+280548.2, WISEA J120604.25+840110.5,
WISEA J221216.27−693121.6, WISEA J235402.79+024014.1)
1. INTRODUCTION
The Wide-field Infrared Survey Explorer (WISE) has
been very successful at identifying the coolest brown
dwarfs in the Solar neighborhood. WISE provided all
sky coverage at four mid-infrared wavelengths centered
at 3.4, 4.6, 12 and 22 µm (W 1, W 2, W 3, and W 4).
This coverage has offered the ideal dataset with which to
identify cool brown dwarfs whose spectral energy distributions peak at mid-infrared wavelengths. Specifically, WISE was designed so that the W 1 band coincides
with a deep water+methane absorption feature and the
W 2 band coincides with a region largely free of opacity. Thus, the WISE W 1 − W 2 color has been especially
useful for identifying late type dwarfs (Kirkpatrick et al.
2011, Cushing et al. 2011, Kirkpatrick et al. 2012, Mace
et al. 2013, and Thompson et al. 2013). Seventeen of
the eighteen spectroscopically confirmed brown dwarfs
a Department of Physics and Astronomy, University of
Toledo, 2801 W. Bancroft St., Toledo, OH 43606, USA;
[email protected]
b Infrared Processing and Analysis Center, MS 100-22, California Institute of Technology, Pasadena, CA 91125, USA
c NASA Exoplanet Science Institute, Mail Code 100-22, California Institute of Technology, 770 South Wilson Ave, Pasadena,
CA 91125, USA
d Department of Physics and Astronomy, UCLA, 430 Portola
Plaza, Box 951547, Los Angeles, CA 90095-1547, USA
e Department of Astronomy, The University of Texas at
Austin, 2515 Speedway, Stop C1400, Austin, TX 78712, USA
f Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA
g Department of Astronomy, University of Virginia, 530 McCormick Road, Charlottesville, VA 22904, USA
h Department of Astronomy & Astrophysics, 525 Davey Lab,
The Pennsylvania State University, University Park, PA 16802,
USA
i School of Physics and Astronomy, Cardiff University, Cardiff
CF24 3AA, UK
with spectral types of Y0 or later were first identified by
WISE (Cushing et al. 2011, Kirkpatrick et al. 2012, Tinney et al. 2012, Kirkpatrick et al. 2013, Cushing et al.
2014a, and Pinfield et al. 2014). The Y0 dwarf WISE
J1217+16B was identified as a companion to a T8.5
dwarf, itself identified by WISE (Liu et al. 2012, Leggett
et al. 2014a). Three objects, WD 0806−661(Luhman
et al. 2011), CFBDSIR J1458+1013B (Liu et al. 2011),
and WISE 0855−0714 (Luhman 2014), likely have effective temperatures similar to (or less than) the above Y
dwarfs, but have yet to be spectroscopically confirmed.
Ground-based follow-up observations at the midinfrared wavelengths where the spectral energy distributions of cold brown dwarfs peak are nearly impossible due
to the high thermal background, forcing follow-up observations to shorter wavelengths where these dwarfs are extremely faint. As a result, obtaining moderate signal-tonoise (S/N) near-infrared spectra for many of the coldest
WISE candidates has only been capable with the Wide
Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST). HST grism spectroscopy of the latest type
brown dwarfs has been invaluable in the study and classification of these objects (Cushing et al. 2011, Kirkpatrick
et al. 2012, Kirkpatrick et al. 2013, and Cushing et al.
2014a).
As effective temperatures cool below 600 K, several
spectroscopic features are predicted by model atmospheres to arise in the Y-band spectral region around
1.07 µm. These include the emergence of ammonia absorption components and the disappearance of optical
alkali resonance lines which have broad wings that are
predicted to extend into the near-infrared (Burrows et
al. 2000, Kirkpatrick et al. 2012). (HST) WFC3 G102
spectroscopy (0.90−1.10 µm) allows us to investigate this
additional wavelength range for such features, as well as
inspect for differences (and similarities) as a function of
arXiv:1502.05362v1 [astro-ph.SR] 18 Feb 2015
High-speed Photometric Observations of ZZ Ceti White Dwarf
Candidates
E. M. Green,1 M.-M. Limoges,2 A. Gianninas,3 P. Bergeron,2 G. Fontaine,2
P. Dufour,2 C. J. O’Malley,1 B. Guvenen,1 L. I. Biddle,1 K. Pearson,1
T. W. Deyoe,1 C. W. Bullivant,1 J. J. Hermes,4 V. Van Grootel,5 and
M. Grosjean5
1 Steward
Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson,
AZ 85721, USA; bgreen@ as.arizona.edu
2 Département
de Physique, Université de Montréal, C. P. 6128, Succ.
Centre-Ville, Montréal, Québec, H3C 3J7, Canada
3 Homer
L. Dodge Department of Physics & Astronomy, University of
Oklahoma, 440 W. Brooks St., Norman, OK 73019, USA
4 Department
of Physics, University of Warwick, Coventry, CV4 7AL, UK
5 Institut
d’Astrophysique et de Géophysique de l’Université de Liège, Allée du
6 Août 17, B-4000 Liège, Belgium
Abstract. We present high-speed photometric observations of ZZ Ceti white dwarf
candidates drawn from the spectroscopic survey of bright DA stars from the Villanova
White Dwarf Catalog by Gianninas et al., and from the recent spectroscopic survey of
white dwarfs within 40 parsecs of the Sun by Limoges et al. We report the discovery of
six new ZZ Ceti pulsators from these surveys, and several photometrically constant DA
white dwarfs, which we then use to refine the location of the ZZ Ceti instability strip.
1.
Introduction
We have recently completed two major spectroscopic surveys of DA stars using the
spectroscopic approach where hydrogen Balmer lines are fitted with the predictions
of detailed model atmospheres. Gianninas et al. (2011) conducted a spectroscopic
survey of over 1300 bright (V < 17.5), hydrogen-rich white dwarfs based largely on
the last published version of the McCook & Sion catalog (McCook & Sion 1999).
The Gianninas et al. sample included 56 known ZZ Ceti stars, 145 photometrically
constant DA stars as well as several white dwarfs whose atmospheric parameters placed
them within or near the empirical boundaries of the ZZ Ceti instability strip (see Figure
35, Gianninas et al. 2011). This includes the ultra-massive ZZ Ceti star, GD 518
(WD 1659+662), discovered by Hermes et al. (2013). More recently, Limoges et
al. (2013, 2014) performed an exhaustive spectroscopic survey of the SUPERBLINK
proper motion database (see Lépine & Shara 2005, and references therein) aimed at
obtaining a complete sample of white dwarfs in the solar neighborhood within 40 pc of
the Sun. Several ZZ Ceti white dwarf candidates were also identified in this survey.
1
Delensing the CMB with the Cosmic Infrared Background
Blake D. Sherwin∗1, 2, 3 and Marcel Schmittfull3
1
Department of Physics, University of California, Berkeley, CA, 94720, USA
Miller Institute for Basic Research in Science, University of California, Berkeley, CA, 94720, USA
3
Berkeley Center for Cosmological Physics, LBNL and University of California, Berkeley, CA, 94720, USA
arXiv:1502.05356v1 [astro-ph.CO] 18 Feb 2015
2
As confusion with lensing B-modes begins to limit experiments that search for primordial B-mode
polarization, robust methods for delensing the CMB polarization sky are becoming increasingly
important. We investigate in detail the possibility of delensing the CMB with the cosmic infrared
background (CIB), emission from dusty star-forming galaxies that is an excellent tracer of the CMB
lensing signal, in order to improve constraints on the tensor-to-scalar ratio r. We find that the maps
of the CIB, such as current Planck satellite maps at 545 GHz, can be used to remove more than half
of the lensing B-mode power. Calculating optimal combinations of different large-scale-structure
tracers for delensing, we find that co-adding CIB data and external arcminute-resolution CMB
lensing reconstruction can lead to significant additional improvements in delensing performance.
We investigate whether measurement uncertainty in the CIB spectra will degrade the delensing
performance if no model of the CIB spectra is assumed, and instead the CIB spectra are marginalized
over, when constraining r. We find that such uncertainty does not significantly affect B-mode surveys
smaller than a few thousand degrees. Even for larger surveys it causes only a moderate reduction
in CIB delensing performance, especially if the surveys have high (arcminute) resolution, which
allows self-calibration of the delensing procedure. Though further work on the impact of foreground
residuals is required, our overall conclusions for delensing with current CIB data are optimistic: this
delensing method can tighten constraints on r by a factor up to ≈ 2.2, and by a factor up to ≈ 4
when combined with external ≈ 3µK-arcmin lensing reconstruction, without requiring the modeling
of CIB properties. CIB delensing is thus a promising method for the upcoming generation of CMB
polarization surveys.
I.
INTRODUCTION
Many inflationary models make a key prediction: the
production of a significant background of gravitational
waves [1–6]. A detection of this stochastic gravitational
wave background would therefore provide particularly
strong evidence for the inflationary paradigm, alongside the observed super-horizon correlation, near scaleinvariance and Gaussianity of the scalar density fluctuations (e.g., [7] and references therein). Such a detection
would also give significant insights into the physics of
inflation. In particular, the ratio of gravitational wave
or tensor power to scalar density power r gives a direct
probe of the energy scale at which inflation occurred.
The currently most promising method to measure this
inflationary gravitational wave background is through
its effects on the cosmic microwave background (CMB).
However, tensor modes can be difficult to constrain accurately via the CMB temperature alone, because scalar
modes also contribute to the temperature, which leads
to cosmic variance and degeneracies. Unlike in temperature, the CMB sky contains two degrees of freedom in
polarization (a magnitude of polarization and a direction); because of this, one can construct a “null” linear
combination from the two degrees of freedom which has
no response to scalar perturbations at linear order, but
still is sourced by inflationary gravitational waves. This
combination, known as the CMB B-mode polarization,
∗
email: [email protected]
therefore allows a constraint on primordial gravitational
waves that is robust and free of leading-order scalar cosmic variance [8–10].
Over the past year, the observational status of primordial B-mode polarization has been a subject of active debate. While the BICEP2 collaboration reported
a measurement of B-mode polarization on large angular
scales [11], subsequent analyses based on Planck polarization measurements suggested that this signal could be
sourced by foreground emission from polarized galactic
dust [12–14]. A recently-released cross-correlation analysis involving data from the BICEP2, Keck Array and
Planck experiments [15] has confirmed that a significant
fraction of the signal seen in BICEP2 was of galactic origin, and has thus replaced any evidence for inflationary
gravitational waves with a new upper limit of r < 0.12
at 95% confidence. For simplicity, and motivated by this
new result, we work under the default assumption that
r = 0 (though cosmic variance from tensor modes could
straightforwardly be incorporated into our analysis if future measurements required this). While foreground polarization appears to be a significant obstacle when analyzing only one frequency channel, component separation
with multifrequency data is a promising path towards removing much of the foreground contamination.
Even in the absence of foregrounds, primordial gravitational waves are not the only source of B-mode polarization. A well-understood, recently observed [11, 15–21]
source of B-mode polarization is the gravitational lensing
of the CMB, the deflection of CMB photons by the distribution of large-scale structure along the photons path.
By remapping the CMB sky, gravitational lensing con-
In Search of Future Earths: Assessing the possibility of finding Earth analogues in the
later stages of their habitable lifetimes
Running title: In Search of Future Earths
Jack T. O'Malley-James: School of Physics and Astronomy, University of St Andrews, North
Haugh, St Andrews, Fife, UK.
Jane S. Greaves: School of Physics and Astronomy, University of St Andrews, North Haugh,
St Andrews, Fife, UK.
John A. Raven: Division of Plant Sciences, University of Dundee at TJHI, The James Hutton
Institute, Invergowrie, Dundee, UK.
Charles S. Cockell: UK Centre for Astrobiology, School of Physics and Astronomy, James
Clerk Maxwell Building, The King's Buildings, University of Edinburgh, Edinburgh, UK.
Corresponding author:
J.T. O'Malley-James
E.mail: [email protected]
1
Earth will become uninhabitable within the next 2-3 billion years as a result of the moving
boundaries of the habitable zone caused by the increasing luminosity of the Sun. Predictions
about the future of habitable conditions on Earth include a decline in species diversity and
habitat extent, ocean loss and changes in the magnitudes of geochemical cycles. However,
testing these predictions on the present-day Earth is difficult. The discovery of a planet that is
a near analogue to the far future Earth (an old-Earth-analogue) could provide a means to test
these predictions. Such a planet would need to have an Earth-like biosphere history, requiring
it to have been in its system's habitable zone for giga-year (Gyr) long periods during the
system's past, and to be approaching the inner-edge of the habitable zone at present. Here we
assess the possibility of finding this very specific type of exoplanet and discuss the benefits
of analysing older Earths in terms of improving our understanding of long-term geological
and bio-geological processes. Finding such a planet in nearby star systems would be ideal,
because it would be close enough to allow for atmospheric characterisation. Hence, as an
illustrative example, G stars within 10 parsecs (pc) of the Sun are assessed as potential oldEarth-analogue hosts. Six of these stars are at appropriate stages of their main sequence
evolution to be good potential hosts. For each of these systems, a hypothetical Earth analogue
is placed at locations within the continuously habitable zone (CHZ) that would allow enough
time for Earth-like biosphere development. Planetary surface temperature evolution over the
host star's main sequence lifetime is assessed using a simple climate model. This is then used
to determine whether the planet would be in the right stage of its late-habitable lifetime to
exhibit detectable biosignatures. The best candidate in terms of the chances of planet
formation in the CHZ and of biosignature detection is 61 Virginis. Predictions from planet
formation studies and biosphere evolution models suggest that only a small fraction (0.36%)
of G stars in the solar neighbourhood could host an old-Earth-analogue. However, if the
development of an Earth-like biosphere is assumed to be rare, requiring a sequence of lowprobability events to occur, then such planets are unlikely to be found in the solar
neighbourhood − although 1000s could be present in the galaxy as a whole.
INTRODUCTION
The fate of Earth, and of any habitable planet, is to eventually become uninhabitable as stellar
evolution pushes the habitable zone boundaries further out over time (Kasting et al., 1993;
Franck et al., 2000; O’Malley-James et al. 2013; Rushby et al., 2013). As an inhabited planet
evolves towards this climate state, the types and abundance of life change (Kasting et al.,
1993; Franck et al., 2000; O’Malley-James et al. 2013; Rushby et al., 2013). This can be
predicted using climate and solar evolution models and by fitting the expected environmental
conditions to the known physical and chemical limits of life. However, the means to test
these predictions do not exist. Inferences about past life and climates on Earth can be
supported or rejected based on geological evidence (Mojzsis et al., 1996; Nisbet & Sleep,
2001; Lepland et al., 2005). However, beyond extrapolations from the evolutionary responses
of species to environmental changes (Thuiller et al., 2008), Earth provides very little evidence
for predictions made about the planet’s future biosphere. The only possible method for testing
2
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOI:10.1002/,
Magnetic field-line lengths inside interplanetary
magnetic flux ropes
1
2
arXiv:1502.05284v1 [astro-ph.SR] 18 Feb 2015
Qiang Hu, Jiong Qiu, and Sam Krucker
1
3
Department of Space Science and
CSPAR, The University of Alabama in
Huntsville, Huntsville, AL, USA.
2
Department of Physics, Montana State
University, Bozeman, MT, USA.
3
Space Science Laboratory, University of
California, Berkeley, CA, USA.
D R A F T
X-2
Abstract.
February 19, 2015, 1:21am
D R A F T
HU ET AL.: MAGNETIC FIELD-LINE LENGTHS IN FLUX ROPES
We report on the detailed and systematic study of field-line
twist and length distributions within magnetic flux ropes embedded in Interplanetary Coronal Mass Ejections (ICMEs). The Grad-Shafranov reconstruction method is utilized together with a constant-twist nonlinear forcefree (Gold-Hoyle) flux rope model to reveal the close relation between the
field-line twist and length in cylindrical flux ropes, based on in-situ Wind
spacecraft measurements. We show that the field-line twist distributions within
interplanetary flux ropes are inconsistent with the Lundquist model. In particular we utilize the unique measurements of magnetic field-line lengths within
selected ICME events as provided by Kahler et al. [2011a] based on energetic
electron burst observations at 1 AU and the associated type III radio emissions detected by the Wind spacecraft. These direct measurements are compared with our model calculations to help assess the flux-rope interpretation
of the embedded magnetic structures. By using the different flux-rope models, we show that the in-situ direct measurements of field-line lengths are consistent with a flux-rope structure with spiral field lines of constant and low
twist, largely different from that of the Lundquist model, especially for relatively large-scale flux ropes.
D R A F T
February 19, 2015, 1:21am
D R A F T
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 05/12/14
TOWARD AN EMPIRICAL THEORY OF PULSAR EMISSION XI. UNDERSTANDING THE ORIENTATIONS
OF PULSAR RADIATION AND SUPERNOVA “KICKS”
Joanna M. Rankin
arXiv:1502.05270v1 [astro-ph.HE] 18 Feb 2015
Physics Department, University of Vermont, Burlington, VT 05405 USA:
[email protected]
Draft version February 19, 2015
ABSTRACT
Two entwined problems have remained unresolved since pulsars were discovered nearly 50 years ago:
the orientation of their polarized emission relative to the emitting magnetic field and the direction of
putative supernova “kicks’ relative to their rotation axes. The rotational orientation of most pulsars
can be inferred only from the (“fiducial’) polarization angle of their radiation, when their beam points
directly at the Earth and the emitting polar fluxtube field is k to the rotation axis. Earlier studies
have been unrevealing owing to the admixture of different types of radiation (core and conal, two
polarization modes), producing both k or ⊥ alignments. In this paper we analyze the some 50 pulsars
having three characteristics: core radiation beams, reliable absolute polarimetry, and accurate proper
motions. The “fiducial’ polarization angle of the core emission, we then find, is usually oriented ⊥ to
the proper-motion direction on the sky. As the primary core emission is polarized ⊥ to the projected
magnetic field in Vela and other pulsars where X-ray imaging reveals the orientation, this shows that
the proper motions usually lie k to the rotation axes on the sky. Two key physical consequences
then follow: first, to the extent that supernova “kicks’ are responsible for pulsar proper motions, they
are mostly k to the rotation axis; and second that most pulsar radiation is heavily processed by the
magnetospheric plasma such that the lowest altitude “parent’ core emission is polarized ⊥ to the
emitting field, propagating as the extraordinary (X) mode.
Keywords: pulsars: general — techniques: polarimetric; emission mechanisms: non-thermal
1. INTRODUCTION
Radio pulsars now contribute importantly to many
fields of physical science, but paradoxically, two fundamental coupled issues have remained unresolved since
they were discovered 47 years ago: the orientation of
their linearly polarized emission relative to the magnetic field in their polar fluxtube emitting regions and
the origin/orientation of their often large space velocities (and thus proper motions) relative to their rotation
axes. Most pulsars are known only by their lighthouselike radio signals, and thus we have no direct means of
determining the orientation of a pulsar’s rotation axis
on the sky, which is crucial both to interpreting the polarization direction of the radiation and the orientation
of the proper motion relative to the spin axis. Pulsars
radiate because highly energetic outward-going charges
are accelerated by the curved dipolar field within their
polar regions, so it is crucial to understand how this radiation is polarized relative to the (projected) B field
orientation on the sky. Figure 1 shows how this field
appears splayed when a pulsar’s beam points directly at
the Earth—the “fiducial” instant—and that the emission
reaching us is associated with that bundle in the plane
of the rotation axis Ω. Clearly, we have no knowledge
of the radial component of a pulsar’s space velocity and
only weak estimates of the radial component of Ω.
At one time it seemed obvious that pulsar radiation
must be polarized k to the projected magnetic field direction. How could it be otherwise given that both the
curving B field and the resulting curvature acceleration
lie in the same plane? Even after the discovery that pulsars emit in two orthogonal (hereafter OPM) polariza-
tion modes [Backer & Rankin (1980); Manchester et al.
(1975)] many assumed in the absence of any direct proof
that the “primary” polarization mode must be k.
This easy presumption was dashed in the new millennium by X-ray imaging of the Vela pulsar [Helfand et al.
(2001); Radhakrishnan & Deshpande (2001)] where arcs
indicated the orientation of the star’s rotation axis Ω
relative to its polarization and proper-motion (PM) directions. Shockingly, the radiation was polarized ⊥
to the projected magnetic field B plane, a circumstance then beautifully confirmed for the radio emission
(Johnston et al. 2005, hereafter Johnston I)—and this
pulsar’s radio emission is almost completely linearly polarized, so there was no OPM ambiguity.
In an earlier paper (Rankin 2007, hereafter Paper I)
we investigated the PPA vs. PM alignments of a number of pulsars, drawing strongly on Johnston I as well
as other sources. Here, a “fiducial” polarization position
angle (PPA) P A0 , at a (“fiducial”) rotational phase representing the magnetic axis longitude, is measured and
referred to infinite frequency as a proxy for the (unseen)
orientation of the rotation axis Ω. These were compared with well determined proper-motion (PM) directions P AV , and the differences Ψ showed strong peaks
at both 0◦ and 90◦ . Given, however, that most of the
pulsars showed strong OPM activity in their profiles, it
was not possible to draw general conclusions about the
polarization orientation with respect to the projected B
direction.
The second coupled key question is how a pulsar’s rotation vector Ω is oriented with respect to its space velocity (of which we can usually measure only their projections on the plane of the sky)? The possibility of a cor-
A composition dependent energy scale and the determination of the cosmic ray primary
mass in the ankle region
A. D. Supanitskya, A. Etchegoyenb, D. Melob , F. Sanchezb
arXiv:1502.05246v1 [astro-ph.HE] 18 Feb 2015
b Instituto
a Instituto de Astronom´ıa y F´ısica del Espacio (IAFE), CONICET-UBA, Buenos Aires, Argentina.
de Tecnolog´ıas en Detecci´on y Astropart´ıculas (CNEA, CONICET, UNSAM), Centro At´omico Constituyentes, San Martin, Buenos Aires, Argentina.
Abstract
At present there are still several open questions about the origin of the ultra high energy cosmic rays. However, great progress
in this area has been made in recent years due to the data collected by the present generation of ground based detectors like the
Pierre Auger Observatory and Telescope Array. In particular, it is believed that the study of the composition of the cosmic rays as
a function of energy can play a fundamental role for the understanding of the origin of the cosmic rays.
The observatories belonging to this generation are composed of arrays of surface detectors and fluorescence telescopes. The
duty cycle of the fluorescence telescopes is ∼ 10 % in contrast with the ∼ 100 % of the surface detectors. Therefore, the energy
calibration of the events observed by the surface detectors is performed by using a calibration curve obtained from a set of high
quality events observed in coincidence by both types of detectors. The advantage of this method is that the reconstructed energy of
the events observed by the surface detectors becomes almost independent of simulations of the showers because just a small part of
the reconstructed energy (the missing energy), obtained from the fluorescence telescopes, comes from simulations. However, the
calibration curve obtained in this way depends on the composition of the cosmic rays, which can introduce biases in composition
analyses when parameters with a strong dependence on primary energy are considered. In this work we develop an analytical
method to study these effects. We consider AMIGA (Auger Muons and Infill for the Ground Array), the low energy extension of
the Pierre Auger Observatory corresponding to the surface detectors, to illustrate the use of the method. In particular, we study
the biases introduced by an energy calibration dependent on composition on the determination of the mean value of the number of
muons, at a given distance to the showers axis, which is one of the parameters most sensitive to primary mass and has an almost
linear dependence with primary energy.
Keywords: Cosmic Rays, Energy Calibration, Chemical Composition
1. Introduction
The cosmic ray energy spectrum extends over more than
eleven orders of magnitude in energy (from below ∼ 109 eV
to above ∼ 1020 eV). It can be approximated by a broken power
law with four spectral features: the knee at a few 1015 eV [1–5],
the ankle at ∼ 4 × 1018 eV [6–12], the cutoff or suppression at
∼ 3 × 1019 eV [9–12], and a second knee at ∼ 1017 eV, recently
reported by the KASCADE-Grande Collaboration [13].
Several experimental techniques are used for the observation of the cosmic rays, depending on the energy range under
consideration. In particular, the direct observation of the primary particles is possible up to energies of the order of ∼ 1015
eV. For larger energies the study of cosmic rays is done by
observing the atmospheric air showers that they generate as a
consequence of their interactions with air molecules in the atmosphere. There are two classes of ground-based detectors,
surface detectors and fluorescence telescopes. The surface detectors observe the lateral distribution of the showers by sampling the secondary particles that reach the Earth’s surface,
whereas the fluorescence telescopes observe the fluorescence
Email address: [email protected] (A. D. Supanitsky)
Preprint submitted to Astroparticle Physics
and Cherenkov photons generated, during the longitudinal development of the showers, as a result of the interaction of the
secondary charged particles with the air molecules [14–17].
Despite great experimental effort done in the last years the
origin of the cosmic rays is still unknown. The observations
used to study their origin comprise: the energy spectrum, the
distribution of the arrival directions, and the composition profile
[18, 19].
Certainly, the detailed study of the composition as a function of energy is of great importance to unveil the origin of
the cosmic rays at all energies (see Ref. [20] for a review on
composition). In particular, it is believed that the composition
information is crucial to find the transition between the galactic and extragalactic components of the cosmic rays (see for
instance Ref. [21]) and to elucidate the origin of the suppression at the highest energies [22]. This feature of the spectrum
could originate as a result of the propagation of the cosmic rays
in the intergalactic medium, or by the end of the efficiency of
the extragalactic sources to accelerate particles at the highest
energies, or by a combination of both effects.
At the highest energies (E & 1015 eV), the composition of
the cosmic rays is studied by using different observable parameters obtained from shower measurements which are very
February 19, 2015
arXiv:1502.05237v1 [astro-ph.SR] 18 Feb 2015
New H−band Stellar Spectral Libraries for the
SDSS–III/APOGEE survey
O. Zamora1,2, D. A. Garc´ıa-Hern´andez1,2 , C. Allende Prieto1,2 , R. Carrera1,2 , L.
Koesterke3 , B. Edvardsson4 , F. Castelli5 , B. Plez6 , D. Bizyaev7,8 , K. Cunha9 , A. E. Garc´ıa
P´erez1,2 , B. Gustafsson4 , J. A. Holtzman8 , J. E. Lawler10 , S. R. Majewski11 , A.
Manchado1,2,12 , Sz. M´esz´aros13 , N. Shane14 , M. Shetrone15 , V. V. Smith16 , G. Zasowski17
–4–
ABSTRACT
The Sloan Digital Sky Survey–III (SDSS–III) Apache Point Observatory
Galactic Evolution Experiment (APOGEE) has obtained high resolution (R ∼
22,500), high signal-to-noise (> 100) spectra in the H−band spectral region
(∼1.5−1.7 µm) for about 146,000 stars in the Milky Way galaxy. We have computed specific spectral libraries with effective temperature (T eff ) ranging from
3500 to 8000 K for the automated chemical analysis of the survey data. The
spectral libraries, used to derive stellar parameters and abundances from the
APOGEE spectra in the SDSS–III data release 12 (DR12), are based on ATLAS9 model atmospheres and the ASSǫT spectral synthesis code. We present
a second set of stellar spectral libraries based on MARCS model atmospheres
and the spectral synthesis code Turbospectrum. The ATLAS9/ASSǫT (T eff =
3500−8000 K) and MARCS/Turbospectrum (T eff = 3500−5500 K) grids of synthetic spectra cover a wide range of metallicity (−2.5 ≤ [M/H] ≤ +0.5 dex),
surface gravity (0 ≤ log g ≤ 5 dex), microturbulence (0.5 ≤ ξ ≤ 8 km s−1 ), carbon (−1 ≤ [C/M] ≤ +1 dex), nitrogen (−1 ≤ [N/M] ≤ +1 dex), and α-elements
(−1 ≤ [α/M] ≤ +1 dex) variations, having thus seven dimensions. We compare
the ATLAS9/ASSǫT and MARCS/Turbospectrum spectral libraries and we apply both of them to the observed H−band spectra of the Sun and the K2 giant
Arcturus, as well as to a selected sample of well-known giant stars observed at
very high-resolution. The new APOGEE synthetic spectral libraries are publicly
available online and can be employed for chemical studies in the H−band using
other high-resolution spectrographs.
Subject headings: astrochemistry – radiative transfer – stars: atmospheres – surveys
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 2/19/04
THE COMPLEX NORTH TRANSITION REGION OF CENTAURUS A: A GALACTIC WIND
Susan G. Neff
NASA’s Goddard Space Flight Center, Laboratory for Observational Cosmology, Mail Code 665, Greenbelt, Maryland, 20771
Jean A. Eilek
arXiv:1502.05236v1 [astro-ph.GA] 18 Feb 2015
Physics Department, New Mexico Tech, Socorro NM 87801 and
National Radio Astronomy Observatory∗† , Socorro NM 87801
Frazer N. Owen
National Radio Astronomy Observatory∗ , Socorro NM 87801
Draft version February 19, 2015
ABSTRACT
We present deep GALEX images of NGC 5128, the parent galaxy of Centaurus A. We detect a
striking “weather ribbon” of Far-UV and Hα emission which extends more than 35 kpc northeast of
the galaxy. This ribbon is associated with a knotty ridge of radio/X-ray emission and is an extension
of the previously known string of optical emission-line filaments. Many phenomena in the region are
too short-lived to have survived transit out from the inner galaxy; something must be driving them
locally. We also detect Far-UV emission from the galaxy’s central dust lane. Combining this with
previous radio and Far-IR measurements, we infer an active starburst in the central galaxy which is
currently forming stars at ∼ 2M yr−1 , and has been doing so for 50-100 Myr. If the wind from this
starburst is enhanced by energy and mass driven out from the AGN, the powerful augmented wind
can be the driver needed for the northern weather system. We argue that both the diverse weather
system, and the enhanced radio emission in the same region, result from the wind’s encounter with
cool gas left by one of the recent merger/encounter events in the history of NGC 5128.
Subject headings: galaxies: active – galaxies: individual (NGC 5128, Centaurus A) – galaxies: jets –
galaxies: starburst – galaxies: winds – radio continuum: galaxies
1. INTRODUCTION
The radio source Centaurus A (“Cen A”) and its parent galaxy, NGC 5128 are the nearest active-galaxy system. At a distance of only 3.8 Mpc, (Harris et al. 2010;
10 ' 1.14kpc), both the radio source Cen A and the galaxy
NGC5128 can be scrutinized with a sensitivity and resolution impossible for other active galaxies.
NGC 5128 is fundamentally a normal elliptical galaxy,
dominated by an old stellar population, with kinematic
signatures typical of other massive ellipticals (e.g., Peng et
al. 2004; Woodley et al. 2010; Rejkuba et al. 2011). However, the galaxy’s optical appearance is dominated by the
unusual, and iconic, dust band (Dunlop 1828, Herschel
1847). The central dust band is the site of vigorous and
ongoing star formation. In this paper, we show that the
extended starburst is strong enough to drive out a wind.
On larger scales, a striking complex of optical emissionline gas, active star formation, cold gas and dust clouds,
and overpressured radio and X-ray knots extends into the
outer reaches of the galaxy, ∼ 10 − 35 kpc to the NE (all
distances are given in projection on the sky). This system
is spatially coincident with the diffuse radio emission we
described in Neff, Eilek, and Owen (2015; Paper 1). We
Electronic address: [email protected]
National Radio Astronomy Observatory is a facility of the
National Science Foundation operated under cooperative agreement
by Associated Universities Inc.
† Adjunct Astronomer at the National Radio Astronomy Observatory.
∗ The
suggest in this paper that both structures are causally
related to an energy-carrying flow in the region.
As in Paper 1, our focus here is the inner ∼ 50 kpc of
the Cen A / NGC 5128 system. We present new GALEX
observations and consider what they, together with new
radio data (Paper 1), reveal about the astrophysics of this
region. In the rest of this section we present an introduction to key aspects of the system. In Table 1 we list
regions and features discussed throughout the paper, and
give approximate distances from the galaxy’s center.
1.1. Large and small-scales
The outer radio lobes — The radio source Cen A was first
identified with the galaxy NGC 5128 by the radio emission from its Outer Lobes (Bolton et al. 1949), which
extend ∼ 600 kpc end-to-end on the sky (e.g., Junkes
et al. 1993; Feain et al. 2011). The outer radio lobes are
also detected as extended γ-ray sources (Abdo et al. 2010),
thought to result from inverse Compton scattering of Cosmic Microwave Background photons. Eilek (2014) argues
that the dynamical age of Cen A’s outer lobes is on the
order of ∼ 1 Gyr (assuming the lobes lie at some finite angle to the sky plane), but that they must have been last
energized no more than ∼ 30 Myr ago to keep shining in
radio and γ-rays. Any energy supplied to the Outer Lobes
by the inner galaxy must, of course, have moved through
the middle (∼ 50-kpc scale) regions which we study in
this paper.
Mon. Not. R. Astron. Soc. 000, 1–?? (2014)
Printed 19 February 2015
(MN LATEX style file v2.2)
On the absence of nova shells
arXiv:1502.05230v1 [astro-ph.SR] 18 Feb 2015
L. Schmidtobreick1, M. Shara2, C. Tappert3, A. Bayo3,4, A. Ederoclite5
1 European
2
3
4
5
Southern Observatory, Casilla 19001, Santiago 19, Chile
Department of Astrophysics, American Museum of Natural History, Central Park West and 79th Street, NY 10024-5192, USA
Instituto de F´ısica y Astronom´ıa, Universidad de Valpara´ıso, Avda. Gran Breta˜
na 1111, Valpara´ıso, Chile
Max-Planck-Institut f¨
ur Astronomie, K¨
onigstuhl 17, D–69117 Heidelberg, Germany
Centro de Estudios de F´ısica del Cosmos de Arag´
on, Plaza San Juan 1, Planta 2, Teruel, E44001, Spain
xxxx
ABSTRACT
We present our wide field Hα+N[II] observations of 15 cataclysmic variables to search
for remnant nova shells. Such shells have been found around other cataclysmic variables that were hitherto not known as novae. Our candidates were selected as objects
in the period regime of high-mass transfer systems that experience - at least occasionally - low mass transfer rates. The fact that we find no indication of a nova shell in
any of these systems allows us to set a lower limit of 13000 years to the recurrence
time of these objects.
Key words: stars: novae, cataclysmic variables
1
INTRODUCTION
A nova eruption in a cataclysmic variable (CV) is a thermonuclear explosion on the surface of the white-dwarf primary once it has accreted a critical mass from its late-type
companion. During the nova eruption, material is ejected
into the interstellar medium, forming an expanding shell
around the CV which can be observed once its angular size
is sufficiently large to be resolvable from the inner binary
(see e.g. Gill & O’Brien 1998). The total ejected mass is estimated between 10−6 and 10−4 M⊙ (Bode & Evans 2008,
other references therein).
In-between nova eruptions the binary is supposed to appear as a ”normal” CV, i.e. its behaviour is dominated by its
current mass-transfer rate and the magnetic field strength
of the white dwarf (Vogt 1989). However, the model of
Shara et al. (1986) predicts that the nova eruption strongly
˙ . After an initial phase of enaffects the mass-transfer rate M
hanced M˙ due to the irradiation of the secondary star by the
˙ is superuption-heated white dwarf (Kovetz et al. 1988), M
posed to decline over the next centuries by several orders of
magnitude (Livio & Shara 1987), possibly even to the point
that M˙ ∼ 0 due to the secondary star losing contact to its
Roche lobe, a scenario that has been termed ”hibernation”.
The recent discovery of ancient nova shells around two lowmass transfer systems, i.e. Z Cam (Shara et al. 2007) and
AT Cnc (Shara et al. 2012) provides strong support for the
idea of the nova - dwarf nova cycle originally proposed by
Vogt (1982). However, even the existence of dwarf novae that
have experienced nova outbursts in the past does not prove
the hibernation model. Individual such cases could equally
c 2014 RAS
be explained by the fact that all CVs can explode as a nova
once sufficient material has been accumulated on the surface
of the white dwarf. Only a careful population study of novae
in comparison to high- and low-mass transfer systems can
yield the answer to this long-debated question.
Most observed old novae do actually show a very
high mass transfer rate (see e.g. Schmidtobreick et al. 2005;
Tappert et al. 2012, 2014) This is expected as the recurrence
time, i.e. the time between two nova outbursts, is supposed
to be smaller for high mass transfer systems and so they are
more likely to be observed during a nova eruption. Still, this
also means that most of the observed old novae are not in
’hibernation’ which could be due to the time-scales and the
relatively short time that has passed after the nova eruption
(less than 100 years in most cases).
The currently known period distribution of novae
shows a significant peak at 3–4 h (Tappert et al. 2013a)
in rough agreement with the theoretical calculations by
Townsley & Bildsten (2005). If one considers all types of
CVs, this range between 3 and 4 h is dominated by
systems with very high M˙ (Townsley & G¨
ansicke 2009;
Schmidtobreick 2013; Rodr´ıguez-Gil et al. 2007). However,
there does exist a small population of dwarf novae in this period range (Ritter & Kolb 2003). Also three old novae that
have been re-discovered only recently (Schmidtobreick et al.
2005; Tappert et al. 2012) present optical spectra that are
much more akin to dwarf nova spectra than to nova-like
stars. At least one of them shows also evidence for dwarfnova like outburst behaviour (V728 Sco; Tappert et al.
˙ CVs in this range be
2013b). Could the presence of low M
a consequence of a previous nova eruption as predicted by
c ESO 2015
Astronomy & Astrophysics manuscript no. hernan˙obispo˙final
February 19, 2015
Analysis of combined radial velocities and activity of BD+20 1790:
evidence supporting the existence of a planetary companion ⋆ , ⋆⋆
M. Hern´an-Obispo1 , M. Tuomi2,3 , M.C. G´alvez-Ortiz4,2 , A. Golovin5 , J. R. Barnes6 , H. R. A. Jones2 , S. R. Kane7 , D.
Pinfield2 , J. S. Jenkins8 , P. Petit9,10 , G. Anglada-Escud´e11 , S. C. Marsden12 , S. Catal´an13 , S. V. Jeffers11 , E. de Castro1 ,
M. Cornide1 , A. Garc´es14 , M. I. Jones15 , N. Gorlova16 , and M. Andreev17
arXiv:1502.05223v1 [astro-ph.EP] 18 Feb 2015
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
Dpto. de Astrof´ısica y Ciencias de la Atm´osfera, Facultad de F´ısica, Universidad Complutense de Madrid, Avda. Complutense s/n,
E-28040, Madrid, Spain
e-mail: [email protected]
Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Hertfordshire AL10 9AB, UK
University of Turku, Tuorla Observatory, Deparment of Physics and Astronomy, V¨ais¨al¨antie 20, Fl-21500, Piikki¨o, Finland
Centro de Astrobiolog´ıa (CSIC-INTA), Ctra. Ajalvir km 4, E-28850 Torrej´on de Ardoz, Madrid, Spain
Main Astronomical Observatory of National Academy of Sciences of Ukraine, Zabolotnogo str., 27, Kiev, 03680, Ukraine
Department of Physical Sciences, The Open University, Walton Hall, Milton Keynes, MK7 6AA, U. K.
Department of Physics & Astronomy, San Francisco State University, 1600 Holloway Avenue, San Francisco, CA 94132, USA
Departmento de Astronom´ıa, Universidad de Chile, Camino del Observatorio 1515, Las Condes, Santiago, Chile, Casilla 36-D
Universit´e de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Plan´etologie, Toulouse, France
CNRS, Institut de Recherche en Astrophysique et Plan´etologie, 14 Avenue Edouard Belin, F-31400 Toulouse, France
Universit¨at G¨ottingen, Institut f¨ur Astrophysik, Friedrich-Hund-Platz 1, 37077 G¨ottingen, Germany
Computational Engineering and Science Research Centre, University of Southern Queensland, Toowoomba, 4350, Australia
Department of Physics, University of Warwick, Coventry CV4 7AL, UK
Institut de Ci´ences de l’Espai (IEEC-CSIC), Facultad de Ci´encies, Campus UAB, 08193, Bellaterra, Spain
Department of Electrical Engineering and Center of Astro-Engineering UC, Pontificia Universidad Cat´olica de Chile, Av. Vicu˜na
Mackenna 4860, 782-0436 Macul, Santiago, Chile
Institute of Astronomy, Katholieke Universiteit Leuven, Celestijnenlaan 200D BUS 2401, 3001 Leuven, Belgium
Terskol Branch of Institute of Astronomy RAS, Kabardino-Balkaria Republic, 361605 Terskol, Russia
recieved – ; accepted –
ABSTRACT
Context. In a previous paper we reported a planetary companion to the young and very active K5Ve star BD+20 1790. We found
that this star has a high level of stellar activity (logR′HK =-3.7) that manifests in a plethora of phenomena (starspots, prominences,
plages, large flares). Based on a careful study of these activity features and a deep discussion and analysis of the effects of the stellar
activity on the radial velocity measurements, we demonstrated that the presence of a planet provided the best explanation for the radial
velocity variations and all the peculiarities of this star. The orbital solution resulted in a close-in massive planet with a period of 7.78
days. However, a paper by Figueira et al. (2010) questioned the evidence for the planetary companion.
Aims. This paper aims to more rigorously assess the nature of the radial velocity measurements with an expanded dataset and new
methods of analysis.
Methods. We have employed Bayesian methods to simultaneously analyse the radial velocity and activity measurements based on a
combined dataset that includes new and previously published observations.
Results. We conclude that the Bayesian analysis and the new activity study support the presence of a planetary companion to BD+20
1790. A new orbital solution is presented, after removing the two main contributions of stellar jitter, one that varies with the photometric period (2.8 days) and another that varies with the synodic period of the star-planet system (4.36 days). We present a new
method to determine these jitter components, considering them as second and third signals in the system. A discussion on possible
star-planet-interaction is included, based on the Bayesian analysis of the activity indices, which indicates that they modulate with the
synodic period. We propose two different sources for flare events in this system: one related to the geometry of the system and the
relative movement of the star and planet, and a second one purely stochastic source that is related to the evolution of stellar surface
active regions. Also, we observe for the first time the magnetic field of the star, from spectropolarimetric data.
Key words. methods: statistical – techniques: radial velocities – stars: activity – planetary systems – stars: individual (BD+20 1790)
1. Introduction
The quest for planets around other stars has become one of the
most successful and productive fields in Astronomy. Since the
surprising discovery of 51 Peg b by Mayor & Queloz in 1995,
the ever increasing rate of newly reported planets 1 has revealed
a large variety of new and strange worlds, in some cases very
different from the ones in our Solar System.
⋆
See Acknowledgements
Tables of the photometry and of the RV are only available at the
CDS.
⋆⋆
1
exoplanet.eu
1
Mon. Not. R. Astron. Soc. 000, 1–16 (2015)
Printed February 19, 2015
(MN LATEX style file v2.2)
arXiv:1502.05221v1 [astro-ph.GA] 18 Feb 2015
The IGIMF and other IMFs in dSphs: the case of
Sagittarius
F. Vincenzo1,2⋆ , F. Matteucci1,2,3, S. Recchi4, F. Calura5, A. McWilliam6,
G.
A. Lanfranchi7
1
Dipartimento di Fisica, Sezione di Astronomia, Universit`
a di Trieste, via G.B. Tiepolo 11, 34100, Trieste, Italy
Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, 34100, Trieste, Italy
3 INFN, Sezione di Trieste, Via Valerio 2, 34100, Trieste, Italy
4 Department of Astrophysics, Vienna University, T¨
urkenschanzstrasse 17, A-1180, Vienna, Austria
5 INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, I-40127, Bologna, Italy
6 The Observatories of the Carnegie Institute of Washington, 813 Santa Barbara Street, Pasadena, CA 91101, USA
7 N´
ucleo de Astrof´ısica Te´
orica, Universidade Cruzeiro do Sul, Rua Galv˜
ao Bueno 868, 01506-000, S˜
ao Paulo, Brazil
2 INAF,
Accepted 2015 February 18. Received 2015 February 16; in original form 2014 October 20
ABSTRACT
We have studied the effects of various initial mass functions (IMFs) on the chemical evolution of the Sagittarius dwarf galaxy (Sgr). In particular, we tested the effects
of the integrated galactic initial mass function (IGIMF) on various predicted abundance patterns. The IGIMF depends on the star formation rate and metallicity and
predicts less massive stars in a regime of low star formation, as it is the case in dwarf
spheroidals. We adopted a detailed chemical evolution model following the evolution
of α-elements, Fe and Eu, and assuming the currently best set of stellar yields. We also
explored different yield prescriptions for the Eu, including production from neutron
star mergers. Although the uncertainties still present in the stellar yields and data prevent us from drawing firm conclusions, our results suggest that the IGIMF applied to
Sgr predicts lower [α/Fe] ratios than classical IMFs and lower [hydrostatic/explosive]
α-element ratios, in qualitative agreement with observations. In our model, the observed high [Eu/O] ratios in Sgr is due to reduced O production, resulting from the
IGIMF mass cutoff of the massive oxygen-producing stars, as well as to the Eu yield
produced in neutron star mergers, a more promising site than core-collapse supernovae,
although many uncertainties are still present in the Eu nucleosynthesis. We find that a
model, similar to our previous calculations, based on the late addition of iron from the
Type Ia supernova time-delay (necessary to reproduce the shape of [X/Fe] vs. [Fe/H]
relations) but also including the reduction of massive stars due to the IGIMF, better
reproduces the observed abundance ratios in Sgr than models without the IGIMF.
Key words: stars: abundances - galaxies: abundances - galaxies: dwarf - galaxies:
evolution - galaxies: formation - Local Group.
1
INTRODUCTION
The Sagittarius (Sgr) dwarf galaxy was the last classical
dwarf spheroidal (dSph) discovered before the advent of
the Sloan Digital Sky Survey (SDSS). Its discovery was
made by Ibata et al. (1994) and it was identified while performing a spectroscopic radial velocity survey of the Galactic bulge stars (Ibata et al. 1997). Its heliocentric distance
(D⊙ = 26 ± 2 kpc, from Simon et al. 2011) makes it the
second closest known satellite galaxy of the Milky Way
⋆
E-mail: [email protected]
c 2015 RAS
(MW) and, because of the strong tidal interaction suffered by the Sgr dSph during its orbit, it has left behind a
well-known stellar stream (Ibata et al. 2001; Majewski et al.
2003; Belokurov et al. 2006), whose chemical characteristics
have been recently studied and compared with the ones of
the Sgr main body and other dSph galaxies by de Boer et al.
(2014). The Sgr dwarf galaxy has been classified as a dwarf
spheroidal because of its very low central surface brightness
(µV = 25.2 ± 0.3 mag arcsec−2 , from Majewski et al. 2003),
its very small total amount of gas (MHI,obs ∼ 104 M⊙ , from
McConnachie 2012) and because of the age and metallicity
of its main stellar population, which dates back to the age of
IFIC/15-10
Does Current Data Prefer a Non-minimally Coupled Inflaton?
Lotfi Boubekeur,1, 2 Elena Giusarma,3 Olga Mena,1 and H´ector Ram´ırez1
1
arXiv:1502.05193v1 [astro-ph.CO] 18 Feb 2015
Instituto de F´ısica Corpuscular (IFIC), CSIC-Universitat de Valencia,
Apartado de Correos 22085, E-46071, Spain.
2
Laboratoire de Physique Math´ematique et Subatomique (LPMS)
Universit´e de Constantine I, Constantine 25000, Algeria.
3
Physics Department and INFN, Universit`
a di Roma “La Sapienza”, Ple Aldo Moro 2, 00185, Rome, Italy
We examine the impact of a non-minimal coupling of the inflaton to the Ricci scalar, 12 ξRφ2 ,
on the inflationary predictions. Such a non-minimal coupling is expected to be present in the
inflaton Lagrangian on fairly general grounds. As a case study, we focus on the simplest inflationary
model governed by the potential V ∝ φ2 , using the latest combined 2015 analysis of Planck and
BICEP2/Keck Array. We find that, for all the data combinations used in this study, a small
positive value of the coupling ξ is favoured at the 2σ level. When considering the cross-correlation
polarization spectra from BICEP2/Keck Array and Planck, a value of r > 0 is found at 95% CL.
PACS numbers: 98.70.Vc, 98.80.Cq, 98.80.Bp
Motivations.— Inflation provides the most theoretically attractive and observationally successful cosmological scenario able to generate the initial conditions of
our universe, while solving the standard cosmological
problems. Despite this remarkable success, the inflationary paradigm is still lacking firm observational confirmation. The picture that emerges from the latest data from
Planck, including also the joint analysis of B-mode polarization measurements from the BICEP2 collaboration [1–
4], is compatible with the inflationary paradigm. According to these observations, structure grows from Gaussian
and adiabatic primordial perturbations. From the theoretical viewpoint, this picture is usually understood as
the dynamics of a single new scalar degree of freedom, the
inflaton, minimally coupled to Einstein gravity. However,
the inflaton φ is expected to have a non-minimal coupling
to the Ricci scalar through the operator 12 ξRφ2 , where ξ
is a dimensionless coupling. Indeed, successful reheating
requires that the inflaton is coupled to the light degrees
of freedom. Such couplings, though weak, will induce
a non-trivial running for ξ. Thus, even starting from
a vanishing value of ξ (away from the conformal fixed
point ξ = −1/6) at some energy scale, a non-trivial nonminimal coupling will be generated radiatively at some
other scale (see e.g. Ref. [5]). Therefore, it is important
to study the impact of such a coupling on the inflationary
predictions, especially in view of the latest Planck 2015
data.
Generically, for successful inflation, the inflaton should
be very weakly coupled∗ . It follows that the magnitude
of ξ is expected to be small. Yet, even with such a suppressed coupling, the inflationary predictions are significantly altered [6–13]. For instance, and as we will see, a
small and positive ξ can enlarge considerably the space of
∗
This requirement is also dictated by the non-detection of large
primordial non-Gaussianities [3] and the soft breaking of the shift
symmetry φ → φ + c, necessary to protect the flatness of the
potential.
phenomenologically acceptable scenarios (See also [14]).
In this letter, we will focus on the simplest inflationary scenario with a potential V ∝ φ2 [15], and a nonzero non-minimal coupling. According to the very recent
Planck 2015 full mission results, the minimally-coupled
version of this scenario (i.e. ξ = 0) is ruled out at more
than 99% confidence level [2, 4], for 50 e-folds of inflation. Thus, before discarding it definitely from the range
of theoretical possibilities, it is worthwhile to explore this
scenario in all generality, given that, as explained earlier,
the presence of non-minimal couplings in the inflaton Lagrangian is quite generic.
Non-minimally coupled Inflaton.— The dynamics of
a non-minimally coupled scalar field φ with a potential
U (φ) is governed, in the Jordan frame, by the following
action†
2
Z
√
MP
ξ
1
S=
d4 x −g
R + Rφ2 − (∂φ)2 − U (φ) ,
2
2
2
(1)
where indices are contracted with the metric gµν , defined
as ds2 = − dt2 + a2 (t) dx2 . Inflation can be conveniently
studied in the Einstein frame, after performing a conforE
mal transformation gµν
= Ω(φ)gµν , with Ω ≡ 1+ξφ2 /MP2
and canonically-normalizing the scalar field. Up to a total derivative, the action takes the familiar form
2
Z
√
MP
1 µν
RE − gE
∂µ ϕ∂ν ϕ − V [φ (ϕ)] ,
S = d4 x −gE
2
2
(2)
where now ϕ is the canonically-normalized inflaton, related to the the original non-minimally coupled scalar
field φ through
2
2
dϕ
1
3 2 Ω0
= + MP
.
(3)
dφ
Ω 2
Ω
†
√
As usual, MP = 1/ 8πGN ' 2.43 × 1018 GeV is the reduced
Planck mass.
arXiv:1502.05190v1 [astro-ph.SR] 18 Feb 2015
DYNAMICS OF ELECTRIC CURRENTS, MAGNETIC FIELD
TOPOLOGY AND HELIOSEISMIC RESPONSE OF A SOLAR
FLARE
I.N. Sharykin1,2 and A.G. Kosovichev1,3,4
Big Bear Solar Observatory, New Jersey Institute of Technology, Big Bear City, CA 92314,
U.S.A
Received
;
accepted
1
Big Bear Solar Observatory
2
Space Research Institute (IKI) of the Russian Academy of Science
3
Stanford University
4
NASA Ames Research Center
–2–
ABSTRACT
The solar flare on July 30, 2011 was of a modest X-ray class (M9.3), but
it made a strong photospheric impact and produced a “sunquake,” observed
with the Helioseismic and Magnetic Imager (HMI) on NASA’s Solar Dynamics
Observatory (SDO). In addition to the helioseismic waves (also observed with
the SDO/AIA instrument), the flare caused a large expanding area of white-light
emission and was accompanied by substantial restructuring of magnetic fields,
leading to the rapid formation of a sunspot structure in the flare region. The
flare produced no significant hard X-ray emission and no coronal mass ejection.
This indicates that the flare energy release was mostly confined to the lower
atmosphere. The absence of significant coronal mass ejection rules out magnetic
rope eruption as a mechanism of helioseismic waves. We discuss the connectivity
of the flare energy release with the electric currents dynamics and show the
potential importance of high-speed plasma flows in the lower solar atmosphere
during the flare energy release.
Subject headings: Sun: flares; Sun: helioseismology; Sun: sunspots; Sun: X-rays
arXiv:1502.05182v1 [astro-ph.SR] 18 Feb 2015
Binarity and the abundance discrepancy problem
in planetary nebulae
Romano L.M. Corradi, Jorge Garc´ıa-Rojas, David Jones, and Pablo Rodr´ıguez–Gil
Instituto de Astrof´ısica de Canarias, E-38200 La Laguna, Tenerife, Spain
Departamento de Astrof´ısica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain
ABSTRACT
The discrepancy between abundances computed using optical recombination lines (ORLs)
and collisionally excited lines (CELs) is a major unresolved problem in nebular astrophysics. We
show here that the largest abundance discrepancies are reached in planetary nebulae with close
binary central stars. This is illustrated by deep spectroscopy of three nebulae with a post commonenvelope (CE) binary star. Abell 46 and Ou5 have O2+ /H+ abundance discrepancy factors larger
than 50, and as high as 300 in the inner regions of Abell 46. Abell 63 has a smaller discrepancy
factor around 10, but still above the typical values in ionized nebulae. Our spectroscopic analysis
supports previous conclusions that, in addition to “standard” hot (Te ∼104 K) gas, a colder
(Te ∼103 K) ionized component that is highly enriched in heavy elements also exists. These
nebulae have low ionized masses, between 10−3 and 10−1 M⊙ depending on the adopted electron
densities and temperatures. Since the much more massive red-giant envelope is expected to be
entirely ejected in the CE phase, the currently observed nebulae would be produced much later,
in post–CE mass loss episodes when the envelope has already dispersed. These observations add
constraints to the abundance discrepancy problem. Possible explanations are revised. Some are
naturally linked to binarity, such as for instance high-metallicity nova ejecta, but it is difficult
at this stage to depict an evolutionary scenario consistent with all the observed properties. The
hypothesis that these nebulae are the result of tidal destruction, accretion and ejection of Jupiterlike planets is also introduced.
Subject headings: planetary nebulae: individual (A 46, A 63, Ou5) – ISM: abundances - binaries: close
– novae, cataclysmic variables – planet–star interactions
1.
recombination lines (ORLs) provide abundance
values that are systematically larger than those
obtained using collisionally excited lines (CELs).
This a long-standing problem in nebular astrophysics, and has obvious implications on the measurement of the chemical content of the Universe,
often done using CELs from emission regions and
the ISM. The abundance discrepancy factor (adf )
between ORLs and CELs is usually between 1.5
and 3 (see e.g. Garc´ıa–Rojas & Esteban 2007; Liu
2012; Esteban et al. 2014), but in PNe it has a significant tail extending to much larger values. The
nebula with the largest known adf (∼70, Liu et al.
2006) is Hf 2–2, which has a close binary central
star (Lutz et al. 1998).
Introduction
This work deals with two apparently unrelated,
main topics in the study of planetary nebulae
(PNe).
The first one is the role of binary evolution,
which is the favoured explanation of the diverse
morphologies displayed by PNe (e.g. Soker 1997;
Balick & Frank 2002), but which may even be the
cause of their mere existence (Moe & De Marco
2001; Soker 2006).
The second issue is the so-called abundance discrepancy problem. It is well known (see e.g. Osterbrock & Ferland 2006) that in photoionized
nebulae – both PNe and H II regions – optical
1
5th Fermi Symposium : Nagoya, Japan : 20-24 Oct, 2014
A strong radio brightening at the jet base of M87 in the period of the
elevated TeV γ-ray state in 2012
K. Hadaa,b , M. Girolettia , M. Kinoc,d , G. Giovanninia,e , F. D’Ammandoa , C. C. Cheungf , M. Beilickeg ,
H. Nagaih , A. Doid , K. Akiyamab,i , M. Honmab,j , K. Niinumak , C. Casadiol , M. Orientia ,
l
´
H. Krawczynski7 , J. L. Gomez
, S. Sawada-Satohb , S. Koyamab,d,i , A. Cesarinim , S. Nakaharan and
M. A. Gurwello
a
INAF Istituto di Radioastronomia, via Gobetti 101, I-40129 Bologna, Italy
Mizusawa VLBI Observatory, National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan
c
Korea Astronomy and Space Science Institute, 776 Daedukdae-ro, Yusong, Daejon 305-348, Korea
d
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1 Yoshinodai, Chuo, Sagamihara 252-5210, Japan
e
Dipartimento di Fisica e Astronomia, Universita` di Bologna, via Ranzani 1, I-40127 Bologna, Italy
f
Space Science Division, Naval Research Laboratory, Washington, DC 20375, USA
g
Physics Department and McDonnell Center for the Space Sciences, Washington University, St.
Louis, MO 63130, USA
h
National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan
i
Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
j
Department of Astronomical Science, The Graduate University for Advanced Studies (SOKENDAI), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
k
Graduate School of Science and Engineering, Yamaguchi
University, 1677-1 Yoshida, Yamaguchi, 753-8512, Japan
l
Instituto de Astrofisica de Andalucia, CSIC, Apartado 3004, 18080 Granada, Spain
m
Department of Physics, University of Trento, I38050, Povo, Trento, Italy
n
Faculty of Science, Kagoshima University, 1-21-35 Korimoto, Kagoshima, Kagoshima 890-0065, Japan
o
Harvard-Smithsonian Center for Astrophysics, Cambridge MA 02138 USA
arXiv:1502.05177v1 [astro-ph.HE] 18 Feb 2015
b
The nearby radio galaxy M87 offers a unique opportunity for exploring the connection between γ-ray production
and jet formation at an unprecedented linear resolution. However, the origin and location of the γ-rays in this
source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and
the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of γ-ray production.
Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI
Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to
October 2012, together with contemporaneous high-energy γ-ray light curves obtained by the Fermi Large
Area Telescope. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS.
We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE
activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure in the radio band.
These results strongly suggest that the TeV γ-ray activity in 2012 originates in the jet base within 0.03 pc
(projected) from the central supermassive black hole.
1. Introduction
The nearby radio galaxy M87 accompanies one of
the best studied AGN jets. Its proximity (16.7 Mpc)
and brightness have enabled detailed studies of this
jet over decades through radio, optical and to Xray at tens of parsec scale resolutions. Furthermore, the inferred very massive black hole (MBH ≃
(3 − 6) × 109 M⊙ ) yields a linear resolution down to
1 milliarcsecond (mas) = 0.08 pc = 140 Schwarzschild
radii (Rs ) (for MBH = 6×109 M⊙ ), making this source
an ideal case to probe the relativistic-jet formation
at an unprecedented compact scale with Very-LongBaseline-Interferometer (VLBI) observations (e.g., Ly
et al. 2007; Kovalev et al. 2007; Hada et al. 2011;
Asada & Nakamura 2012; Doeleman et al. 2012;
Hada et al. 2013). M87 is now widely known to
eConf C141020.1
show γ-ray emission up to the very-high-energy (VHE;
E > 100 GeV) regime, where this source often exhibits active ﬂaring episodes. The location and the
physical processes of such emission have been a matter of debate over the past years, and there are two
candidate sites which can be responsible for the VHE
γ-ray production. One is a very active knot HST-1
which is located at more than 100 pc from the nucleus (Stawarz et al. 2006; Cheung et al. 2007; Harris
et al. 2009). This argument is based on the famous
VHE ﬂare event in 2005, where HST-1 underwent a
large radio-to-X-ray outburst jointly with a VHE ﬂare.
In contrast, the other candidate is the core/jet base,
which is very close to the central black hole. This
argument is based on the VHE event in 2008, where
the core/VHE showed a remarkable correlation in the
light curves. There was another VHE event in 2010,
1
FAST RADIO BURSTS: COLLISIONS BETWEEN NEUTRON
STARS AND ASTEROIDS/COMETS
arXiv:1502.05171v1 [astro-ph.HE] 18 Feb 2015
J. J. Geng1,2 , and Y. F. Huang1,2
ABSTRACT
Fast radio bursts (FRBs) are newly discovered radio transient sources. Their
high dispersion measures indicate a cosmological origin. But due to the lack of
observational data in other wavelengths, their progenitors still remain unclear.
Here we suggest the collisions between neutron stars and asteroids/comets as
a promising mechanism for FRBs. During the impact process, a hot plasma
fireball will form after the material of the small body penetrates into the neutron
star surface. The ionized matter inside the fireball will then expand along the
magnetic field lines. Coherent radiation from the thin shell at the top of the
fireball will account for the observed FRBs. Our scenario can reasonably explain
the main features of FRBs, such as their durations, luminosities, and the event
rate. We argue that for a single neutron star, FRBs are not likely to happen
repeatedly in a forseeable time span since such impacts are of low probability.
We predict that faint X-ray afterglows should be associated with FRBs, which
may be detected by future rapid follow-up observations of FRBs.
Subject headings: pulsars: general — radio continuum: general — stars: neutron
— minor planets, asteroids: general
1.
INTRODUCTION
Recently, the discovery of a number of fast radio bursts (FRBs) was reported (Lorimer
et al. 2007; Thornton et al. 2013). Typically, they are single radio pulses with flux densities
Sν ∼ a few Jy and durations δt ∼ a few ms at frequency νFRB ∼ 1 GHz. There have
been no counterparts observed in other wavelengths till now (Petroff et al. 2014), maybe
due to the lack of rapid, multiwavelength follow-up after the bursts. According to their
1
2
School of Astronomy and Space Science, Nanjing University, Nanjing 210046, China; [email protected]
Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education,
Nanjing 210046, China
February 19, 2015, accepted for publication in The Astrophysical Journal Letters
Preprint typeset using LATEX style emulateapj v. 5/2/11
THE CONFINED X-CLASS FLARES OF SOLAR ACTIVE REGION 2192
J. K. Thalmann1 , Y. Su1 , M. Temmer1 and A. M. Veronig1
arXiv:1502.05157v1 [astro-ph.SR] 18 Feb 2015
February 19, 2015, accepted for publication in The Astrophysical Journal Letters
ABSTRACT
The unusually large NOAA active region 2192, observed in October 2014, was outstanding in its productivity of major two-ribbon flares without coronal mass ejections. On a large scale, a predominantly
north-south oriented magnetic system of arcade fields served as a strong, also lateral, confinement for
a series of large two-ribbon flares originating from the core of the active region. The large initial separation of the flare ribbons, together with an almost absent growth in ribbon separation, suggests a
confined reconnection site high up in the corona. Based on a detailed analysis of the confined X1.6 flare
on October 22, we show how exceptional the flaring of this active region was. We provide evidence for
repeated energy release, indicating that the same magnetic field structures were repeatedly involved
in magnetic reconnection. We find that a large number of electrons was accelerated to non-thermal
energies, revealing a steep power law spectrum, but that only a small fraction was accelerated to high
energies. The total non-thermal energy in electrons derived (on the order of 1025 J) is considerably
higher than that in eruptive flares of class X1, and corresponds to about 10% of the excess magnetic
energy present in the active-region corona.
Subject headings: Sun: photosphere — Sun: atmosphere — Sun: magnetic topology — Sun: activity
— Sun: flares — Sun: X-rays, gamma rays
1. INTRODUCTION
Coronal mass ejections (CMEs) and flares are interpreted to be different manifestations of a sudden instability and the associated release of magnetic energy in the
solar corona. In general, they can occur independently of
each other. Their association rate, however, is strongly
increasing with the strength of the event. As can be inferred from Figure 1 of Yashiro et al. (2006), in about
10%, 40%, and 75% of GOES class C1-, M1-, and X1flares, respectively, a CME association is found. Flares
≥X2.5 have an association rate >90%. Sometimes, however, the Sun shows striking deviations from this trend.
On 2014 Oct 17, active region (AR) NOAA 2192 appeared on the east limb of the Sun and developed into
the largest AR since NOAA 6368 in Nov 1990. In particular, the large size of NOAA 2192 was unexpected, as
it occurred in unusually weak solar cycle 24. During its
passage across the visible solar disk, between Oct 17 and
30, it produced six X- and 30 M-class flares, as well as
numerous smaller events. The GOES soft X-ray (SXR)
flux of the six largest flares peaked on Oct 19 05:03 UT
(X1.1), Oct 22 14:28 UT (X1.6), Oct 24 22:41 UT (X3.1),
Oct 25 17:08 UT (X1.0), Oct 26 10:56 UT (X2.0), and
Oct 27 14:47 UT (X2.0). The highly exceptional aspect
of the flaring activity was the lack of eruptive events:
none of the X-flares was accompanied by a CME.
In this letter, we investigate NOAA 2192 in the period
Oct 22–24, regarding its productivity of a series of large
(≥M5) though confined flares and a single eruptive M4.0
flare. During this period, the AR was located within
roughly 25◦ from disk center so that foreshortening effects were minimal. Additionally, we analyze in detail
the X1.6 flare on Oct 22 which, in contrast to the other
X-flares during the considered period, was well covered
[email protected]
1 Institute of Physics/IGAM, University of Graz, Universit¨
atsplatz 5/II, 8010 Graz, Austria
also by hard X-ray (HXR) data.
2. DATA AND METHODS
We use data from the Atmospheric Imaging Assembly
(AIA; Lemen et al. 2012) on board the Solar Dynamics Observatory (SDO; Pesnell et al. 2012). In particular, 1700 ˚
A (sampling the photosphere at the temperature minimum) and 1600 ˚
A (picturing photospheric plus
transition region emission) filtergrams were used for the
analysis of flare ribbons. For the clear distinction of the
low-atmosphere imprint of the four events under study,
we use 1700 ˚
A data. In this way, we avoid a contamination of the signal due to ejected material during the
eruptive M4.0 flare (which, in the line-of-sight, overlaps
with the actual flare ribbons that we aim to track). For
detailed analysis of the flare ribbons and the energy deposited by non-thermal electrons during the X1.6 flare,
we use 1600 ˚
A images.
Short-term brightenings not related to flaring activity
were removed by applying a 3-min running-median filter
to the image sequences (at a 1-min cadence). These filtered images were used to track the location and time
evolution of flare pixels. For the identification of flare
pixels, we use the 99-percentile intensity of the entire
series of filtered images as a threshold for detection. Importantly, the 99-percentile determines only the brightest
pixels in a series of images in a particular wavelength due
to its definition based on the relative occurrence of intensity values. Effects of blooming and saturated pixels
around the flare peak time were minimized by demanding that a flare pixel has to be identified in at least five
consecutive images.
The above data sets are complemented by
SDO/Helioseismic and Magnetic Imager (HMI; Schou
et al. 2012) magnetic field data.
The large-scale
coronal magnetic field environment around NOAA 2192
is retrieved via the potential field source surface (PFSS)
Mon. Not. R. Astron. Soc. 000, 000–000 (0000)
Printed 19 February 2015
(MN LATEX style file v2.2)
Measuring the equation of state of the high−z intergalactic medium
using curvature statistics
Hamsa Padmanabhan1⋆ , R. Srianand1 †, T. Roy Choudhury2‡
1
arXiv:1502.05140v1 [astro-ph.CO] 18 Feb 2015
2
Inter-University Centre for Astronomy and Astrophysics, Pune 411007, India
National Centre for Radio Astrophysics, Tata Institute of Fundamental Research, Pune 411007, India
ABSTRACT
Using hydrodynamical simulations, we explore the use of the mean and percentiles of the
curvature distribution function to recover the equation of state of the high-z (2 < z < 4)
intergalactic medium (IGM). We find that the mean and percentiles of the absolute curvature
distribution exhibit tight correlation with the temperatures measured at respective characteris¯ i ’s at each redshift. Hence, they provide nearly independent probes of the
tic overdensities ∆
same underlying temperature-density distribution, and can in principle be used to simultaneously recover both parameters T0 and γ of the IGM effective equation of state. We quantify
the associated errors in the recovered parameters T0 and γ from the intrinsic scatter in the
characteristic overdensities and the uncertainties in the curvature measurement.
Key words: dark ages, reionization, first stars - intergalactic medium - quasars : absorption
lines
1 INTRODUCTION
The Lyman-α absorption lines seen in the spectra of highredshift quasars arise predominantly from the low to moderate overdensity intergalactic medium (IGM; Cen et al. 1994;
Zhang, Anninos & Norman 1995; Petitjean, Mueket & Kates
1995; Hernquist et al. 1996; Miralda-Escud´e et al. 1996). The
balance between photoionization heating and adiabatic cooling
in the IGM leads to a mean temperature-density relation that
is well-approximated (Hui & Gnedin 1997) by a power-law,
T = T0 ∆γ−1 for overdensities ∆ ≤ 10. The parameters of this
“equation of state” depend upon the reionization history of the
universe and their predicted values at different epochs vary in
different reionization scenarios, so long as either hydrogen or
helium reionization happened relatively recently. Otherwise, the
temperature-density relation is predominantly set by the shape of
the ionizing ultraviolet background and the adiabatic expansion
or collapse of large-scale structure (Hui & Gnedin 1997). Nonradiative mechanisms of energy injection like, e.g. blazar heating
(Chang, Broderick & Pfrommer 2012; Puchwein et al. 2012) may
also distort the T − ∆ relation. Hence, it is important to measure
the temperature-density relation at various redshifts in order to
constrain the epoch and evolution of reionization and the properties
of the ionizing and heating sources.
Detailed studies of different properties of the IGM have
contributed to understanding the thermal history of the IGM,
by using observed data together with the results of numeri-
⋆
Electronic address: [email protected]
† Electronic address: [email protected]
‡ Electronic address: [email protected]
c 0000 RAS
cal simulations. The methods include (a) using the Doppler
b-parameter-column density (b − N ) or b-distribution of
the Lyman-α absorption lines originating from the IGM
(Haehnelt & Steinmetz 1998; Schaye et al. 1999, 2000;
McDonald et al. 2001; Choudhury, Srianand & Padmanabhan
2001; Ricotti, Gnedin & Shull 2000; Bryan & Machacek 2000;
Bolton et al. 2010, 2012), (b) using the small-scale power spectrum of the Lyman-α forest (Theuns, Schaye & Haehnelt 2000;
Zaldarriaga, Hui & Tegmark 2001), and (c) using the wavelet
decomposition of the Lyman-α lines (Theuns & Zaroubi 2000;
Zaldarriaga 2002; Theuns et al. 2002a,b; Lidz et al. 2010). The
latter two methods do not require Voigt profile decomposition of
the spectral lines. Some evidence for a peak in the intergalactic
medium temperature around z ∼ 3, together with a nearly isothermal profile (i.e. γ = 1) signifying the end of He II reionization, has
been reported (Schaye et al. 2000; Theuns et al. 2002b; Lidz et al.
2010), however, e.g. McDonald et al. (2001) do not find evidence
supporting this claim. Interestingly, the parameters (T0 and γ)
derived using the above methods typically have large uncertainties
(& 30%).
Recently, the temperature of the IGM over redshifts 4.5 to
∼ 2.8 has been measured to a high precision (. 10%) by
Becker et al. (2011) (and extended upto z ∼ 1.5 by Boera et al.
2014) using the curvature statistic, which also avoids the fitting of individual spectral lines and can be used to detect
the additional heating effect in high-redshift quasar near-zones
(Padmanabhan, Choudhury & Srianand 2014, hereafter Paper I).
The curvature is normally sensitive to both T0 and γ, however, it
¯ it becomes more
is found that at a characteristic overdensity ∆,
¯ and fairly independent of γ. The value of ∆
¯ insensitive to T (∆)
creases with decreasing redshifts, going towards non-linear densi-
Astronomy & Astrophysics manuscript no. extQSO_paper_ArXiV
February 19, 2015
c
ESO
2015
Where is the fuzz? Undetected Lyman α nebulae around QSOs at
z∼2.3∗
Edmund Christian Herenz1 , Lutz Wisotzki1 , Martin Roth1 , and Friedrich Anders1
Leibniz-Institut für Astrophysik Potsdam (AIP)
An der Sternware 16, 14482 Potsdam e-mail: [email protected]
arXiv:1502.05132v1 [astro-ph.GA] 18 Feb 2015
Received . . . ; accepted: 17 February 2014
ABSTRACT
We observed a small sample of 5 radio-quiet QSOs with integral field spectroscopy to search for possible extended emission in the
Lyα line. We subtracted the QSO point sources using a simple PSF self-calibration technique that takes advantage of the simultaneous
availability of spatial and spectral information. In 4 of the 5 objects we find no significant traces of extended Lyα emission beyond
the contribution of the QSO nuclei itself, while in UM 247 there is evidence for a weak and spatially quite compact excess in the Lyα
line at several kpc outside the nucleus. For all objects in our sample we estimated detection limits for extended, smoothly distributed
Lyα emission by adding fake nebulosities into the datacubes and trying to recover them after PSF subtraction. Our observations are
consistent with other studies showing that giant Lyα nebulae such as those found recently around some quasars are very rare. Lyα
fuzz around typical radio-quiet QSOs is fainter, less extended and is therefore much harder to detect. The faintness of these structures
is consistent with the idea that radio-quiet QSOs typically reside in dark matter haloes of modest masses.
Key words. galaxies: quasars: emission lines – galaxies: high-redshift
galaxies: quasars: individual: UM 247, Q 0027+0103, Q 0256-0003, Q 0308+0129, Q 2243+0141,
1. Introduction
The circum-galactic medium (CGM), defined as gas outside of
the main stellar body of galaxies but still within the virial radii
of their dark-matter haloes, is of crucial importance in galaxy
evolution. It may act as a reservoir for fuelling star formation in
the galaxy, and it is also subject to feedback processes that expell material from the galaxy. If violent enough, this feedback
may in turn heat up the CGM and thus prevent it from contributing further to the formation of stars. Observing the CGM at high
redshifts will hence provide relevant insights about galaxy formation.
One observational approach to study the ‘cold’ gas phase
4
of the CGM (T <
∼ 10 K) at high redshift uses absorption
signatures against background sources. This has provided important statistical constraints on several properties of the CGM
(e.g. Adelberger et al. 2005; Hennawi et al. 2006; Steidel et al.
2010). However, the spatial distribution of the CGM in individual galaxy cannot be captured with this method.
An alternative approach is to map the CGM in the Lyman
α emission line. Several mechanisms have been identified that
should lead to Lyα emission from the CGM: Cooling of infalling
gravitationally heated gas (e.g., Haiman et al. 2000), cooling
following superwind-driven shocks (e.g., Taniguchi & Shioya
2000) and – possibly most important for our investigation – Lyα
fluorescence induced by exposure to UV radiation. While the
metagalactic UV background alone is predicted to produce only
a very faint glow that is probably out of reach for the current generation of optical instruments (e.g., Kollmeier et al. 2010), Lyα
∗
Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, operated jointly by the MaxPlanck Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC)
fluorescence caused by the much stronger UV radiation from a
QSO should boost the emission into the detectable regime (Rees
1988; Haiman & Rees 2001; Kollmeier et al. 2010). Searching
for Lyα signatures of the CGM around luminous QSOs is the
topic of the present study.
Haiman & Rees (2001) estimated that a z ∼ 3 QSO harbouring a 5 × 1011 M halo should be surrounded by Lyα fuzz
extending radially outwards to 200 –300 , at a surface brightness of
∼ 5 × 10−17 erg s−1 cm−2 arcsec−2 . At least in this model, Lyα
fuzz is predicted as a generic property of high-z QSOs, where
the surface brightness of this fuzz depends only on the mass of
the halo. More recent theoretical work suggests, however, that
a substantial fraction of gas is accreted within filamentary cold
streams (Dekel et al. 2009; Faucher-Giguère & Kereš 2011; Rosdahl & Blaizot 2012, and references therein). If these streams are
optically thick to ionizing radiation, they will develop a highly
ionized skin in the presence of a quasar. This skin then acts like
a mirror converting up to two thirds of the incident ionizing radiation into Lyα photons. In this scenario the expected surface
brightness then depends on the ionizing photon flux produced by
the quasar and the projected spatial configuration of the streams
(Kollmeier et al. 2010; Hennawi & Prochaska 2013).
Revealing extended Lyα structures around QSOs requires
a proper subtraction of the PSF-broadened nuclear component,
which will outshine – even under good seeing conditions – the
expected CGM signal close to the QSO. Observationally, this
makes the detection of circum-QSO Lyα fuzz much harder than
searching for Lyα ‘blobs’, now routinely found in large-area
narrow-band surveys (Steidel et al. 2000; Matsuda et al. 2004;
Saito et al. 2006; Ouchi et al. 2009; Matsuda et al. 2011; Erb
et al. 2011; Prescott et al. 2012, 2013). Many of these ‘blobs’
have no obvious central source of ionizing photons (e.g., NilsArticle number, page 1 of 12
To appear in ApJ Letters
Preprint typeset using LATEX style emulateapj v. 01/23/15
THE FIRST SCIENCE RESULTS FROM SPHERE: DISPROVING THE PREDICTED BROWN DWARF
AROUND V471 TAU*
A. Hardy1,6 , M.R. Schreiber1,6 , S.G. Parsons 1 , C. Caceres1,6 , G. Retamales1,6 Z. Wahhaj2 , D. Mawet2 , H.
Canovas1,6 , L. Cieza5,6 , T.R. Marsh3 , M.C.P. Bours3 , V.S. Dhillon,4 , A. Bayo1,6
arXiv:1502.05116v1 [astro-ph.EP] 18 Feb 2015
To appear in ApJ Letters
ABSTRACT
Variations of eclipse arrival times have recently been detected in several post common envelope binaries
consisting of a white dwarf and a main sequence companion star. The generally favoured explanation
for these timing variations is the gravitational pull of one or more circumbinary substellar objects
periodically moving the centre of mass of the host binary. Using the new extreme-AO instrument
SPHERE, we image the prototype eclipsing post-common envelope binary V471 Tau in search of
the brown dwarf that is believed to be responsible for variations in its eclipse arrival times. We
report that an unprecedented contrast of ∆mH = 12.1 at a separation of 260 mas was achieved, but
resulted in a non-detection. This implies that there is no brown dwarf present in the system unless
it is three magnitudes fainter than predicted by evolutionary track models, and provides damaging
evidence against the circumbinary interpretation of eclipse timing variations. In the case of V471 Tau,
a more consistent explanation is offered with the Applegate mechanism, in which these variations are
prescribed to changes in the quadrupole moment within the main-sequence star.
Keywords: stars: individual(V471 Tau) — planet-star interactions — binaries: eclipsing — brown
dwarfs — white dwarfs — binaries: close
1. INTRODUCTION
Circumbinary substellar objects, although anticipated
for a long time, have only recently been identified around
main-sequence binary stars (Doyle et al. 2011). Long
before this however, claims have been made for circumbinary substellar objects around close white dwarfmain sequence binaries. Initially consisting of a main
sequence binary with separations of order ∼1 AU
(Zorotovic & Schreiber 2013), these systems are believed
to have been through a spectacular phase of binary star
evolution to explain their current close separation. When
the more massive star of the binary evolves off the main
sequence, it fills its Roche-lobe and generates dynamically unstable mass transfer onto the secondary star. As
the time scale for dynamically unstable mass transfer
is much shorter than the thermal time scale of the secondary, the latter cannot adjust its structure fast enough
to incorporate the overflowing mass. Instead, a common
envelope of material forms around both the secondary
star and the core of the giant. Drag forces between
the envelope and the central binary then cause the envelope to be expelled at the expense of orbital energy and
angular momentum of the binary (e.g. Webbink 1984;
* Based on observations collected at the European Southern
Observatory, Chile, program ID 60.A-9355(A)
1 Departamento de F´
ısica y Astronom´ıa, Universidad Valpara´ıso, Avenida Gran Breta˜
na 1111, Valpara´ıso, Chile.
2 European Southern Observatory, Casilla 19001, Vitacura,
Santiago, Chile.
3 Department of Physics, Gibbet Hill Road, University of Warwick, Coventry, CV4 7AL, UK
4 Department of Physics and Astronomy, University of
Sheffield, Sheffield S3 7RH, UK
5 Nucleo de Astronomia, Universidad Diego Portales, Av.
Ej´
ercito 441, Santiago, Chile.
6 Millennium Nucleus ‘Protoplanetary Disks in ALMA Early
Science’, Universidad de Valpara´ıso, Avenida Gran Breta˜
na
1111, Valpara´ıso, Chile
Zorotovic et al. 2010; Ivanova et al. 2013). The emerging
white dwarf-main sequence binaries contain separations
of just a few solar radii, and are known as post common
envelope binaries (PCEBs) (Nebot G´omez-Mor´an et al.
2011).
Shortly after the discovery of the first PCEB it was
realised that it displays variations in its eclipse arrival times. Today, similar variations are seen in almost all eclipsing PCEBs with long enough coverage
(Parsons et al. 2010b; Zorotovic & Schreiber 2013), for
which the most common hypothesis is the presence of
a circumbinary object, typically a brown dwarf or multiple giant planets. In this scenario, the gravitational
pull of the circumbinary objects periodically move the
center of mass of the host binary stars, thereby changing the light travel time of the eclipse signal to Earth
(Irwin 1959). Indeed, the planetary model employed to
explain the eclipse timing variations (ETVs) seen in the
PCEB NN Ser (Beuermann et al. 2010) successfully predicted new eclipse arrival times (Beuermann et al. 2013;
Marsh et al. 2014), providing support to the circumbinary interpretation but raising questions regarding the
formation of these third objects. Zorotovic & Schreiber
(2013) favour a scenario in which the circumbinary objects form as a consequence of the common envelope evolution, in a so-called second generation scenario. This is
based on the finding that nearly all PCEBs with long
enough coverage show ETVs, yet only a small fraction of
main sequence binaries seem to host circumbinary substellar objects. Indeed, Schleicher & Dreizler (2014) were
able to develop a model in which a second generation
protoplanetary disk forms during common envelope evolution and produces giant planets through the disk instability model. In contrast, Bear & Soker (2014) prefer
the first generation scenario, in which the objects form
at a similar time to their main-sequence hosts, and survive the common-envelope phase. They claim that if a
STEP Group at USTC
http://space.ustc.edu.cn/dreams/
February 19, 2015, 1:22am
Investigating plasma motion of magnetic clouds at 1 AU through a
velocity-modified cylindrical force-free flux rope model
arXiv:1502.05112v1 [astro-ph.SR] 18 Feb 2015
Yuming Wang,1,2,∗ Zhenjun Zhou,1 Chenglong Shen,1,2 Rui Liu,1,3 and S. Wang1
1
CAS Key Laboratory of Geospace Environment, Department of Geophysics and Planetary Sciences, University of Science and
Technology of China, Hefei, Anhui 230026, China
2
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China,
Hefei, Anhui 230026, China
3
Mengcheng National Geophysical Observatory, School of Earth and Space Sciences, University of Science and Technology of China,
Hefei, China
∗
Corresponding Author, Contact: [email protected]
Contents
1 Introduction
2 Model and Events
2.1 Velocity-modified cylindrical force-free
rope model . . . . . . . . . . . . . . . .
2.1.1 Model description . . . . . . . .
2.1.2 Parameters . . . . . . . . . . . .
2.1.3 Evaluation of the goodness-of-fit
2.2 Events and model testing . . . . . . . .
that the expansion speed is correlated with the radial propagation speed and most MCs did not expand self-similarly at 1 AU. In
our statistics, about 62%/17% of MCs underwent a under/overexpansion at 1 AU and the expansion rate is about 0.6 on average.
Third, most interestingly, we find that a significant poloidal motion did exist in some MCs. Three speculations about the cause
of the poloidal motion are therefore proposed. These findings advance our understanding of the MC’s properties at 1 AU as well
as the dynamic evolution of CMEs from the Sun to interplanetary
space.
1
2
flux
. . .
. . .
. . .
. . .
. . .
3 Effects of velocity on the fitting results
2
2
3
4
5
1
7
4 Statistical properties of plasma motion
4.1 Linear propagating motion . . . . . . . . . . .
4.2 Expanding motion . . . . . . . . . . . . . . .
4.3 Poloidal motion . . . . . . . . . . . . . . . . .
11
11
11
13
5 conclusions and discussion
14
A Derivation of poloidal speed vϕ
16
B Additional tables
17
Introduction
Since first identified by Burlaga et al. in 1981, magnetic
clouds (MCs) have been studied extensively in the past
decades. They are large-scale organized magnetic structures
in interplanetary space, developed from coronal mass ejections (CMEs), and play an important role in understanding
the evolution of CMEs from the Sun to the heliosphere and
the associated geoeffectiveness.
The current knowledge of MCs are mostly from in-situ
one-dimensional observations, and various MC fitting models have been developed to reconstruct the global picture of
MCs in two or three dimensions. It is now believed that an
MC is a loop-like magnetic flux rope with two ends rooting
on the Sun [e.g., Burlaga et al., 1981; Larson et al., 1997;
Janvier et al., 2013]. The modeling efforts mainly focus on
two aspects. One is to reconstruct a realistic geometry and
magnetic field configuration. In past decades, MC fitting
models have been developed from cylindrically symmetrical
force-free flux ropes [e.g., Goldstein, 1983; Marubashi , 1986;
Burlaga, 1988; Lepping et al., 1990; Kumar and Rust , 1996]
gradually to asymmetrically cylindrical (non-)force-free flux
ropes [e.g., Mulligan and Russell , 2001; Hu and Sonnerup,
2002;
Hidalgo et al.,
2002a,b;
Cid et al.,
2002;
Vandas and Romashets , 2003] and torus-shaped flux ropes
[e.g., Romashets and Vandas , 2003; Marubashi and Lepping,
2007; Hidalgo and Nieves-Chinchilla, 2012]. Some comparisons among various MC fitting models could be found
Abstract. Magnetic clouds (MCs) are the interplanetary counterparts of coronal mass ejections (CMEs), and usually modeled
by a flux rope. By assuming the quasi-steady evolution and selfsimilar expansion, we introduce three types of global motion into
a cylindrical force-free flux rope model, and developed a new
velocity-modified model for MCs. The three types of the global
motion are the linear propagating motion away from the Sun,
the expanding and the poloidal motion with respect to the axis
of the MC. The model is applied to 72 MCs observed by Wind
spacecraft to investigate the properties of the plasma motion of
MCs. First, we find that some MCs had a significant propagation
velocity perpendicular to the radial direction, suggesting the direct evidence of the CME’s deflected propagation and/or rotation
in interplanetary space. Second, we confirm the previous results
1
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 01/23/15
DYNAMICAL STABILITY OF IMAGED PLANETARY SYSTEMS IN FORMATION:
APPLICATION TO HL TAU
D. Tamayo1,2,3 , A. H. M. J. Triaud1,3,4 , K. Menou1,4 and H. Rein1
arXiv:1502.05099v1 [astro-ph.EP] 18 Feb 2015
Draft version February 19, 2015
ABSTRACT
We present a general and simple framework for understanding the dynamical stability of planets
embedded in a protoplanetary nebula over typical disk lifetimes, and provide estimates for the maximum allowable planetary masses. We collect these easily evaluated dynamical constraints into a
workflow that can help guide the design and interpretation of new observational campaigns and numerical simulations of gap opening in such systems. We argue that the locations of resonances should
be significantly shifted from integer period ratios in massive disks like HL Tau, and that theoretical
uncertainties in the exact shift, together with observational errors, imply a large uncertainty in the
dynamical state and stability in such disks. This renders our results largely insensitive to an improved
determination of the gaps’ orbital radii, and presents an important barrier to using systems like HL
Tau as a proxy for the initial conditions following planet formation. An important observational
avenue to breaking this degeneracy is to search for eccentric gaps, which could implicate resonantly
interacting planets. Unfortunately, massive disks like HL Tau should induce swift pericenter precession that would smear out any such eccentric features of planetary origin. This strongly motivates
pushing toward more typical, less massive disks. For a nominal non-resonant model of the HL Tau
system with five planets, we find a maximum mass for the outer three bodies of approximately 2
Neptune masses. In a resonant configuration, these planets can reach at least the mass of Saturn.
The inner two planets’ masses are unconstrained by dynamical stability arguments.
Subject headings: celestial mechanics—planets and satellites: dynamical evolution and stability
1. INTRODUCTION
Given its wavelength coverage, sensitivity and resolution, the Atacama Large Millimeter/Submillimeter Array (ALMA) represents a unique observatory for studying not only the extent and masses of protoplanetary
disks, but also their structure. The latter provides the
possibility of detecting forming planets as they gravitationally sculpt the surrounding disk. From October
24-31 2014, as part of its science verification process,
the ALMA team pointed 25-30 antennas at the young
star HL Tau (J2000 04:31:38.45 +18:13:59.0), measuring
the continuum emission in Band 6 (211-275 GHz), with
baselines from 15 m to 15 km5 . As of the time of this
writing, no scientific publication has appeared analysing
these data; however, the ALMA project has publicly released an image of the system (see Fig. 1) that reveals
several concentric dust gaps, perhaps revealing forming
planets massive enough to sculpt the nearby disk material.
This striking discovery during the science verification
phase suggests that ALMA may find similar features in
several other young systems. In this paper, we [email protected]
1 Department of Physical & Environmental Sciences, University of Toronto at Scarborough, Toronto, Ontario M1C 1A4,
Canada
2 Canadian Institute for Theoretical Astrophysics, 60 St.
George St, University of Toronto, Toronto, Ontario M5S 3H8,
Canada
3 Centre for Planetary Sciences Fellow
4 Department of Astronomy & Astrophysics, University of
Toronto, Toronto, Ontario M5S 3H4, Canada
5
https://public.nrao.edu/news/pressreleases/
planet-formation-alma
fore generally investigate the constraints one can extract
on potential planetary masses from the requirement that
observed systems be dynamically stable. We then apply
these constraints to the HL Tau system using approximate gap separations from the publicly released image.
As argued in Appendix A, our analysis should be largely
insensitive to an improved determination of the orbital
radii.
We begin in Sec. 2 with a list of assumptions that we
make in our analysis. Section 3 then investigates stability in non-resonant systems and applies our thresholds
to numerical simulations of the HL Tau system, finding
that resonances may be important in shaping the system’s stability. Section 4 discusses the role of resonances,
and provides simple criteria for evaluating whether resonances must be considered for stability. Section 5 outlines a procedure for numerically investigating resonant
systems, applying it to HL Tau. Section 6.1 provides a
summarized workflow for how one might analyze a potential planetary system embedded in a gas disk, making
reference to more detailed discussion in the main text.
Our conclusions for the HL Tau system can be found in
Sec. 6.2.
2. ASSUMPTIONS
By specializing to protoplanetary disks observed by
ALMA, we can restrict the relevant range of scales, simplifying our analysis. For any pair of putative planets
(throughout the paper, subscripts 1 and 2 to the inner
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 5/2/11
CHANDRA AND XMM-NEWTON OBSERVATIONS OF THE BIMODAL PLANCK SZ-DETECTED CLUSTER
PLCKG345.40-39.34 (A3716) WITH HIGH AND LOW ENTROPY SUBCLUSTER CORES
arXiv:1502.05088v1 [astro-ph.CO] 17 Feb 2015
Felipe Andrade-Santos1 , Christine Jones1 , William R. Forman1 , Stephen S. Murray1,2 , Ralph P. Kraft1 ,
Alexey Vikhlinin1 , Reinout J. van Weeren1 , Paul E. J. Nulsen1 , Laurence P. David1 , William A. Dawson3 ,
Monique Arnaud4 , Etienne Pointecouteau5,6 , Gabriel W. Pratt4 , & Jean-Baptiste Melin7
1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
of Physics and Astronomy, The Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
3 Lawrence Livermore National Laboratory, P.O. Box 808 L- 210, Livermore, CA, 94551, USA
4 Laboratoire AIM, IRFU/Service d’Astrophysique CEA/DSM CNRS Universit´
e Paris Diderot, Bˆ
at. 709, CEA-Saclay, 91191
Gif-sur-Yvette Cedex, France
5 Universit´
e de Toulouse, UPS-OMP, IRAP, F-31028 Toulouse Cedex 4, France
6 CNRS, IRAP, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse Cedex 4, France
7 DSM/Irfu/SPP, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France
Draft version February 19, 2015
2 Department
ABSTRACT
We present results from Chandra, XMM-Newton, and ROSAT observations of the Planck SZdetected cluster A3716 (PLCKG345.40-39.34 – G345). We show that G345 is, in fact, two subclusters
separated on the sky by 400 kpc. We measure the subclusters’ gas temperatures (∼ 2–3 keV), total
(∼ 1–2 ×1014 M⊙ ) and gas (∼ 1–2 ×1013 M⊙ ) masses, gas mass fraction within r500 , entropy profiles,
and X-ray luminosities (∼ 1043 erg s−1 ). Using the gas density and temperature profiles for both
subclusters, we show that there is good (0.8σ) agreement between the expected Sunyaev-Zel’dovich
signal predicted from the X-ray data and that measured from the Planck mission, and better agreement within 0.6σ when we re-computed the Planck value assuming a two component cluster model,
with relative amplitudes fixed based on the X-ray data. Dynamical analysis shows that the two galaxy
subclusters are very likely (> 97% probability) gravitationally bound, and in the most likely scenario,
the subclusters will undergo core passage in 500 ± 200 Myr. The northern subcluster is centrally
peaked and has a low entropy core, while the southern subcluster has a high central entropy. The
high central entropy in the southern subcluster can be explained either by the mergers of several
groups, as suggested by the presence of five giant ellipticals or by AGN energy injection, as suggested
by the presence of a strong radio source in one of its massive elliptical galaxies, or by a combination
of both processes.
Subject headings: galaxy clusters: general — cosmology: large-structure formation
1. INTRODUCTION
Clusters of galaxies are the largest gravitationally
bound structures in the Universe to have reached virial
equilibrium. In the standard ΛCDM cosmology, massive halos dominated by dark matter assemble by accretion of smaller groups and clusters. Under the
influence of gravity, uncollapsed matter and smaller
collapsed halos fall into larger halos and, occasionally, halos of comparable mass merge with one another. Observations of substructures in clusters of galaxies (e.g., Jones & Forman 1984, 1999; Mohr et al. 1995;
Buote & Tsai 1996; Jeltema et al. 2005; Lagan´a et al.
2010; Andrade-Santos et al. 2012, 2013) and the growth
of structure (Vikhlinin et al. 2009) strongly suggest
that clusters formed recently (Forman & Jones 1982;
Richstone et al. 1992). To better understand the evolution of the Universe, it is important to identify and
characterize these large structures.
The first catalog of 189 Sunyaev-Zel’dovich (SZ) clusters detected by the Planck mission was released in
early 2011 (Planck Collaboration et al. 2011). Through
a Chandra XVP (X-ray Visionary Program – PI: Christine Jones) and HRC Guaranteed Time Observations
(PI: Stephen S. Murray), which comprise the ChandraPlanck Legacy Program for Massive Clusters of Galax-
ies1 , we are obtaining Chandra exposures sufficient to
collect at least 10,000 source counts for each of the ESZ
Planck clusters at z ≤ 0.35. PLCKG345.40-39.34 (hereafter G345 - RA = 20:52:16.8, DEC = -52:49:30.7) is a
nearby cluster (z = 0.045). G345 was the first cluster observed as part of the Chandra XVP, which revealed it to
be a double cluster (we distinguish the north and south
subclusters as G345N and G345S, respectively).
Historically, Abell et al. (1989) had catalogued a single cluster of galaxies (A3716 – cross mark in Figure 1)
located ∼ 6′ northwest of G345S and ∼ 6′ southwest of
G345N.
In Section 2 of this paper we present the Chandra,
XMM-Newton, and ROSAT observations and data reduction. In Section 3 we describe the X-ray spatial and
spectral analysis. The total mass and gas mass for each
subcluster are presented in Section 4, and the entropy
profiles of both subclusters are computed in Section 5. In
Section 6 we compute the expected Sunyaev-Zel’dovich
signal from the X-ray data, as well as a comparison between the X-ray emission and the Planck reconstructed
Sunyaev-Zel’dovich map. Finally, we present a dynamical model for the two subclusters in Section 7 and the
conclusions in Section 8.
1 hea-www.cfa.harvard.edu/CHANDRA PLANCK CLUSTERS/
Mon. Not. R. Astron. Soc. 000, 000–000 (0000)
Printed 19 February 2015
(MN LATEX style file v2.2)
Near-IR imaging of T Cha: evidence for scattered-light
disk structures at solar system scales
arXiv:1502.05084v1 [astro-ph.EP] 17 Feb 2015
A.
Cheetham,1? N. Hu´elamo,2 S. Lacour,3 , I. de Gregorio-Monsalvo4,5 and P. Tuthill1
1
2
3
4
5
Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia
Centro de Astrobiolog´ıa (INTA-CSIC); ESAC Campus, P.O. Box 78, E-28691 Villanueva de la Ca˜
nada, Spain
Observatoire de Paris, LESIA/CNRS UMR, 5 place Jules Janssen, Meudon, France
Joint ALMA Observatory (JAO), Alonso de Cordova 3107, Vitacura, Santiago de Chile
European Southern Observatory, Garching bei Munchen, D-85748 Germany
Received 19 January 2015; accepted 15 February 2015
ABSTRACT
T Chamaeleontis is a young star surrounded by a transitional disk, and a plausible
candidate for ongoing planet formation. Recently, a substellar companion candidate
was reported within the disk gap of this star. However, its existence remains controversial, with the counter-hypothesis that light from a high inclination disk may also
be consistent with the observed data. The aim of this work is to investigate the origin
of the observed closure phase signal to determine if it is best explained by a compact
companion. We observed T Cha in the L0 and Ks filters with sparse aperture masking, with 7 datasets covering a period of 3 years. A consistent closure phase signal is
recovered in all L0 and Ks datasets. Data were fit with a companion model and an
inclined circumstellar disk model based on known disk parameters: both were shown
to provide an adequate fit. However, the absence of expected relative motion for an
orbiting body over the 3-year time baseline spanned by the observations rules out the
companion model. Applying image reconstruction techniques to each dataset reveals
a stationary structure consistent with forward scattering from the near edge of an
inclined disk.
Key words: Techniques: high angular resolution, interferometric – Stars: circumstellar material, planetary systems – Stars: individual: T Cha
1
INTRODUCTION
Transitional disks represent an important phase in the evolution of circumstellar disks; a stepping stone between a gas
rich protoplanetary disk and a rocky, gas poor debris disk.
The characteristic deficit of mid-IR emission that defines this class is explained by a lack of dust at intermediate
distances from the central star, leading to the appearance
of an annular gap in their disks. These holes or gaps can
be created by dust clearing due to the presence of a planetary system, causing transition disks to be prime targets
for young planetary systems in formation. However, alternative explanations exist, such as the presence of a close
binary companion (e.g. CoKu Tau 4, Ireland & Kraus 2008)
or grain growth (Dullemond & Dominik 2005).
In the last few years, several substellar companion candidates have been reported within the annular gaps of transition disks. Among these are T Cha, HD 142527 and LkCa
15 (Hu´elamo et al. 2011; Biller et al. 2012; Kraus & Ireland
?
Email: [email protected]
c 0000 RAS
2012), all based on the Sparse Aperture Masking (SAM)
technique (Tuthill et al. 2006). In this paper, we investigate
the nature of the companion candidate reported for T Cha.
T Cha is a G8 star at an estimated distance of 108±9 pc
(Torres et al. 2008). The star is surrounded by an inner
disk (Olofsson et al. 2013) and an outer disk separated by
a dust gap of ∼20 AU (Hu´elamo et al. 2015). The highly
inclined nature of its outer disk has led to the suggestion
that forward scattering may produce closure phases similar
to those expected from the reported companion candidate
(Olofsson et al. 2013). To investigate this claim, we have
obtained 4 epochs of SAM data for T Cha covering 3 years,
to look for evidence of orbital motion that would reveal the
nature of this object.
Several multiwavelength studies of the T Cha system
have been performed, yielding precise constraints on the
disk and stellar parameters (Cieza et al. 2011; Olofsson et al.
2011, 2013; Hu´elamo et al. 2015). Based on these constraints,
we adopt a mass of 1.5 M , temperature of 5400 K and stellar radius of 1.4 R . The outer disk around T Cha has been
recently imaged by ALMA (Hu´elamo et al. 2015), showing
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 12/16/11
arXiv:1502.05072v1 [astro-ph.CO] 17 Feb 2015
NEW LIMITS ON POLARIZED POWER SPECTRA AT 126 AND 164 MHZ: RELEVANCE TO EPOCH OF
REIONIZATION MEASUREMENTS
David F. Moore1 , James E. Aguirre1 , Aaron R. Parsons2,8 , Zaki S. Ali2 , Richard F. Bradley3,4,5 , Chris L.
Carilli6,7 , David R. DeBoer8 , Matthew R. Dexter8 , Nicole E. Gugliucci5 , Daniel C. Jacobs9 , Pat Klima4 ,
Adrian Liu2,10 , David H. E. MacMahon8 , Jason R. Manley11 , Jonathan C. Pober12 , Irina I. Stefan7 , and William
P. Walbrugh11
Draft version February 19, 2015
ABSTRACT
Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation
power spectrum of highly redshifted 21 cm Hi emission from the epoch of reionization, yet observational
constraints on the level of polarized emission are poor. Using the Donald C. Backer Precision Array
for Probing the Epoch of Reionization (PAPER), we present the first limits on the power spectra of
all four Stokes parameters in two frequency bands, centered at 126 MHz (z = 10.3) and 164 MHz
(z = 7.66). This data comes from from a three-month observing campaign of a 32-antenna deployment,
for which unpolarized power spectrum results have been reported at z = 7.7 (Parsons et al. 2014) and
7.5 < z < 10.5 (Jacobs et al. 2014). The power spectra in this paper are processed in the same
way, and show no definitive detection of polarized power. The limits are sufficiently low that we
are able to show that the excess unpolarized power reported in those works is not due to leakage of
Faraday-rotated polarized emission. Building upon the Moore et al. (2013) simulations of polarized
point sources, we further argue that our upper limits and previous observations imply that the mean
polarization fraction of point sources at these frequencies is ∼ 2 × 10−3, roughly an order of magnitude
lower than that observed for point sources at 1.4 GHz.
1. INTRODUCTION
The study of the polarized sky at meter wavelengths
has become critical to the characterization of foregrounds
for the 21 cm Epoch of Reionization (EoR) power spectrum. Smooth-spectrum sources occupy a well-defined
“wedge” in the cylindrical (kk , k⊥ ) space of the EoR
power spectrum (e.g., Morales et al. 2012; Pober et al.
2013), but it has been understood for some time (e.g.,
Geil et al. 2011; Pen et al. 2009) that Faraday-rotated
polarized emission can contaminate the measurement of
unpolarized EoR emission. When mapped into the power
spectrum, this Faraday-rotated polarized emission generates power which mimics the high kk emission of the
EoR, scattering power into the otherwise-clean EoR window. Moore et al. (2013, hereafter M13) simulated the
effects of the low-level forest of weak, polarized point
sources leaking into the full 3-D EoR power spectrum.
M13 used the best measurements available at the time,
notably the existing (unpolarized) point source surveys
(Hales et al. 1988; Cohen et al. 2004) and the few polar1 Dept. of Physics and Astronomy, U. Pennsylvania, Philadelphia, PA
2 Astronomy Dept., U. California, Berkeley, CA
3 Dept. of Electrical and Computer Engineering, U. Virginia,
Charlottesville, VA
4 National Radio Astronomy Obs., Charlottesville, VA
5 Dept. of Astronomy, U. Virginia, Charlottesville, VA
6 National Radio Astronomy Obs., Socorro, NM
7 Cavendish Lab., Cambridge, UK
8 Radio Astronomy Lab., U. California, Berkeley, CA
9 School of Earth and Space Exploration, Arizona State U.,
Tempe, AZ
10 Berkeley Center for Cosmological Physics, U. California,
Berkeley, CA
11 Square Kilometer Array, South Africa Project, Cape Town,
South Africa
12 Dept. of Physics, U. Washington, Seattle, WA
ization measurements below 200 MHz (Pen et al. 2009;
Bernardi et al. 2009), and found that the level of emission could be problematic for any experiment with ∼ 1%
leakage from Stokes Q → I. However, this result is
highly dependent on the actual level of polarized emission present at these frequencies.
Recent new measurements of the polarized sky at low
frequencies have begun to reveal some features of both
diffuse and point source emission. Bernardi et al. (2013,
hereafter B13) conducted a 2400 square degree survey
using the Murchison Widefield Array (MWA) at 189
MHz. They detected one polarized point source, PMN
J0351-2744, with a polarization fraction of 2% and polarized flux of 320 mJy. At somewhat higher frequencies (350 MHz), Gieß¨
ubel et al. (2013) studied sources
behind M31 and found that the fractional polarization
decreased by a factor of 0.14 relative to that measured
at 1.4 GHz, with typical values of a few percent at 350
MHz. A study of the depolarization of point sources by
Farnes et al. (2014) showed a systematic trend for depolarization of steep-spectrum point sources as frequency
decreased, resulting in very low polarization fractions below 300 MHz. In diffuse emission, both Bernardi et al.
(2013) and Jeli´c et al. (2014, hereafter Je14) report significant amounts of weakly polarized emission with rotation measure (RM) |Φ| ≤ 25 rad m−2 , which they attributed to Galactic synchrotron. It is worth noting that
the reported rotation measures for point sources vary
over the range of values reported in the all-sky RM map
of Oppermann et al. (2012), whereas the diffuse emission
typically shows low RM. This may be due to a subtle
difference between the Faraday depths of diffuse emission and point sources. Since we expect diffuse emission
mostly to be generated within our Galaxy, we expect
it to be emitting from within the magnetized, ionized
Submitted to ApJ
Preprint typeset using LATEX style emulateapj v. 5/2/11
PROPER MOTIONS AND TRAJECTORIES FOR 16 EXTREME RUNAWAY AND HYPERVELOCITY STARS1
Warren R. Brown2 , Jay Anderson3 , Oleg Y. Gnedin4 , Howard E. Bond3,5 , Margaret J. Geller2 , and Scott J.
Kenyon2
arXiv:1502.05069v1 [astro-ph.SR] 17 Feb 2015
2 Smithsonian Astrophysical Observatory, 60 Garden St, Cambridge, MA 02138
3 Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218
4 Department of Astronomy, University of Michigan, Ann Arbor, MI 48109
5 Department
of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802
Submitted to ApJ
ABSTRACT
We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars,
mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically
consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center
“hypervelocity star” origin. The trajectory of HE 0437–5439 is consistent with both Milky Way and
Large Magellanic Cloud origins. A Galactic center origin is excluded at 3σ confidence for two of
the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways
and provide evidence for ∼500 km s−1 ejections from the disk. We also measure a significant proper
motion for the unbound sdO star US 708. Its 1,000 km s−1 motion is in some tension with proposed
supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the
future, we expect Gaia will better constrain the origin of these remarkable unbound stars.
Subject headings: stars: early-type — kinematics and dynamics — Galaxy: halo — Galaxy: kinematics and dynamics — Magellanic Clouds
1. INTRODUCTION
The Milky Way is a gravitationally bound system
of a few 1011 stars, but it also hosts some unbound
stars. Where these unbound stars come from is an open
question. It is very difficult to explain unbound mainsequence stars by supernova explosions in close binary
systems (Blaauw 1961) or by dynamical encounters between binaries (Poveda et al. 1967), the mechanisms that
explain disk “runaway” stars. Simulations show that
stars cannot be launched at speeds greater than their binary orbital velocity, and this velocity cannot exceed the
escape velocity from a star’s surface, since that would require the stars to be orbiting inside one another (Leonard
1991). The escape velocity from the surface of mainsequence stars is comparable to the 500 – 600 km s−1
local Galactic escape velocity (Piffl et al. 2014). Only
dynamical interactions with an object more massive and
compact than a star can easily explain unbound mainsequence stars. Hills (1988) predicted that three-body
exchange interactions between stars and the central massive black hole (MBH) will inevitably unbind stars from
the Galaxy, and he called such objects “hypervelocity
stars” (HVSs). We therefore refer to unbound mainsequence stars in this paper as HVSs.
The first HVS was discovered by Brown et al. (2005),
and a couple dozen B-type HVSs are now known in
the Milky Way halo (Edelmann et al. 2005; Brown et al.
2006a,b, 2007a,b; Heber et al. 2008; Brown et al. 2009,
2012b, 2014; Zheng et al. 2014). These stars have radial
velocities up to +700 km s−1 ; no star moving towards us
[email protected];[email protected];[email protected]
1 Based on observations with the NASA/ESA Hubble Space
Telescope obtained at the Space Telescope Science Institute,
and from the data archive at STScI, which is operated by the
Association of Universities for Research in Astronomy, Inc.,
under NASA contract NAS5-26555.
with a comparable radial velocity has ever been observed.
To date, detailed spectroscopic analyses of the B-type
HVSs find that they are main-sequence B stars at 10 –
100 kpc distances (Edelmann et al. 2005; Przybilla et al.
2008a,b; L´
opez-Morales & Bonanos 2008; Brown et al.
2012a, 2013). The short lifetimes of B stars require that
these unbound stars were ejected from a region with recent star formation, such as the Galactic disk or the
Galactic center.
Full three-dimensional trajectories for the HVSs are a
crucial test of their origin. Measuring the trajectories for
the HVSs, however, requires accurate distances and absolute proper motions. Distance estimates to individual
stars are determined by comparing either spectroscopic
Teff and log g or broadband colors to stellar evolution
tracks. Typical precision is 15% (Brown et al. 2014);
the accuracy depends on the choice of stellar evolution
tracks (i.e. assumptions about metallicity). Proper motions are equally difficult to measure, because HVSs are
distant and they should be on radial trajectories. Expected HVS proper motions are typically . 1 mas yr−1
and thus cannot be measured with ground based surveys. The only unbound star with a statistically significant proper motion measurement to date is the runaway
B star HD 271791, a star that was ejected in the direction of Galactic rotation from the outer disk (Heber et al.
2008). HIP 60350 is another example of a high velocity B
star with a significant proper motion measurement, but
it is unbound at less than 1σ significance (Irrgang et al.
2010). Here, we present Hubble Space Telescope (HST)
proper motion measurements for a sample of 16 stars
with extreme radial velocities, 12 of which are unbound
to the Milky Way in radial velocity alone.
Proper motions are important because they can discriminate between Galactic center and Galactic disk origins. The star US 708 (hereafter HVS 2), for example, is
Mon. Not. R. Astron. Soc. 000, 1–7 (2014)
Printed 19 February 2015
(MN LATEX style file v2.2)
arXiv:1502.05068v1 [astro-ph.HE] 17 Feb 2015
X-ray Emission from Charge Exchange in the Cygnus Loop
SNR
Shawn
R. Roberts1⋆ , Q. Daniel Wang1
1
Department of Astronomy, University of Massachusetts, Amherst, MA 01002, USA
Submitted to MNRAS, September, 2014
ABSTRACT
The Cygnus Loop has been the focus of substantial debate concerning the contribution of charge exchange (CX) to supernova remnant (SNR) X-ray emission. We
take advantage of a distinct feature of CX, enhanced Kα forbidden line emission, and
employ the energy centroid of the OVII Kα triplet as a diagnostic. Based on X-ray
spectra extracted from an extensive set of Suzaku observations, we measure the energy
centroid shifts of the triplet on and off the shock rim of the remnant. We find that
enhanced forbidden to resonance line emission exists throughout much of the rim and
this enhancement azimuthally correlates with non-radiative Hα filaments, a tracer of
strong neutral-plasma interaction in the optical. We also show that alternative mechanisms cannot explain the enhancement observed. These results demonstrate the need
to model the CX contribution to the X-ray emission of SNRs, particularly for shocks
propagating in a partially neutral medium. Such modelling may be critically important to the correct measurements of the ionization, thermal, and chemical properties
of SNRs.
Key words: atomic processes; ISM: supernova remnants; ISM: abundances; ISM:
individual: Cygnus Loop; X-rays: ISM
1
INTRODUCTION
X-ray astrophysics is layered with a host of interesting phenomena. As spectral capabilities continue to grow within
the X-ray regime, our ability to discern the importance of
these varied phenomena has led to many interesting surprises. One X-ray emission mechanism whose importance
we are increasingly becoming aware of is charge exchange
or charge transfer (CX; e.g., Lallement 2004). While CX reactions have been studied in an astrophysical context for
decades, their importance for X-ray astrophysics wasn’t realized until comet observations in 1996 (Lisse et al. 1996;
Cravens 1997). Since this cometary epiphany, CX has been
proposed and studied in a wide range of X-ray sources, from
planetary atmospheres (Dennerl 2002; Dennerl et al. 2006;
Kharchenko et al. 2008) through individual star-forming
complexes (Townsley et al. 2011) and star-forming galaxies
(Tsuru et al. 2007; Liu et al. 2011, 2012) to the cores of clusters (Fabian et al. 2011). More recently, the first evidence for
X-ray emitting CX in a supernova remnant (SNR) has been
suggested (Katsuda et al. 2011).
CX is a special case of inelastic collision in which one or
more electrons is exchanged between an ion and typically a
neutral atom, most likely hydrogen. As the colliding particles approach each other, their potentials overlap, allowing
the photonless exchange and subsequent cascade of the captured electron. This is an extremely efficient process with
a relatively large reaction cross section (about 5 orders of
⋆
E-mail: [email protected]
magnitude greater than for electron collisional excitation).
If the capturing ion is sufficiently ionized (e.g. OVIII) then
the cascading is expected to result in X-ray emission. Thus,
in principle, X-ray emitting CX can be an important process at any astrophysical site where heavily ionized species
interact with neutrals, making the thin post-shock region of
SNRs a very promising site for CX.
While CX has been used to explain the broad component of Hα emission in many SNRs for more than 30 years
(Chevalier et al. 1980; Ghavamian et al. 2001), X-ray emission produced by CX has not yet been directly apparent.
This is at least partly due to the lack of a non-dispersive
high-resolution spectrometer. We are instead confined to Xray CCD data for our studies of diffuse emission. At the
spectral resolution of such data, the pure line emission of
CX can mimic the continuum + line emission of a thermal
plasma, making it difficult to disentangle the relative contributions of the different emission mechanisms. Thus, without
accounting for the contribution from CX, the modelling of
the properties of the shocked plasma (e.g. chemical abundances) can potentially be misleading.
To delve into the importance of CX for SNR emission
we need to be able to spatially and spectrally resolve the
nuances of the plasma evolution in nearby remnants, such
as the Cygnus Loop SNR. This remnant has an angular diameter of ≈ 3◦ and suffers relatively little absorption. Such
a close SNR provides a unique laboratory for the study of
shock physics and non-equilibrium conditions. These spatial
and spectral studies have been done to some extent in the
X-ray regime with Suzaku and XMM-Newton, for which ob-
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 2/19/04
THE COMPLEX NORTH TRANSITION REGION OF CENTAURUS A: RADIO STRUCTURE
Susan G. Neff
NASA’s Goddard Space Flight Center, Laboratory for Observational Cosmology, Mail Code 665, Greenbelt, Maryland, 20771
Jean A. Eilek
arXiv:1502.05066v1 [astro-ph.GA] 17 Feb 2015
Physics Department, New Mexico Tech, Socorro NM 87801 and
National Radio Astronomy Observatory∗ † , Socorro NM 87801
Frazer N. Owen
National Radio Astronomy Observatory∗ , P. O. Box O, Socorro NM 87801
Draft version February 19, 2015
ABSTRACT
We present deep radio images of the inner ∼ 50 kpc of Centaurus A, taken with the Karl G. Jansky
Very Large Array (VLA) at 90 cm. We focus on the Transition Regions between the inner galaxy
– including the active nucleus, inner radio lobes and star-forming disk – and the outer radio lobes.
We detect previously unknown extended emission around the Inner Lobes, including radio emission
from the star-forming disk. We find that the radio-loud part of the North Transition Region, known
as the North Middle Lobe, is significantly overpressured relative to the surrounding ISM. We see
no evidence for a collimated flow from the Active Galactic Nucleus (AGN) through this region. Our
images show that the structure identified by Morganti et al. (1999) as a possible large-scale jet appears
to be part of a narrow ridge of emission within the broader, diffuse, radio-loud region. This knotty
radio ridge is coincident with other striking phenomena: compact X-ray knots, ionized gas filaments,
and streams of young stars. Several short-lived phenomena in the North Transition Region, as well as
the frequent re-energization required by the Outer Lobes, suggest that energy must be flowing through
both Transition Regions at the present epoch. We suggest that the energy flow is in the form of a
galactic wind.
Subject headings: galaxies: active – galaxies: individual (NGC 5128, Centaurus A) – galaxies: jets –
radio continuum: galaxies
1. INTRODUCTION
The radio source Centaurus A (“Cen A”) and its parent
galaxy, NGC 5128, form an important active-galaxy system which can reveal the interplay between the history
of the parent galaxy and the development of the radio
source it creates. Only 3.8 Mpc away, (Harris et al. 2010;
10 ' 1.14 kpc ), Cen A/NGC 5128 can be studied with a
level of detail unavailable in other systems. Interestingly
– or perhaps because the system is so close – both the
stellar galaxy and the radio source show details that we
do not normally expect.
NGC 5128 is fundamentally a normal elliptical galaxy,
dominated by an old stellar population with kinematic signatures typical of other massive ellipticals. The galaxy’s
optical appearance is, however, somewhat unusual, being dominated by the iconic central dust band, which is
thought to be the result of a merger with a smaller gas-rich
system. Other remnants of the merger may include cold
HI still orbiting the galaxy and a ribbon of young stars
and emission-line gas clouds (e.g., Graham 1998, Mould et
al. 2000, Struve et al. 2010). In a companion paper (Neff,
Electronic address: [email protected]
National Radio Astronomy Observatory is a facility of the
National Science Foundation operated under cooperative agreement
by Associated Universities Inc.
† Adjunct Astronomer at the National Radio Astronomy Observatory.
∗ The
Eilek, & Owen 2014, Paper 2), we show that ribbon extends to ∼ 35 kpc from the galaxy, is coincident with a
complex, knotty structure seen in radio and X-rays, and
is probably energized by a galactic wind moving through
the middle regions of the radio source.
Our focus in this paper is the radio structure of the
inner ∼ 50 kpc of Cen A. We present new VLA 90 cm
observations, and discuss their implications for the astrophysics of the region and the overall energy flow in the
system. In Section 2, we review relevant background information on the Cen A / NGC 5128 system and present
a first look at our results. Section 3 describes our radio
observations and imaging methods. Section 4 presents
our 90 cm image, and discusses its limitations. Section 5
compares our image to the 20 cm results from Morganti
et al. , 1999 (“M99”). In Sections 6 and 7, we present and
analyze various components of the Inner and Transition
regions. In Section 8, we discuss the need for energy flow
through the Transition Regions and argue that those regions are dominated by a wind rather than jets. Section
9 summarizes our results.
2. CEN A: SETTING THE STAGE
We begin with the three spatial scales of the system –
small, large and intermediate, as summarized in Figure 1
and Table 1 – and continue with an introduction to key
IRFU-15-01
Prospects for Annihilating Dark Matter in the inner Galactic halo
by the Cherenkov Telescope Array
Valentin Lefranc∗ and Emmanuel Moulin†
CEA Irfu, Service de Physique des Particules, Centre de Saclay, Gif-sur-Yvette 91190, France
Paolo Panci‡
Institut d’Astrophysique de Paris, UMR 7095 CNRS,
Universit´e Pierre et Marie Curie, 98 bis Boulevard Arago, Paris 75014, France
arXiv:1502.05064v1 [astro-ph.HE] 17 Feb 2015
Joseph Silk§
Institut d’Astrophysique de Paris, UMR 7095 CNRS,
Universit´e Pierre et Marie Curie, 98 bis Boulevard Arago, Paris 75014, France
Department of Physics and Astronomy, 3701 San Martin Drive,
The Johns Hopkins University, Baltimore MD 21218, USA and
BIPAC, 1 Keble Road, University of Oxford, Oxford OX1 3RH UK
We compute the sensitivity to dark matter annihilations for the forthcoming large Cherenkov
Telescope Array (CTA) in several primary channels and over a range of dark matter masses from
30 GeV up to 80 TeV. For all channels, we include inverse Compton scattering of e± by dark
matter annihilations on the ambient photon background, which yields substantial contributions to
the overall γ-ray flux. We improve the analysis over previous work by: i) implementing a spectral
and morphological analysis of the γ-ray emission; ii) taking into account the most up-to-date cosmic
ray background obtained from a full CTA Monte Carlo simulation and a description of the diffuse
astrophysical emission; and iii) including the systematic uncertainties in the rich observational CTA
datasets. We find that our spectral and morphological analysis improves the CTA sensitivity by
roughly a factor 2. For the hadronic channels, CTA will be able to probe thermal dark matter
candidates over a broad range of masses if the systematic uncertainties in the datasets will be
controlled better than the percent level. For the leptonic modes, the CTA sensitivity will be well
below the thermal value of the annihilation cross-section. In this case, even with larger systematics,
thermal dark matter candidates up to masses of a few TeV will be easily studied.
PACS numbers: 95.35.+d, 95.30.Cq, 98.35.-a, 95.85.Pw
I.
INTRODUCTION
A substantial body of astrophysical data (such as galaxy rotation curves and weak lensing observations, Large Scale
Structure surveys and the precise data from CMB measurements [1, 2]) have revealed the crucial gravitational role of a
new kind of particle, dubbed as Dark Matter (DM). This provides one of the most compelling hints for beyond Standard
Model (SM) physics. Nevertheless, in order to identify the microphysics nature of this particle, we need information
on its mass and interaction characteristics. One common strategy dubbed as indirect detection of DM particles aims
at detecting the secondary stable SM products of a DM annihilation or decay in our Galaxy and beyond, on top of the
astrophysical background. For popular DM candidates, belonging to the class of Weakly Interacting Massive Particles
(WIMPs), these SM products are injected in the realm of investigation of very high energy astrophysics. In particular,
the current and upcoming ground-based telescope, are very appropriate for probing the properties of WIMPs in the
mass domain of few 100 GeV up to 100 TeV, by looking at the very high energy γ-rays from the annihilation of such
particles.
ˇ
Over the last several years, in order to identify the nature of DM, the Imaging Atmospheric (or Air) Cerenkov
Telescopes (IACTs), such as H.E.S.S. [3], MAGIC [4] and VERITAS [5], have conducted rich observational programs
in the inner region of the Milky Way halo [6–8], nearby satellite galaxies [9–18] and galaxy clusters [19–21]. So far,
∗ Electronic
address:
address:
‡ Electronic address:
§ Electronic address:
† Electronic
[email protected]
[email protected]
[email protected]
[email protected]
Mon. Not. R. Astron. Soc. 000, 1–7 (2014)
Printed 19 February 2015
(MN LATEX style file v2.2)
arXiv:1502.05054v1 [astro-ph.HE] 17 Feb 2015
Chandra Identification of Two AGN Discovered by INTEGRAL
John A. Tomsick1⋆ , Roman Krivonos1,2, Farid Rahoui3,4, Marco Ajello5 , Jerome Rodriguez6,
Nicolas
Barriere1, Arash Bodaghee7, and Sylvain Chaty6,8
1
Space Sciences Laboratory, 7 Gauss Way, University of California, Berkeley, CA 94720-7450, USA
Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia
3 European Southern Observatory, Karl Schwarzschild-Strasse 2, 85748 Garching bei Munchen, Germany
4 Department of Astronomy, Harvard University, 60 Garden Street, Cambridge, MA 02138, USA
5 Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA
6 Laboratorie AIM (UMR-E 9005 CEA/DSM-CNRS-Université Paris Diderot), Irfu/Service d’Astrophysique, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France
7 Georgia College & State University, CBX 82, Milledgeville, GA 31061, USA
8 Institut Universitaire de France, 103 Boulevard Saint-Michel, 75005 Paris, France
2 Space
19 February 2015
ABSTRACT
Here, we report on observations of two hard X-ray sources that were originally discovered with the INTEGRAL satellite: IGR J04059+5416 and IGR J08297–4250. We use the
Chandra X-ray Observatory to localize the sources and then archival near-IR images to identify the counterparts. Both sources have counterparts in the catalog of extended 2 Micron
All-Sky Survey sources, and the counterpart to IGR J04059+5416 has been previously identified as a galaxy. Thus, we place IGR J04059+5416 in the class of Active Galactic Nuclei
(AGN), and we suggest that IGR J08297–4250 is also an AGN. If this identification is correct, the near-IR images suggest that the host galaxy of IGR J08297–4250 may be merging with a smaller nearby galaxy. For IGR J04059+5416, the 0.3–86 keV spectrum from
Chandra and INTEGRAL is consistent with an absorbed power-law with a column density
22
−2
of NH = (3.1+2.0
and a photon index of Γ = 1.4 ± 0.7, and we suggest that it
−1.5 ) × 10 cm
is a Seyfert galaxy. For IGR J08297–4250, the photon index is similar, Γ = 1.5 ± 0.8, but the
23
−2
source is highly absorbed (NH = (6.1+10.1
).
−4.3 ) × 10 cm
Key words: X-rays: general, X-rays: galaxies, galaxies: active, galaxies: Seyfert, stars: individual(IGR J04059+5416, IGR J08297–4250)
1 INTRODUCTION
The hard X-ray imaging by the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) satellite (Winkler et al. 2003) has
led to the discovery of a large number of new or previously poorly
studied “IGR” sources. The most recent published catalogs of 20–
100 keV sources detected by the Imager on-Board the INTEGRAL
Satellite (IBIS; Ubertini et al. 2003) include more than 700 sources
for the whole sky (Bird et al. 2010) and 402 sources within 17.5◦ of
the Galactic plane1 (Krivonos et al. 2012). The current INTEGRAL
source catalog2 includes 954 sources detected by IBIS. The nature
of ≈200 of the IGR sources is still unknown according to the most
up-to-date list3 .
While INTEGRAL has found many new sources, it only localizes them to 1′ –5′ , which is not nearly adequate for finding
⋆
1
2
3
E-mail: [email protected] (JAT)
see http://hea.iki.rssi.ru/integral
see http://isdc.unige.ch/integral/catalog/38/catalog.html
see http://irfu.cea.fr/Sap/IGR-Sources/
c 2014 RAS
optical/IR counterparts. Short exposures of IGR sources with the
Chandra X-ray Observatory allow for a major advance in understanding the nature of these sources by providing sub-arcsecond
positions, leading to unique optical/IR counterparts, as well as
0.3–10 keV spectra that can be used to measure column densities and continuum shapes (see Tomsick et al. 2006, 2008, 2009,
2012; Bodaghee et al. 2012, for examples of our previous work on
Chandra follow-up of IGR sources). Other groups have been carrying out similar investigations (e.g., Fiocchi et al. 2010; Ratti et al.
2010; Paizis et al. 2011, 2012; Nowak et al. 2012; Karasev et al.
2012).
Prior to this study, little was known about IGR J04059+5416
and IGR J08297–4250. They were discovered in the Krivonos et al.
(2012) hard X-ray survey, which includes 9-years of INTEGRAL
data. They are listed as being unidentified in the Krivonos et al.
(2012) catalog, and there have not been any more publications on
these sources. A search of the SIMBAD database does not suggest
any likely counterparts. The Chandra observations that we made
to identify the sources are described in Section 2. Then, the results
Mon. Not. R. Astron. Soc. 000, 000–000 (0000)
Printed February 19, 2015
(MN LATEX style file v2.2)
arXiv:1502.05053v1 [astro-ph.GA] 17 Feb 2015
The rise and fall of stellar disks across the peak of cosmic
star formation history: mergers versus smooth accretion
C. Welker1,2⋆ , Y. Dubois1, J. Devriendt2,3, C. Pichon1,4, S. Kaviraj5 and S. Peirani1
1
2
3
4
5
CNRS and UPMC Univ. Paris 06, UMR 7095, Institut d’Astrophysique de Paris, 98 bis Boulevard Arago, F-75014 Paris, France
Sub-department of Astrophysics, University of Oxford, Keble Road, Oxford OX1 3RH
Observatoire de Lyon, UMR 5574, 9 avenue Charles Andr´
e, Saint Genis Laval 69561, France
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA
Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts, AL10 9AB
Accepted . Received ; in original form
ABSTRACT
Building galaxy merger trees from a state-of-the-art cosmological hydrodynamics simulation, Horizon-AGN, we perform a statistical study of how mergers and smooth accretion drive galaxy morphologic properties above z > 1. More specifically, we investigate
how stellar densities, effective radii and shape parameters derived from the inertia tensor depend on mergers of different mass ratios. We find strong evidence that smooth
accretion tends to flatten small galaxies over cosmic time, leading to the formation of
disks. On the other hand, mergers, and not only the major ones, exhibit a propensity
to puff up and destroy stellar disks, confirming the origin of elliptical galaxies. We also
find that elliptical galaxies are more susceptible to grow in size through mergers than
disc galaxies with a size-mass evolution r0.5 ∝ Ms1.2 instead of r0.5 ∝ Ms−0.5 − M 0.5
depending on the merger mass ratio. The gas content drive the size-mass evolution due
to merger with a faster size growth for gas-poor galaxies r0.5 ∝ Ms2 than for gas-rich
galaxies r0.5 ∝ Ms .
Key words: galaxies: formation — galaxies: evolution — galaxies: interactions —
galaxies: kinematics and dynamics — methods: numerical
1
INTRODUCTION
While our present understanding of galaxy evolution derives mainly from the nearby (z < 1) Universe, the bulk
of today’s stellar mass formed around the broad peak of
cosmic star formation history at z ∼ 2 (e.g. Madau et al.
1998; Hopkins & Beacom 2006). Although it represents a
significant epoch in the evolution of the observable Universe, the properties of galaxies remain largely unexplored at
this epoch, as it has only recently become accessible by current observational facilities (CANDELS, GOODS, Herschel,
ALMA, Chen et al. 2014; Lin et al. 2010; Cooper et al.
2012; L´
opez-Sanjuan et al. 2013).
As a result, what constitutes arguably the most important aspect of hierarchical galaxy formation and evolution
is still being debated. To what extent mergers, as opposed
to secular evolution driven by (cold) gas inflows, explain
the diversity of galaxies? The significance of mergers, considered a cornerstone of the bottom-up growth of galaxies,
has been heavily debated in recent work. They are certainly
⋆
E-mail: [email protected]
capable of inducing star formation, black hole growth and
morphological transformations (e.g. Springel et al. 2005),
but it is not obvious that, at z > 1, mergers drive the
evolution of galaxy properties like stellar mass, size and
morphology (Shankar et al. 2004; Law 2009; Kaviraj et al.
2013) considering the steady input of accreted streams of
gas (Kereˇs et al. 2005; Ocvirk et al. 2008; Dekel et al. 2009)
and the gas-rich nature of galaxies (Tacconi et al. 2010;
Santini et al. 2014).
Observational studies suggest that large fractions of
star-forming galaxies around z ∼ 2 are not mergers
but show kinematics and visual morphologies that are
more consistent with systems dominated by turbulent
discs (e.g. Forster Schreiber et al. 2006; Shapiro et al. 2008;
Genzel & Burkert 2008; Mancini et al. 2011; Kaviraj et al.
2013). In addition, many primordial spheroids that are forming the bulk of their stellar mass at z ∼ 2 do not show
the tidal features that would be expected from recent major mergers (Kaviraj et al. 2013). Notwithstanding these advances, large statistical studies remain difficult, both due to
the fact that observational efforts rely on pencil-beam surveys that are susceptible to low-number statistics and cosmic
S UBMITTED TO A P JL
Preprint typeset using LATEX style emulateapj v. 5/2/11
EVERY INTERACTING DOUBLE WHITE DWARF BINARY MAY MERGE
K EN J. S HEN 1
arXiv:1502.05052v1 [astro-ph.SR] 17 Feb 2015
Submitted to ApJL
ABSTRACT
Interacting double white dwarf binaries can give rise to a wide variety of astrophysical outcomes ranging
from faint thermonuclear and Type Ia supernovae to the formation of neutron stars and stably accreting AM
Canum Venaticorum systems. One key factor affecting the final outcome is whether mass transfer remains
dynamically stable or instead diverges, leading to the tidal disruption of the donor and the merger of the binary.
It is typically thought that for low ratios of the donor mass to the accretor mass, mass transfer remains stable,
especially if accretion occurs via a disk. In this Letter, we examine disk-accreting binaries with extremely low
mass ratios and find that the initial phase of hydrogen-rich mass transfer leads to a classical nova-like outburst
on the accretor. Dynamical friction within the expanding nova shell shrinks the orbit and causes the mass
transfer rate to increase dramatically above the accretor’s Eddington limit, resulting in a binary merger. While
further calculations are necessary to confirm this outcome for the entire range of binaries previously thought to
be dynamically stable, it appears likely that most, if not all, interacting double white dwarf binaries will merge
during the course of their evolution.
Subject headings: binaries: close— novae, cataclysmic variables— nuclear reactions, nucleosynthesis,
abundances— supernovae: general— white dwarfs
1. INTRODUCTION
Double white dwarf (WD) binaries are formed via common
envelope evolution in binary systems. The emission of gravitational waves in the closest double WD binaries shrinks their
orbits and causes mass transfer between the WDs within a
Hubble timescale. The resulting interaction leads to a broad
range of outcomes including explosive transient phenomena
and the formation of exotic classes of stars.
Previous work has found that the dynamical stability of
mass transfer in these binaries, and thus their ultimate fates,
depends on the mass ratio of the donor mass to the accretor mass, q = M2 /M1 , and on the timescale for angular momentum feedback from the accretor to the orbit (Marsh et al.
2004). However, systems are generally assumed to be dynamically stable if accretion occurs via a disk, as torques on
the extended lever arm of the disk should yield efficient angular momentum feedback. In this Letter, we challenge this
assumption by showing that even disk-accreting double WD
binaries with extremely low-mass He WD donors will merge
unstably.
The overlying non-degenerate H-rich layers surrounding
these low-mass WDs alters the initial phase of mass transfer compared to simplified pure He WD donors (D’Antona
et al. 2006; Kaplan et al. 2012) and leads to classical novalike outbursts (Shen et al. 2013). Here, we consider the impact of the expanding nova shell on the binary orbit. We calculate the evolution of two fiducial double WD binaries with
masses 0.1 + 0.5 and 0.1 + 1.0 M and find that these binaries,
which were previously thought to form stable mass transfer
AM Canum Venaticorum systems, instead merge due to dynamical friction within the accretor’s expanded outer layers.
In Section 2, we estimate the shrinking of the orbital separation due to mass ejection via dynamical friction within the
common envelope formed by the nova shell. In Section 3,
we provide a prescription for how the donor’s mass transfer
1 Department of Astronomy and Theoretical Astrophysics Center, University of California, Berkeley, CA 94720, USA.
rate responds to the decrease in orbital separation. In Section
4, we describe our fiducial calculations, and we consider the
implications of our findings in Section 5.
2. COMMON ENVELOPE ENERGETICS
We begin by estimating the change in the binary separation,
a, if some or all of the energy to eject the expanding nova
shell comes from the orbit via dynamical friction of the two
WDs with the common envelope. This is at best a rough approximation (see Ivanova et al. 2013 and references therein);
however, we use it here for simplicity and defer analysis with
alternative common envelope physics to future work.
The specific energy required to eject material to infinity
from the Roche surface of the accreting WD is
"
2 #
GMt 1/(1 + q) q/(1 + q) 1
1
ΦL1 = −
+
−
− xL1
a
1 − xL1
xL1
2 1+q
GMt
f (q),
(1)
a
where the total mass of the system is Mt = M1 + M2 , and xL1
is the distance from the center of the donor to the L1 Lagrange point in units of the binary separation. The function
f (q) varies from 1.4 − 1.9 for q = 0.1 − 0.5. While this specific energy changes slightly as mass is being ejected from
the system, it will remain correct to first order in the ratio of
the ejecta mass, Mej , to the total mass and thus is sufficient for
our purposes of estimation.
Assuming that a fraction αCE of the energy to eject Mej
comes from the binding energy of the orbit, we obtain the
change in binary separation following the common envelope
event,
Mej
da
2(1 + q) f (q)
=−
αCE
+1 .
(2)
a
M1
q
which is again correct to first order in Mej /Mt . Note that
because of the factor of 2/q, the ratio da/a can be more than
an order of magnitude larger than Mej /M1 .
≡−
Submitted to Astrophysical Journal, February 19, 2015
Preprint typeset using LATEX style emulateapj v. 01/23/15
THE ROTATIONAL BEHAVIOR OF KEPLER STARS WITH PLANETS
´ n1,4,† , I. C. Lea
˜ o1 , J. P. Bravo1 , D. B. de Freitas1 , C. E. Ferreira Lopes1,3 , S. Alves2, 5 , M.
F. Paz-Chincho
Catelan2,4, 6 , B. L. Canto Martins1 , and J. R. De Medeiros1
1 Departamento
de F´ısica Te´
orica e Experimental, Universidade Federal do Rio Grande do Norte, Campus Universit´
ario, Natal, RN,
Brazil
de Astrof´ısica, Pontificia Universidad Cat´
olica de Chile, Av. Vicu˜
na Mackenna 4860, Macul, Santiago, Chile
3 SUPA Wide–Field Astronomy Unit, Institute for Astronomy, School of Physics and Astronomy, University of Edinburgh, Royal
Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK and
4 Millennium Institute of Astrophysics (MAS), Santiago, Chile
Submitted to Astrophysical Journal, February 19, 2015
arXiv:1502.05051v1 [astro-ph.SR] 17 Feb 2015
2 Instituto
ABSTRACT
We analyzed the host stars of the present sample of confirmed planets detected by Kepler and Kepler
Objects of Interest (KOI) to compute new photometric rotation periods and to study the behavior
of their angular momentum. Lomb–Scargle periodograms and wavelet maps were computed for 3, 807
stars. For 540 of these stars, we were able to detect rotational modulation of the light curves at
a significance level of greater than 99%. For 63 of these 540 stars, no rotation measurements were
previously available in the literature. According to the published masses and evolutionary tracks of the
stars in this sample, the sample is composed of M– to F–type stars (with masses of 0.48–1.53 M⊙ ) with
rotation periods that span a range of 2 to 89 days. These periods exhibit an excellent agreement with
previously reported (for the stars for which such values are available), and the observed rotational
period distribution strongly agrees with theoretical predictions. Furthermore, for the 540 sources
considered here, the stellar angular momentum provides an important test of Kraft’s relation based
on the photometric rotation periods. Finally, this study directly contributes in a direct approach to
our understanding of how angular momentum is distributed between the host star and its (detected)
planetary system; the role of angular momentum exchange in such systems is an unavoidable piece of
the stellar rotation puzzle.
Keywords: — stars:— general — Hertzsprung-Russell and C-M diagrams — stars: kinematics and
dynamics — stars: rotation
1. INTRODUCTION
The study of stellar rotation is important for our understanding of many kinematic properties of stars, including the origin and evolution of the rotation itself;
more generally, the origin and evolution of stellar angular momentum; and the role of rotation in stellar
magnetism and mixing of chemical elements. Stars are
formed from rotating molecular cloud cores and preserve only a very small fraction of the initial angular
momentum of these cores (e.g., Palla 2002; Lamm et al.
2005). Indeed, the initial angular momentum of the
initial molecular cloud is approximately eight orders of
magnitude greater than that of the stars that eventually form from such a cloud core (e.g., Bodenheimer
1989). With the discovery of extra–solar planets (e.g.,
Mayor & Queloz 1995), new questions concerning rotation have begun to arise, including questions regarding the large–scale effects of the presence of planetary
companions on stellar rotation (Pont 2009; Alves et al.
2010; Lanza 2010) and on the heating of stellar coronae and chromospheres (e.g., Poppenhaeger et al. 2010;
Canto Martins et al. 2011; Lanza 2012). Admittedly, the
study of the rotation of planet–hosting stars is hampered by one major factor: for the great majority of
5 CAPES Foundation, Ministry of Education of Brazil,
Bras´ılia, DF, Brazil
6 Centro de Astro–Ingenier´
ıa, Ponticia Universidad Cat´
olica
de Chile, Santiago, Chile
† Correspondence author: [email protected]
stars with planets only the projected rotational velocity v sin i is available, and all such bodies identified,
to date are essentially slow–to–very–moderate rotators
(Alves et al. 2010). Despite this constraint, Alves et
al. have noted that stars without any detected planets
demonstrate a clear angular momentum deficit compared
with stars with planets.
The recent space–missions CoRoT (Baglin et al. 2007)
and Kepler (Koch et al. 2010) are providing a new window into the study of the behavior of stellar rotation,
including the rotation of stars with planets, by virtue
of their high–precision photometry, which yields excellent quality light curves (LCs). In addition to their
high precision, these missions provide nearly uninterrupted photometric measurements of unprecedented duration and cadence, allowing for refined analysis of the
features involved in the behavior of these temporal series, as in de Freitas et al. (2013b). In recent works,
McQuillan et al. (2013a,b, 2014), Nielsen et al. (2013),
Walkowicz & Basri (2013), and Reinhold et al. (2013)
computed the rotation periods (Prot ) for a large fraction
of the Kepler targets, including the Kepler Objects of Interest (hereafter KOI). McQuillan et al. (2013a,b, 2014)
analyzed the first 11 quarters of Kepler data by computing Fourier periodograms and autocorrelation functions
(ACF), and validated their results through visual inspection. McQuillan et al. (2013a) computed the rotation periods of 1, 570 M dwarfs and identified a possible change
in the slope of the period–mass relation at approximately
0.55 M⊙ . McQuillan et al. (2013b) measured the periods
Mon. Not. R. Astron. Soc. 000, 1–18 (0000)
Printed 19 February 2015
(MN LATEX style file v2.2)
arXiv:1502.05050v1 [astro-ph.GA] 17 Feb 2015
The LMC geometry and outer stellar populations
from early DES data
Eduardo Balbinot1,2,3? , B. X. Santiago2,3 , L. Girardi4 , A. Pieres2,3 , L. N. da Costa3,5 ,
M. A. G. Maia3,5 , R. A. Gruendl6,7 A. R. Walker8 , B. Yanny9 , A. Drlica-Wagner9 ,
A. Benoit-Levy10 , T. M. C. Abbott8 , S. S. Allam9,11 , J. Annis9 , J. P. Bernstein12 ,
R. A. Bernstein13 , E. Bertin14 , D. Brooks10 , E. Buckley-Geer9 , A. Carnero Rosell3,5 ,
C. E. Cunha15 , D. L. DePoy16 , S. Desai17,18 , H. T. Diehl9 , P. Doel10 , J. Estrada9 ,
A. E. Evrard19,20,14 , A. Fausti Neto3 , D. A. Finley9 , B. Flaugher9 , J. A. Frieman9,18 ,
D. Gruen17,21 , K. Honscheid22 , D. James8 , K. Kuehn23 , N. Kuropatkin9 , O. Lahav10 ,
M. March24 , J. L. Marshall16 , C. Miller19,20 , R. Miquel25,26 R. Ogando3,5 , J. Peoples9 ,
A. Plazas27 , V. Scarpine9 , M. Schubnell12 , I. Sevilla-Noarbe28 , R. C. Smith8 ,
M. Soares-Santos9 , E. Suchyta22 , M. E. C. Swanson7 , G. Tarle19 , D. L. Tucker9 ,
R. Wechsler29,18 , J. Zuntz30
(Affiliations can be found after the references)
19 February 2015
ABSTRACT
The Dark Energy Camera has captured a large set of images as part of Science Verification (SV) for the Dark Energy Survey. The SV footprint covers a large portion of the
outer Large Magellanic Cloud (LMC), providing photometry 1.5 magnitudes fainter
than the main sequence turn-off of the oldest LMC stellar population. We derive geometrical and structural parameters for various stellar populations in the LMC disk. For
the distribution of all LMC stars, we find an inclination of i = −38.14◦ ± 0.08◦ (near
side in the North) and a position angle for the line of nodes of θ0 = 129.51◦ ±0.17◦ . We
find that stars younger than ∼ 4 Gyr are more centrally concentrated than older stars.
Fitting a projected exponential disk shows that the scale radius of the old populations
is R>4Gyr = 1.41 ± 0.01 kpc, while the younger population has R<4Gyr = 0.72 ± 0.01
kpc. However, the spatial distribution of the younger population deviates significantly
from the projected exponential disk model. The distribution of old stars suggests a
large truncation radius of Rt = 13.5 ± 0.8 kpc. If this truncation is dominated by the
tidal field of the Galaxy, we find that the LMC is ' 24+9
−6 times less massive than the
encircled Galactic mass. By measuring the Red Clump peak magnitude and comparing
with the best-fit LMC disk model, we find that the LMC disk is warped and thicker
in the outer regions north of the LMC centre. Our findings may either be interpreted
as a warped and flared disk in the LMC outskirts, or as evidence of a spheroidal halo
component.
Key words: galaxies: Magellanic Clouds; galaxies: stellar content; stars: statistics
1
INTRODUCTION
The Milky Way (MW) satellite system offers a variety of examples of dwarf galaxies. Most of its members are essentially
?
e-mail: [email protected]
c 0000 RAS
gas-free and contain mainly old stars (McConnachie 2012).
The evolution of these systems is closely related to the formation of the Galaxy and the process of mass assembly of
the large scale structures in the Universe (Klypin et al. 1999;
Moore et al. 1999; Stewart et al. 2008). On the other hand,
the Large and Small Magellanic Clouds (LMC and SMC,
Mon. Not. R. Astron. Soc. 000, 1–11 (2002)
Printed February 19, 2015
(MN LATEX style file v2.2)
arXiv:1502.05046v1 [astro-ph.HE] 17 Feb 2015
Estimating the fossil disc mass during supermassive black hole
mergers: the importance of torque implementation
M.
Tazzari1,2,3⋆ and G. Lodato3
1
European Southern Observatory (ESO), Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany
Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany
3 Dipartimento di Fisica, Università degli Studi di Milano, Via Celoria 16, Milano, I-20133, Italy
2 Excellence
ABSTRACT
In this paper, we revisit the issue of estimating the “fossil” disc mass in the circumprimary
disc, during the merger of a supermassive black hole binary. As the binary orbital decay speeds
up due to the emission of gravitational waves, the gas in the circumprimary disc might be
forced to accrete rapidly and could in principle provide a significant electromagnetic counterpart to the gravitational wave emission. Since the luminosity of such flare is proportional
to the gaseous mass in the circumprimary disc, estimating such mass accurately is important.
Previous investigations of this issue have produced contradictory results, with some authors
estimating super-Eddington flares and large disc mass, while others suggesting that the “fossil” disc mass is very low, even less than a Jupiter mass. Here, we perform simple 1D calculations to show that such very low estimates of the disc mass are an artifact of the specific
implementation of the tidal torque in 1D models. In particular, for moderate mass ratios of
the binary, the usual formula for the torque used in 1D models significantly overestimates the
width of the gap induced by the secondary and this artificially leads to a very small leftover
circumprimary disc. Using a modified torque, calibrated to reproduce the correct gap width
as estimated by 3D models, leads to fossil disc masses of the order of one solar mass. The
rapid accretion of the whole circumprimary disc would produce peak luminosities of the order of 1-20 times the Eddington luminosity. Even if a significant fraction of the gas escapes
accretion by flowing out the secondary orbit during the merger (an effect not included in
our calculations), we would still predict close to Eddington luminosities that might be easily
detected.
Key words: accretion, accretion discs – black hole physics – hydrodynamics – gravitational
waves – galaxies: formation
1 INTRODUCTION
The merger of two supermassive black holes (SMBH) is estimated
to be one of the most intense events of gravitational wave (GW)
emission and the ability to detect one such event is of paramount
importance not only to test General Relativity directly but also to
provide constraints on galaxy formation.
There is now strong evidence for the presence of SMBHs with
masses between 106 and 109 M⊙ in the nuclei of the local Universe
galaxies (Kormendy & Richstone 1995; Magorrian et al. 1998; Ferrarese & Ford 2005). In a hierarchical scenario for the galaxy evolution, nearby galaxies are the outcome of several mergers between
progenitors of smaller mass. If each of the progenitors of a galaxy
merger contains a SMBH in its centre the resulting galaxy will naturally host a pair of SMBHs.
It is expected that dynamical friction between the binary and
⋆
E-mail: [email protected].
the crowded field of stars in which they are embedded can remove
angular momentum from the binary and reduce its orbital separation down to parsec scales (Begelman, Blandford & Rees 1980),
at which the process is expected to stall due to the depletion of
the loss cone (Milosavljevi´c & Merritt 2001). Galaxy scale 3D hydrodynamical simulations (Escala et al. 2005; Cuadra et al. 2009;
Dotti et al. 2006) have shown that further shrinkage of the binary
due to angular momentum loss to the gaseous background can lead
to orbital separations of the order of 0.1 pc. At separations of the
order of 0.001 pc, the main driver of orbital decay is the emission of
GW. The evolution of the system between 0.1 and 0.001 pc is much
harder to study. 3D and 2D simulations are appropriate to study the
dynamics of binary-disc interaction on short timescales, and are
thus well suited to study the final stages of the orbital decay, before
and during the GW-driven phase (D’Orazio, Haiman & MacFadyen
2013; Farris et al. 2014). The overall evolution of the system during
the long disc-driven decay is best followed using simple 1D diffusion models for the disc (Armitage & Natarajan 2002; Lodato et al.
Mon. Not. R. Astron. Soc. 000, 1–?? (2014)
Printed 19 February 2015
(MN LATEX style file v2.2)
arXiv:1502.05045v1 [astro-ph.CO] 17 Feb 2015
Group-galaxy correlations in redshift space as a probe of
the growth of structure
F.
G. Mohammad?1,2 , S. de la Torre3 , D. Bianchi4,1 , L. Guzzo1 , J. A. Peacock5
1
INAF-Osservatorio astronomico di Brera, via Emilio Bianchi 46, I-23807 Merate (LC), Italy
dell’Insubria, Dipartimento di Scienza e Alta Tecnologia, Via Valleggio 11, I22100 Como, Italy
3 Aix-Marseille Universit´
e, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388, Marseille, France
4 ICG - Institute of Cosmology and Gravitation, University of Portsmouth, Dennis Sciama Building, Burnaby Road,
Portsmouth, PO1 3FX, United Kingdom
5 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, United Kingdom
2 Universit`
a
19 February 2015
ABSTRACT
We investigate the use of the cross-correlation between galaxies and galaxy groups to
measure redshift-space distortions (RSD) and thus probe the growth rate of cosmological structure. This is compared to the classical approach based on using galaxy
auto-correlation. We make use of realistic simulated galaxy catalogues that have been
constructed by populating simulated dark matter haloes with galaxies through halo
occupation prescriptions. We adapt the classical RSD dispersion model to the case of
the group-galaxy cross-correlation function and estimate the RSD parameter β by fitting both the full anisotropic correlation function ξ s (rp , π) and its multipole moments.
In addition, we define a modified version of the latter statistics by truncating the multipole moments to exclude strongly non-linear distortions at small transverse scales.
We fit these three observable quantities in our set of simulated galaxy catalogues and
estimate statistical and systematic errors on β for the case of galaxy-galaxy, groupgroup, and group-galaxy correlation functions. When ignoring off-diagonal elements
of the covariance matrix in the fitting, the truncated multipole moments of the groupgalaxy cross-correlation function provide the most accurate estimate, with systematic
errors below 3% when fitting transverse scales larger than 10 h−1 Mpc. When including
the full covariance matrix, however, the three observables perform more similarly and
are more stable with respect to the included scales. Group auto-correlation provides
marginally smaller systematic errors, followed by group-galaxy cross-correlation and
galaxy auto-correlation. Although statistical errors are generally larger for groups, the
use of group-galaxy cross-correlation can potentially allow the reduction of systematics
while using simple linear or dispersion models.
Key words: Cosmology: large-scale structure of Universe – Galaxies: statistics.
1
INTRODUCTION
Cosmological observations over the past 15 years have established that the Universe is undergoing a late-time phase
of accelerated expansion (e.g. Riess et al. 1998; Perlmutter
et al. 1999). This is most naturally explained with the presence of a ‘dark energy’, a fluid with negative pressure filling the entire Universe, which is currently indistinguishable
from Einstein’s cosmological constant. Alternatively, however, one may consider reproducing observations by modifying the very nature of the gravitational equations of General
Relativity (e.g. Carroll et al. 2004; Joyce et al. 2014). Geo-
?
E-mail: [email protected]
metrical probes such as the Cosmic Microwave Background
(CMB), Baryon Acoustic Oscillations (BAO) and type 1a
Supernovae (SN1a) constrain the expansion history H(z),
which however can be equally well fitted by both scenarios.
This degeneracy between dark energy and a modification
of standard gravity can only be lifted by dynamical probes
looking at the growth of structure inside the Universe, which
is directly linked to the underlying theory of gravity (Joyce
et al. 2014).
The growth of cosmological structure induces galaxy peculiar velocities, i.e. coherent flows of galaxies towards matter overdensities. When redshifts are used to map galaxy
positions, we are sensitive to such peculiar velocities, whose
line-of-sight component combines with the cosmological red-
Draft version February 19, 2015
Preprint typeset using LATEX style emulateapj v. 5/2/11
CONNECTING DARK MATTER HALOS WITH
THE GALAXY CENTER AND THE SUPERMASSIVE BLACK HOLE
´
´ n1 and Andy D. Goulding1,2
Akos
Bogda
arXiv:1502.05043v1 [astro-ph.GA] 17 Feb 2015
1 Harvard
Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA; [email protected] and
2 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA
Draft version February 19, 2015
ABSTRACT
Observational studies of nearby galaxies have demonstrated correlations between the mass of the
central supermassive black holes (BHs) and properties of the host galaxies, notably the stellar bulge
mass or central stellar velocity dispersion. Motivated by these correlations, the theoretical paradigm
has emerged, in which BHs and bulges co-evolve. However, this picture was challenged by observational
and theoretical studies, which hinted that the fundamental connection may be between BHs and dark
matter halos, and not necessarily with their host galaxies. Based on a study of 3130 elliptical galaxies
– selected from the Sloan Digital and ROSAT All Sky Surveys – we demonstrate that the central
stellar velocity dispersion exhibits a significantly tighter correlation with the total gravitating mass,
traced by the X-ray luminosity of the hot gas, than with the stellar mass. This hints that the central
stellar velocity dispersion, and hence the central gravitational potential, may be the fundamental
property of elliptical galaxies that is most tightly connected to the larger-scale dark matter halo.
Furthermore, using the central stellar velocity dispersion as a surrogate for the BH mass, we find that
in elliptical galaxies the inferred BH mass and inferred total gravitating mass within the virial radius
1.6+0.6
1.8+0.7
(or within five effective radii) can be expressed as MBH ∝ Mtot −0.4 (or MBH ∝ M5reff−0.6 ). These
results are consistent with a picture in which the BH mass is directly set by the central stellar velocity
dispersion, which, in turn, is determined by the total gravitating mass of the system.
Subject headings: galaxies: elliptical and lenticular, cD, – galaxies: evolution, – galaxies: halos, –
X-rays: galaxies, – X-rays: ISM
1. INTRODUCTION
Supermassive black holes (BHs), located at the center
of every massive galaxy, are believed to have a profound
effect on the evolution of their host galaxies. Rapidly
growing BHs shine as active galactic nuclei (AGNs), and
can release copious amounts of energy into their surroundings, which can heat and expel their gas supply,
thereby quenching star formation and the BH growth
(Silk & Rees 1998; King 2003; Wyithe & Loeb 2003;
Hopkins et al. 2006). Energetic feedback from BHs was
suggested to play a pivotal role in building the presentday observed galaxy luminosity function (Croton et al.
2006; Faber et al. 2007), and hence BHs are a key component of modern galaxy formation theories.
Observational studies of nearby galaxies have pointed
out the existence of scaling relations between the BH
mass (MBH ) and various properties of the host galaxies, in particular with the stellar bulge mass, Mbulge
(Magorrian et al. 1998; H¨
aring & Rix 2004); or central stellar velocity dispersion, σc (Gebhardt et al. 2000;
Ferrarese & Merritt 2000; G¨
ultekin et al. 2009). To explain the existence of these relatively tight scaling relations, a widely accepted theoretical paradigm has been
established, in which BHs grow in a self-regulated manner and co-evolve with their host bulges from early
epochs to the present day (Silk & Rees 1998; Fabian
1999; King 2003; Hopkins et al. 2006). However, several observational and theoretical studies revealed correlations between the BH mass and the total gravitating mass, thereby hinting that dark matter halos may
govern the BH growth (Ferrarese 2002; Bandara et al.
2009; Booth & Schaye 2010, 2011; Bogd´an et al. 2012).
Follow-up studies demonstrated that in late-type or
bulgeless galaxies the dark matter halo and the BH
mass do not correlate tightly (Kormendy & Bender 2011;
Sun et al. 2013). The absence of a tight relation in
these galaxies may be due to low-mass or disk-dominated
galaxies being less efficient at growing BHs, and dynamical effects that increase the scatter in the MBH − Mtot relation for low-mass systems (Volonteri et al. 2011). However, observationally it is not well established, whether
dark matter halos play a major role in driving the BH
growth in elliptical galaxies. Therefore, to fully explore
the connection between the BH mass and total halo mass,
it is essential to perform a dedicated study focusing on
elliptical galaxies.
The main difficulty in conclusively determining
whether the growth of BHs in elliptical galaxies is driven
by stellar bulges or their large-scale dark matter halos
arises from the demanding nature of BH and halo mass
measurements. Specifically, dynamical BH mass measurements have only been obtained for ∼ 80 galaxies,
with about half of these BHs being located in elliptical
galaxies (G¨
ultekin et al. 2009; McConnell & Ma 2013).
Of these systems, only a handful have precisely measured
halo masses (Humphrey et al. 2006), which is insufficient
to probe the importance of dark matter halos. Additionally, the local sample of galaxies with available BH and
halo mass measurements does not offer a representative
sample because galaxies residing in clusters are overrepresented and include several brightest cluster galaxies.
Therefore, instead of using individual galaxies, we may
Why CLEAN when you can PURIFY? A new approach for
next-generation radio-interferometric imaging
arXiv:1502.05037v1 [astro-ph.IM] 17 Feb 2015
Rafael E. Carrillo∗ , Jason D. McEwen† , and Yves Wiaux‡
∗ Institute of Electrical Engineering, Ecole Polytechnique F´ed´erale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.
† Mullard Space Science Laboratory, University College London (UCL), Holmbury St Mary, Surrey RH5 6NT, UK.
‡ Institute of Sensors, Signals, and Systems, Heriot-Watt University, Edinburgh EH14 4AS, UK.
Abstract—In recent works, sparse models and convex optimization
techniques have been applied to radio-interferometric (RI) imaging
showing the potential to outperform state-of-the-art imaging algorithms
in the field. In this talk, I will review our latest contributions in RI
imaging, which leverage the versatility of convex optimization to both
handle realistic continuous visibilities and offer a highly parallelizable
structure paving the way to high-dimensional data scalability. Firstly,
I will review our recently proposed average sparsity approach, SARA,
which relies on the observation that natural images exhibit strong average
sparsity over multiple coherent bases. Secondly, I will discuss efficient
implementations of SARA, and sparse regularization problems in general,
for large-scale imaging problems in a new toolbox dubbed PURIFY.
The advent of next-generation radio telescopes, such as the new
LOw Frequency ARray (LOFAR), the recently upgraded Karl G.
Jansky Very Large Array (VLA) and the future Square Kilometer
Array (SKA), has posed several challenges for image reconstruction
and the design of data processing systems [1]. The new telescopes
will achieve much higher dynamic range than current instruments at
a higher angular resolution. Also, these telescopes will acquire a
massive amount of data, thus posing large-scale inverse problems in
the perspective of image reconstruction. These challenges have triggered an intense research in the community to reformulate imaging
and calibration techniques for radio interferometry (RI).
The RI measurement equation can be discretized as y = Φx + n,
where y ∈ CM denotes the vector of measured visibilities, Φ ∈
CM ×N is a discretization of the measurement operator and n ∈ CM
represents the observation noise. In [2] we propose an imaging
algorithm dubbed sparsity averaging reweighted analysis (SARA)
based on average sparsity over multiple bases, showing superior
reconstruction qualities relative to state-of-the-art imaging methods
in the field. A sparsity dictionary composed of a concatenation
of q coherent bases, Ψ = [Ψ1 , Ψ2 , . . . , Ψq ], is used and average
sparsity is promoted through the minimization of an analysis ℓ0 prior,
¯ 0 , where Ψ† denotes the adjoint operator of Ψ [3].
kΨ† xk
SARA adopts a reweighted ℓ1 minimization scheme to promote
¯ 0 . The algorithm replaces
average sparsity through the prior kΨ† xk
the ℓ0 norm by a weighted ℓ1 norm and solves a sequence of weighted
ℓ1 problems where the weights are essentially the inverse of the values
of the solution of the previous problem [2], [3]. The weighted ℓ1
problem is defined as:
¯ 1 subject to ky − Φxk
¯ 2 ≤ ǫ,
min kWΨ† xk
N
¯
x∈R
+
(1)
where W ∈ RD×D denotes the diagonal matrix with positive weights,
N
RN
and ǫ is an upper bound on
+ denotes the positive orthant in R
the ℓ2 norm of the noise, which can be accurately estimated. Hence,
we focus our attention on solving problem (1) efficiently, especially
for large-scale data problems, i.e. when the number of visibilities
is very large (M ≫ N ). In this case, we propose to split the data
vector y and the measurement operator into R blocks in the following
manner:
T T
y = [y1T , . . . , yR
] and Φ = [ΦT1 , . . . , ΦTR ]T ,
(2)
where each yi is modelled as yi = Φi x + ni and ni denotes the
noise vector. With this partition the optimization problem in (1) can
be reformulated as
¯ 2 ≤ ǫi , i = 1, . . . , R, (3)
¯ 1 subject to kyi − Φi xk
min kWΨ† xk
N
¯
x∈R
+
where each ǫi is an appropriate bound for the ℓ2 norm of ni .
In [4] we propose a general algorithmic framework based on the
simultaneous-direction method of multipliers (SDMM) to solve (3).
The proposed framework offers a parallel implementation structure
that decomposes the original problem into several small simple
problems, hence allowing implementation in multicore architectures
or in computer clusters, or on graphics processing units. These
implementations provide both flexibility in memory requirements
and a significant gain in terms of speed, thus enabling scalability to large-scale problems. A beta version of an SDMM-based
imaging software written in C and dubbed PURIFY was released
that handles various sparsity priors, including SARA, thus providing
a new powerful framework for RI imaging (toolbox available at
http://basp-group.github.io/purify/). Even though this beta version of
PURIFY is not parallelized yet, we discuss in detail the extraordinary parallel and distributed optimization potential of SDMM to be
exploited in future versions. We also discuss other possible research
avenues for big data scalability. One possibility is to incorporate
ideas from stochastic gradient methods into proximal splitting and
augmented Lagrangian methods. The key idea is to use only one data
block yi , or a subset of blocks, at each iteration of the reconstruction.
By doing so, the computational complexity per iteration will be
reduced. Thus, the total processing time of the algorithm will also be
reduced if the convergence rate of the original problem is preserved.
See for example [5] and references therein for first theoretical results.
R EFERENCES
[1] S. Wijnholds, A.-J. van der Veen, F. De Stefani, E. La Rosa, and
A. Farina, “Signal processing challenges for radio astronomical arrays,” in
Proceedings of the IEEE International Conference on Acoustics, Speech
and Signal Processing (ICASSP), 2014, pp. 5382–5386.
[2] R. E. Carrillo, J. D. McEwen, and Y. Wiaux, “Sparsity averaging
reweighted analysis (SARA): a novel algorithm for radio-interferometric
imaging,” MNRAS, vol. 426, no. 2, pp. 1223–1234, 2012.
[3] R. E. Carrillo, J. D. McEwen, D. V. D. Ville, J.-P. Thiran, and Y. Wiaux,
“Sparsity averaging for compressive imaging,” IEEE Signal Process.
Letters, vol. 20, no. 6, pp. 591–594, 2013.
[4] R. E. Carrillo, J. D. McEwen, and Y. Wiaux, “PURIFY: a new approach
to radio-interferometric imaging,” MNRAS, vol. 439, no. 4, pp. 3591–
3604, 2014.
[5] S. Azadi and S. Sra, “Towards an optimal stochastic alternating direction
method of multipliers,” in Proceedings of the 31st International Conference on Machine Learning, 2014, pp. 620–628.