Table of contents

Volume 702

Number 2, 2009 September 10

Previous issue Next issue

825

, , , , and

We perform an autocorrelation study of the Auger data with the aim to constrain the number density ns of ultrahigh energy cosmic ray (UHECR) sources, estimating at the same time the effect on ns of the systematic energy scale uncertainty and of the distribution of UHECR. The use of global analysis has the advantage that no biases are introduced, either in ns or in the related error bar, by the a priori choice of a single angular scale. The case of continuous, uniformly distributed sources is nominally disfavored at 99% CL and the fit improves if the sources follow the large-scale structure of matter in the universe. The best-fit values obtained for the number density of proton sources are within a factor ∼2 around ns ≃ 1 × 10−4Mpc−3 and depend mainly on the Auger energy calibration scale, with lower densities being preferred if the current scale is correct. The data show no significant small-scale clustering on scales smaller than a few degrees. This might be interpreted as a signature of magnetic smearing of comparable size, comparable with the indication of a ≈3° magnetic deflection coming from cross-correlation results. The effects of some approximations done on the above results are also discussed.

833

, , , , , , and

We present sensitive near-infrared (NIR) VLT ISAAC spectroscopic observations of the z = 6.08 quasar SDSS J030331.40-001912.9. This quasi-stellar object (QSO) is more than a magnitude fainter than other QSOs at z ∼ 6 for which NIR spectroscopy has been obtained to date and is therefore presumably more representative of the QSO population at the end of cosmic reionization. Combining rest-frame UV continuum luminosity with the width measurements of the Mg ii and C iv lines, we derive a black hole mass of 2+1.0−0.5 × 108M, the lowest mass observed for z ∼ 6 QSOs to date, and derive an Eddington ratio of 1.6+0.4−0.6, among the highest value derived for QSOs at any redshift. The Spitzer 24 μm nondetection of this QSO does not leave space for a significant hot dust component in its optical/NIR spectral energy distribution, in common with one other faint QSO at z = 6, but in contrast to more than 20 more z = 6 QSOs and all known lower redshift QSOs with sufficiently deep multiwavelength photometry. We conclude that we have found evidence for differences in the intrinsic properties of at least one z ∼ 6 QSO as compared to the lower redshift population.

838

, , , , , , , , , et al

Spectra in the extreme ultraviolet range from 107 to 353 Å emitted from Fe ions in various ionization stages have been observed at the Heidelberg electron beam ion trap (EBIT) with a flat-field grating spectrometer. A series of transition lines and their intensities have been analyzed and compared with collisional-radiative simulations. The present collisional-radiative model reproduces well the relative line intensities and facilitates line identification of ions produced in the EBIT. The polarization effect on the line intensities resulting from nonthermal unidirectional electron impact was explored and found to be significant (up to 24%) for a few transition lines. Based upon the observed line intensities, relative charge state distributions (CSD) of ions were determined, which peaked at Fe23+ tailing toward lower charge states. Another simulation on ion charge distributions including the ionization and electron capture processes generated CSDs which are in general agreement with the measurements. By observing intensity ratios of specific lines from levels collisionally populated directly from the ground state and those starting from the metastable levels of Fe xxi, Fe x and other ionic states, the effective electron densities were extracted and found to depend on the ionic charge. Furthermore, it was found that the overlap of the ion cloud with the electron beam estimated from the effective electron densities strongly depends on the charge state of the ion considered, i.e. under the same EBIT conditions, higher charge ions show less expansion in the radial direction.

851

, , , , , and

Quasar outflows may play important role in the evolution of its host galaxy and central black hole, and are most often studied in absorption lines. In this paper, we present a detailed study of multiple outflows in the obscured ultraluminous infrared quasar Q1321+058. The outflows reveal themselves in the complex optical and ultraviolet (UV) emission-line spectrum, with a broad component blueshifted by 1650 km s−1 and a narrow component by 360 km s−1, respectively. The higher velocity component shows ever strong N iii] (N iii]/C iii] = 3.8 ± 0.3 and N iii]/C iv = 0.53) and strong Si iii] (Si iii]/C iii] ≃ 1), in addition to strong [O iii]λ5007 and [Ne iii]λ3869 emission. A comparison of these line ratios with photoionization models suggests an overabundance of N and Si relative to C. The abundance pattern is consistent with a fast chemical enriching process associated with a recent starburst, triggered by a recent galaxy merger. The outflow extends to several tens to hundred parsecs from the quasar, and covers only a very small sky. We find that the outflow with line emitting gas is energetically insufficient to remove the interstellar medium of the host galaxy, but total kinetic energy may be much larger than suggested by the emission lines. The velocity range and the column density suggest that the outflow might be part of the low-ionization broad absorption line region as seen in a small class of quasars. The optical and UV continuum is starlight dominated and can be modeled with a young-aged (1 Myr) plus an intermediate-aged (∼0.5–1 Gyr) stellar populations, suggesting a fast building of the stellar mass in the host galaxy, consistent with the starburst-type metal abundances inferred from the high-velocity outflow spectrum. The broadband spectral energy distribution shows that it is an obscured quasar with its bulk emission in the middle infrared. The star formation rate, independently estimated from UV, far-infrared, and emission-line luminosity, is much lower than that is required for the co-evolution of the black hole and its host spheroid.

862

and

White-light observations of interplanetary disturbances have been dominated by interplanetary coronal mass ejections (ICMEs). This is because the other type of disturbance, the corotating interaction region (CIR), has proved difficult to detect using white-light imagers. Recently, a number of papers have appeared presenting CIR observations using the Solar Terrestrial Relations Observatory (STEREO) Heliospheric Imagers (HIs), but have mostly only focused on a single spacecraft and imager. In this paper, we present observations of a single CIR that was observed by all three current white-light heliospheric imagers (SMEI and both STEREO HIs), as well as the in situ instruments on both STEREO satellites and ACE. We begin with a discussion of the geometry of the CIR structure, and show how the apparent leading edge structure is expected to change as it corotates relative to the observer. We use these calculations to predict elongation–time profiles for CIRs of different speeds for each of the imagers, and also to predict the arrival times at the in situ instruments. We show that although all three measured different parts, they combine to produce a self-consistent picture of the CIR. Finally, we offer some thoughts on why CIRs have proved so difficult to detect in white-light heliospheric images.

871

and

We investigate the dynamical effects of dark matter subhalos on the structure and evolution of a galactic disk, using the semi-analytic method that includes approximated and empirical relations as achieved in detailed numerical simulations of the cold dark matter model. We calculate the upper limit for the size of a galactic disk at a specific redshift z, based on the orbital properties of subhalos characterized by their pericentric distances from the center of a host halo. We find that this possibly largest size of a disk as determined by the smallest pericentric distances of subhalos shows the characteristic properties, which are basically in agreement with an observed galactic disk at low and high z. Namely, it is found that a massive disk can have a larger size than a less massive one, because of its stability against the destruction effect of subhalos. Also, with fixed mass, the size of a galactic disk at low z can be larger than that at high z, reflecting the orbital evolution of subhalos with respect to a host halo. These results suggest that the presence and structure of a galactic disk may be dynamically limited by the interaction with dark matter substructures, especially at high z.

880

, , , , and

High-quality oscillator strengths are needed to convert the amount of interstellar absorption into an accurate abundance from data acquired with the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer. In order to help clarify the appropriate values to use for ultraviolet transitions in Cu ii, we report lifetimes from our beam-foil experiments. The Cu ii results for the 1358.773 Å line (3d101S0–3d94p1Po1) provide further evidence for a short lifetime for the upper level of interest. Additional results for transitions that can repopulate this upper level indicate that the data are not affected by cascades. Our measured lifetime, the most precisely determined to date, is in excellent agreement with the most recent experimental results and is consistent with theoretical values. When combined with the theoretical branching fractions of Donnelly et al. and Dong & Fritzsche, we obtain a value of 0.273 ± 0.028 for the oscillator strength, which is in very good agreement with the theoretical results of Donnelly et al. and the most recent recommended value given by Morton.

884

, , , , and

