A Study of CO Emission in High Redshift QSOs Using the Owens Valley Millimeter Array

Searches for CO emission in high-redshift objects have traditionally suffered from the accuracy of optically-derived redshifts due to lack of bandwidth in correlators at radio observatories. This problem has motivated the creation of the new COBRA continuum correlator, with 4 GHz available bandwidth, at the Owens Valley Radio Observatory Millimeter Array. Presented here are the first scientific results from COBRA. We report detections of redshifted CO(J=3-2) emission in the QSOs SMM J04135+10277 and VCV J140955.5+562827, as well as a probable detection in RX J0911.4+0551. At redshifts of z=2.846, z=2.585, and z=2.796, we find integrated CO flux densities of 5.4 Jy km/s, 2.4 Jy km/s, and 2.9 Jy km/s for SMM J04135+10277, VCV J140955.5+562827, and RX J0911.4+0551, respectively, over linewidths of Delta(V_{FWHM}) ~ 350 km/s. These measurements, when corrected for gravitational lensing, correspond to molecular gas masses of order M(H_2) ~ 10^{9.6-11.1} solar masses, and are consistent with previous CO observations of high-redshift QSOs. We also report 3-sigma upper limits on CO(3-2) emission in the QSO LBQS 0018-0220 of 1.3 Jy km/s. We do not detect significant 3mm continuum emission from any of the QSOs, with the exception of a tentative (3-sigma) detection in RX J0911.4+0551 of S_{3mm}=0.92 mJy/beam.Comment: 18 pages, 5 figures, 2 tables, accepted to ApJ. Changes made for version 2: citations added, 2 objects added to Table 2 and Figure

A CO Survey of Gravitationally Lensed Quasars with the IRAM Interferometer

We present the results of a CO survey of gravitationally lensed quasars, conducted with the Plateau de Bure Interferometer over the last three years. Among the 18 objects surveyed, one was detected in CO line emission, while six were detected in the continuum at 3mm and three in the continuum at 1mm. The low CO detection rate may at least in part be due to uncertainties in the redshifts derived from quasar broad emission lines. The detected CO source, the z=3.2 radio quiet quasar MG0751+2716, is quite strong in the CO(4-3) line and in the millimeter/submillimeter continuum, the latter being emission from cool dust. The integrated CO line flux is 5.96 +- 0.45 Jy.km/s, and the total molecular gas mass is estimated to be in the range M(H_2) = 1.6-3.1 X 10^9 solar masses.Comment: 5 pages, 2 figures, uses aa.cls and psfig.st

Molecular Gas in the Lensed Lyman Break Galaxy cB58

We have used the IRAM Plateau de Bure Interferometer to map CO(3-2) emission from the gravitationally lensed Lyman break galaxy MS1512-cB58. This is the first detection of a molecular emission line in any Lyman break system; its integrated intensity implies a total molecular gas mass of 6.6e9 Msun, while its width implies a dynamical mass of 1.0e10 csc^2i Msun (for a flat Lambda=0.7 cosmology). These estimates are in excellent concordance with nearly all parameters of the system measured at other wavelengths, and yield a consistent picture of past and future star formation with no obvious discrepancies requiring explanation by differential lensing. In particular, we find that the age and remaining lifetime of the current episode of star formation are likely to be similar; the surface densities of star formation and molecular gas mass are related by a Schmidt law; and the fraction of baryonic mass already converted into stars is sufficient to account for the observed enrichment of the interstellar medium to 0.4 Zsun. Barring substantial gas inflow or a major merger, the stars forming in the current episode will have mass and coevality at z=0 similar to those of a spiral bulge. Assuming cB58 is a typical Lyman break galaxy apart from its magnification, its global parameters suggest that the prescriptions for star formation used in some semi-analytic models of galaxy evolution require moderate revision, although the general prediction that gas mass fraction should increase with redshift is validated. [abridged]Comment: 41 pages, 6 figures, accepted by Ap

Radio Observations of Infrared Luminous High Redshift QSOs

