sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Haas, M. R.

Normalized to: Haas, M.

170 article(s) in total. 893 co-authors, from 1 to 50 common article(s). Median position in authors list is 5,0.

[1]  [pdf] - 2124861
Time-delay measurement of MgII broad line response for the highly-accreting quasar HE 0413-4031: Implications for the MgII-based radius-luminosity relation
Comments: 39 pages (23 pages - Main text, 16 pages - Appendix), 21 figures, 14 tables; accepted for publication in the Astrophysical Journal
Submitted: 2020-05-18
We present the monitoring of the AGN continuum and MgII broad line emission for the quasar HE 0413-4031 ($z=1.38$) based on the six-year monitoring by the South African Large Telescope (SALT). We managed to estimate a time-delay of $302.6^{+28.7}_{-33.1}$ days in the rest frame of the source using seven different methods: interpolated cross-correlation function (ICCF), discrete correlation function (DCF), $z$-transformed DCF, JAVELIN, two estimators of data regularity (Von Neumann, Bartels), and $\chi^2$ method. This time-delay is below the value expected from the standard radius-luminosity relation. However, based on the monochromatic luminosity of the source and the SED modelling, we interpret this departure as the shortening of the time-delay due to the higher accretion rate of the source, with the inferred Eddington ratio of $\sim 0.4$. The MgII line luminosity of HE 0413-4031 responds to the continuum variability as $L_{\rm line}\propto L_{\rm cont}^{0.43\pm 0.10}$, which is consistent with the light-travel distance of the location of MgII emission at $R_{\rm out} \sim 10^{18}\,{\rm cm}$. Using the data of 10 other quasars, we confirm the radius-luminosity relation for broad MgII line, which was previously determined for broad H$\beta$ line for lower-redshift sources. In addition, we detect a general departure of higher-accreting quasars from this relation in analogy to H$\beta$ sample. After the accretion-rate correction of the light-travel distance, the MgII-based radius-luminosity relation has a small scatter of only $0.10$ dex.
[2]  [pdf] - 2082424
Dust Reverberation of 3C273: torus structure and lag -- luminosity relation
Comments: 24 pages, 19 Figures, accepted for publication in the Astronomical Journal
Submitted: 2020-04-21
We monitored the z=0.158 quasar 3C273 between 2015 and 2019 in the optical (BVrz) and near-infrared (JHK) with the aim to perform dust reverberation mapping. Accounting for host galaxy and accretion disk contributions, we obtained pure dust light curves in JHK. Cross correlations between the V-band and the dust light curves yield an average rest-frame delay for the hot dust of tau_cent ~ 410 days. This is a factor 2 shorter than expected from the the dust ring radius R_x ~ 900 light days reported from interferometric studies. The dust covering factor (CF) is about 8%, much smaller than predicted from the half covering angle of 45 degree found for active galactic nuclei (AGN). We analyse the asymmetric shape of the correlation functions and explore whether an inclined bi-conical bowl-shaped dust torus geometry could bring these findings (tau_cent, R_x and CF) into a consistent picture. The hot varying dust emission originates from the edge of the bowl rim with a small covering angle 40 < theta < 45 degree, and we see only the near side of the bi-conus. Such a dust gloriole with R_x = 900 +/- 200 ld and an inclination i = 12 degree matches the data remarkably well. Comparing the results of 3C273 with literature for less luminous AGN, we find a lag-luminosity relation tau propto L^alpha with alpha between 0.33 and 0.40, flatter than the widely adopted relation with alpha = 0.5. We address several explanations for the new lag - luminosity relation.
[3]  [pdf] - 2026180
Spectral evolution and radial dust transport in the prototype young eruptive system EX Lup
Comments: 12 pages, 4 figures, accepted for publication in ApJ
Submitted: 2019-10-29
EX Lup is the prototype of a class of pre-main sequence eruptive stars defined by their repetitive outbursts lasting several months. In 2008 January-September EX Lup underwent its historically largest outburst, brightening by about 4 magnitudes in visual light. In previous studies we discovered on-going silicate crystal formation in the inner disk during the outburst, but also noticed that the measured crystallinity fraction started decreasing after the source returned to the quiescent phase. Here we present new observations of the 10 $\mu$m silicate feature, obtained with the MIDI and VISIR instruments at Paranal Observatory. The observations demonstrate that within five years practically all crystalline forsterite disappeared from the surface of the inner disk. We reconstruct this process by presenting a series of parametric axisymmetric radiative transfer models of an expanding dust cloud that transports the crystals from the terrestrial zone to outer disk regions where comets are supposed to form. Possibly the early Sun also experienced similar flare-ups, and the forming planetesimals might have incorporated crystalline silicate material produced by such outbursts. Finally, we discuss how far the location of the dust cloud could be constrained by future JWST observations.
[4]  [pdf] - 2025865
The New EXor Outburst of ESO-H-alpha~99 observed by Gaia ATLAS and TESS
Comments: Accepted by The Astronomical Journal
Submitted: 2019-09-23
We report photometry and spectroscopy of the outburst of the young stellar object ESO-Halpha 99. The outburst was first noticed in Gaia alert Gaia18dvc and later by ATLAS. We have established the outburst light curve with archival ATLAS ``Orange'' filter photometry, Gaia data, new V-band photometry, and J, H, and K_s photometry from IRIS and UKIRT. The brightness has fluctuated several times near the light curve maximum. The TESS satellite observed ESO-Halpha 99 with high cadence during one of these minor minima and found brightness fluctuations on timescales of days and hours. Imaging with UKIRT shows the outline of an outflow cavity, and we find one knot of H_2 1-0 S(1) emission, now named MHO 1520, on the symmetry axis of this nebula, indicating recent collimated outflow activity from ESO-Halpha 99. Its pre-outburst SED shows a flat FIR spectrum, confirming its early evolutionary state and its similarity to other deeply embedded objects in the broader EXor class. The pre-outburst luminosity is 34 L_sun, a much higher luminosity than typical EXors, indicating that ESO-Halpha 99 may be a star of intermediate mass. Infrared and optical spectroscopy show a rich emission line spectrum, including HI lines, strong red CaII emission, as well as infrared CO bandhead emission, all characteristic EXors in the broadest sense. Comparison of the present spectra with an optical spectrum obtained in 1993, presumably in the quiescent state of the object, shows that during the present outburst the continuum component of the spectrum has increased notably more than the emission lines. The Halpha equivalent width during the outburst is down to one half of its 1993 level, and shock-excited emission lines are much less prominent.
[5]  [pdf] - 1924965
Time delay measurement of Mg II line in CTS C30.10 with SALT
Comments: submitted to ApJ; comments welcome
Submitted: 2019-01-28
We report 6 yr monitoring of a distant bright quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope (SALT). We measured the rest-frame time-lag of $562\pm 2$ days between the continuum variations and the response of the Mg II emission line, using the Javelin approach. More conservative approach, based on five different methods, imply the time delay of $564^{+109}_{-71}$ days. This time delay, combined with other available measurements of Mg II line delay, mostly for lower redshift sources, shows that the Mg II line reverberation implies a radius-luminosity relation very similar to the one based on a more frequently studied H$\beta$ line.
[6]  [pdf] - 1790934
Broad-line region structure and line profile variations in the changing look AGN HE1136-2304
Comments: 21 pages, 33 figures
Submitted: 2018-08-22, last modified: 2018-11-27
A strong X-ray outburst was detected in HE1136-2304 in 2014. Accompanying optical spectra revealed that the spectral type has changed from a nearly Seyfert 2 type (1.95), classified by spectra taken 10 and 20 years ago, to a Seyfert 1.5 in our most recent observations. We seek to investigate a detailed spectroscopic campaign on the spectroscopic properties and spectral variability behavior of this changing look AGN and compare this to other variable Seyfert galaxies. We carried out a detailed spectroscopic variability campaign of HE1136-2304 with the 10 m Southern African Large Telescope (SALT) between 2014 December and 2015 July. The broad-line region (BLR) of HE1136-2304 is stratified with respect to the distance of the line-emitting regions. The integrated emission line intensities of Halpha, Hbeta, HeI 5876, and HeII 4686 originate at distances of 15.0 (+4.2,-3.8), 7.5 (+4.6,-5.7), 7.3 (+2.8,-4.4), and 3.0 (+5.3,-3.7) light days with respect to the optical continuum at 4570AA. The variability amplitudes of the integrated emission lines are a function of distance to the ionizing continuum source as well. We derived a central black hole mass of 3.8 (+-3.1) 10exp(7) M_solar based on the line widths and distances of the BLR. The outer line wings of all BLR lines respond much faster to continuum variations indicating a Keplerian disk component for the BLR. The response in the outer wings is about two light days shorter than the response of the adjacent continuum flux with respect to the ionizing continuum flux. The vertical BLR structure in HE1136-2304 confirms a general trend that the emission lines of narrow line active galactic nuclei (AGNs) originate at larger distances from the midplane in comparison to AGNs showing broader emission lines. Otherwise, the variability behavior of this changing look AGN is similar to that of other AGN.
[7]  [pdf] - 1791001
Photodynamical analysis of the triply eclipsing hierarchical triple system EPIC 249432662
Comments: 20 pages, accepted for publication in MNRAS
Submitted: 2018-09-12, last modified: 2018-11-22
Using Campaign 15 data from the K2 mission, we have discovered a triply-eclipsing triple star system: EPIC 249432662. The inner eclipsing binary system has a period of 8.23 days, with shallow $\sim$3% eclipses. During the entire 80-day campaign, there is also a single eclipse event of a third-body in the system that reaches a depth of nearly 50% and has a total duration of 1.7 days, longer than for any previously known third-body eclipse involving unevolved stars. The binary eclipses exhibit clear eclipse timing variations. A combination of photodynamical modeling of the lightcurve, as well as seven follow-up radial velocity measurements, has led to a prediction of the subsequent eclipses of the third star with a period of 188 days. A campaign of follow-up ground-based photometry was able to capture the subsequent pair of third-body events as well as two further 8-day eclipses. A combined photo-spectro-dynamical analysis then leads to the determination of many of the system parameters. The 8-day binary consists of a pair of M stars, while most of the system light is from a K star around which the pair of M stars orbits.
[8]  [pdf] - 1759444
The broad line region and dust torus structure of the Seyfert 1 WPVS48
Comments: 5 pages, 3 figures, Accepted for publication in Proceedings of Science. Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe, 9-13 April 2018, Padova, Italy
Submitted: 2018-10-02
Optical and near-mid-infrared reverberation mapping data obtained at Universit\"{a}tssternwarte Bochum in Chile and with the Spitzer Space Telescope allow us to explore the geometry of both the H$\alpha$ BLR and the dust torus for the nearby Seyfert 1 galaxy WPVS\,48. On average, the H$\alpha$ variations lag the blue AGN continuum by about 18 days, while the dust emission variations lag by 70 days in the J+K and by 90 days in the L+M bands. The IR echoes are sharp, while the H$\alpha$ echo is smeared. This together favours a bowl-shaped toroidal geometry where the dust sublimation radius is defined by a bowl surface, which is virtually aligned with a single iso-delay surface, thus leading to the sharp IR echoes. The BLR clouds, however, are located inside the bowl and spread over a range of iso-delay surfaces, leading to a smeared echo.
[9]  [pdf] - 1767677
ALMA resolves the stellar birth explosion in distant quasar 3C298
Comments: Astrophysical Journal Letters, in press
Submitted: 2018-09-22
Galaxies are believed to experience star formation and black hole driven nuclear activity symbiotically. The symbiosis may be more extreme in the distant universe, as far-infrared photometry with the Herschel Space Observatory has found many cases of ultra-luminous cool dust emission in z>1 radio galaxies and quasars, which could have its origin in the central black hole activity, or in extreme starbursts. We here present strong evidence for an extreme circumnuclear starburst in the z=1.439 quasar 3C298. Our unparalleled 0.18 arcsecond resolution ALMA image at rest-frame 410micrometer wavelength shows that the ~40K dust in its host galaxy resides in an asymmetric circumnuclear structure. The morphology of this structure implies a starburst origin and a symbiotic physical relation with the AGN driven radio source. The symbiosis is likely to be a general property of distant massive galaxies.
[10]  [pdf] - 1795758
Simultaneous H alpha and dust reverberation mapping of 3C120: Testing the bowl-shaped torus geometry
Comments: 17 pages, 18 figures, 3 tables, accepted for publication in in section 4. Extragalactic astronomy of Astronomy and Astrophysics
Submitted: 2018-08-10, last modified: 2018-08-23
At the Universitaetssternwarte Bochum near Cerro Armazones we have monitored the Seyfert-1 galaxy 3C 120 between September 2014 and March 2015 in BVRI and a narrow band filter covering the redshifted H alpha line; in addition we obtained a single con-temporary spectrum with FAST at Mt. Hopkins. Compared to earlier epochs 3C 120 is about a factor of three brighter, allowing us to study the shape of the broad line region (BLR) and the dust torus in a high luminosity phase. The analysis of the light curves yields that the dust echo is rather sharp and symmetric in contrast to the more complex broad H alpha BLR echo. We investigate how far this supports an optically thick bowl-shaped BLR and dust torus geometry as proposed by Kawaguchi & Mori (2010) and Goad et al. (2012). The comparison with several parameterizations of these models supports the following geometry: the BLR clouds lie inside the bowl closely above the bowl rim, up to a half covering angle 0 deg < theta < 40 deg (measured against the equatorial plane). Then the BLR is spread over many isodelay surfaces, yielding a smeared and structured echo as observed. Furthermore, if the BLR clouds shield the bottom of the bowl rim against radiation from the nucleus, the hot dust emission comes essentially from the top edge of the bowl (40 deg < theta < 45 deg). Then, for small inclinations as for 3C120, the top dust edge forms a ring which largely coincides with a narrow range of isodelay surfaces, yielding the observed sharp dust echo. The scale height of the BLR increases with radial distance from the black hole. This leads to luminosity dependent foreshortening effects of the lag. We discuss implications and possible corrections of the foreshortening for the black hole mass determination and consequences for the lag (size) - luminosity relationships and the difference to interferometric torus sizes.
[11]  [pdf] - 1747959
Optical and UV surface brightness of translucent dark nebulae: Dust albedo, radiation field and fluorescence emission by H2
Comments: 30 pages, 21 figures, Astronomy and Astrophysics in press
Submitted: 2018-06-16
Photometry of the nebulae LDN1780, LDN1642 and LBN406 is used to derive scattering properties of dust and to investigate the presence of UV fluorescence emission by molecular hydrogen and the extended red emission (ERE). We used multi-wavelength optical photometry and imaging at ground-based telescopes and archival imaging and spectroscopic UV data from the spaceborn GALEX and SPEAR/FIMS instruments. We used Monte Carlo RT and both observational data and synthetic models for the ISRF in the solar neighbourhood. The line-of-sight extinctions through the clouds have been determined using near infrared excesses of background stars and the 200/250um far infrared emission by dust measured using the ISO and Herschel space observatories. The optical surface brightness of the target clouds can be explained in terms of scattered light. The dust albedo ranges from 0.58 at 3500A to 0.72 at 7500A. The SED of LDN1780 is explained in terms of optical depth and background scattered light effects instead of ERE. The FUV surface brightness of LDN 1780 cannot be explained by scattered light only. In LDN1780 H2 fluorescent emission in the wavelength range 1400A-1700A has been detected and analysed. Our albedo values agree with the predictions of the dust model of Weingartner and Draine and with the THEMIS CMM model for evolved core-mantle grains. The H2 fluorescent emission in LDN1780 shows a pronounced dichotomy with a preference for its southern side where enhanced illumination impinges from the Sco OB2 association and the O star zeta Oph. A good correlation is found between the H2 fluorescence and a previously mapped 21-cm excess emission. The H2 fluorescence emission in LDN1780 has been modelled using a PDR code; the resulting values for H2 column density and the total gas density are consistent with the estimates derived from CO observations and optical extinction along the line of sight.
[12]  [pdf] - 1674933
Kepler Data Validation I -- Architecture, Diagnostic Tests, and Data Products for Vetting Transiting Planet Candidates
Comments: 84 pages, 26 figures, 2 tables. Accepted for publication in PASP on 9 March 2018. Published in PASP on 24 April 2018
Submitted: 2018-03-12, last modified: 2018-04-24
The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for ~200,000 individual stellar targets in at least one observing quarter in the primary Kepler Mission. Light curves for target stars are extracted in the Kepler Science Data Processing Pipeline, and are searched for transiting planet signatures. A Threshold Crossing Event is generated in the transit search for targets where the transit detection threshold is exceeded and transit consistency checks are satisfied. These targets are subjected to further scrutiny in the Data Validation (DV) component of the Pipeline. Transiting planet candidates are characterized in DV, and light curves are searched for additional planets after transit signatures are modeled and removed. A suite of diagnostic tests is performed on all candidates to aid in discrimination between genuine transiting planets and instrumental or astrophysical false positives. Data products are generated per target and planet candidate to document and display transiting planet model fit and diagnostic test results. These products are exported to the Exoplanet Archive at the NASA Exoplanet Science Institute, and are available to the community. We describe the DV architecture and diagnostic tests, and provide a brief overview of the data products. Transiting planet modeling and the search for multiple planets on individual targets are described in a companion paper. The final revision of the Kepler Pipeline code base is available to the general public through GitHub. The Kepler Pipeline has also been modified to support the TESS Mission which will commence in 2018.
[13]  [pdf] - 1767401
Long-term optical, UV, and X-ray continuum variations in the changing-look AGN HE 1136-2304
Comments: 25 pages, 21 figures, Astronomy & Astrophysics in press
Submitted: 2018-04-16
A strong outburst in the X-ray continuum and a change of its Seyfert spectral type was detected in HE 1136-2304 in 2014. The spectral type changed from nearly Seyfert 2 type (1.95) to Seyfert 1.5 type in comparison to previous observations taken ten to twenty years before. In a subsequent variability campaign we wanted to investigate whether this outburst was a single event or whether the variability pattern following the outburst was similar to those seen in other variable Seyfert galaxies. In addition to a SALT spectral variability campaign, we carried out optical continuum as well as X-ray and UV (Swift) monitoring studies from 2014 to 2017. HE 1136-2304 strongly varied on timescales of days to months from 2014 to 2017. No systematic trends were found in the variability behavior following the outburst in 2014. A general decrease in flux would have been expected for a tidal disruption event. This could not be confirmed. More likely the flux variations are connected to irregular fluctuations in the accretion rate. The strongest variability amplitudes have been found in the X-ray regime: HE 1136-2304 varied by a factor of eight during 2015. The amplitudes of the continuum variability (from the UV to the optical) systematically decreased with wavelength following a power law F_var = a ${\times}$ {\lambda}^-c with c = 0.84. There is a trend that the B-band continuum shows a delay of three light days with respect to the variable X-ray flux. The Seyfert type 1.5 did not change despite the strong continuum variations for the period between 2014 and 2017.
[14]  [pdf] - 1670495
Planetary Candidates Observed by Kepler. VIII. A Fully Automated Catalog With Measured Completeness and Reliability Based on Data Release 25
Comments: 61 pages, 23 Figures, 9 Tables, Accepted to The Astrophysical Journal Supplement Series
Submitted: 2017-10-18, last modified: 2018-03-04
We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs, Twicken et al. 2016). The Robovetter also vetted simulated data sets and measured how well it was able to separate TCEs caused by noise from those caused by low signal-to-noise transits. We discusses the Robovetter and the metrics it uses to sort TCEs. For orbital periods less than 100 days the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates between 200 and 500 days around FGK dwarf stars, the Robovetter is 76.7% complete and the catalog is 50.5% reliable. The KOI catalog, the transit fits and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive.
[15]  [pdf] - 1589381
Galaxy overdensities around 3C radio galaxies and quasars at $1<z<2.5$ revealed by Spitzer $3.6/4.5\mu$m and Pan-STARRS
Submitted: 2017-10-06
Luminous radio sources are thought to lie in galaxy clusters or proto-clusters. The complete sample of 64 high-redshift 3C sources at $1<z<2.5$ has been mapped with the Spitzer Space Telescope. The IRAC $3.6\mu$m and $4.5\mu$m 5-$\sigma$ detection limit of $4\mu$Jy (22.4 AB mag) allows us to search for the brightest candidate cluster member galaxies associated with the 3C sources. To remove the contamination of foreground stars and galaxies along the lines of sight toward the 3C sources we apply color cuts: removed sources satisfy either the IRAC1/2 cut $[3.6] - [4.5] < -0.3$ or the Pan-STARRS/IRAC cut $i - [4.5] < 0.5$ if detected by Pan-STARRS. For both selection methods, about half of the 3C radio sources show significant overdensities ($>3\sigma$) within 30" ($\sim$250 kpc) projected distance from the radio source compared to the surrounding galaxy densities measured in the $50''-120''$ annulus. The Pan-STARRS/IRAC cut reveals higher average overdensities than the IRAC1/2 cut and also a decline of overdensities at $z>1.4$. To infer the nature of the cluster members, we rerun the analysis using a stronger IRAC1/2 cut $[3.6] - [4.5] < -0.1$ which removes $1<z<1.4$ passive ellipticals but not star-forming galaxies. For the strong cut, the overdensities, on average, completely disappear at $1<z<1.4$. We therefore suggest that the 4.5$\mu$m detected cluster member galaxies are mainly passive ellipticals.
[16]  [pdf] - 1584134
Starburst-driven superwinds in quasar host galaxies
Comments: To appear in the Astrophysical Journal Letters
Submitted: 2017-06-01
During five decades astronomers have been puzzled by the presence of strong absorption features including metal lines, observed in the optical and ultraviolet spectra of quasars, signalling in- and outflowing gas winds with relative velocities up to several thousands of km/sec. In particular the location of these winds - close to the quasar, further out in its host galaxy, or in its direct environment - and the possible impact on their surroundings have been issues of intense discussion and uncertainty. Using our Herschel Space Observatory data, we report a tendency for this so-called associated metal absorption to occur along with prodigious star formation in the quasar host galaxy, indicating that the two phenomena are likely to be interrelated, that the gas winds likely occur on the kiloparsec scale and would then have a strong impact on the interstellar medium of the galaxy. This correlation moreover would imply that the unusually high cold dust luminosities in these quasars are connected with ongoing star formation. Given that we find no correlation with the AGN strength, the wind feedback which we establish in these radio-loud objects is most likely associated with their host star formation rather than with their black hole accretion.
[17]  [pdf] - 1546376
Time-Series Analysis of Broadband Photometry of Neptune from K2
Comments: 31 pages, 6 figure, accepted for publication in AJ
Submitted: 2017-02-08
We report here on our search for excess power in photometry of Neptune collected by the K2 mission that may be due to intrinsic global oscillations of the planet Neptune. To conduct this search, we developed new methods to correct for instrumental effects such as intrapixel variability and gain variations. We then extracted and analyzed the time-series photometry of Neptune from 49 days of nearly continuous broadband photometry of the planet. We find no evidence of global oscillations and place an upper limit of $\sim$5 ppm at 1000 \uhz\ for the detection of a coherent signal. With an observed cadence of 1-minute and point-to-point scatter less than 0.01\%, the photometric signal is dominated by reflected light from the Sun, which is in turn modulated by atmospheric variability of Neptune at the 2\% level. A change in flux is also observed due to the increasing distance between Neptune and the K2 spacecraft, and solar variability with convection-driven solar p modes present.
[18]  [pdf] - 1534015
Brightness variations of the FUor-type eruptive star V346 Nor
Comments: 7 pages, 4 figures, appendices. Accepted for publication in A&A Letters
Submitted: 2017-01-01
Decades after the beginning of its FU Orionis-type outburst, V346 Nor unexpectedly underwent a fading event of $\Delta{}K$ = 4.6 mag around 2010. We obtained near-infrared observations and re-analysed data from the VISTA/VVV survey to outline the brightness evolution. In our VLT/NaCO images, we discovered a halo of scattered light around V346~Nor with a size of about 0.04 arcsec (30 au). The VISTA data outlined a well-defined minimum in the light curve at late 2010/early 2011, and tentatively revealed a small-amplitude periodic modulation of 58 days. Our latest data points from 2016 demonstrate that the source is still brightening but has not reached the 2008 level yet. We used a simple accretion disk model with varying accretion rate and line-of-sight extinction to reproduce the observed near-infrared magnitudes and colors. We found that before 2008, the flux changes of V346 Nor were caused by a correlated change of extinction and accretion rate, while the minimum around 2010 was mostly due to decreasing accretion. The source reached a maximal accretion rate of ${\approx}10^{-4} M_{\odot}$ yr$^{-1}$ in 1992. A combination of accretion and extinction changes was already invoked in the literature to interpret the flux variations of certain embedded young eruptive stars.
[19]  [pdf] - 1515721
Detection of Potential Transit Signals in 17 Quarters of Kepler Data: Results of the Final Kepler Mission Transiting Planet Search (DR25)
Comments: 60 pages, 20 figures, 3 tables. Submitted to AAS Journals on 4/19/16. Resubmitted following referee report on 6/23/16. Accepted for publication in AJ on 6/30/16. Published in AJ on 11/11/16
Submitted: 2016-04-20, last modified: 2016-11-16
We present results of the final Kepler Data Processing Pipeline search for transiting planet signals in the full 17-quarter primary mission data set. The search includes a total of 198,709 stellar targets, of which 112,046 were observed in all 17 quarters and 86,663 in fewer than 17 quarters. We report on 17,230 targets for which at least one transit signature is identified that meets the specified detection criteria: periodicity, minimum of three observed transit events, detection statistic (i.e., signal-to-noise ratio) in excess of the search threshold, and passing grade on three statistical transit consistency tests. Light curves for which a transit signal is identified are iteratively searched for additional signatures after a limb-darkened transiting planet model is fitted to the data and transit events are removed. The search for additional planets adds 16,802 transit signals for a total of 34,032; this far exceeds the number of transit signatures identified in prior pipeline runs. There was a strategic emphasis on completeness over reliability for the final Kepler transit search. A comparison of the transit signals against a set of 3402 well-established, high-quality Kepler Objects of Interest yields a recovery rate of 99.8%. The high recovery rate must be weighed against a large number of false-alarm detections. We examine characteristics of the planet population implied by the transiting planet model fits with an emphasis on detections that would represent small planets orbiting in the habitable zone of their host stars.
[20]  [pdf] - 1494915
The Astropy Problem
Muna, Demitri; Alexander, Michael; Allen, Alice; Ashley, Richard; Asmus, Daniel; Azzollini, Ruyman; Bannister, Michele; Beaton, Rachael; Benson, Andrew; Berriman, G. Bruce; Bilicki, Maciej; Boyce, Peter; Bridge, Joanna; Cami, Jan; Cangi, Eryn; Chen, Xian; Christiny, Nicholas; Clark, Christopher; Collins, Michelle; Comparat, Johan; Cook, Neil; Croton, Darren; Davids, Isak Delberth; Depagne, Éric; Donor, John; Santos, Leonardo A. dos; Douglas, Stephanie; Du, Alan; Durbin, Meredith; Erb, Dawn; Faes, Daniel; Fernández-Trincado, J. G.; Foley, Anthony; Fotopoulou, Sotiria; Frimann, Søren; Frinchaboy, Peter; Garcia-Dias, Rafael; Gawryszczak, Artur; George, Elizabeth; Gonzalez, Sebastian; Gordon, Karl; Gorgone, Nicholas; Gosmeyer, Catherine; Grasha, Katie; Greenfield, Perry; Grellmann, Rebekka; Guillochon, James; Gurwell, Mark; Haas, Marcel; Hagen, Alex; Haggard, Daryl; Haines, Tim; Hall, Patrick; Hellwing, Wojciech; Herenz, Edmund Christian; Hinton, Samuel; Hlozek, Renee; Hoffman, John; Holman, Derek; Holwerda, Benne Willem; Horton, Anthony; Hummels, Cameron; Jacobs, Daniel; Jensen, Jens Juel; Jones, David; Karick, Arna; Kelley, Luke; Kenworthy, Matthew; Kitchener, Ben; Klaes, Dominik; Kohn, Saul; Konorski, Piotr; Krawczyk, Coleman; Kuehn, Kyler; Kuutma, Teet; Lam, Michael T.; Lane, Richard; Liske, Jochen; Lopez-Camara, Diego; Mack, Katherine; Mangham, Sam; Mao, Qingqing; Marsh, David J. E.; Mateu, Cecilia; Maurin, Loïc; McCormac, James; Momcheva, Ivelina; Monteiro, Hektor; Mueller, Michael; Munoz, Roberto; Naidu, Rohan; Nelson, Nicholas; Nitschelm, Christian; North, Chris; Nunez-Iglesias, Juan; Ogaz, Sara; Owen, Russell; Parejko, John; Patrício, Vera; Pepper, Joshua; Perrin, Marshall; Pickering, Timothy; Piscionere, Jennifer; Pogge, Richard; Poleski, Radek; Pourtsidou, Alkistis; Price-Whelan, Adrian M.; Rawls, Meredith L.; Read, Shaun; Rees, Glen; Rein, Hanno; Rice, Thomas; Riemer-Sørensen, Signe; Rusomarov, Naum; Sanchez, Sebastian F.; Santander-García, Miguel; Sarid, Gal; Schoenell, William; Scholz, Aleks; Schuhmann, Robert L.; Schuster, William; Scicluna, Peter; Seidel, Marja; Shao, Lijing; Sharma, Pranav; Shulevski, Aleksandar; Shupe, David; Sifón, Cristóbal; Simmons, Brooke; Sinha, Manodeep; Skillen, Ian; Soergel, Bjoern; Spriggs, Thomas; Srinivasan, Sundar; Stevens, Abigail; Streicher, Ole; Suchyta, Eric; Tan, Joshua; Telford, O. Grace; Thomas, Romain; Tonini, Chiara; Tremblay, Grant; Tuttle, Sarah; Urrutia, Tanya; Vaughan, Sam; Verdugo, Miguel; Wagner, Alexander; Walawender, Josh; Wetzel, Andrew; Willett, Kyle; Williams, Peter K. G.; Yang, Guang; Zhu, Guangtun; Zonca, Andrea
Submitted: 2016-10-10
The Astropy Project ( is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.
[21]  [pdf] - 1411241
The K2 Ecliptic Plane Input Catalog (EPIC) and Stellar Classifications of 138,600 Targets in Campaigns 1-8
Comments: 19 pages, 12 figures, 5 tables, accepted for publication in ApJS. An electronic version of Table 5 is available as an ancillary file (sidebar on the right), and source codes are available at and v3: minor text changes and updated uncertainties in Table 5; v4: minor text changes to match published version
Submitted: 2015-12-08, last modified: 2016-05-19
The K2 Mission uses the Kepler spacecraft to obtain high-precision photometry over ~80 day campaigns in the ecliptic plane. The Ecliptic Plane Input Catalog (EPIC) provides coordinates, photometry and kinematics based on a federation of all-sky catalogs to support target selection and target management for the K2 mission. We describe the construction of the EPIC, as well as modifications and shortcomings of the catalog. Kepler magnitudes (Kp) are shown to be accurate to ~0.1 mag for the Kepler field, and the EPIC is typically complete to Kp~17 (Kp~19 for campaigns covered by SDSS). We furthermore classify 138,600 targets in Campaigns 1-8 (~88% of the full target sample) using colors, proper motions, spectroscopy, parallaxes, and galactic population synthesis models, with typical uncertainties for G-type stars of ~3% in Teff, ~0.3 dex in log(g), ~40% in radius, ~10% in mass, and ~40% in distance. Our results show that stars targeted by K2 are dominated by K-M dwarfs (~41% of all selected targets), F-G dwarfs (~36%) and K giants (~21%), consistent with key K2 science programs to search for transiting exoplanets and galactic archeology studies using oscillating red giants. However, we find a significant variation of the fraction of cool dwarfs with galactic latitude, indicating a target selection bias due to interstellar reddening and the increased contamination by giant stars near the galactic plane. We discuss possible systematic errors in the derived stellar properties, and differences to published classifications for K2 exoplanet host stars. The EPIC is hosted at the Mikulski Archive for Space Telescopes (MAST):
[22]  [pdf] - 1483340
Measuring Transit Signal Recovery in the Kepler Pipeline. III. Completeness of the Q1-Q17 DR24 Planet Candidate Catalogue, with Important Caveats for Occurrence Rate Calculations
Comments: 8 pages, 5 figures, full electronic version of Table 1 available at the NASA Exoplanet Archive; accepted by ApJ May 2nd, 2016
Submitted: 2016-05-18
With each new version of the Kepler pipeline and resulting planet candidate catalogue, an updated measurement of the underlying planet population can only be recovered with an corresponding measurement of the Kepler pipeline detection efficiency. Here, we present measurements of the sensitivity of the pipeline (version 9.2) used to generate the Q1-Q17 DR24 planet candidate catalog (Coughlin et al. 2016). We measure this by injecting simulated transiting planets into the pixel-level data of 159,013 targets across the entire Kepler focal plane, and examining the recovery rate. Unlike previous versions of the Kepler pipeline, we find a strong period dependence in the measured detection efficiency, with longer (>40 day) periods having a significantly lower detectability than shorter periods, introduced in part by an incorrectly implemented veto. Consequently, the sensitivity of the 9.2 pipeline cannot be cast as a simple one-dimensional function of the signal strength of the candidate planet signal as was possible for previous versions of the pipeline. We report on the implications for occurrence rate calculations based on the Q1-Q17 DR24 planet candidate catalog and offer important caveats and recommendations for performing such calculations. As before, we make available the entire table of injected planet parameters and whether they were recovered by the pipeline, enabling readers to derive the pipeline detection sensitivity in the planet and/or stellar parameter space of their choice.
[23]  [pdf] - 1443954
The bimodal initial mass function in the Orion Nebula Cloud
Comments: Accepted at MNRAS, 12 pages, 13 figures, 3 tables
Submitted: 2016-05-11
Due to its youth, proximity and richness the Orion Nebula Cloud (ONC) is an ideal testbed to obtain a comprehensive view on the Initial Mass Function (IMF) down to the planetary mass regime. Using the HAWK-I camera at the VLT, we have obtained an unprecedented deep and wide near-infrared JHK mosaic of the ONC (90% completeness at K~19.0mag, 22'x28). Applying the most recent isochrones and accounting for the contamination of background stars and galaxies, we find that ONC's IMF is bimodal with distinct peaks at about 0.25 and 0.025 M_sun separated by a pronounced dip at the hydrogen burning limit (0.08 M_sun), with a depth of about a factor 2-3 below the log-normal distribution. Apart from ~920 low-mass stars (M < 1.4 M_sun) the IMF contains ~760 brown dwarf (BD) candidates and ~160 isolated planetary mass object (IPMO) candidates with M > 0.005 M_sun, hence about ten times more substellar candidates than known before. The substellar IMF peak at 0.025 M_sun could be caused by BDs and IPMOs which have been ejected from multiple systems during the early star-formation process or from circumstellar disks.
[24]  [pdf] - 1411581
False positive probabilties for all Kepler Objects of Interest: 1284 newly validated planets and 428 likely false positives
Comments: 20 pages, 8 figures. Published in ApJ. Instructions to reproduce results can be found at
Submitted: 2016-05-09
We present astrophysical false positive probability calculations for every Kepler Object of Interest (KOI)---the first large-scale demonstration of a fully automated transiting planet validation procedure. Out of 7056 KOIs, we determine that 1935 have probabilities <1% to be astrophysical false positives, and thus may be considered validated planets. 1284 of these have not yet been validated or confirmed by other methods. In addition, we identify 428 KOIs likely to be false positives that have not yet been identified as such, though some of these may be a result of unidentified transit timing variations. A side product of these calculations is full stellar property posterior samplings for every host star, modeled as single, binary, and triple systems. These calculations use 'vespa', a publicly available Python package able to be easily applied to any transiting exoplanet candidate.
[25]  [pdf] - 1521160
Campaign 9 of the $K2$ Mission: Observational Parameters, Scientific Drivers, and Community Involvement for a Simultaneous Space- and Ground-based Microlensing Survey
Henderson, Calen B.; Poleski, Radosław; Penny, Matthew; Street, Rachel A.; Bennett, David P.; Hogg, David W.; Gaudi, B. Scott; Zhu, W.; Barclay, T.; Barentsen, G.; Howell, S. B.; Mullally, F.; Udalski, A.; Szymański, M. K.; Skowron, J.; Mróz, P.; Kozłowski, S.; Wyrzykowski, Ł.; Pietrukowicz, P.; Soszyński, I.; Ulaczyk, K.; Pawlak, M.; Sumi, T.; Abe, F.; Asakura, Y.; Barry, R. K.; Bhattacharya, A.; Bond, I. A.; Donachie, M.; Freeman, M.; Fukui, A.; Hirao, Y.; Itow, Y.; Koshimoto, N.; Li, M. C. A.; Ling, C. H.; Masuda, K.; Matsubara, Y.; Muraki, Y.; Nagakane, M.; Ohnishi, K.; Oyokawa, H.; Rattenbury, N.; Saito, To.; Sharan, A.; Sullivan, D. J.; Tristram, P. J.; Yonehara, A.; Bachelet, E.; Bramich, D. M.; Cassan, A.; Dominik, M.; Jaimes, R. Figuera; Horne, K.; Hundertmark, M.; Mao, S.; Ranc, C.; Schmidt, R.; Snodgrass, C.; Steele, I. A.; Tsapras, Y.; Wambsganss, J.; Bozza, V.; Burgdorf, M. J.; Jørgensen, U. G.; Novati, S. Calchi; Ciceri, S.; D'Ago, G.; Evans, D. F.; Hessman, F. V.; Hinse, T. C.; Husser, T. -O.; Mancini, L.; Popovas, A.; Rabus, M.; Rahvar, S.; Scarpetta, G.; Skottfelt, J.; Southworth, J.; Unda-Sanzana, E.; Bryson, S. T.; Caldwell, D. A.; Haas, M. R.; Larson, K.; McCalmont, K.; Packard, M.; Peterson, C.; Putnam, D.; Reedy, L.; Ross, S.; Van Cleve, J. E.; Akeson, R.; Batista, V.; Beaulieu, J. -P.; Beichman, C. A.; Bryden, G.; Ciardi, D.; Cole, A.; Coutures, C.; Foreman-Mackey, D.; Fouqué, P.; Friedmann, M.; Gelino, C.; Kaspi, S.; Kerins, E.; Korhonen, H.; Lang, D.; Lee, C. -H.; Lineweaver, C. H.; Maoz, D.; Marquette, J. -B.; Mogavero, F.; Morales, J. C.; Nataf, D.; Pogge, R. W.; Santerne, A.; Shvartzvald, Y.; Suzuki, D.; Tamura, M.; Tisserand, P.; Wang, D.
Comments: 19 pages, 11 figures, 3 tables; submitted to PASP
Submitted: 2015-12-30, last modified: 2016-03-07
$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $\pi_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$.
[26]  [pdf] - 1395742
Star formation in 3CR radio galaxies and quasars at z < 1
Comments: Submitted to the Astronomical Journal on 2015-07-28, accepted on 2016-02-22, not all figures are included in the preprint due to size limitations but available from
Submitted: 2016-02-24, last modified: 2016-02-26
Using the Herschel Space Observatory we have observed a representative sample of 87 powerful 3CR sources at redshift z < 1. The far-infrared (FIR, 70-500 micron) photometry is combined with mid-infrared (MIR) photometry from the Wide-Field Infrared Survey Explorer (WISE) and catalogued data to analyse the complete spectral energy distributions (SEDs) of each object from optical to radio wavelength. To disentangle the contributions of different components, the SEDs are fitted with a set of templates to derive the luminosities of host galaxy starlight, dust torus emission powered by active galactic nuclei (AGN) and cool dust heated by stars. The level of emission from relativistic jets is also estimated, in order to isolate the thermal host galaxy contribution. The new data are in line with the orientation-based unification of high-excitation radio-loud AGN, in that the dust torus becomes optically thin longwards of 30 micron. The low excitation radio galaxies and the MIR weak sources represent MIR- and FIR-faint AGN population different from the high-excitation MIR-bright objects; it remains an open question whether they are at a later evolutionary state or an intrinsically different population. The derived luminosities for host starlight and dust heated by star formation are converted to stellar masses and star formation rates (SFR). The host-normalized SFR of the bulk of the 3CR sources is low when compared to other galaxy populations at the same epoch. Estimates of the dust mass yield a 1--100 times lower dust/stellar mass ratio than for the Milky Way, indicating that these 3CR hosts have very low levels of interstellar matter explaining the low level of star formation. Less than 10% of the 3CR sources show levels of star formation above those of the main sequence of star forming galaxies.
[27]  [pdf] - 1415034
Planetary Candidates Observed by Kepler. VII. The First Fully Uniform Catalog Based on The Entire 48 Month Dataset (Q1-Q17 DR24)
Comments: Accepted to the Astrophysical Journal Supplement Series. 30 pages, 9 figures, 7 tables. We make the DR24 robovetter decision code publicly available at, with input and output examples provided using the same data as contained in the full paper's tables
Submitted: 2015-12-18, last modified: 2016-02-18
We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness of individual objects, we find that our robotic vetting is overall comparable to, and in most cases is superior to, the human vetting procedures employed by past catalogs. This catalog is the first to utilize artificial transit injection to evaluate the performance of our vetting procedures and quantify potential biases, which are essential for accurate computation of planetary occurrence rates. With respect to the cumulative Kepler Object of Interest (KOI) catalog, we designate 1,478 new KOIs, of which 402 are dispositioned as planet candidates (PCs). Also, 237 KOIs dispositioned as false positives (FPs) in previous Kepler catalogs have their disposition changed to PC and 118 PCs have their disposition changed to FP. This brings the total number of known KOIs to 8,826 and PCs to 4,696. We compare the Q1-Q17 DR24 KOI catalog to previous KOI catalogs, as well as ancillary Kepler catalogs, finding good agreement between them. We highlight new PCs that are both potentially rocky and potentially in the habitable zone of their host stars, many of which orbit solar-type stars. This work represents significant progress in accurately determining the fraction of Earth-size planets in the habitable zone of Sun-like stars. The full catalog is publicly available at the NASA Exoplanet Archive.
[28]  [pdf] - 1426245
Identifying False Alarms in the Kepler Planet Candidate Catalog
Comments: 6 pages. Accepted for publication in PASP
Submitted: 2016-02-09
We present a new automated method to identify instrumental features masquerading as small, long period planets in the \kepler\ planet candidate catalog. These systematics, mistakenly identified as planet transits, can have a strong impact on occurrence rate calculations because they cluster in a region of parameter space where Kepler's sensitivity to planets is poor. We compare individual transit-like events to a variety of models of real transits and systematic events, and use a Bayesian Information Criterion to evaluate the likelihood that each event is real. We describe our technique and test its performance on simulated data. Results from this technique are incorporated in the \kepler\ Q1-17 DR24 planet candidate catalog of \citet{Coughlin15}.
[29]  [pdf] - 1298241
Light curves of the latest FUor: Indication of a close binary
Submitted: 2015-10-22
We monitored the recent FUor 2MASS J06593158-0405277 (V960 Mon) since November 2009 at various observatories and multiple wavelengths. After the outburst by nearly 2.9 mag in $r$ around September 2014 the brightness gently fades until April 2015 by nearly 1 mag in $U$ and 0.5 mag in $z$. Thereafter the brightness at $\lambda>5000 \AA$ was constant until June 2015 while the shortest wavelengths ($U, B$) indicate a new rise, similar to that seen for the FUor V2493 Cyg (HBC722). Our near-infrared (NIR) monitoring between December 2014 and April 2015 shows a smaller outburst amplitude ($\sim$2 mag) and a smaller (0.2 $-$ 0.3 mag) post-outburst brightness decline. Optical and NIR color-magnitude diagrams indicate that the brightness decline is caused by growing extinction. The post-outburst light curves are modulated by an oscillating color-neutral pattern with a period of about 17 days and an amplitude declining from $\sim$0.08 mag in October 2014 to $\sim$0.04 mag in May 2015. The properties of the oscillating pattern lead us to suggest the presence of a close binary with eccentric orbit.
[30]  [pdf] - 1292310
A Machine Learning Technique to Identify Transit Shaped Signals
Comments: Accepted for Publication to The Astrophysical Journal
Submitted: 2015-08-31
We describe a new metric that uses machine learning to determine if a periodic signal found in a photometric time series appears to be shaped like the signature of a transiting exoplanet. This metric uses dimensionality reduction and k-nearest neighbors to determine whether a given signal is sufficiently similar to known transits in the same data set. This metric is being used by the Kepler Robovetter to determine which signals should be part of the Q1-Q17 DR24 catalog of planetary candidates. The Kepler Mission reports roughly 20,000 potential transiting signals with each run of its pipeline, yet only a few thousand appear sufficiently transit shaped to be part of the catalog. The other signals tend to be variable stars and instrumental noise. With this metric we are able to remove more than 90% of the non-transiting signals while retaining more than 99% of the known planet candidates. When tested with injected transits, less than 1% are lost. This metric will enable the Kepler mission and future missions looking for transiting planets to rapidly and consistently find the best planetary candidates for follow-up and cataloging.
[31]  [pdf] - 1277170
The stability of the optical flux variation gradient for 3C120
Comments: 8 pages, 6 figures. Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics v2: Language-Editor Version
Submitted: 2015-08-07, last modified: 2015-08-27
New $B$- and $V$-band monitoring in 2014 $-$ 2015 reveals that the Seyfert 1 Galaxy, 3C120, has brightened by a magnitude of $1.4$, compared to our campaign that took place in 2009 $-$ 2010. This allowed us to check for the debated luminosity and time-dependent color variations claimed for SDSS quasars. For our 3C120 data, we find that the $B/V$ flux ratio of the variable component in the bright epoch is indistinguishable from the faint one. We do not find any color variability on different timescales ranging from about $1$ to $1800$ days. We suggest that the luminosity and time-dependent color variability is an artifact caused by analyzing the data in magnitudes instead of fluxes. The flux variation gradients of both epochs yield consistent estimates of the host galaxy contribution to our 7.5" aperture. These results confirm that the optical flux variation gradient method works well for Seyfert galaxies.
[32]  [pdf] - 1265869
Ground-based transit observations of the HAT-P-18, HAT-P-19, HAT-P-27/WASP-40 and WASP-21 systems
Submitted: 2015-08-25
As part of our ongoing effort to investigate transit timing variations (TTVs) of known exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs due to the known properties of their host systems based on the respective discovery papers. During the past three years 46 transit observations were carried out, mostly using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to refine the systems orbital parameters. In all cases we found no hints for significant TTVs, or changes in the system parameters inclination, fractional stellar radius and planet to star radius ratio. However, comparing our results with those available in the literature shows that we can confirm the already published values.
[33]  [pdf] - 1259988
The Bochum Survey of the Southern Galactic Disk: II. Follow-up measurements and multi-filter photometry for 1323 square degrees monitored in 2010-2015
Submitted: 2015-08-07
This paper is the second in a series describing the southern Galactic Disk Survey (GDS) performed at the Universit\"atssternwarte Bochum near Cerro Armazones in Chile. Haas et al. (2012, Paper I) presented the survey design and a the characteristics of the observations and data. They identified ~2200 variable stars in an area of 50 square degrees with more than 50 observations in 2011. Here we present the first complete version of the GDS covering all 268 fields with 1323 square degrees along the Galactic disk including revised data from Paper I. The individual fields were observed up to 272 times and comprise a maximum time span between September 2010 and May 2015. We detect a total of 64151 variable sources, which are presented in a catalog including some of their properties and their light curves. A comparison with the International Variable Star Index (VSX) and All Sky Automated Survey (ASAS) indicates that 56794 of these sources are previously unknown variables. Furthermore, we present UBVr'i'z' photometry for all sources within the GDS, resulting in a new multi-color catalog of nearly 16x10^6 sources detected in at least one filter. Both the GDS and the near-infrared VISTA Variables in the Via Lactea survey (VVV) complement each other in the overlap area of about 300 square degrees enabling future comparison studies.
[34]  [pdf] - 1277107
Measuring Transit Signal Recovery in the Kepler Pipeline II: Detection Efficiency as Calculated in One Year of Data
Comments: 13 pages, 7 figures, 1 electronic table, accepted by ApJ
Submitted: 2015-07-17
The Kepler planet sample can only be used to reconstruct the underlying planet occurrence rate if the detection efficiency of the Kepler pipeline is known, here we present the results of a second experiment aimed at characterising this detection efficiency. We inject simulated transiting planet signals into the pixel data of ~10,000 targets, spanning one year of observations, and process the pixels as normal. We compare the set of detections made by the pipeline with the expectation from the set of simulated planets, and construct a sensitivity curve of signal recovery as a function of the signal-to-noise of the simulated transit signal train. The sensitivity curve does not meet the hypothetical maximum detection efficiency, however it is not as pessimistic as some of the published estimates of the detection efficiency. For the FGK stars in our sample, the sensitivity curve is well fit by a gamma function with the coefficients a = 4.35 and b = 1.05. We also find that the pipeline algorithms recover the depths and periods of the injected signals with very high fidelity, especially for periods longer than 10 days. We perform a simplified occurrence rate calculation using the measured detection efficiency compared to previous assumptions of the detection efficiency found in the literature to demonstrate the systematic error introduced into the resulting occurrence rates. The discrepancies in the calculated occurrence rates may go some way towards reconciling some of the inconsistencies found in the literature.
[35]  [pdf] - 1264048
Terrestrial Planet Occurrence Rates for the Kepler GK Dwarf Sample
Comments: 19 Pages, 17 Figures, Submitted ApJ. Python source to support Kepler pipeline completeness estimates available at
Submitted: 2015-06-12
We measure planet occurrence rates using the planet candidates discovered by the Q1-Q16 Kepler pipeline search. This study examines planet occurrence rates for the Kepler GK dwarf target sample for planet radii, 0.75<Rp<2.5 Rearth, and orbital periods, 50<Porb<300 days, with an emphasis on a thorough exploration and identification of the most important sources of systematic uncertainties. Integrating over this parameter space, we measure an occurrence rate of F=0.77 planets per star, with an allowed range of 0.3<F<1.9. The allowed range takes into account both statistical and systematic uncertainties, and values of F beyond the allowed range are significantly in disagreement with our analysis. We generally find higher planet occurrence rates and a steeper increase in planet occurrence rates towards small planets than previous studies of the Kepler GK dwarf sample. Through extrapolation, we find that the one year orbital period terrestrial planet occurrence rate, zeta_1=0.1, with an allowed range of 0.01<zeta_1<2, where zeta_1 is defined as the number of planets per star within 20% of the Rp and Porb of Earth. For G dwarf hosts, the zeta_1 parameter space is a subset of the larger eta_earth parameter space, thus zeta_1 places a lower limit on eta_earth for G dwarf hosts. From our analysis, we identify the leading sources of systematics impacting Kepler occurrence rate determinations as: reliability of the planet candidate sample, planet radii, pipeline completeness, and stellar parameters.
[36]  [pdf] - 1154942
The unification of powerful quasars and radio galaxies and their relation to other massive galaxies
Comments: 12 pages, 2 figures, 1 table. Accepted for publication in ApJ Letters
Submitted: 2015-05-13
The unification model for powerful radio galaxies and radio-loud quasars postulates that these objects are intrinsically the same but viewed along different angles. Herschel Space Observatory data permit the assessment of that model in the far-infrared spectral window. We analyze photometry from Spitzer and Herschel for the distant 3CR hosts, and find that radio galaxies and quasars have different mid-infrared, but indistinguishable far-infrared colors. Both these properties, the former being orientation dependent and the latter orientation invariant, are in line with expectations from the unification model. Adding powerful radio-quiet active galaxies and typical massive star-forming galaxies to the analysis, we demonstrate that infrared colors not only provide an orientation indicator, but can also distinguish active from star-forming galaxies.
[37]  [pdf] - 1182990
Circumstellar Disks revealed by $H$/$K$ Flux Variation Gradients
Comments: 5 pages, 5 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2015-04-28
The variability of young stellar objects (YSO) changes their brightness and color preventing a proper classification in traditional color-color and color magnitude diagrams. We have explored the feasibility of the flux variation gradient (FVG) method for YSOs, using $H$ and $K$ band monitoring data of the star forming region RCW\,38 obtained at the University Observatory Bochum in Chile. Simultaneous multi-epoch flux measurements follow a linear relation $F_{H}=\alpha + \beta \cdot F_{K}$ for almost all YSOs with large variability amplitude. The slope $\beta$ gives the mean $HK$ color temperature $T_{var}$ of the varying component. Because $T_{var}$ is hotter than the dust sublimation temperature, we have tentatively assigned it to stellar variations. If the gradient does not meet the origin of the flux-flux diagram, an additional non- or less-varying component may be required. If the variability amplitude is larger at the shorter wavelength, e.g. $\alpha < 0$, this component is cooler than the star (e.g. a circumstellar disk); vice versa, if $\alpha > 0$, the component is hotter like a scattering halo or even a companion star. We here present examples of two YSOs, where the $HK$ FVG implies the presence of a circumstellar disk; this finding is consistent with additional data at $J$ and $L$. One YSO shows a clear $K$-band excess in the $JHK$ color-color diagram, while the significance of a $K$-excess in the other YSO depends on the measurement epoch. Disentangling the contributions of star and disk it turns out that the two YSOs have huge variability amplitudes ($\sim 3-5$\,mag). The $HK$ FVG analysis is a powerful complementary tool to analyze the varying components of YSOs and worth further exploration of monitoring data at other wavelengths.
[38]  [pdf] - 1449779
The broad-line region and dust torus size of the Seyfert 1 galaxy PGC50427
Comments: 11 pages, 23 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2015-02-24
We present the results of a three years monitoring campaigns of the $z = 0.024$ type-1 active galactic nucleus (AGN) PGC50427. Through the use of Photometric Reverberation Mapping with broad and narrow band filters, we determine the size of the broad-line emitting region by measuring the time delay between the variability of the continuum and the H$\alpha$ emission line. The H$\alpha$ emission line responds to blue continuum variations with an average rest frame lag of $19.0 \pm 1.23$ days. Using single epoch spectroscopy we determined a broad-line H$\alpha$ velocity width of 1020 km s$^{-1}$ and in combination with the rest frame lag and adoption a geometric scaling factor $f = 5.5$, we calculate a black hole mass of $M_{BH} \sim 17 \times 10^{6} M_{\odot}$. Using the flux variation gradient method, we separate the host galaxy contribution from that of the AGN to calculate the rest frame 5100\AA~ luminosity at the time of our monitoring campaign. The rest frame lag and the host-subtracted luminosity permit us to derive the position of PGC50427 in the BLR size -- AGN luminosity diagram, which is remarkably close to the theoretically expected relation of $R \propto L^{0.5}$. The simultaneous optical and NIR ($J$ and $K_{s}$) observations allow us to determine the size of the dust torus through the use of dust reverberation mapping method. We find that the hot dust emission ($\sim 1800K$) lags the optical variations with an average rest frame lag of $46.2 \pm 2.60$ days. The dust reverberation radius and the nuclear NIR luminosity permit us to derive the position of PGC50427 on the known $\tau - M{V}$ diagram. The simultaneus observations for the broad-line region and dust thermal emission demonstrate that the innermost dust torus is located outside the BLR in PGC50427, supporting the unified scheme for AGNs. (Abstract shortened, see the manuscript.)
[39]  [pdf] - 1224352
Planetary Candidates Observed by Kepler VI: Planet Sample from Q1-Q16 (47 Months)
Comments: 18 pages, to be published in the Astrophysical Journal Supplement Series
Submitted: 2015-02-06
\We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive.
[40]  [pdf] - 1224258
The progenitor of the FUor-type young eruptive star 2MASS J06593158-0405277
Comments: 16 pages, 3 figures, 2 tables, accepted for publication in ApJ Letters
Submitted: 2015-01-30
Only a dozen confirmed FU Orionis-type young outbursting stars (FUors) are known today; this explains the interest in the recent FUor eruption of 2MASS J06593158-0405277. Its outburst and expected decline will be subject to numerous studies in the future. Almost equally important for the understanding of the eruption mechanism, however, is the physical characterization of the FUor's precursor. Here we analyze unpublished archival data and summarize - and partly revise - all relevant photometry from optical to submillimeter wavelengths. Our analysis implies that the FUor is possibly associated with eight T Tauri star candidates and a strong Class 0 source. Adopting a distance of 450 pc for the FUor, we derive a quiescent bolometric luminosity and temperature of L_bol = 4.8 L_Sun and T_bol = 1190 K, typical for young Class II sources. The central star has a temperature of T_eff = 4000 K, a mass of 0.75 M_Sun, and an age of about 6 x 10^5 yr. The SED implies a circumstellar mass of 0.01 - 0.06 M_Sun, and the system is surrounded by a faint infrared nebulosity. Our results provide an almost complete picture of a FUor progenitor, supporting the interpretation of future post-outburst studies.
[41]  [pdf] - 1224231
Planetary Candidates Observed by Kepler V: Planet Sample from Q1-Q12 (36 Months)
Comments: Accepted for publication, ApJS
Submitted: 2015-01-28, last modified: 2015-01-29
The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive.
[42]  [pdf] - 942390
Star formation in z>1 3CR host galaxies as seen by Herschel
Comments: 30 pages, 13 figures, 4 tables. Accepted for publication in A&A
Submitted: 2015-01-07
We present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z>1 3CR sources, from the Herschel GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG-strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation-based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation (negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether their supermassive black holes are accreting or not.
[43]  [pdf] - 1179586
Modelling photometric reverberation data -- a disk-like broad-line region and a potentially larger black hole mass for 3C120
Comments: 9 pages, 11 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2013-09-30, last modified: 2014-07-02
We consider photometric reverberation mapping, where the nuclear continuum variations are monitored via a broad-band filter and the echo of emission line clouds of the broad line region (BLR) is measured with a suitable narrow-band (NB) filter. We investigate how an incomplete emission-line coverage by the NB filter influences the BLR size determination. This includes two basic cases: 1) a symmetric cut of the blue and red part of the line wings, and 2) the filter positioned asymmetrically to the line centre so that essentially a complete half of the emission line is contained in the NB filter. Under the assumption that the BLR size is dominated by circular Keplerian orbits, we find that symmetric cutting of line wings may lead to overestimating the BLR size by less than 5%. The case of asymmetric half-line coverage, similar as for our data of the Seyfert 1 galaxy 3C120, yields the BLR size with a bias of less than 1%. Our results suggest that any BLR size bias due to narrow-band line cut in photometric reverberation mapping is small and in most cases negligible. We used well sampled photometric reverberation mapping light curves with sharp variation features in both the continuum and the Hbeta light curves to determine the geometry type of the Hbeta BLR for 3C120. Modelling of the light curve, under the assumption that the BLR is essentially virialised, argues against a spherical geometry and favours a nearly face-on disk-like geometry with inclination i = 10 +/- 4 deg and extension from 22 to 28 light days. The low inclination may lead to a larger black hole mass than the derived when using the average geometry scaling factor f=5.5. We discuss deviations of Seyfert 1 galaxies from the M_BH - sigma relation.
[44]  [pdf] - 1210100
3C 220.3: a radio galaxy lensing a submillimeter galaxy
Comments: 17 pages, 14 Figures, accepted for publication in ApJ
Submitted: 2014-06-11
Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.
[45]  [pdf] - 1209552
Precision asteroseismology of the pulsating white dwarf GD 1212 using a two-wheel-controlled Kepler spacecraft
Comments: 8 pages, 4 figures, accepted for publication in The Astrophysical Journal
Submitted: 2014-05-14
We present a preliminary analysis of the cool pulsating white dwarf GD 1212, enabled by more than 11.5 days of space-based photometry obtained during an engineering test of the two-reaction-wheel-controlled Kepler spacecraft. We detect at least 19 independent pulsation modes, ranging from 828.2-1220.8 s, and at least 17 nonlinear combination frequencies of those independent pulsations. Our longest uninterrupted light curve, 9.0 days in length, evidences coherent difference frequencies at periods inaccessible from the ground, up to 14.5 hr, the longest-period signals ever detected in a pulsating white dwarf. These results mark some of the first science to come from a two-wheel-controlled Kepler spacecraft, proving the capability for unprecedented discoveries afforded by extending Kepler observations to the ecliptic.
[46]  [pdf] - 844879
Broad-line region structure and kinematics in the radio galaxy 3C 120
Comments: 18 pages, 25 figures, Astronomy & Astrophysics in press
Submitted: 2014-05-07
Broad emission lines originate in the surroundings of supermassive black holes in the centers of active galactic nuclei (AGN). One method to investigate the extent, structure, and kinematics of the BLR is to study the continuum and line profile variability in AGN. We selected the radio-loud Seyfert 1 galaxy 3C 120 as a target for this study. We took spectra with a high signal-to-noise ratio of 3C 120 with the 9.2m Hobby-Eberly Telescope between Sept. 2008 and March 2009. In parallel, we photometrically monitored the continuum flux at the Wise observatory. We analyzed the continuum and line profile variations in detail (1D and 2D reverberation mapping) and modeled the geometry of the line-emitting regions based on the line profiles. We show that the BLR in 3C 120 is stratified with respect to the distance of the line-emitting regions from the center with respect to the line widths (FWHM) of the rms profiles and with respect to the variability amplitude of the emission lines. The emission line wings of H{\alpha} and H{\beta} respond much faster than their central region. This is explained by accretion disk models. In addition, these lines show a stronger response in the red wings. However, the velocity-delay maps of the helium lines show a stronger response in the blue wing. Furthermore, the HeII{\lambda}4686 line responds faster in the blue wing in contradiction to observations made one and a half years later when the galaxy was in a lower state. The faster response in the blue wing is an indication for central outflow motions when this galaxy was in a bright state during our observations. The vertical BLR structure in 3C 120 coincides with that of other AGN. We confirm the general trend: the emission lines of narrow line AGN originate at larger distances from the midplane than AGN with broader emission lines.
[47]  [pdf] - 1173339
Damped Ly{\alpha} Absorption Systems in Semi-Analytic Models with Multiphase Gas
Comments: 27 pages, 15 figures, submitted to MNRAS
Submitted: 2013-08-12, last modified: 2014-04-23
We investigate the properties of damped Ly{\alpha} absorption systems (DLAs) in semi-analytic models of galaxy formation, including partitioning of cold gas in galactic discs into atomic, molecular, and ionized phases with a molecular gas-based star formation recipe. We investigate two approaches for partitioning gas into these constituents: a pressure-based and a metallicity-based recipe. We identify DLAs by passing lines of sight through our simulations to compute HI column densities. We find that models with "standard" gas radial profiles - where the average specific angular momentum of the gas disc is equal to that of the host dark matter halo - fail to reproduce the observed column density distribution of DLAs. These models also fail to reproduce the distribution of velocity widths {\Delta}v, overproducing low {\Delta}v relative to high {\Delta}v systems. Models with "extended" radial gas profiles - corresponding to gas discs with higher specific angular momentum - are able to reproduce quite well the column density distribution of absorbers over the column density range 19 < log NHI < 22.5 in the redshift range 2 < z < 3.5. The model with pressure-based gas partitioning also reproduces the observed line density of DLAs, HI gas density, and {\Delta}v distribution at z < 3 remarkably well. However all of the models investigated here underproduce DLAs and the HI gas density at z > 3. If this is the case, the flatness in the number of DLAs and HI gas density over the redshift interval 0 < z < 5 may be due to a cosmic coincidence where the majority of DLAs at z > 3 arise from intergalactic gas in filaments while those at z < 3 arise predominantly in galactic discs. We further investigate the dependence of DLA metallicity on redshift and {\Delta}v, and find reasonably good agreement with the observations, particularly when including the effects of metallicity gradients (abbrv.).
[48]  [pdf] - 1202585
Contamination in the Kepler Field. Identification of 685 KOIs as False Positives Via Ephemeris Matching Based On Q1-Q12 Data
Comments: Accepted to the Astronomical Journal on February 21, 2014. 13 pages, 7 figures, 2 tables. The tables are shortened in the text, but the full .tex file for each table is available as arXiv ancillary files. (See "Ancillary File" links below "Download" on the right side of the page.)
Submitted: 2014-01-06, last modified: 2014-04-14
The Kepler mission has to date found almost 6,000 planetary transit-like signals, utilizing three years of data for over 170,000 stars at extremely high photometric precision. Due to its design, contamination from eclipsing binaries, variable stars, and other transiting planets results in a significant number of these signals being false positives. This directly affects the determination of the occurrence rate of Earth-like planets in our Galaxy, as well as other planet population statistics. In order to detect as many of these false positives as possible, we perform ephemeris matching among all transiting planet, eclipsing binary, and variable star sources. We find that 685 Kepler Objects of Interest - 12% of all those analyzed - are false positives as a result of contamination, due to 409 unique parent sources. Of these, 118 have not previously been identified by other methods. We estimate that ~35% of KOIs are false positives due to contamination, when performing a first-order correction for observational bias. Comparing single-planet candidate KOIs to multi-planet candidate KOIs, we find an observed false positive fraction due to contamination of 16% and 2.4% respectively, bolstering the existing evidence that multi-planet KOIs are significantly less likely to be false positives. We also analyze the parameter distributions of the ephemeris matches and derive a simple model for the most common type of contamination in the Kepler field. We find that the ephemeris matching technique is able to identify low signal-to-noise false positives that are difficult to identify with other vetting techniques. We expect false positive KOIs to become more frequent when analyzing more quarters of Kepler data, and note that many of them will not be able to be identified based on Kepler data alone.
[49]  [pdf] - 1208796
Rapidly growing black holes and host galaxies in the distant Universe from the Herschel Radio Galaxy Evolution Project
Comments: 37 pages, accepted to A&A
Submitted: 2014-04-03
We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.
[50]  [pdf] - 1208751
Stellar Signatures of AGN Jet Triggered Star Formation
Submitted: 2014-04-01
To investigate feedback between relativistic jets emanating from Active Galactic Nuclei (AGN) and the stellar population of the host galaxy, we analyze the long-term evolution of the galaxy-scale simulations by Gaibler et al. (2012) of jets in massive, gas-rich galaxies at z ~ 2 - 3 and of stars formed in the host galaxies. We find strong, jet-induced differences in the resulting stellar populations of galaxies that host relativistic jets and galaxies that do not, including correlations in stellar locations, velocities, and ages. Jets are found to generate distributions of increased radial and vertical velocities that persist long enough to effectively extend the stellar structure of the host. The jets cause the formation of bow shocks that move out through the disk, generating rings of star formation within the disk. The bow shock often accelerates pockets of gas in which stars form, yielding populations of stars with significant radial and vertical velocities, some of which have large enough velocities to escape the galaxy. These stellar population signatures can serve to identify past jet activity as well as jet-induced star formation.
[51]  [pdf] - 1203483
The K2 Mission: Characterization and Early results
Comments: 25 pages, 11 figures, Accepted to PASP
Submitted: 2014-02-20, last modified: 2014-03-07
The K2 mission will make use of the Kepler spacecraft and its assets to expand upon Kepler's groundbreaking discoveries in the fields of exoplanets and astrophysics through new and exciting observations. K2 will use an innovative way of operating the spacecraft to observe target fields along the ecliptic for the next 2-3 years. Early science commissioning observations have shown an estimated photometric precision near 400 ppm in a single 30 minute observation, and a 6-hour photometric precision of 80 ppm (both at V=12). The K2 mission offers long-term, simultaneous optical observation of thousands of objects at a precision far better than is achievable from ground-based telescopes. Ecliptic fields will be observed for approximately 75-days enabling a unique exoplanet survey which fills the gaps in duration and sensitivity between the Kepler and TESS missions, and offers pre-launch exoplanet target identification for JWST transit spectroscopy. Astrophysics observations with K2 will include studies of young open clusters, bright stars, galaxies, supernovae, and asteroseismology.
[52]  [pdf] - 1203585
Validation of Kepler's Multiple Planet Candidates. III: Light Curve Analysis & Announcement of Hundreds of New Multi-planet Systems
Comments: 138 pages, 8 Figures, 5 Tables. Accepted for publications in the Astrophysical Journal
Submitted: 2014-02-26
The Kepler mission has discovered over 2500 exoplanet candidates in the first two years of spacecraft data, with approximately 40% of them in candidate multi-planet systems. The high rate of multiplicity combined with the low rate of identified false-positives indicates that the multiplanet systems contain very few false-positive signals due to other systems not gravitationally bound to the target star (Lissauer, J. J., et al., 2012, ApJ 750, 131). False positives in the multi- planet systems are identified and removed, leaving behind a residual population of candidate multi-planet transiting systems expected to have a false-positive rate less than 1%. We present a sample of 340 planetary systems that contain 851 planets that are validated to substantially better than the 99% confidence level; the vast majority of these have not been previously verified as planets. We expect ~2 unidentified false-positives making our sample of planet very reliable. We present fundamental planetary properties of our sample based on a comprehensive analysis of Kepler light curves and ground-based spectroscopy and high-resolution imaging. Since we do not require spectroscopy or high-resolution imaging for validation, some of our derived parameters for a planetary system may be systematically incorrect due to dilution from light due to additional stars in the photometric aperture. None the less, our result nearly doubles the number of verified exoplanets.
[53]  [pdf] - 1203535
Spectral Energy Distributions of QSOs at z>5: common AGN-heated dust and occasionally strong star-formation
Comments: 23 pages, 15 figures, 7 tables; accepted for publication in ApJ; a version of the paper including all online-only panels of Figure A.1 can be found at
Submitted: 2014-02-24
We present spectral energy distributions (SEDs) of 69 QSOs at z>5, covering a rest frame wavelength range of 0.1mu to ~80mu, and centered on new Spitzer and Herschel observations. The detection rate of the QSOs with Spitzer is very high (97% at lambda_rest ~< 4mu), but drops towards the Herschel bands with 30% detected in PACS (rest frame mid-infrared) and 15% additionally in the SPIRE (rest frame far-infrared; FIR). We perform multi-component SED fits for Herschel-detected objects and confirm that to match the observed SEDs, a clumpy torus model needs to be complemented by a hot (~1300K) component and, in cases with prominent FIR emission, also by a cold (~50K) component. In the FIR detected cases the luminosity of the cold component is on the order of 10^13 L_sun which is likely heated by star formation. From the SED fits we also determine that the AGN dust-to-accretion disk luminosity ratio declines with UV/optical luminosity. Emission from hot (~1300K) dust is common in our sample, showing that nuclear dust is ubiquitous in luminous QSOs out to redshift 6. However, about 15% of the objects appear under-luminous in the near infrared compared to their optical emission and seem to be deficient in (but not devoid of) hot dust. Within our full sample, the QSOs detected with Herschel are found at the high luminosity end in L_UV/opt and L_NIR and show low equivalent widths (EWs) in H_alpha and in Ly_alpha. In the distribution of H_alpha EWs, as determined from the Spitzer photometry, the high-redshift QSOs show little difference to low redshift AGN.
[54]  [pdf] - 772442
Masses, Radii, and Orbits of Small Kepler Planets: The Transition from Gaseous to Rocky Planets
Marcy, Geoffrey W.; Isaacson, Howard; Howard, Andrew W.; Rowe, Jason F.; Jenkins, Jon M.; Bryson, Stephen T.; Latham, David W.; Howell, Steve B.; Gautier, Thomas N.; Batalha, Natalie M.; Rogers, Leslie A.; Ciardi, David; Fischer, Debra A.; Gilliland, Ronald L.; Kjeldsen, Hans; Christensen-Dalsgaard, Jørgen; Huber, Daniel; Chaplin, William J.; Basu, Sarbani; Buchhave, Lars A.; Quinn, Samuel N.; Borucki, William J.; Koch, David G.; Hunter, Roger; Caldwell, Douglas A.; Van Cleve, Jeffrey; Kolbl, Rea; Weiss, Lauren M.; Petigura, Erik; Seager, Sara; Morton, Timothy; Johnson, John Asher; Ballard, Sarah; Burke, Chris; Cochran, William D.; Endl, Michael; MacQueen, Phillip; Everett, Mark E.; Lissauer, Jack J.; Ford, Eric B.; Torres, Guillermo; Fressin, Francois; Brown, Timothy M.; Steffen, Jason H.; Charbonneau, David; Basri, Gibor S.; Sasselov, Dimitar D.; Winn, Joshua; Sanchis-Ojeda, Roberto; Christiansen, Jessie; Adams, Elisabeth; Henze, Christopher; Dupree, Andrea; Fabrycky, Daniel C.; Fortney, Jonathan J.; Tarter, Jill; Holman, Matthew J.; Tenenbaum, Peter; Shporer, Avi; Lucas, Philip W.; Welsh, William F.; Orosz, Jerome A.; Bedding, T. R.; Campante, T. L.; Davies, G. R.; Elsworth, Y.; Handberg, R.; Hekker, S.; Karoff, C.; Kawaler, S. D.; Lund, M. N.; Lundkvist, M.; Metcalfe, T. S.; Miglio, A.; Aguirre, V. Silva; Stello, D.; White, T. R.; Boss, Alan; Devore, Edna; Gould, Alan; Prsa, Andrej; Agol, Eric; Barclay, Thomas; Coughlin, Jeff; Brugamyer, Erik; Mullally, Fergal; Quintana, Elisa V.; Still, Martin; hompson, Susan E.; Morrison, David; Twicken, Joseph D.; Désert, Jean-Michel; Carter, Josh; Crepp, Justin R.; Hébrard, Guillaume; Santerne, Alexandre; Moutou, Claire; Sobeck, Charlie; Hudgins, Douglas; Haas, Michael R.; Robertson, Paul; Lillo-Box, Jorge; Barrado, David
Comments: 94 pages, 55 figures, 25 tables. Accepted by ApJS
Submitted: 2014-01-16
We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).
[55]  [pdf] - 769891
Dust reverberation-mapping of the Seyfert 1 galaxy WPVS48
Comments: 4 pages, 6 figures, Published in Astronomy and Astrophysics
Submitted: 2014-01-13
Using robotic telescopes of the Universitatssternwarte Bochum near Cerro Armazones in Chile, we monitored the z=0.0377 Seyfert 1 galaxy WPVS48 (2MASX J09594263-3112581) in the optical (B and R) and near-infrared (NIR, J and Ks) with a cadence of two days. The light curves show unprecedented variability details. The NIR variation features of WPVS48 are consistent with the corresponding optical variations, but the features appear sharper in the NIR than in the optical, suggesting that the optical photons undergo multiple scatterings. The J and Ks emission, tracing the hot (1600 K) dust echo, lags the B and R variations by on average 64 +/- 4 days and 71 +/- 5 days, respectively (restframe). WPVS48 lies on the known tau-M_V relationship. However, the observed lag is about three times shorter than expected from the dust sublimation radius r_sub inferred from the optical-UV luminosity, and explanations for this common discrepancy are searched for. The sharp NIR echos argue for a face-on torus geometry and allow us to put forward two potential scenarios: 1) as previously proposed, in the equatorial plane of the accretion disk the inner region of the torus is flattened and may come closer to the accretion disk. 2) The dust torus with inner radius r_sub is geometrically and optically thick, so that the observer only sees the facing rim of the torus wall, which lies closer to the observer than the torus equatorial plane and therefore leads to an observed foreshortened lag. Both scenarios are able to explain the factor three discrepancy between tau and r_sub. Longer-wavelength dust reverberation data might enable one to distinguish between the scenarios.
[56]  [pdf] - 1180422
Detection of Potential Transit Signals in Sixteen Quarters of Kepler Mission Data
Comments: Accepted by ApJ Supplement
Submitted: 2013-11-01, last modified: 2014-01-07
We present the results of a search for potential transit signals in four years of photometry data acquired by the Kepler Mission. The targets of the search include 111,800 stars which were observed for the entire interval and 85,522 stars which were observed for a subset of the interval. We found that 9,743 targets contained at least one signal consistent with the signature of a transiting or eclipsing object, where the criteria for detection are periodicity of the detected transits, adequate signal-to-noise ratio, and acceptance by a number of tests which reject false positive detections. When targets that had produced a signal were searched repeatedly, an additional 6,542 signals were detected on 3,223 target stars, for a total of 16,285 potential detections. Comparison of the set of detected signals with a set of known and vetted transit events in the Kepler field of view shows that the recovery rate for these signals is 96.9%. The ensemble properties of the detected signals are reviewed.
[57]  [pdf] - 1202361
Planetary Candidates Observed by Kepler IV: Planet Sample From Q1-Q8 (22 Months)
Comments: 12 pages, 8 figures, Accepted ApJ Supplement
Submitted: 2013-12-18
We provide updates to the Kepler planet candidate sample based upon nearly two years of high-precision photometry (i.e., Q1-Q8). From an initial list of nearly 13,400 Threshold Crossing Events (TCEs), 480 new host stars are identified from their flux time series as consistent with hosting transiting planets. Potential transit signals are subjected to further analysis using the pixel-level data, which allows background eclipsing binaries to be identified through small image position shifts during transit. We also re-evaluate Kepler Objects of Interest (KOI) 1-1609, which were identified early in the mission, using substantially more data to test for background false positives and to find additional multiple systems. Combining the new and previous KOI samples, we provide updated parameters for 2,738 Kepler planet candidates distributed across 2,017 host stars. From the combined Kepler planet candidates, 472 are new from the Q1-Q8 data examined in this study. The new Kepler planet candidates represent ~40% of the sample with Rp~1 Rearth and represent ~40% of the low equilibrium temperature (Teq<300 K) sample. We review the known biases in the current sample of Kepler planet candidates relevant to evaluating planet population statistics with the current Kepler planet candidate sample.
[58]  [pdf] - 1157757
Physical properties of simulated galaxy populations at z=2 - II. Effects of cosmology, reionization and ISM physics
Comments: 11 pages, 6 figures, accepted by MNRAS
Submitted: 2012-11-13, last modified: 2013-08-06
We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on the assumed star formation law, the equation of state imposed on the unresolved interstellar medium, the stellar initial mass function, the reionization history, and the assumed cosmology. This work complements that of Paper I, where we studied the effects of varying models for galactic winds driven by star formation and AGN. The normalisation of the matter power spectrum strongly affects the galaxy mass function, but has a relatively small effect on the physical properties of galaxies residing in haloes of a fixed mass. Reionization suppresses the stellar masses and gas fractions of low-mass galaxies, but by z = 2 the results are insensitive to the timing of reionization. The stellar initial mass function mainly determines the physical properties of galaxies through its effect on the efficiency of the feedback, while changes in the recycled mass and metal fractions play a smaller role. If we use a recipe for star formation that reproduces the observed star formation law independently of the assumed equation of state of the unresolved ISM, then the latter is unimportant. The star formation law, i.e. the gas consumption time scale as a function of surface density, determines the mass of dense, star-forming gas in galaxies, but affects neither the star formation rate nor the stellar mass. This can be understood in terms of self-regulation: the gas fraction adjusts until the outflow rate balances the inflow rate.
[59]  [pdf] - 1157527
Physical properties of simulated galaxy populations at z=2 - I. Effect of metal-line cooling and feedback from star formation and AGN
Comments: 25 pages, 15 figures, accepted by MNRAS
Submitted: 2012-11-05, last modified: 2013-08-06
We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range.
[60]  [pdf] - 694662
The substellar content of the ONC
Comments: 3 pages, 3 figures, Poster summary to appear in the proceedings of the conference 'Brown dwarfs come of age', May 20-24 2013, Memorie della Societa Astronomica Italiana
Submitted: 2013-07-17
The Substellar Initial Mass Function (SIMF) of many star-forming regions is still poorly known but the detailed knowledge of its shape will help to distinguish among the substellar formation theories. The Orion Nebula Cluster (ONC) is one of the most extensively studied star forming regions. We here present deep, wide-field JHK observations of the ONC taken with HAWK-I@VLT. These observations extend the IMF into the brown dwarf and free-floating planetary mass regime with unprecedented sensitivity. To obtain a clean sample of ONC members, we exclude potential background sources with the help of CO extinction maps. Masses are assigned by means of evolutionary tracks in the H vs. J-H Color-Magnitude Diagram (CMD). Besides the well known stellar peak at ~0.25 M_sun we find a pronounced second peak at ~0.04 M_sun in the SIMF and indications for a third rise in the free-floating planetary mass regime.
[61]  [pdf] - 1172232
Eclipsing high-mass binaries I. Light curves and system parameters for CPD-518946, PISMIS24-1 and HD319702
Comments: 7 pages, 4 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2013-06-23
We present first results of a comprehensive photometric O-star survey performed with a robotic twin refractor at the Universit\"atssternwarte Bochum located near Cerro Armazones in Chile. For three high-mass stars, namely Pismis24-1, CPD-518946 and HD319702, we determined the period through the Lafler-Kinman algorithm and model the light curves within the framework of the Roche geometry. For Pismis24-1, a previously known eclipsing binary, we provide first light curves and determined a photometric period of 2.36 days together with an orbital inclination of 61.8 degrees. The best-fitting model solution to the light curves suggest a detached configuration. With a primary temperature of T1 = 42520K we obtain the temperature of the secondary component as T2 = 41500K. CPD-518946 is another known eclipsing binary for which we present a revised photometric period of 1.96 days with an orbital inclination of 58.4 degrees. The system has likely a semi-detached configuration and a mass ratio q = M1/M2 = 2.8. If we adopt a primary temperature of T1 = 34550K we obtain T2 = 21500K for the secondary component. HD319702 is a newly discovered eclipsing binary member of the young open cluster NGC6334. The system shows well-defined eclipses favouring a detached configuration with a period of 2.0 days and an orbital inclination of 67.5 degrees. Combining our photometric result with the primary spectral type O8 III(f) (T1 = 34000K) we derive a temperature of T2 = 25200K for the secondary component.
[62]  [pdf] - 1164940
Identification of Background False Positives from Kepler Data
Comments: 72 pages, 38 figures. To appear in PASP
Submitted: 2013-02-28, last modified: 2013-06-11
The Kepler Mission was launched on March 6, 2009 to perform a photometric survey of more than 100,000 dwarf stars to search for Earth-size planets with the transit technique. The reliability of the resulting planetary candidate list relies on the ability to identify and remove false positives. Major sources of astrophysical false positives are planetary transits and stellar eclipses on background stars. We describe several new techniques for the identification of background transit sources that are separated from their target stars, indicating an astrophysical false positive. These techniques use only Kepler photometric data. We describe the concepts and construction of these techniques in detail as well as their performance and relative merits.
[63]  [pdf] - 1171791
Revealing the heavily obscured AGN population of High Redshift 3CRR Sources with Chandra X-ray Observations
Comments: 42 pages, 11 figures, 3 tables; ApJ, in press
Submitted: 2013-06-04
Chandra observations of a complete, flux-limited sample of 38 high-redshift (1<z<2), low-frequency selected (and so unbiased in orientation) 3CRR radio sources are reported. The sample includes 21 quasars (= broad line radio galaxies) and 17 narrow line radio galaxies (NLRGs) with matched 178 MHz radio luminosity (log L_R ~ 44-45). The quasars have high radio core-fraction, high X-ray luminosities (log L_X ~ 45-46) and soft X-ray hardness ratios (HR ~ -0.5) indicating low obscuration. The NLRGs have lower core-fraction, lower apparent X-ray luminosities (log L_X ~ 43-45) and mostly hard X-ray hardness ratios (HR>0) indicating obscuration (log N_H ~ 22-24 cm^-2). These properties and the correlation between obscuration and radio core-fraction are consistent with orientation-dependent obscuration as in Unification models. About half the NLRGs have soft X-ray hardness ratios and/or high [OIII] emission line to X-ray luminosity ratio suggesting obscuration by Compton thick (CT) material so that scattered nuclear or extended X-ray emission dominates (as in NGC1068). The ratios of unobscured to Compton-thin (10^{22}< N_H(int) < 1.5 x 10^{24} cm^-2) to CT (N_H(int) > 1.5 x 10^{24} cm^-2) is 2.5:1.4:1 in this high luminosity, radio-selected sample. The obscured fraction is 0.5, higher than is typically reported for AGN at comparable luminosities from multi-wavelength surveys (0.1-0.3). Assuming random nuclear orientation, the unobscured half-opening angle of the disk/wind/torus structure is ~ 60deg and the obscuring material covers 30deg of which ~ 12deg is Compton thick. The multi-wavelength properties reveal that many NLRGs have intrinsic absorption 10-1000x higher than indicated by their X-ray hardness ratios, and their true L_X values are ~10--100x larger than the hardness-ratio absorption corrections would indicate.
[64]  [pdf] - 667654
A sub-Mercury-sized exoplanet
Comments: Accepted and published in Nature (2013 Feb 28). This is the submitted version of paper, merged with the Supplementary Information
Submitted: 2013-05-23
Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.
[65]  [pdf] - 1171434
Complete infrared spectral energy distributions of mm detected quasars at z>5
Comments: 15 pages, 7 figures, accepted for publication in ApJ
Submitted: 2013-05-17
We present Herschel far-infrared (FIR) photometry of eleven quasars at redshift z>5 that have previously been detected at 1.2mm. We perform full spectral energy distribution (SED) fits over the wavelength range lambda_rest ~0.1-400mu for those objects with good Herschel detections. These fits reveal the need for an additional far-infrared (FIR) component besides the emission from a dusty AGN-powered torus. This additional FIR component has temperatures of T_FIR ~ 40-60K with luminosities of L_(8-1000mu) ~ 10^13 L_sun (accounting for 25-60% of the bolometric FIR luminosity). If the FIR dust emission is due to star formation it would suggest star formation rates in excess of 1000 solar masses per year. We show that at long wavelengths (lambda_rest > 50mu) the contribution of the AGN-powered torus emission is negligible. This explains how previous FIR studies of high-redshift quasars that relied on single component fits to (ground-based) observations at lambda_obs > 350mu reached T_FIR and L_FIR values similar to our complete SED fits. Stacking the Herschel data of four individually undetected sources reveals a significant average signal in the PACS bands but not in SPIRE. The average SED of sources with individual Herschel detections shows a striking surplus in near- and mid-infrared emission when compared to common AGN templates. The comparison between two average SEDs (sources with and without individual Herschel detections) matched in the UV/optical indicates that for these objects the strength of the MIR emission may correlate with the strength of the FIR emission.
[66]  [pdf] - 657404
Kepler-62: A five-planet system with planets of 1.4 and 1.6 Earth radii in the Habitable Zone
Comments: Published online 18 April 2013 in Science Express
Submitted: 2013-04-27
We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.
[67]  [pdf] - 1166054
A super-Earth-sized planet orbiting in or near the habitable zone around Sun-like star
Comments: Accepted for publication in the Astrophysical Journal
Submitted: 2013-04-17
We present the discovery of a super-earth-sized planet in or near the habitable zone of a sun-like star. The host is Kepler-69, a 13.7 mag G4V-type star. We detect two periodic sets of transit signals in the three-year flux time series of Kepler-69, obtained with the Kepler spacecraft. Using the very high precision Kepler photometry, and follow-up observations, our confidence that these signals represent planetary transits is >99.1%. The inner planet, Kepler-69b, has a radius of 2.24+/-0.4 Rearth and orbits the host star every 13.7 days. The outer planet, Kepler-69c, is a super-Earth-size object with a radius of 1.7+/-0.3 Rearth and an orbital period of 242.5 days. Assuming an Earth-like Bond albedo, Kepler-69c has an equilibrium temperature of 299 +/- 19 K, which places the planet close to the habitable zone around the host star. This is the smallest planet found by Kepler to be orbiting in or near habitable zone of a Sun-like star and represents an important step on the path to finding the first true Earth analog.
[68]  [pdf] - 650488
Detection of Potential Transit Signals in the First Twelve Quarters of Kepler Mission Data
Comments: 29 pages, 14 figures, accepted for publication in ApJ Supp
Submitted: 2012-12-12, last modified: 2013-04-09
We present the results of a search for potential transit signals in the first three years of photometry data acquired by the Kepler Mission. The targets of the search include 112,321 targets which were observed over the full interval and an additional 79,992 targets which were observed for a subset of the full interval. From this set of targets we find a total of 11,087 targets which contain at least one signal which meets the Kepler detection criteria: those criteria are periodicity of the signal, an acceptable signal-to-noise ratio, and three tests which reject false positives. Each target containing at least one detected signal is then searched repeatedly for additional signals, which represent multi-planet systems of transiting planets. When targets with multiple detections are considered, a total of 18,406 potential transiting planet signals are found in the Kepler Mission dataset. The detected signals are dominated by events with relatively low signal-to-noise ratios and by events with relatively short periods. The distribution of estimated transit depths appears to peak in the range between 20 and 30 parts per million, with a few detections down to fewer than 10 parts per million. The detections exhibit signal-to-noise ratios from 7.1 sigma, which is the lower cut-off for detections, to over 10,000 sigma, and periods ranging from 0.5 days, which is the shortest period searched, to 525 days, which is the upper limit of achievable periods given the length of the data set and the requirement that all detections include at least 3 transits. The detected signals are compared to a set of known transit events in the Kepler field of view, many of which were identified by alternative methods; the comparison shows that the current search recovery rate for targets with known transit events is 98.3%.
[69]  [pdf] - 1164544
Fundamental Properties of Kepler Planet-Candidate Host Stars using Asteroseismology
Comments: 19 pages, 10 figures, 4 tables; accepted for publication in ApJ; machine-readable versions of tables 1-3 are available as ancillary files or in the source code; v2: minor changes to match published version
Submitted: 2013-02-11, last modified: 2013-03-26
We have used asteroseismology to determine fundamental properties for 66 Kepler planet-candidate host stars, with typical uncertainties of 3% and 7% in radius and mass, respectively. The results include new asteroseismic solutions for four host stars with confirmed planets (Kepler-4, Kepler-14, Kepler-23 and Kepler-25) and increase the total number of Kepler host stars with asteroseismic solutions to 77. A comparison with stellar properties in the planet-candidate catalog by Batalha et al. shows that radii for subgiants and giants obtained from spectroscopic follow-up are systematically too low by up to a factor of 1.5, while the properties for unevolved stars are in good agreement. We furthermore apply asteroseismology to confirm that a large majority of cool main-sequence hosts are indeed dwarfs and not misclassified giants. Using the revised stellar properties, we recalculate the radii for 107 planet candidates in our sample, and comment on candidates for which the radii change from a previously giant-planet/brown-dwarf/stellar regime to a sub-Jupiter size, or vice versa. A comparison of stellar densities from asteroseismology with densities derived from transit models in Batalha et al. assuming circular orbits shows significant disagreement for more than half of the sample due to systematics in the modeled impact parameters, or due to planet candidates which may be in eccentric orbits. Finally, we investigate tentative correlations between host-star masses and planet candidate radii, orbital periods, and multiplicity, but caution that these results may be influenced by the small sample size and detection biases.
[70]  [pdf] - 1165512
The extended narrow-line region of two type-I quasi-stellar objects
Comments: 7 pages, 5 figures, accepted for publication in ApJ
Submitted: 2013-03-25
We investigate the narrow-line region (NLR) of two radio-quiet QSOs, PG1012+008 and PG1307+085, using high signal-to-noise spatially resolved long-slit spectra obtained with FORS1 at the Very Large Telescope. Although the emission is dominated by the point-spread function of the nuclear source, we are able to detect extended NLR emission out to several kpc scales in both QSOs by subtracting the scaled central spectrum from outer spectra. In contrast to the nuclear spectrum, which shows a prominent blue wing and a broad line profile of the [O III] line, the extended emission reveals no clear signs of large scale outflows. Exploiting the wide wavelength range, we determine the radial change of the gas properties in the NLR, i.e., gas temperature, density, and ionization parameter, and compare them with those of Seyfert galaxies and type-II QSOs. The QSOs have higher nuclear temperature and lower electron density than Seyferts, but show no significant difference compared to type-II QSOs, while the ionization parameter decreases with radial distance, similar to the case of Seyfert galaxies. For PG1012+008, we determine the stellar velocity dispersion of the host galaxy. Combined with the black hole mass, we find that the luminous radio-quiet QSO follows the local M_BH-sigma* relation of active galactic nuclei.
[71]  [pdf] - 639919
Photometric reverberation mapping of 3C120
Comments: 11 pages, 11 figures, Published in Astronomy and Astrophysics
Submitted: 2013-03-14
We present the results of a five month monitoring campaign of the local active galactic nuclei (AGN) 3C120. Observations with a median sampling of two days were conducted with the robotic 15cm telescope VYSOS-6 located near Cerro Armazones in Chile. Broad band (B,V) and narrow band (NB) filters were used in order to measure fluxes of the AGN and the H_beta broad line region (BLR) emission line. The NB flux is constituted by about 50% continuum and 50% H_beta emission line. To disentangle line and continuum flux, a synthetic H_beta light curve was created by subtracting a scaled V-band light curve from the NB light curve. Here we show that the H_beta emission line responds to continuum variations with a rest frame lag of 23.6 +/- 1.69 days. We estimate a virial mass of the central black hole M_BH = 57 +/- 27 * 10^6 solar masses, by combining the obtained lag with the velocity dispersion of a single contemporaneous spectrum. Using the flux variation gradient (FVG) method, we determined the host galaxy subtracted rest frame 5100A luminosity at the time of our monitoring campaign with an uncertainty of 10% (L_AGN = 6.94 +/- 0.71* 10^43 ergs^-1). Compared with recent spectroscopic reverberation results, 3C120 shifts in the R_BLR - L_AGN diagram remarkably close to the theoretically expected relation of R-L^0.5. Our results demonstrate the performance of photometric AGN reverberation mapping, in particular for efficiently determining the BLR size and the AGN luminosity
[72]  [pdf] - 639948
Size and disk-like shape of the broad-line region of ESO399-IG20
Comments: 9 pages, 13 figures, Published in Astronomy and Astrophysics
Submitted: 2013-03-14
We present photometric reverberation mapping of the narrow-line Seyfert 1 galaxy ESO399-IG20 performed with the robotic 15 cm telescope VYSOS-6 at the Cerro Armazones Observatory. Through the combination of broad- and narrow-band filters we determine the size of the broad-line emitting region (BLR) by measuring the time delay between the variability of the continuum and the H_alpha emission line. We use the flux variation gradient method to separate the host galaxy contribution from that of the active galactic nucleus (AGN), and to calculate the 5100A luminosity L_AGN of the AGN. Both measurements permit us to derive the position of ESO399-IG20 in the BLR size -- AGN luminosity R_BLR-L_AGN^0.5 diagram. We infer the basic geometry of the BLR through modelling of the light curves. The pronounced sharp variability patterns in both the continuum and the emission line light curves allow us to reject a spherical BLR geometry. The light curves are best fitted by a disk-like BLR seen nearly face-on with an inclination angle of 6+/-3 degrees and with an extension from 16 to 20 light days.
[73]  [pdf] - 1164952
Measuring Transit Signal Recovery in the Kepler Pipeline I: Individual Events
Comments: 40 pages, 18 figures, submitted to ApJ
Submitted: 2013-03-01
The Kepler Mission was designed to measure the frequency of Earth-size planets in the habitable zone of Sun-like stars. A crucial component for recovering the underlying planet population from a sample of detected planets is understanding the completeness of that sample - what fraction of the planets that could have been discovered in a given data set were actually detected. Here we outline the information required to determine the sample completeness, and describe an experiment to address a specific aspect of that question, which is the issue of transit signal recovery. We investigate the extent to which the Kepler pipeline preserves individual transit signals by injecting simulated transits into the pixel-level data, processing the modified pixels through the pipeline, and comparing the measured transit signal-to-noise ratio (SNR) to that expected without perturbation by the pipeline. We inject simulated transit signals across the full focal plane for a set of observations of length 89 days. On average, we find that the SNR of the injected signal is recovered at MS = 0.9973(+/-0.0012)xBS-0.0151(+/-0.0049), where MS is the measured SNR and BS is the baseline, or expected, SNR. The 1{\sigma} width of the distribution around this correlation is +/-2.64%. We discuss the pipeline processes that cause the measured SNR to deviate significantly from the baseline SNR; these are primarily the handling of data adjacent to spacecraft re-pointings and the removal of harmonics prior to the measurement of the SNR. Finally we outline the further work required to characterise the completeness of the Kepler pipeline.
[74]  [pdf] - 1164538
Kepler-68: Three Planets, One With a Density Between That of Earth and Ice Giants
Comments: 32 pages, 13 figures, Accepted to ApJ
Submitted: 2013-02-11
NASA's Kepler Mission has revealed two transiting planets orbiting Kepler-68. Follow-up Doppler measurements have established the mass of the innermost planet and revealed a third jovian-mass planet orbiting beyond the two transiting planets. Kepler-68b, in a 5.4 day orbit has mass 8.3 +/- 2.3 Earth, radius 2.31 +/- 0.07 Earth radii, and a density of 3.32 +/- 0.92 (cgs), giving Kepler-68b a density intermediate between that of the ice giants and Earth. Kepler-68c is Earth-sized with a radius of 0.953 Earth and transits on a 9.6 day orbit; validation of Kepler-68c posed unique challenges. Kepler-68d has an orbital period of 580 +/- 15 days and minimum mass of Msin(i) = 0.947 Jupiter. Power spectra of the Kepler photometry at 1-minute cadence exhibit a rich and strong set of asteroseismic pulsation modes enabling detailed analysis of the stellar interior. Spectroscopy of the star coupled with asteroseismic modeling of the multiple pulsation modes yield precise measurements of stellar properties, notably Teff = 5793 +/- 74 K, M = 1.079 +/- 0.051 Msun, R = 1.243 +/- 0.019 Rsun, and density 0.7903 +/- 0.0054 (cgs), all measured with fractional uncertainties of only a few percent. Models of Kepler-68b suggest it is likely composed of rock and water, or has a H and He envelope to yield its density of about 3 (cgs).
[75]  [pdf] - 1152390
The Herschel view of the environment of the radio galaxy 4C+41.17 at z = 3.8
Comments: 17 pages, accepted for publication in the Monthly Notices of the Royal Astronomical Society
Submitted: 2012-10-23
We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.
[76]  [pdf] - 1151608
Jet and Torus Orientations in High Redshift Radio Galaxies
Comments: 18 pages, 12 figures, 6 tables
Submitted: 2012-09-24
We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.
[77]  [pdf] - 1151163
Jet interactions with a giant molecular cloud in the Galactic centre and ejection of hypervelocity stars
Comments: 4 pages, 1 figure. Astronomy and Astrophysics Letters, in press
Submitted: 2012-09-06
The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue that these stars, {as well as many more abundant bound OB stars in the innermost few parsecs,} were generated by the interactions of an AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium.
[78]  [pdf] - 1151059
Extreme host galaxy growth in powerful early-epoch radio galaxies
Comments: To appear in Astrophysical Journal Letters
Submitted: 2012-09-03
During the first half of the universe's age, a heyday of star-formation must have occurred because many massive galaxies are in place after that epoch in cosmic history. Our observations with the revolutionary Herschel Space Observatory reveal vigorous optically obscured star-formation in the ultra-massive hosts of many powerful high-redshift 3C quasars and radio galaxies. This symbiotic occurrence of star-formation and black hole driven activity is in marked contrast to recent results dealing with Herschel observations of X-ray selected active galaxies. Three archetypal radio galaxies, at redshifts 1.132,1.575, and 2.474 are presented here, with inferred star-formation rates of hundreds of solar masses per year. A series of spectacular coeval AGN/starburst events may have formed these ultra-massive galaxies and their massive central black holes during their relatively short lifetimes.
[79]  [pdf] - 1124380
Rapid Coeval Black Hole and Host Galaxy Growth in MRC 1138-262: The Hungry Spider
Comments: 11 pages, 2 figures, accepted for publication in ApJ
Submitted: 2012-06-25
We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.
[80]  [pdf] - 1124273
Kepler-36: A Pair of Planets with Neighboring Orbits and Dissimilar Densities
Comments: Accepted for publication in Science. Published online on June 21, 2012. Main Text and supplemental information included in a single merged file, 69 pages. Attachments to the supplemental material are available for free on Science website
Submitted: 2012-06-20
In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.
[81]  [pdf] - 1124200
Gas-rich mergers and feedback are ubiqitous amongst starbursting radio galaxies, as revealed by JVLA, IRAM PdBI and Herschel
Comments: 12 pages, 6 colour figures; accepted for publication in MNRAS
Submitted: 2012-06-18
We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.
[82]  [pdf] - 1123581
Toward Understanding the Origin of Turbulence in Molecular Clouds: Small Scale Structures as Units of Dynamical Multi-Phase Interstellar Medium
Comments: 12 pages, 9 figures. To be published in the Astrophysical Journal
Submitted: 2012-05-22
In order to investigate the origin of the interstellar turbulence, detailed observations in the CO J=1--0 and 3--2 lines have been carried out in an interacting region of a molecular cloud with an HII region. As a result, several 1,000 to 10,000 AU scale cloudlets with small velocity dispersion are detected, whose systemic velocities have a relatively large scatter of a few km/s. It is suggested that the cloud is composed of small-scale dense and cold structures and their overlapping effect makes it appear to be a turbulent entity as a whole. This picture strongly supports the two-phase model of turbulent medium driven by thermal instability proposed previously. On the surface of the present cloud, the turbulence is likely to be driven by thermal instability following ionization shock compression and UV irradiation. Those small scale structures with line width of ~ 0.6 km/s have a relatively high CO line ratio of J=3--2 to 1--0, 1 < R(3-2/1-0) < 2. The large velocity gradient analysis implies that the 0.6 km/s width component cloudlets have an average density of 10^{3-4} cm^{-3}, which is relatively high at cloud edges, but their masses are only < 0.05 M_{sun}.
[83]  [pdf] - 501954
The Transiting Circumbinary Planets Kepler-34 and Kepler-35
Comments: Accepted and published in Nature (2012 Jan 26). This is the submitted version of paper, merged with the Supplementary Information; 56 pages total with 20 figures
Submitted: 2012-04-17
Most Sun-like stars in the Galaxy reside in gravitationally-bound pairs of stars called "binary stars". While long anticipated, the existence of a "circumbinary planet" orbiting such a pair of normal stars was not definitively established until the discovery of Kepler-16. Incontrovertible evidence was provided by the miniature eclipses ("transits") of the stars by the planet. However, questions remain about the prevalence of circumbinary planets and their range of orbital and physical properties. Here we present two additional transiting circumbinary planets, Kepler-34 and Kepler-35. Each is a low-density gas giant planet on an orbit closely aligned with that of its parent stars. Kepler-34 orbits two Sun-like stars every 289 days, while Kepler-35 orbits a pair of smaller stars (89% and 81% of the Sun's mass) every 131 days. Due to the orbital motion of the stars, the planets experience large multi-periodic variations in incident stellar radiation. The observed rate of circumbinary planets implies > ~1% of close binary stars have giant planets in nearly coplanar orbits, yielding a Galactic population of at least several million.
[84]  [pdf] - 1093259
Transit Timing Observations from Kepler: IV. Confirmation of 4 Multiple Planet Systems by Simple Physical Models
Comments: In the proofs process, corrections were made to tables -- most crucially, the timing data for Kepler-30b and the depths and radii of planets in Kepler-31 and 32
Submitted: 2012-01-25, last modified: 2012-04-02
Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets' masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems' architectures, even in cases for which high-precision Doppler follow-up is impractical.
[85]  [pdf] - 1116907
Planetary Candidates Observed by Kepler, III: Analysis of the First 16 Months of Data
Comments: Submitted to ApJS. Machine-readable tables are available at,, and the NASA Exoplanet Archive
Submitted: 2012-02-27
New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.
[86]  [pdf] - 1093257
Transit Timing Observations from Kepler: II. Confirmation of Two Multiplanet Systems via a Non-parametric Correlation Analysis
Comments: 23 pages, 8 figures, 4 tables, 1 electronic table, accepted to ApJ
Submitted: 2012-01-25
We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.
[87]  [pdf] - 967262
Detection of Potential Transit Signals in the First Three Quarters of Kepler Mission Data
Submitted: 2012-01-04, last modified: 2012-01-18
We present the results of a search for potential transit signals in the first three quarters of photometry data acquired by the Kepler Mission. The targets of the search include 151,722 stars which were observed over the full interval and an additional 19,132 stars which were observed for only 1 or 2 quarters. From this set of targets we find a total of 5,392 detections which meet the Kepler detection criteria: those criteria are periodicity of the signal, an acceptable signal-to-noise ratio, and a composition test which rejects spurious detections which contain non-physical combinations of events. The detected signals are dominated by events with relatively low signal-to-noise ratio and by events with relatively short periods. The distribution of estimated transit depths appears to peak in the range between 40 and 100 parts per million, with a few detections down to fewer than 10 parts per million. The detected signals are compared to a set of known transit events in the Kepler field of view which were derived by a different method using a longer data interval; the comparison shows that the current search correctly identified 88.1% of the known events. A tabulation of the detected transit signals, examples which illustrate the analysis and detection process, a discussion of future plans and open, potentially fruitful, areas of further research are included.
[88]  [pdf] - 1092280
Kepler-21b: A 1.6REarth Planet Transiting the Bright Oscillating F Subgiant Star HD 179070
Comments: Accepted to ApJ
Submitted: 2011-12-09
We present Kepler observations of the bright (V=8.3), oscillating star HD 179070. The observations show transit-like events which reveal that the star is orbited every 2.8 days by a small, 1.6 R_Earth object. Seismic studies of HD 179070 using short cadence Kepler observations show that HD 179070 has a frequencypower spectrum consistent with solar-like oscillations that are acoustic p-modes. Asteroseismic analysis provides robust values for the mass and radius of HD 179070, 1.34{\pm}0.06 M{\circ} and 1.86{\pm}0.04 R{\circ} respectively, as well as yielding an age of 2.84{\pm}0.34 Gyr for this F5 subgiant. Together with ground-based follow-up observations, analysis of the Kepler light curves and image data, and blend scenario models, we conservatively show at the >99.7% confidence level (3{\sigma}) that the transit event is caused by a 1.64{\pm}0.04 R_Earth exoplanet in a 2.785755{\pm}0.000032 day orbit. The exoplanet is only 0.04 AU away from the star and our spectroscopic observations provide an upper limit to its mass of ~10 M_Earth (2-{\sigma}). HD 179070 is the brightest exoplanet host star yet discovered by Kepler.
[89]  [pdf] - 1092227
Kepler-22b: A 2.4 Earth-radius Planet in the Habitable Zone of a Sun-like Star
Comments: Accepted to ApJ
Submitted: 2011-12-07
A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.
[90]  [pdf] - 1083059
The hot-Jupiter Kepler-17b: discovery, obliquity from stroboscopic starspots, and atmospheric characterization
Comments: 21 pages, 13 figures, Accepted for publication in ApJ on October 14, 2011
Submitted: 2011-07-28, last modified: 2011-10-26
This paper reports the discovery and characterization of the transiting hot giant exoplanet Kepler-17b. The planet has an orbital period of 1.486 days, and radial velocity measurements from the Hobby-Eberly Telescope (HET) show a Doppler signal of 420+/-15 m.s-1. From a transit-based estimate of the host star's mean density, combined with an estimate of the stellar effective temperature T_eff=5630+/-100 K from high-resolution spectra, we infer a stellar host mass of 1.061+/-0.067 M_sun and a stellar radius of 1.019+/-0.033 R_jup. We estimate the planet mass and radius to be Mp=2.450+/-0.114 M_jup and Rp=1.312+/-0.018 R_jup and a planet density near 1.35 The host star is active, with dark spots that are frequently occulted by the planet. The continuous monitoring of the star reveals a stellar rotation period of 11.89 days, 8 times the the planet's orbital period; this period ratio produces stroboscopic effects on the occulted starspots. The temporal pattern of these spot-crossing events shows that the planet's orbit is prograde and the star's obliquity is smaller than 15 deg. We detected planetary occultations of Kepler-17b with both the Kepler and Spitzer Space Telescopes. We use these observations to constrain the eccentricity, e, and find that it is consistent with a circular orbit (e<0.0011). The brightness temperatures of the planet the infrared bandpasses are T_3.6um=1880+/-100 K and T4.5um=1770+/-150 K. We measure the optical geometric albedo A_g in the Kepler bandpass and find A_g = 0.10+/-0.02. The observations are best described by atmospheric models for which most of the incident energy is re-radiated away from the day side.
[91]  [pdf] - 423238
Observational Constraints, Stellar Models, and Kepler Data for theta Cyg, the Brightest Star Observable by Kepler
Comments: Proceedings submission for Resolving the Future of Long-Baseline Interferometry, Socorro, NM, March 28-31, 2011, 10 pages
Submitted: 2011-10-10
The V=4.48 F4 main-sequence star theta Cyg is the brightest star observable in the Kepler spacecraft field-of-view. Short-cadence (58.8 s) photometric data were obtained by Kepler during 2010 June-September. Preliminary analysis shows solar-like oscillations in the frequency range 1200- 2500 microHz. To interpret these data and to motivate further observations, we use observational constraints from the literature to construct stellar evolution and pulsation models for this star. We compare the observed large frequency separation of the solar-like oscillations with the model predictions, and discuss the prospects for gamma Doradus-like g-mode pulsations, given the observational constraints. We discuss the value of angular diameter measurements from optical interferometry for constraining stellar properties and the implications for asteroseismology.
[92]  [pdf] - 1084622
Kepler 18-b, c, and d: A System Of Three Planets Confirmed by Transit Timing Variations, Lightcurve Validation, Spitzer Photometry and Radial Velocity Measurements
Comments: ApJS in press
Submitted: 2011-10-04
We report the detection of three transiting planets around a Sunlike star, which we designate Kepler-18. The transit signals were detected in photometric data from the Kepler satellite, and were confirmed to arise from planets using a combination of large transit-timing variations, radial-velocity variations, Warm-Spitzer observations, and statistical analysis of false-positive probabilities. The Kepler-18 star has a mass of 0.97M_sun, radius 1.1R_sun, effective temperature 5345K, and iron abundance [Fe/H]= +0.19. The planets have orbital periods of approximately 3.5, 7.6 and 14.9 days. The innermost planet "b" is a "super-Earth" with mass 6.9 \pm 3.4M_earth, radius 2.00 \pm 0.10R_earth, and mean density 4.9 \pm 2.4 g cm^-3. The two outer planets "c" and "d" are both low-density Neptune-mass planets. Kepler-18c has a mass of 17.3 \pm 1.9M_earth, radius 5.49 \pm 0.26R_earth, and mean density 0.59 \pm 0.07 g cm^-3, while Kepler-18d has a mass of 16.4 \pm 1.4M_earth, radius 6.98 \pm 0.33R_earth, and mean density 0.27 \pm 0.03 g cm^-3. Kepler-18c and Kepler-18d have orbital periods near a 2:1 mean-motion resonance, leading to large and readily detected transit timing variations.
[93]  [pdf] - 1084255
Chandra X-ray Observations of the redshift 1.53, radio-loud quasar: 3C 270.1
Comments: 20 pages, 7 figures, accepted by ApJ
Submitted: 2011-09-20, last modified: 2011-09-21
Chandra X-ray observations of the high redshift (z =1.532) radio-loud quasar 3C270.1 in 2008 February show the nucleus to have a power-law spectrum, Gamma = 1.66 +/- 0.08, typical of a radio-loud quasar, and a marginally-detected Fe Kalpha emission line. The data also reveal extended X-ray emission, about half of which is associated with the radio emission from this source. The southern emission is co-spatial with the radio lobe and peaks at the position of the double radio hotspot. Modeling this hotspot including Spitzer upper limits rules out synchrotron emission from a single power-law population of electrons, favoring inverse-Compton emission with a field of ~11nT, roughly a third of the equipartition value. The northern emission is concentrated close to the location of a 40 deg. bend where the radio jet is presumed to encounter external material. It can be explained by inverse Compton emission involving Cosmic Microwave Background photons with a field of ~3nT, roughly a factor of nine below the equipartition value. The remaining, more diffuse X-ray emission is harder (HR=-0.09 +/- 0.22). With only 22.8+/-5.6 counts, the spectral form cannot be constrained. Assuming thermal emission with a temperature of 4 keV yields an estimate for the luminosity of 1.8E44 erg/s, consistent with the luminosity-temperature relation of lower-redshift clusters. However deeper Chandra X-ray observations are required to delineate the spatial distribution, and better constrain the spectrum of the diffuse emission to verify that we have detected X-ray emission from a high-redshift cluster.
[94]  [pdf] - 518887
Disentangling galaxy environment and host halo mass
Comments: 15 pages, 9 figures, 2 tables. Accepted by MNRAS
Submitted: 2011-03-02, last modified: 2011-09-20
[Abridged] The properties of observed galaxies and dark matter haloes in simulations depend on their environment. The term environment has been used to describe a wide variety of measures that may or may not correlate with each other. Popular measures of environment include the distance to the N'th nearest neighbour, the number density of objects within some distance, or the mass of the host dark matter halo. We use results from the Millennium simulation and a semi-analytic model for galaxy formation to quantify the relations between environment and halo mass. We show that the environmental parameters used in the observational literature are in effect measures of halo mass, even if they are measured for a fixed stellar mass. The strongest correlation between environment and halo mass arises when the number of objects is counted out to a distance of 1.5-2 times the virial radius of the host halo and when the galaxies/haloes are required to be relatively bright/massive. For observational studies the virial radius is not easily determined, but the number of neighbours out to 1-2 Mpc/h gives a similarly strong correlation. For the distance to the N'th nearest neighbour the correlation with halo mass is nearly as strong provided N>2. We demonstrate that this environmental parameter becomes insensitive to halo mass if it is constructed from dimensionless quantities. This can be achieved by scaling the minimum luminosity/mass of neighbours to that of the object in question and by dividing the distance to a length scale associated with either the neighbour or the galaxy under consideration. We show how such a halo mass independent environmental parameter can be defined for observational and numerical studies. The results presented here will help future studies to disentangle the effects of halo mass and external environment on the properties of galaxies and dark matter haloes.
[95]  [pdf] - 1084125
Kepler-16: A Transiting Circumbinary Planet
Comments: Science, in press; for supplemental material see
Submitted: 2011-09-15
We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.
[96]  [pdf] - 1083943
Photometric AGN reverberation mapping - an efficient tool for BLR sizes, black hole masses and host-subtracted AGN luminosities
Comments: 12 pages, 15 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2011-09-08
Photometric reverberation mapping employs a wide bandpass to measure the AGN continuum variations and a suitable band, usually a narrow band (NB), to trace the echo of an emission line in the broad line region (BLR). The narrow band catches both the emission line and the underlying continuum, and one needs to extract the pure emission line light curve. We performed a test on two local AGNs, PG0003+199 (=Mrk335) and Ark120, observing well-sampled broad- (B, V) and narrow-band light curves with the robotic 15cm telescope VYSOS-6 on Cerro Armazones, Chile. In PG0003+199, H_alpha dominates the flux in the NB by 85%, allowing us to measure the time lag of H_alpha against B without the need to correct for the continuum contribution. In Ark120, H_beta contributes only 50% to the flux in the NB. The cross correlation of the B and NB light curves shows two distinct peaks of similar strength, one at lag zero from the autocorrelated continuum and one from the emission line at tau_cent = 47.5 +/- 3.4 days. We constructed a synthetic H_beta light curve, by subtracting a scaled V light curve, which traces the continuum, from the NB light curve. The cross correlation of this synthetic H_beta light curve with the B light curve shows only one major peak at tau_cent = 48.0 +/- 3.3 days, while the peak from the autocorrelated continuum at lag zero is absent. We conclude that, as long as the emission line contributes at least 50% to the bandpass, the pure emission line light curve can be reconstructed from photometric monitoring data so that the time lag can be measured. For both objects the lags we find are consistent with spectroscopic reverberation results. While the dense sampling (median 2 days) enables us to determine tau_cent with small (10%) formal errors, we caution that gaps in the light curves may lead to much larger systematic uncertainties. (Abstract shortened, see the manuscript.)
[97]  [pdf] - 1083916
Discovery and Atmospheric Characterization of Giant Planet Kepler-12b: An Inflated Radius Outlier
Comments: Revised for ApJ
Submitted: 2011-09-07
We report the discovery of planet Kepler-12b (KOI-20), which at 1.695\pm0.030 RJ is among the handful of planets with super-inflated radii above 1.65 RJ. Orbiting its slightly evolved G0 host with a 4.438-day period, this 0.431\pm0.041 MJ planet is the least-irradiated within this largest-planet-radius group, which has important implications for planetary physics. The planet's inflated radius and low mass lead to a very low density of 0.111\pm0.010 g cm-3. We detect the occultation of the planet at a significance of 3.7{\sigma} in the Kepler bandpass. This yields a geometric albedo of 0.14\pm0.04; the planetary flux is due to a combination of scattered light and emitted thermal flux. We use multiple observations with Warm Spitzer to detect the occultation at 7{\sigma} and 4{\sigma} in the 3.6 and 4.5 {\mu}m bandpasses, respectively. The occultation photometry timing is consistent with a circular orbit, at e < 0.01 (1{\sigma}), and e < 0.09 (3{\sigma}). The occultation detections across the three bands favor an atmospheric model with no dayside temperature inversion. The Kepler occultation detection provides significant leverage, but conclusions regarding temperature structure are preliminary, given our ignorance of opacity sources at optical wavelengths in hot Jupiter atmospheres. If Kepler-12b and HD 209458b, which intercept similar incident stellar fluxes, have the same heavy element masses, the interior energy source needed to explain the large radius of Kepler-12b is three times larger than that of HD 209458b. This may suggest that more than one radius-inflation mechanism is at work for Kepler-12b, or that it is less heavy-element rich than other transiting planets.
[98]  [pdf] - 1083909
The Kepler-19 System: A Transiting 2.2 R_Earth Planet and a Second Planet Detected via Transit Timing Variations
Comments: 50 pages, 15 figures, accepted for publication in ApJ
Submitted: 2011-09-07
We present the discovery of the Kepler-19 planetary system, which we first identified from a 9.3-day periodic transit signal in the Kepler photometry. From high-resolution spectroscopy of the star, we find a stellar effective temperature Teff=5541 \pm 60 K, a metallicity [Fe/H]=-0.13 \pm 0.06, and a surface gravity log(g)=4.59 \pm 0.10. We combine the estimate of Teff and [Fe/H] with an estimate of the stellar density derived from the photometric light curve to deduce a stellar mass of M_star = 0.936 \pm 0.040 M_Sun and a stellar radius of R_star = 0.850 \pm 0.018 R_Sun. We rule out the possibility that the transits result from an astrophysical false positive by first identifying the subset of stellar blends that reproduce the precise shape of the light curve. We conclude that the planetary scenario is more than three orders of magnitude more likely than a blend. The blend scenario is independently disfavored by the achromaticity of the transit: we measure a transit depth with Spitzer at 4.5 {\mu}m of 547+113-110 ppm, consistent with the depth measured in the Kepler optical bandpass of 567\pm6 ppm. We determine a physical radius of the planet Kepler-19b of R_p = 2.209 \pm 0.048 R_Earth. From radial-velocity observations of the star, we find an upper limit on the planet mass of 20.3 M_Earth, corresponding to a maximum density of 10.4 g cm^-3. We report a significant sinusoidal deviation of the transit times from a predicted linear ephemeris, which we conclude is due to an additional perturbing body in the system. We cannot uniquely determine the orbital parameters of the perturber, as various dynamical mechanisms match the amplitude, period, and shape of the transit timing signal and satisfy the host star's radial velocity limits. However, the perturber in these mechanisms has period <160 days and mass <6 M_Jup, confirming its planetary nature as Kepler-19c. [Abridged]
[99]  [pdf] - 384954
The First Kepler Mission Planet Confirmed With The Hobby-Eberly Telescope: Kepler-15b, a Hot Jupiter Enriched In Heavy Elements
Comments: preprint, 33 pages, 12 figures, submitted to ApJ
Submitted: 2011-07-13
We report the discovery of Kepler-15b, a new transiting exoplanet detected by NASA's Kepler mission. The transit signal with a period of 4.94 days was detected in the quarter 1 (Q1) Kepler photometry. For the first time, we have used the High-Resolution-Spectrograph (HRS) at the Hobby-Eberly Telescope (HET) to determine the mass of a Kepler planet via precise radial velocity (RV) measurements. The 24 HET/HRS radial velocities (RV) and 6 additional measurements from the FIES spectrograph at the Nordic Optical Telescope (NOT) reveal a Doppler signal with the same period and phase as the transit ephemeris. We used one HET/HRS spectrum of Kepler-15 taken without the iodine cell to determine accurate stellar parameters. The host star is a metal-rich ([Fe/H]=0.36+/-0.07) G-type main sequence star with T_eff=5515+/-124 K. The amplitude of the RV-orbit yields a mass of the planet of 0.66+/-0.1 M_Jup. The planet has a radius of 0.96+/-0.06 R_Jup and a mean bulk density of 0.9+/-0.2 g/cm^3. The planetary radius resides on the lower envelope for transiting planets with similar mass and irradiation level. This suggests significant enrichment of the planet with heavy elements. We estimate a heavy element mass of 30-40 M_Earth within Kepler-15b.
[100]  [pdf] - 1077592
Kepler-14b: A massive hot Jupiter transiting an F star in a close visual binary
Comments: 10 pages, 8 figures, submitted to ApJ
Submitted: 2011-06-27
We present the discovery of a hot Jupiter transiting an F star in a close visual (0.3" sky projected angular separation) binary system. The dilution of the host star's light by the nearly equal magnitude stellar companion (~ 0.5 magnitudes fainter) significantly affects the derived planetary parameters, and if left uncorrected, leads to an underestimate of the radius and mass of the planet by 10% and 60%, respectively. Other published exoplanets, which have not been observed with high-resolution imaging, could similarly have unresolved stellar companions and thus have incorrectly derived planetary parameters. Kepler-14b (KOI-98) has a period of P = 6.790 days and correcting for the dilution, has a mass of Mp = 8.40 +0.19-0.18 MJ and a radius of Rp = 1.136 +0.073-0.054 RJ, yielding a mean density of rho = 7.1 +- 1.1 g cm-3.
[101]  [pdf] - 522125
The correlation structure of dark matter halo properties
Comments: 5 pages, 2 figures. Accepted by MNRAS Letters. Minor changes
Submitted: 2011-03-28, last modified: 2011-05-25
We investigate the correlation between nine different dark matter halo properties using a rank correlation analysis and a Principal Component Analysis for a sample of haloes spanning five orders of magnitude in mass. We consider mass and dimensionless measures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and environment, where the latter is defined in a way that makes it insensitive to mass. We find that concentration is the most fundamental property. Except for environment, all parameters are strongly correlated with concentration. Concentration, age, substructure, mass, sphericity and relaxedness can be considered a single family of parameters, albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more independent, although spin does correlate strongly with substructure and both spin and triaxiality correlate substantially with concentration. Although mass sets the scale of a halo, all other properties are more sensitive to concentration.
[102]  [pdf] - 1076360
Quantifying the anisotropy in the infrared emission of powerful AGN
Comments: 6 pages, 2 figures; accepted by ApJ
Submitted: 2011-05-02
We use restframe near- and mid-IR data of an isotropically selected sample of quasars and radio galaxies at 1.0 \leq z \leq 1.4, which have been published previously, to study the wavelength-dependent anisotropy of the IR emission. For that we build average SEDs of the quasar subsample (= type 1 AGN) and radio galaxies (= type 2 AGN) from ~1-17 {\mu}m and plot the ratio of both average samples. From 2 to 8 {\mu}m restframe wavelength the ratio gradually decreases from 20 to 2 with values around 3 in the 10{\mu}m silicate feature. Longward of 12{\mu}m the ratio decreases further and shows some high degree of isotropy at 15 {\mu}m (ratio ~1.4). The results are consistent with upper limits derived from the X-ray/mid-IR correlation of local Seyfert galaxies. We find that the anisotropy in our high-luminosity radio-loud sample is smaller than in radio-quiet lower-luminosity AGN which may be interpreted in the framework of a receding torus model with luminosity-dependent obscuration properties. It is also shown that the relatively small degree of anisotropy is consistent with clumpy torus models.
[103]  [pdf] - 1053392
The Kepler Cluster Study: Stellar Rotation in NGC6811
Comments: 13 pages (aastex 12pt preprint style), 4 figures, 1 table, accepted for publication in ApJ Letters
Submitted: 2011-04-14
We present rotation periods for 71 single dwarf members of the open cluster NGC6811 determined using photometry from NASA's Kepler Mission. The results are the first from The Kepler Cluster Study which combine Kepler's photometry with ground-based spectroscopy for cluster membership and binarity. The rotation periods delineate a tight sequence in the NGC6811 color-period diagram from ~1 day at mid-F to ~11 days at early-K spectral type. This result extends to ~1 Gyr similar prior results in the ~600 Myr Hyades and Praesepe clusters, suggesting that rotation periods for cool dwarf stars delineate a well-defined surface in the 3-dimensional space of color (mass), rotation, and age. It implies that reliable ages can be derived for field dwarf stars with measured colors and rotation periods, and it promises to enable further understanding of various aspects of stellar rotation and activity for cool stars.
[104]  [pdf] - 1052692
Planet Occurrence within 0.25 AU of Solar-type Stars from Kepler
Comments: Submitted to ApJ, 22 pages, 10 figures
Submitted: 2011-03-13
We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]
[105]  [pdf] - 2096108
Characteristics of planetary candidates observed by Kepler, II: Analysis of the first four months of data
Comments: 106 pages, 15 figures, contains tables of candidates. Submitted to ApJ
Submitted: 2011-02-02, last modified: 2011-03-12
On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.
[106]  [pdf] - 1041881
The rates and modes of gas accretion on to galaxies and their gaseous haloes
Comments: Accepted for publication in MNRAS, 23 pages and 18 figures. Revised version: minor changes
Submitted: 2010-11-10, last modified: 2011-02-18
(Abridged) We study the rate at which gas accretes onto galaxies and haloes and investigate whether the accreted gas was shocked to high temperatures before reaching a galaxy. For this purpose we use a suite of large cosmological, hydrodynamical simulations from the OWLS project. We improve on previous work by considering a wider range of halo masses and redshifts, by distinguishing accretion onto haloes and galaxies, by including important feedback processes, and by comparing simulations with different physics. The specific rate of gas accretion onto haloes is, like that for dark matter, only weakly dependent on halo mass. For halo masses Mhalo>>10^11 Msun it is relatively insensitive to feedback processes. In contrast, accretion rates onto galaxies are determined by radiative cooling and by outflows driven by supernovae and active galactic nuclei. Galactic winds increase the halo mass at which the central galaxies grow the fastest by about two orders of magnitude to Mhalo~10^12 Msun. Gas accretion is bimodal, with maximum past temperatures either of order the virial temperature or <~10^5 K. The fraction of gas accreted on to haloes in the hot mode is insensitive to feedback and metal-line cooling. It increases with decreasing redshift, but is mostly determined by halo mass, increasing gradually from less than 10% for ~10^11 Msun to greater than 90% at 10^13 Msun. In contrast, for accretion onto galaxies the cold mode is always significant and the relative contributions of the two accretion modes are more sensitive to feedback and metal-line cooling. The majority of stars present in any mass halo at any redshift were formed from gas accreted in the cold mode, although the hot mode contributes typically over 10% for Mhalo>~10^11 Msun. Galaxies, but not necessarily their gaseous haloes, are predominantly fed by gas that did not experience an accretion shock when it entered the host halo.
[107]  [pdf] - 1052059
Megamaser detection and nuclear obscuration in Seyfert galaxies
Submitted: 2011-02-11
We revisit the relation between H2O maser detection rate and nuclear obscuration for a sample of 114 Seyfert galaxies, drawn from the CfA, 12um and IRAS F25/F60 catalogs. These sources have mid-infrared spectra from the Spitzer Space Telescope and they are searched for X-ray and [O III], 5007Angstrom fluxes from the literature. We use the strength of the [O IV], 25.9um emission line as tracer for the intrinsic AGN strength. After normalization by [O IV] the observed X-ray flux provides information about X-ray absorption. The distribution of X-ray / [O IV] flux ratios is significantly different for masers and non-masers: The maser detected Seyfert-2s (Sy 1.8-2.0) populate a distinct X-ray / [O IV] range which is, on average, about a factor four lower than the range of Seyfert-2 non-masers and about a factor of ten lower than the range of Seyfert-1s (Sy 1.0-1.5). Non-masers are almost equally distributed over the entire X-ray / [O IV] range. This provides evidence that high nuclear obscuration plays a crucial role for the probability of maser detection. Furthermore, after normalization with [O IV], we find a similar but weaker trend for the distribution of the maser detection rate with the absorption of the 7um dust continuum. This suggests that the obscuration of the 7 um continuum occurs on larger spatial scales than that of the X-rays. Hence, in the AGN unified model, at moderate deviation from edge-on, the 7um dust absorption may occur without proportionate X-ray absorption. The absorption of [O III] appears unrelated to maser detections. The failure to detect masers in obscured AGN is most likely due to insufficient observational sensitivity.
[108]  [pdf] - 318272
KEPLER's First Rocky Planet: Kepler-10b
Comments: Accepted, Astrophysical Journal, November 25, 2010; Eexpected publication date: February 20, 2011
Submitted: 2011-02-03
NASA's Kepler Mission uses transit photometry to determine the frequency of earth-size planets in or near the habitable zone of Sun-like stars. The mission reached a milestone toward meeting that goal: the discovery of its first rocky planet, Kepler-10b. Two distinct sets of transit events were detected: 1) a 152 +/- 4 ppm dimming lasting 1.811 +/- 0.024 hours with ephemeris T[BJD]=2454964.57375+N*0.837495 days and 2) a 376 +/- 9 ppm dimming lasting 6.86 +/- 0.07 hours with ephemeris T[BJD]=2454971.6761+N*45.29485 days. Statistical tests on the photometric and pixel flux time series established the viability of the planet candidates triggering ground-based follow-up observations. Forty precision Doppler measurements were used to confirm that the short-period transit event is due to a planetary companion. The parent star is bright enough for asteroseismic analysis. Photometry was collected at 1-minute cadence for >4 months from which we detected 19 distinct pulsation frequencies. Modeling the frequencies resulted in precise knowledge of the fundamental stellar properties. Kepler-10 is a relatively old (11.9 +/- 4.5 Gyr) but otherwise Sun-like Main Sequence star with Teff=5627 +/- 44 K, Mstar=0.895 +/- 0.060 Msun, and Rstar=1.056 +/- 0.021 Rsun. Physical models simultaneously fit to the transit light curves and the precision Doppler measurements yielded tight constraints on the properties of Kepler-10b that speak to its rocky composition: Mpl=4.56 +/- 1.29 Mearth, Rpl=1.416 +/- 0.036 Rearth, and density=8.8 +/- 2.9 gcc. Kepler-10b is the smallest transiting exoplanet discovered to date.
[109]  [pdf] - 300715
KOI-126: A Triply-Eclipsing Hierarchical Triple with Two Low-Mass Stars
Comments: To be published in Science on 2/4/2011. Announced at Jan. 2011 AAS meeting and made available on Science Express. Includes Supporting Online Material
Submitted: 2011-02-02
The Kepler spacecraft has been monitoring the light from 150,000 stars in its primary quest to detect transiting exoplanets. Here we report on the detection of an eclipsing stellar hierarchical triple, identified in the Kepler photometry. KOI-126 (A,(B, C)), is composed of a low-mass binary (masses M_B = 0.2413+/-0.0030 M_Sun, M_C = 0.2127+/-0.0026 M_Sun; radii R_B = 0.2543+/-0.0014 R_Sun, R_C = 0.2318+/-0.0013 R_Sun; orbital period P_1 = 1.76713+/-0.00019 days) on an eccentric orbit about a third star (mass M_A = 1.347+/-0.032 M_Sun; radius R_A = 2.0254+/-0.0098 R_Sun; period of orbit around the low-mass binary P_2 = 33.9214+/-0.0013 days; eccentricity of that orbit e_2 = 0.3043+/-0.0024). The low-mass pair probe the poorly sampled fully-convective stellar domain offering a crucial benchmark for theoretical stellar models.
[110]  [pdf] - 1051826
A Closely-Packed System of Low-Mass, Low-Density Planets Transiting Kepler-11
Comments: published in Nature
Submitted: 2011-02-01
When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.
[111]  [pdf] - 270069
Five Kepler target stars that show multiple transiting exoplanet candidates
Comments: Accepted to ApJ
Submitted: 2010-06-14, last modified: 2010-11-09
We present and discuss five candidate exoplanetary systems identified with the Kepler spacecraft. These five systems show transits from multiple exoplanet candidates. Should these objects prove to be planetary in nature, then these five systems open new opportunities for the field of exoplanets and provide new insights into the formation and dynamical evolution of planetary systems. We discuss the methods used to identify multiple transiting objects from the Kepler photometry as well as the false-positive rejection methods that have been applied to these data. One system shows transits from three distinct objects while the remaining four systems show transits from two objects. Three systems have planet candidates that are near mean motion commensurabilities---two near 2:1 and one just outside 5:2. We discuss the implications that multitransiting systems have on the distribution of orbital inclinations in planetary systems, and hence their dynamical histories; as well as their likely masses and chemical compositions. A Monte Carlo study indicates that, with additional data, most of these systems should exhibit detectable transit timing variations (TTV) due to gravitational interactions---though none are apparent in these data. We also discuss new challenges that arise in TTV analyses due to the presence of more than two planets in a system.
[112]  [pdf] - 1033441
Mid-infrared spectroscopy of high-redshift 3CRR sources
Comments: 16 pages, 8 Figures, accepted for ApJ
Submitted: 2010-07-01
Using the Spitzer Space Telescope, we have obtained rest frame 9-16mu spectra of 11 quasars and 9 radio galaxies from the 3CRR catalog at redshifts 1.0<z<1.4. This complete flux-limited 178MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). The mean radio galaxy spectrum shows a clear silicate absorption feature (tau_9.7mu = 1.1) whereas the mean quasar spectrum shows silicates in emission. The mean radio galaxy spectrum matches a dust-absorbed mean quasar spectrum in both shape and overall flux level. The data for individual objects conform to these results. The trend of the silicate depth to increase with decreasing core fraction of the radio source further supports that for this sample, orientation is the main driver for the difference between radio galaxies and quasars, as predicted by AGN unification. However, comparing our high-z sample with lower redshift 3CRR objects reveals that the absorption of the high-z radio galaxy MIR continuum is lower than expected from a scaled up version of lower luminosity sources, and we discuss some effects that may explain these trends.
[113]  [pdf] - 1032785
Herschel/PACS far-infrared photometry of two z>4 quasars
Comments: 4 pages, 3 figures, accepted for publication in the A&A special issue on Herschel
Submitted: 2010-05-27
We present Herschel far-infrared (FIR) observations of two sub-mm bright quasars at high redshift: SDSS J1148+5251 (z=6.42) and BR 1202-0725 (z=4.69) obtained with the PACS instrument. Both objects are detected in the PACS photometric bands. The Herschel measurements provide additional data points that constrain the FIR spectral energy distributions (SEDs) of both sources, and they emphasise a broad range of dust temperatures in these objects. For lambda_rest ~< 20mu, the two SEDs are very similar to the average SEDs of quasars at low redshift. In the FIR, however, both quasars show excess emission compared to low-z QSO templates, most likely from cold dust powered by vigorous star formation in the QSO host galaxies. For SDSS J1148+5251 we detect another object at 160mu with a distance of ~10 arcseconds from the QSO. Although no physical connection between the quasar and this object can be shown with the available data, it could potentially confuse low-resolution measurements, thus resulting in an overestimate of the FIR luminosity of the z=6.42 quasar.
[114]  [pdf] - 1032443
Tracing the sites of obscured star formation in the Antennae galaxies with Herschel-PACS
Comments: 4 pages, 4 figures, 2 tables (A&A Herschel special issue)
Submitted: 2010-05-13
FIR imaging of interacting galaxies allows locating even hidden sites of star formation and measuring of the relative strength of nuclear and extra-nuclear star formation. We want to resolve the star-forming sites in the nearby system of the Antennae. Thanks to the unprecedented sharpness and depth of the PACS camera onboard ESA's Herschel Space Observatory, it is possible for the first time to achieve a complete assessment of individual star-forming knots in the FIR with scan maps at 70, 100, and 160 um. We used clump extraction photometry and SED diagnostics to derive the properties related to star-forming activity. The PACS 70, 100, and 160 um maps trace the knotty structure of the most recent star formation along an arc between the two nuclei in the overlap area. The resolution of the starburst knots and additional multi-wavelength data allow their individual star formation history and state to be analysed. In particular, the brightest knot in the mid-infrared (K1), east of the southern nucleus, exhibits the highest activity by far in terms of dust heating and star formation rate, efficiency, and density. With only 2 kpc in diameter, this area has a 10-1000 um luminosity, which is as high as that of our Milky Way. It shows the highest deficiency in radio emission in the radio-to-FIR luminosity ratio and a lack of X-ray emission, classifying it as a very young complex. The brightest 100 and 160 um emission region (K2), which is close to the collision front and consists of 3 knots, also shows a high star formation density and efficiency and lack of X-ray emission in its most obscured part, but an excess in the radio-to-FIR luminosity ratio. This suggests a young stage, too, but different conditions in its interstellar medium. Our results provide important checkpoints for numerical simulations of interacting galaxies when modelling the star formation and stellar feedback.
[115]  [pdf] - 1018961
The Kepler Pixel Response Function
Comments: 10 pages, 5 figures, accepted by ApJ Letters. Version accepted for publication.
Submitted: 2010-01-02, last modified: 2010-03-24
Kepler seeks to detect sequences of transits of Earth-size exoplanets orbiting Solar-like stars. Such transit signals are on the order of 100 ppm. The high photometric precision demanded by Kepler requires detailed knowledge of how the Kepler pixels respond to starlight during a nominal observation. This information is provided by the Kepler pixel response function (PRF), defined as the composite of Kepler's optical point spread function, integrated spacecraft pointing jitter during a nominal cadence and other systematic effects. To provide sub-pixel resolution, the PRF is represented as a piecewise-continuous polynomial on a sub-pixel mesh. This continuous representation allows the prediction of a star's flux value on any pixel given the star's pixel position. The advantages and difficulties of this polynomial representation are discussed, including characterization of spatial variation in the PRF and the smoothing of discontinuities between sub-pixel polynomial patches. On-orbit super-resolution measurements of the PRF across the Kepler field of view are described. Two uses of the PRF are presented: the selection of pixels for each star that maximizes the photometric signal to noise ratio for that star, and PRF-fitted centroids which provide robust and accurate stellar positions on the CCD, primarily used for attitude and plate scale tracking. Good knowledge of the PRF has been a critical component for the successful collection of high-precision photometry by Kepler.
[116]  [pdf] - 1018957
Kepler Mission Design, Realized Photometric Performance, and Early Science
Comments: 16 pages, 5 figures, 1 table 26 Jan revision replaced Subject headings with keywords from approved list
Submitted: 2010-01-01, last modified: 2010-01-26
The Kepler Mission, launched on Mar 6, 2009 was designed with the explicit capability to detect Earth-size planets in the habitable zone of solar-like stars using the transit photometry method. Results from just forty-three days of data along with ground-based follow-up observations have identified five new transiting planets with measurements of their masses, radii, and orbital periods. Many aspects of stellar astrophysics also benefit from the unique, precise, extended and nearly continuous data set for a large number and variety of stars. Early results for classical variables and eclipsing stars show great promise. To fully understand the methodology, processes and eventually the results from the mission, we present the underlying rationale that ultimately led to the flight and ground system designs used to achieve the exquisite photometric performance. As an example of the initial photometric results, we present variability measurements that can be used to distinguish dwarf stars from red giants.
[117]  [pdf] - 518880
Variations in Integrated Galactic Initial Mass Functions due to Sampling Method and Cluster Mass Function
Comments: Resubmitted to A&A, 14 pages, 9 Figures
Submitted: 2010-01-13
[abridged] Stars are thought to be formed predominantly in clusters. The clusters are formed following a cluster initial mass function (CMF) similar to the stellar initial mass function (IMF). Both the IMF and the CMF favour low-mass objects. The numerous low-mass clusters will lack high mass stars. If the integrated galactic initial mass function originates from stars formed in clusters, the IGIMF could be steeper than the IMF. We investigate how well constrained this steepening is and how it depends on the choice of sampling method and CMF. We compare analytic sampling to several implementations of random sampling of the IMF, and different CMFs. We implement different IGIMFs into GALEV to obtain colours and metallicities for galaxies. Choosing different ways of sampling the IMF results in different IGIMFs. Depending on the lower cluster mass limit and the slope of the cluster mass function, the steepening varies between very strong and negligible. We find the size of the effect is continuous as a function of the power-law slope of the CMF, if the CMF extends to masses smaller than the maximum stellarmass. The number of O-stars detected by GAIA might help in judging on the importance of the IGIMF effect. The impact of different IGIMFs on integrated galaxy photometry is small, within the intrinsic scatter of observed galaxies. Observations of gas fractions and metallicities could rule out at least the most extreme sampling methods. As we still do not understand the details of star formation, one sampling method cannot be favoured over another. Also, the CMF at very low cluster masses is not well constrained observationally. These uncertainties need to be taken into account when using an IGIMF, with severe implications for galaxy evolution models and interpretations of galaxy observations.
[118]  [pdf] - 1018972
Discovery and Rossiter-McLaughlin Effect of Exoplanet Kepler-8b
Comments: 26 pages, 8 figures, 2 tables; In preparation for submission to the Astrophysical Journal
Submitted: 2010-01-04
We report the discovery and the Rossiter-McLaughlin effect of Kepler-8b, a transiting planet identified by the NASA Kepler Mission. Kepler photometry and Keck-HIRES radial velocities yield the radius and mass of the planet around this F8IV subgiant host star. The planet has a radius RP = 1.419 RJ and a mass, MP = 0.60 MJ, yielding a density of 0.26 g cm^-3, among the lowest density planets known. The orbital period is P = 3.523 days and orbital semima jor axis is 0.0483+0.0006/-0.0012 AU. The star has a large rotational v sin i of 10.5 +/- 0.7 km s^-1 and is relatively faint (V = 13.89 mag), both properties deleterious to precise Doppler measurements. The velocities are indeed noisy, with scatter of 30 m s^-1, but exhibit a period and phase consistent with the planet implied by the photometry. We securely detect the Rossiter-McLaughlin effect, confirming the planet's existence and establishing its orbit as prograde. We measure an inclination between the projected planetary orbital axis and the projected stellar rotation axis of lambda = -26.9 +/- 4.6 deg, indicating a moderate inclination of the planetary orbit. Rossiter-McLaughlin measurements of a large sample of transiting planets from Kepler will provide a statistically robust measure of the true distribution of spin-orbit orientations for hot jupiters in general.
[119]  [pdf] - 1018974
Kepler Science Operations
Comments: Recently submitted to ApJL
Submitted: 2010-01-03
Kepler's primary mission is a search for earth-size exoplanets in the habitable zone of late-type stars using the transit method. To effectively accomplish this mission, Kepler orbits the Sun and stares nearly continuously at one field-of-view which was carefully selected to provide an appropriate density of target stars. The data transmission rates, operational cycles, and target management requirements implied by this mission design have been optimized and integrated into a comprehensive plan for science operations. The commissioning phase completed all critical tasks and accomplished all objectives within a week of the pre-launch plan. Since starting science, the nominal data collection timeline has been interrupted by two safemode events, several losses of fine point, and some small pointing adjustments. The most important anomalies are understood and mitigated, so Kepler's technical performance metrics have improved significantly over this period and the prognosis for mission success is excellent. The Kepler data archive is established and hosting data for the science team, guest observers, and public. The first data sets to become publicly available include the monthly full-frame images, dropped targets, and individual sources as they are published. Data are released through the archive on a quarterly basis; the Kepler Results Catalog will be released annually starting in 2011.
[120]  [pdf] - 1018963
Selection, Prioritization, and Characteristics of Kepler Target Stars
Comments: Submitted to Astrophysical Journal Letters
Submitted: 2010-01-03
The Kepler Mission began its 3.5-year photometric monitoring campaign in May 2009 on a select group of approximately 150,000 stars. The stars were chosen from the ~half million in the field of view that are brighter than 16th magnitude. The selection criteria are quantitative metrics designed to optimize the scientific yield of the mission with regards to the detection of Earth-size planets in the habitable zone. This yields more than 90,000 G-type stars on or close to the Main Sequence, >20,000 of which are brighter than 14th magnitude. At the temperature extremes, the sample includes approximately 3,000 M-type dwarfs and a small sample of O and B-type MS stars <200. Small numbers of giants are included in the sample which contains ~5,000 stars with surface gravities log(g) < 3.5. We present a brief summary of the selection process and the stellar populations it yields in terms of surface gravity, effective temperature, and apparent magnitude. In addition to the primary, statistically-derived target set, several ancillary target lists were manually generated to enhance the science of the mission, examples being: known eclipsing binaries, open cluster members, and high proper-motion stars.
[121]  [pdf] - 1018950
Instrument Performance in Kepler's First Months
Comments: 5 pages, 2 figures, 2 tables, Astrophysical Journal Letters, accepted
Submitted: 2010-01-01
The Kepler Mission relies on precise differential photometry to detect the 80 parts per million (ppm) signal from an Earth-Sun equivalent transit. Such precision requires superb instrument stability on time scales up to ~2 days and systematic error removal to better than 20 ppm. To this end, the spacecraft and photometer underwent 67 days of commissioning, which included several data sets taken to characterize the photometer performance. Because Kepler has no shutter, we took a series of dark images prior to the dust cover ejection, from which we measured the bias levels, dark current, and read noise. These basic detector properties are essentially unchanged from ground-based tests, indicating that the photometer is working as expected. Several image artifacts have proven more complex than when observed during ground testing, as a result of their interactions with starlight and the greater thermal stability in flight, which causes the temperature-dependent artifact variations to be on the timescales of transits. Because of Kepler's unprecedented sensitivity and stability, we have also seen several unexpected systematics that affect photometric precision. We are using the first 43 days of science data to characterize these effects and to develop detection and mitigation methods that will be implemented in the calibration pipeline. Based on early testing, we expect to attain Kepler's planned photometric precision over 80%-90% of the field of view.
[122]  [pdf] - 1018953
Initial Characteristics of Kepler Long Cadence Data For Detecting Transiting Planets
Comments: 9 pages, 5 figures
Submitted: 2010-01-01
The Kepler Mission seeks to detect Earth-size planets transiting solar-like stars in its ~115 deg^2 field of view over the course of its 3.5 year primary mission by monitoring the brightness of each of ~156,000 Long Cadence stellar targets with a time resolution of 29.4 minutes. We discuss the photometric precision achieved on timescales relevant to transit detection for data obtained in the 33.5-day long Quarter 1 (Q1) observations that ended 2009 June 15. The lower envelope of the photometric precision obtained at various timescales is consistent with expected random noise sources, indicating that Kepler has the capability to fulfill its mission. The Kepler light curves exhibit high precision over a large dynamic range, which will surely permit their use for a large variety of investigations in addition to finding and characterizing planets. We discuss the temporal characteristics of both the raw flux time series and the systematic error-corrected flux time series produced by the Kepler Science Pipeline, and give examples illustrating Kepler's large dynamic range and the variety of light curves obtained from the Q1 observations.
[123]  [pdf] - 1018954
Overview of the Kepler Science Processing Pipeline
Comments: 8 pages, 3 figures
Submitted: 2010-01-01
The Kepler Mission Science Operations Center (SOC) performs several critical functions including managing the ~156,000 target stars, associated target tables, science data compression tables and parameters, as well as processing the raw photometric data downlinked from the spacecraft each month. The raw data are first calibrated at the pixel level to correct for bias, smear induced by a shutterless readout, and other detector and electronic effects. A background sky flux is estimated from ~4500 pixels on each of the 84 CCD readout channels, and simple aperture photometry is performed on an optimal aperture for each star. Ancillary engineering data and diagnostic information extracted from the science data are used to remove systematic errors in the flux time series that are correlated with these data prior to searching for signatures of transiting planets with a wavelet-based, adaptive matched filter. Stars with signatures exceeding 7.1 sigma are subjected to a suite of statistical tests including an examination of each star's centroid motion to reject false positives caused by background eclipsing binaries. Physical parameters for each planetary candidate are fitted to the transit signature, and signatures of additional transiting planets are sought in the residual light curve. The pipeline is operational, finding planetary signatures and providing robust eliminations of false positives.
[124]  [pdf] - 1018946
Initial Characteristics of Kepler Short Cadence Data
Comments: 5 pages, 4 figures, ApJ Letters in press
Submitted: 2009-12-31
The Kepler Mission offers two options for observations -- either Long Cadence (LC) used for the bulk of core mission science, or Short Cadence (SC) which is used for applications such as asteroseismology of solar-like stars and transit timing measurements of exoplanets where the 1-minute sampling is critical. We discuss the characteristics of SC data obtained in the 33.5-day long Quarter 1 (Q1) observations with Kepler which completed on 15 June 2009. The truly excellent time series precisions are nearly Poisson limited at 11th magnitude providing per-point measurement errors of 200 parts-per-million per minute. For extremely saturated stars near 7th magnitude precisions of 40 ppm are reached, while for background limited measurements at 17th magnitude precisions of 7 mmag are maintained. We note the presence of two additive artifacts, one that generates regularly spaced peaks in frequency, and one that involves additive offsets in the time domain inversely proportional to stellar brightness. The difference between LC and SC sampling is illustrated for transit observations of TrES-2.
[125]  [pdf] - 341102
The physics driving the cosmic star formation history
Comments: Accepted for publication in MNRAS, 27 pages and 18 figures. Revised version: minor changes
Submitted: 2009-09-28, last modified: 2009-11-24
We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.
[126]  [pdf] - 1017682
Galactic consequences of clustered star formation
Comments: 4 pages, 2 figures, to appear in online version of Proceedings of IAU S266, a two page version will appear in the Proceedings of IAU S262
Submitted: 2009-09-24
If all stars form in clusters and both the stars and the clusters follow a power law distribution which favours the creation of low mass objects, then the numerous low mass clusters will be deficient in high mass stars. Therefore, the mass function of stars, integrated over the whole galaxy (the Integrated Galactic Initial Mass Function, IGIMF) will be steeper at the high mass end than the underlying IMF of the stars. We show how the steepness of the IGIMF depends on the sampling method and on the assumptions made for the star cluster mass function. We also investigate the O-star content, integrated photometry and chemical enrichment of galaxies that result from several IGIMFs, as compared to more standard IMFs.
[127]  [pdf] - 1017448
Polycyclic aromatic hydrocarbon selected galaxies
Comments: 14 pages, 11 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2009-09-09
Polycyclic aromatic hydrocarbon (PAH) emission has been found in both starbursts and modestly starforming galaxies, but the relation between starforming activity and PAH luminosity is still a matter of debate. The different correlation degrees could be caused by the variety of optical and far-infrared sample selection criteria. In order to obtain a census of the typical properties of PAH emitting galaxies, we here study moderately distant galaxies which have been selected by their PAH emission. Combining the ISOCAM Parallel Survey at 6.7 micron with 2MASS we have colour-selected a sample of 120 candidates for strong PAH emission. We obtained optical and mid-infrared spectra of 75 and 19 sources, respectively, and analysed IRAS-ADDSCANs and available Spitzer 3.6-160 micron photometry. The Spitzer spectra exhibit clear PAH features and corroborate that our photometric selection criteria trace the PAH emission of galaxies fairly well. The optical spectra show emission lines, at median redshift z~0.1, as well as Hdelta and CaII absorption, indicating ongoing starformation as well as post-starbursts. The mid- and far-infrared spectral energy distributions (SEDs) provide evidence for a broad range of far-infrared (FIR) luminosities, but in general the dust is colder than for starburst galaxies like M82. For most galaxies the monocromatic luminosity peaks at about equal height at optical, 6.7 micron (PAH) and FIR wavelengths. In about 15% of the sources the FIR luminosity exceeds the optical and PAH energy output by about a factor 5-10 despite the cool dust temperature; in these galaxies a large dust mass of 10^8 - 10^9 M_sun is inferred.
[128]  [pdf] - 24516
The Origin of Extended Disk Galaxies at z=2
Comments: 5 pages, 4 figures. Accepted for publication in MNRAS Letters. Minor changes to match published version
Submitted: 2009-05-21, last modified: 2009-07-20
Galaxy formation models typically assume that the size and rotation speed of galaxy disks are largely dictated by the mass, concentration, and spin of their surrounding dark matter haloes. Equally important, however, are the fraction of baryons in the halo that collect into the central galaxy, as well as the net angular momentum that they are able to retain during its assembly process. We explore the latter using a set of four large cosmological N-body/gasdynamical simulations drawn from the OWLS (OverWhelmingly Large Simulations) project. These runs differ only in their implementation of feedback from supernovae. We find that, when expressed as fractions of their virial values, galaxy mass and net angular momentum are tightly correlated. Galaxy mass fractions, m_d=M_gal/M_vir, depend strongly on feedback, but only weakly on halo mass or spin over the halo mass range explored here (M_vir>1e11 h^{-1}M_sun). The angular momentum of a galaxy, j_d=J_gal/J_vir, correlates with m_d in a manner that is insensitive to feedback and that deviates strongly from the simple j_d = m_d assumption often adopted in semi-analytic models of galaxy formation. The m_d-j_d correlation implies that, in a given halo, galaxy disk size is maximal when the central galaxy makes up a substantial fraction (~20%-30%) of all baryons within the virial radius. At z=2, such systems may host gaseous disks with radial scale lengths as large as those reported for star-forming disks by the SINS survey, even in moderately massive haloes of average spin. Extended disks at z=2 may thus signal the presence of systems where galaxy formation has been particularly efficient, rather than the existence of haloes with unusually high spin parameter.
[129]  [pdf] - 315452
Clustering of red galaxies around the z=1.53 quasar 3C270.1
Comments: 12 pages, 7 figures, accepted by ApJ
Submitted: 2009-01-05
In the paradigm of hierarchical galaxy formation, luminous radio galaxies mark mass assembly peaks that should contain clusters of galaxies. Observations of the z=1.53 quasar 3C270.1 with the Spitzer Space Telescope at 3.6-24 micron and with the 6.5-m MMT in the z'- and Y-bands allow detection of potential cluster members via photometric redshifts. Compared with nearby control fields, there is an excess of 11 extremely red objects (EROs) at 1.33 < z_phot < 1.73, consistent with a proto-cluster around the quasar. The spectral energy distributions (SEDs) of 3/4 of the EROs are better fitted with passive elliptical galaxies than withdust-reddened starbursts, and of four sources well-detected on an archival HST snapshot image, all have undisturbed morphologies. However, one ERO, not covered by the HST image, is a double source with 0.8" separation on the z' image and a marginal (2sigma) 24 micron detection indicating a dust-enshrouded starburst. The EROs are more luminous than L* (H = -23.6 AB mag at z=1.5).
[130]  [pdf] - 14843
Near- and mid-infrared photometry of high-redshift 3CR sources
Comments: 11 pages, 5 figures, accepted by ApJ
Submitted: 2008-07-24
Using the Spitzer Space Telescope, we have obtained 3.6--24 micron photometry of 38 radio galaxies and 24 quasars from the 3CR catalog at redshift 1<z<2.5. This 178 MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). Quasar and radio galaxy subsamples matched in isotropic radio luminosity are compared. The quasars all have similar spectral energy distributions (SEDs), nearly constant in nu F_nu through the rest 1.6-10 micron range, consistent with a centrally heated dust distribution which outshines the host galaxy contribution. The radio galaxy SEDs show larger dispersion, but the mean radio galaxy SED declines from rest 1.6 to 3 micron and then rises from 3 to 8 micron. The radio galaxies are on average a factor 3-10 less luminous in this spectral range than the quasars. These characteristics are consistent with composite emission from a heavily reddened AGN plus starlight from the host galaxy. The mid-infrared colors and radio to mid-infrared spectral slopes of individual galaxies are also consistent with this picture. Individual galaxies show different amounts of extinction and host galaxy starlight, consistent with the orientation-dependent unified scheme.
[131]  [pdf] - 518879
ACS imaging of star clusters in M51 II. The luminosity function and mass function across the disk
Comments: 15 pages, 9 figures, accepted for publication in Astronomy and Astrophysics, replaced equation 3
Submitted: 2008-06-23, last modified: 2008-07-23
Whether or not there exists a physical upper mass limit for star clusters is still unclear. HST/ACS data for the rich cluster population in the interacting galaxy M51 enables us to investigate this in more detail. We investigate whether the cluster luminosity function (LF) in M51 shows evidence for an upper limit to the mass function. The variations of the LF parameters with position on the disk are addressed. We determine the cluster LF for all clusters in M51 falling within our selection criteria, as well as for several subsets of the sample. In that way we can determine the properties of the cluster population as a function of galactocentric distance and background intensity. By comparing observed and simulated LFs we can constrain the underlying cluster initial mass function and/or cluster disruption parameters. A physical upper mass limit for star clusters will appear as a bend dividing two power law parts in the LF, if the cluster sample is large enough to sample the full range of cluster masses. The location of the bend in the LF is indicative of the value of the upper mass limit. The slopes of the power laws are an interplay between upper masses, disruption and fading. The LF of the cluster population of M51 is better described by a double power law than by a single power law. We show that the cluster initial mass function is likely to be truncated at the high mass end. We conclude from the variation of the LF parameters with galactocentric distance that both the upper mass limit and the cluster disruption parameters are likely to be a function of position in the galactic disk. At higher galactocentric distances the maximum mass is lower, cluster disruption slower, or both.
[132]  [pdf] - 14576
A near-infrared survey of the entire R Corona Australis cloud
Comments: 11 pages, 12 Figures, accepted by Astronomy and Astrophysics
Submitted: 2008-07-16
To understand low- to intermediate-mass star-formation in the nearby R CrA molecular cloud, we try to identify the stellar content that is accessible with near-infrared observations. We obtained a JHK band mosaic of 10 x 60 arcmin square covering the entire RCrA molecular cloud with unprecedented sensitivity. We present a catalogue of about 3500 near-infrared sources fainter than the saturation limit K = 10 mag, reaching K = 18mag. We analysed the extended sources by inspecting their morphology and point sources by means of colour-colour and colour-magnitude diagrams. Additionally, we compared the extinction inferred from the NIR data with the line-of-sight dust emission at 1.2 mm. Sources towards high dust emission but relatively low H-K show a projected mm-exces; these sources are either immediately surrounded by cold circumstellar material or, if too red to be a true foreground object, they are embedded in the front layer of the 1.2 mm emitting dust cloud. In both cases they are most likely associated with the cloud.
[133]  [pdf] - 12828
Nuclear activity in nearby galaxies: Mid-infrared imaging with the VLT
Comments: accepted by A&A: 14 pages, 5 Figures;
Submitted: 2008-05-21
Dust enshrouded activity can ideally be studied by mid-infrared (MIR) observations. In order to explore the AGN versus star forming origin of the nuclear MIR emission of galaxies, observations of high spatial resolution are required. Here we report on 11.3mic. observations with VISIR at the VLT, reaching 0.35'' spatial resolution (FWHM). During the scientific verification of VISIR we have observed a sample of 36 nearby galaxies having a variety of optically classified nuclear activity: 17 black hole driven active galactic nuclei (AGN), 10 starbursts (SBs) and 9 quiet spirals. 16/17 AGN are detected and unresolved, 5/10 SBs are detected and resolved with structured emission up to a few arcsec, while for 5/10 SB and all 9 quiet nuclei low upper limits are provided. The morphology of the resolved SB nuclei follows that seen at radio frequencies. The compactness of AGN and the extent of the SB nuclei is consistent with predictions from radiative transfer models and with MIR spectra of lower spatial resolution. We explore the nuclear MIR surface brightness as a quantitative measure. While AGN and SB cannot be distinguished with MIR data from 4m class telescopes, our data provide evidence that, up to a distance of 100 Mpc, AGN and SB can well be separated by means of MIR surface brightness when using 8m class telescopes.
[134]  [pdf] - 12635
Mid-Infrared Spectra of High Redshift (z>2) Radio Galaxies
Comments: 6 pages, 3 figures, includes emulateapj.cls, in press ApJL
Submitted: 2008-05-14
We present the first mid-infrared Spitzer/Infrared Spectrograph (IRS) observations of powerful radio galaxies at z>2. These radio galaxies, 4C +23.56 (z=2.48) and 6C J1908+7220 (z=3.53), both show strong mid-infrared continua, but with 6C J1908+7220 also showing strong PAH emission at rest-frame 6.2 and 7.7um. In 4C+23.56 we see no obvious PAH features above the continuum. The PAH emission in 6C J1908+7220 is the amongst the most distant observed to date and implies that there is a large instantaneous star formation rate (SFR). This is consistent with the strong detection of 6C J1908+7220 at far-IR and sub-mm wavelengths, indicative of large amounts of cold dust, ~10^9Msun. Powerful radio galaxies at lower redshifts tend to have weak or undetectable PAH features and typically have lower far-IR luminosities. In addition, 4C 23.56 shows moderate silicate absorption as seen in less luminous radio galaxies, indicating tau_{9.7um}=0.3+/-0.05. This feature is shifted out of the observed wavelength range for 6C J1908+7220. The correlation of strong PAH features with large amounts of cold dust, despite the presence of a powerful AGN, is in agreement with other recent results and implies that star formation at high redshift is, in some cases at least, associated with powerful, obscured AGN.
[135]  [pdf] - 11507
Spitzer Observations of M33 and the Hot Star, H II Region Connection
Comments: Accepted MNRAS: 47 pages, 12 figures, 6 tables, abstract truncated
Submitted: 2008-04-04
We have observed emission lines of [S IV] 10.51, H(7-6) 12.37, [Ne II] 12.81, [Ne III] 15.56, and [S III] 18.71 um in a number of extragalactic H II regions with the Spitzer Space Telescope. A previous paper presented our data and analysis for the substantially face-on spiral galaxy M83. Here we report our results for the local group spiral galaxy M33. The nebulae selected cover a wide range of galactocentric radii (R_G). The observations were made with the Infrared Spectrograph with the short wavelength, high resolution module. The above set of five lines is observed cospatially, thus permitting a reliable comparison of the fluxes. From the measured fluxes, we determine the ionic abundance ratios including Ne++/Ne+, S3+/S++, and S++/Ne+ and find that there is a correlation of increasingly higher ionization with larger R_G. By sampling the dominant ionization states of Ne (Ne+, Ne++) and S (S++, S3+) for H II regions, we can estimate the Ne/H, S/H, and Ne/S ratios. We find from linear least-squares fits that there is a decrease in metallicity with increasing R_G: d log (Ne/H)/dR_G = -0.058+-0.014 and d log (S/H)/dR_G = -0.052+-0.021 dex kpc-1. There is no apparent variation in the Ne/S ratio with R_G. Unlike our previous similar study of M83, where we conjectured that this ratio was an upper limit, for M33 the derived ratios are likely a robust indication of Ne/S. This occurs because the H II regions have lower metallicity and higher ionization than those in M83. Both Ne and S are primary elements produced in alpha-chain reactions, following C and O burning in stars, making their yields depend very little on the stellar metallicity. Thus, it is expected that Ne/S remains relatively constant throughout a galaxy. The median (average) Ne/S ratio derived for H II regions in M33 is 16.3 (16.9), just slightly higher than
[136]  [pdf] - 4217
The reddest ISO-2MASS quasar
Comments: 8 pages, 9 figures
Submitted: 2007-08-22
In the course of the NIR/MIR AGN search combining the 6.7 mu ISOCAM Parallel Survey and 2MASS we have discovered 24 type-1 quasars about a third of which are too red to be discriminated by optical/UV search techniques. Here we report on a detailed case study of the reddest type-1 quasar of our sample (J2341) at redshift z=0.236 with M_K=-25.8 and J-K=1.95. We performed spectroscopy in the optical with VLT/FORS1 and in the MIR with Spitzer as well as NIR imaging with ISPI at CTIO. The optical and NIR observations reveal a star forming emission-line galaxy at the same redshift as the quasar with a projected linear separation of 1.8 arcsec (6.7 kpc). The quasar and its companion are embedded in diffuse extended continuum emission. Compared with its companion the quasar exhibits redder optical-NIR colours, which we attribute to hot nuclear dust. The MIR spectrum shows only few emission lines superimposed on a power-law spectral energy distribution. However, the lack of strong FIR emission suggests that our potentially interacting object contains much less gas and dust and is in a stage different from dust reddened ULIRG-AGN like Mrk 231. The optical spectrum shows signatures for reddening in the emission-lines and no post-starburst stellar population is detected in the host galaxy of the quasar. The optical continuum emission of the active nucleus appears absorbed and diluted. Even the combination of absorption and host dilution is not able to match J2341 with standard quasar templates. While the BLR shows only a rather moderate absorption of E_(B-V)=0.3, the continuum shorter than 4500 AA requires strong obscuration with E_(B-V)=0.7, exceeding the constraints from the low upper limit on the 9.7 mu silicate absorption. This leads us to conclude that the continuum of J2341 is intrinsically redder than that of typical quasars.
[137]  [pdf] - 3214
VISIR/VLT mid-infrared imaging of Seyfert nuclei: Nuclear dust emission and the Seyfert-2 dichotomy
Comments: 11 pages, 6 figures, accepted by A&A
Submitted: 2007-07-18
Half of the Seyfert-2 galaxies escaped detection of broad lines in their polarised spectra observed so far. Some authors have suspected that these non-HBLRs contain real Sy2 nuclei without intrinsic broad line region hidden behind a dust torus. If this were true, then their nuclear structure would fundamentally differ from that of Sy2s with polarised broad lines: in particular, they would not be explained by orientation-based AGN unification. Further arguments for two physically different Sy2 populations have been derived from the warm and cool IRAS F25/F60 ratios. These ratios, however, refer to the entire host galaxies and are unsuitable to conclusively establish the absence of a nuclear dust torus. Instead, a study of the Seyfert-2 dichotomy should be performed on the basis of nuclear properties only. Here we present the first comparison between [OIII] 5007A and mid-infrared imaging at matching spatial resolution. Exploring the Seyfert-2 dichotomy we find that the distributions of nuclear mid-infrared/[OIII] luminosity ratios are indistinguishable for Sy1s and Sy2s with and without detected polarised broad lines and irrespective of having warm or cool IRAS F25/F60 ratios. We find no evidence for the existence of a population of real Sy2s with a deficit of nuclear dust emission. Our results suggest 1) that all Seyfert nuclei possess the same physical structure including the putative dust torus and 2) that the cool IRAS colours are caused by a low contrast of AGN to host galaxy. Then the Seyfert-2 dichotomy is explained in part by unification of non-HBLRs with narrow-line Sy1s and to a larger rate by observational biases caused by a low AGN/host contrast and/or an unfavourable scattering geometry.
[138]  [pdf] - 811
ACS imaging of star clusters in M51. I. Identification and radius distribution
Comments: 17 pages, 23 figures. Accepted for publication in A&A
Submitted: 2007-04-26
We use HST/ACS observations of the spiral galaxy M51 in F435W, F555W and F814W to select a large sample of star clusters with accurate effective radius measurements in an area covering the complete disc of M51. We present the dataset and study the radius distribution and relations between radius, colour, arm/interarm region, galactocentric distance, mass and age. We select a sample of 7698 (F435W), 6846 (F555W) and 5024 (F814W) slightly resolved clusters and derive their effective radii by fitting the spatial profiles with analytical models convolved with the point spread function. The radii of 1284 clusters are studied in detail. We find cluster radii between 0.5 and ~10 pc, and one exceptionally large cluster candidate with a radius of 21.6 pc. The median radius is 2.1 pc. We find 70 clusters in our sample which have colours consistent with being old GC candidates and we find 6 new "faint fuzzy" clusters in, or projected onto, the disc of M51. The radius distribution can not be fitted with a power law, but a log-normal distribution provides a reasonable fit to the data. This indicates that shortly after the formation of the clusters from a fractal gas, their radii have changed in a non-uniform way. We find an increase in radius with colour as well as a higher fraction of redder clusters in the interarm regions, suggesting that clusters in spiral arms are more compact. We find a correlation between radius and galactocentric distance which is considerably weaker than the observed correlation for old Milky Way GCs. We find weak relations between cluster luminosity and radius, but we do not observe a correlation between cluster mass and radius.
[139]  [pdf] - 90026
Spitzer Observations of M83 and the Hot Star, H II Region Connection
Comments: Accepted MNRAS: 31 pages, 7 figures, 3 tables
Submitted: 2007-03-09
We have undertaken a program to observe emission lines of SIV 10.5, NeII 12.8, NeIII 15.6, & SIII 18.7 um in a number of extragalactic HII regions with the Spitzer Space Telescope. We report our results for the nearly face-on spiral galaxy M83. The nebulae selected cover a wide range of galactocentric radii (R_G). The observations were made with the Infrared Spectrograph in the short wavelength, high dispersion configuration. The above set of 4 lines is observed cospatially, thus permitting a reliable comparison of the fluxes. From the measured fluxes, we determine the ionic abundance ratios including Ne++/Ne+, S3+/S++, and S++/Ne+ and find that there is a correlation of increasingly higher ionization with larger R_G. By sampling the dominant ionization states of Ne and S for HII regions, Ne/S ~ (Ne+ + Ne++)/(S++ + S3+). Our findings of ratios that exceed the benchmark Orion value are more likely due to other effects than a true gradient in Ne/S. Both Ne and S are primary elements produced in alpha- chain reactions. It is expected that Ne/S remains relatively constant throughout a galaxy. This type of observation and method of analysis has the potential for accurate measurements of Ne/S, particularly for HII regions with lower metallicity & higher ionization than those here, such as those in M33. Our observations may also be used to test the predicted ionizing spectral energy distribution of various stellar atmosphere models. We compare the fractional ionization ratios <Ne++>/<S++> & <Ne++>/<S3+> vs. <S3+>/<S++> with predictions made from our photoionization models using several state-of-the-art stellar atmosphere model grids. A second paper using Spitzer observations of HII regions in the galaxy M33 will follow.
[140]  [pdf] - 88969
Narrow-line AGN in the ISO-2MASS Survey
Comments: Accepted for publication in A&A
Submitted: 2007-01-29
Context: A long-standing challenge of observational AGN research is to find type 2 quasars, the luminous analogues of Seyfert-2 galaxies. Aims: We search for luminous narrow-line type 2 AGN, characterise their properties, and compare them with broad-line type 1 AGN. Methods: Combining the ISOCAM parallel survey at 6.7 mu with 2MASS, we have selected AGN via near-mid-infrared colours caused by the hot nuclear dust emission. We performed spectroscopy in the optical and, for a subset of the sample, also in the mid-infrared with Spitzer. Results: We find nine type 2 AGN at redshift 0.1<z<0.5, three of them have even quasar-like [OIII] luminosities. At the given redshift and luminosity range the number of type 2 AGN is at least as high as that of type 1s. At z>0.5 we did not find type 2 AGN, probably because the hottest dust emission, still covered by the NIR filters, is obscured. The optical spectra of the type 2 host galaxies show young and old stellar populations. Only one object is an ultraluminous infrared galaxy with starburst. The 5-38 mu spectra of the two type 2 sources observed show a strong continuum with PAH emission in one case and silicate absorption in the other case. Conclusions: The near-mid-infrared selection is a successful strategy to find luminous type 2 AGN at low z. The objects exhibit a large range of properties so that it is difficult to infer details by means of popular SED fitting with simple average templates.
[141]  [pdf] - 88008
Unification of 3CR Radio Galaxies and Quasars
Comments: 2 pages, Conference proceedings to appear in "The Central Engine of Active Galactic Nuclei", ed. L. C. Ho and J.-M. Wang (San Francisco: ASP)
Submitted: 2006-12-21
We have observed seven powerful FR2 radiogalaxies and seven quasars with the Spitzer IRS. Both samples are comparable in both, redshift range and isotropic 178 Hz luminosity. Both samples are found to have similar distributions in the luminosity ratios of Mid-IR high- and low-excitation lines ([NeV]/[NeII]), and of Mid-IR high-excitation line to radio power ratio ([NeV]/P_178MHz). However, the MIR/FIR ratio is generally higher for quasars. We further observed Silicate features at 10 and 18micron in emission. In our sample only quasars show emission features, while silicate absorption is seen only in the radio galaxies. These observations are all in agreement with unification schemes that explain both groups as the same class of objects seen under different orientation angles.
[142]  [pdf] - 87243
Properties of the Narrow-Line Region in Seyfert Galaxies
Comments: 4 pages, 4 figures, Conference proceedings to appear in "The Central Engine of Active Galactic Nuclei", ed. L. C. Ho and J.-M. Wang (San Francisco: ASP)
Submitted: 2006-11-28
We study the narrow-line region (NLR) of six Seyfert-1 and six Seyfert-2 galaxies by means of spatially resolved optical spectroscopy and photoionization modelling. From spatially resolved spectral diagnostics, we find a transition between the AGN-excited NLR and the surrounding star-forming regions, allowing us to determine the NLR size independent of stellar contamination. CLOUDY photoionization models show that the observed transition represents a true difference in ionization source and cannot be explained by variations of physical parameters. The electron density and ionization parameter decrease with radius indicating that the NLR is photoionized by the central source only. The velocity field suggests a disky NLR gas distribution.
[143]  [pdf] - 84610
Variation of the cluster luminosity function across the disk of M51
Comments: 2 pages, no figures. To appear in "Mass loss from stars and the evolution of stellar clusters". Proc. of a workshop held in honor of Henny Lamers, Lunteren, The Netherlands, May 29 - June 1, 2006. Editors Alex de Koter, Linda Smith and Rens Waters (San Francisco: ASP)
Submitted: 2006-09-01
We study the luminosity function (LF) of the star clusters in M51. Comparing the observed LF with the LF resulting from artificial cluster populations suggests that there exists an upper mass limit for clusters and that this limit and/or the cluster disruption varies with galactocentric distance.
[144]  [pdf] - 83839
Size and properties of the narrow-line region in Seyfert-1 galaxies from spatially-resolved optical spectroscopy
Comments: 25 pages, 9 figures, accepted for publication in A&A
Submitted: 2006-07-27
Spatially resolved emission-line spectroscopy is a powerful tool to determine the physical conditions in the narrow-line region (NLR) of active galactic nuclei (AGNs). We recently used optical long-slit spectroscopy to study the NLRs of a sample of six Seyfert-2 galaxies. We have shown that such an approach, in comparison to the commonly used [OIII] narrow-band imaging alone, allows us to probe the size of the NLR in terms of AGN photoionisation. Moreover, several physical parameters of the NLR can be directly accessed. We here apply the same methods to study the NLR of six Seyfert-1 galaxies and compare our results to those of Seyfert-2 galaxies. We employ diagnostically valuable emission-line ratios to determine the physical properties of the NLR, including the core values and radial dependencies of density, ionisation parameter, and reddening. Tracking the radial change of emission-line ratios in diagnostic diagrams allows us to measure the transition between AGN-like and HII-like line excitation, and thus we are able to measure the size of the NLR. In the diagnostic diagrams, we find a transition between line ratios falling in the AGN regime and those typical for HII regions in two Seyfert-1 galaxies, thus determining the size of the NLR. The central electron temperature and ionisation parameter are in general higher in type-1 Seyferts than in type 2s. In almost all cases, both electron density and ionisation parameter decrease with radius. In general, the decrease is faster in Seyfert-1 galaxies than in type 2s. In several objects, the gaseous velocity distribution is characteristic for rotational motion in an (inclined) emission-line disk in the centre. We give estimates of the black hole masses. We discuss our findings in detail for each object.
[145]  [pdf] - 82807
Size and properties of the narrow-line region in Seyfert-2 galaxies from spatially-resolved optical spectroscopy
Comments: 23 pages, 41 figures, accepted for publication in A&A
Submitted: 2006-06-14
While [OIII] narrow-band imaging is commonly used to measure the size of the narrow-line regions (NLRs) in active galactic nuclei (AGNs), it can be contaminated by emission from surrounding starbursts. Recently, we have shown that long-slit spectroscopy provides a valuable alternative approach to probe the size in terms of AGN photoionisation. Moreover, several parameters of the NLR can be directly accessed. We here apply the same methods developed and described for the Seyfert-2 galaxy NGC1386 to study the NLR of five other Seyfert-2 galaxies by using high-sensitivity spatially-resolved optical spectroscopy obtained at the VLT and the NTT. We probe the AGN-photoionisation of the NLR and thus, its ``real'' size using diagnostic line-ratio diagrams.We derive physical properties of the NLR such as reddening, ionisation parameter, electron density, and velocity as a function of distance from the nucleus. For NGC5643, the diagnostic diagrams unveil a similar transition between line ratios falling in the AGN regime and those typical for HII regions as found for NGC1386, thus determining the size of the NLR. For the other four objects, all measured line ratios fall in the AGN regime. In almost all cases, both electron density and ionisation parameter decrease with radius. Deviations from this general behaviour (such as a secondary peak) seen in both the ionisation parameter and electron density can be interpreted as signs of shocks from the interaction of a radio jet and the NLR gas. In several objects, the gaseous velocity distribution is characteristic for rotational motion in an (inclined) emission-line disk in the centre. We compare our results to those of NGC1386 and show that the latter can be considered as prototypical also for this larger sample. We discuss our findings in detail for each object.
[146]  [pdf] - 1233781
The radii of thousands of star clusters in M51 with HST/ACS
Comments: 4 pages, 2 figures, contribution to "Globular Clusters: Guides to Galaxies", March 6th-10th, 2006
Submitted: 2006-05-01
We exploit the superb resolution of the new HST/ACS mosaic image of M51 to select a large sample of young (< 1 Gyr) star clusters in the spiral disk, based on their sizes. The image covers the entire spiral disk in B, V, I and H_alpha, at a resolution of 2 pc per pixel. The surface density distribution of 4357 resolved clusters shows that the clusters are more correlated with clouds than with stars, and we find a hint of enhanced cluster formation at the corotation radius. The radius distribution of a sample of 769 clusters with more accurate radii suggests that young star clusters have a preferred effective radius of ~3 pc, which is similar to the preferred radius of the much older GCs. However, in contrast to the GCs, the young clusters in M51 do not show a relation between radius and galactocentric distance. This means that the clusters did not form in tidal equilibrium with their host galaxy, nor that their radius is related to the ambient pressure.
[147]  [pdf] - 1233779
The Maximum Mass of Star Clusters
Comments: 4 pages, contribution to "Globular Clusters: Guides to Galaxies", March 6th-10th, 2006
Submitted: 2006-04-27
When an universal untruncated star cluster initial mass function (CIMF) described by a power-law distribution is assumed, the mass of the most massive star cluster in a galaxy (M_max) is the result of the size-of-sample (SoS) effect. This implies a dependence of M_max on the total number of star clusters (N). The SoS effect also implies that M_max within a cluster population increases with equal logarithmic intervals of age. This is because the number of clusters formed in logarithmic age intervals increases (assuming a constant cluster formation rate). This effect has been observed in the SMC and LMC. Based on the maximum pressure (P_int) inside molecular clouds, it has been suggested that a physical maximum mass (M_max[phys]) should exist. The theory predicts that M_max[phys] should be observable, i.e. lower than M_max that follows from statistical arguments, in big galaxies with a high star formation rate. We compare the SoS relations in the SMC and LMC with the ones in M51 and model the integrated cluster luminosity function (CLF) for two cases: 1) M_max is determined by the SoS effect and 2) M_max=M_max[phys]=constant. The observed CLF of M51 and the comparison of the SoS relations with the SMC and LMC both suggest that there exists a M_max[phys] of 5*10^5 M_sun in M51. The CLF of M51 looks very similar to the one observed in the ``Antennae'' galaxies. A direct comparison with our model suggests that there M_max[phys]=2*10^6 M_sun.
[148]  [pdf] - 78498
Observational evidence for a truncation of the star cluster initial mass function at the high mass end
Comments: 4 pages, 2 figures. A&A Letters, accepted
Submitted: 2005-12-12
We present the luminosity function (LF) of star clusters in M51 based on HST/ACS observations taken as part of the Hubble Heritage project. The clusters are selected based on their size and with the resulting 5990 clusters we present one of the largest cluster samples of a single galaxy. We find that the LF can be approximated with a double power-law distribution with a break around M_V = -8.9. On the bright side the index of the power-law distribution is steeper (a = 2.75) than on the faint-side (a = 1.93), similar to what was found earlier for the ``Antennae'' galaxies. The location of the bend, however, occurs about 1.6 mag fainter in M51. We confront the observed LF with the model for the evolution of integrated properties of cluster populations of Gieles et al., which predicts that a truncated cluster initial mass function would result in a bend in, and a double power-law behaviour of, the integrated LF. The combination of the large field-of view and the high star cluster formation rate of M51 make it possible to detect such a bend in the LF. Hence, we conclude that there exists a fundamental upper limit to the mass of star clusters in M51. Assuming a power-law cluster initial mass function with exponentional cut-off of the form NdM ~ M^-b * exp(-M/M_C)dM, we find that M_C = 10^5 M_sun. A direct comparison with the LF of the ``Antennae'' suggests that there M_C = 4*10^5 M_sun.
[149]  [pdf] - 77637
Size and Properties of the NLR in the Seyfert-2 Galaxy NGC 1386
Comments: 14 pages, 20 figures, accepted for publication in A&A
Submitted: 2005-11-09
We study the narrow-line region (NLR) of the Seyfert-2 galaxy NGC 1386 by means of long-slit spectroscopy obtained with FORS1 at the VLT. We use the galaxy itself for subtracting the stellar template, applying reddening corrections to fit the stellar template to the spectra of the NLR. The continuum gets steadily redder towards the nucleus. The spatial distribution of the reddening derived from the Balmer decrement differs from the continuum reddening, indicating dust within the NLR with a varying column density along the line of sight. Using spatially resolved spectral diagnostics, we find a transition between central line ratios falling into the AGN regime and outer ones in the HII-region regime, occuring at a radius of ~6 arcsec (310 pc) in all three diagnostic diagrams. Applying CLOUDY photoionisation models, we show that the observed distinction between HII-like and AGN-like ratios in NGC 1386 represents a true difference in ionisation source and cannot be explained by variations of physical parameters such as ionisation parameter, electron density or metallicity. We interpret it as a real border between the NLR, i.e. the central AGN-photoionised region and surrounding HII regions. We derive surface brightness, electron density, and ionisation parameter as a function of distance from the nucleus. Both the electron density and the ionisation parameter decrease with radius. We discuss the consequences of these observations for the interpretation of the empirical NLR size-luminosity relation. In the outer part of the NLR, we find evidence for shocks, resulting in a secondary peak of the electron-density and ionisation-parameter distribution north of the nucleus. We compare the NLR velocity curve with the stellar one and discuss the differences.
[150]  [pdf] - 75922
Spitzer IRS spectroscopy of 3CR radio galaxies and quasars: Testing the unified schemes
Comments: 5 pages, 3 figures, accepted as Letter for A&A
Submitted: 2005-09-13
With the Spitzer Space Telescope Infrared Spectrograph (IRS) we have observed seven powerful FR2 radio galaxies and seven quasars. Both samples, the galaxies and the quasars, are comparable in isotropic 178 MHz luminosity (10^26.5 W/Hz < P 178MHz < 10^29.5 W/Hz) and in redshift range (0.05 < z < 1.5). We find for both samples similar distributions in the luminosity ratios of high- to low-excitation lines ([NeV] 14.3 mu / [NeII] 12.8 mu) and of high-excitation line to radio power ([NeV] 24.3 \mu / P 178MHz). This solves the long debate about the apparent difference of quasars and radio galaxies in favor of the orientation-dependent unified schemes. Furthermore, the luminosity ratio [OIII] 500.7 nm / [OIV] 25.9 mu of most galaxies is by a factor of ten lower than that of the quasars. This suggests that the optical emission from the central NLR is essentially absorbed (A_V} > 3) in the powerful FR2 galaxies and that the [OIII] 500.7 nm luminosity does not serve as isotropic tracer for testing the unified schemes.
[151]  [pdf] - 75666
MAMBO 1.25mm observations of 3CR radio galaxies and quasars at z=1.5: On the debate of the unified schemes
Comments: 8 pages, 4 Figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2005-09-05
In order to explore the nature of the 850 micron flux difference between powerful radio galaxies and steep radio-spectrum quasars at z=1.5 reported by Willott et al. (2002), we have observed 9 sources from their sample of 11 quasars at 1.25 mm. For 7 sources the 1.25 mm fluxes are much brighter than one would expect from a purely thermal dust model fitted to the submm data, providing evidence for the synchrotron nature of the observed 1.25 mm radiation. If we extrapolate a power-law synchrotron spectrum to shorter wavelengths, then for 6 of the 9 sources also the 850 micron fluxes are dominated by synchrotron radiation. We discuss how far the (sub)-millimetre data can be interpreted in accordance with the orientation-dependent unified schemes for powerful radio galaxies and quasars. In this case the results challenge the reported evidence for the receding torus model and for the evolutionary trend of a declining dust luminosity with increasing projected size of the radio lobes.
[152]  [pdf] - 74489
The ISO-2MASS AGN survey: On the type-1 sources
Comments: 5 pages, 5 figures, accepted for publication as A&A Letter
Submitted: 2005-07-13
We combined the ISOCAM Parallel Mode Survey at 6.7 micron (LW2 filter) with the Two Micron All Sky Survey in order to obtain a powerful tool to search for AGN independent of dust extinction. Using moderate colour criteria H-K>0.5 and K-LW2>2.7 we have selected a sample of 77 AGN candidates in an effective area of about 10 square degrees. By means of optical spectroscopy we find 24 (30%) type-1 QSOs at redshifts 0.1<z<2.3; nine of them have z>0.8. About one third of the ISO-2MASS QSOs show so red optical colours, that they are missed in optical and UV AGN surveys like SDSS, 2DF, or HES. With a surface density of about 2 deg^(-2) down to R<18 mag the ISO-2MASS QSOs outnumber the 1.35 deg^(-2) of the SDSS quasar survey by 50%; we find a combined optical-IR QSO surface density of 2.7 deg^(-2). Since only two of the ISO-2MASS QSOs have also J-K>2, the inclusion of the ISO mid-infrared photometry significantly extends the capabilities of the pure 2MASS red AGN survey. We suggest that the newly found red AGN resemble young members of the quasar population, and that quasars spend much of their lifetime in a dust enshrouded phase.
[153]  [pdf] - 74334
Obscured Activity: AGN, Quasars, Starbursts and ULIGs observed by the Infrared Space Observatory
Comments: Accepted for publication in 'ISO science legacy - a compact review of ISO major achievements', Space Science Reviews - dedicated ISO issue. To be published by Springer in 2005. 62 pages (low resolution figures version). Higher resolution PDFs available from or
Submitted: 2005-07-07
Some of the most active galaxies in the Universe are obscured by large quantities of dust and emit a substantial fraction of their bolometric luminosity in the infrared. Observations of these infrared luminous galaxies with the Infrared Space Observatory (ISO) have provided a relatively unabsorbed view to the sources fuelling this active emission. The improved sensitivity, spatial resolution and spectroscopic capability of ISO over its predecessor Infrared Astronomical Satellite (IRAS), has enabled significant advances in the understanding of the infrared properties of active galaxies. ISO surveyed a wide range of active galaxies which, in the context of this review, includes those powered by intense bursts of star-formation as well as those containing a dominant active galactic nucleus (AGN). Mid infrared imaging resolved for the first time the dust enshrouded nuclei in many nearby galaxies, while a new era in infrared spectroscopy was opened by probing a wealth of atomic, ionic and molecular lines as well as broad band features in the mid and far infrared. This was particularly useful since it resulted in the understanding of the power production, excitation and fuelling mechanisms in the nuclei of active galaxies including the intriguing but so far elusive ultraluminous infrared galaxies. Detailed studies of various classes of AGN and quasars greatly improved our understanding of the unification scenario. Far-infrared imaging and photometry also revealed the presence of a new very cold dust component in galaxies and furthered our knowledge of the far-infrared properties of faint starbursts, ULIGs and quasars. We summarise almost nine years of key results based upon ISO data spanning the full range of luminosity and type of active galaxies.
[154]  [pdf] - 72360
Discovery of 10 micron silicate emission in quasars. Evidence of the AGN unification scheme
Comments: accepted A&A Letters: 2 Fig., 4 pages
Submitted: 2005-04-12
According to the unified scheme, AGN are surrounded by a dust-torus, and the observed diversity of AGN properties results from the different orientations relative to our line of sight. The strong resonance of silicate dust at 10 micron is therefore, as expected, seen in absorption towards many type-2 AGN. In type-1 AGN, it should be seen in emission because the hot inner surface of the dust torus becomes visible. However, this has not been observed so far, thus challenging the unification scheme or leading to exotic modifications of the dust-torus model. Here we report the discovery of the 10 micron silicate feature in emission in two luminous quasars with the Infrared Spectrograph of the Spitzer Space Telescope.
[155]  [pdf] - 71324
VPMS J1342+2840 - an unusual quasar from the variability and proper motion survey
Comments: 4 pages, 3 figures, accepted for publication in A&A Letters
Submitted: 2005-02-25
We report the discovery of the highly peculiar, radio-loud quasar VPMS J1342+2840 (z ~ 1.3) from the variability and proper motion survey. We present spectroscopic, imaging and photometric observations. The unusual spectrum shows a strong depression of the continuum over a wide wavelength range in the blue part without the typical structures of broad absorption line (BAL) troughs. The image of the quasar is unresolved and there is no evidence for a foreground object on the line of sight. The broad-band spectral energy distribution is not consistent with obvious dust reddening with the standard SMC extinction curve. The downturn of the continuum flux of VPMS J1342+2840 at short wavelengths can be caused by dust reddening only if the reddening curve is steeper then the SMC curve in the ultraviolet and is very flat at longer wavelengths. Alternatively, the dominant spectral features can be explained by low-ionization BALs forming unusually wide, overlapping absorption troughs.
[156]  [pdf] - 69878
The dust SED in the dwarf galaxy NGC 1569: Indications for an altered dust composition?
Comments: 4 pages, conference proceedings paper, "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu & R.J. Tuffs, AIP Conf. Ser., in press
Submitted: 2004-12-17
We discuss the interpretation of the dust SED from the mid-infrared to the millimeter range of NGC 1569. The model developed by D\'esert et al. (1990) including three dust components (Polyaromatic Hydrocarbons, Very Small Grains and big grains) can explain the data using a realistic interstellar radiation field and adopting an enhanced abundance of VSGs. A simple three-temperature model is also able to reproduce the data but requires a very low dust temperature which is considered to be unlikely in this low-metallicity starburst galaxy. The high abundance of Very Small Grains might be due to large grain destruction in supernova shocks. This possibility is supported by ISO data showing that the emission at 14.3 $\mu$m, tracing VSGs, is enhanced with respect to the emission at 6.7 $\mu$m and 850 $\mu$m in regions of high star formation.
[157]  [pdf] - 67800
Dust in 3CR radio galaxies: On the FR1-FR2 difference
Comments: 4 pages 3 figures Letter accepted by A&A
Submitted: 2004-09-29
We compare three 3CR samples of 11 FR1 galaxies, 17 FR2 galaxies and 18 lobe-dominated quasars contained in the ISO Data Archive. In contrast to the powerful FR2 galaxies with edge-brightened lobes, the low radio power FR1 galaxies in our sample do not exhibit any high MIR or FIR dust luminosity, which is typical for a buried, intrinsically more luminous AGN. This consolidates the fact already inferred from optical studies that their AGNs have only a relatively low luminosity. Also the FR1 galaxies show a high FIR/MIR luminosity ratio, compared to quasars, suggesting that their FIR luminosity is substantially powered by the interstellar radiation field (ISRF) of the giant elliptical hosts. Finally, we discuss the FR1-FR2 morphological dichotomy. FR1 galaxies do not have more interstellar matter (ISM) than FR2s as traced -- on the large scale -- by the cool FIR emitting dust and -- in the nuclear region -- by the warm MIR emitting dust. Due to the lack of central gas we suggest that the black holes of our FR1 galaxies are fed at a lower accretion rate than those of the FR2 galaxies.
[158]  [pdf] - 65268
The ISOPHOT-MAMBO survey of 3CR radio sources: Further evidence for the unified schemes
Comments: 16 pages, 7 figures, 3 tables, accepted by Astronomy & Astrophysics
Submitted: 2004-06-04
We present the complete set of ISOPHOT observations of 3CR radio galaxies and quasars, which are contained in the ISO Data Archive, providing 75 mid- and far-infrared spectral energy distributions (SEDs) between 5 and 200 micron. For 28 sources they are supplemented with MAMBO 1.2 mm observations and for 15 sources with new submillimetre data from the SCUBA archive. We check the orientation-dependent unified scheme, in which the powerful FR2 narrow line galaxies are quasars viewed at high inclination, so that their nuclei are hidden behind a dust torus intercepting the optical-ultraviolet AGN radiation and reemitting it in the infrared. We find that (1) both the quasars and the galaxies show a high mid- to far-infrared luminosity ratio typical for powerful AGNs and (2) -- when matched in 178 MHz luminosity -- both show the same ratio of isotropic far-infrared to isotropic 178 MHz lobe power. Therefore, from our large sample investigated here we find strong evidence for the orientation-dependent unification of the powerful FR2 galaxies with the quasars.
[159]  [pdf] - 65168
Observations of Cold Dust in Nearby Elliptical Galaxies
Comments: 22 pages, 9 figures, tentatively scheduled for the ApJ 10 September 2004, v612 issue
Submitted: 2004-06-01
Spectral energy distribution (SED) analyses, that include new millimeter (mm) to far-infrared (FIR) observations, obtained with continuum instruments on the Nobeyama and James Clerk Maxwell Telescopes and Infrared Space Observatory (ISO), are presented for seven, nearby (< 45Mpc) FIR-bright elliptical galaxies. These are analyzed together with archival FIR and short-wave radio data obtained from the NASA/IPAC Extragalactic Database (NED). The radio to infrared SEDs are best fit by power-law plus greybody models of dust residing in the central galactic regions within 2.4kpc diameter and with temperatures between about 21 and 28K, emissivity index approximately 2, and masses from about 0.16 to 1.9 million solar masses. The emissivity index is consistent with dust constituting amorphous silicate and carbonaceous grains previously modeled for stellar-heated dust observed in the Galaxy and other nearby extragalactic sources. Using updated dust absorption coefficients for this type of dust, dust masses are estimated that are similar to those determined from earlier FIR data alone, even though the latter results implied hotter dust temperatures. Fluxes and masses that are consistent with the new FIR and submm data are estimated for dust cooler than 20K within the central galactic regions. Tighter physical constraints for such cold, diffuse dust (if it exists) with low surface brightness, will need sensitive FIR-submm observations with the Spitzer Space Telescope, SCUBA2, or ALMA.
[160]  [pdf] - 64312
On the Measurement of Elemental Abundance Ratios in Inner Galaxy H II Regions
Comments: 29 pages, 6 figures, ApJ in press
Submitted: 2004-04-20
Although variations in elemental abundance ratios in the Milky Way Galaxy certainly exist, details remain uncertain, particularly in the inner Galaxy, where stars and H II regions in the Galactic plane are obscured optically. In this paper we revisit two previously studied, inner Galaxy H II regions: G333.6-0.2 and W43. We observed three new positions in G333.6-0.2 with the Kuiper Airborne Observatory and reobserved the central position with the Infrared Space Observatory's Long Wavelength Spectrometer in far-infrared lines of S++, N++, N+, and O++. We also added the N+ lines at 122 and 205 micron to the suite of lines measured in W43 by Simpson et al. (1995). The measured electron densities range from ~40 to over 4000 cm-3 in a single H II region, indicating that abundance analyses must consider density variations, since the critical densities of the observed lines range from 40 to 9000 cm-3. We propose a method to handle density variations and make new estimates of the S/H and N/H abundance ratios. We find that our sulfur abundance estimates for G333.6-0.2 and W43 agree with the S/H abundance ratios expected for the S/H abundance gradient previously reported by Simpson et al. The estimated N/H, S/H, and N/S ratios are the most reliable because of their small corrections for unseen ionization states (< 10%). We compute models of the two H II regions to estimate corrections for the other unseen ionization states. We discuss these predictions and conclude that only a few of the latest models adequately reproduce H II region observations, including the well-known, relatively-large observed Ne++/O++ ratios in low- and moderate-excitation H II regions.
[161]  [pdf] - 64214
Mid-infrared selection of AGN
Comments: 4 pages, 2 figures, accepted for publication as Letter in Astronomy & Astrophysics
Submitted: 2004-04-15
Since a large fraction of active galactic nuclei (AGN) is missed in common UV-excess surveys and is even hard to find in radio, near-IR and X-ray surveys, we have used a new AGN selection technique which is expected to be not affected by extinction. Within the scientific verification of the ISOCAM Parallel Survey at 6.7 micron we have discovered objects with exceptional mid-infrared (MIR) emission. They are essentially not detected on IRAS-ADDSCANs and only very few of them show up in the NVSS and FIRST radio surveys. Various colour criteria of the 6.7 micron data with 2MASS and optical wavebands show that the sources reach more extreme IR colours than the sources in the Hubble Deep Field-South and the ELAIS survey. The comparison with known object types suggests that we have found AGN with a pronounced MIR emission, probably due to circum-nuclear dust. First results from optical spectroscopy of ten candidates corroborate this interpretation showing four AGN, two reddened LINER and four extremely reddened emission-line galaxies with MIR/FIR flux ratios higher than for known pure starburst galaxies. The results will make a significant contribution to the debate on the entire AGN population.
[162]  [pdf] - 63948
ISOCAM survey and dust models of 3CR radio galaxies and quasars
Comments: Accepted by A&A, 30pages, 10 Figures, available at:
Submitted: 2004-04-02
We present a survey of all 3CR sources imaged with ISOCAM onboard the Infrared Space Observatory (ISO). For the source, we present spatially integrated mid--infrared (MIR, 5 - 18mic.) fluxes measured from newly calibrated ISOCAM images. In total, we detected 68 AGN of the 3CR catalogue, at redshifts z < 2.5. The one with the highest redshift is 4C+72.26 at z = 3.53. ISOCAM data are combined with other photometric measurements to construct the spectral energy distribution (SED) from optical to radio wavelengths. In order to describe dust emission we apply new radiative transfer models. By varying three parameters, luminosity, effective size and extinction, we obtain a fit to the SED for our objects. In the MIR the hot dust component is mainly due to small grains and PAHs. In the models, a type 1 AGN is represented by a compact dust distribution, the dust is therefore very warm and emission of PAHs is weak because of photo--destruction. In AGNs of type 2, the dust is relatively colder but PAH bands are strong.
[163]  [pdf] - 57639
The spiral galaxy M33 mapped in the FIR by ISOPHOT: A spatially resolved study of the warm and cold dust
Comments: 10 pages, 9 figures; the full resolution images of Fig.1 can be obtained by anonymous ftp from /pub/iso-arch
Submitted: 2003-06-26
The Sc galaxy M33 has been mapped with ISOPHOT in the far-infrared, at 60, 100, and 170mue. The spatial resolution of these FIR maps allows the separation of spiral arms and interarm regions and the isolation of a large number of star-forming regions. The spectral energy distribution in the FIR indicates a superposition of two components, a warm one originating from dust at ~45K, and a cold one, at ~16K. The warm component is concentrated towards the spiral arms and the star-forming regions, and is likely heated by the UV radiation from OB stars. The cold component is more smoothly distributed over the disk, and heated by the diffuse interstellar radiation. For the about 60 star-forming regions detected the H-alpha/FIR flux ratio increases significantly with the distance from the galaxy center, probably due to decreasing extinction. An anti-correlation of F_Ha/F_60 with F_170 suggests the intrinsic extinction to be related to the cold dust surface brightness according to A_V/S_170~0.03mag/MJy*sr. For the total galaxy the star formation rate (SFR) derived from the FIR is in agreement with that derived from the de-extincted H-alpha emission. For individual star-forming regions, a consistency between SFRs derived from the optical and from the FIR requires only a fraction of the UV radiation to be absorbed locally. The individual star-forming regions also show a local radio-FIR correlation. This local correlation is, however, due to quite different components than to those that lead to the well-known global radio-FIR correlation for entire galaxies.
[164]  [pdf] - 55815
Dust Temperatures in the Infrared Space Observatory Atlas of Bright Spiral Galaxies
Comments: AJ, 2003, in press
Submitted: 2003-03-26
We examine far-infrared and submillimeter spectral energy distributions for galaxies in the Infrared Space Observatory Atlas of Bright Spiral Galaxies. For the 71 galaxies where we had complete 60-180 micron data, we fit blackbodies with lambda^-1 emissivities and average temperatures of 31 K or lambda^-2 emissivities and average temperatures of 22 K. Except for high temperatures determined in some early-type galaxies, the temperatures show no dependence on any galaxy characteristic. For the 60-850 micron range in eight galaxies, we fit blackbodies with lambda^-1, lambda-2, and lambda^-beta (with beta variable) emissivities to the data. The best results were with the lambda^-beta emissivities, where the temperatures were ~30 K and the emissivity coefficient beta ranged from 0.9 to 1.9. These results produced gas to dust ratios that ranged from 150 to 580, which were consistent with the ratio for the Milky Way and which exhibited relatively little dispersion compared to fits with fixed emissivities.
[165]  [pdf] - 51694
A very young star forming region detected by the ISOPHOT Serendipity Survey
Comments: 14 pages, 14 figures, 7 tables. Accepted for publication in Astronomy & Astrophysics
Submitted: 2002-09-15
We present a multi-wavelength study of the star forming region ISOSS J 20298+3559, which was identified by a cross-correlation of cold compact sources from the 170 micron ISOPHOT Serendipity Survey (ISOSS) database coinciding with objects detected by the MSX, 2MASS and IRAS infrared surveys. ISOSS J 20298+3559 is associated with a massive dark cloud complex (M ~ 760 M$_{\odot}$) and located in the Cygnus X giant molecular cloud. We derive a distance of 1800 pc on the basis of optical extinction data. The low average dust temperature (T ~ 16 K) and large mass (M ~ 120 M$_{\odot}$) of the dense inner part of the cloud, which has not been dispersed, indicates a recent begin of star formation. The youth of the region is supported by the early evolutionary stage of several pre- and protostellar objects discovered across the regio n: I) Two candidate Class 0 objects with masses of 8 and 3.5 M$_{\odot}$, II) a gravitationally bound, cold (T ~ 12 K) and dense (n(H$_{2}$) \~ 2 x 10$^{5}$ cm$^{-3}$) cloud core with a mass of 50 M$_{\odot}$ and III) a Herbig B2 star with a mass of 6.5 M$_{\odot}$ and a bolometric luminosity of 2200 L$_{\odot}$, showing evidence for ongoing accretion and a stellar age of less than 40000 years. The dereddened SED of the Herbig star is well reproduced by an accretion disc + star model. The externally heated cold cloud core is a good candidate for a massive pre-protostellar object. The star formation efficiency in the central cloud region is about 14 %.
[166]  [pdf] - 1232908
Star Formation in the ISO Atlas of Spiral Galaxies
Comments: Accepted by The Astronomical Journal, 20 pages, 14 figures
Submitted: 2002-06-17, last modified: 2002-06-18
We investigate star formation along the Hubble sequence using the ISO Atlas of Spiral Galaxies. Using mid-infrared and far-infrared flux densities normalized by K-band flux densities as indicators of recent star formation, we find several trends. First, star formation activity is stronger in late-type (Sc - Scd) spirals than in early-type (Sa - Sab) spirals. This trend is seen both in nuclear and disk activity. These results confirm several previous optical studies of star formation along the Hubble sequence but conflict with the conclusions of most of the previous studies using IRAS data, and we discuss why this might be so. Second, star formation is significantly more extended in later-type spirals than in early-type spirals. We suggest that these trends in star formation are a result of differences in the gas content and its distribution along the Hubble sequence, and it is these differences that promote star formation in late-type spiral galaxies. We also search for trends in nuclear star formation related to the presence of a bar or nuclear activity. The nuclear star formation activity is not significantly different between barred and unbarred galaxies. We do find that star formation activity appears to be inhibited in LINERs and transition objects compared to HII galaxies. The mean star formation rate in the sample is 1.4 Msun/yr based on global far-infrared fluxes. Combining these data with CO data gives a mean gas consumption time of 6.4 x 10^8 yr, which is ~5 times lower than the values found in other studies. Finally, we find excellent support for the Schmidt Law in the correlation between molecular gas masses and recent star formation in this sample of spiral galaxies.
[167]  [pdf] - 44948
Determination of limits on disc masses around six pulsars at 15 and 90 microns
Comments: 7 pages, 2 figures, revised version, accepted for publication in A&A
Submitted: 2001-09-25, last modified: 2002-03-21
We have searched for evidence of emission at 15 microns with ISOCAM and at 90 microns with ISOPHOT from dust orbiting six nearby pulsars, both in binaries and in isolation, located at distances between about 100 to 1000 pc. No emission was detected at any of the pulsar positions, and for the nearest pulsar J0108-1431 the 3 sigma upper limits on the flux density is about 66 mJy at 15 microns and 22.5 mJy at 90 microns. Upper limits on the masses of circumpulsar dust are inferred at a given temperature using a simple modelling of the radiated flux; they are compared to upper limits of orbiting mass obtained with the dust heating model of Foster & Fisher (1996). These results suggest that it is unlikely that any of them have sufficiently massive, circumpulsar discs out of which planets may form in the future.
[168]  [pdf] - 45271
Infrared to millimetre photometry of ultra-luminous IR galaxies: new evidence favouring a 3-stage dust model
Comments: 23 pages, 11 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2001-10-09
Infrared to millimetre spectral energy distributions have been obtained for 41 bright ultra-luminous infrared galaxies. The observations were carried out with ISOPHOT between 10 and 200 micron and supplemented for 16 sources with SCUBA at 450 and 850 micron and with SEST at 1.3 mm. In addition, seven sources were observed at 1.2 and 2.2 $\mu$m with the 2.2 m telescope on Calar Alto. These new SEDs represent the most complete set of infrared photometric templates obtained so far on ULIRGs in the local universe.
[169]  [pdf] - 41044
Dust emission from 3C radio galaxies and quasars: New ISO observations favour the unified scheme
Comments: 10 pages, 6 figures, Latex
Submitted: 2001-02-20
In order to test the unified scheme for luminous radio galaxies and quasars we observed 10 galaxy/quasar pairs from the 3CR catalogue with ISOPHOT at infrared wavelengths between 5 and 180 micron. Each pair was selected such that both the 178 MHz luminosity and the redshift match as close as possible between the radio galaxy and the quasar in order to minimize effects of cosmic evolution. 13 of the 20 sources were detected in at least one waveband. 12 sources show clear evidence of a thermal bump at FIR wavelength, while in the remaining 7 sources the upper limits are still compatible with the presence of luminous dust emission. In agreement with the predictions of the unified scheme, the quasars and galaxies in our sample cannot be distinguished by their observed mid- and far-infrared properties. This is in contrast to the IRAS results which indicated that radio galaxies radiate significantly less mid- to far-infrared emission than quasars. However, the IRAS samples are dominated by low-redshift sources (z < 0.5), while our sample contains several of the most luminous radio galaxies at redshift z ~ 1. The latter have already been suspected to contain a hidden quasar for other reasons, e.g. an extended emission line region aligned with the radio axis. From the ratio between FIR luminosity emitted by dust and the radio power at 178 MHz, we conclude that the radio galaxy/quasar unification might be perfectly valid for the most luminous 3C sources at high redshift (z > 0.8). At lower redshifts (z < 0.5), however, some of the FRII radio galaxies contain active nuclei which emit less UV-optical continuum than the quasars of similar radio power. (abrigded)
[170]  [pdf] - 38336
The bright 175um knots of the Andromeda Galaxy
Comments: 9 pages. 12 figures, accepted by A&A
Submitted: 2000-09-29
Discrete far-infrared (FIR) sources of M31 are identified in the ISO 175um map and characterized via their FIR colours, luminosities and masses. With a mean size of 800pc, they probably represent several clouds in chance projection or giant cloud complexes. The 175um data point provides crucial information in addition to the IRAS 60 and 100um data: At least two modified Planck curves with temperatures of about 40K and 15-21K are necessary to fit the SEDs of the knots. We distinguish between three types of knots - cold, medium, warm - in order to recognize trends. Comparisons with radio and optical tracers show that - statistically - the cold knots can be identified well with CO and HI radio sources and thus might represent mainly molecular cloud complexes. The warm knots coincide with known HII regions and supernova remnants. The medium knots might contain a balanced mixture of molecular clouds and HII regions. The cold knots have a considerable luminosity and their discovery raises the question of hidden star formation. Though the optically dark dust lanes in M31 generally match the FIR ring, surprisingly we do not find a convincing coincidence of our knots with individual dark clouds which might therefore show mainly foreground dust features. The ratio of FIR luminosity to dust mass, L/M, is used to measure the energy content of the dust. The knots have a clear L/M excess over the rest of M31, providing evidence that they are powered by star formation in addition to the interstellar radiation field. The L/M ratio of the warm knots is comparable to that of Galactic HII regions like M42 or NGC2024, while that of the cold knots still reaches values like in the average Orion complex. Thus both the warm and the cold knots are interpreted as containing large cloud complexes with considerable ongoing star formation.