We use Newtonian N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (≲0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The Newtonian N-body simulations include a dense cluster of 10 M stellar-mass black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20%–30% of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with SBHs.

890

, , , , and

The coalescence of a massive black hole (MBH) binary leads to the gravitational-wave recoil of the system and its ejection from the galaxy core. We have carried out N-body simulations of the motion of a MBH = 3.7 × 106M MBH remnant in the "Via Lactea I" simulation, a Milky Way-sized dark matter halo. The black hole receives a recoil velocity of Vkick = 80, 120, 200, 300, and 400 km s−1 at redshift 1.5, and its orbit is followed for over 1 Gyr within a "live" host halo, subject only to gravity and dynamical friction against the dark matter background. We show that, owing to asphericities in the dark matter potential, the orbit of the MBH is highly nonradial, resulting in a significantly increased decay timescale compared to a spherical halo. The simulations are used to construct a semi-analytic model of the motion of the MBH in a time-varying triaxial Navarro–Frenk–White dark matter halo plus a spherical stellar bulge, where the dynamical friction force is calculated directly from the velocity dispersion tensor. Such a model should offer a realistic picture of the dynamics of kicked MBHs in situations where gas drag, friction by disk stars, and the flattening of the central cusp by the returning black hole are all negligible effects. We find that MBHs ejected with initial recoil velocities Vkick ≳ 500 km s−1 do not return to the host center within a Hubble time. In a Milky Way-sized galaxy, a recoiling hole carrying a gaseous disk of initial mass ∼MBH may shine as a quasar for a substantial fraction of its "wandering" phase. The long decay timescales of kicked MBHs predicted by this study may thus be favorable to the detection of off-nuclear quasar activity.

901

and

We present a three-dimensional model of the density distribution of a coronal mass ejection (CME) from 2008 April 26. This CME was observed by the two spacecraft composing the Solar Terrestrial Relations Observatory (STEREO), which tracked the CME from near the Sun, into the interplanetary medium (IPM), and all the way to 1 AU. The CME was directed toward STEREO-B and hit that spacecraft on 2008 April 29. The STEREO images of the CME show an internal structure that can be interpreted as having a flux rope shape. The two different perspectives on the event provided by the two STEREO spacecraft allow us to make a particularly strong argument for the flux rope interpretation, and the STEREO data also allow us to study the evolution of the flux rope in the IPM. The flux rope is oriented close to the ecliptic plane, but with the western leg tilted northwards by about 20°. This implies an orientation roughly perpendicular to the neutral line of the active region at the event's point of origin, apparently an unusual geometry given that previous analyses have found that CME flux ropes are usually, but not always, oriented parallel to the neutral lines of their source regions. The CME model also consists of a front out ahead of the flux rope, possibly a shock launched by the flux rope driver. The model density distribution is reasonably successful at reproducing the CME appearance both close to the Sun in coronagraphic images, and far from the Sun in images of the IPM from STEREO's heliospheric imagers. This suggests that self-similar expansion is a reasonable first-order approximation for this particular CME, and also indicates that the flux rope's orientation does not change much during its journey through the IPM.

911

, , , , and

This work presents the results of a Chandra study of 21 broad absorption line (BAL) radio-loud quasars (RLQs). We conducted a Chandra snapshot survey of 12 bright BAL RLQs selected from Sloan Digital Sky Survey Data/Faint Images of the Radio Sky data and possessing a wide range of radio and C iv absorption properties. Optical spectra were obtained nearly contemporaneously with the Hobby–Eberly Telescope; no strong flux or BAL variability was seen between epochs. In addition to the snapshot targets, we include in our sample nine additional BAL RLQs possessing archival Chandra coverage. We compare the properties of (predominantly high-ionization) BAL RLQs to those of non-BAL RLQs as well as to BAL radio-quiet quasars (RQQs) and non-BAL RQQs for context. All 12 snapshots and 8/9 archival BAL RLQs are detected, with observed X-ray luminosities less than those of non-BAL RLQs having comparable optical/UV luminosities by typical factors of 4.1–8.5. (BAL RLQs are also X-ray weak by typical factors of 2.0–4.5 relative to non-BAL RLQs having both comparable optical/UV and radio luminosities.) However, BAL RLQs are not as X-ray weak relative to non-BAL RLQs as are BAL RQQs relative to non-BAL RQQs. While some BAL RLQs have harder X-ray spectra than typical non-BAL RLQs, some have hardness ratios consistent with those of non-BAL RLQs, and there does not appear to be a correlation between X-ray weakness and spectral hardness, in contrast to the situation for BAL RQQs. RLQs are expected to have X-ray continuum contributions from both accretion-disk corona and small-scale jet emission. While the entire X-ray continuum in BAL RLQs cannot be obscured to the same degree as in BAL RQQs, we calculate that the jet is likely partially covered in many BAL RLQs. We comment briefly on implications for geometries and source ages in BAL RLQs.

929

, , , and

We present a multi-wavelength photometric and spectroscopic study of a newly discovered candidate cluster [DBS2003] 45. Our H, Ks photometry confirms that [DBS2003] 45 is a cluster. An average visual extinction AV ∼  7.1  ±   0.5 is needed to fit the cluster sequence with a model isochrone. Low-resolution spectroscopy indicates that half a dozen early B- and at least one late O-type giant stars are present in the cluster. We estimate the age of the cluster to be between 5 and 8 Myr based on spectroscopic analysis. Assuming an age of 6 Myr, we fit the observed mass function with a power law, N(M) ∝ M−Γ, and find an index Γ ∼ 1.27  ±   0.15, which is consistent with the Salpeter value. We estimate that the total cluster mass is around 103M by integrating the derived mass function between 0.5 and 45 M. Both mid-infrared and radio wavelength observations show that a bubble filled with ionized gas is associated with the cluster. The total ionizing photon flux estimated from radio continuum measurements is consistent with the number of hot stars we detected. Infrared bright point sources along the rim of the bubble suggest that there is triggered star formation at the periphery of the H ii region.

940

, , and

In the spectra of 139 early-type Large Magellanic Cloud (LMC) stars observed with Far Ultraviolet Spectroscopic Explorer and with deep radio Parkes H i 21 cm observations along with those stars, we search for and analyze the absorption and emission from high-velocity gas at +90 ⩽ vLSR ⩽ +175 km s−1. The H i column density of the high-velocity clouds (HVCs) along these sightlines ranges from <1018.4 to 1019.2 cm−2. The incidence of the HVC metal absorption is 70%, significantly higher than the H i emission occurrence of 32%. We find that the mean metallicity of the HVC is [O i/H i] = −0.51 ± 0.120.16. There is no strong evidence for a large variation in the HVC metallicity, implying that these HVCs have a similar origin and are part of the same complex. The mean and scatter of the HVC metallicities are more consistent with the present-day LMC oxygen abundance than that of the Small Magellanic Cloud or the MW. We find that on average [Si ii/O i] = +0.48 ± 0.150.25 and [Fe ii/O i] = +0.33 ± 0.140.21, implying that the HVC complex is dominantly ionized. The HVC complex has a multiphase structure with neutral (O i, Fe ii), weakly ionized (Fe ii, N ii), and highly ionized (O vi) components, and has evidence of dust but no molecules. All the observed properties of the HVC can be explained by an energetic outflow from the LMC. This is the first example of a large (>106M) HVC complex that is linked to stellar feedback occurring in a dwarf spiral galaxy.

955

, , , , and

Differential Ks-band luminosity functions (LFs) are presented for a complete sample of 1613 nearby bright galaxies segregated by visible morphology. The LF for late-type spirals follows a power law that rises toward low luminosities whereas the LFs for ellipticals, lenticulars, and bulge-dominated spirals are peaked and decline toward both higher and lower luminosities. Each morphological type (E, S0, S0/a–Sab, Sb–Sbc, Sc–Scd) contributes approximately equally to the overall Ks-band luminosity density of galaxies in the local universe. Type averaged bulge/disk ratios are used to subtract the disk component leading to the prediction that the Ks-band LF for bulges is bimodal with ellipticals dominating the high luminosity peak, comprising 60% of the bulge luminosity density in the local universe with the remaining 40% contributed by lenticulars and the bulges of spirals. Overall, bulges contribute 30% of the galaxy luminosity density at Ks in the local universe with spiral disks making up the remainder. If bulge luminosities indicate central black hole (BH) masses, then our results predict that the BH mass function is also bimodal.

970

, , , , and

We present a new technique for determining the quantity and composition of dust in astrophysical environments using <6 keV X-rays. We argue that high-resolution X-ray spectra as enabled by the Chandra and XMM-Newton gratings should be considered a powerful and viable new resource for delving into a relatively unexplored regime for directly determining dust properties: composition, quantity, and distribution. We present initial cross section measurements of astrophysically likely iron-based dust candidates taken at the Lawrence Berkeley National Laboratory Advanced Light Source synchrotron beamline, as an illustrative tool for the formulation of our technique for determining the quantity and composition of interstellar dust with X-rays. (Cross sections for the materials presented here will be made available for astrophysical modeling in the near future.) Focused at the 700 eV Fe Liii and Lii photoelectric edges, we discuss a technique for modeling dust properties in the soft X-rays using L-edge data to complement K-edge X-ray absorption fine structure analysis techniques discussed by Lee & Ravel. The paper is intended to be a techniques paper of interest and useful to both condensed matter experimentalists and astrophysicists. For the experimentalists, we offer a new prescription for normalizing relatively low signal-to-noise ratio L-edge cross section measurements. For astrophysics interests, we discuss the use of X-ray absorption spectra for determining dust composition in cold and ionized astrophysical environments and a new method for determining species-specific gas and dust ratios. Possible astrophysical applications of interest, including relevance to Sagittarius A*, are offered. Prospects for improving on this work in future X-ray missions with higher throughput and spectral resolution are also presented in the context of spectral resolution goals for gratings and calorimeters, for proposed and planned missions such as Astro-H and the International X-ray Observatory.

980

, , , , and

We present a regularized maximum likelihood weak-lensing reconstruction of the Deep Lens Survey F2 field (4 deg2). High signal-to-noise ratio peaks in our lensing significance map appear to be associated with possible projected filamentary structures. The largest apparent structure extends for over a degree in the field and has contributions from known optical clusters at three redshifts (z ∼ 0.3, 0.43, 0.5). Noise in weak-lensing reconstructions is known to potentially cause "false positives"; we use Monte Carlo techniques to estimate the contamination in our sample, and find that 10%–25% of the peaks are expected to be false detections. For significant lensing peaks, we estimate the total signal-to-noise ratio of detection using a method that accounts for pixel-to-pixel correlations in our reconstruction. We also report the detection of a candidate relative underdensity in the F2 field with a total signal-to-noise ratio of ∼5.5.

989

, , , , , , , , and

We present two exoplanets detected at Keck Observatory. HD 179079 is a G5 subgiant that hosts a hot Neptune planet with M sin i = 27.5 M in a 14.48 days, low-eccentricity orbit. The stellar reflex velocity induced by this planet has a semiamplitude of K = 6.6 m s−1. HD 73534 is a G5 subgiant with a Jupiter-like planet of M sin i = 1.1 MJup and K = 16 m s−1 in a nearly circular 4.85 yr orbit. Both stars are chromospherically inactive and metal-rich. We discuss a known, classical bias in measuring eccentricities for orbits with velocity semiamplitudes, K, comparable to the radial velocity uncertainties. For exoplanets with periods longer than 10 days, the observed exoplanet eccentricity distribution is nearly flat for large amplitude systems (K > 80 m s−1), but rises linearly toward low eccentricity for lower amplitude systems (K > 20 m s−1).

998

, , , , , , and

Solar energetic particle (SEP) events are traditionally classified as "impulsive" or "gradual." It is now widely accepted that in gradual SEP events, particles are accelerated at coronal mass ejection-driven (CME-driven) shocks. In many of these large SEP events, particle spectra exhibit double power law or exponential rollover features, with the break energy or rollover energy ordered as (Q/A)α, with Q being the ion charge in e and A the ion mass in units of proton mass mp. This Q/A dependence of the spectral breaks provides an opportunity to study the underlying acceleration mechanism. In this paper, we examine how the Q/A dependence may depend on shock geometry. Using the nonlinear guiding center theory, we show that α  ∼   1/5 for a quasi-perpendicular shock. Such a weak Q/A dependence is in contrast to the quasi-parallel shock case where α can reach 2. This difference in α reflects the difference of the underlying parallel and perpendicular diffusion coefficients κ|| and κ. We also examine the Q/A dependence of the break energy for the most general oblique shock case. Our analysis offers a possible way to remotely examine the geometry of a CME-driven shock when it is close to the Sun, where the acceleration of particle to high energies occurs.

1005

, , , and

We employ a high-resolution ΛCDM N-body simulation to present merger rate predictions for dark matter (DM) halos and investigate how common merger-related observables for galaxies—such as close pair counts, starburst counts, and the morphologically disturbed fraction—likely scale with luminosity, stellar mass, merger mass ratio, and redshift from z = 0 to z = 4. We investigate both rate at which subhalos first enter the virial radius of a larger halo (the "infall rate"), and the rate at which subhalos become destroyed, losing 90% of the mass they had at infall (the "destruction rate"). For both merger rate definitions, we provide a simple "universal" fitting formula that describes our derived merger rates for DM halos a function of dark halo mass, merger mass ratio, and redshift, and go on to predict galaxy merger rates using number density matching to associate halos with galaxies. For example, we find that the instantaneous (destruction) merger rate of m/M > 0.3 mass-ratio events into typical LfL* galaxies follows the simple relation dN/dt ≃ 0.03(1 + f) Gyr−1(1 + z)2.1. Despite the rapid increase in merger rate with redshift, only a small fraction of >0.4 L* high-redshift galaxies (∼3% at z = 2) should have experienced a major merger (m/M > 0.3) in the very recent past (t < 100 Myr). This suggests that short-lived, merger-induced bursts of star formation should not contribute significantly to the global star formation rate at early times, in agreement with several observational indications. In contrast, a fairly high fraction (∼20%) of those z = 2 galaxies should have experienced a morphologically transformative merger within a virial dynamical time (∼500 Myr at z = 2). We compare our results to observational merger rate estimates from both morphological indicators and pair-fraction-based determinations between z = 0and2 and show that they are consistent with our predictions. However, we emphasize that great care must be made in these comparisons because the predicted observables depend very sensitively on galaxy luminosity, redshift, overall mass ratio, and uncertain relaxation timescales for merger remnants. We show that the majority of bright galaxies at z = 3 should have undergone a major merger (>0.3) in the previous 700 Myr and conclude that mergers almost certainly play an important role in delivering baryons and influencing the kinematic properties of Lyman break galaxies (LBGs).

1016

, , , and

Many models for coronal loops have difficulty explaining the observed EUV brightness of the transition region, which is often significantly less than theoretical models predict. This discrepancy has been addressed by a variety of approaches including filling factors and time-dependent heating, with varying degrees of success. Here, we focus on an effect that has been ignored so far: the absorption of EUV light with wavelengths below 912 Å by the resonance continua of neutral hydrogen and helium. Such absorption is expected to occur in the low-lying transition region of hot, active region loops that is colocated with cool chromospheric features and called "moss" as a result of the reticulated appearance resulting from the absorption. We use cotemporal and cospatial spectroheliograms obtained with the Solar and Heliospheric Observatory/SUMER and Hinode/EIS of Fe xii 1242 Å, 195 Å, and 186.88 Å, and compare the density determination from the 186/195 Å line ratio to that resulting from the 195/1242 Å line ratio. We find that while coronal loops have compatible density values from these two line pairs, upper transition region moss has conflicting density determinations. This discrepancy can be resolved by taking into account significant absorption of 195 Å emission caused by the chromospheric inclusions in the moss. We find that the amount of absorption is generally of the order of a factor of 2. We compare to numerical models and show that the observed effect is well reproduced by three-dimensional radiative MHD models of the transition region and corona. We use STEREO A/B data of the same active region and find that increased angles between line of sight and local vertical cause additional absorption. Our determination of the amount of chromospheric absorption of TR emission can be used to better constrain coronal heating models.

1025

, , , and

We have conducted millimeter-wave observations of deuterated species of various carbon-chain molecules toward a low-mass star-forming region, L1527, which shows extraordinary richness of carbon-chain molecules in a vicinity of the protostar (Warm Carbon Chain Chemistry; WCCC). We have detected the spectral lines of l-C3D, C4D, C4HD, DC3N, DC5N, and c-C3HD, where l-C3D and C4HD are detected for the first time in space. The deuterium fractionation ratios are found to be moderate (2% to 7%), although they tend to be higher than those in the starless core, TMC-1. The upper limit to the [CH2DOH]/[CH3OH] ratio is also as low as 3%. Therefore, high deuterium fractionation ratios reported for hot corino sources are not seen in L1527. The observed ratios mean that the depletion of CO onto dust grains had not proceeded far in L1527, compared to the hot corino case. This would be consistent with a short timescale of the starless core phase, as suggested for the possible origin of WCCC.

1036

, , , and

We have re-analyzed continuum and recombination lines radio data available in the literature in order to derive the luminosity function (LF) of Galactic H ii regions. The study is performed by considering the first and fourth Galactic quadrants independently. We estimate the completeness level of the sample in the fourth quadrant at 5 Jy, and the one in the first quadrant at 2 Jy. We show that the two samples (fourth or first quadrant) include, as well as giant and supergiant H ii regions, a significant number of subgiant sources. The LF is obtained, in each Galactic quadrant, with a generalized Schmidt's estimator using an effective volume derived from the observed spatial distribution of the considered H ii regions. The re-analysis also takes advantage of recently published ancillary absorption data allowing to solve the distance ambiguity for several objects. A single power-law fit to the LFs retrieves a slope equal to −2.23 ± 0.07 (fourth quadrant) and to −1.85 ± 0.11 (first quadrant). We also find marginal evidence of a luminosity break at Lknee = 1023.45 erg s−1 Hz−1 for the LF in the fourth quadrant. We convert radio luminosities into equivalent Hα and Lyman continuum luminosities to facilitate comparisons with extragalactic studies. We obtain an average total H ii regions Lyman continuum luminosity of 0.89 ± 0.23 × 1053 s−1, corresponding to 30% of the total ionizing luminosity of the Galaxy.

1048

, , and

Umbral dots (UDs) were observed in a stable sunspot in NOAA 10944 by the Hinode Solar Optical Telescope on 2007 March 1. The observation program consisted of blue continuum images and spectropolarimetric profiles of Fe i 630 nm line. An automatic detection algorithm for UDs was applied to the 2 hr continuous blue continuum images, and using the obtained data, the lifetime, size, and proper motion of UDs were calculated. The magnetic structure of the sunspot was derived through the inversion of the spectropolarimetric profiles. We calculated the correlations between UD's parameters (size, lifetime, occurrence rate, proper motion) and magnetic fields (field strength, inclination, azimuth), and obtained the following results. (1) Both the lifetime and size of UDs are almost constant regardless of the magnetic field strength at their emergence site. (2) The speed of UDs increases as the field inclination angle at their emergence site gets larger. (3) The direction of movement of UDs is nearly parallel to the direction of the horizontal component of magnetic field in the region with strongly inclined field, while UDs in the region with weakly inclined field show virtually no proper motion. Our results describe the basic properties of magnetoconvection in sunspots. We will discuss our results in comparison to recent magnetohydrodynamic simulations by Schüssler & Vögler and Rempel et al.

1058

, , , , , , , and

We investigate the formation of the stellar halos of four simulated disk galaxies using high-resolution, cosmological SPH + N-body simulations. These simulations include a self-consistent treatment of all the major physical processes involved in galaxy formation. The simulated galaxies presented here each have a total mass of ∼1012M, but span a range of merger histories. These simulations allow us to study the competing importance of in situ star formation (stars formed in the primary galaxy) and accretion of stars from subhalos in the building of stellar halos in a ΛCDM universe. All four simulated galaxies are surrounded by a stellar halo, whose inner regions (r < 20 kpc) contain both accreted stars, and an in situ stellar population. The outer regions of the galaxies' halos were assembled through pure accretion and disruption of satellites. Most of the in situ halo stars formed at high redshift out of smoothly accreted cold gas in the inner 1 kpc of the galaxies' potential wells, possibly as part of their primordial disks. These stars were displaced from their central locations into the halos through a succession of major mergers. We find that the two galaxies with recently quiescent merger histories have a higher fraction of in situ stars (∼20%–50%) in their inner halos than the two galaxies with many recent mergers (∼5%–10% in situ fraction). Observational studies concentrating on stellar populations in the inner halo of the Milky Way will be the most affected by the presence of in situ stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic feature of galaxy formation.

1068

, , and

We have studied the sensitivity of s-process nucleosynthesis in massive stars to ±2σ variations in the rates of the triple-α and 12C(α, γ)16O reactions. We simulated the evolution of massive stars from H burning through Fe-core collapse, followed by a supernova explosion. We found that the production factors of s-process nuclides between 58Fe and 96Zr change strongly with changes in the He burning reaction rates; using the Lodders solar abundances rather than those of Anders and Grevesse reduces s-process nucleosynthesis; later burning phases beyond core He burning and shell C burning have a significant effect on post-explosive production factors. We also discuss the implications of the uncertainties in the helium burning rates for evidence of a new primary neutron capture process (LEPP) in massive stars.

1078

and

With the development of one-dimensional stellar evolution codes including rotation and the increasing number of observational data for stars of various evolutionary stages, it becomes more and more possible to follow the evolution of the rotation profile and angular momentum distribution in stars. In this context, understanding the interplay between rotation and convection in the very extended envelopes of giant stars is very important considering that all low- and intermediate-mass stars become red giants after the central hydrogen burning phase. In this paper, we analyze the interplay between rotation and convection in the envelope of red giant stars using three-dimensional numerical experiments. We make use of the Anelastic Spherical Harmonics code to simulate the inner 50% of the envelope of a low-mass star on the red giant branch. We discuss the organization and dynamics of convection, and put a special emphasis on the distribution of angular momentum in such a rotating extended envelope. To do so, we explore two directions of the parameter space, namely, the bulk rotation rate and the Reynolds number with a series of four simulations. We find that turbulent convection in red giant stars is dynamically rich, and that it is particularly sensitive to the rotation rate of the star. Reynolds stresses and meridional circulation establish various differential rotation profiles (either cylindrical or shellular) depending on the convective Rossby number of the simulations, but they all agree that the radial shear is large. Temperature fluctuations are found to be large and in the slowly rotating cases, a dominant ℓ = 1 temperature dipole influences the convective motions. Both baroclinic effects and turbulent advection are strong in all cases and mostly oppose one another.

1098

, , , , , and

We present deep Gemini GMOS optical spectroscopy of nine luminous quasars at redshifts z ∼ 0.5, drawn from the Sloan Digital Sky Survey type 2 quasar sample. Our targets were selected to have high intrinsic luminosities (MV < −26 mag) as indicated by the [O iii] λ5007 Å emission-line luminosity (L[O iii]). Our sample has a median black hole mass of ∼108.8M inferred assuming the local MBH–σ* relation and a median Eddington ratio of ∼0.7, using stellar velocity dispersions σ* measured from the G band. We estimate the contamination of the stellar continuum from scattered quasar light based on the strength of broad Hβ, and provide an empirical calibration of the contamination as a function of L[O iii]; the scattered-light fraction is ∼30% of L5100 for objects with L[O iii] = 109.5L. Population synthesis indicates that young poststarburst populations (<0.1 Gyr) are prevalent in luminous type 2 quasars, in addition to a relatively old population (>1 Gyr) which dominates the stellar mass. Broad emission complexes around He ii λ4686 Å with luminosities up to 108.3L are unambiguously detected in three out of the nine targets, indicative of Wolf–Rayet (WR) populations. Population synthesis shows that ∼5 Myr poststarburst populations contribute substantially to the luminosities (>50% of L5100) of all three objects with WR detections. We find two objects with double cores and four with close companions. Our results may suggest that luminous type 2 quasars trace an early stage of galaxy interaction, perhaps responsible for both the quasar and the starburst activity.

1118

, , , and

We measured the titanium (Ti) isotope composition, i.e., 50Ti/47Ti, 48Ti/47Ti, and 46Ti/47Ti, in five calcium-rich–aluminum-rich refractory inclusions (CAIs) from the oxidized CV3 chondrite Allende and in two CAIs from the reduced CV3 chondrite Efremovka. Our data indicate that CAIs are enriched in 50Ti/47Ti and 46Ti/47Ti and are slightly depleted in 48Ti/47Ti compared to normal Ti defined by ordinary chondrites, eucrites, ureilites, mesosiderites, Earth, Moon, and Mars. Some CAIs have an additional 50Ti excess of ∼8ε relative to bulk carbonaceous chondrites, which are enriched in 50Ti by ∼2ε relative to terrestrial values, leading to a total excess of ∼10ε. This additional 50Ti excess is correlated with nucleosynthetic anomalies found in 62Ni and 96Zr, all indicating an origin from a neutron-rich stellar source. Bulk carbonaceous chondrites show a similar trend, however, the extent of the anomalies is either less than or similar to the smallest anomalies seen in CAIs. Mass balance calculations suggest that bulk Allende Ti possibly consists of a mixture of at least two Ti components, anomalous Ti located in CAIs and a normal component possibly for matrix and chondrules. This argues for a heterogeneous distribution of Ti isotopes in the solar system. The finding that anomalous Ti is concentrated in CAIs suggests that CAIs formed in a specific region of the solar system and were, after their formation, not homogeneously redistributed within the solar system. Combining the CAI data with improved model predictions for early solar system irradiation effects indicates that a local production scenario for the relatively short lived radionuclides can be excluded, because the production of, e.g., 10Be, 26Al, and 41Ca, would result in a significant collateral shift in Ti isotopes, which is not seen in the measured data.

1127

, , , , , , , , and

We present subarcsecond resolution mid-infrared (mid-IR) photometry in the wavelength range from 8 to 20 μm of 18 Seyfert galaxies, reporting high spatial resolution nuclear fluxes for the entire sample. We construct spectral energy distributions (SEDs) that the active galactic nucleus (AGN) dominates, relatively uncontaminated by starlight, adding near-IR measurements from the literature at similar angular resolution. We find that the IR SEDs of intermediate-type Seyferts are flatter and present higher 10 to 18 μm ratios than those of Seyfert 2 galaxies. We fit the individual SEDs with clumpy dusty torus models using the in-house-developed BayesClumpy tool. We find that the clumpy models reproduce the high spatial resolution measurements. Regardless of the Seyfert type, even with high spatial resolution data, near- to mid-IR SED fitting poorly constrains the radial extent of the torus. For the Seyfert 2 galaxies, we find that edge-on geometries are more probable than face-on views, with a number of clouds along equatorial rays of N0 = 5–15. The 10 μm silicate feature is generally modeled in shallow absorption. For the intermediate-type Seyferts, N0 and the inclination angle of the torus are lower than those of the Seyfert 2 nuclei, with the silicate feature appearing in weak emission or absent. The columns of material responsible for the X-ray absorption are larger than those inferred from the model fits for most of the galaxies, which is consistent with X-ray absorbing gas being located within the dust sublimation radius, whereas the mid-IR flux arises from an area farther from the accretion disk. The fits yield both the bolometric luminosity of the intrinsic AGN and the torus-integrated luminosity, from which we derive the reprocessing efficiency of the torus. In the models, the outer radial extent of the torus scales with the AGN luminosity, and we find the tori to be confined to scales less than 5 pc.

1150

, , , and

The solar tachocline at the bottom of the convection zone is an important region for the dynamics of the Sun and the solar dynamo. In this region, the sound speed inferred by global helioseismology exhibits a bump of approximately 0.4% relative to the standard solar model. Global helioseismology does not provide any information on possible latitudinal variations or asymmetries between the northern and southern hemisphere. Here, we develop a time–distance helioseismology technique, including surface- and deep-focusing measurement schemes and a combination of both, for two-dimensional tomographic imaging of the solar tachocline that infers radial and latitudinal variations in the sound speed. We test the technique using artificial solar oscillation data obtained from numerical simulations. The technique successfully recovers major features of the simplified tachocline models. The technique is then applied to SOHO/MDI medium-ℓ data and provides for the first time a full two-dimensional sound-speed perturbation image of the solar tachocline. The one-dimensional radial profile obtained by latitudinal averaging of the image is in good agreement with the previous global helioseismology result. It is found that the amplitude of the sound-speed perturbation at the tachocline varies with latitude, but it is not clear whether this is in part or fully an effect of instrumental distortion. Our initial results demonstrate that time–distance helioseismology can be used to probe the deep interior structure of the Sun, including the solar tachocline.

1157

, , , , , , , , , et al

Recent observations have revealed that some Type Ia supernovae exhibit narrow, time-variable Na i D absorption features. The origin of the absorbing material is controversial, but it may suggest the presence of circumstellar gas in the progenitor system prior to the explosion, with significant implications for the nature of the supernova (SN) progenitors. We present the third detection of such variable absorption, based on six epochs of high-resolution spectroscopy of the Type Ia supernova SN 2007le from the Keck I Telescope and the Hobby–Eberly Telescope. The data span a time frame of approximately three months, from 5 days before maximum light to 90 days after maximum. We find that one component of the Na i D absorption lines strengthened significantly with time, indicating a total column density increase of ∼2.5 × 1012 cm−2. The data limit the typical timescale for the variability to be more than 2 days but less than 10 days. The changes appear to be most prominent after maximum light rather than at earlier times when the ultraviolet flux from the SN peaks. As with SN 2006X, we detect no change in the Ca ii H and K absorption lines over the same time period, rendering line-of-sight effects improbable and suggesting a circumstellar origin for the absorbing material. Unlike the previous two supernovae exhibiting variable absorption, SN 2007le is not highly reddened (EBV = 0.27 mag), also pointing toward circumstellar rather than interstellar absorption. Photoionization calculations show that the data are consistent with a dense (107 cm−3) cloud or clouds of gas located ∼0.1 pc (3 × 1017 cm) from the explosion. These results broadly support the single-degenerate scenario previously proposed to explain the variable absorption, with mass loss from a nondegenerate companion star responsible for providing the circumstellar gas. We also present possible evidence for narrow Hα emission associated with the SN, which will require deep imaging and spectroscopy at late times to confirm.

1171

and

The existence of a shallow decay phase in the early X-ray afterglows of gamma-ray bursts is a common feature. Here we investigate the possibility that this is connected to the formation of a highly magnetized millisecond pulsar, pumping energy into the fireball on timescales longer than the prompt emission. In this scenario, the nascent neutron star could undergo a secular bar-mode instability, leading to gravitational wave losses which would affect the neutron star spin-down. In this case, nearby gamma-ray bursts with isotropic energies of the order of 1050 ergs would produce a detectable gravitational wave signal emitted in association with an observed X-ray light-curve plateau, over relatively long timescales of minutes to about an hour. The peak amplitude of the gravitational wave signal would be delayed with respect to the gamma-ray burst trigger, offering gravitational wave interferometers such as the advanced LIGO and Virgo the challenging possibility of catching its signature on the fly.

1179

and

LS I +61°303 is a γ-ray binary with periodic radio outbursts coincident with the orbital period of P = 26.5 days. The origin of the radio emission is unclear, it could be due either to a jet, as in microquasars, or to the shock boundary between the Be star and a possible pulsar wind. We here analyze the radio spectral index over 6.7 years from Green Bank Interferometer data at 2.2 GHz and 8.3 GHz. We find two new characteristics in the radio emission. The first characteristic is that the periodic outbursts indeed consist of two consecutive outbursts; the first outburst is optically thick, whereas the second outburst is optically thin. The spectrum of LS I +61°303 is well reproduced by the shock-in-jet model commonly used in the context of microquasars and active galactic nuclei: the optically thin spectrum is due to shocks caused by relativistic plasma ("transient jet") traveling through a preexisting much slower steady flow ("steady jet"). This steady flow is responsible for the preceding optically thick spectrum. The second characteristic we find is that the observed spectral evolution, from optically thick to optically thin emission, occurs twice during the orbital period. We observed this occurrence at the orbital phase of the main 26.5 d outburst and also at an earlier phase, shifted by ΔΦ∼ 0.3 (i.e., almost 8 d before). We show that this result qualitatively and quantitatively agrees with the two-peak accretion/ejection model proposed in the past for LS I +61°303. We conclude that the radio emission in LS I +61°303 originates from a jet and suggest that the variable TeV emission comes from the usual Compton losses expected as an important by-product in the shock-in-jet theory.

1190

and

We performed a series of high-resolution (up to 10243) direct numerical simulations of hydro and magnetohydrodynamic (MHD) turbulence. Our simulations correspond to the "strong" MHD turbulence regime that cannot be treated perturbatively. We found that for simulations with normal viscosity the slopes for energy spectra of MHD are similar to ones in hydro, although slightly more shallower. However, for simulations with hyperviscosity the slopes were very different, for instance, the slopes for hydro simulations showed a pronounced and well defined bottleneck effect, while the MHD slopes were relatively much less affected. We believe that this is indicative of MHD strong turbulence being less local than the Kolmogorov turbulence. This calls for revision of MHD strong turbulence models that assume local "as-in-hydro case" cascading. Nonlocality of MHD turbulence casts doubt on numerical determination of the slopes with currently available (5123–10243) numerical resolutions, including simulations with normal viscosity. We also measure various so-called alignment effects and discuss their influence on the turbulent cascade.

1199

, , , and

The dynamical state of galaxy groups at intermediate redshifts can provide information about the growth of structure in the universe. We examine three goodness-of-fit tests, the Anderson–Darling (A–D), Kolmogorov, and χ2 tests, in order to determine which statistical tool is best able to distinguish between groups that are relaxed and those that are dynamically complex. We perform Monte Carlo simulations of these three tests and show that the χ2 test is profoundly unreliable for groups with fewer than 30 members. Power studies of the Kolmogorov and A–D tests are conducted to test their robustness for various sample sizes. We then apply these tests to a sample of the second Canadian Network for Observational Cosmology Redshift Survey (CNOC2) galaxy groups and find that the A–D test is far more reliable and powerful at detecting real departures from an underlying Gaussian distribution than the more commonly used χ2 and Kolmogorov tests. We use this statistic to classify a sample of the CNOC2 groups and find that 34 of 106 groups are inconsistent with an underlying Gaussian velocity distribution, and thus do not appear relaxed. In addition, we compute velocity dispersion profiles (VDPs) for all groups with more than 20 members and compare the overall features of the Gaussian and non-Gaussian groups, finding that the VDPs of the non-Gaussian groups are distinct from those classified as Gaussian.

1211

and

We perform time-resolved spectroscopy on the prompt emission in gamma-ray bursts (GRBs) and identify a thermal, photospheric component peaking at a temperature of a few hundreds keV. This peak does not necessarily coincide with the broad-band (keV–GeV) power peak. We show that this thermal component exhibits a characteristic temporal behavior. We study a sample of 56 long bursts, all strong enough to allow time-resolved spectroscopy. We analyze the evolution of both the temperature and flux of the thermal component in 49 individual time-resolved pulses, for which the temporal coverage is sufficient, and find that the temperature is nearly constant during the first few seconds, after which it decays as a power law with a sample-averaged index of −0.68. The thermal flux first rises with an averaged power-law index of 0.63 after which it decays with an averaged index of −2. The break times are the same to within errors. We find that the ratio of the observed to the emergent thermal flux typically exhibits a monotoneous power-law increase during the entire pulse as well as during complex bursts. Thermal photons carry a significant fraction (∼30% to more than 50%) of the prompt emission energy (in the observed 25–1900 keV energy band), thereby significantly contributing to the high radiative efficiency. Finally, we show here that the thermal emission can be used to study the properties of the photosphere, hence the physical parameters of the GRB fireball.

1230

, , and

We have re-analyzed the NRAO VLA Sky Survey (NVSS) data to derive rotation measures (RMs) toward 37,543 polarized radio sources. The resulting catalog of RM values covers the sky area north of declination −40° with an average density of more than one RM per square degree. We present an image of the median RM over 82% of the sky with a resolution of 8° and a typical error of ±1–2 rad m−2. The image shows large-scale structures in RM that extend to very high Galactic latitudes. A simple analysis of the RM structure at high Galactic latitudes is used to derive properties of the Galactic halo magnetic field in the solar neighborhood. We find the component of the local field perpendicular to the plane (the z-component) equal to +0.30 μG for z < 0 and −0.14 μG for z>0. The reversal of sign across the Galactic plane is consistent with a quadrupole field geometry for the poloidal component of the halo field. The halo magnetic field component parallel to the disk is also found to be antisymmetric and generally consistent with a toroidal field, with strength +0.83 μG for z < 0 and −0.39 μG for z>0. We have identified five regions of the sky where the foreground median RM is consistently less than 1 rad m−2 over several degrees. These holes in the foreground RM will be useful for future studies of possible small-scale fluctuations in cosmic magnetic field structures. In addition to allowing measurement of RMs toward polarized sources, the new analysis of the NVSS data removes the effects of bandwidth depolarization for |RM| ≳ 100 rad m−2 inherent in the original NVSS source catalog. This new catalog of RMs and polarized flux densities is available online, and will be a valuable resource for further studies of the Galactic magnetic field and magnetoionic medium, and extragalactic magnetic fields.

1237

We define and analyze the photometric orbit (PhO) of an extrasolar planet observed in reflected light. In our definition, the PhO is a Keplerian entity with six parameters: semimajor axis, eccentricity, mean anomaly at some particular time, argument of periastron, inclination angle, and effective radius, which is the square root of the geometric albedo times the planetary radius. Preliminarily, we assume a Lambertian phase function. We study in detail the case of short-period giant planets (SPGPs) and observational parameters relevant to the Kepler mission: 20 ppm photometry with normal errors, 6.5 hr cadence, and three-year duration. We define a relevant "planetary population of interest" in terms of probability distributions of the PhO parameters. We perform Monte Carlo experiments to estimate the ability to detect planets and to recover PhO parameters from light curves. We calibrate the completeness of a periodogram search technique, and find structure caused by degeneracy. We recover full orbital solutions from synthetic Kepler data sets and estimate the median errors in recovered PhO parameters. We treat in depth a case of a Jupiter body-double. For the stated assumptions, we find that Kepler should obtain orbital solutions for many of the 100–760 SPGP that Jenkins & Doyle estimate Kepler will discover. Because most or all of these discoveries will be followed up by ground-based radial velocity observations, the estimates of inclination angle from the PhO may enable the calculation of true companion masses: Kepler photometry may break the "msin i" degeneracy. PhO observations may be difficult. There is uncertainty about how low the albedos of SPGPs actually are, about their phase functions, and about a possible noise floor due to systematic errors from instrumental and stellar sources. Nevertheless, simple detection of SPGPs in reflected light should be robust in the regime of Kepler photometry, and estimates of all six orbital parameters may be feasible in at least a subset of cases.

1250

, , , and

We compare assembly of dark matter (DM) halos with and without baryons from identical initial conditions, within the context of cosmological evolution in the ΛCDM WMAP3 Universe (baryons+DM, hereafter BDM model, and pure DM, PDM model). In representative PDM and BDM models, we find that baryons contribute decisively to the evolution of the central region, leading to an isothermal DM cusp, and thereafter to a flat DM density core—the result of heating by dynamical friction of the DM+baryon substructure during a quiescent evolution epoch. This process ablates the cold gas from an embedded disk, cutting the star formation rate by a factor of 10, and heats up the spheroidal gas and stellar components, triggering their expansion. The substructure is more resilient to the tidal disruption in the presence of baryons. The disk which formed from inside-out as gas-dominated is transformed into an intermediate Hubble type by z ∼ 2 and to an early type by z ∼ 0.5, based on its gas contents and spheroidal-to-disk stellar mass ratio. We find that only a relatively small ∼20% fraction of DM particles in PDM and BDM models are bound within the radius of maximal circular velocity in the halo, slightly less so within halo characteristic radii—most of the DM particles perform larger radial excursions. The DM particles are unbound to the cusp region. We also find that the fraction of baryons within the halo virial radius somewhat increases during the major mergers and decreases during the minor mergers. The net effect appears to be negligible—an apparent result of our choice of feedback from stellar evolution. Furthermore, we find that the DM halos are only partially relaxed beyond their virialization. While the substructure is being tidally disrupted, mixing of its debris in the halo is not efficient and becomes even less so with z. The phase–space correlations (streamers) formed after z ∼ 1 will survive largely to the present time—an important implication for embedded disk evolution.

1268

, , and

In the lenticular galaxy NGC 1023 a third population of globular clusters (GCs), called faint fuzzies (FFs), was discovered next to the blue and red GC populations by Larsen & Brodie. While these FFs have colors comparable to the red population, the new population is fainter, larger (Reff>7 pc) and, most importantly, shows clear signs of corotation with the galactic disk of NGC 1023. We present N-body simulations verifying the hypothesis that these disk-associated FFs are related to the young massive cluster complexes (CCs) observed by Bastian et al. in M51, who discovered a mass–radius relation for these CCs. Our models have an initial configuration based on the observations from M51 and are placed on various orbits in a galactic potential derived for NGC 1023. All computations end up with a stable object containing 10%–60% of the initial CC mass after an integration time of 5 Gyr. A conversion to visual magnitudes demonstrates that the resulting objects cover exactly the observed range for FFs. Moreover, the simulated objects show projected half-mass radii between 3.6 and 13.4 pc, in good agreement with the observed FF sizes. We conclude that objects like the young massive CCs in M51 are likely progenitors of the FFs observed in NGC 1023.

1275

, , , and

This paper addresses the challenge of understanding the typical star formation histories of red-sequence galaxies, using linestrength indices and mass-to-light ratios as complementary constraints on their stellar age distribution. We first construct simple parametric models of the star formation history that bracket a range of scenarios, and fit these models to the linestrength indices of low-redshift cluster red-sequence galaxies. For giant galaxies, we confirm the downsizing trend, i.e., the stellar populations are younger, on average, for lower σ galaxies. We find, however, that this trend flattens or reverses at σ ≲ 70 km s−1. We then compare predicted stellar mass-to-light ratios with dynamical mass-to-light ratios derived from the fundamental plane (FP), or by the SAURON group. For galaxies with σ ∼ 70 km s−1, models with a late "frosting" of young stars and models with exponential star formation histories have stellar mass-to-light ratios that are larger than observed dynamical mass-to-light ratios by factors of 1.7 and 1.4, respectively, and so are rejected. The single stellar population (SSP) model is consistent with the FP, and requires a modest amount of dark matter (between 20% and 30%) to account for the difference between stellar and dynamical mass-to-light ratios. A model in which star formation was "quenched" at intermediate ages is also consistent with the observations, although in this case less dark matter is required for low mass galaxies. We also find that the contribution of stellar populations to the "tilt" of the fundamental plane is highly dependent on the assumed star formation history: for the SSP model, the tilt of the FP is driven primarily by stellar-population effects. For a quenched model, two-thirds of the tilt is due to stellar populations and only one-third is due to dark matter or non-homology.

1297

, , , , , and

We discuss the results stemming from observations of the white-light and [Fe xiv] emission corona during the total eclipse of the Sun of 2008 August 1, in Mongolia (Altaj region) and in Russia (Akademgorodok, Novosibirsk, Siberia). Corresponding to the current extreme solar minimum, the white-light corona, visible up to 20 solar radii, was of a transient type with well pronounced helmet streamers situated above a chain of prominences at position angles 48°, 130°, 241°, and 322°. A variety of coronal holes, filled with a number of thin polar plumes, were seen around the poles. Furthering an original method of image processing, stars up to 12 mag, a Kreutz-group comet (C/2008 O1) and a coronal mass ejection (CME) were also detected, with the smallest resolvable structures being of, and at some places even less than, 1 arcsec. Differences, presumably motions, in the corona and prominences are seen even with the 19 minutes time difference between our sites. In addition to the high-resolution coronal images, which show the continuum corona (K-corona) that results from electron scattering of photospheric light, images of the overlapping green-emission-line (530.3 nm, [Fe xiv]) corona were obtained with the help of two narrow-passband filters (centered on the line itself and for the continuum in the vicinity of 529.1 nm, respectively), each with an FWHM of 0.15 nm. Through solar observations, on whose scheduling and details we consulted, with the Solar and Heliospheric Observatory, Hinode's XRT and SOT, Transition Region and Coronal Explorer, and STEREO, as well as Wilcox Solar Observatory and Solar and Heliospheric Observatory/Michelson Doppler Imager magnetograms, we set our eclipse observations in the context of the current unusually low and prolonged solar minimum.

1309

, , , and

We show that the period clustering of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), their X-ray luminosities, ages, and statistics can be explained with fall back disks with large initial specific angular momentum. The disk evolution models are developed by comparison to self-similar analytical models. The initial disk mass and angular momentum set the viscous timescale. An efficient torque, with (1 − ω2*) dependence on the fastness parameter ω*, leads to period clustering in the observed AXP–SGR period range under a wide range of initial conditions. The timescale t0 for the early evolution of the fall back disk, and the final stages of fall back disk evolution, when the disk becomes passive, are the crucial determinants of the evolution. The disk becomes passive at temperatures around 100 K, which provides a natural cutoff for the X-ray luminosity and defines the end of evolution in the observable AXP and SGR phase. This low value for the minimum temperature for active disk turbulence indicates that the fall back disks are active up to a large radius, ≳1012 cm. We find that transient AXPs and SGRs are likely to be older than their persistent cousins. A fall back disk with mass transfer rates corresponding to the low quiescent X-ray luminosities of the transient sources in early evolutionary phases would have a relatively lower initial mass, such that the mass-flow rate in the disk is not sufficient for the inner disk to penetrate into the light cylinder of the young neutron star, making mass accretion onto the neutron star impossible. The transient AXP phase therefore must start later. The model results imply that the transient AXP/SGRs, although older, are likely to be similar in number to persistent sources. This is because the X-ray luminosities of AXPs and SGRs are found to decrease faster at the end of their evolution, and the X-ray luminosities of transient AXP and SGRs in quiescence lie in the luminosity range of X-ray cutoff phase. Taking the range of quiescent X-ray luminosities of transient AXPs and SGRs ∼1033–1034 erg s−1, our simulations imply that the duration of the cutoff phase is comparable to the lifetime in the persistent phase for a large range of initial conditions.

