Brought to you by:

Table of contents

Volume 741

Number 1, 2011 November 1

Previous issue Next issue

L1

, , , , , , , , , et al

Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here, we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star's rotation period is 35.1 ± 1.0 days, and its projected obliquity with respect to the stellar binary orbit is 1fdg6 ± 2fdg4. Therefore, the three largest sources of angular momentum—the stellar orbit, the planetary orbit, and the primary's rotation—are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the "pseudosynchronous" period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2–4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

L2

, , and

The observed population of Hot Jupiters displays a stunning variety of physical properties, including a wide range of densities and core sizes for a given planetary mass. Motivated by the observational sample, this Letter studies the accretion of rocky planets by Hot Jupiters, after the Jovian planets have finished their principal migration epoch and become parked in ∼4 day orbits. In this scenario, rocky planets form later and then migrate inward due to torques from the remaining circumstellar disk, which also damps the orbital eccentricity. This mechanism thus represents one possible channel for increasing the core masses and metallicities of Hot Jupiters. This Letter determines probabilities for the possible end states for the rocky planet: collisions with the Jovian planets, accretion onto the star, ejection from the system, and long-term survival of both planets. These probabilities depend on the mass of the Jovian planet and its starting orbital eccentricity, as well as the eccentricity damping rate for the rocky planet. Since these systems are highly chaotic, a large ensemble (N ∼ 103) of simulations with effectively equivalent starting conditions is required. Planetary collisions are common when the eccentricity damping rate is sufficiently low, but are rare otherwise. For systems that experience planetary collisions, this work determines the distributions of impact velocities—both speeds and impact parameters—for the collisions. These velocity distributions help determine the consequences of the impacts, e.g., where energy and heavy elements are deposited within the giant planets.

L3

, , and

We investigate the pulsation driving mechanism responsible for the long-period photometric variations observed in LS IV-14°116, a subdwarf B star showing a He-enriched atmospheric composition. To this end, we perform detailed nonadiabatic pulsation computations over fully evolutionary post-He-core-flash stellar structure models, appropriate for hot subdwarf stars at evolutionary phases previous to the He-core burning stage. We found that the variability of LS IV-14°116 can be attributed to non-radial g-mode pulsations excited by the epsilon-mechanism acting in the He-burning shells that appear before the star settles in the He-core burning stage. Even more interestingly, our results show that LS IV-14°116 could be the first known pulsating star in which the epsilon-mechanism of mode excitation is operating. Last but not the least, we find that the period range of destabilized modes is sensitive to the exact location of the burning shell, something that might help in distinguishing between the different evolutionary scenarios proposed for the formation of this star.

L4

, , , , , , and

We present results from an 87 ks Suzaku observation of the canonical low-excitation radio galaxy (LERG) NGC 6251. We have previously suggested that LERGs violate conventional active galactic nucleus unification schemes: they may lack an obscuring torus and are likely to accrete in a radiatively inefficient manner, with almost all of the energy released by the accretion process being channeled into powerful jets. We model the 0.5–20 keV Suzaku spectrum with a single power law of photon index Γ = 1.82+0.04− 0.05, together with two collisionally ionized plasma models whose parameters are consistent with the known galaxy- and group-scale thermal emission. Our observations confirm that there are no signatures of obscured, accretion-related X-ray emission in NGC 6251, and we show that the luminosity of any such component must be substantially sub-Eddington in nature.

L5

, , , , , , , , , et al

We report on the extremely intense and fast gamma-ray flare above 100 MeV detected by AGILE from the Crab Nebula in mid-April 2011. This event is the fourth of a sequence of reported major gamma-ray flares produced by the Crab Nebula in the period 2007/mid-2011. These events are attributed to strong radiative and plasma instabilities in the inner Crab Nebula, and their properties are crucial for theoretical studies of fast and efficient particle acceleration up to 1015 eV. Here we study the very rapid flux and spectral evolution of the event that on 2011 April 16 reached the record-high peak flux of F = (26 ± 5) × 10−6 photons cm−2 s−1 with a rise-time timescale that we determine to be in the range 6–10 hr. The peak flaring gamma-ray spectrum reaches a distinct maximum near 500 MeV with no substantial emission above 1 GeV. The very rapid rise time and overall evolution of the Crab Nebula flare strongly constrain the acceleration mechanisms and challenge MHD models. We briefly discuss the theoretical implications of our observations.

