We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present ALMA and ACA [CI]$^{3}P_{1}-^{3}P_{0}$ ([CI](1-0)) observations of NGC6240, which we combine with ALMA CO(2-1) and IRAM PdBI CO(1-0) data to study the physical properties of the massive molecular (H$_2$) outflow. We discover that the receding and approaching sides of the H$_2$ outflow, aligned east-west, exceed 10 kpc in their total extent. High resolution ($0.24"$) [CI](1-0) line images surprisingly reveal that the outflow emission peaks between the two AGNs, rather than on either of the two, and that it dominates the velocity field in this nuclear region. We combine the [CI](1-0) and CO(1-0) data to constrain the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$) in the outflow, which is on average $2.1\pm1.2~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$. We estimate that $60\pm20$ % of the total H$_2$ gas reservoir of NGC6240 is entrained in the outflow, for a resulting mass-loss rate of $\dot{M}_{\rm out}=2500\pm1200~M_{\odot}~yr^{-1}$ $\equiv 50\pm30$ SFR. This energetics rules out a solely star formation-driven wind, but the puzzling morphology challenges a classic radiative-mode AGN feedback scenario. For the quiescent gas we compute $\langle\alpha_{\rm CO}\rangle = 3.2\pm1.8~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$, which is at least twice the value commonly employed for (U)LIRGs. We observe a tentative trend of increasing $r_{21}\equiv L^{\prime}_{\rm CO(2-1)}/L^{\prime}_{\rm CO(1-0)}$ ratios with velocity dispersion and measure $r_{21}>1$ in the outflow, whereas $r_{21}\simeq1$ in the quiescent gas. We propose that molecular outflows are the location of the warmer, strongly unbound phase that partially reduces the opacity of the CO lines in (U)LIRGs, hence driving down their global $\alpha_{\rm CO}$ and increasing their $r_{21}$ values.

All measurements of cosmic star formation must assume an initial distribution of stellar masses -- the stellar initial mass function -- in order to extrapolate from the star-formation rate measured for typically rare, massive stars (> 8 Msun) to the total star-formation rate across the full stellar mass spectrum. The shape of the stellar initial mass function in various galaxy populations underpins our understanding of the formation and evolution of galaxies across cosmic time. Classical determinations of the stellar initial mass function in local galaxies are traditionally made at ultraviolet, optical and near-infrared wavelengths, which cannot be probed in dust-obscured galaxies, especially in distant starbursts, whose apparent star-formation rates are hundreds to thousands of times higher than in our Milky Way, selected at submillimetre (rest-frame far-infrared) wavelengths. The 13C/18O abundance ratio in the cold molecular gas -- which can be probed via the rotational transitions of the 13CO and C18O isotopologues -- is a very sensitive index of the stellar initial mass function, with its determination immune to the pernicious effects of dust. Here we report observations of 13CO and C18O emission for a sample of four dust-enshrouded starbursts at redshifts of approximately two to three, and find unambiguous evidence for a top-heavy stellar initial mass function in all of them. A low 13CO/C18O ratio for all our targets -- alongside a well-tested, detailed chemical evolution model benchmarked on the Milky Way -- implies that there are considerably more massive stars in starburst events than in ordinary star-forming spiral galaxies. This can bring these extraordinary starbursts closer to the `main sequence' of star-forming galaxies, though such main-sequence galaxies may not be immune to changes in initial stellar mass function, depending upon their star-formation densities.

In this work we extend the work on the recently discovered role of Cosmic Rays (CRs) in regulating the average CO/$\rm H_2$ abundance ratio in molecular clouds (and thus their CO line visibility) in starburst galaxies, and find that it can lead to a CO-poor/CI-rich $\rm H_2 $ gas phase even in environments with Galactic or in only modestly enhanced CR backgrounds expected in ordinary star-forming galaxies. Furthermore, the same CR-driven astro-chemistry raises the possibility of a widespread phase transition of molecular gas towards a CO-poor/CI-rich phase in: a) molecular gas outflows found in star-forming galaxies, b) active galactic nuclei (AGNs), and c) near synchrotron-emitting radio jets and the radio-loud cores of powerful radio galaxies. For main sequence galaxies we find that CRs can render some of their molecular gas mass CO-invisible, compounding the effects of low metallicities. Imaging the two fine structure lines of atomic carbon with resolution high enough to search beyond the CI/CO-bright line regions associated with central starbursts can reveal such a CO-poor/CI-rich molecular gas phase, provided that relative brightness sensitivity levels of $T_b$(CI $1-0$)/$T_b$(CO $J=1-0$)$\sim $0.15 are reached. The capability to search for such gas in the Galaxy is now at hand with the new high-frequency survey telescope HEAT deployed in Antarctica and future ones to be deployed in Dome A. ALMA can search for such gas in star-forming spiral disks, galactic molecular gas outflows and the CR-intense galactic and circumgalactic gas-rich environments of radio-loud objects.

Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atomic Carbon in high-z galaxies, at previously inaccessible frequencies from the ground. We have targeted three gravitationally lensed galaxies at redshift of about 3 and conducted a comparative study of the observed high-J CO/CI ~ratios with well-studied nearby galaxies. Atomic Carbon (CI(2-1)) was detected in one of the three targets and marginally in a second, while in all three targets the $J=7\to6$ CO line is detected. The CO(7-6)/CI(2-1), CO(7-6)/CO(1-0) line ratios and the CO(7-6)/(far-IR continuum) luminosity ratio are compared to those of nearby objects. A large excitation status in the ISM of these high-z objects is seen, unless differential lensing unevenly boosts the CO line fluxes from the warm and dense gas more than the CO(1-0), CI(2-1), tracing a more widely distributed cold gas phase. We provide estimates of total molecular gas masses derived from the atomic Carbon and the Carbon monoxide CO(1-0), which within the uncertainties turn out to be equal.

We use state-of-the-art chemical models to track the cosmic evolution of the CNO isotopes in the interstellar medium (ISM) of galaxies, yielding powerful constraints on their stellar initial mass function (IMF). We re-assess the relative roles of massive stars, asymptotic giant branch (AGB) stars and novae in the production of rare isotopes such as 13C, 15N, 17O and 18O, along with 12C, 14N and 16O. The CNO isotope yields of super-AGB stars, novae and fast-rotating massive stars are included. Having reproduced the available isotope enrichment data in the solar neighbourhood, and across the Galaxy, and having assessed the sensitivity of our models to the remaining uncertainties, e.g. nova yields and star-formation history, we show that we can meaningfully constrain the stellar IMF in galaxies using C, O and N isotope abundance ratios. In starburst galaxies, where data for multiple isotopologue lines are available, we find compelling new evidence for a top-heavy stellar IMF, with profound implications for their star-formation rates and efficiencies, perhaps also their stellar masses. Neither chemical fractionation nor selective photodissociation can significantly perturb globally-averaged isotopologue abundance ratios away from the corresponding isotope ones, as both these processes will typically affect only small mass fractions of molecular clouds in galaxies. Thus the Atacama Large Millimetre Array now stands ready to probe the stellar IMF, and even the ages of specific starburst events in star-forming galaxies across cosmic time unaffected by the dust obscuration effects that plague optical/near-infrared studies.

We explore the effects of the expected higher cosmic ray (CR) ionization rates $\zeta_{\rm CR}$ on the abundances of carbon monoxide (CO), atomic carbon (C), and ionized carbon (C$^+$) in the H$_2$ clouds of star-forming galaxies. The study of Bisbas et al. (2015) is expanded by: a) using realistic inhomogeneous Giant Molecular Cloud (GMC) structures, b) a detailed chemical analysis behind the CR-induced destruction of CO, and c) exploring the thermal state of CR-irradiated molecular gas. CRs permeating the interstellar medium with $\zeta_{\rm CR}$$\gtrsim 10\times$(Galactic) are found to significantly reduce the [CO]/[H$_2$] abundance ratios throughout the mass of a GMC. CO rotational line imaging will then show much clumpier structures than the actual ones. For $\zeta_{\rm CR}$$\gtrsim 100\times$(Galactic) this bias becomes severe, limiting the utility of CO lines for recovering structural and dynamical characteristics of H$_2$-rich galaxies throughout the Universe, including many of the so-called Main Sequence (MS) galaxies where the bulk of cosmic star formation occurs. Both C$^+$ and C abundances increase with rising $\zeta_{\rm CR}$, with C remaining the most abundant of the two throughout H$_2$ clouds, when $\zeta_{\rm CR}\sim (1-100)\times$(Galactic). C$^+$ starts to dominate for $\zeta_{\rm CR}$$\gtrsim 10^3\times$(Galactic). The thermal state of the gas in the inner and denser regions of GMCs is invariant with $T_{\rm gas}\sim 10\,{\rm K}$ for $\zeta_{\rm CR}\sim (1-10)\times$(Galactic). For $\zeta_{\rm CR}$$\sim 10^3\times$(Galactic) this is no longer the case and $T_{\rm gas}\sim 30-50\,{\rm K}$ are reached. Finally we identify OH as the key species whose $T_{\rm gas}-$sensitive abundance could mitigate the destruction of CO at high temperatures.

HE0450-2958, an interacting quasar-starburst galaxy pair at $z$ = 0.285, is one of the best known examples of strong star formation activity in the presence of a quasar-driven jet. We present new multi-band JVLA-imaging covering 1 to 6 GHz and reaching an angular resolution of up to $0{}_{{}^.}^{"}6$ (a 6-fold improvement over existing radio data). We confirm the previous detection of a spatially extended radio component around the quasar indicating that there is on-going star formation activity in the quasar host galaxy. For the first time, we directly detect a jet-like bipolar outflow from the quasar aligned with its companion star-forming galaxy (SFG) and several blobs of ionized gas in its vicinity identified in previous studies. Within the companion SFG we find evidence for a flattening of the synchrotron spectral index towards the point of intersection with the jet axis, further suggesting that the outflow may actually be impacting its interstellar medium (ISM). We discuss two possible mechanisms that could have triggered the starburst in the companion SFG: a wet-dry merger with the quasar and jet-induced star formation. While triggering through interaction-driven gas dynamics cannot be excluded with current data, our new observations make HE0450-2958 a strong candidate for jet-induced star formation, and one of the rare links between local systems (like Minkowski's Object or Centaurus A) and the high-z regime where radio-optical alignments suggest that this phenomenon could be more common.

Images of dust continuum and carbon monoxide (CO) line emission are powerful tools for deducing structural characteristics of galaxies, such as disk sizes, H$_2$ gas velocity fields and enclosed H$_2$ and dynamical masses. We report on a fundamental constraint set by the cosmic microwave background (CMB) on the observed structural and dynamical characteristics of galaxies, as deduced from dust continuum and CO-line imaging at high redshifts. As the CMB temperature rises in the distant Universe, the ensuing thermal equilibrium between the CMB and the cold dust and H$_2$ gas progressively erases all spatial and spectral contrasts between their brightness distributions and the CMB. For high-redshift galaxies, this strongly biases the recoverable H$_2$ gas and dust mass distributions, scale lengths, gas velocity fields and dynamical mass estimates. This limitation is unique to mm/submm wavelengths and unlike its known effect on the global dust continuum and molecular line emission of galaxies, it cannot be addressed simply. We nevertheless identify a unique signature of CMB-affected continuum brightness distributions, namely an increasing rather than diminishing contrast between such brightness distributions and the CMB when the cold dust in distant galaxies is imaged at frequencies beyond the Raleigh-Jeans limit. For the molecular gas tracers, the same effect makes the atomic carbon (CI) lines maintain a larger contrast than the CO lines against the CMB.

We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests.

We use a high-resolution hydrodynamic simulation that tracks the non-equilibrium abundance of molecular hydrogen within a massive high-redshift galaxy to produce mock Atacama Large Millimeter Array (ALMA) maps of the fine-structure lines of atomic carbon, CI 1-0 and CI 2-1. Inspired by recent observational and theoretical work, we assume that CI is thoroughly mixed within giant molecular clouds and demonstrate that its emission is an excellent proxy for H2. Nearly all of the H2 associated with the galaxy can be detected at redshifts z<4 using a compact interferometric configuration with a large synthesized beam (that does not resolve the target galaxy) in less than 4 h of integration time. Low-resolution imaging of the \CI lines (in which the target galaxy is resolved into three to four beams) will detect ~80 per cent of the H2 in less than 12 h of aperture synthesis. In this case, the resulting data cube also provides the crucial information necessary for determining the dynamical state of the galaxy. We conclude that ALMA observations of the CI 1-0 and 2-1 emission are well-suited for extending the interval of cosmic look-back time over which the H2 distributions, the dynamical masses, and the Tully-Fisher relation of galaxies can be robustly probed.

We report mapping observations of the $^{12}$CO $J=3-2$, $4-3$, $6-5$, and $7-6$ transitions and the CI 492 GHz transition toward the central 40$''\times$40$''$ region of the Circinus galaxy, using the Atacama Pathfinder EXperiment (APEX) telescope. We also detected $^{13}$CO $J=3-2$ at the central position of Circinus. These observations are to date the highest CO transitions reported in Circinus. With large velocity gradient (LVG) modeling and likelihood analysis we try to obtain density, temperature, and column density of the molecular gas in three regions: the nuclear region ($D< 18''\sim$ 360 pc), the entire central 45$''$ ($D<45''\sim$ 900pc) region, and the star-forming (S-F) ring (18$''<D<45''$). In the nuclear region, we can fit the CO excitation with a single excitation component, yielding an average condition of $n_{\rm H_2} \sim 10^{3.2}$cm$^{-3}$, $T_{\rm kin}$ $\sim$ 200 K, and d$v$/d$r \sim$ 3 km s$^{-1}$pc$^{-1}$. In the entire 45$''$ region, two excitation components are needed with $n_{\rm H_2}$ $\sim$ 10$^{4.2}$ and 10$^{3.0}$ cm$^{-3}$, $T_{\rm kin}\sim$ 60 K and 30 K, and $M_{\rm H_2}\sim$ $2.3\times 10^7$ M$_\odot$ and $6.6 \times 10^7$ M$_\odot$, respectively. The gas excitation in the S-F ring can also be fitted with two LVG components, after subtracting the CO fluxes in the 18$''$ region. The S-F ring region contributes 80\% of the molecular mass in the 45$''$ region. For the 45$''$ region, we find a conversion factor of $N({\rm H_2})/I_{\rm CO\ 1-0}$ = $0.37\times 10^{20}$ cm$^{-2} ({\rm K\ km\ s}^{-1})^{-1}$, about 1/5 of the Galactic disk value. The luminosity ratios of CI and $^{12}$CO $J=3-2$ ($R_{\rm CI/CO\ J=3-2}$) in Circinus basically follow a linear trend. The average $R_{\rm CI/CO\ J=3-2}$ in Circinus is $\sim$ 0.2, lying at an intermediate value between non-AGN nuclear region and high-redshift galaxies.

We used the SPIRE/FTS instrument aboard the Herschel Space Observatory (HSO) to obtain the Spectral Line Energy Distributions (SLEDs) of CO from J=4-3 to J=13-12 of Arp 193 and NGC 6240, two classical merger/starbursts selected from our molecular line survey of local Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11} L_{sol}). The high-J CO SLEDs are then combined with ground-based low-J CO, {13}CO, HCN, HCO+, CS line data and used to probe the thermal and dynamical states of their large molecular gas reservoirs. We find the two CO SLEDs strongly diverging from J=4-3 onwards, with NGC6240 having a much higher CO line excitation than Arp193, despite their similar low-J CO SLEDs and L_{FIR}/L_{CO,1-0}, L_{HCN}/L_{CO} (J=1-0) ratios (proxies of star formation efficiency and dense gas mass fraction). In Arp193, one of the three most extreme starbursts in the local Universe, the molecular SLEDs indicate a small amount ~(5-15)% of dense gas (n>=10^{4}cm^{-3}) unlike NGC6240 where most of the molecular gas (~(60-70)%) is dense n~(10^4-10^5)cm^{-3}. Strong star-formation feedback can drive this disparity in their dense gas mass fractions, and also induce extreme thermal and dynamical states for the molecular gas.In NGC6240, and to a lesser degree in Arp193, we find large molecular gas masses whose thermal states cannot be maintained by FUV photons from Photon Dominated Regions (PDRs). We argue that this may happen often in metal-rich merger/starbursts, strongly altering the initial conditions of star formation. ALMA can now directly probe these conditions across cosmic epoch, and even probe their deeply dust-enshrouded outcome, the stellar IMF averaged over galactic evolution.

Using the IRAM 30-m telescope, CN and CO isotopologues have been measured toward the central regions of the nearby starburst galaxy NGC253 and the prototypical ultraluminous infrared galaxy Mrk231. In NGC253, the 12C/13C ratio is 40+-10. Assuming that the ratio also holds for the CO emitting gas, this yields 16O/18O = 145+-36 and 16O/17O = 1290+-365 and a 32S/34S ratio close to that measured for the local interstellar medium (20-25). No indication for vibrationally excited CN is found. Peak line intensity ratios between NGC253 and Mrk231 are ~100 for 12C16O and 12C18O J=1-0, while the ratio for 13C16O J=1-0 is ~250. This and similar 13CO and C18O line intensities in the J=1-0 and 2-1 transitions of Mrk231 suggest 12C/13C ~ 100 and 16O/18O ~ 100, in agreement with values obtained for the less evolved ultraluminous merger Arp220. Also accounting for other extragalactic data, 12C/13C ratios appear to vary over a full order of magnitude, from >100 in ultraluminous high redshift galaxies to ~100 in more local such galaxies to ~40 in weaker starbursts not undergoing a large scale merger to 25 in the Central Molecular Zone of the Milky Way. With 12C being predominantly synthesized in massive stars, while 13C is mostly ejected by longer lived lower mass stars at later times, this is qualitatively consistent with our results of decreasing carbon isotope ratios with time and rising metallicity. It is emphasized, however, that both infall of poorly processed material, initiating a nuclear starburst, as well as the ejecta from newly formed massive stars (in particular in case of a top-heavy stellar initial mass function) can raise the carbon isotope ratio for a limited amount of time.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

Results from a large, multi-J CO, {13}CO, and HCN line survey of Luminous Infrared Galaxies (L_{IR}>=10^{10} L_{\odot}) in the local Universe (z<=0.1), complemented by CO J=4--3 up to J=13--12 observations from the Herschel Space Observatory (HSO), paints a new picture for the average conditions of the molecular gas of the most luminous of these galaxies with turbulence and/or large cosmic ray (CR) energy densities U_{CR} rather than far-UV/optical photons from star-forming sites as the dominant heating sources. Especially in ULIRGs (L_{IR}>10^{12} L_{\odot}) the Photon Dominated Regions (PDRs) can encompass at most \sim few% of their molecular gas mass while the large U_{CR} and the strong turbulence in these merger/starbursts, can volumetrically heat much of their molecular gas to T_{kin}\sim(100-200)K, unhindered by the high dust extinctions. Moreover the strong supersonic turbulence in ULIRGs relocates much of their molecular gas at much higher average densities than in isolated spirals. This renders low-J CO lines incapable of constraining the properties of the bulk of the molecular gas in ULIRGs, with substantial and systematic underestimates of its mass possible when only such lines are used. A comparative study of multi-J HCN lines and CO SLEDs from J=1--0 up to J=13--12 of NGC 6240 and Arp 193 offers a clear example of two merger/starbursts whose similar low-J CO SLEDs, and L_{IR}/L_{CO,1-0}, L_{HCN, 1-0}/L_{CO,1-0} ratios, yield no indications about their strongly diverging CO SLEDs beyond J=4--3, and ultimately the different physical conditions in their molecular ISM. The much larger sensitivity of ALMA and its excellent site in the Atacama desert now allows the observations necessary to ....

We make use of our 'minimal' cold interstellar medium (ISM) emission line model that predicts the molecular and atomic line emission per unit dense, star-forming gas mass (Geach & Papadopoulos 2012; Paper I) to examine the utility of key line ratios in surveys of the so-called star formation 'mode' as traced by xi_SF = M_dense(H_2)/M_total(H_2). We argue that xi_SF and its proxies provide very sensitive, extinction-free discriminators of rapid starburst/merger-driven versus secular quiescent/disk-like stellar mass assembly, with the most promising diagnostic to be applied in the near-future being CO(4-3)/[CI](1-0). These lines are accessible across nearly the full range 0<z<2 (thus covering the bulk of galaxy evolution) with the Atacama Large Millimeter Array. In addition to their diagnostic power, another advantage of this combination is the similar observed frequencies (Delta nu_0 ~ 30 GHz) of the lines, resulting in nearly spatially-matched beams for a fixed aperture, thus mitigating the effects of resolution/morphology bias in the interpretation of galaxy-averaged line ratios. Finally we discuss the capability of deep blind redshift surveys with the high frequency component of the Square Kilometer Arrray (SKA) in discovering H_2-rich galaxies with very low xi_SF values. These could be the progenitors of starburst galaxies seen prior to the onset of star formation; such galaxies could be a class of extreme outliers from local (gas surface density)-(star formation rate) scaling laws, which would exclude them from current star formatation or stellar mass selected samples. Our conservative model suggests that SKA could detect such systems residing at z~3 at a rate of 20-200 per hour.

We present a study using the Karl G. Jansky Very Large Array of 12CO J=1-0 emission in three strongly lensed submillimetre-selected galaxies (SMMJ16359, SMMJ14009 and SMMJ02399) at z=2.5-2.9. These galaxies span L(IR) = 10^11 - 10^13 Lsun, offering an opportunity to compare the interstellar medium of LIRGs and ULIRGs at high redshift. We estimate molecular gas masses in the range (2-40) x 10^9 Msun using a method that assumes canonical underlying brightness temperature ratios for star-forming and non-star-forming gas phases and a maximal star-formation efficiency. A more simplistic method using X(CO) = 0.8 yields gas masses twice as high. The observed CO(3-2)/CO(1-0) brightness temperature ratio for SMMJ14009, r(3-2)/(1-0) = (0.95 \pm 0.12), is indicative of warm star-forming gas, possibly influenced by the central AGN. We search for 12CO(1-0) emission in the Lyman-break galaxy, A2218 #384, located at z=2.517 in the same field as SMMJ16359, and assign a 3-sigma gas mass limit of <6 x 10^8 Msun. We use rest-frame 115-GHz free-free flux densities in SMMJ14009 and SMMJ02399 - measurements tied directly to the photionisation rate of massive stars and made possible by the VLA's bandwidth - to estimate star-formation rates of 400-600 Msun/yr and to estimate the fraction of L(IR) due to the AGN.

Cosmic rays (CRs) control the thermal, ionization and chemical state of the dense H_2 gas regions that otherwise remain shielded from far-UV and optical stellar radiation propagating through the dusty ISM of galaxies. It is in such CR-dominated regions (CRDRs) rather than Photon-dominated regions (PDRs) of H_2 clouds where the star formation initial conditions are set, making CRs the ultimate star-formation feedback factor in galaxies, able to operate even in their most deeply dust-enshrouded environments. CR-controlled star formation initial conditions naturally set the stage for a near-invariant stellar Initial Mass Function (IMF) in galaxies as long as their average CR energy density U_{CR} permeating their molecular ISM remains within a factor of ~10 of its Galactic value. Nevertheless, in the extreme environments of the compact starbursts found in merging galaxies, where U_{CR}\sim(few)x10^{3}U_{CR,Gal}, CRs dramatically alter the initial conditions of star formation. In the resulting extreme CRDRs H_2 cloud fragmentation will produce far fewer low mass (<8 M_{sol}) stars, yielding a top-heavy stellar IMF. This will be a generic feature of CR-controlled star-formation initial conditions, lending a physical base for a bimodal IMF during galaxy formation, with a top-heavy one for compact merger-induced starbursts, and an ordinary IMF preserved for star formation in isolated gas-rich disks. In this scheme the integrated galactic IMFs (IGIMF) are expected to be strong functions of the star formation history of galaxies.

We predict the space density of molecular gas reservoirs in the Universe, and place a lower limit on the number counts of carbon monoxide (CO), hydrogen cyanide (HCN) molecular and [CII] atomic emission lines in blind redshift surveys in the submillimeter-centimeter spectral regime. Our model uses: (a) recently available HCN Spectral Line Energy Distributions (SLEDs) of local Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun), (b) a value for epsilon=SFR/M_dense(H_2) provided by new developments in the study of star formation feedback on the interstellar medium and (c) a model for the evolution of the infrared luminosity density. Minimal 'emergent' CO SLEDs from the dense gas reservoirs expected in all star-forming systems in the Universe are then computed from the HCN SLEDs since warm, HCN-bright gas will necessarily be CO-bright, with the dense star-forming gas phase setting an obvious minimum to the total molecular gas mass of any star-forming galaxy. We include [CII] as the most important of the far-infrared cooling lines. Optimal blind surveys with the Atacama Large Millimeter Array (ALMA) could potentially detect very distant (z~10-12) [CII] emitters in the >ULIRG galaxy class at a rate of ~0.1-1 per hour (although this prediction is strongly dependent on the star formation and enrichment history at this early epoch), whereas the (high-frequency) Square Kilometer Array (SKA) will be capable of blindly detecting z>3 low-J CO emitters at a rate of ~40-70 per hour. The [CII] line holds special promise for the detection of metal-poor systems with extensive reservoirs of CO-dark molecular gas where detection rates with ALMA can reach up to 2-7 per hour in Bands 4-6.

We report results from a large molecular line survey of Luminous Infrared Galaxies (L_{IR} >= 10^{11} L_sol) in the local Universe (z<=0.1), conducted during the last decade with the James Clerk Maxwell Telescope (JCMT) and the IRAM 30-m telescope. This work presents the CO and {13}CO line data for 36 galaxies, further augmented by multi-J total CO luminosities available for other IR-bright galaxies from the literature. This yields a sample of N=70 galaxies with the star-formation (SF) powered fraction of their IR luminosities spanning L_{IR} (10^{10}-2x10^{12}) L_sol and a wide range of morphologies. Simple comparisons of their available CO Spectral Line Energy Distributions (SLEDs) with local ones, as well as radiative transfer models discern a surprisingly wide range of average ISM conditions, with most of the surprises found in the high-excitation regime. These take the form of global CO SLEDs dominated by a very warm (T_{kin}>=100 K) and dense (n>=10^4 cm^{-3}) gas phase, involving galaxy-sized (~(few)x10^9 M_sol) gas mass reservoirs under conditions that would otherwise amount only ~1% of mass per typical SF molecular cloud in the Galaxy. Some of the highest excitation CO SLEDs are found in the so-called Ultra Luminous Infrared Galaxies and seem irreducible to ensembles of ordinary SF-powered regions. Highly supersonic turbulence and high cosmic ray (CR) energy densities rather than far-UV/optical photons or SNR-induced shocks from individual SF sites can globally warm the large amounts of dense gas found in these merger-driven starbursts and easily power their extraordinary CO line excitation.....

In this work we conclude the analysis of our CO line survey of Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11}L_{sol}) in the local Universe (Paper\,I), by focusing on the influence of their average ISM properties on the total molecular gas mass estimates via the so-called X_{co}=M(H_2)/L_{co,1-0} factor. One-phase radiative transfer models of the global CO Spectral Line Energy Distributions (SLEDs) yield an X_{co} distribution with: <X_{co}>\sim(0.6+/-0.2) M_{sol}(K km s^{-1} pc^2)^{-1} over a significant range of average gas densities, temperatures and dynamical states. The latter emerges as the most important parameter in determining X_{co}, with unbound states yielding low values and self-gravitating states the highest ones. Nevertheless in many (U)LIRGs where available higher-J CO lines (J=3--2, 4--3, and/or J=6--5) or HCN line data from the literature allow a separate assessment of the gas mass at high densities (>=10^{4} cm^{-3}) rather than a simple one-phase analysis we find that {\it near-Galactic X_{co} (3-6)\, M_sol\,(K\,km^{-1}\,pc^2)^{-1} values become possible.} We further show that in the highly turbulent molecular gas in ULIRGs a high-density component will be common and can be massive enough for its high X_{co} to dominate the average value for the entire galaxy. ......... ...this may have thus resulted to systematic underestimates of molecular gas mass in ULIRGs.

We have used the IRAM Plateau de Bure Interferometer and the Expanded Very Large Array to obtain a high resolution map of the CO(6-5) and CO(1-0) emission in the lensed, star-forming galaxy SMMJ2135-0102 at z=2.32. The kinematics of the gas are well described by a model of a rotationally-supported disk with an inclination-corrected rotation speed, v_rot = 320+/-25km/s, a ratio of rotational- to dispersion- support of v/sigma=3.5+/-0.2 and a dynamical mass of 6.0+/-0.5x10^10Mo within a radius of 2.5kpc. The disk has a Toomre parameter, Q=0.50+/-0.15, suggesting the gas will rapidly fragment into massive clumps on scales of L_J ~ 400pc. We identify star-forming regions on these scales and show that they are 10x denser than those in quiescent environments in local galaxies, and significantly offset from the local molecular cloud scaling relations (Larson's relations). The large offset compared to local molecular cloud linewidth-size scaling relations imply that supersonic turbulence should remain dominant on scales ~100x smaller than in the kinematically quiescent ISM of the Milky Way, while the molecular gas in SMMJ2135 is expected to be ~50x denser than that in the Milky Way on all scales. This is most likely due to the high external hydrostatic pressure we measure for the interstellar medium (ISM), P_tot/kB ~ (2+/-1)x10^7K/cm3. In such highly turbulent ISM, the subsonic regions of gravitational collapse (and star-formation) will be characterised by much higher critical densities, n_crit>=10^8/cm3, a factor ~1000x more than the quiescent ISM of the Milky Way.

We report the detection of the CO J=1-0 emission line in [H89]1821+643, one of the most optically luminous QSOs in the local Universe and a template ULIRG-to-QSO transition object, located in a rich cool-core cluster at z=0.297. The CO emission is likely to be extended, highly asymmetric with respect to the center of the host elliptical where the QSO resides, and corresponds to a molecular gas mass of ~ 8.0 x 10^{9} M_sun. The dynamical mass enclosed by the CO line emission region could amount up to ~ 1.7x10^{12} M_sun (80% of the total mass of the elliptical host). The bulk of the CO emission is located at ~9 kpc south-east from the nuclei position, and close to a faint optical structure, suggesting that the CO emission could either represent a gas-rich companion galaxy merging with the elliptical host or a tail-like structure reminiscent of a previous interaction. We argue that the first scenario is more likely given the large masses implied for the CO source, which would imply a highly asymmetrical elliptical host. The close alignment between the CO emission major axis and the radio-plume suggests a possible role of the latter excitation of the ambient gas reservoir. The stacking technique was used to search for CO emission and 3mm continuum emission from galaxies in the surrounding cluster, however no detection was found either toward individual galaxies or the stacked ensemble of galaxies, with a 3 sigma limit of < 1.1x10^9 M_sun for the molecular gas.

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CRDRs) characterize the ISM of galaxies during events of high-density star formation, fundamentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray energy densities in such systems (U_{CR} (few)x(10^3-10^4) U_{CR,Gal}) will indeed raise the minimum temperature of this phase (where the initial conditions of star formation are set) from ~10K (as in the Milky Way) to (50-100)K. Moreover in such extreme CRDRs the gas temperature remains fully decoupled from that of the dust, with T_{kin} >> T_{dust}, even at high densities (n(H_2)~10^5--10^6 cm^{-3}), quite unlike CRDRs in the Milky Way where T_k T_{dust} when n(H_2) >= 10^5 cm^{-3}. These dramatically different star formation initial conditions will: a) boost the Jeans mass of UV-shielded gas regions by factors of ~10--100 with respect to those in quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point of the effective equation of state (EOS) of molecular gas. Both these effects occur across the entire density range of typical molecular clouds, and may represent {\it a new paradigm for all high-density star formation in the Universe}, with cosmic rays as the key driving mechanism, operating efficiently even in the high dust extinction environments of extreme starbursts...

We report the results of a pilot study with the EVLA of 12CO J=1-0 emission from four SMGs at z=2.2-2.5, each with an existing detection of CO J=3-2. Using the EVLA's most compact configuration we detect strong, broad J=1-0 line emission from all of our targets. The median line width ratio, sigma(1-0)/sigma(3-2) = 1.15 +/- 0.06, suggests that the J=1-0 is more spatially extended than the J=3-2 emission, a situation confirmed by our maps which reveal velocity structure in several cases and typical sizes of ~16 kpc FWHM. The median Tb ratio is r(3-2/1-0) = 0.55 +/- 0.05, noting that our value may be biased high because of the J=3-2-based sample selection. Naively, this suggests gas masses ~2x higher than estimates made using higher-J transitions of CO, with the discrepency due to the difference in assumed Tb ratio. We also estimate masses using the 12CO J=1-0 line and the observed global Tb ratios, assuming standard underlying Tb ratios as well as a limiting SFE, i.e. without calling upon X(CO). Using this new method, we find a median molecular gas mass of (2.5 +/- 0.8) x 10^10 Msun, with a plausible range stretching 3x higher. Even larger masses cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical mass for our sample is (2.3 +/- 1.4) x 10^11 Msun. We examine the Schmidt-Kennicutt relation for all the distant galaxy populations for which CO J=1-0 or J=2-1 data are available, finding small systematic differences. These have previously been interpreted as evidence for different modes of star formation, but we argue that these differences are to be expected, given the still considerable uncertainties. Finally, we discuss the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for breaking them, and the future prospects for imaging and studying cold, quiescent molecular gas at high redshifts [abridged].

The high-density star formation typical of the merger/starburst events that power the large IR luminosities of Ultra Luminous Infrared Galaxies (ULIRGs) (L_{IR}>10^{12}Lsol) throughout the Universe results to extraordinarily high cosmic ray (CR) energy densities of U_CR~(few)x(10^3--10^4)U_{CR,Gal} permeating their interstellar medium (ISM), a direct consequence of the large supernovae remnants (SNRs) number densities in such systems. Unlike far-UV photons emanating from their numerous star forming sites, these large CR energy densities in ULIRGs will volumetrically heat and raise the ionization fraction of dense (n>10^4 cm^{-3}) UV-shielded gas cores throughout their compact star-forming volumes. Such conditions can turn most of the large molecular gas masses found in such systems and their high redshift counterparts (M(H2)~10^9-10^10 M_{sol}) into giant CR-dominated Regions (CRDRs) rather than ensembles of Photon-dominated Regions (PDRs) which dominate in less IR-luminous systems where star formation and molecular gas distributions are much more extended. The molecular gas in CRDRs will have a {\it minimum} temperature of T_{kin}~(80--160)K, and very high ionization fractions of x(e)>10^{-6} throughout its UV-shielded dense cores, which in turn will {\it fundamentally alter the initial conditions for star formation in such systems.}. Observational tests of CRDRs can be provided by ......

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

We utilize detailed time-varying models of the coupled evolution of stars and the HI, H_2, and CO-bright H_2 gas phases in galaxy-sized numerical simulations to explore the evolution of gas-rich and/or metal-poor systems, expected to be numerous in the Early Universe. The inclusion of the CO-bright H_2 gas phase, and the realistic rendering of star formation as an H_2-regulated process (and the new feedback processes that this entails) allows the most realistic tracking of strongly evolving galaxies, and much better comparison with observations. We find that while galaxies eventually settle into states conforming to Schmidt-Kennicutt (S-K) relations, significant and systematic deviations of their star formation rates (SFRs) from the latter occur, especially pronounced and prolonged for ... ...This indicates potentially serious limitations of (S-K)-type relations as reliable sub-grid elements of star formation physics in simulations of structure formation in the Early Universe. We anticipate that galaxies with marked deviations from the S-K relations will be found at high redshifts as unbiased inventories of total gas mass become possible with ALMA and the EVLA.

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

Luminous infrared galaxies ($L_{\rm{IR}}>10^{11} L_{\odot}$) are often associated with interacting galactic systems and are thought to be powered by merger--induced starbursts and/or dust--enshrouded AGN. In such systems, the evolution of the dense, star forming molecular gas as a function of merger separation is of particular interest. Here, we present observations of the CO(3-2) emission from a sample of luminous infrared galaxy mergers that span a range of galaxy-galaxy separations. The excitation of the molecular gas is studied by examining the CO(3-2)/CO(1-0) line ratio, $r_{31}$, as a function of merger extent. We find these line ratios, $r_{31}$, to be consistent with kinetic temperatures of $T_k$=(30--50) K and gas densities of $n_{\rm{H}_2}=10^3 \rm{cm}^{-3}$. We also find weak correlations between $r_{31}$ and both merger progression and star formation efficiency ($L_{\rm{fIR}} / L_{\rm{CO(1-0)}}$). These correlations show a tendency for gas excitation to increase as the merger progresses and the star formation efficiency rises. To conclude, we calculate the contributions of the CO(3-2) line to the 850 $\mu$m fluxes measured with SCUBA, which are seen to be significant ($\sim$24%).

We report on new sensitive CO J=6-5 line observations of several luminous infrared Galaxies (LIRGs: L$_{\rm IR}$(8-1000$\mu $m$)\ga 10^{11}$ L$_{\odot}$), 36% (8/22) of them ULIRGs (L$_{\rm IR}$$>10^{12}$ L$_{\odot}$), and two powerful local AGN: the optically luminous QSO PG 1119+120, and the powerful radio galaxy 3C 293 using the James Clerk Maxwell Telescope (JCMT) on Mauna Kea in Hawaii. We combine these observations with existing low-J CO data and dust emission Spectral Energy Distributions (SEDs) in the far-infrared - submillimetre from the literature to constrain the properties of the star-forming ISM in these systems. We then build the first {\it local} CO Spectral Line Energy Distributions (SLEDs) for the {\it global} molecular gas reservoirs that reach up to high J-levels. These CO SLEDs are neither biased by strong lensing (which affects many of those constructed for high-redshift galaxies), nor suffer from undersampling of CO-bright regions (as most current high-J CO observations of nearby extended systems do). We find: ...

We report new single dish CO J=6--5 line observations for the archetypal Ultra Luminous Infrared Galaxy (ULIRG) Arp 220 with the James Clerk Maxwell Telescope atop Mauna Kea in Hawaii. The J=6--5 line is found to be faint, with brightness temperature ratios (6--5)/(1--0), (6--5)/(3--2) of R_{65/10} = 0.080+/-0.017 and R_{65/32} = 0.082+/-0.019, suggesting very low excitation conditions that cannot be reconciled with the warm and very dense molecular gas present in one of the most extreme starbursts in the local Universe. We find that an optically thick dust continuum, with tau(f>350 GHz)>=1 for the bulk of the warm dust and gas in Arp 220, submerges this line to an almost black body curve, reducing its flux, and {\it affecting its CO Spectral Line Energy Distribution (SLED) at high frequencies}. This also resolves the CII line deficiency in this object, first observed by {\it ISO}: the near absence of that line is a dust optical depth effect, not a dense Photodissociation Region (PDR) phenomenon. Finally we briefly comment on the possibility of such extreme ISM states in other ULIRGs in the distant Universe, and its consequences for the diagnostic utility of high frequency molecular and atomic ISM lines in such systems. In the case of Arp 220 we anticipate that the now spaceborne Herschel Space Observatory will find faint high-J CO lines at f>=690 GHz that would appear as sub-thermally excited with respect to the low-J ones as a result of the effects of dust absorption.

We present new multi-wavelength observations of the first submm-selected galaxy, SMM J02399-0136 at z=2.8. These observations include mapping of the CO J=1-0 emission using elements of the Expanded VLA, as well as high-resolution 1.4-GHz imaging and optical/IR data from the VLA, HST, Spitzer and Keck. Together these new data provide fundamental insights into the mass and distribution of stars, gas and dust within this archetypal SMG. The CO J=1-0 emission, with its minimal excitation and density requirements, traces the bulk of the metal-rich molecular gas, and reveals a mass of ~10^11 M(sun), extending over ~5" (~25 kpc in the source plane), although there is tentative evidence that it may be significantly larger. Our data suggest that three or more distinct structures are encompassed by this molecular gas reservoir, including the BAL quasar from which the redshift of the SMG was initially determined. In particular, the new rest-frame near-IR observations identify a massive, obscured, starburst coincident with a previously known Ly-alpha cloud. This starburst dominates the far-IR emission from the system and requires a re-assessment of previous claims that the gas reservoir resides in a massive, extended disk around the BAL QSO. Instead it appears that SMM J02399-0136 comprises a merger between a far-IR-luminous, but highly obscured starburst, the BAL QSO host and a faint third component. Our findings suggest that this SMG and its immediate environment mark a vast and complex galactic nursery and that detailed studies of other SMGs are likely to uncover a similarly rich diversity of properties.

We use time-varying models of the coupled evolution of the HI, H_2 gas phases and stars in galaxy-sized numerical simulations to: a) test for the emergence of the Kennicutt-Schmidt (K-S) and the H_2-pressure relation, b) explore a realistic H_2-regulated star formation recipe which brings forth a neglected and potentially significant SF-regulating factor, and c) go beyond typical galactic environments (for which these galactic empirical relations are deduced) to explore the early evolution of very gas-rich galaxies. In this work we model low mass galaxies ($M_{\rm baryon} \le 10^9 \msun$), while incorporating an independent treatment of CO formation and destruction, the most important tracer molecule of H2 in galaxies, along with that for the H2 gas itself. We find that both the K-S and the H_2-pressure empirical relations can robustly emerge in galaxies after a dynamic equilibrium sets in between the various ISM states, the stellar component and its feedback. (abridged)

We present a new observational study of the CO(1-0) line emission as an H2 gas mass tracer under extreme conditions in extragalactic environments. Our approach is to study the full neutral interstellar medium (H2, HI and dust) of two galaxies whose bulk interstellar medium (ISM) resides in environments that mark (and bracket) the excitation extremes of the ISM conditions found in infrared luminous galaxies, the starburst NGC3310 and the quiescent spiral NGC157. Our study maintains a robust statistical notion of the so-called X factor (i.e. a large ensemble of clouds is involved) while exploring its dependency on the very different average ISM conditions prevailing within these two systems. These are constrained by fully-sampled CO(3-2) and CO(1-0) observations, at a matched beam resolution of Half Power Beam Width 15", obtained with the JCMT the Nobeyama 45-m telescope, combined with sensitive 850 and 450 micron dust emission and HI interferometric images which allow a complete view of all the neutral ISM components. We found an X factor varying by a factor of 5 within the spiral galaxy NGC157 and about 2 times lower than the Galactic value in NGC3310. In addition, the dust emission spectrum in NGC3310 shows a pronounced submm "excess". We explain and fit this by an enhanced abundance of Very Small Grains (VSGs). The alternative approach to interpret this excess as very cold dust would yield a cold dust temperature of 5-11 K with a correspondingly very low gas-to-dust mass ratio of 5-43, very unlikely for this starburst galaxy with low metallicity.

We present new high resolution numerical simulations of the ISM in a central R ~32 parsecs region around a supermassive black hole at a galactic center. Three-dimensional hydrodynamic modeling of the ISM (Wada & Norman 2002) with the nuclear starburst now includes tracking of the formation of molecular hydrogen out of the neutral hydrogen phase. In a quasi equilibrium state, mass fraction of H2 is about 0.4 (total H2 mass is ~1.5 10^6 Msun) of the total gas mass for the uniform far UV (FUV) with G_0 = 10. The gas forms an inhomogeneous disk, whose scale-height becomes larger in the outer region. H2 forms a thin nuclear disk in the inner ~ 5 pc, which is surrounded by molecular clouds swelled up toward h < 10 pc. The velocity field of the disk is highly turbulent in the torus region, whose velocity dispersion is ~ 20 km/s on average. Average supernova rate (SNR) of ~ 5 10^-5/yr is large enough to energize these structures. Gas column densities toward the nucleus larger than 10^22 cm^-2 are observed if the viewing angle is smaller than \theta_v ~ 50 deg from the edge-on. However, the column densities are distributed over almost two orders of magnitude around the average for any given viewing angle due to the clumpy nature of the torus. For a stronger FUV (G_0 =100), the total H2 mass in an equillibrium is only slightly smaller (~ 0.35). Finally the morphology and kinematics of the circumnuclear molecular gas disks emerging from our models is similar to that revealed by recent near infrared observations using VLTI/Keck.

We report single-dish multi-transition measurements of the 12^CO, HCN, and HCO^+ molecular line emission as well as HNC J=1-0 and HNCO in the two ultraluminous infra-red galaxies Arp220 and NGC6240. Using this new molecular line inventory, in conjunction with existing data in the literature, we compiled the most extensive molecular line data sets to date for such galaxies. The many rotational transitions, with their different excitation requirements, allow the study of the molecular gas over a wide range of different densities and temperatures with significant redundancy, and thus allow good constraints on the properties of the dense gas in these two systems. The mass (~(1-2) x 10^10 Msun) of dense gas (>10^5-6 cm^-3) found accounts for the bulk of their molecular gas mass, and is consistent with most of their IR luminosities powered by intense star bursts while self-regulated by O,B star cluster radiative pressure onto the star-forming dense molecular gas. The highly excited HCN transitions trace a gas phase ~(10-100)x denser than that of the sub-thermally excited HCO^+ lines (for both galaxies). These two phases are consistent with an underlying density-size power law found for Galactic GMCs (but with a steeper exponent), with HCN lines tracing denser and more compact regions than HCO^+. Whether this is true in IR-luminous, star forming galaxies in general remains to be seen, and underlines the need for observations of molecular transitions with high critical densities for a sample of bright (U)LIRGs in the local Universe -- a task for which the HI-FI instrument on board Herschel is ideally suited to do.

We used the mm/sub-mm receivers on the James Clerk Maxwell Telescope (JCMT) to observe the CO J=3--2, 2--1 lines in five local, optically powerful AGN and the J=4--3 line in 3C 293 (a powerful radio galaxy). Luminous CO J=3--2 emission and high CO (3--2)/(1--0) intensity ratios are found in all objects, indicating highly excited molecular gas. In 3C 293 an exceptionally bright CO J=4--3 line is found which cannot be easily explained given its quiescent star-forming environment and low AGN X-ray luminosity. In this object shocks emanating from a well-known interaction of a powerful jet with a dense ISM may be responsible for the high excitation of its molecular gas on galaxy-wide scales. Star formation can readily account for the gas excitation in the rest of the objects, although high X-ray AGN luminosities can also contribute significantly in two cases. Measuring and eventually imaging CO line ratios in local luminous QSO hosts can be done by a partially completed ALMA during its early phases of commissioning, promising a sensitive probe of starburst versus AGN activity in obscured environments at high linear resolutions.

We report on interferometric imaging of the CO J=1--0 and J=3--2 line emission from the controversial QSO/galaxy pair HE 0450--2958. {\it The detected CO J=1--0 line emission is found associated with the disturbed companion galaxy not the luminous QSO,} and implies $\rm M_{gal}(H_2)\sim (1-2)\times 10^{10} M_{\odot}$, which is $\ga 30% $ of the dynamical mass in its CO-luminous region. Fueled by this large gas reservoir this galaxy is the site of an intense starburst with $\rm SFR\sim 370 M_{\odot} yr^{-1}$, placing it firmly on the upper gas-rich/star-forming end of Ultra Luminous Infrared Galaxies (ULIRGs, $\rm L_{IR}>10^{12} L_{\odot}$). This makes HE 0450--2958 the first case of extreme starburst and powerful QSO activity, intimately linked (triggered by a strong interaction) but not coincident. The lack of CO emission towards the QSO itself renews the controversy regarding its host galaxy by making a gas-rich spiral (the typical host of Narrow Line Seyfert~1 AGNs) less likely. Finally, given that HE 0450--2958 and similar IR-warm QSOs are considered typical ULIRG$\to $(optically bright QSO) transition candidates, our results raise the possibility that some may simply be {\it gas-rich/gas-poor (e.g. spiral/elliptical) galaxy interactions} which ``activate'' an optically bright unobscured QSO in the gas-poor galaxy, and a starburst in the gas-rich one. We argue that such interactions may have gone largely unnoticed even in the local Universe because the combination of tools necessary to disentagle the progenitors (high resolution and S/N optical {\it and} CO imaging) became available only recently.

We report on detections of the high-excitation CO J=6-5, J=4-3 lines in Mrk231, a prototypical Ultra Luminous Infrared Galaxy (ULIRG) and Seyfert 1 QSO. These observations are combined with CO J=3-2, HCN J=4-3 (this work), and CO J=2-1, J=1-0, 13CO J=2-1, HCN J=1-0 measurements taken from the literature to provide better constraints on the properties of the molecular gas in an extreme starburst/QSO in the local Universe. We find that the CO J=4-3 and J=6-5 transitions trace a different gas phase from that dominating the lower three CO transitions, with n(H_2) ~ (1-3)x10^4 cm-3 and Tk ~ (40-70) K. This phase is responsible for the luminous HCN emission, and contains most of the H2 gas mass of this galaxy. The total CO line cooling emanating from this dense phase is found similar to that of the [CII] line at 158 micron, suggesting a very different thermal balance to that seen in lower IR-luminosity galaxies, and one likely dominated by dense photon-dominated regions. Our dense "sampling" of the CO rotational ladder and the HCN lines enables us to produce well-constrained Spectral Line Energy Distributions (SLEDs) for the dense molecular gas in Mrk231 and compare them to those of high redshift starbursts, many of which have SLEDs that may be affected by strong lensing. Finally, we use our local molecular line excitation template to assess the capabilities of future cm and mm/sub-mm arrays in detecting CO and HCN transitions in similar systems throughout the local and distant universe.

We present Australia Telescope Compact Array radio continuum observations of the quasar/galaxy system HE0450-2958. An asymetric triple linear morphology is observed, with the central radio component coincident with the quasar core and a second radio component associated with a companion galaxy at a projected distance of 7kpc from the quasar. The system obeys the far-infrared to radio continuum correlation, implying the radio emission is energetically dominated by star formation activity. However, there is undoubtedly some contribution to the overall radio emission from a low-luminosity AGN core and a pair of radio lobes. Long baseline radio interferometric observations of the quasar core place a 3sigma upper limit of 0.6mJy at 1400MHz on the AGN contribution to the quasar's radio emission; less than 30% of the total. The remaining 70% of the radio emission from the quasar is associated with star formation activity and provides the first direct evidence for the quasar's host galaxy. A re-anlaysis of the VLT spectroscopic data shows extended emission line regions aligned with the radio axis and extended on scales of ~20kpc. This is interpreted as evidence for jet-cloud interactions, similar to those observed in radio galaxies and Seyferts. The emission lines in the companion galaxy are consistent with radiative shocks and its spatial association with the eastern radio lobe implies large-scale jet-induced star formation has played a role in this galaxy's evolution.

A large CO, HCN multi-transition survey of 30 Luminous Infrared Galaxies ($\rm L_{IR}>10^{11} L_{\odot}$) is nearing completion with the James Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii), and the IRAM 30-meter telescope at Pico Veleta (Spain). The CO J=1--0, 2--1, 3--2, 4--3,6--5, $ ^{13}$CO J=2--1, HCN J=1--0, 3--2, 4--3 observations, resulting from $\sim 250$ hours of JCMT, $\sim 100$ hours of 30-m observing time and data from the literature constitute {\it the largest extragalactic molecular line survey to date}, and can be used to address a wide range of issues and eventually allow the construction of reliable Spectral Line Energy Distributions (SLEDs) for the molecular gas in local starbursts. First results suggest that: a) HCN and HCO$^+$ J=1--0 line luminosities can be poor mass estimators of dense molecular gas ($\rm n\geq 10^4 cm^{-3}$) unless their excitation is accounted for, b) CO cooling of such gas in ULIRGs may be comparable to that of the CII line at $\rm 158 \mu m$, and c) low excitation of the {\it global} molecular gas reservoir remains possible in such systems. In such cases the expected low CO $\rm J+1\to J$ line luminosities for $\rm J+1\geq 4$ can lead to a strong bias against their detection from ULIRGs at high redshifts.

It has been recently argued that the HCN J=1--0 line emission may not be an unbiased tracer of dense molecular gas ($\rm n\ga 10^4 cm^{-3}$) in Luminous Infrared Galaxies (LIRGs: $\rm L_{FIR}> 10^{11} L_{\odot}$) and HCO$^+$ J=1--0 may constitute a better tracer instead (Graci\'a-Carpio et al. 2006), casting doubt into earlier claims supporting the former as a good tracer of such gas (Gao & Solomon 2004; Wu et al. 2006). In this paper new sensitive HCN J=4--3 observations of four such galaxies are presented, revealing a surprisingly wide excitation range for their dense gas phase that may render the J=1--0 transition from either species a poor proxy of its mass. Moreover the well-known sensitivity of the HCO$^+$ abundance on the ionization degree of the molecular gas (an important issue omitted from the ongoing discussion about the relative merits of HCN and HCO$^+$ as dense gas tracers) may severely reduce the HCO$^+$ abundance in the star-forming and highly turbulent molecular gas found in LIRGs, while HCN remains abundant. This may result to the decreasing HCO$^+$/HCN J=1--0 line ratio with increasing IR luminosity found in LIRGs, and casts doubts on the HCO$^+$ rather than the HCN as a good dense molecular gas tracer. Multi-transition observations of both molecules are needed to identify the best such tracer, its relation to ongoing star formation, and constrain what may be a considerable range of dense gas properties in such galaxies.

Ultra Luminous Infra-Red Galaxies (ULIRGs) -- gas-rich mergers harboring the most extreme star-forming conditions encountered in the local Universe -- are thought to re-enact the galaxy formation processes we are only barely able to glimpse in the distant Universe. Combining new single-dish molecular line observations of 12CO, 13CO, HCO+, HCN, and CS towards the two ULIRGs Arp220 and NGC6240 with existing data in the literature, we have compiled the most extensive molecular line data-sets to date of such galaxies. The data allow us to put strong constraints on the properties of the dense star forming gas in these two systems, and compare the relative usefulness of CS, HCN and HCO+ as tracers of dense gas. In addition, we have build molecular line templates based on our observations, and demonstrate that Herschel/HI-FI will be able to detect the high-J transitions of most of the above molecules in a large sample of ULIRGs out to z<=0.5, assuming Arp220 and NGC6240 are representative of the ULIRG population at these redshifts.

We report deep K-band (18-27GHz) observations with the 100-m Green Bank Telescope of HCN(1-0) line emission towards the two submillimeter-selected galaxies (SMGs) SMMJ02399-0136 (z=2.81) and SMMJ16359+6612 (z=2.52). For both sources we have obtained spectra with channel-to-channel rms noise of <=0.5mJy, resulting in velocity-integrated line fluxes better than < 0.1 Jy km/s, although we do not detect either source. Such sensitive observations -- aided by gravitational lensing of the sources -- permit us to put upper limits of L_HCN(1-0) < 2x10^10 K km/s pc^2 on the intrinsic HCN(1-0) line luminosities of the two SMGs. The far-infrared (FIR) luminosities for all three SMGs with sensitive HCN(1-0) observations to date are found to be consistent with the tight FIR-HCN luminosity correlation observed in Galactic molecular clouds, quiescent spirals and (ultra) luminous infrared galaxies in the local Universe. Thus, the observed HCN luminosities remain in accordance with what is expected from the universal star formation efficiency per dense molecular gas mass implied by the aforementioned correlation, and more sensitive observations with today's large aperture radio telescopes hold the promise of detecting HCN(1-0) emission in similar objects in the distant Universe.

We present models of the evolution of the gaseous and stellar content of galaxies incorporating the formation of H_2 out of HI gas as part of such a model. We do so by formulating a subgrid model for gas clouds that uses well-known cloud scaling relations and solves for the HI-H_2 balance set by the H_2 formation on dust grains and its FUV-induced photodissociation by the temporally and spatially varying interstellar radiation field. This allows the seamless tracking of the evolution of the H_2 gas phase, its precursor Cold Neutral Medium (CNM) HI gas, simultaneously with the star formation. Our most important findings are: a) a significant dependence of the HI-H_2 transition and the resultant H_2 gas mass on the ambient metallicity and the H_2 formation rate, b) the important influence of the characteristic star formation timescale (regulating the ambient FUV radiation field) on the equilibrium H_2 gas mass and c) the possibility of a diffuse H_2 gas phase. Finally, we implement and briefly explore a novel approach of using the ambient H_2 gas mass fraction as a criterion for the onset of star formation. (abridged)

Preliminary results from a sensitive survey of the CO J=1--0, 2--1, 3--2, 4--3, 6--5, HCN J=1--0, 3--2, 4--3, and CI J=1--0 lines of a sample of 30 Ultraluminous Infrared Galaxies (ULIRGs) are presented. These reveal a tandalizing picture of the physical conditions of the molecular gas in these extraordinary galaxies ($\rm L_{FIR}>10^{12} L_{\odot}$), with a diffuse phase dominating the low-J CO lines and a much denser and warmer phase dominating the CO 4-3 and 6-5 and all the HCN lines. The CI J=1--0 emission was found to be a robust tracer of their total molecular gas mass under a large range of physical conditions, a potent alternative to the much weaker emission from the $ ^{13}$CO isotopologue, and especially promising as an H$_2$ tracer for similar objects at high redshifts.

We report sensitive imaging observations of the CO J=1-0 line emission in the powerful high-redshift radio galaxy 4C 41.17 (z=3.8) with the NRAO Very Large Array (VLA), conducted in order to detect the large concomitant H_2 gas reservoir recently unveiled in this system by De Breuck et al (2005) via the emission of the high excitation J=4-3 line. Our observations fail to detect the J=1-0 line but yield sensitive lower limits on the R_43=(4-3)/(1-0) brightness temperature ratio of R_43 ~ 0.55 - >1.0 for the bulk of the H_2 gas mass. Such high ratios are typical of the high-excitation molecular gas phase ``fueling'' the star formation in local starbursts, but quite unlike these objects, much of the molecular gas in 4C 41.17 seems to be in such a state, and thus participating in the observed starburst episode. The widely observed and unique association of highly excited molecular gas with star forming sites allows CO line emission with large (high-J)/(low-J) intensity ratios to serve as an excellent ``marker'' of the spatial distribution of star formation in distant dust-obscured starbursts, unaffected by extinction.

Sensitive new observations of the fine structure line $^3$P$_2$$\to $ $^3$P$_1$ (J=2--1) of the neutral atomic carbon CI ($\nu_{rest}\sim 809$ GHz) in the strongly lensed Ultra Luminous Infrared Galaxy (ULIRG) IRAS F10214+4724 at z=2.3 with the mm/sub-mm telescope James Clerk Maxwel (JCMT) are presented. These do not confirm the presence of emission from this line at the flux levels or angular extent previously reported in the literature. The new 2$\sigma $ upper limits are: $\rm S_{CI}\la 7 Jy km s^{-1}$ (central position), and $\rm < S_{CI} > \la 8.5 Jy km s^{-1}$ (average over the two $\rm [\delta (RA), \delta (Dec)]=[0'',\pm 10'']$ positions). A CI emission assumed fully concomitant with the bulk of H$_2$ and confined entirely within the strongly lensed object yields an upper limit of $\rm M_{CI}(H_2)\la 1.5 \times 10^{10} M_{\odot}$, compatible with the reported CO-derived H$_2$ gas mass, within the uncertainties of the two methods. A comparison with the recent detection of the $^3$P$_1$$\to $ $^3$P$_0$ (J=1--0) line in this galaxy by Weiss et al. (2004) is made and the large discrepancy with the previous CI measurements is briefly discussed.

We present new sensitive wide-band measurements of the fine structure line 3^P_1 -> 3^P_0 (J=1-0, 492GHz) of neutral atomic carbon (CI) in the two typical Ultra Luminous Infrared Galaxies NGC6240 and Arp220. We then use them along with several other CI measurements in similar objects found in the literature to estimate their global molecular gas content under the assumption of a full CI-H_2 concomitance. We find excellent agreement between the H_2 gas mass estimated with this method and the standard methods using 12^CO. This may provide a new way to measure H_2 gas mass in galaxies, and one which may be very valuable in ULIRGs since in such systems the bright 12^CO emission is known to systematically overestimate the gas mass while their 13^CO emission is usually very weak. At redshifts z>=1 the CI J=1-0 line shifts to much more favorable atmospheric windows and can become a viable alternative tracer of the H_2 gas fueling starburst events in the distant Universe.

I examine the possibility that the Cold Neutral Medium of the interstellar medium in galaxies contains a molecular gas phase that may represent a significant and even the dominant amount of its mass in metal-poor regions. In spiral galaxies such regions are found at large galactocentric distances where diffuse H$_2$ gas will be untraceable through its feeble $ ^{12}$CO J=1--0 emission, the very lack of it being also responsible its higher kinetic temperatures ($\rm T_k\sim 60-100$ K). Sensitive sub-mm imaging of spiral galaxies has demonstrated the existence of dust well inside their HI gas distribution while recent observational work suggests a high H$_2$ formation rate from HI association onto grains. The latter is indeed a critical unknown and a high value can easily compensate for the reduction of the available grain surface in the metal-poor environments giving rise to a CO-poor, diffuse H$_2$ phase that may then contribute significantly to the mass and pressure of the interstellar medium in such environments.

We report on the detection of the lowest CO J = 1 - 0 transition in the powerful high-redshift radio galaxy 4C60.07 at z = 3.79. The CO emission is distributed in two spatially and kinematically distinct components as was previously known from the observations of the higher excitation CO J = 4 - 3 line. The total molecular gas mass in 4C60.07 inferred from the CO J = 1 - 0 emission is M(H_2) = 1.3 x 10^11 M_sun, sufficient to fuel the inferred star-formation rate of ~1600 M_sun yr^-1 for 10^8 yrs. From our high-resolution CO J = 1 - 0 VLA maps we find the dynamical mass of 4C60.07 to be comparable to that of a giant elliptical at the present time. A significant fraction of the mass is in the form of molecular gas suggesting that 4C60.07 is in an early state of its evolution. The merging nature of 4C60.07 along with its large dynamical mass imply that this system is a giant elliptical caught in its formative stages.

We report the first detection of the lowest CO transition in a sub-millimetre bright galaxy and extremely red object (ERO) at z = 1.44 using the Very Large Array. The total J = 1 - 0 line luminosity of ERO J164502+4626.4 is (7+-1) x 10^{10} K km s^{-1} pc^2, which yields a total molecular gas mass of ~6 x 10^{10} Msun. We also present a map of the 850-um continuum emission obtained using SCUBA, from which we infer a far-IR luminosity and dust mass of L_FIR ~ 9 x 10^{12} Lsun and M_d ~ 9 x 10^{8} Msun. We find tentative evidence that the CO and sub-mm dust emission is extended over several tens of kpc. If confirmed by high-resolution imaging, this implies that ERO J164502+4626.4 is not simply a high redshift counterpart of a typical Ultra Luminous Infrared Galaxy (ULIRG).

We present results from axisymmetric stellar core collapse simulations in general relativity. Our hydrodynamics code has proved robust and accurate enough to allow for a detailed analysis of the global dynamics of the collapse. Contrary to traditional approaches based on the 3+1 formulation of the gravitational field equations, our framework uses a foliation based on a family of outgoing light cones, emanating from a regular center, and terminating at future null infinity. Such a coordinate system is well adapted to the study of interesting dynamical spacetimes in relativistic astrophysics such as stellar core collapse and neutron star formation. Perhaps most importantly this procedure allows for the unambiguous extraction of gravitational waves at future null infinity without any approximation, along with the commonly used quadrupole formalism for the gravitational wave extraction. Our results concerning the gravitational wave signals show noticeable disagreement when those are extracted by computing the Bondi news at future null infinity on the one hand and by using the quadrupole formula on the other hand. We have strong indication that for our setup the quadrupole formula on the null cone does not lead to physical gravitational wave signals. The Bondi gravitational wave signals extracted at infinity show typical oscillation frequencies of about 0.5 kHz.

There is now strong evidence suggesting that the 12CO J = 1-0 transition, widely used to trace H2 gas, significantly underestimates its mass in metal-poor regions. In spiral disks such regions are found in large galactocentric distances where we show that any unaccounted H2 gas phase is likely to be diffuse (~5-20 cm^-3) and warmer (T(kin) ~ 50-100 K) than the cool (T(kin) ~ 15-20 K) CO-luminous one. Moreover we find that a high value of the H2 formation rate on grains, suggested by recent observational work, can compensate for the reduction of the available grain surface in the metal-poor part of typical galactic disks and thus enhance this CO-poor H2 component which may be contributing significantly to the mass and pressure of spiral disks beyond their optical radius.

The ultraluminous broad absorption line quasar APM08279+5255 is one of the most luminous systems known. Here, we present an analysis of its nuclear CO(1-0) emission. Its extended distribution suggests that the gravitational lens in this system is highly elliptical, probably a highly inclined disk. The quasar core, however, lies in the vicinity of naked cusp, indicating that APM08279+5255 is truly the only odd-image gravitational lens. This source is the second system for which the gravitational lens can be used to study structure on sub-kpc scales in the molecular gas associated with the AGN host galaxy. The observations and lens model require CO distributed on a scale of $\sim 400$ pc. Using this scale, we find that the molecular gas mass makes a significant, and perhaps dominant, contribution to the total mass within a couple hundred parsecs of the nucleus of APM08279+5255.

We examine the physical conditions of the molecular gas in the extremely red object (ERO), HR10 at z=1.44, as constrained by the recent detection of CO J+1 -> J, J=1, 4 lines. The line ratios of such widely spaced CO transitions are extremely sensitive probes of gas excitation; in the case of the HR10, they provide a rare opportunity to study the state of the molecular gas in a gas-rich starburst at a cosmologically significant distance. Contrary to earlier claims, we find the CO J=5-4 transition to be of low excitation. The brightness temperature ratio of (5-4)/(2-1) ~ 0.15. This has serious consequences for estimates of H2 gas mass and the detectability of high-J CO lines from similar objects at high redshifts (z>~3).

We present new numerical algorithms for the coupled Einstein-perfect fluid system in axisymmetry. Our framework uses a foliation based on a family of light cones, emanating from a regular center, and terminating at future null infinity. This coordinate system is well adapted to the study of the dynamical spacetimes associated with isolated relativistic compact objects such as neutron stars. In particular, the approach allows the unambiguous extraction of gravitational waves at future null infinity and avoids spurious outer boundary reflections. The code can accurately maintain long-term stability of polytropic equilibrium models of relativistic stars. We demonstrate global energy conservation in a strongly perturbed neutron star spacetime, for which the total energy radiated away by gravitational waves corresponds to a significant fraction of the Bondi mass. As a first application we present results in the study of pulsations of axisymmetric relativistic stars, extracting the frequencies of the different fluid modes in fully relativistic evolutions of the Einstein-perfect fluid system and making a first comparison between the gravitational news function and the predicted wave using the approximations of the quadrupole formula.

The non-linear behaviour of oscillation modes in compact stars is a topic of considerable current interest. Accurate numerical studies of such phenomena are likely to require powerful new approaches to both fluid and spacetime computations. We propose that a key ingredient of such methods will be the non-linear evolution of deviations from the background stationary equilibrium star. We investigate the feasibility of this approach by applying it to non-linear radial oscillations of a neutron star, and explore numerically various non-linear features of this problem, for a large range of amplitudes. Quadratic and higher order mode coupling and non-linear transfer of energy is demonstrated and analysed in detail.

We study the interaction of massless scalar fields with self-gravitating neutron stars by means of fully dynamic numerical simulations of the Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to spherical symmetry and the neutron stars are approximated by relativistic polytropes. Studying the nonlinear dynamics of isolated neutron stars is very effectively performed within the characteristic formulation of general relativity, in which the spacetime is foliated by a family of outgoing light cones. We are able to compactify the entire spacetime on a computational grid and simultaneously impose natural radiative boundary conditions and extract accurate radiative signals. We study the transfer of energy from the scalar field to the fluid star. We find, in particular, that depending on the compactness of the neutron star model, the scalar wave forces the neutron star either to oscillate in its radial modes of pulsation or to undergo gravitational collapse to a black hole on a dynamical timescale. The radiative signal, read off at future null infinity, shows quasi-normal oscillations before the setting of a late time power-law tail.

We present the results of numerical simulations of the spherically symmetric gravitational collapse of supermassive stars (SMS). The collapse is studied using a general relativistic hydrodynamics code. The coupled system of Einstein and fluid equations is solved employing observer time coordinates, by foliating the spacetime by means of outgoing null hypersurfaces. The code contains an equation of state which includes effects due to radiation, electrons and baryons, and detailed microphysics to account for electron-positron pairs. In addition energy losses by thermal neutrino emission are included. We are able to follow the collapse of SMS from the onset of instability up to the point of black hole formation. Several SMS with masses in the range $5\times 10^5 M_{\odot}- 10^9 M_{\odot}$ are simulated. In all models an apparent horizon forms initially, enclosing the innermost 25% of the stellar mass. From the computed neutrino luminosities, estimates of the energy deposition by $\nu\bar{\nu}$-annihilation are obtained. Only a small fraction of this energy is deposited near the surface of the star, where, as proposed recently by Fuller & Shi (1998), it could cause the ultrarelativistic flow believed to be responsible for $\gamma$-ray bursts. Our simulations show that for collapsing SMS with masses larger than $5\times 10^5 M_{\odot}$ the energy deposition is at least two orders of magnitude too small to explain the energetics of observed long-duration bursts at cosmological redshifts. In addition, in the absence of rotational effects the energy is deposited in a region containing most of the stellar mass. Therefore relativistic ejection of matter is impossible.

Molecular hydrogen is an important component of galaxies because it fuels star formation and accretion onto AGN, the two processes that generate the large infrared luminosities of gas-rich galaxies. Observations of spectral-line emission from the tracer molecule CO are used to probe the properties of this gas. But the lines that have been studied in the local Universe, mostly the lower rotational transitions of J = 1-0 and J = 2-1, have hitherto been unobservable in high-redshift galaxies. Instead, higher transitions have been used, although the densities and temperatures required to excite these higher transitions may not be reached by much of the gas. As a result, past observations may have underestimated the total amount of molecular gas by a substantial amount. Here we report the discovery of large amounts of low-excitation molecular gas around the infrared-luminous quasar, APM 08279+5255 at z = 3.91, using the two lowest excitation lines of 12CO (J = 1-0 and J = 2-1). The maps confirm the presence of hot and dense gas near the nucleus, and reveal an extended reservoir of molecular gas with low excitation that is 10 to 100 times more massive than the gas traced by higher-excitation observations. This raises the possibility that significant amounts of low-excitation molecular gas may lurk in the environments of high-redshift (z > 3) galaxies.

We present tests and results of a new axisymmetric, fully general relativistic code capable of solving the coupled Einstein-matter system for a perfect fluid matter field. Our implementation is based on the Bondi metric, by which the spacetime is foliated with a family of outgoing light cones. We use high-resolution shock-capturing schemes to solve the fluid equations. The code can accurately maintain long-term stability of a spherically symmetric, relativistic, polytropic equilibrium model of a neutron star. In axisymmetry, we demonstrate global energy conservation of a perturbed neutron star in a compactified spacetime, for which the total energy radiated away by gravitational waves corresponds to a significant fraction of the Bondi mass.

Black holes are superb sources of gravitational wave signals, for example when they are born in stellar collapse. We explore the subtleties that may emerge if mass accretion events increase significantly the mass of the black hole during its gravitational wave emission. We find the familiar damped-oscillatory radiative decay but now both decay rate and frequencies are modulated by the mass accretion rate. Any appreciable increase in the horizon mass during emission reflects on the instantaneous signal frequency, which shows a prominent negative branch in the dot[f](f) evolution diagram. The features of the frequency evolution pattern reveal key properties of the accretion event, such as the total accreted mass and the accretion rate. For slow accretion rates the frequency evolution follows verbatim the accretion rate, as expected from dimensional arguments. In view of the possibility of detection of black hole ``ringing'' by the upcoming gravitational wave experiments, the deciphering of the late time frequency dynamics may provide direct insight into otherwise obscured aspects of the black hole birth process.

We consider exact solutions for static black holes localized on a three-brane in five-dimensional gravity in the Randall-Sundrum scenario. We show that the Reissner-Nordstrom metric is an exact solution of the effective Einstein equations on the brane, re-interpreted as a black hole without electric charge, but with instead a tidal 'charge' arising via gravitational effects from the fifth dimension. The tidal correction to the Schwarzschild potential is negative, which is impossible in general relativity, and in this case only one horizon is admitted, located outside the Schwarzschild horizon. The solution satisfies a closed system of equations on the brane, and describes the strong-gravity regime. Current observations do not strongly constrain the tidal charge, and significant tidal corrections could in principle arise in the strong-gravity regime and for primordial black holes.

We study the stability of three-dimensional numerical evolutions of the Einstein equations, comparing the standard ADM formulation to variations on a family of formulations that separate out the conformal and traceless parts of the system. We develop an implementation of the conformal-traceless (CT) approach that has improved stability properties in evolving weak and strong gravitational fields, and for both vacuum and spacetimes with active coupling to matter sources. Cases studied include weak and strong gravitational wave packets, black holes, boson stars and neutron stars. We show under what conditions the CT approach gives better results in 3D numerical evolutions compared to the ADM formulation. In particular, we show that our implementation of the CT approach gives more long term stable evolutions than ADM in all the cases studied, but is less accurate in the short term for the range of resolutions used in our 3D simulations.

We present sensitive sub-mm imaging of the Seyfert 1 galaxy NGC 7469 at 850 $\mu $m and 450 $\mu $m with the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) and ^{12}CO J=3-2 line observations of its central starbursting region. The global dust spectrum, as constrained by the new set of sub-mm data and available 1.30 mm and IRAS 100 $\mu $m, 60 $\mu $m data reveals a dominant warm dust component with a temperature of $\rm T_{\rm d} \sim 35$ K and a global molecular gas-to-dust ratio $\rm M(H_2)/M_{\rm d}\sim 600$. Such high values are typical for IR-bright spirals and in order to reconcile them with the significantly lower ratio of ~ 100 obtained for the Milky Way a cold dust reservoir, inconspicuous at FIR wavelengths, is usually postulated. We present evidence that, rather than unaccounted cold dust mass, this high ratio suggests an overestimate of $\rm M(H_2)$ from its associated ^{12}CO J=1-0 line luminosity by a factor of ~ 5 when a Milky Way value for this conversion is used in starburst nuclei environments.

Formulations of Eulerian general relativistic ideal hydrodynamics in conservation form are analyzed in some detail, with particular emphasis to geometric source terms. Simple linear transformations of the equations are introduced and the associated equivalence class is exploited for the optimization of such sources. A significant reduction of their complexity is readily possible in generic spacetimes. The local characteristic structure of the standard member of the equivalence class is analyzed for a general equation of state (EOS). This extends previous results restricted to the polytropic case. The properties of all other members of the class, in particular specialized forms employing Killing symmetries, are derivable from the standard form. Special classes of EOS are identified for both spacelike and null foliations, which lead to explicit inversion of the state vector and computational savings. The entire approach is equally applicable to spacelike or lightlike foliations and presents a complete proposal for numerical relativistic hydrodynamics on stationary or dynamic geometries.

We have studied numerically the so-called Bondi-Hoyle (wind) accretion on to a rotating (Kerr) black hole in general relativity. We have used the Kerr-Schild form of the Kerr metric, free of coordinate singularities at the black hole horizon. The `test-fluid' approximation has been adopted, assuming no dynamical evolution of the gravitational field. We have used a recent formulation of the general relativistic hydrodynamic equations which casts them into a first-order hyperbolic system of conservation laws. Our studies have been performed using a Godunov-type scheme based on Marquina's flux-formula. We find that regardless of the value of the black hole spin the final accretion pattern is always stable, leading to constant accretion rates of mass and momentum. The flow is characterized by a strong tail shock which is increasingly wrapped around the central black hole as the hole angular momentum increases. The rotation induced asymmetry in the pressure field implies that besides the well known drag, the black hole will experience also a lift normal to the flow direction.

We report the detection of sub-mm emission from dust at 850 microns and of the 12CO J=4-3 line in the two distant powerful radio galaxies 4C 60.07 (z=3.79) and 6C 1909+722 (z=3.53). In the case of 4C 60.07 the dust emission is also detected at 1.25 mm. The estimated molecular gas masses are large, of the order of ~(0.5-1)x10^{11} Solar. The large FIR luminosities (L_fir ~ 10^{13} Solar) suggest that we are witnessing two major starburst phenomena, while the observed large velocity widths (FWHM > 500 km/sec) are characteristic of mergers. In the case of 4C 60.07 the CO emission extends over ~30 kpc and spans a velocity range of >1000 km/sec. It consists of two distinct features with FWHM of >= 550 km/sec and ~150 km/sec and line centers separated by >=700 km/sec The least massive of these components is probably very gas-rich with potentially >=60% of its dynamical mass in the form of molecular gas. The extraordinary morphology of the CO emission in this object suggests that it is not just a scaled-up version of a local Ultra Luminous Infrared Galaxy, and it may be a formative stage of the elliptical host of the residing radio-loud AGN. Finally we briefly explore the effects of the wide range of gas excitation conditions expected for starburst environments on the luminosity of high-J CO lines. We conclude that in unlensed objects, CO (J+1-->J), J+1>3 lines can be significantly weak with respect to CO J=1-0 and this can hinder their detection even in the presence of substantial molecular gas masses.

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

Characteristic methods show excellent promise in the evolution of single black hole spacetimes. The effective coupling with matter fields may help the numerical exploration of important astrophysical systems such as neutron star black hole binaries. To this end we investigate formalisms for numerical relativistic hydrodynamics which can be adaptable to null (characteristic) foliations of the spacetime. The feasibility of the procedure is demonstrated with one-dimensional results on the evolution of self-gravitating matter accreting onto a dynamical black hole.

We introduce a formulation of Eulerian general relativistic hydrodynamics which is applicable for (perfect) fluid data prescribed on either spacelike or null hypersurfaces. Simple explicit expressions for the characteristic speeds and fields are derived in the general case. A complete implementation of the formalism is developed in the case of spherical symmetry. The algorithm is tested in a number of different situations, predisposing for a range of possible applications. We consider the Riemann problem for a polytropic gas, with initial data given on a retarded/advanced time slice of Minkowski spacetime. We compute perfect fluid accretion onto a Schwarzschild black hole spacetime using ingoing null Eddington-Finkelstein coordinates. Tests of fluid evolution on dynamic background include constant density and TOV stars sliced along the radial null cones. Finally, we consider the accretion of self-gravitating matter onto a central black hole and the ensuing increase in the mass of the black hole horizon.

We present fully sampled ^{12}$CO, ^{13}CO J=2-1, 3-2 maps of the inner $\sim 1'\times 1'$ region of NGC 1068. We combine these measurements with an existing interferometric map of ^{12}CO J=1-0 that includes single dish data and thus contains all the flux present. This allows a reliable estimate of the ^{12}CO (J=3-2)/(J=1-0) ratio at the highest angular resolution currently possible and the use of this sensitive line ratio to probe the physical condition of the molecular gas. We also present two measurements of the faint C ^{18}O J=2--1 emission which confirm earlier measurements of a high $C ^{18}O/ ^{13}CO$ intensity ratio in this galaxy. A simple two-phase model for the gas can account for all the observed line ratios if C ^{18}O emission and part of the ^{13}CO emission arise in a dense spatially concentrated component. The ^{12}CO emission originates from a warmer, diffuse gas phase.The dense gas phase contains the bulk of the molecular gas mass while the diffuse phase may not be virialized leading to an overestimate of molecular gas mass when deduced from the luminosity of the ^{12}CO J=1--0 line and a standard galactic conversion factor.

We report on sensitive sub-mm imaging observations of the prototype Seyfert~2/starburst galaxy NGC 1068 at 850 $\mu $m and 450 $\mu $m using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). We find clear evidence of dust emission associated with the extended HI component which together with the very faint $ ^{12}$CO J=1--0 emission give a gas-to-dust ratio of $\rm M_{\rm gas}/M_{\rm dust} \sim 70-150$. This contrasts with the larger ratio $\rm M_{\rm gas}/M_{\rm dust}\sim 330$ estimated within a galactocentric radius of $r\leq 1.36$ kpc, where the gas is mostly molecular and starburst activity occurs. The large gas-to-dust ratio found for the starburst region is attributed to a systematic overestimate of the molecular gas mass in starburst environments when the luminosity of the $ ^{12}$CO J=1--0 line and a standard galactic conversion factor is used. On the other hand sub-mm imaging proves to be a more powerful tool than conventional CO imaging for revealing the properties of the diffuse $\rm H_2$ that coexists with HI. This molecular gas phase is characterized by low densities ($\rm n(H_2)<10^3$ cm$ ^{-3}$), very faint emission from sub-thermally excited CO, and contains more mass than HI, namely $\rm M(H_2)/M(HI)\sim 5$.