1321

, , , and

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of 12C16O in galaxies at redshift z = 0 to z = 10. This prediction relies on a recently presented simulation of the molecular cold gas content in ∼3 × 107 evolving galaxies based on the Millennium Simulation. We combine this simulation with a model for the conversion between molecular mass and CO-line intensities, which incorporates the following mechanisms: (1) molecular gas is heated by the cosmic microwave background (CMB), starbursts (SBs), and active galactic nuclei (AGNs); (2) molecular clouds in dense or inclined galaxies can overlap; (3) compact gas can attain a smooth distribution in the densest part of disks; (4) CO luminosities scale with metallicity changes between galaxies; and (5) CO luminosities are always detected against the CMB. We analyze the relative importance of these effects and predict the cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies (i.e., no AGN, no massive SB) at high z ≳ 7 in CO emission will be dramatically hindered by the weak contrast against the CMB, in contradiction to earlier claims that CMB heating will ease the detection of high-redshift CO. The full simulation of extragalactic CO lines and the predicted CO-LFs at any redshift can be accessed online (http://s-cubed.physics.ox.ac.uk/, go to "S3-SAX") and they should be useful for the modeling of CO-line surveys with future telescopes, such as the Atacama Large Millimeter/submillimeter Array or the Square Kilometre Array.

1336

, , and

Based on photometric data of the central parts of eight globular clusters and one open cluster presented by An and his collaborators, we select red horizontal branch (RHB) stars in the (gr)0g0 diagram and make a statistical study of the distributions of their colors and absolute magnitudes in the SDSS ugriz system. Meanwhile, absolute magnitudes in the Johnson VRI system are calculated through the translation formulae between gri and VRI in the literature. The calibrations of absolute magnitude as functions of metallicity and age are established by linear regressions of the data. It is found that metallicity coefficients in these calibrations decrease, while age coefficients increase, from the blue u filter to the red z filter. The calibration of Mi = 0.06[Fe/H] + 0.040t + 0.03 has the smallest scatter of 0.04 mag, and thus i is the best filter in the ugriz system when RHB stars are used for distance indicators. The comparison of the MI calibration from our data with that from red clump stars indicates that the previous suggestion that the I filter is better than the V filter in distance determination may not be true because of its significant dependence on age.

1343

, , , and

The issue whether Moreton waves are flare-ignited or coronal mass ejection (CME)-driven, or a combination of both, is still a matter of debate. We develop an analytical model describing the evolution of a large-amplitude coronal wave emitted by the expansion of a circular source surface in order to mimic the evolution of a Moreton wave. The model results are confronted with observations of a strong Moreton wave observed in association with the X3.8/3B flare/CME event from 2005 January 17. Using different input parameters for the expansion of the source region, either derived from the real CME observations (assuming that the upward moving CME drives the wave), or synthetically generated scenarios (expanding flare region, lateral expansion of the CME flanks), we calculate the kinematics of the associated Moreton wave signature. Those model input parameters are determined which fit the observed Moreton wave kinematics best. Using the measured kinematics of the upward moving CME as the model input, we are not able to reproduce the observed Moreton wave kinematics. The observations of the Moreton wave can be reproduced only by applying a strong and impulsive acceleration for the source region expansion acting in a piston mechanism scenario. Based on these results we propose that the expansion of the flaring region or the lateral expansion of the CME flanks is more likely the driver of the Moreton wave than the upward moving CME front.

1353

, , , , , , , , , et al

We present the first results from a high sampling rate, multimonth reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from telescopes around the world. The primary goal of this campaign was to obtain either new or improved Hβ reverberation lag measurements for several relatively low luminosity active galactic nuclei (AGNs). We feature results for NGC 4051 here because, until now, this object has been a significant outlier from AGN scaling relationships, e.g., it was previously a ∼2–3σ outlier on the relationship between the broad-line region (BLR) radius and the optical continuum luminosity—the RBLRL relationship. Our new measurements of the lag time between variations in the continuum and Hβ emission line made from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of RBLR = 1.87+0.54−0.50 light days and black hole mass of MBH = (1.73+0.55−0.52) × 106M. This radius is consistent with that expected from the RBLRL relationship, based on the present luminosity of NGC 4051 and the most current calibration of the relation by Bentz et al.. We also present a preliminary look at velocity-resolved Hβ light curves and time delay measurements, although we are unable to reconstruct an unambiguous velocity-resolved reverberation signal.

1367

and

We present a uniform X-ray spectral analysis of eight type-1 active galactic nuclei that have been previously observed with relativistically broadened iron emission lines. Utilizing data from the XMM-Newton European Photon Imaging Camera (EPIC-pn) we carefully model the spectral continuum, taking complex intrinsic absorption and emission into account. We then proceed to model the broad Fe Kα feature in each source with two different accretion disk emission line codes, as well as a self-consistent, ionized accretion disk spectrum convolved with relativistic smearing from the inner disk. Comparing the results, we show that relativistic blurring of the disk emission is required to explain the spectrum in most sources, even when one models the full reflection spectrum from the photoionized disk.

1387

, , , and

We investigate the secular evolution of the orbital semimajor axis and eccentricity due to mass transfer in eccentric binaries, allowing for both mass and angular momentum loss from the system. Adopting a delta function mass transfer rate at the periastron of the binary orbit, we find that, depending on the initial binary properties at the onset of mass transfer, the orbital semimajor axis and eccentricity can either increase or decrease at a rate linearly proportional to the magnitude of the mass transfer rate at periastron. The range of initial binary mass ratios and eccentricities that lead to increasing orbital semimajor axes and eccentricities broadens with increasing degrees of mass loss from the system and narrows with increasing orbital angular momentum loss from the binary. Comparison with tidal evolution timescales shows that the usual assumption of rapid circularization at the onset of mass transfer in eccentric binaries is not justified, irrespective of the degree of systemic mass and angular momentum loss. This work extends our previous results for conservative mass transfer in eccentric binaries and can be incorporated into binary evolution and population synthesis codes to model non-conservative mass transfer in eccentric binaries.

1393

, , , , , , , , , et al

We use very deep near-infrared (NIR) imaging data obtained in MOIRCS Deep Survey (MODS) to investigate the evolution of the galaxy stellar mass function back to z ∼ 3. The MODS data reach J = 24.2, H = 23.1, and K = 23.1 (5σ, Vega magnitude) over 103 arcmin2 (wide) and J = 25.1, H = 23.7, and K = 24.1 over 28 arcmin2 (deep) in the GOODS-North region. The wide and very deep NIR data allow us to measure the number density of galaxies down to low stellar mass (109–1010M) even at high redshift with high statistical accuracy. The normalization of the mass function decreases with redshift, and the integrated stellar mass density becomes ∼8%–18% of the local value at z ∼ 2 and ∼4%–9% at z ∼ 3, which are consistent with results of previous studies in general fields. Furthermore, we found that the low-mass slope becomes steeper with redshift from α ∼ −1.3 at z ∼ 1 to α ∼ −1.6 at z ∼ 3 and that the evolution of the number density of low-mass (109–1010M) galaxies is weaker than that of M* (∼1011M) galaxies. This indicates that the contribution of low-mass galaxies to the total stellar mass density has been significant at high redshift. The steepening of the low-mass slope with redshift is an opposite trend expected from the stellar mass dependence of the specific star formation rate reported in previous studies. The present result suggests that the hierarchical merging process overwhelmed the effect of the stellar mass growth by star formation and was very important for the stellar mass assembly of these galaxies at 1 ≲ z ≲ 3.

1413

, , and

We examine the radius evolution of close in giant planets with a planet evolution model that couples the orbital–tidal and thermal evolution. For 45 transiting systems, we compute a large grid of cooling/contraction paths forward in time, starting from a large phase space of initial semimajor axes and eccentricities. Given observational constraints at the current time for a given planet (semimajor axis, eccentricity, and system age), we find possible evolutionary paths that match these constraints, and compare the calculated radii to observations. We find that tidal evolution has two effects. First, planets start their evolution at larger semimajor axis, allowing them to contract more efficiently at earlier times. Second, tidal heating can significantly inflate the radius when the orbit is being circularized, but this effect on the radius is short-lived thereafter. Often circularization of the orbit is proceeded by a long period while the semimajor axis slowly decreases. Some systems with previously unexplained large radii that we can reproduce with our coupled model are HAT-P-7, HAT-P-9, WASP-10, and XO-4. This increases the number of planets for which we can match the radius from 24 (of 45) to as many as 35 for our standard case, but for some of these systems we are required to be viewing them at a special time around the era of current radius inflation. This is a concern for the viability of tidal inflation as a general mechanism to explain most inflated radii. Also, large initial eccentricities would have to be common. We also investigate the evolution of models that have a floor on the eccentricity, as may be due to a perturber. In this scenario, we match the extremely large radius of WASP-12b. This work may cast some doubt on our ability to accurately determine the interior heavy element enrichment of normal, noninflated close in planets, because of our dearth of knowledge about these planets' previous orbital–tidal histories. Finally, we find that the end state of most close in planetary systems is disruption of the planet as it moves ever closer to its parent star.

1428

and

Recently, Hennebelle and Chabrier derived an analytical theory for the mass spectrum of non-self-gravitating clumps associated with overdensities in molecular clouds and for the initial mass function (IMF) of gravitationally bound prestellar cores, as produced by the turbulent collapse of the cloud. In this companion paper, we examine the effects of the nonisothermality of the flow, of the turbulence forcing and of local fluctuation of the velocity dispersion, on the mass function. In particular, we investigate the influence of a polytropic equation of state (eos) and of the effective adiabatic exponent γ and find that it has a drastic influence on the low-mass part of the IMF. We also consider a barotropic eos (i.e., a piecewise polytropic eos) that mimics the thermal behavior of the molecular gas and compare the prediction of our theory with the results of numerical simulations and with the observationally derived IMF, for cloud parameters which satisfy Larson-type relations. We find that for clouds whose density is, at all scales, almost an order of magnitude larger than the density inferred for the CO clumps in the Galaxy, a good agreement is obtained between the theory and the observed IMF, suggesting that star formation preferentially occurs in high-density environments. We derive an analytical expression for the IMF which generalizes the expression previously obtained for the isothermal case. This easy-to-implement analytical IMF should serve as a template to compare observational or numerical results with the theory.

1443

and

We report the observations of the magnetohydrodynamic (MHD) waves propagating along magnetic flux tubes in the solar photosphere. We identified 20 isolated strong peaks (8 peaks for pores and 12 peaks for intergranular magnetic structure) in the power spectra of the line-of-sight (LOS) magnetic flux, the LOS velocity, and the intensity for 14 different magnetic concentrations. The observation is performed with the spectro-polarimeter of the Solar Optical Telescope aboard the Hinode satellite. The oscillation periods are located in 3–6 minutes for the pores and in 4–9 minutes for the intergranular magnetic elements. These peaks correspond to the magnetic, the velocity, and the intensity fluctuation in time domain with root-mean-square amplitudes of 4–17 G (0.3%–1.2%), 0.03–0.12 km s−1, and 0.1%–1%, respectively. Phase differences between the LOS magnetic flux (ϕB), the LOS velocity (ϕv), the intensities of the line core (ϕI,core), and the continuum intensity (ϕI,cont) have striking concentrations at around −90° for ϕB − ϕv and ϕv − ϕI,core, around 180° for ϕI,core − ϕB, and around 10° for ϕI,core − ϕI,cont. Here, for example, ϕB − ϕv ∼ −90° means that the velocity leads the magnetic field by a quarter of cycle. The observed phase relation between the magnetic and the photometric intensity fluctuations would not be consistent with that caused by the opacity effect, if the magnetic field strength decreases with height along the oblique LOS. We suggest that the observed fluctuations are due to longitudinal (sausage-mode) and/or transverse (kink-mode) MHD waves. The observed phase relation between the fluctuations in the magnetic flux and the velocity is consistent with the superposition of the ascending wave and the descending wave reflected at chromosphere/corona boundary (standing wave). Even with such reflected waves, the residual upward Poynting flux is estimated to be 2.7 × 106 erg cm−2 s−1 for a case of the kink wave. Seismology of the magnetic flux tubes is possible to obtain various physical parameters from the observed period and amplitude of the oscillations.

1458

, , , and

We present a study of the classification of z ∼ 1 extremely red objects (EROs), using a combination of Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS), Spitzer Infrared Array Camera (IRAC), and ground-based images of the COSMOS field. Our sample includes ∼5300 EROs with iKs ⩾ 2.45 (AB, equivalently IKs = 4 in Vega) and Ks ⩽ 21.1 (AB). For EROs in our sample, we compute, using the ACS F814W images, their concentration, asymmetry, as well as their Gini coefficient and the second moment of the brightest 20% of their light. Using those morphology parameters and the Spitzer IRAC [3.6] − [8.0] color, the spectral energy distribution (SED) fitting method, we classify EROs into two classes: old galaxies (OGs) and young, dusty starburst galaxies (DGs). We found that the fraction of OGs and DGs in our sample is similar, about 48% of EROs in our sample are OGs, and 52% of them are DGs. To reduce the redundancy of these three different classification methods, we performed a principal component analysis on the measurements of EROs, and find that morphology parameters and SEDs are efficient in segregating OGs and DGs. The [3.6] − [8.0] color, which depends on reddening, redshift, and photometric accuracy, is difficult to separate EROs around the discriminating line between starburst and elliptical. We investigate the dependence of the fraction of EROs on their observational properties, and the results suggest that DGs become increasingly important at fainter magnitudes, redder colors, and higher redshifts. The clustering of the entire EROs, DGs, and OGs was estimated by calculating their correlation function, and we find that the clustering of EROs is much stronger than that of full K-limited samples of galaxies; the clustering amplitude of OGs is a factor of ∼2 larger than that of DGs.