We present Very Large Array (VLA) observations at 1.4 GHz and 5 GHz of a sample of 12 Quasi-stellar Objects (QSOs) at z = 3.99 to 4.46. The sources were selected as the brightest sources at 250 GHz from the recent survey of Omont et al. (2001). We detect seven sources at 1.4 GHz with flux densities, S_{1.4} > 50 microJy. These centimeter (cm) wavelength observations imply that the millimeter (mm) emission is most likely thermal dust emission. The radio-through-optical spectral energy distributions for these sources are within the broad range defined by lower redshift, lower optical luminosity QSOs. For two sources the radio continuum luminosities and morphologies indicate steep spectrum, radio loud emission from a jet-driven radio source. For the remaining 10 sources the 1.4 GHz flux densities, or limits, are consistent with those expected for active star forming galaxies. If the radio emission is powered by star formation in these systems, then the implied star formation rates are of order 1e3 M_solar/year. We discuss the angular sizes and spatial distributions of the radio emitting regions, and we consider briefly these results in the context of co-eval black hole and stellar bulge formation in galaxies.Comment: to appear in the A

Sensitive Radio Observations of High Redshift Dusty QSOs

We present sensitive radio continuum imaging at 1.4 GHz and 4.9 GHz of seven high redshift QSOs selected for having a 240 GHz continuum detection, which is thought to be thermal dust emission. We detect radio continuum emission from four of the sources: BRI 0952-0115, BR 1202-0725, LBQS 1230+1627B, and BRI 1335-0417. The radio source in BR 1202-0725 is resolved into two components, coincident with the double mm and CO sources. We compare the results at 1.4 GHz and 240 GHz to empirical and semi-analytic spectral models based on star forming galaxies at low redshift. The radio-to-submm spectral energy distribution for BR 1202-0725, LBQS 1230+1627B, and BRI 1335-0417 are consistent with that expected for a massive starburst galaxy, with implied massive star formation rates of order 1000 solar masses per year (without correcting for possible amplification by gravitational lensing). The radio-to-submm spectral energy distribution for BRI 0952-0115 suggests a low-luminosity radio jet source driven by the AGN.Comment: 12 pages, Latex emulateapj format, including 1 table and 3 figures. The Astrophysical Journal, to appear in the January 2000 issu

A Compact Starburst Core in the Dusty Lyman Break Galaxy Westphal-MD11

Using the IRAM Plateau de Bure Interferometer, we have searched for CO(3-2) emission from the dusty Lyman break galaxy Westphal-MD11 at z = 2.98. Our sensitive upper limit is surprisingly low relative to the system's 850 um flux density and implies a far-IR/CO luminosity ratio as elevated as those seen in local ultraluminous mergers. We conclude that the observed dust emission must originate in a compact structure radiating near its blackbody limit and that a relatively modest molecular gas reservoir must be fuelling an intense nuclear starburst (and/or deeply buried active nucleus) that may have been triggered by a major merger. In this regard, Westphal-MD11 contrasts strikingly with the lensed Lyman break galaxy MS1512-cB58, which is being observed apparently midway through an extended episode of more quiescent disk star formation.Comment: 5 pages, 1 figure (emulateapj), accepted by ApJ

Molecular line radiative transfer in protoplanetary disks: Monte Carlo simulations versus approximate methods

We analyze the line radiative transfer in protoplanetary disks using several approximate methods and a well-tested Accelerated Monte Carlo code. A low-mass flaring disk model with uniform as well as stratified molecular abundances is adopted. Radiative transfer in low and high rotational lines of CO, C18O, HCO+, DCO+, HCN, CS, and H2CO is simulated. The corresponding excitation temperatures, synthetic spectra, and channel maps are derived and compared to the results of the Monte Carlo calculations. A simple scheme that describes the conditions of the line excitation for a chosen molecular transition is elaborated. We find that the simple LTE approach can safely be applied for the low molecular transitions only, while it significantly overestimates the intensities of the upper lines. In contrast, the Full Escape Probability (FEP) approximation can safely be used for the upper transitions (J_{\rm up} \ga 3) but it is not appropriate for the lowest transitions because of the maser effect. In general, the molecular lines in protoplanetary disks are partly subthermally excited and require more sophisticated approximate line radiative transfer methods. We analyze a number of approximate methods, namely, LVG, VEP (Vertical Escape Probability) and VOR (Vertical One Ray) and discuss their algorithms in detail. In addition, two modifications to the canonical Monte Carlo algorithm that allow a significant speed up of the line radiative transfer modeling in rotating configurations by a factor of 10--50 are described.Comment: 47 pages, 12 figures, accepted for publication in Ap