L6

, , , and

Results of modeling the heliosphere are compared with observations of the Voyager spacecraft and the IBEX mission simultaneously. The MHD solutions are tested against observational data for different strengths and orientations of the local interstellar magnetic field (LIMF) used in the simulations for asymmetric magnetized solar wind flow. We show that the model reproduces approximately the position of the IBEX ribbon and the termination shock crossing distance for Voyager 2, when the LIMF vector lies in the proximity of the hydrogen deflection plane with the inclination angle to the local interstellar flow equal to 39°  ±  9° and its magnitude is 2.4 ± 0.3 μG. In ecliptic coordinates this solution corresponds to the LIMF vector pointing from (longitude, latitude) = (227° ± 7°, 35° ± 7°).

L7

and

We present Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) observations of EUV cyclones in the quiet Sun. These cyclones are rooted in the rotating network magnetic fields (RNFs). Such cyclones can last several to more than 10 hr and, at the later phase, they are found to be associated with EUV brightenings (microflares) and even EUV waves. SDO Helioseismic and Magnetic Imager (HMI) observations show a ubiquitous presence of RNFs. Using HMI line-of-sight magnetograms on 2010 July 8, we find 388 RNFs in an area of 800 × 980 arcsec2 near the disk center where no active region is present. The sense of rotation shows a weak hemisphere preference. The unsigned magnetic flux of the RNFs is about 4.0 × 1021 Mx, or 78% of the total network flux. These observational phenomena at small scale reported in this Letter are consistent with those at large scale in active regions. The ubiquitous RNFs and EUV cyclones over the quiet Sun may suggest an effective way to heat the corona.

L8

, , , and

We present a spectroscopic analysis of Very Large Telescope/X-Shooter observations of six O-type stars in the low-metallicity (Z ∼ 1/7 Z) galaxies IC 1613, WLM, and NGC 3109. The stellar and wind parameters of these sources allow us, for the first time, to probe the mass loss versus metallicity dependence of stellar winds at metallicities below that of the Small Magellanic Cloud (at Z ∼ 1/5 Z) by means of a modified wind momentum versus luminosity diagram. The wind strengths that we obtain for the objects in WLM and NGC 3109 are unexpectedly high and do not agree with theoretical predictions. The objects in IC 1613 tend toward a higher than expected mass-loss rate, but remain consistent with predictions within their error bars. We discuss potential systematic uncertainties in the mass-loss determinations to explain our results. However, if further scrutinization of these findings point towards an intrinsic cause for this unexpected sub-SMC mass-loss behavior, implications would include a higher than anticipated number of Wolf–Rayet stars and Ib/Ic supernovae in low-metallicity environments, but a reduced number of long-duration gamma-ray bursts produced through a single-star evolutionary channel.

L9

, , , , , and

In the interstellar medium (ISM), an important channel of water formation is the reaction of atoms on the surface of dust grains. Here, we report on a laboratory study of the formation of water via the O+D reaction network. While prior studies were done on ices, as appropriate to the formation of water in dense clouds, we explored how water formation occurs on bare surfaces, i.e., in conditions mimicking the transition from diffuse to dense clouds (Av ∼ 1–5). Reaction products were detected during deposition and afterward when the sample is brought to a high temperature. We quantified the formation of water and intermediary products, such as D2O2, over a range of surface temperatures (15–25 K). The detection of OD on the surface signals the importance of this reactant in the overall scheme of water formation in the ISM.

L10

, , , and

Titan is the only satellite that possesses a thick atmosphere, composed mainly of N2 and CH4. However, its origin and evolution remain largely unknown. Knowledge of the acquirement of a N2 atmosphere on Titan would provide insights into nitrogen evolution in planetary atmospheres as well as the formation of satellite systems around gas giants. Previous studies have proposed that the atmospheric N2 would have been converted from NH3 via shock heating by accreting satellitesimals in the highly reducing proto-atmosphere composed of NH3 and CH4. Nevertheless, the validity of this mechanism strongly depends on both the composition of the proto-atmosphere and kinetics of shock chemistry. Here, we show that a CO2-rich oxidizing proto-atmosphere is necessary to form N2 from NH3 efficiently by atmospheric shock heating. Efficient shock production of N2 is inhibited in a reducing proto-atmosphere composed of NH3 and CH4, because CH4 plays as the coolant gas owing to its large heat capacity. Our calculations show that the amount of N2 produced in a CO2-rich proto-atmosphere could have reached ∼20 times that on the present Titan. Although further quantitative analysis are required (especially, the occurrence of catalytic reactions), our results imply that the chemical composition of satellitesimals that formed the Saturnian system is required to be oxidizing if the current atmospheric N2 is derived from the shock heating in the proto-atmosphere during accretion. This supports the formation of regular satellites in an actively supplied circumplanetary disk using CO2-rich materials originated from the solar nebula at the final stage of gas giant formation.