1472

and

Column densities for H i, Al iii, Si iv, C iv, and O vi toward 109 stars and 30 extragalactic objects have been assembled to study the extensions of these species away from the Galactic plane into the Galactic halo. H i and Al iii mostly trace the warm neutral and warm ionized medium, respectively, while Si iv, C iv, and O vi trace a combination of warm photoionized and collisionally ionized plasmas. The much larger object sample compared to previous studies allows us to consider and correct for the effects of the sample bias that has affected earlier but smaller surveys of the gas distributions. We find that Si iv and C iv have similar exponential scale heights of 3.2(+1.0, −0.6) and 3.6(+1.0, −0.8) kpc. The scale height of O vi is marginally smaller with h = 2.6 ± 0.6 kpc. The transition temperature gas is ∼3 times more extended than the warm ionized medium traced by Al iii with h = 0.90(+0.62, −0.33) kpc and ∼12 times more extended than the warm neutral medium traced by H i with h = 0.24 ± 0.06 kpc. There is a factor of 2 decrease in the dispersion of the log of the column density ratios for transition temperature gas for lines of sight in the Galactic disk compared to extragalactic lines of sight through the entire halo. The observations are compared to the predictions of the various models for the production of the transition temperature gas in the halo. The appendix presents a revision to the electron scale height of Gaensler et al.'s 2008 study based on electron dispersion measures.