In our program of studying numerically the so-called Bondi-Hoyle accretion in the fully relativistic regime, we present here first results concerning the evolution of matter accreting supersonically onto a rotating (Kerr) black hole. These computations generalize previous results where the non-rotating (Schwarzschild) case was extensively considered. We parametrize our initial data by the asymptotic conditions for the fluid and explore the dependence of the solution on the angular momentum of the black hole. Towards quantifying the robustness of our numerical results, we use two different geometrical foliations of the black hole spacetime, the standard form of the Kerr metric in Boyer-Lindquist coordinates as well as its Kerr-Schild form, which is free of coordinate singularities at the black hole horizon. We demonstrate some important advantages of using such horizon adapted coordinate systems. Our numerical study indicates that regardless of the value of the black hole spin the final accretion pattern is always stable, leading to constant accretion rates of mass and momentum. The flow is characterized by a strong tail shock, which, unlike the Schwarzschild case, is increasingly wrapped around the central black hole as the hole angular momentum increases. The rotation induced asymmetry in the pressure field implies that besides the well known drag, the black hole will experience also a lift normal to the flow direction. This situation exhibits some analogies with the Magnus effect of classical fluid dynamics.

We present a new geometrical approach to the study of accretion flows onto rotating (Kerr) black holes. Instead of Boyer-Lindquist coordinates, the standard choice in all existing numerical simulations in the literature, we employ the simplest example of a horizon adapted coordinate system, the Kerr-Schild coordinates. This choice eliminates boundary ambiguities and unphysical divergent behavior at the event horizon. Computations of Bondi-Hoyle accretion onto extreme Kerr black holes, performed here for the first time, demonstrate the key advantages of this procedure. We argue it offers the best approach to the numerical study of the, observationally, increasingly more accesible relativistic inner region around black holes.