L11

, , , , , , , , and

We present new constraints on the ratio of black hole (BH) mass to total galaxy stellar mass at 0.3 < z < 0.9 for a sample of 32 type-1 active galactic nuclei (AGNs) from the XMM-COSMOS survey covering the range MBH ∼ 107.2 − 8.7M. Virial MBH estimates based on Hβ are available from the COSMOS Magellan/IMACS survey. We use high-resolution Hubble Space Telescope (HST) imaging to decompose the light of each type-1 AGN and host galaxy, and employ a specially built mass-to-light ratio to estimate the stellar masses (M*). The MBHM* ratio shows a zero offset with respect to the local relation for galactic bulge masses, and we also find no evolution in the mass ratio MBH/M*∝(1 + z)0.02 ± 0.34 up to z ∼ 0.9. Interestingly, at the high-MBH end there is a positive offset from the z = 0 relation, which can be fully explained by a mass function bias with a cosmic scatter of σμ = 0.3, reaffirming that the intrinsic distribution is consistent with zero evolution. From our results we conclude that since z ∼ 0.9 no substantial addition of stellar mass is required: the decline in star formation rates and merger activity at z < 1 support this scenario. Nevertheless, given that a significant fraction of these galaxies show a disk component, their bulges are indeed undermassive. This is a direct indication that for the last 7 Gyr the only essential mechanism required for these galaxies to obey the z = 0 relation is a redistribution of stellar mass to the bulge, likely driven by secular processes, i.e., internal instabilities and minor merging.

L12

and

Among the most explored directions in the study of dense stellar systems is the investigation of the effects of the retention of supernova remnants, especially that of the massive stellar remnant black holes (BHs), in star clusters. By virtue of their eventual high central concentration, these stellar mass BHs potentially invoke a wide variety of physical phenomena, the most important ones being emission of gravitational waves (GWs), formation of X-ray binaries, and modification of the dynamical evolution of the cluster. Here we propose, for the first time, that rapid removal of stars from the outer parts of a cluster by the strong tidal field in the inner region of our Galaxy can unveil its BH sub-cluster, which appears as a star cluster that is gravitationally bound by an invisible mass. We study the formation and properties of such systems through direct N-body computations and estimate that they can be present in significant numbers in the inner region of the Milky Way. We call such objects "dark star clusters" (DSCs) as they appear dimmer than normal star clusters of similar mass and they comprise a predicted, new class of entities. The finding of DSCs will robustly cross-check BH retention; they will not only constrain the uncertain natal kicks of BHs, thereby the widely debated theoretical models of BH formation, but will also pinpoint star clusters as potential sites for GW emission for forthcoming ground-based detectors such as the Advanced LIGO. Finally, we also discuss the relevance of DSCs for the nature of IRS 13E.

L13

, , , , and

PSR J1734−3333 is a radio pulsar rotating with a period P = 1.17 s and slowing down with a period derivative $\dot{P}=2.28\times 10^{-12}$, the third largest among rotation-powered pulsars. These properties are midway between those of normal rotation-powered pulsars and magnetars, two populations of neutron stars that are notably different in their emission properties. Here we report on the measurement of the second period derivative of the rotation of PSR J1734−3333 and calculate a braking index n = 0.9 ± 0.2. This value is well below 3, the value expected for an electromagnetic braking due to a constant magnetic dipole, and indicates that this pulsar may soon have the rotational properties of a magnetar. While there are several mechanisms that could lead to such a low braking index, we discuss this observation, together with the properties exhibited by some other high-$\dot{P}$ rotation-powered pulsars, and interpret it as evidence of a possible evolutionary route for magnetars through a radio-pulsar phase, supporting a unified description of the two classes of the object.

L14

, , , , , and