1490

, , , , and

Collisions between dust aggregates are the key to understand the formation of planetesimals because the collision inevitably takes place in protoplanetary disks. To clarify whether or not dust aggregates can grow through their mutual collisions at relatively high velocities, we carry out more than 4000 runs of three-dimensional numerical simulations of collisions between icy equal-mass clusters formed under ballistic particle–cluster aggregation (BPCA) as well as those of ballistic cluster–cluster aggregation, including offset collisions with various values of the impact parameter. Since our BPCA clusters have a fractal dimension of 3 and a relatively compact structure, their results enable us to determine the criteria for growth and disruption of compressed aggregates at their collisions in protoplanetary disks. The results show that ice dust aggregates are able to grow at collisions with velocities up to 50 m s−1, in spite of their initial structures and impact parameters. We also find that the mass of ejecta relative to the total mass of colliding aggregates decreases with increasing the size of the aggregates. These results demonstrate the feasibility of growth and survival for dust aggregates through their mutual collisions with relatively high velocities in protoplanetary disks.

1502

It is widely believed that lenticular (S0) galaxies were initially spirals from which the gas has been removed by interactions with hot cluster gas, or by ram pressure stripping of cool gas from spirals that are orbiting within rich clusters of galaxies. However, problems with this interpretation are that (1) some lenticulars, such as NGC 3115, are isolated field galaxies rather than cluster members. (2) The distribution of flattening values of S0 galaxies in clusters, in groups, and in the field are statistically indistinguishable. This is surprising because one might have expected most of the progenitors of field S0 galaxies to have been flattened late-type galaxies, whereas lenticulars in clusters are thought to have mostly been derived from bulge-dominated early-type galaxies. (3) It should be hardest for ram pressure to strip massive luminous galaxies with deep potential wells. However, no statistically significant differences are seen between the luminosity distributions of early-type Shapley–Ames galaxies in clusters, groups, and in the field. (4) Finally both ram pressure stripping and evaporation by hot intracluster gas would be most efficient in rich clusters. However, the small number of available data in the Shapley–Ames sample appears to show no statistically significant differences between the relative frequencies of dust-poor S01 and dust-rich S03 galaxies in clusters, groups, and in the field. It is tentatively concluded that ram pressure stripping and heating by intracluster gas, may not be the only evolutionary channels that lead to the formation of lenticular galaxies. It is speculated that gas starvation, or gas ejection by active nuclei, may have played a major role in the formation of a significant fraction of all S0 galaxies.