Structure of W3(OH) from Very High Spectral Resolution Observations of 5 Centimeter OH Masers

Recent studies of methanol and ground-state OH masers at very high spectral resolution have shed new light on small-scale maser processes. The nearby source W3(OH), which contains numerous bright masers in several different transitions, provides an excellent laboratory for high spectral resolution techniques. We present a model of W3(OH) based on EVN observations of the rotationally-excited 6030 and 6035 MHz OH masers taken at 0.024 km/s spectral resolution. The 6.0 GHz masers are becoming brighter with time and show evidence for tangential proper motions. We confirm the existence of a region of magnetic field oriented toward the observer to the southeast and find another such region to the northeast in W3(OH), near the champagne flow. The 6.0 GHz masers trace the inner edge of a counterclockwise rotating torus feature. Masers at 6030 MHz are usually a factor of a few weaker than at 6035 MHz but trace the same material. Velocity gradients of nearby Zeeman components are much more closely correlated than in the ground state, likely due to the smaller spatial separation between Zeeman components. Hydroxyl maser peaks at very long baseline interferometric resolution appear to have structure on scales both smaller than that resolvable as well as on larger scales.Comment: 21 pages using emulateapj.cls including 16 figures and 2 tables, accepted to Ap

Generation of potential/surface density pairs in flat disks Power law distributions

We report a simple method to generate potential/surface density pairs in flat axially symmetric finite size disks. Potential/surface density pairs consist of a ``homogeneous'' pair (a closed form expression) corresponding to a uniform disk, and a ``residual'' pair. This residual component is converted into an infinite series of integrals over the radial extent of the disk. For a certain class of surface density distributions (like power laws of the radius), this series is fully analytical. The extraction of the homogeneous pair is equivalent to a convergence acceleration technique, in a matematical sense. In the case of power law distributions, the convergence rate of the residual series is shown to be cubic inside the source. As a consequence, very accurate potential values are obtained by low order truncation of the series. At zero order, relative errors on potential values do not exceed a few percent typically, and scale with the order N of truncation as 1/N**3. This method is superior to the classical multipole expansion whose very slow convergence is often critical for most practical applications.Comment: Accepted for publication in Astronomy & Astrophysics 7 pages, 8 figures, F90-code available at http://www.obs.u-bordeaux1.fr/radio/JMHure/intro2applawd.htm

Rotating molecular outflows: the young T Tauri star in CB26

The disk-outflow connection is thought to play a key role in extracting excess angular momentum from a forming proto-star. Though jet rotation has been observed in a few objects, no rotation of molecular outflows has been unambiguously reported so far. We report new millimeter-interferometric observations of the edge-on T Tauri star - disk system in the isolated Bok globule CB26. The aim of these observations was to study the disk-outflow relation in this 1Myr old low-mass young stellar object. The IRAM PdBI array was used to observe 12CO(2-1) at 1.3mm in two configurations, resulting in spectral line maps with 1.5 arcsec resolution. We use an empirical parameterized steady-state outflow model combined with 2-D line radiative transfer calculations and chi^2-minimization in parameter space to derive a best-fit model and constrain parameters of the outflow. The data reveal a previously undiscovered collimated bipolar molecular outflow of total length ~2000 AU, escaping perpendicular to the plane of the disk. We find peculiar kinematic signatures that suggest the outflow is rotating with the same orientation as the disk. However, we could not ultimately exclude jet precession or two misaligned flows as possible origin of the observed peculiar velocity field. There is indirect indication that the embedded driving source is a binary system, which, together with the youth of the source, could provide the clue to the observed kinematic features of the outflow. CB26 is so far the most promising source to study the rotation of a molecular outflow. Assuming that the outflow is rotating, we compute and compare masses, mass flux, angular momenta, and angular momentum flux of disk and outflow and derive disk dispersal timescales of 0.5...1 Myr, comparable to the age of the system.Comment: 14 pages, 6 figures, to appear in Astronomy & Astrophysic