In this Letter, we introduce a technique for finding resonance radii in a disk galaxy. We use a two-dimensional velocity field in Hα emission obtained with Fabry–Perot interferometry, derive the classical rotation curve, and subtract it off, leaving a residual velocity map. As the streaming motions should reverse sign at corotation, we detect these reversals and plot them in a histogram against galactocentric radius, excluding points where the amplitude of the reversal is smaller than the measurement uncertainty. The histograms show well-defined peaks which we assume to occur at resonance radii, identifying corotations as the most prominent peaks corresponding to the relevant morphological features of the galaxy (notably bars and spiral arm systems). We compare our results with published measurements on the same galaxies using other methods and different types of data.

L15

, , , and

We present radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C303. From these data we derive the magnetoplasma and electrodynamic parameters of this 50 kpc long jet. For one component of this jet we obtain for the first time a direct determination of a galactic-scale electric current (∼3 × 1018 A), and its direction—positive away from the active galactic nucleus. Our analysis strongly supports a model where the jet energy flow is mainly electromagnetic.

L16

, , , , , , , , , and

We report the discovery of the first identified pulsating DA white dwarf, WD J1916+3938 (Kepler ID 4552982), in the field of the Kepler mission. This ZZ Ceti star was first identified through ground-based, time-series photometry, and follow-up spectroscopy confirms that it is a hydrogen-atmosphere white dwarf with Teff = 11,129 ± 115 K and log g = 8.34 ± 0.06, placing it within the empirical ZZ Ceti instability strip. The object shows up to 0.5% amplitude variability at several periods between 800 and 1450 s. Extended Kepler observations of WD J1916+3938 could yield the best light curve, to date, of any pulsating white dwarf, allowing us to directly study the interior of an evolved object representative of the fate of the majority of stars in our Galaxy.

L17

, , , , and

There are few observational constraints on how the escape of ionizing photons from starburst galaxies depends on galactic parameters. Here we report on the first major detection of an ionization cone in NGC 5253, a nearby starburst galaxy. This high-excitation feature is identified by mapping the emission-line ratios in the galaxy using [S iii] λ9069, [S ii] λ6716, and Hα narrowband images from the Maryland-Magellan Tunable Filter at Las Campanas Observatory. The ionization cone appears optically thin, which suggests the escape of ionizing photons. The cone morphology is narrow with an estimated solid angle covering just 3% of 4π steradians, and the young, massive clusters of the nuclear starburst can easily generate the radiation required to ionize the cone. Although less likely, we cannot rule out the possibility of an obscured active galactic nucleus source. An echelle spectrum along the minor axis shows complex kinematics that are consistent with outflow activity. The narrow morphology of the ionization cone supports the scenario that an orientation bias contributes to the difficulty in detecting Lyman continuum emission from starbursts and Lyman break galaxies.

L18

, , , , and

We report on the follow-up XMM-Newton observation of the planet-hosting star HD 189733 we obtained in 2011 April. We observe a flare just after the secondary transit of the hot Jupiter. This event shares the same phase and many of the characteristics of the flare we observed in 2009. We suggest that a systematic interaction between planet and stellar magnetic fields when the planet passes close to active regions on the star can lead to periodic variability phased with planetary motion. By means of high-resolution X-ray spectroscopy with the Reflection Grating Spectrometer on board XMM-Newton, we determine that the corona of this star is unusually dense.

L19

, , , and

The question is addressed as to what extent incompressible magnetohydrodynamics can describe random magnetic and velocity fluctuations measured in the solar wind. It is demonstrated that distributions of spectral indices for the velocity, magnetic field, and total energy obtained from high-resolution numerical simulations of magnetohydrodynamic turbulence are qualitatively and quantitatively similar to solar wind observations at 1 AU. Both simulations and observations show that in the inertial range the magnetic field spectrum Eb is steeper than the velocity spectrum Ev with EbEv and that the magnitude of the residual energy ER = EvEb decreases nearly following a k−2 scaling.

L20

, , , , , , , and