1507

, , , , , , , , , and

We have used Spitzer/Infrared Array Camera (IRAC) to conduct a photometric monitoring program of the IC1396A dark globule in order to study the mid-IR (3.6–8 μm) variability of the heavily embedded young stellar objects (YSOs) present in that area. We obtained light curves covering a 14 day timespan with a twice daily cadence for 69 YSOs, and continuous light curves with approximately 12 s cadence over 7 hr for 38 YSOs. Typical accuracies for our relative photometry were 1%–2% for the long timespan data and a few millimagnitude, corresponding to less than 0.5%, for the 7 hr continuous "staring-mode" data. More than half of the YSOs showed detectable variability, with amplitudes from ∼0.05 mag to ∼0.2 mag. About 30% of the YSOs showed quasi-sinusoidal light-curve shapes with apparent periods from 5 to 12 days and light-curve amplitudes approximately independent of wavelength over the IRAC bandpasses. We have constructed models which simulate the time-dependent spectral energy distributions of Class I and II YSOs in order to attempt to explain these light curves. Based on these models, the apparently periodic light curves are best explained by YSO models where one or two high-latitude photospheric spots heat the inner wall of the circumstellar disk, and where we view the disk at fairly large inclination angle. Disk inhomogeneities, such as increasing the height where the accretion funnel flows to the stellar hot spot, enhances the light-curve modulations. The other YSOs in our sample show a range of light-curve shapes, some of which are probably due to varying accretion rate or disk shadowing events. One star, IC1396A-47, shows a 3.5 hr periodic light curve; this object may be a PMS Delta Scuti star.

1530

, , , , , , , and

We investigate a peculiar feature at the hottest, blue end of the horizontal branch of Galactic globular cluster ω Centauri, using the high-precision and nearly complete catalog that has been constructed from a survey taken with the Advanced Camera for Survey on board the Hubble Space Telescope, that covers the inner 10 × 10 arcmin. It is a densely populated clump of stars with an almost vertical structure in the F435W-(F435W−F625W) plane, that we termed "blue clump." A comparison with theoretical models leads to the conclusion that this feature must necessarily harbor either hot flasher stars or canonical He-rich stars—progeny of the blue main sequence (MS) subpopulation observed in this cluster—or a mixture of both types, plus possibly a component from the normal-He population hosted by the cluster. A strong constraint coming from theory is that the mass of the objects in the "blue clump" has to be very finely tuned, with a spread of at most only ∼0.03 M. By comparing observed and theoretical star counts along both the H- and He-burning stages we find that at least 15% of the expected He-rich horizontal branch stars are missing from the color–magnitude diagram. This missing population could be the progeny of red giants that failed to ignite central He-burning and have produced He-core white dwarfs (WDs). Our conclusion supports the scenario recently suggested by Calamida et al. for explaining the observed ratio of WDs to MS stars in ω Centauri.

1536

, , , and

Using a 2.5D time-dependent numerical code we recently developed, we solve the full compressible Navier–Stokes equations to determine the structure of the boundary layer (BL) between the white dwarf (WD) and the accretion disk in nonmagnetic cataclysmic variable systems. In this preliminary work, our numerical approach does not include radiation. In the energy equation, we either take the dissipation function (Φ) into account or we assume that the energy dissipated by viscous processes is instantly radiated away (Φ = 0). For a slowly rotating nonmagnetized accreting WD, the accretion disk extends all the way to the stellar surface. There, the matter impacts and spreads toward the poles as new matter continuously piles up behind it. We carry out numerical simulations for different values of the alpha-viscosity parameter (α), corresponding to different mass accretion rates. In the high viscosity cases (α = 0.1), the spreading BL sets off a gravity wave in the surface matter. The accretion flow moves supersonically over the cusp making it susceptible to the rapid development of gravity wave and/or Kelvin–Helmholtz shearing instabilities. This BL is optically thick and extends more than 30° to either side of the disk plane after only 3/4 of a Keplerian rotation period (tK = 19 s). In the low viscosity cases (α = 0.001), the spreading BL does not set off gravity waves and it is optically thin.

1553

, , and

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self-consistently combined Fokker–Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a ∼10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a nonthermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

1567

and

We present results of an analysis of the J, H, and Ks Two Micron All Sky Survey (2MASS) images of 139 spiral edge-on galaxies selected from the Revised Flat Galaxies Catalog. The basic structural parameters scale length (h), scale height (z0), and central surface brightness of the stellar disks (μ0) are determined for all selected galaxies in the near-infrared (NIR) bands. The mean relative ratios of the scale heights of the thin stellar disks in the J:H:Ks bands are 1.16:1.08:1.00, respectively. Comparing the scale heights obtained from the NIR bands for the same objects, we estimate the scale heights of the thin stellar disks corrected for the internal extinction. We find that the extinction-corrected scale height is, on average, 11% smaller than that in the K band. Using the extinction-corrected structural parameters, we find that the dark-to-luminous mass ratio is, on average, 1.3 for the galaxies in our sample within the framework of a simplified galactic model. The relative thicknesses of the stellar disks z0/h correlates with their face-on central surface brightnesses obtained from the 2MASS images. We also find that the scale height of the stellar disks shows no systematic growth with radius in most of our galaxies.

1575

, , and

We examine the case for quark-novae (QNe) as possible sources for the reionization and early metal enrichment of the universe. QNe are predicted to arise from the explosive collapse (and conversion) of sufficiently massive neutron stars into quark stars (QSs). A QN can occur over a range of timescales following the supernova (SN) event. For QNe that arise days to weeks after the SNe, we show that dual shock that arises as the QN ejecta encounter the SN ejecta can produce enough photons to reionize hydrogen in most of the intergalactic medium (IGM) by z ∼ 6. Such events can explain the large optical depth τe ∼ 0.1 as measured by WMAP, if the clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a normal initial mass function for the oldest stars can be reconciled with a large optical depth as well as the mean metallicity of the early IGM post reionization. We find that QN also make a contribution to r-process element abundances for atomic numbers A ⩾ 130. We predict that the main cosmological signatures of QNe are the gamma-ray bursts that announce their birth. These will be clustered at redshifts in the range z ∼ 7–8 in our model.

1584

, , , , , , , and

We present single-dish 350 μm dust continuum polarimetry as well as HCN and HCO+J = 4 → 3 rotational emission spectra obtained on NGC 1333 IRAS 4. The polarimetry indicates a uniform field morphology over a 20'' radius from the peak continuum flux of IRAS 4A, in agreement with models of magnetically supported cloud collapse. The field morphology around IRAS 4B appears to be quite distinct, however, with indications of depolarization observed toward the peak flux of this source. Inverse P Cygni profiles are observed in the HCN J = 4 → 3 line spectra toward IRAS 4A, providing a clear indication of infall gas motions. Taken together, the evidence gathered here appears to support the scenario that IRAS 4A is a cloud core in a critical state of support against gravitational collapse.