We develop and calibrate a new method for estimating the gravitational radiation emitted by complex motions of matter sources in the vicinity of black holes. We compute numerically the linearized curvature perturbations induced by matter fields evolving in fixed black hole backgrounds, whose evolution we obtain using the equations of relativistic hydrodynamics. The current implementation of the proposal concerns non-rotating holes and axisymmetric hydrodynamical motions. As first applications we study i) dust shells falling onto the black hole isotropically from finite distance, ii) initially spherical layers of material falling onto a moving black hole, and iii) anisotropic collapse of shells. We focus on the dependence of the total gravitational wave energy emission on the flow parameters, in particular shell thickness, velocity and degree of anisotropy. The gradual excitation of the black hole quasi-normal mode frequency by sufficiently compact shells is demonstrated and discussed. A new prescription for generating physically reasonable initial data is discussed, along with a range of technical issues relevant to numerical relativity.

Despite the fact that the Schwarzschild and Kerr solutions for the Einstein equations, when written in standard Schwarzschild and Boyer-Lindquist coordinates, present coordinate singularities, all numerical studies of accretion flows onto collapsed objects have been widely using them over the years. This approach introduces conceptual and practical complications in places where a smooth solution should be guaranteed, i.e., at the gravitational radius. In the present paper, we propose an alternative way of solving the general relativistic hydrodynamic equations in background (fixed) black hole spacetimes. We identify classes of coordinates in which the (possibly rotating) black hole metric is free of coordinate singularities at the horizon, independent of time, and admits a spacelike decomposition. In the spherically symmetric, non-rotating case, we re-derive exact solutions for dust and perfect fluid accretion in Eddington-Finkelstein coordinates, and compare with numerical hydrodynamic integrations. We perform representative axisymmetric computations. These demonstrations suggest that the use of those coordinate systems carries significant improvements over the standard approach, especially for higher dimensional studies.