At high redshift, the universe is so young that core-collapse supernovae (SNe) are suspected to be the dominant source of dust production. However, some observations indicate that the dust production by SNe is an inefficient process, casting doubts on the existence of abundant SNe-dust in the early universe. Recently, Perley et al. reported that the afterglow of GRB 071025—an unusually red gamma-ray burst (GRB) at z ∼ 5—shows evidence for SNe-produced dust. Since this is perhaps the only high-redshift GRB exhibiting compelling evidence for SNe-dust but the result could easily be affected by small systematics in photometry, we re-examined the extinction properties of GRB 071025 using our own optical/near-infrared data at a different epoch. In addition, we tested SNe-dust models with different progenitor masses and dust destruction efficiencies to constrain the dust formation mechanisms. By searching for the best-fit model of the afterglow spectral energy distribution, we confirm the previous claim that the dust in GRB 071025 is most likely to originate from SNe. We also find that the SNe-dust model of 13 or 25 M without dust destruction fits the extinction property of GRB 071025 best, while pair-instability SNe models with a 170 M progenitor poorly fit the data. Our results indicate that, at least in some systems at high redshift, SNe with intermediate initial masses within 10–30 M were the main contributors for the dust enrichment, and the dust destruction effect due to reverse shock was negligible.

L21

, , and

Coronal bright fronts (CBFs) are large-scale wavefronts that propagate through the solar corona at hundreds of kilometers per second. While their kinematics have been studied in detail, many questions remain regarding the temporal evolution of their amplitude and pulse width. Here, contemporaneous high cadence, multi-thermal observations of the solar corona from the Solar Dynamic Observatory (SDO) and Solar TErrestrial RElations Observatory (STEREO) spacecraft are used to determine the kinematics and expansion rate of a CBF wavefront observed on 2010 August 14. The CBF was found to have a lower initial velocity with weaker deceleration in STEREO observations compared to SDO observations (∼340 km s−1 and −72 m s−2 as opposed to ∼410 km s−1 and −279 m s−2). The CBF kinematics from SDO were found to be highly passband-dependent, with an initial velocity ranging from 379 ± 12 km s−1 to 460 ± 28 km s−1 and acceleration ranging from −128 ± 28 m s−2 to −431 ± 86 m s−2 in the 335 Å and 304 Å passbands, respectively. These kinematics were used to estimate a quiet coronal magnetic field strength range of ∼1–2 G. Significant pulse broadening was also observed, with expansion rates of ∼130 km s−1 (STEREO) and ∼220 km s−1 (SDO). By treating the CBF as a linear superposition of sinusoidal waves within a Gaussian envelope, the resulting dispersion rate of the pulse was found to be ∼8–13 Mm2 s−1. These results are indicative of a fast-mode magnetoacoustic wave pulse propagating through an inhomogeneous medium.

L22

and

The observed similarities between the mass function of prestellar cores (CMF) and the stellar initial mass function (IMF) have led to the suggestion that the IMF is already largely determined in the gas phase. However, theoretical arguments show that the CMF may differ significantly from the IMF. In this Letter, we study the relation between the CMF and the IMF, as predicted by the IMF model of Padoan and Nordlund. We show that (1) the observed mass of prestellar cores is on average a few times smaller than that of the stellar systems they generate; (2) the CMF rises monotonically with decreasing mass, with a noticeable change in slope at approximately 3–5 M, depending on mean density; (3) the selection of cores with masses larger than half their Bonnor–Ebert mass yields a CMF approximately consistent with the system IMF, rescaled in mass by the same factor as our model IMF, and therefore suitable to estimate the local efficiency of star formation, and to study the dependence of the IMF peak on cloud properties; and (4) only one in five pre-brown-dwarf core candidates is a true progenitor to a brown dwarf.

L23

and

We present preliminary though statistically significant evidence that shows that multiplanetary systems that exhibit a 2/1 period commensurability are in general younger than multiplanetary systems without commensurabilities, or even systems with other commensurabilities. An immediate possible conclusion is that the 2/1 mean-motion resonance in planetary systems tends to be disrupted after typically a few Gyrs.

L24

, , , , , , , , , et al

Strange-looking dust cloud around asteroid (596) Scheila was discovered on 2010 December 11.44–11.47. Unlike normal cometary tails, it consisted of three tails and faded within two months. We constructed a model to reproduce the morphology of the dust cloud based on the laboratory measurement of high-velocity impacts and the dust dynamics. As a result, we succeeded in reproducing the peculiar dust cloud by an impact-driven ejecta plume consisting of an impact cone and downrange plume. Assuming an impact angle of 45°, our model suggests that a decameter-sized asteroid collided with (596) Scheila from the direction of (αim, δim) = (60°, −40°) in J2000 coordinates on 2010 December 3. The maximum ejection velocity of the dust particles exceeded 100 m s−1. Our results suggest that the surface of (596) Scheila consists of materials with low tensile strength.