1593

, , , and

We present a follow-up analysis of the unique magnetic luminosity-variable carbon-atmosphere white dwarf SDSS J142625.71+575218.3. This includes the results of some 106.4 hr of integrated light photometry which have revealed, among other things, the presence of a new periodicity at 319.720 s which is not harmonically related to the dominant oscillation (417.707 s) previously known in that star. Using our photometry and available spectroscopy, we consider the suggestion made by Montgomery et al. that the luminosity variations in SDSS J142625.71+575218.3 may not be caused by pulsational instabilities, but rather by photometric activity in a carbon-transferring analog of AM CVn. This includes a detailed search for possible radial velocity variations due to rapid orbital motion on the basis of Multiple Mirror Telescope spectroscopy. At the end of the exercise, we unequivocally rule out the interacting binary hypothesis and conclude instead that, indeed, the luminosity variations are caused by g-mode pulsations as in other pulsating white dwarfs. This is in line with the preferred possibility put forward by Montgomery et al.

1604

, , , and

We analyze measured proton and electron temperatures in the high-speed solar wind in order to calculate the separate rates of heat deposition for protons and electrons. When comparing with other regions of the heliosphere, the fast solar wind has the lowest density and the least frequent Coulomb collisions. This makes the fast wind an optimal testing ground for studies of collisionless kinetic processes associated with the dissipation of plasma turbulence. Data from the Helios and Ulysses plasma instruments were collected to determine mean radial trends in the temperatures and the electron heat conduction flux between 0.29 and 5.4 AU. The derived heating rates apply specifically for these mean plasma properties and not for the full range of measured values around the mean. We found that the protons receive about 60% of the total plasma heating in the inner heliosphere, and that this fraction increases to approximately 80% by the orbit of Jupiter. A major factor affecting the uncertainty in this fraction is the uncertainty in the measured radial gradient of the electron heat conduction flux. The empirically derived partitioning of heat between protons and electrons is in rough agreement with theoretical predictions from a model of linear Vlasov wave damping. For a modeled power spectrum consisting only of Alfvénic fluctuations, the best agreement was found for a distribution of wavenumber vectors that evolves toward isotropy as distance increases.

1615

, , , and

We present the results of a high angular resolution Very Large Array (VLA) Class I 44 GHz and Class II 6.7 GHz CH3OH maser survey of a sample of ∼20 massive young stellar object (MYSO) outflow candidates selected on the basis of extended 4.5 μm emission in Spitzer Galactic Legacy Infrared Mid-Plane Survey Extraordinaire images. These 4.5 μm selected candidates are referred to as extended green objects (EGOs), for the common coding of this band as green in three-color Infrared Array Camera images. The detection rate of 6.7 GHz Class II CH3OH masers, which are associated exclusively with massive YSOs, toward EGOs is ≳64%—nearly double the detection rate of surveys using other MYSO selection criteria. The detection rate of Class I 44 GHz CH3OH masers, which trace molecular outflows, is ∼89% toward EGOs associated with 6.7 GHz CH3OH masers. The two types of CH3OH masers exhibit different spatial distributions: 6.7 GHz masers are centrally concentrated and usually coincide with 24 μm emission, while 44 GHz masers are widely distributed and generally trace diffuse 4.5 μm features. We also present results of a complementary James Clerk Maxwell Telescope (JCMT) single-pointing molecular line survey of EGOs in the outflow tracers HCO+(3–2) and SiO(5–4). The HCO+ line profiles and high SiO detection rate (90%) are indicative of the presence of active outflows. No 44 GHz continuum emission is detected at the 5 mJy beam−1 (5σ) level toward 95% of EGOs surveyed, excluding bright ultracompact H ii regions as powering sources for the 4.5 μm outflows. The results of our surveys constitute strong evidence that EGOs are young, massive YSOs, with active outflows, presumably powered by ongoing accretion.

1648

and

The black hole (BH) X-ray transient XTE J1550 − 564 has undergone a strong outburst in 1998 and two relativistic X-ray jets have been detected years later with the Chandra X-ray observatory; the eastern jet was found previously to have decelerated after its first detection. Here, we report a full analysis of the evolution of the western jet; significant deceleration is also detected in the western side. Our analysis indicates that there is a cavity outside the central source and the jets first traveled with constant velocity and then were slowed down by the interactions between the jets and the interstellar medium. The best-fitted radius of the cavity is ∼0.31 pc on the eastern side and ∼0.44 pc on the western side, and the densities also show asymmetry, of ∼0.034 cm−3 on the east to ∼0.12 cm−3 on the west. The best-fitted magnetic fields on both sides are ∼0.5 mG. Similar analysis is also applied to another microquasar system, H 1743 − 322, and a large-scale low-density region is also found. Based on these results and the comparison with other microquasar systems, we suggest a generic scenario for microquasar jets, classifying the observed jets into three main categories, with different jet morphologies (and sizes) corresponding to different scales of vacuous environments surrounding them. We also suggest that either continuous jets or accretion disk winds, or both may be responsible for creating these cavities. Therefore X-ray jets from microquasars provide us with a promising method of probing the environment of accreting BHs.

1662

and

The differential migration of two planets due to planet–disk interaction can result in capture into the 2:1 eccentricity-type mean-motion resonances. Both the sequence of 2:1 eccentricity resonances that the system is driven through by continued migration and the possibility of a subsequent capture into the 4:2 inclination resonances are sensitive to the migration rate within the range expected for type II migration due to planet–disk interaction. If the migration rate is fast, the resonant pair can evolve into a family of 2:1 eccentricity resonances different from those found by Lee. This new family has outer orbital eccentricity e2 ≳ 0.4–0.5, asymmetric librations of both eccentricity resonance variables, and orbits that intersect if they are exactly coplanar. Although this family exists for an inner-to-outer planet mass ratio m1/m2 ≳ 0.2, it is possible to evolve into this family by fast migration only for m1/m2 ≳ 2. Thommes and Lissauer have found that a capture into the 4:2 inclination resonances is possible only for m1/m2 ≲ 2. We show that this capture is also possible for m1/m2 ≳ 2 if the migration rate is slightly slower than that adopted by Thommes and Lissauer. There is significant theoretical uncertainty in both the sign and the magnitude of the net effect of planet–disk interaction on the orbital eccentricity of a planet. If the eccentricity is damped on a timescale comparable to or shorter than the migration timescale, e2 may not be able to reach the values needed to enter either the new 2:1 eccentricity resonances or the 4:2 inclination resonances. Thus, if future observations of extrasolar planetary systems were to reveal certain combinations of mass ratio and resonant configuration, they would place a constraint on the strength of eccentricity damping during migration, as well as on the rate of the migration itself.

1673

, , , , , , and

The Rossi X-ray Timing Explorer has observed five outbursts from the transient 2.5 ms accretion-powered pulsar SAX J1808.4−3658 during 1998–2008. We present a pulse timing study of the most recent outburst and compare it with the previous timing solutions. The spin frequency of the source continues to decrease at a rate of (−5.5 ± 1.2) × 10−18 Hz s−1, which is consistent with the previously determined spin derivative. The spin down occurs mostly during quiescence, and is most likely due to the magnetic dipole torque from a B = 1.5 × 108 G dipolar field at the neutron star surface. We also find that the 2 hr binary orbital period is increasing at a rate of (3.80 ± 0.06) × 10−12 s s−1, also consistent with previous measurements. It remains uncertain whether this orbital change reflects secular evolution or short-term variability.

1679

and

X-ray monitoring observations were performed with the Swift observatory of the ultraluminous X-ray sources Holmberg IX X-1, NGC 5408 X-1, and NGC 4395 X-2 and also of the nuclear X-ray source in NGC 4395. Holmberg IX X-1 remains in the hard X-ray spectral state as its flux varies by a factor of 7 up to a (isotropic) luminosity of 2.8 × 1040 erg s−1. This behavior may suggest an unusually massive compact object. We find excess power at periods near 60 days and 28 days in the X-ray emission from Holmberg IX X-1. Additional monitoring is required to test the significance of these signals. NGC 5408 X-1 and NGC 4395 X-2 appear to remain in the soft spectral state found by Chandra and XMM with little variation in spectral hardness even as the luminosity changes by a factor of 9. We found an outburst from the nuclear source in NGC 4395 reaching an X-ray luminosity of 9 × 1040 erg s−1, several times higher than any previously reported.

L101

Disk galaxies are in hydrostatic equilibrium along their vertical axis. The pressure allowing for this configuration consists of thermal, turbulent, magnetic, and cosmic-ray components. For the Milky Way the thermal pressure contributes ∼10% of the total pressure near the plane, with this fraction dropping toward higher altitudes. Out of the rest, magnetic fields contribute ∼1/3 of the pressure to distances of ∼3 kpc above the disk plane. In this Letter, we attempt to extrapolate these local values to high-redshift, rapidly accreting, rapidly star-forming disk galaxies and study the effect of the extra pressure sources on the accretion of gas onto the galaxies. In particular, magnetic field tension may convert a smooth cold-flow accretion to clumpy, irregular star formation regions and rates. The infalling gas accumulates on the edge of the magnetic fields, supported by magnetic tension. When the mass of the infalling gas exceeds some threshold mass, its gravitational force cannot be balanced by magnetic tension anymore, and it falls toward the disk's plane, rapidly making stars. Simplified estimations of this threshold mass are consistent with clumpy star formation observed in SINS, UDF, GOODS, and GEMS surveys. We discuss the shortcomings of pure hydrodynamic codes in simulating the accretion of cold flows into galaxies, and emphasize the need for magnetohydrodynamic simulations.

L105

, , , , , , , , , and

Images of Titan's clouds, possible over the past 10 years, indicate primarily discrete convective methane clouds near the south and north poles and an immense stratiform cloud, likely composed of ethane, around the north pole. Here we present spectral images from Cassini's Visual Mapping Infrared Spectrometer that reveal the increasing presence of clouds in Titan's tropical atmosphere. Radiative transfer analyses indicate similarities between summer polar and tropical methane clouds. Like their southern counterparts, tropical clouds consist of particles exceeding 5 μm. They display discrete structures suggestive of convective cumuli. They prevail at a specific latitude band between 8°–20° S, indicative of a circulation origin and the beginning of a circulation turnover. Yet, unlike the high latitude clouds that often reach 45 km altitude, these discrete tropical clouds, so far, remain capped to altitudes below 26 km. Such low convective clouds are consistent with the highly stable atmospheric conditions measured at the Huygens landing site. Their characteristics suggest that Titan's tropical atmosphere has a dry climate unlike the south polar atmosphere, and despite the numerous washes that carve the tropical landscape.

L110

, , , , , , , , , et al

We describe and present initial results of a weak lensing survey of nearby (z ≲ 0.1) galaxy clusters in the Sloan Digital Sky Survey (SDSS). In this first study, galaxy clusters are selected from the SDSS spectroscopic galaxy cluster catalogs of Miller et al. and Berlind et al. We report a total of seven individual low-redshift cluster weak lensing measurements that include A2048, A1767, A2244, A1066, A2199, and two clusters specifically identified with the C4 algorithm. Our program of weak lensing of nearby galaxy clusters in the SDSS will eventually reach ∼200 clusters, making it the largest weak lensing survey of individual galaxy clusters to date.

L114

, , , , , , , , , and

We analyze a 100 ks Chandra observation of the powerful radio galaxy, 4C 60.07 at z = 3.79. We identify extended X-ray emission with LX ∼ 1045 erg s−1 across a ∼90 kpc region around the radio galaxy. The energetics of this X-ray halo and its morphological similarity to the radio emission from the galaxy suggest that it arises from inverse Compton (IC) scattering, by relativistic electrons in the radio jets, of cosmic microwave background photons and potentially far-infrared photons from the dusty starbursts around this galaxy. The X-ray emission has a similar extent and morphology to the Lyα halo around the galaxy, suggesting that it may be ionizing this halo. Indeed, we find that the GHz-radio and X-ray and Lyα luminosities of the halo around 4C 60.07 are identical to those of 4C 41.17 (also at z = 3.8) implying that these three components are linked by a single physical process. This is only the second example of highly extended IC emission known at z>3, but it underlines the potential importance of IC emission in the formation of the most massive galaxies at high redshifts. In addition, we detect two X-ray luminous active galactic nuclei (AGNs) within ∼30 kpc of the radio galaxy. These two companion AGNs imply that the radio and starburst activity in the radio galaxy is triggered through multiple mergers of massive progenitors on a short timescale, ≲100 Myr. These discoveries demonstrate the wealth of information which sensitive X-ray observations can yield into the formation of massive galaxies at high redshifts.

L119

We report on the discovery of the young, nearby, brown dwarf 2MASS J0041353−562112. The object has a spectral type of M7.5; it shows Li absorption and signatures of accretion, which implies that it still has a disk and suggests an age below 10 Myr. The space motion vector and position on the sky indicate that the brown dwarf is probably a member of the ∼20 Myr old Tuc-Hor association, or that it may be an ejected member of the ∼12 Myr old β Pic association; both would imply that 2MASS J0041353−562112 may in fact be older than 10 Myr. No accreting star or brown dwarf was previously known in these associations. Assuming an age of 10 Myr, the brown dwarf has a mass of about 30 MJup and is located at 35 pc distance. The newly discovered object is the closest accreting brown dwarf known. Its membership to an association older than 10 Myr implies that either disks in brown dwarfs can survive as long as in more massive stars, perhaps even longer, or that star formation in Tuc-Hor or β Pic occurred more recently than previously thought. The history and evolution of this object can provide new fundamental insight into the formation process of stars, brown dwarfs, and planets.

L124

, , , , , and

A new interstellar molecule, HSCN (thiocyanic acid), an energetic isomer of the well-known species HNCS, has been detected toward Sgr B2(N) with the Arizona Radio Observatory 12 m telescope. Eight rotational transitions in the Ka = 0 ladder were observed in the 2 mm and 3 mm bands. Five consecutive transitions in the 3 mm band are unblended, but three in the 2 mm band are partially masked by lines of other molecules. The peak intensity of all eight transitions are well described by a rotational temperature that is in very good agreement with that of many other molecules in this source. The line width and radial velocity of HSCN match closely with those of the ground state isomer HNCS (isothiocyanic acid), HNCO (isocyanic acid), and HOCN (cyanic acid); preliminary maps indicate that all four molecules are similarly distributed in Sgr B2. Although HSCN is calculated to lie over 3000 K higher in energy than HNCS, its column density of 1.3 × 1013 cm−2 in Sgr B2(N) is only three times lower than that of HNCS. The fractional abundances of HSCN and HNCS relative to H2 are 4.5 × 10−12 and 1.1 × 10−11. By analogy with the isomeric pair HCN and HNC, these two sulfur-bearing isomers are plausibly formed from a common cation precursor.

L128

, , , , , , , , , and

We report the detection of several molecular gas-phase and ice absorption features in three photometrically selected young stellar object (YSO) candidates in the central 280 pc of the Milky Way. Our spectra, obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, reveal gas-phase absorption from CO2 (15.0 μm), C2H2 (13.7 μm), and HCN (14.0 μm). We attribute this absorption to warm, dense gas in massive YSOs. We also detect strong and broad 15 μm CO2 ice absorption features, with a remarkable double-peaked structure. The prominent long-wavelength peak is due to CH3OH-rich ice grains, and is similar to those found in other known massive YSOs. Our IRS observations demonstrate the youth of these objects, and provide the first spectroscopic identification of massive YSOs in the Galactic Center.

L133

, , and

The presence of fine structures in sunspot vector magnetic fields has been confirmed from Hinode as well as other earlier observations. We studied 43 sunspots based on the data sets taken from ASP/DLSP, Hinode (SOT/SP), and SVM (USO). In this Letter, (1) we introduce the concept of signed shear angle (SSA) for sunspots and establish its importance for non-force-free fields. (2) We find that the sign of global α (force-free parameter) is well correlated with that of the global SSA and the photospheric chirality of sunspots. (3) Local α patches of opposite signs are present in the umbra of each sunspot. The amplitude of the spatial variation of local α in the umbra is typically of the order of the global α of the sunspot. (4) We find that the local α is distributed as alternately positive and negative filaments in the penumbra. The amplitude of azimuthal variation of the local α in the penumbra is approximately an order of magnitude larger than that in the umbra. The contributions of the local positive and negative currents and α in the penumbra cancel each other giving almost no contribution for their global values for the whole sunspot. (5) Arc-like structures (partial rings) with a sign opposite to that of the dominant sign of α of the umbral region are seen at the umbral–penumbral boundaries of some sunspots. (6) Most of the sunspots studied belong to the minimum epoch of the 23rd solar cycle and do not follow the so-called hemispheric helicity rule.

L138

, , and

We address the origin of the patchy dark and bright emission structure, known as "moss," observed by TRACE extreme ultraviolet observations of the solar disk. Here we propose an explanation based on turbulent, patchy heat conduction from the corona into the transition region. Computer simulations demonstrate that magnetic turbulence in coronal loops develops a flux rope structure with current sheets near the flux rope boundaries. Localized heating due to current sheet activity such as magnetic reconnection is followed by heat conduction along turbulent magnetic field lines. The field line trajectories tend to remain near the flux rope boundaries, resulting in selective heating of the plasma in the transition region. This can explain the network of bright regions in the observed moss morphology.

L143

, , , , , , , , , et al

We propose a physically motivated and self-consistent prescription for the modeling of transient neutron star low-mass X-ray binary (LMXB) properties, such as duty cycle (DC), outburst duration, and recurrence time. We apply this prescription to the population synthesis models of field LMXBs presented by Fragos et al., and compare the transient LMXB population to the Chandra X-ray survey of the two elliptical galaxies NGC 3379 and NGC 4278, which revealed several transient sources. We are able to exclude models with a constant DC for all transient systems, while models with a variable DC based on the properties of each system are consistent with the observed transient populations. We predict that the majority of the observed transient sources in these two galaxies are LMXBs with red giant donors. Finally, our comparison suggests that transient LMXBs are very rare in globular clusters (GCs), and thus the number of identified transient LMXBs may be used as a tracer of the relative contribution of field and GC LMXB populations.

L148

, , , and

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in active galactic nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1%. This lends support to the hypothesis that the universe is seeded by B fields that are generated in AGNs via this mechanism and subsequently injected into intergalactic space by the jet outflows.

L153

, , , , , , , , , et al

We present new results from BRAVA, a large-scale radial velocity survey of the Galactic bulge, using M giant stars selected from the Two Micron All Sky Survey catalog as targets for the Cerro Tololo Inter-American Observatory 4 m Hydra multi-object spectrograph. The purpose of this survey is to construct a new generation of self-consistent bar models that conform to these observations. We report the dynamics for fields at the edge of the Galactic bulge at latitudes b = −8° and compare to the dynamics at b = −4°. We find that the rotation curve V(r) is the same at b = −8° as at b = −4°. That is, the Galactic boxy bulge rotates cylindrically, as do boxy bulges of other galaxies. The summed line-of-sight velocity distribution at b = −8° is Gaussian, and the binned longitude–velocity plot shows no evidence for either a (disk) population with cold dynamics or for a (classical bulge) population with hot dynamics. The observed kinematics are well modeled by an edge-on N-body bar, in agreement with published structural evidence. Our kinematic observations indicate that the Galactic bulge is a prototypical product of secular evolution in galaxy disks, in contrast with stellar population results that are most easily understood if major mergers were the dominant formation process.

L158

, , , , , , and

A new and more comprehensive model of charge-exchange induced X-ray emission, due to ions precipitating into the Jovian atmosphere near the poles, has been used to analyze spectral observations made by the Chandra X-ray Observatory. The model includes for the first time carbon ions, in addition to the oxygen and sulfur ions previously considered, in order to account for possible ion origins from both the solar wind and the Jovian magnetosphere. By comparing the model spectra with newly reprocessed Chandra observations, we conclude that carbon ion emission provides a negligible contribution, suggesting that solar wind ions are not responsible for the observed polar X-rays. In addition, results of the model fits to observations support the previously estimated seeding kinetic energies of the precipitating ions (∼0.7–2 MeV u−1), but infer a different relative sulfur-to-oxygen abundance ratio for these Chandra observations.

L163

and

Disk fragmentation resulting from the gravitational instability has been proposed as an efficient mechanism for forming giant planets. We use the planet Fomalhaut b, the triple-planetary system HR 8799, and the potential protoplanet associated with HL Tau to test the viability of this mechanism. We choose the above systems since they harbor planets with masses and orbital characteristics favored by the fragmentation mechanism. We do not claim that these planets must have formed as the result of fragmentation, rather the reverse: if planets can form from disk fragmentation, then these systems are consistent with what we should expect to see. We use the orbital characteristics of these recently discovered planets, along with a new technique to more accurately determine the disk cooling times, to place both lower and upper limits on the disk surface density—and thus mass—required to form these objects by disk fragmentation. Our cooling times are over an order of magnitude shorter than those of Rafikov, which makes disk fragmentation more feasible for these objects. We find that the required mass interior to the planet's orbital radius is ∼0.1 M for Fomalhaut b, the protoplanet orbiting HL Tau, and the outermost planet of HR 8799. The two inner planets of HR 8799 probably could not have formed in situ by disk fragmentation.

L168

, , and

We employ tomographic observations of a small region of plage to study the propagation of waves from the solar photosphere to the chromosphere using a Fourier phase-difference analysis. Our results show the expected vertical propagation for waves with periods of 3 minutes. Waves with 5 minute periods, i.e., above the acoustic cutoff period, are found to propagate only at the periphery of the plage, and only in the direction in which the field can be reasonably expected to expand. We conclude that field inclination is critically important in the leakage of p-mode oscillations from the photosphere into the chromosphere.

L172

, , , , , , and

We have identified in an acid resistant residue of the carbonaceous chondrite Murchison a large number (458) of highly refractory metal nuggets (RMNs) that once were most likely hosted by Ca,Al-rich inclusions (CAIs). While osmium isotopic ratios of two randomly selected particles rule out a presolar origin, the bulk chemistry of 88 particles with sizes in the submicron range determined by energy dispersive X-ray (EDX) spectroscopy shows striking agreement with predictions of single-phase equilibrium condensation calculations. Both chemical composition and morphology strongly favor a condensation origin. Particularly important is the presence of structurally incompatible elements in particles with a single-crystal structure, which also suggests the absence of secondary alteration. The metal particles represent the most pristine early solar system material found so far and allow estimation of the cooling rate of the gaseous environment from which the first solids formed by condensation. The resulting value of 0.5 K yr−1 is at least 4 orders of magnitude lower than the cooling rate of molten CAIs. It is thus possible, for the first time, to see through the complex structure of most CAIs and infer the thermal history of the gaseous reservoir from which their components formed by condensation.

L177

, , and

We report the discovery of three pulsars whose large dispersion measures (DMs) and angular proximity to Sgr A* indicate the existence of a Galactic center population of neutron stars. The relatively long periods (0.98–1.48 s) most likely reflect strong selection against short-period pulsars from radio-wave scattering at the observation frequency of 2 GHz used in our survey with the Green Bank Telescope. One object (PSR J1746−2850I) has a characteristic spindown age of only 13 kyr along with a high surface magnetic field ∼4 × 1013 G. It and a second object found in the same telescope pointing, PSR J1746−2850II (which has the highest known DM among pulsars), may have originated from recent star formation in the Arches or Quintuplet clusters given their angular locations. Along with a third object, PSR J1745−2910, and two similar high-dispersion, long-period pulsars reported by Johnston et al., the five objects found so far are 10–15 arcmin from Sgr A*, consistent with there being a large pulsar population in the Galactic center, most of whose members are undetectable in relatively low-frequency surveys because of pulse broadening from the same scattering volume that angularly broadens Sgr A* and OH/IR masers.

L182

, , and

Magnetic fields are dragged in from the interstellar medium during the gravitational collapse that forms star/disk systems. Consideration of mean field magnetohydrodynamics in these disks shows that magnetic effects produce sub-Keplerian rotation curves and truncate the inner disk. This Letter explores the ramifications of these predicted disk properties for the migration of extrasolar planets. Sub-Keplerian flow in gaseous disks drives a new migration mechanism for embedded planets and modifies the gap-opening processes for larger planets. This sub-Keplerian migration mechanism dominates over Type I migration for sufficiently small planets (mP ≲ 1 M) and/or close orbits (r ≲ 1 AU). Although the inclusion of sub-Keplerian torques shortens the total migration time by only a moderate amount, the mass accreted by migrating planetary cores is significantly reduced. Truncation of the inner disk edge (for typical system parameters) naturally explains final planetary orbits with periods P ∼ 4 days. Planets with shorter periods, P ∼ 2 days, can be explained by migration during FU-Orionis outbursts, when the mass accretion rate is high and the disk edge moves inward. Finally, the midplane density is greatly increased at the inner truncation point of the disk (the X-point); this enhancement, in conjunction with continuing flow of gas and solids through the region, supports the in situ formation of giant planets.

L187

, , , and

We present unambiguous evidence for a parsec scale wind in the broad-line radio galaxy 3C 382, the first radio-loud active galactic nucleus, with RL = log10(f5 GHz/f4400)>1, whereby an outflow has been measured with X-ray grating spectroscopy. A 118 ks Chandra grating (HETG) observation of 3C 382 has revealed the presence of several high ionization absorption lines in the soft X-ray band, from Fe, Ne, Mg, and Si. The absorption lines are blueshifted with respect to the systemic velocity of 3C 382 by −840 ± 60 km s−1 and are resolved by Chandra with a velocity width of σ = 340 ± 70 km s−1. The outflow appears to originate from a single zone of gas of column density NH = 1.3 × 1021 cm−2 and ionization parameter log(ξ/erg cm s−1) = 2.45. From the above measurements we calculate that the outflow is observed on parsec scales, within the likely range from 10to1000 pc, i.e., consistent with an origin in the narrow-line region.