We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $\alpha_{\rm CO} = 4.1$ and $\alpha_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments.

We present extremely deep upper limits on the radio emission from 4U~1957+11, an X-ray binary that is generally believed to be a persistently accreting black hole that is almost always in the soft state. We discuss a more comprehensive search for Type I bursts than in past work, revealing a stringent upper limit on the burst rate, bolstering the case for a black hole accretor. The lack of detection of this source at the 1.07 $\mu$Jy/beam noise level indicates jet suppression that is stronger than expected even in the most extreme thin disk models for radio jet production -- the radio power here is 1500--3700 times lower than the extrapolation of the hard state radio/X-ray correlation, with the uncertainties depending primarily on the poorly constrained source distance. We also discuss the location and velocity of the source and show that it must have either formed in the halo or with a strong asymmetric natal kick.

We present submillimeter spectra of the (proto-)super star cluster (SSC) candidates in the starbursting center of the nearby galaxy NGC 253 identified by Leroy et al. (2018). The 2.5pc resolution of our ALMA cycle 3 observations approach the size of the SSCs and allows the study of physical and chemical properties of the molecular gas in these sources. In the 14 SSC sources and in the frequency ranges 342.0-345.8 GHz and 353.9-357.7 GHz we detect 55 lines belonging to 19 different chemical species. The SSCs differ significantly in chemical complexity, with the richest clusters showing 19 species and the least complex showing 4 species. We detect HCN isotopologues and isomers (H$^{13}$CN, HC$^{15}$N, H$^{15}$NC), abundant HC$_3$N, SO and S$^{18}$O, SO$_2$, and H$_2$CS. The gas ratios CO/HCN, CO/HCO$^+$ are low, ~1-10, implying high dense gas fractions in the SSCs. Line ratio analyses suggests chemistry consistent with photon-dominated regions and mechanical heating. None of the SSCs near the galaxy center show line ratios that imply an X-ray dominated region, suggesting that heating by any (still unknown) AGN does not play a major role. The gas temperatures are high in most sources, with an average rotational temperature of ~130 K in SO$_2$. The widespread existence of vibrationally excited HCN and HC$_3$N transitions implies strong IR radiation fields, potentially trapped by a greenhouse effect due to high continuum opacities.

We report on high angular resolution polarimetric observations of the nearby radio galaxy M87 using the Very Long Baseline Array at 24 GHz ($\lambda=$1.3 cm) and 43 GHz ($\lambda=$7 mm) in 2017-2018. New images of the linear polarization substructure in the nuclear region are presented, characterized by a two-component pattern of polarized intensity and smooth rotation of the polarization plane around the 43 GHz core. From a comparison with an analogous dataset from 2007, we find that this global polarization pattern remains stable on a time interval of 11 yr, while showing smaller month-scale variability. We discuss the possible Faraday rotation toward the M87 nucleus at centimeter to millimeter wavelengths. These results can be interpreted in a scenario where the observed polarimetric pattern is associated with the magnetic structure in the confining magnetohydrodynamic wind, which also serves as the source of the observed Faraday rotation.

We present observations of CO(3-2) and $^{13}$CO(3-2) emission near the supernebula in the dwarf galaxy NGC 5253, which contains one of the best examples of a potential globular cluster in formation. The 0.3" resolution images reveal an unusual molecular cloud, "Cloud D1", coincident with the radio-infrared supernebula. The ~6-pc diameter cloud has a linewidth, $\Delta$ v = 21.7 km/s, that reflects only the gravitational potential of the star cluster residing within it. The corresponding virial mass is 2.5 x 10$^5$ M$_\odot$. The cluster appears to have a top-heavy initial mass function, with $M_{low}$~1-2 M$_\odot$. Cloud D1 is optically thin in CO(3-2) probably because the gas is hot. Molecular gas mass is very uncertain but constitutes < 35% of the dynamical mass within the cloud boundaries. In spite of the presence of an estimated ~1500-2000 O stars within the small cloud, the CO appears relatively undisturbed. We propose that Cloud D1 consists of molecular clumps or cores, possibly star-forming, orbiting with more evolved stars in the core of the giant cluster.

We present 3, 15, and 33 GHz imaging towards galaxy nuclei and extranuclear star-forming regions using the Karl G. Jansky Very Large Array as part of the Star Formation in Radio Survey. With $3-33$ GHz radio spectra, we measured the spectral indices and corresponding thermal (free-free) emission fractions for a sample of 335 discrete regions having significant detections in at least two radio bands. After removing 14 likely background galaxies, we find that the median thermal fraction at 33 GHz is $92 \pm 0.8\%$ with a median absolute deviation of $11\%$, when a two-component power-law model is adopted to fit the radio spectrum. Limiting the sample to 238 sources that are confidently identified as star-forming regions, and not affected by potential AGN contamination (i.e., having galactocentric radii $r_{\rm G} \geq 250$ pc), results in a median thermal fraction of $93 \pm 0.8 \%$ with a median absolute deviation of $10\%$. We further measure the thermal fraction at 33 GHz for 163 regions identified at 7" resolution to be $94 \pm 0.8 \%$ with a median absolute deviation of $8\%$. Together, these results confirm that free-free emission dominates the radio spectra of star-forming regions on scales up to $\sim$500 pc in normal star-forming galaxies. We additionally find a factor of $\sim$1.6 increase in the scatter of the measured spectral index and thermal fraction distributions as a function of decreasing galactocentric radius. This trend is likely reflective of the continuous star-formation activity occurring in the galaxy centers, resulting a larger contribution of diffuse nonthermal emission relative to star-forming regions in the disk.

Centimeter-wave transitions are important counterparts to the rotational mm-wave transitions usually observed to study gas in star-forming regions. However, given their relative weakness, these transitions have historically been neglected. We present Australia Telescope Compact Array 4cm- and 15mm-band spectral line observations of nine nearby star-forming galaxies in the H75 array configuration. Thirteen different molecular lines are detected across the sample from OH, NH$_3$, H$_2$O, H$_2$CO, and c-C$_3$H$_2$, as well as 18 radio recombination lines (RRLs) in NGC 253. Excited OH $^2\Pi_{3/2}$ absorption is detected towards NGC 253 (J=5/2), NGC 4945 (J=9/2), and Circinus (J=9/2); the latter two represent only the third and fourth extragalactic J=9/2 detections. These lines in Circinus suggest rotation temperatures in excess of 2000 K, and thus it is likely that the populations of OH rotational states are not governed by a Boltzmann distribution. Circinus's OH lines are blueshifted from the systemic velocity by ~35 km s$^{-1}$, while NGC 4945's are redshifted by ~100 km s$^{-1}$. NGC 4945's OH absorption likely indicates infall onto the nucleus. The NH$_3$ (1,1) through (6,6) lines in NGC 4945 display a superposition of emission and absorption similar to that seen in other dense gas tracers. Strong (3,3) emission points towards maser activity. The relative NH$_3$ absorption strengths in NGC 4945 show similar anomalies as in previous studies of Arp 220 (weak (1,1) and strong (5,5) absorption). A trend towards higher LTE electron temperatures with increasing RRL frequency is present in NGC 253, likely indicative of stimulated emission within the nuclear region.

We present 0.15'' (~2.5pc) resolution ALMA CO(3-2) observations of the starbursting center in NGC253. Together with archival ALMA CO(1-0) and CO(2-1) data we decompose the emission into a disk and non-disk component. We find ~7-16% of the CO luminosity to be associated with the non-disk component ($1.2-4.2 \times 10^7$ K km s$^{-1}$ pc$^2$). The total molecular gas mass in the center of NGC253 is $\sim 3.6 \times 10^8$ M$_\odot$ with $\sim 0.5 \times 10^8$ M$_\odot$ (~15%) in the non-disk component. These measurements are consistent across independent mass estimates through three CO transitions. The high-resolution CO(3-2) observations allow us to identify the molecular outflow within the non-disk gas. Using a starburst conversion factor, we estimate the deprojected molecular mass outflow rate, kinetic energy and momentum in the starburst of NGC253. The deprojected molecular mass outflow rate is in the range ~14-39 M$_\odot$ yr$^{-1}$ with an uncertainty of 0.4dex. The large spread arises due to different interpretations of the kinematics of the observed gas while the errors are due to unknown geometry. The majority of this outflow rate is contributed by distinct outflows perpendicular to the disk, with a significant contribution by diffuse molecular gas. This results in a mass loading factor $\eta = \dot{M}_\mathrm{out} / \dot{M}_\mathrm{SFR}$ in the range $\eta \sim 8-20$ for gas ejected out to ~300pc. We find the kinetic energy of the outflow to be $\sim 2.5-4.5 \times 10^{54}$ erg and ~0.8dex typical error which is ~0.1% of the total or ~8% of the kinetic energy supplied by the starburst. The outflow momentum is $4.8-8.7 \times 10^8$ M$_\odot$ km s$^{-1}$ (~0.5dex error) or ~2.5-4% of the kinetic momentum released into the ISM by feedback. The unknown outflow geometry and launching sites are the primary source of uncertainty in this study.

There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.

We describe a next major frontier in observational studies of galaxy evolution and star formation: linking the physical conditions in the cold, star-forming interstellar medium to host galaxy and local environment. The integrated gas content of galaxies has been surveyed extensively over the last decades. The link between environment and cold gas density, turbulence, excitation, dynamical state, and chemical makeup remain far less well understood. We know that these properties do vary dramatically and theoretical work posits a strong connection between the state of the gas, its ability to form stars, and the impact of stellar feedback. A next major step in the field will be to use sensitive cm-, mm-, and submm-wave spectroscopy and high resolution spectroscopic imaging to survey the state of cold gas across the whole local galaxy population. Such observations have pushed the capabilities of the current generation of telescopes. We highlight three critical elements for progress in the next decade: (1) robust support and aggressive development of ALMA, (2) the deployment of very large heterodyne receiver arrays on single dish telescopes, and (3) development of a new interferometric array that dramatically improves on current capabilities at cm- and mm-wavelengths (~ 1-115 GHz).

To trace chemical evolution in a prototypical starburst environment, we spatially resolve the carbon and oxygen isotope ratios across the central molecular zone (full size ~$\sim 600$ pc) in the nearby starburst galaxy NGC 253. We imaged the emission of the optically thin isotopologues $^{13}$CO, C$^{18}$O, C$^{17}$O, $^{13}$C$^{18}$O at a spatial resolution $\sim50$ pc, comparable to the typical size of giant molecular associations. Optical depth effects and contamination of $^{13}$C$^{18}$O by C$_4$H is discussed and accounted for to derive column densities. This is the first extragalactic detection of the double isotopologue $^{13}$C$^{18}$O. Derived isotopic ratios $^{12}$C/$^{13}$C$\sim21\pm6$, $^{16}$O/$^{18}$O$\sim130\pm40$, and $^{18}$O/$^{17}$O$\sim4.5\pm0.8$ differ from the generally adopted values in the nuclei of galaxies. The molecular clouds in the central region of NGC 253 show similar rare isotope enrichment to those within the central molecular zone of the Milky way. This enrichment is attributed to stellar nucleosynthesis. Measured isotopic ratios suggest an enhancement of $^{18}$O as compared to our Galactic center, which we attribute to an extra $^{18}$O injection from massive stars. Our observations show evidence for mixing of distinct gas components with different degrees of processing. We observe an extra molecular component of highly processed gas on top of the already proposed less processed gas being transported to the central region of NGC 253. Such multicomponent nature and optical depth effects may hinder the use of isotopic ratios based on spatially unresolved line to infer the star formation history and/or initial stellar mass function properties in the nuclei of galaxies.

Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larger galaxy and surrounding cluster. This White Paper points out what advances can be achieved in the field by new technologies, concentrating on the zone where jets are accelerated to relativistic speeds and collimated. The ngVLA and new space VLBI missions will give higher angular resolution, sensitivity and fidelity in the radio, penetrating this zone for additional objects and allowing us to resolve fundamental questions over the physics of jet acceleration and collimation. Interferometry in other bands would allow us to probe directly flaring components. We also emphasize the need for polarimetry, which is essential to revealing the role and configuration of magnetic fields.

PAN is a scientific instrument suitable for deep space and interplanetary missions. It can precisely measure and monitor the flux, composition, and direction of highly penetrating particles ($> \sim$100 MeV/nucleon) in deep space, over at least one full solar cycle (~11 years). The science program of PAN is multi- and cross-disciplinary, covering cosmic ray physics, solar physics, space weather and space travel. PAN will fill an observation gap of galactic cosmic rays in the GeV region, and provide precise information of the spectrum, composition and emission time of energetic particle originated from the Sun. The precise measurement and monitoring of the energetic particles is also a unique contribution to space weather studies. PAN will map the flux and composition of penetrating particles, which cannot be shielded effectively, precisely and continuously, providing valuable input for the assessment of the related health risk, and for the development of an adequate mitigation strategy. PAN has the potential to become a standard on-board instrument for deep space human travel. PAN is based on the proven detection principle of a magnetic spectrometer, but with novel layout and detection concept. It will adopt advanced particle detection technologies and industrial processes optimized for deep space application. The device will require limited mass (~20 kg) and power (~20 W) budget. Dipole magnet sectors built from high field permanent magnet Halbach arrays, instrumented in a modular fashion with high resolution silicon strip detectors, allow to reach an energy resolution better than 10\% for nuclei from H to Fe at 1 GeV/n.

The nuclei of most normal galaxies contain supermassive black holes, which can accrete gas through a disk and become active. These Active Galactic Nuclei, AGN, can form jets which are observed on scales from AU to Mpc and from meter wavelengths to TeV gamma ray energies. High resolution radio imaging and multi-wavelength/messenger campaigns are elucidating the conditions under which this happens. Evidence is presented that: AGN jets are formed when the black hole spins and the accretion disk is strongly magnetized, perhaps on account of gas accretingat high latitude beyond the black hole sphere of influence; AGN jets are collimated close to the black hole by magnetic stress associated with a disk wind; higher power jets can emerge from their galactic nuclei in a relativistic, supersonic and proton-dominated state and they terminate in strong, hot spot shocks; lower power jets are degraded to buoyant plumes and bubbles; jets may accelerate protons to EeV energies which contribute to the cosmic ray spectrum and which may initiate pair cascades that can efficiently radiate synchrotron gamma rays; jets were far more common when the universe was a few billion years old and black holes and massive galaxies were growing rapidly; jets can have a major influence on their environments, stimulating and limiting the growth of galaxies. The observational prospects for securing our understanding of AGN jets are bright.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

We test models of starburst driven outflows using observations of the 22.2 GHz H$_2$O and 36.2 GHz class I CH$_3$OH maser lines. We have observed the starburst galaxy NGC 253 using the Karl G. Jansky Very Large Array. We present evidence for entrainment of star-forming dense-molecular gas in the outflow of NGC 253. We also show that H$_2$O masers are associated with forming super star clusters and not with supernova remnants. We detect four new 36 GHz CH$_3$OH masers in the central kpc and show possible evidence for a star-formation origin of two class I CH$_3$OH masers. Such high resolution observations are essential for understanding the origin of these masers.

Outflows from AGB stars enrich the Galactic environment with metals and inject mechanical energy into the ISM. Radio spectroscopy can recover both properties through observations of molecular lines. We present results from SWAG: "Survey of Water and Ammonia in the Galactic Center". The survey covers the entire Central Molecular Zone (CMZ), the inner 3.35deg x 0.9deg (~480 x 130pc) of the Milky Way that contains ~5x10^7 Mo of molecular gas. Although our survey primarily targets the CMZ, we observe across the entire sightline through the Milky Way. AGB stars are revealed by their signature of double peaked 22 GHz water maser lines. They are distinguished by their spectral signatures and their luminosities, which reach up to 10^-7 Lo. Higher luminosities are usually associated with Young Stellar Objects located in CMZ star forming regions. We detect a population of ~600 new water masers that can likely be associated with AGB outflows.

NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.

Gas density is widely believed to play a governing role in star formation. However, the exact role of density in setting the star formation rate remains debated. We also lack a general theory that explains how the gas density distribution in galaxies is set. The primary factor preventing the resolution of these issues is the limited number of observations of the gas density distribution across diverse environments. Centimeter- and millimeter-wave spectroscopy offer the most promising way forward in this field, but the key density-sensitive transitions are faint compared to the capabilities of current telescopes. In this chapter, we describe how a next-generation Very Large Array (ngVLA) represents the natural next step forward in this sensitivity-limited field. Such a facility would provide a crucial link between the `Milky Way' and `Extragalactic' views of star formation and dramatically advance our understanding of the drive and role of gas density in galaxies, building on current results from ALMA, NOEMA, the Green Bank Telescope, and other current facilities working in this area.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO(1-0) and 12CO(2-1) in the central 40" (680 pc) of the nuclear starburst galaxy NGC 253, including its molecular outflow. We measure the ratio of brightness temperature for CO(2-1)/CO(1-0), r_21, in the central starburst and outflow-related features. We discuss how r_21 can be used to constrain the optical depth of the CO emission, which impacts the inferred mass of the outflow and consequently the molecular mass outflow rate. We find r_21 less than or equal to 1 throughout, consistent with a majority of the CO emission being optically-thick in the outflow, as it is in the starburst. This suggests that the molecular outflow mass is 3-6 times larger than the lower limit reported for optically thin CO emission from warm molecular gas. The implied molecular mass outflow rate is 25-50 solar masses per year, assuming that conversion factor for the outflowing gas is similar to our best estimates for the bulk of the starburst. This is a factor of 9-19 times larger than the star formation rate in NGC 253. We see tentative evidence for an extended, diffuse CO(2-1) component.

We have obtained a smooth time series for the Electric Vector Position Angle (EVPA) of the blazar OJ 287 at centimeter wavelengths, by making $\pm n\pi$ adjustments to archival values from 1974 to 2016. The data display rotation reversals in which the EVPA rotates counter-clockwise (CCW) for 180 deg and then rotates clockwise (CW) by a similar amount. The time scale of the rotations is a few weeks to a year, and the scale for a double rotation, including the reversal, is one to three years. We have seen four of these events in 40 years. A model consisting of two successive outbursts in polarized flux density, with EVPAs counter-rotating, superposed on a steady polarized jet, can explain many of the details of the observations. Polarization images support this interpretation. The model can also help to explain similar events seen at optical wavelengths. The outbursts needed for the model can be generated by the super-magnetosonic jet model of Nakamura et al. (2010) and Nakamura and Meier (2014), which requires a strong helical magnetic field. This model produces forward and reverse pairs of fast and slow MHD waves, and the plasma inside the two fast/slow pairs rotates around the jet axis, but in opposite directions.

The Survey of Water and Ammonia in Nearby galaxies (SWAN) studies atomic and molecular species across the nuclei of four star forming galaxies: NGC\,253, IC\,342, NGC\,6946, and NGC\,2146. As part of this survey, we present Karl G. Jansky Very Large Array (VLA) molecular line observations of three galaxies: IC\,342, NGC\,6946 and NGC\,2146. NGC\,253 is covered in a previous paper. These galaxies were chosen to span an order of magnitude in star formation rates and to select a variety of galaxy types. We target the metastable transitions of ammonia NH$_{3}$(1,1) to (5,5), the 22\,GHz water (H$_2$O) ($6_{16}-5_{23}$) transition, and the 36.1\,GHz methanol (CH$_3$OH) ($4_{-1}-3_{0}$) transition. {We use the NH$_{3}$\ metastable lines to perform thermometry of the dense molecular gas.} We show evidence for uniform heating across the central kpc of IC\,342 with two temperature components for the molecular gas, similar to NGC 253,} of 27\,K and 308\,K, and that the dense molecular gas in NGC\,2146 has a temperature $<$86 K. We identify two new water masers in IC\,342, and one new water maser in each of NGC\,6946 and NGC\,2146. The two galaxies NGC\,253 and NGC\,2146, with the most vigorous star formation, host H$_2$O kilomasers. Lastly, we detect the first 36\,GHz CH$_3$OH\ masers in IC\,342 and NGC\,6946. For the four external galaxies the total CH$_3$OH\ luminosity in each galaxy suggests a correlation with galactic star formation rate, whereas the morphology of the emission is similar to that of HNCO, a weak shock tracer.

We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.

Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While $^{12}$CO has been targeted extensively, isotopologues such as $^{13}$CO have the advantage of being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new dataset of $^{13}$CO(1-0) observations with the IRAM 30-m telescope of the full discs of 9 nearby spiral galaxies from the EMPIRE survey at a spatial resolution of $\sim$1.5kpc. $^{13}$CO(1-0) is mapped out to $0.7-1r_{25}$ and detected at high signal-to-noise throughout our maps. We analyse the $^{12}$CO(1-0)-to-$^{13}$CO(1-0) ratio ($\Re$) as a function of galactocentric radius and other parameters such as the $^{12}$CO(2-1)-to-$^{12}$CO(1-0) intensity ratio, the 70-to-160$\mu$m flux density ratio, the star-formation rate surface density, the star-formation efficiency, and the CO-to-H$_2$ conversion factor. We find that $\Re$ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favored explanations for the observed $\Re$ variations, while abundance changes may also be at play. We calculate a spatially-resolved $^{13}$CO(1-0)-to-H$_2$ conversion factor and find an average value of $1.0\times10^{21}$ cm$^{-2}$ (K.km/s)$^{-1}$ over our sample with a standard deviation of a factor of 2. We find that $^{13}$CO(1-0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than $^{12}$CO(1-0) in the galaxy centres where the fraction of dense gas is larger.

Centaurus A, with its gas-rich elliptical host galaxy, NGC 5128, is the nearest radio galaxy at a distance of 3.8 Mpc. Its proximity allows us to study the interaction between an active galactic nucleus, radio jets, and molecular gas in great detail. We present ALMA observations of low J transitions of three CO isotopologues, HCN, HCO$^{+}$, HNC, CN, and CCH toward the inner projected 500 pc of NGC 5128. Our observations resolve physical sizes down to 40 pc. By observing multiple chemical probes, we determine the physical and chemical conditions of the nuclear interstellar medium of NGC 5128. This region contains molecular arms associated with the dust lanes and a circumnuclear disk (CND) interior to the molecular arms. The CND is approximately 400 pc by 200 pc and appears to be chemically distinct from the molecular arms. It is dominated by dense gas tracers while the molecular arms are dominated by $^{12}$CO and its rare isotopologues. The CND has a higher temperature, elevated CN/HCN and HCN/HNC intensity ratios, and much weaker $^{13}$CO and C$^{18}$O emission than the molecular arms. This suggests an influence from the AGN on the CND molecular gas. There is also absorption against the AGN with a low velocity complex near the systemic velocity and a high velocity complex shifted by about 60 km s$^{-1}$. We find similar chemical properties between the CND in emission and both the low and high velocity absorption complexes implying that both likely originate from the CND. If the HV complex does originate in the CND, then that gas would correspond to gas falling toward the supermassive black hole.

We present 33 GHz imaging for 112 pointings towards galaxy nuclei and extranuclear star-forming regions at $\approx$2" resolution using the Karl G. Jansky Very Large Array (VLA) as part of the Star Formation in Radio Survey. A comparison with 33 GHz Robert C. Byrd Green Bank Telescope single-dish observations indicates that the interferometric VLA observations recover $78\pm4 %$ of the total flux density over 25" regions ($\approx$ kpc-scales) among all fields. On these scales, the emission being resolved out is most likely diffuse non-thermal synchrotron emission. Consequently, on the $\approx30-300$ pc scales sampled by our VLA observations, the bulk of the 33 GHz emission is recovered and primarily powered by free-free emission from discrete HII regions, making it an excellent tracer of massive star formation. Of the 225 discrete regions used for aperture photometry, 162 are extranuclear (i.e., having galactocentric radii $r_{\rm G} \geq 250$ pc) and detected at $>3\sigma$ significance at 33 GHz and in H$\alpha$. Assuming a typical 33 GHz thermal fraction of 90 %, the ratio of optically-thin 33 GHz-to-uncorrected H$\alpha$ star formation rates indicate a median extinction value on $\approx30-300$ pc scales of $A_{\rm H\alpha} \approx 1.26\pm0.09$ mag with an associated median absolute deviation of 0.87 mag. We find that 10 % of these sources are "highly embedded" (i.e., $A_{\rm H\alpha}\gtrsim3.3$ mag), suggesting that on average HII regions remain embedded for $\lesssim1$ Myr. Finally, we find the median 33 GHz continuum-to-H$\alpha$ line flux ratio to be statistically larger within $r_{\rm G}<250$ pc relative the outer-disk regions by a factor of $1.82\pm0.39$, while the ratio of 33 GHz-to-24 $\mu$m flux densities are lower by a factor of $0.45\pm0.08$, which may suggest increased extinction in the central regions.

The Galactic Center contains large amounts of molecular and ionized gas as well as a plethora of energetic objects. Water masers are an extinction-insensitive probe for star formation and thus ideal for studies of star formation stages in this highly obscured region. With the Australia Telescope Compact Array, we observed 22 GHz water masers in the entire Central Molecular Zone with sub-parsec resolution as part of the large SWAG survey: ``Survey of Water and Ammonia in the Galactic Center''. We detect of order 600 22 GHz masers with isotropic luminosities down to ~10^-7 Lo. Masers with luminosities of >~10^-6 Lo are likely associated with young stellar objects. They appear to be close to molecular gas streamers and may be due to star formation events that are triggered at pericenter passages near Sgr A*. Weaker masers are more widely distributed and frequently show double line features, a tell-tale sign for an origin in evolved star envelopes.

We report ALMA observations of $^{12}$CO(3-2) and $^{13}$CO(3-2) in the gas-poor dwarf galaxy NGC 5253. These 0.3"(5.5 pc) resolution images reveal small, dense molecular gas clouds that are located in kinematically distinct, extended filaments. Some of the filaments appear to be falling into the galaxy and may be fueling its current star formation. The most intense CO(3-2) emission comes from the central $\sim$100 pc region centered on the luminous radio-infrared HII region known as the supernebula. The CO(3-2) clumps within the starburst region are anti-correlated with H$\alpha$ on $\sim$5 pc scales, but are well-correlated with radio free-free emission. Cloud D1, which enshrouds the supernebula, has a high $^{12}$CO/$^{13}$CO ratio, as does another cloud within the central 100 pc starburst region, possibly because the clouds are hot. CO(3-2) emission alone does not allow determination of cloud masses as molecular gas temperature and column density are degenerate at the observed brightness, unless combined with other lines such as $^{13}$CO.

The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.

The paradigm in which magnetic fields play a crucial role in launching/collimating outflows in many astrophysical objects continues to gain support. However, semi-analytical models including the effect of magnetic fields on the dynamics and morphology of jets are still missing due to the intrinsic difficulties in integrating the equations describing a collimated, relativistic flow in the presence of gravity. Only few solutions have been found so far, due to the highly nonlinear character of the equations together with the need to blindly search for singularities. These numerical problems prevented a full exploration of the parameter space. We present a new integration scheme to solve r-self-similar, stationary, axisymmetric magnetohydrodynamics equations describing collimated, relativistic outflows crossing smoothly all the singular points (Alfven point and modified slow/fast points). For the first time, we are able to integrate from the disk mid-plane to downstream of the modified fast point. We discuss an ensemble of jet solutions, emphasising trends and features that can be compared to observables. We present, for the first time with a semi-analytical MHD model, solutions showing counter-rotation of the jet for a substantial fraction of its extent. We find diverse jet configurations with bulk Lorentz factors up to 10 and potential sites for recollimation between $10^3-10^7$ gravitational radii. Such extended coverage of the intervals of quantities, such as magnetic-to-thermal energy ratios at the base or the heights/widths of the recollimation region, makes our solutions suitable for application to many different systems where jets are launched.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs.

We present Very Large Array observations of the 33 GHz radio continuum emission from 22 local ultraluminous and luminous infrared (IR) galaxies (U/LIRGs). These observations have spatial (angular) resolutions of 30--720 pc (0.07"-0.67") in a part of the spectrum that is likely to be optically thin. This allows us to estimate the size of the energetically dominant regions. We find half-light radii from 30 pc to 1.7 kpc. The 33 GHz flux density correlates well with the IR emission, and we take these sizes as indicative of the size of the region that produces most of the energy. Combining our 33 GHz sizes with unresolved measurements, we estimate the IR luminosity and star formation rate per area, and the molecular gas surface and volume densities. These quantities span a wide range (4 dex) and include some of the highest values measured for any galaxy (e.g., $\mathrm{\Sigma_{SFR}^{33GHz} \leq 10^{4.1} M_{\odot} yr^{-1} kpc^{-2}}$). At least $13$ sources appear Compton thick ($\mathrm{N_{H}^{33GHz} \geq 10^{24} cm^{-2}}$). Consistent with previous work, contrasting these data with observations of normal disk galaxies suggests a nonlinear and likely multi-valued relation between SFR and molecular gas surface density, though this result depends on the adopted CO-to-H$_{2}$ conversion factor and the assumption that our 33 GHz sizes apply to the gas. 11 sources appear to exceed the luminosity surface density predicted for starbursts supported by radiation pressure and supernovae feedback, however we note the need for more detailed observations of the inner disk structure. U/LIRGs with higher surface brightness exhibit stronger [{\sc Cii}] 158$\mu$m deficits, consistent with the suggestion that high energy densities drive this phenomenon.

We present Karl G Jansky Very Large Array molecular line observations of the nearby starburst galaxy NGC 253, from SWAN: "Survey of Water and Ammonia in Nearby galaxies". SWAN is a molecular line survey at centimeter wavelengths designed to reveal the physical conditions of star forming gas over a range of star forming galaxies. NGC 253 has been observed in four 1GHz bands from 21 to 36 GHz at 6" ($\sim100$pc) spatial and 3.5 km s$^{-1}$ spectral resolution. In total we detect 19 transitions from seven molecular and atomic species. We have targeted the metastable inversion transitions of ammonia (NH$_{3}$) from (1,1) to (5,5) and the (9,9) line, the 22.2 GHz water (H$_2$O) ($6_{16}-5_{23}$) maser, and the 36.1 GHz methanol (CH$_3$OH) ($4_{-1}-3_{0}$) maser. Utilizing NH$_{3}$ as a thermometer, we present evidence for uniform heating over the central kpc of NGC 253. The molecular gas is best described by a two kinetic temperature model with a warm 130K and a cooler 57K component. A comparison of these observations with previous ALMA results suggests that the molecular gas is not heated in photon dominated regions or shocks. It is possible that the gas is heated by turbulence or cosmic rays. In the galaxy center we find evidence for NH$_{3}$(3,3) masers. Furthermore we present velocities and luminosities of three water maser features related to the nuclear starburst. We partially resolve CH$_3$OH masers seen at the edges of the bright molecular emission, which coincides with expanding molecular superbubbles. This suggests that the masers are pumped by weak shocks in the bubble surfaces.

We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).

Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions.

We present a detailed study of a molecular outflow feature in the nearby starburst galaxy NGC 253 using ALMA. We find that this feature is clearly associated with the edge of NGC 253's prominent ionized outflow, has a projected length of ~300 pc, with a width of ~50 pc and a velocity dispersion of ~40 km s^-1, consistent with an ejection from the disk about 1 Myr ago. The kinematics of the molecular gas in this feature can be interpreted (albeit not uniquely) as accelerating at a rate of 1 km s^-1 pc^-1. In this scenario, the gas is approaching escape velocity at the last measured point. Strikingly, bright tracers of dense molecular gas (HCN, CN, HCO+, CS) are also detected in the molecular outflow: We measure an HCN(1-0)/CO(1-0) line ratio of ~1/10 in the outflow, similar to that in the central starburst region of NGC 253 and other starburst galaxies. By contrast, the HCN/CO line ratio in the NGC 253 disk is significantly lower (~1/30), similar to other nearby galaxy disks. This strongly suggests that the streamer gas originates from the starburst, and that its physical state does not change significantly over timescales of ~1 Myr during its entrainment in the outflow. Simple calculations indicate that radiation pressure is not the main mechanism for driving the outflow. The presence of such dense material in molecular outflows needs to be accounted for in simulations of galactic outflows.

We use the IRAM Large Program EMPIRE and new high-resolution ALMA data to measure 13CO(1-0)/C18O(1-0) intensity ratios across nine nearby spiral galaxies. These isotopologues of CO are typically optically thin across most of the area in galaxy disks, and this ratio allows us to gauge their relative abundance due to chemistry or stellar nucleosynthesis effects. Resolved 13CO/C18O gradients across normal galaxies have been rare due to the faintness of these lines. We find a mean 13CO/C18O ratio of 6.0$\pm$0.9 for the central regions of our galaxies. This agrees well with results in the Milky Way, but differs from results for starburst galaxies (3.4$\pm$0.9) and ultraluminous infrared galaxies (1.1$\pm$0.4). In our sample, the 13CO/C18O ratio consistently increases with increasing galactocentric radius and decreases with increasing star formation rate surface density. These trends qualitatively agree with expectations for carbon and oxygen isotopic abundance variations due to stellar nucleosynthesis, with a possible effect of fractionation.

High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be \tau ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $\eta$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{\eta_{12}}{\eta_{13}}\,\tau_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.

We present subarcsecond images of submillimeter CO and continuum emission from a local galaxy forming massive star clusters: the blue compact dwarf galaxy II Zw 40. At $\sim$0.4" resolution (20 pc), the CO(3-2), CO(1-0), 3mm and 870${\mu}$m continuum maps illustrate star formation on the scales of individual molecular clouds. Dust contributes about a third of the 870${\mu}$m continuum emission, with free-free accounting for the rest. On these scales, there is not a good correspondence between gas, dust, and free-free emission. Dust continuum is enhanced toward the star-forming region as compared to the CO emission. We suggest that an unexpectedly low and spatially variable gas-to-dust ratio is the result of massive clusters of the starburst.

SWAG ("Survey of Water and Ammonia in the Galactic Center") is a multi-line interferometric survey toward the Center of the Milky Way conducted with the Australia Telescope Compact Array. The survey region spans the entire ~400pc Central Molecular Zone and comprises ~42 spectral lines at pc spatial and sub-km/s spectral resolution. In addition, we deeply map continuum intensity, spectral index, and polarization at the frequencies where synchrotron, free-free, and thermal dust sources emit. The observed spectral lines include many transitions of ammonia, which we use to construct maps of molecular gas temperature, opacity and gas formation temperature (see poster by Nico Krieger et al., this volume). Water masers pinpoint the sites of active star formation and other lines are good tracers for density, radiation field, shocks, and ionization. This extremely rich survey forms a perfect basis to construct maps of the physical parameters of the gas in this extreme environment.

We observe Arp 220, the nearest Ultra-Luminous Infrared Galaxy (ULIRG), over 4 GHz in the K and Ka bands. We provide constraints for the kinematics,morphology, and identify molecular species on scales resolving both nuclei (0.6" or 230 pc). We detect multiple molecular species, including hydroxyl in both cores. We tentatively detect H2O at 21.84 GHz in both nuclei, indicating the likely presence of maser emission. The observed frequency range also contains metastable ammonia transitions from (J,K) = (1,1) to (5,5), as well as the (9,9) inversion line, which, together are a well-known thermometer of dense molecular gas. Furthermore, the non-metastable (4,2) and (10,9) and possibly the (3,1) lines are also detected. We apply a standard temperature analysis to Arp 220. However, the analysis is complicated in that standard LTE assumptions do not hold. There are indications that a substantial fraction of ammonia could be in the non-metastable transitions as opposed to only the metastable ones. Thus, the non-metastable transitions could be essential to constraining the temperature. We compare all of these data to ALMA observations of this source, confirming the outflow previously observed by other tracers in both nuclei.

We present ALMA observations of the inner 1' (1.2 kpc) of the Circinus galaxy, the nearest Seyfert. We target CO (1-0) in the region associated with a well-known multiphase outflow driven by the central active galactic nucleus (AGN). While the geometry of Circinus and its outflow make disentangling the latter difficult, we see indications of outflowing molecular gas at velocities consistent with the ionized outflow. We constrain the mass of the outflowing molecular gas to be 1.5e5 -5.1e6 solar masses, yielding a molecular outflow rate of 0.35-12.3 solar masses per year. The values within this range are comparable to the star formation rate in Circinus, indicating that the outflow indeed regulates star formation to some degree. The molecular outflow in Circinus is considerably lower in mass and energetics than previously-studied AGN-driven outflows, especially given its high ratio of AGN luminosity to bolometric luminosity. The molecular outflow in Circinus is, however, consistent with some trends put forth in Cicone et al. (2014), including a linear relation between kinetic power and AGN luminosity, as well as its momentum rate vs. bolometric luminosity (although the latter places Circinus among the starburst galaxies in that sample). We detect additional molecular species including CN and C17O.

AGN feedback is invoked as one of the most relevant mechanisms that shape the evolution of galaxies. Our goal is to understand the interplay between AGN feedback and the interstellar medium in M51, a nearby spiral galaxy with a modest AGN and a kpc-scale radio jet expanding through the disc of the galaxy. For that purpose, we combine molecular gas observations in the CO(1-0) and HCN(1-0) lines from the Plateau de Bure interferometer with archival radio, X-ray, and optical data. We show that there is a significant scarcity of CO emission in the ionisation cone, while molecular gas emission tends to accumulate towards the edges of the cone. The distribution and kinematics of CO and HCN line emission reveal AGN feedback effects out to r~500pc, covering the whole extent of the radio jet, with complex kinematics in the molecular gas which displays strong local variations. We propose that this is the result of the almost coplanar jet pushing on molecular gas in different directions as it expands; the effects are more pronounced in HCN than in CO emission, probably as the result of radiative shocks. Following previous interpretation of the redshifted molecular line in the central 5" as caused by a molecular outflow, we estimate the outflow rates to be Mdot_H2~0.9Msun/yr and Mdot_dense~0.6Msun/yr, which are comparable to the molecular inflow rates (~1Msun/yr); gas inflow and AGN feedback could be mutually regulated processes. The agreement with findings in other nearby radio galaxies suggests that this is not an isolated case, and probably the paradigm of AGN feedback through radio jets, at least for galaxies hosting low-luminosity active nuclei.

Compiling data from literature and the ALMA archive, we show enhanced HCN(4-3)/HCO$^+$(4-3) and/or HCN(4-3)/CS(7-6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN-enhancement). The number of sample galaxies is significantly increased from our previous work. We expect this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN-to-HCO$^+$ and HCN-to-CS in AGNs as compared to SB galaxies by a factor of a few to even $>$ 10 is a plausible explanation for the submillimeter HCN-enhancement. However, a counter argument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we could not fully discriminate these two scenarios at this moment due to insufficient amount of multi-transition, multi-species data, the former scenario equivalently claims for abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest it is difficult with conventional gas phase X-ray dominated region (XDR) ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral-neutral reactions that are efficiently activated at high temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.

This white paper discusses how a "next-generation" Very Large Array (ngVLA) operating in the frequency range 1-116 GHz could be a groundbreaking tool to study the detailed astrophysics of the "matter cycle" in the Milky Way and other galaxies. If optimized for high brightness sensitivity, the ngVLA would bring detailed microwave spectroscopy and modeling of the full radio spectral energy distribution into regular use as survey tools at resolutions of 0.1- 1 arcseconds. This wavelength range includes powerful diagnostics of density, excitation, and chemistry in the cold ISM, as well as multiple tracers of the rate of recent star formation, the magnetic field, shocks, and properties of the ionized ISM. We highlight design considerations that would make this facility revolutionary in this area, the foremost of which is a large amount of collecting area on ~km-length baselines. We also emphasize the strong case for harnessing the large proposed collecting area of the ngVLA for very long baseline applications as part of the core design. This would allow measurements of the three dimensional space motions of galaxies to beyond the Local Group and mapping of the Milky Way out to the far side of the disk. Finally, we discuss the gains from the proposed combination of very high resolution and sensitivity to thermal emission, which include observing the feeding of black holes and resolving forming protoclusters.

We report Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter Array (ALMA) spectroscopy in the redshifted molecular oxygen (O$_2$) 56.265~GHz and 424.763~GHz transitions from the $z=0.88582$ gravitational lens towards PKS\,1830$-$21. The ALMA non-detection of O$_2$ 424.763~GHz absorption yields the $3\sigma$ upper limit $N({\rm O}_2) \leq 5.8 \times 10^{17}$~cm$^{-2}$ on the O$_2$ column density, assuming that the O$_2$ level populations are thermalized at the gas kinetic temperature of 80~K. The VLA spectrum shows absorption by the CH$_3$CHO 56.185~GHz and 56.265~GHz lines, with the latter strongly blended with the O$_2$ 56.265~GHz line. Since the two CH$_3$CHO lines have the same equilibrium strength, we used the known CH$_3$CHO 56.185~GHz line profile to subtract out the CH$_3$CHO 56.265~GHz feature from the VLA spectrum, and then carried out a search for O$_2$ 56.265~GHz absorption in the residual spectrum. The non-detection of redshifted O$_2$ 56.265~GHz absorption in the CH$_3$CHO-subtracted VLA spectrum yields $N({\rm O}_2) \leq 2.3 \times 10^{17}$~cm$^{-2}$. Our $3\sigma$ limits on the O$_2$ abundance relative to H$_2$ are then $X({\rm O}_2) \leq 9.1 \times 10^{-6}$ (VLA) and $X({\rm O}_2) \leq 2.3 \times 10^{-5}$ (ALMA). These are $5-15$ times lower than the best previous constraint on the O$_2$ abundance in an external galaxy. The low O$_2$ abundance in the $z= 0.88582$ absorber may arise due to its high neutral carbon abundance and the fact that its molecular clouds appear to be diffuse or translucent clouds with low number density and high kinetic temperature.

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $\mu$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

We present Karl G. Jansky Very Large Array (VLA) Ka band (33 GHz) and Atacama Large Millimeter Array (ALMA) Band 3 (94.5 GHz) continuum images covering the nucleus and two extranuclear star-forming regions within the nearby galaxy NGC 3627 (M 66), observed as part of the Star Formation in Radio Survey (SFRS). Both images achieve an angular resolution of $\lesssim$2\arcsec, allowing us to map the radio spectral indices and estimate thermal radio fractions at a linear resolution of $\lesssim$90 pc at the distance of NGC 3627. The thermal fraction at 33 GHz reaches unity at and around the peaks of each HII region; we additionally observed the spectral index between 33 and 94.5 GHz to become both increasingly negative and positive away from the peaks of the HII regions, indicating an increase of non-thermal extended emission from diffusing cosmic-ray electrons and the possible presence of cold dust, respectively. While the ALMA observations were optimized for collecting continuum data, they also detected line emission from the $J=1\rightarrow0$ transitions of HCN and HCO$^{+}$. The peaks of dense molecular gas traced by these two spectral lines are spatially offset from the peaks of the 33 and 94.5 GHz continuum emission for the case of the extranuclear star-forming regions, indicating that our data reach an angular resolution at which one can spatially distinguish sites of recent star formation from the sites of future star formation. Finally, we find trends of decreasing dense gas fraction and velocity dispersion with increasing star formation efficiency among the three regions observed, indicating that the dynamical state of the dense gas, rather than its abundance, plays a more significant role in the star formation process.

We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.

A local dwarf galaxy, NGC 5253, has a young super star cluster that may provide an example of highly efficient star formation. Here we report the detection and imaging, with the Submillimeter Array, of the J= 3-2 rotational transition of CO at the location of the massive cluster associated with the supernebula. The gas cloud is hot, dense, quiescent, and extremely dusty. Its gas-to-dust ratio is lower than the Galactic value, which we attribute to dust enrichment by Wolf-Rayet stars within the embedded star cluster. Its star formation efficiency exceeds 50%, ten times higher than clouds in the Milky Way: this cloud is a factory of stars and soot. We suggest that high efficiency results from the force-feeding of star formation by a streamer of gas falling into the galaxy.

We study the kinematics of ridge lines on the pc-scale jet of the active galactic nucleus BL Lac. We show that the ridge lines display transverse patterns that move superluminally downstream, and that the moving patterns are analogous to waves on a whip. Their apparent speeds $\beta_\mathrm{app}$ (units of $c$) range from 3.9 to 13.5, corresponding to $\beta_\mathrm{wave}^\mathrm{gal}= 0.981 - 0.998$ in the galaxy frame. We show that the magnetic field in the jet is well-ordered with a strong transverse component, and assume that it is helical and that the transverse patterns are Alfv\'en waves propagating downstream on the longitudinal component of the magnetic field. The wave-induced transverse speed of the jet is non-relativistic ($\beta_\mathrm{tr}^\mathrm{gal} \lesssim 0.09$). In 2010 the wave activity subsided and the jet then displayed a mild wiggle that had a complex oscillatory behaviour. The Alfv\'en waves appear to be excited by changes in the position angle of the recollimation shock, in analogy to exciting a wave on a whip by shaking the handle. A simple model of the system with plasma sound speed $\beta_\mathrm{s}=0.3$ and apparent speed of a slow MHD wave $\beta_\mathrm{app,S}=4$ yields Lorentz factor of the beam $\Gamma_\mathrm{beam} \sim 4.5$, pitch angle of the helix (in the beam frame) $\alpha\sim 67^\circ$, Alfv\'en speed $\beta_\mathrm{A}\sim 0.64$, and magnetosonic Mach number $M_\mathrm{ms}\sim 4.7$. This describes a plasma in which the magnetic field is dominant and in a rather tight helix, and Alfv\'en waves are responsible for the moving transverse patterns.

We present spatially resolved ($\sim$50 pc) imaging of molecular gas species in the central kiloparsec of the nearby starburst galaxy NGC 253, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). A total of 50 molecular lines are detected over a 13 GHz bandwidth imaged in the 3 mm band. Unambiguous identifications are assigned for 27 lines. Based on the measured high CO/C$^{17}$O isotopic line ratio ($\gtrsim$350), we show that $^{12}$CO(1-0) has moderate optical depths. A comparison of the HCN and HCO$^{+}$ with their $^{13}$C-substituted isotopologues shows that the HCN(1-0) and HCO$^{+}$(1-0) lines have optical depths at least comparable to CO(1-0). H$^{13}$CN/H$^{13}$CO$^{+}$ (and H$^{13}$CN/HN$^{13}$C) line ratios provide tighter constraints on dense gas properties in this starburst. SiO has elevated abundances across the nucleus. HNCO has the most distinctive morphology of all the bright lines, with its global luminosity dominated by the outer parts of the central region. The dramatic variation seen in the HNCO/SiO line ratio suggests that some of the chemical signatures of shocked gas are being erased in the presence of dominating central radiation fields (traced by C$_{2}$H and CN). High density molecular gas tracers (including HCN, HCO$^+$, and CN) are detected at the base of the molecular outflow. We also detect hydrogen $\beta$ recombination lines that, like their $\alpha$ counterparts, show compact, centrally peaked morphologies, distinct from the molecular gas tracers. A number of sulfur based species are mapped (CS, SO, NS, C$_{2}$S, H$_{2}$CS and CH$_{3}$SH) and have morphologies similar to SiO.

We use ALMA to derive the mass, length, and time scales associated with the nuclear starburst in NGC 253. This region forms ~2 M_sun/yr of stars and resembles other starbursts in scaling relations, with star formation consuming the gas reservoir 10 times faster than in galaxy disks. We present observations of CO, the high effective density transitions HCN(1-0), HCO+(1-0), CS(2-1), and their isotopologues. We identify ten clouds that appear as peaks in line emission and enhancements in the HCN-to-CO ratio. These clouds are massive (~10^7 M_sun) structures with sizes (~30 pc) similar to GMCs in other systems. Compared to disk galaxy GMCs, they show high line widths (~20-40 km/s) given their size, with implied Mach numbers ~90. The clouds also show high surface (~6,000 M_sun/pc^2) and volume densities (n_H2~2,000 cm^-3). Given these, self-gravity can explain the line widths. This short free fall time (~0.7 Myr) helps explain the more efficient star formation in NGC 253. We also consider the starburst region as a whole. The geometry is confused by the high inclination, but simple models support a non-axisymmetric, bar-like geometry with a compact, clumpy region of high gas density embedded in an extended CO distribution. Even for the whole region, the surface density still exceeds that of a disk galaxy GMC. The orbital time (~10 Myr), disk free fall time (<~ 3 Myr), and disk crossing time (<~ 3 Myr) are each much shorter than in a normal spiral galaxy disk. Some but not all aspects of the structure correspond to predictions from assuming vertical dynamical equilibrium or a marginally stable rotating disk. Finally, the CO-to-H2 conversion factor implied by our cloud calculations is approximately Galactic, contrasting with results showing a low value for the whole starburst region. The contrast provides resolved support for the idea of mixed molecular ISM phases in starburst galaxies.

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{\Sigma_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{\Sigma_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

NGC1266 is a nearby lenticular galaxy that harbors a massive outflow of molecular gas powered by the mechanical energy of an active galactic nucleus (AGN). It has been speculated that such outflows hinder star formation (SF) in their host galaxies, providing a form of feedback to the process of galaxy formation. Previous studies, however, indicated that only jets from extremely rare, high power quasars or radio galaxies could impart significant feedback on their hosts. Here we present detailed observations of the gas and dust continuum of NGC1266 at millimeter wavelengths. Our observations show that molecular gas is being driven out of the nuclear region at $\dot{M}_{\rm out} \approx 110 M_\odot$ yr$^{-1}$, of which the vast majority cannot escape the nucleus. Only 2 $M_\odot$ yr$^{-1}$ is actually capable of escaping the galaxy. Most of the molecular gas that remains is very inefficient at forming stars. The far-infrared emission is dominated by an ultra-compact ($\lesssim50$pc) source that could either be powered by an AGN or by an ultra-compact starburst. The ratio of the SF surface density ($\Sigma_{\rm SFR}$) to the gas surface density ($\Sigma_{\rm H_2}$) indicates that SF is suppressed by a factor of $\approx 50$ compared to normal star-forming galaxies if all gas is forming stars, and $\approx$150 for the outskirt (98%) dense molecular gas if the central region is is powered by an ultra-compact starburst. The AGN-driven bulk outflow could account for this extreme suppression by hindering the fragmentation and gravitational collapse necessary to form stars through a process of turbulent injection. This result suggests that even relatively common, low-power AGNs are able to alter the evolution of their host galaxies as their black holes grow onto the M-$\sigma$ relation.

The nearby Sy 1 galaxy NGC 1097 represents an ideal laboratory to explore the molecular chemistry in the presence and surroundings of an active galactic nucleus. Exploring the distribution of different molecular species allows us to understand the physical processes affecting the ISM both in the AGN vicinity as well as in the outer star forming molecular ring. We carried out 3 mm ALMA observations of HCN, HCO+, CCH, CS, HNCO, SiO, HC3N, and SO as well as the 13C isotopologues. All species were imaged over the central 2 kpc (~30") of the galaxy at a resolution of ~2.2"x1.5 (150 pc x 100 pc). HCO+ and CS appear to be slightly enhanced in the star forming ring. CCH, showing the largest variations across NGC 1097, is suggested to be a good tracer of both obscured and early stage star formation. HNCO, SiO and HC3N are significantly enhanced in the inner circumnuclear disk surrounding the AGN. Differences in the molecular abundances are observed between the star forming ring and the inner circumnuclear disk. We conclude that the HCN/HCO+ and HCN/CS differences observed between AGN dominated and starburst galaxies are not due to a HCN enhancement due to X-rays, but rather this enhancement is produced by shocked material at distances of 200 pc from the AGN. Additionally we claim the lower HCN/CS to be a combination of a small under-abundance of CS in AGNs together with excitation effects, where a high dense gas component (~10^6 cm^-3) may be more prominent in SB galaxies. However the most promising are the differences found among the dense gas tracers which, at our modest spatial resolution, seem to outline the physical structure of the molecular disk around the AGN. In this picture, HNCO probes the well shielded gas in the disk, surrounding the dense material moderately exposed to X-ray radiation traced by HC3N. Finally SiO might be the innermost molecule in the disk structure.

We present a survey of 3 mm molecular lines in IRAS 04296+2923, one of the brightest known molecular-line emitting galaxies, and one of the closest LIRGs. Data are from the Owens Valley and CARMA millimeter interferometers. Species detected at ~<4" resolution include C^18O, HCN, HCO+, HNC, CN, CH_3OH and, tentatively, HNCO. Along with existing CO, ^13CO and radio continuum data, these lines constrain the chemical properties of the inner disk. Dense molecular gas in the nucleus fuels a star formation rate ~>10 M_sun/yr and is traced by lines of HCN, HCO+, HNC, and CN. A correlation between HCN and star formation rate is observed on sub-kpc scales, consistent with global relations. Toward the nucleus, CN abundances are similar to those of HCN, indicating emission comes from a collection (~40-50) of moderate visual extinction, photon-dominated region clouds. The CO isotopic line ratios are unusual: CO(1-0)/^13CO(1-0) and CO(1-0)/C^18O(1-0) line ratios are large toward the starburst, as is commonly observed in LIRGs, but farther out in the disk these ratios are remarkably low (~<3). ^13CO/C^18O abundance ratios are lower than in Galactic clouds, possibly because the C^18O is enriched by massive star ejecta from the starburst. ^13CO is underabundant relative to CO. Extended emission from CH_3OH indicates that dynamical shocks pervade both the nucleus and the inner disk. The unusual CO isotopologue ratios, the CO/HCN intensity ratio versus L_IR, the HCN/CN abundance ratio and the gas consumption time versus inflow rate, all indicate that the starburst in IRAS 04296+2923 is in an early stage of development.

We present parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

Parsec-scale VLBA images of BL Lac at 15 GHz show that the jet contains a permanent quasi-stationary emission feature 0.26 mas (0.34 pc projected) from the core, along with numerous moving features. In projection, the tracks of the moving features cluster around an axis at position angle -166.6 deg that connects the core with the standing feature. The moving features appear to emanate from the standing feature in a manner strikingly similar to the results of numerical 2-D relativistic magneto-hydrodynamic (RMHD) simulations in which moving shocks are generated at a recollimation shock. Because of this, and the close analogy to the jet feature HST-1 in M87, we identify the standing feature in BL Lac as a recollimation shock. We assume that the magnetic field dominates the dynamics in the jet, and that the field is predominantly toroidal. From this we suggest that the moving features are compressions established by slow and fast mode magneto-acoustic MHD waves. We illustrate the situation with a simple model in which the slowest moving feature is a slow-mode wave, and the fastest feature is a fast-mode wave. In the model the beam has Lorentz factor about 3.5 in the frame of the host galaxy, and the fast mode wave has Lorentz factor about 1.6 in the frame of the beam. This gives a maximum apparent speed for the moving features 10c. In this model the Lorentz factor of the pattern in the galaxy frame is approximately 3 times larger than that of the beam itself.

We describe a new paradigm for understanding both relativistic motions and particle acceleration in the M87 jet: a magnetically dominated relativistic flow that naturally produces four relativistic magnetohydrodynamic (MHD) shocks (forward/reverse fast and slow modes). We apply this model to a set of optical super- and subluminal motions discovered by Biretta and coworkers with the {\em Hubble Space Telescope} during 1994 -- 1998. The model concept consists of ejection of a {\em single} relativistic Poynting jet, which possesses a coherent helical (poloidal + toroidal) magnetic component, at the remarkably flaring point HST-1. We are able to reproduce quantitatively proper motions of components seen in the {\em optical} observations of HST-1 with the same model we used previously to describe similar features in radio VLBI observations in 2005 -- 2006. This indicates that the quad relativistic MHD shock model can be applied generally to recurring pairs of super/subluminal knots ejected from the upstream edge of the HST-1 complex as observed from radio to optical wavelengths, with forward/reverse fast-mode MHD shocks then responsible for observed moving features. Moreover, we identify such intrinsic properties as the shock compression ratio, degree of magnetization, and magnetic obliquity and show that they are suitable to mediate diffusive shock acceleration of relativistic particles via the first-order Fermi process. We suggest that relativistic MHD shocks in Poynting-flux dominated helical jets may play a role in explaining observed emission and proper motions in many AGNs.

We present maps of a large number of dense molecular gas tracers across the Central Molecular Zone of our Galaxy. The data were taken with the CSIRO/CASS Mopra telescope in Large Projects in the 1.3cm, 7mm, and 3mm wavelength regimes. Here, we focus on the brightness of the shock tracers SiO and HNCO, molecules that are liberated from dust grains under strong (SiO) and weak (HNCO) shocks. The shocks may have occurred when the gas enters the bar regions and the shock differences could be due to differences in the moving cloud mass. Based on tracers of ionizing photons, it is unlikely that the morphological differences are due to selective photo-dissociation of the molecules. We also observe direct heating of molecular gas in strongly shocked zones, with a high SiO/HNCO ratios, where temperatures are determined from the transitions of ammonia. Strong shocks appear to be the most efficient heating source of molecular gas, apart from high energy emission emitted by the central supermassive black hole Sgr A* and the processes within the extreme star formation region Sgr B2.

We present a large-scale, interferometric survey of ammonia (1,1) and (2,2) toward the Galactic Center observed with the Australia Telescope Compact Array (ATCA). The survey covers Delta l ~1degree (~150pc) at an assumed distance of 8.5 kpc) and Delta b ~0.2degree (~30pc) which spans the region between the supermassive black hole SgrA* and the massive star forming region SgrB2. The resolution is ~20'' (~0.8pc) and emission at scales >~2' (>~3.2pc) is filtered out due to missing interferometric short spacings. Consequently, the data represent the denser, compact clouds and disregards the large scale, diffuse gas. Many of the clumps align with the 100 pc dust ring and mostly anti-correlate with 1.2cm continuum emission. We present a kinetic temperature map of the dense gas. The temperature distribution peaks at ~38K with a width at half maximum between 18K and 61K (measurements sensitive within Tkin~10-80K). Larger clumps are on average warmer than smaller clumps which suggests internal heating sources. Our observations indicate that the circumnuclear disk ~1.5 pc around SgrA* is supplied with gas by the 20km/s molecular cloud. This gas is substantially cooler than gas ~3-15pc away from SgrA*. We find a strong temperature gradient across SgrB2. Ammonia column densities correlate well with SCUBA 850um fluxes, but the relation is shifted from the origin, which may indicate a requirement for a minimum amount of dust to form and shield ammonia. Around the Arches and Quintuplet clusters we find shell morphologies with UV-influenced gas in their centers, followed by ammonia and radio continuum layers.

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

We present a new, semi-analytic formalism to model the acceleration and collimation of relativistic jets in a gravitational potential. The gravitational energy density includes the kinetic, thermal, and electromagnetic mass contributions. The solutions are close to self-similar throughout the integration, from very close to the black hole to the region where gravity is unimportant. The field lines are tied to the conditions very close to the central object and eventually overcollimate, possibly leading to a collimation shock. This collimation shock could provide the conditions for diffusive shock acceleration, leading to the observed electron populations with a power-law energy distribution in jets. We provide the derivation, a detailed analysis of a solution, and describe the effects the parameters have on the properties of the solutions, such as the Lorentz factor and location of the collimation shock. We also discuss the deviations from self-similarity. By comparing the new gravity term with the gravity term obtained from a non-relativistic formalism in a previous work, we show they are equivalent in the non-relativistic limit. This equivalence shows the approach taken in that work is valid and allows us to comment on its limitations.

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

The OVRO 40-m telescope has been monitoring the 15 GHz radio flux density of over 1200 blazars since 2008. The 15 GHz light curve of the flat spectrum radio quasar J1359+4011 shows a strong and persistent quasi-periodic oscillation. The time-scale of the oscillation varies between 120 and 150 days over a $\sim4$ year time span. We interpret this as the active galactic nucleus mass-scaled analog of low-frequency quasi-periodic oscillations from Galactic microquasars, or as evidence of modulation of the accretion flow by thermal instabilites in the "inner" accretion disc.

We present a 3-dimensional data cube of the K band continuum and the Br-gamma, H2 S(0) and S(1) lines within the central 18.5" x 13.8" (520 pc x 390 pc) region of NGC 6946. Data were obtained using OSIRIS, a near-infrared Integral Field Unit at Keck Observatory, with Laser Guide Star Adaptive Optics. The 0.3" resolution allows us to investigate the stellar bulge and the forming star clusters in the nuclear region on 10 parsec scales. We detect giant HII regions associated with massive young star clusters in the nuclear spiral/ring (R ~ 30pc) and in the principal shocks along the nuclear bar. Comparisons of the Br-gamma fluxes with Pa-alpha line emission and radio continuum indicate A_K ~ 3, A_V ~ 25 for the nuclear star forming regions. The most luminous HII regions are restricted to within 70pc of the center, despite the presence of high gas columns at larger radii (R ~ 200pc). H2 emission is restricted to clouds within R ~ 60pc of the center, and in this sense resembles the distribution of the HCN line emission. We propose that gas-assisted migration of the young star clusters is contributing to the buildup of the nuclear bar and nuclear star cluster (R < 30pc) in this galaxy.

We present the first 100 pc scale view of the dense molecular gas in the central ~ 1.3 kpc region of the type-1 Seyfert NGC 1097 traced by HCN (J=4-3) and HCO+ (J=4-3) lines afforded with ALMA band 7. This galaxy shows significant HCN enhancement with respect to HCO+ and CO in the low-J transitions, which seems to be a common characteristic in AGN environments. Using the ALMA data, we study the characteristics of the dense gas around this AGN and search for the mechanism of HCN enhancement. We find a high HCN (J=4-3) to HCO+ (J=4-3) line ratio in the nucleus. The upper limit of the brightness temperature ratio of HCN (v2=1^{1f}, J=4-3) to HCN (J=4-3) is 0.08, which indicates that IR pumping does not significantly affect the pure rotational population in this nucleus. We also find a higher HCN (J=4-3) to CS (J=7-6) line ratio in NGC 1097 than in starburst galaxies, which is more than 12.7 on the brightness temperature scale. Combined from similar observations from other galaxies, we tentatively suggest that this ratio appears to be higher in AGN-host galaxies than in pure starburst ones similar to the widely used HCN to HCO+ ratio. LTE and non-LTE modeling of the observed HCN and HCO+ lines using J=4-3 and 1-0 data from ALMA, and J=3-2 data from SMA, reveals a high HCN to HCO+ abundance ratio (5 < [HCN]/[HCO+] < 20: non-LTE analysis) in the nucleus, and that the high-J lines (J=4-3 and 3-2) are emitted from dense (10^{4.5} < n_H2 [/cc] < 10^6), hot (70 < Tkin [K] < 550) regions. Finally we propose that the high temperature chemistry is more plausible to explain the observed enhanced HCN emission in NGC 1097 than the pure gas phase PDR/XDR chemistry.

We report the detection of a 22GHz water maser line in the nearest (D~3.8Mpc) radio galaxy Centaurus A using the Australia Telescope Compact Array (ATCA). The line is centered at a velocity of ~960kms-1, which is redshifted by about 400kms-1 from the systemic velocity. Such an offset, as well as the width of ~120kms-1, could be consistent with either a nuclear maser arising from an accretion disk of the central supermassive black hole, or for a jet maser that is emitted from the material that is shocked near the base of the jet in Centaurus\,A. The best spatial resolution of our ATCA data constrains the origin of the maser feature within <3pc from the supermassive black hole. The maser exhibits a luminosity of ~1Lo, which classifies it as a kilomaser, and appears to be variable on timescales of months. A kilomaser can also be emitted by shocked gas in star forming regions. Given the small projected distance from the core, the large offset from systemic velocity, as well as the smoothness of the line feature, we conclude that a jet maser line emitted by shocked gas around the base of the AGN is the most likely explanation. For this scenario we can infer a minimum density of the radio jet of ~>10cm-3, which indicates substantial mass entrainment of surrounding gas into the propagating jet material.

We present a kinematic analysis of the dense molecular gas in the central 200 parsecs of the nearby galaxy NGC1097, based on Cycle 0 observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA). We use the HCN(4-3) line to trace the densest interstellar molecular gas, and quantify its kinematics, and estimate an inflow rate for the molecular gas. We find a striking similarity between the ALMA kinematic data and the analytic spiral inflow model that we have previously constructed based on ionized gas velocity fields on larger scales. We are able to follow dense gas streaming down to 40 pc distance from the supermassive black hole in this Seyfert 1 galaxy. In order to fulfill marginal stability, we deduce that the dense gas is confined to a very thin disc, and we derive a dense gas inflow rate of 0.09 Msun/yr at 40 pc radius. Combined with previous values from the Ha and CO gas, we calculate a combined molecular and ionized gas inflow rate of 0.2 Msun/yr at 40 pc distance from the central supermassive black hole of NGC1097.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

We present a new, approximate method for modelling the acceleration and collimation of relativistic jets in the presence of gravity. This method is self-similar throughout the computational domain where gravitational effects are negligible and, where significant, self-similar within a flux tube. These solutions are applicable to jets launched from a small region (e.g., near the inner edge of an accretion disk). As implied by earlier work, the flow can converge onto the rotation axis, potentially creating a collimation shock. In this first version of the method, we derive the gravitational contribution to the relativistic equations by analogy with non-relativistic flow. This approach captures the relativistic kinetic gravitational mass of the flowing plasma, but not that due to internal thermal and magnetic energies. A more sophisticated treatment, derived from the basic general relativistic magnetohydrodynamical equations, is currently being developed. Here we present an initial exploration of parameter space, describing the effects the model parameters have on flow solutions and the location of the collimation shock. These results provide the groundwork for new, semi-analytic models of relativistic jets which can constrain conditions near the black hole by fitting the jet break seen increasingly in X-ray binaries.

We have carried out a systematic search with Spitzer Warm Mission and archival data for infrared emission from the hotspots in radio lobes that have been described by Hardcastle et al. (2004). These hotspots have been detected with both radio and X-ray observations, but an observation at an intermediate frequency in the infrared can be critical to distinguish between competing models for particle acceleration and radiation processes in these objects. Between the archival and warm mission data, we report detections of 18 hotspots; the archival data generally include detections at all four IRAC bands, the Warm Mission data only at 3.6 microns. Using a theoretical formalism adopted from Godfrey et al. (2009), we fit both archival and warm mission spectral energy distributions [SEDs] - including radio, X-ray, and optical data from Hardcastle as well as the Spitzer data - with a synchrotron self-Compton [SSC] model, in which the X-rays are produced by Compton scattering of the radio frequency photons by the energetic electrons which radiate them. With one exception, a SSC model requires that the magnetic field be less or much less than the equipartition value which minimizes total energy and has comparable amounts of energy in the magnetic field and in the energetic particles. This conclusion agrees with those of comparable recent studies of hotspots, and with the analysis presented by Hardcastle et al. (2004). We also show that the infrared data rule out the simplest synchrotron only models for the SEDs. We briefly discuss the implications of these results and of alternate interpretations of the data.

We present 2" - 10" imaging of eleven transitions from nine molecular species across the nuclear bar in Maffei 2. The data were obtained with the BIMA and OVRO interferometers. The ten detected transitions are compared with existing CO isotopologues, HCN, CS and millimeter continuum data. Dramatic spatial variations among the mapped species are observed across the nuclear bar. A principle component analysis is performed to characterize correlations between the transitions, star formation and molecular column density. The analysis reveals that HCN, HNC, HCO+ and 3 mm continuum are tightly correlated, indicating a direct connection to massive star formation. We find two main morphologically distinct chemical groups, CH3OH, SiO and HNCO comprising the grain chemistry molecules, versus HCN, HNC, HCO+ and C2H, molecules strong in the presence of star formation. The grain chemistry molecules, HNCO, CH3OH and SiO, trace hydrodynamical bar shocks. The near constancy of the HNCO/CH3OH, SiO/CH3OH and SiO/HNCO ratios argue that shock properties are uniform across the nucleus. HCN/HCO+, HCN/HNC, HCN/CS and HCN/CO ratios are explained primarily by variations in density. High HCO+/N2H+ ratios are correlated with the C2H line, suggesting that this ratio may be a powerful new dense photon-dominated region (PDR) probe in external galaxies. C2H reveals a molecular outflow along the minor axis. The morphology and kinematics of the outflow are consistent with an outflow age of 6-7 Myrs.

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices \alpha \approx -0.67 (where F_\nu \propto \nu^\alpha) with \pm 0.15 (1\sigma) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

We present the first images of the J=5-4 and J=16-15 lines of the dense gas tracer, cyanoacetylene, HC_3N, in an external galaxy. The central 200 pc of the nearby star-forming spiral galaxy, IC 342, was mapped using the VLA and the Plateau de Bure Interferometer. HC_3N(5-4) line emission is found across the nuclear mini-spiral, but is very weak towards the starburst site, the location of the strongest mid-IR and radio emission. The J=16-15 and 10-9 lines are also faint near the large HII region complex, but are brighter relative to the 5-4 line, consistent with higher excitation. The brightest HC_3N emission is located in the northern arm of the nuclear minispiral, 100 pc away from the radio/IR source to the southwest of the nucleus. This location appears less affected by ultraviolet radiation, and may represent a more embedded, earlier stage of star formation. HC_3N excitation temperatures are consistent with those determined from C^{18}O; the gas is dense, 10^{4-5}/cc, and cool, T_K ~< 40 K. So as to not violate limits on the total H_2 mass determined from C^{18}O, at least two dense components are required to model IC 342's giant molecular clouds. These observations suggest that HC_3N(5-4) is an excellent probe of the dense, quiescent gas in galaxies. The high excitation combined with faint emission towards the dense molecular gas at the starburst indicates that it currently lacks large masses of very dense gas. We propose a scenario where the starburst is being caught in the act of dispersing or destroying its dense gas in the presence of the large HII region. This explains the high star formation efficiency seen in the dense component. The little remaining dense gas appears to be in pressure equilibrium with the starburst HII region.

Observations of relativistic jets from black holes systems suggest that particle acceleration often occurs at fixed locations within the flow. These sites could be associated with critical points that allow the formation of standing shock regions, such as the magnetosonic modified fast point. Using the self-similar formulation of special relativistic magnetohydrodynamics by Vlahakis & K\"onigl, we derive a new class of flow solutions that are both relativistic and cross the modified fast point at a finite height. Our solutions span a range of Lorentz factors up to at least 10, appropriate for most jets in X-ray binaries and active galactic nuclei, and a range in injected particle internal energy. A broad range of solutions exists, which will allow the eventual matching of these scale-free models to physical boundary conditions in the analysis of observed sources.

We present a first look at the local LIRG, IRAS04296+2923. This barred spiral, overlooked because of its location in the Galactic plane, is among the half dozen closest LIRGs. More IR-luminous than either M82 or the Antennae, it may be the best local example of a nuclear starburst caused by bar-mediated secular evolution. We present Palomar J and Pa beta images, VLA maps from 20-1.3cm, a Keck LWS image at 11.7mic and OVRO CO(1-0) and ^13CO(1-0), and 2.7 mm continuum images. The J-band image shows a symmetric barred spiral. Two bright, compact mid-IR/radio sources in the nucleus comprise a starburst that is equivalent to 10^5 O7 stars, probably a pair of young super star clusters separated by 30pc. The nuclear starburst is forming stars at the rate of ~12Msun/yr, half of the total star formation rate for the galaxy of ~25Msun/yr. IRAS04296 is bright in CO, and among the most gas-rich galaxies in the local universe. The CO luminosity of the inner half kpc is equivalent to that of the entire Milky Way. While the most intense CO emission extends over a 15"(2 kpc) region, the nuclear starburst is confined to ~1-2"(150-250 pc) of the dynamical center. From ^13CO, we find that the CO conversion factor in the nucleus is higher than the Galactic value by a factor 3-4, typical of gas-rich spiral nuclei. The nuclear star formation efficiency is M_gas/SFR^nuc = 2.7x10^-8 yr^-1, corresponding to gas consumption timescale, tau_SF^nuc~4x10^7 yrs. The star formation efficiency is ten times lower in the disk, tau_SF^disk~3.3x10^8 yrs. The low absolute star formation efficiency in the disk implies that the molecular gas is not completely consumed before it drifts into the nucleus, and is capable of fueling a sustained nuclear starburst. IRAS04296 is beginning a 100Myr period as a LIRG, during which it will turn much of its 6x10^9Msun of molecular gas into a nuclear cluster of stars. (abridged)

This is the first in a series of papers that introduces a new paradigm for understanding the jet in M87: a collimated relativistic flow in which strong magnetic fields play a dominant dynamical role. Here wefocus on the flow downstream of HST-1 - an essentially stationary flaring feature that ejects trails of superluminal components. We propose that these components are quad relativistic magnetohydrodynamic shock fronts (forward/reverse fast and slow modes) in a narrow jet with a helically twisted magnetic structure. And we demonstrate the properties of such shocks with simple one-dimensional numerical simulations. Quasi-periodic ejections of similar component trails may be responsible for the M87 jet substructures observed further downstream on 100 - 1,000 pc scales. This new paradigm requires the assimilation of some new concepts into the astrophysical jet community, particularly the behavior of slow/fast-mode waves/shocks and of current-driven helical kink instabilities. However, the prospects of these ideas applying to a large number of other jet systems may make this worth the effort.

The Nuclear Spectroscopic Telescope Array (NuSTAR) is a NASA Small Explorer mission that will carry the first focusing hard X-ray (5 -- 80 keV) telescope to orbit. NuSTAR will offer a factor 50 -- 100 sensitivity improvement compared to previous collimated or coded mask imagers that have operated in this energy band. In addition, NuSTAR provides sub-arcminute imaging with good spectral resolution over a 12-arcminute field of view. After launch, NuSTAR will carry out a two-year primary science mission that focuses on four key programs: studying the evolution of massive black holes through surveys carried out in fields with excellent multiwavelength coverage, understanding the population of compact objects and the nature of the massive black hole in the center of the Milky Way, constraining explosion dynamics and nucleosynthesis in supernovae, and probing the nature of particle acceleration in relativistic jets in active galactic nuclei. A number of additional observations will be included in the primary mission, and a guest observer program will be proposed for an extended mission to expand the range of scientific targets. The payload consists of two co-aligned depth-graded multilayer coated grazing incidence optics focused onto solid state CdZnTe pixel detectors. To be launched in early 2012 on a Pegasus rocket into a low-inclination Earth orbit. Data will be publicly available at GSFC's High Energy Astrophysics Science Archive Research Center (HEASARC) following validation at the science operations center located at Caltech.

The fundamental role played by black holes in our study of microquasars, gamma ray bursts, and the outflows from active galactic nuclei requires an appreciation for, and at times some in-depth analysis of, curved spacetime. We highlight misconceptions surrounding the notion of coordinate transformation in general relativity as applied to metrics for rotating black holes that are beginning to increasingly appear in the literature. We emphasize that there is no coordinate transformation that can turn the metric of a rotating spacetime into that for a Schwarzschild spacetime, or more generally, that no coordinate transformation exists that can diagonalize the metric for a rotating spacetime. We caution against the notion of "local" coordinate transformation, which is often incorrectly associated with a global analysis of the spacetime.

Explorers have made breakthroughs in many fields of astrophysics. The science from both these missions contributed to three Nobel Prizes - Giacconi (2002), Mather, and Smoot (2006). Explorers have: marked the definitive beginning of precision cosmology, discovered that short gamma-ray bursts are caused by compact star mergers and have measured metalicity to redshifts z>6. NASA Explorers do cutting-edge science that cannot be done by facility-class instruments. The Explorer program provides a rapid response to changing science and technology, to enable cutting-edge science at moderate cost. Explorers also enable innovation, and engage & train scientists, managers and engineers, adding human capital to NASA and the nation. The astrophysics Explorer launch rate now being achieved is 1 per 3 years, and budget projections are in the $150M/year range for the next five years. A newly Vigorous Explorer Program should be created to: 1. Reach the long-stated goal of annual astrophysics launches; 2. Find additional launch options for Explorers and actively encourage cost savings in launchers and spacecraft, such as new commercial vehicles and innovative partnerships. 3. Mitigate risk via stronger technical development and sub-orbital programs, and through longer, more thorough, Phase A programs, potentially reducing the need for a 30% contingency; 4. Strive to protect the funding for missions that have reached Phase B, to prevent significant launch slips and cancellations, with a goal of 4 to 5 years from Phase B to launch; 5. Review the project management procedures and requirements to seek cost reductions, including the risk management strategy and the review and reporting process; 6. Review and possibly modify the cost caps for all Explorer classes to optimize scientific returns per dollar. [ABRIDGED]

We present 7mm Very Large Array continuum images of the starburst galaxy M82. On arcsecond scales, two-thirds of the 7mm continuum consists of free-free emission from HII regions. In the subarcsecond resolution map, we identify 14 compact sources, including 9 bright HII regions with N_Lyc > 10^51 sec^-1. Four of the HII regions have rising spectra, implying emission measures > 10^8 cm^-6 pc. Except for one compact source with peculiar features, all other compact radio sources are found in dust lanes and do not have optical or near-infrared continuum counterparts. Four regions of extended, high brightness (EM > 10^7 cm-6 pc) radio emission are found in our high resolution map, including some as large as ~2", or 30 pc, representing either associations of small HII regions, or sheetlike structures of denser gas. The good correlation between 7 mm emission and Spitzer IRAC 8 micron continuum-removed PAH feature suggests that PAH emission may track the recently formed OB stars. We find an excellent correlation between molecular gas and star formation, particularly dense gas traced by HCN, down to the ~ 45 pc scale in M82. We also find star formation efficiencies (SFEs) of 1-10% on the same scale, based on CO maps. The highest SFE are found in regions with the highest dense gas fractions.

(Abridged) We present one of the first physically-motivated two-dimensional general relativistic magnetohydrodynamic (GRMHD) numerical simulations of a radiatively-cooled black-hole accretion disk. The fiducial simulation combines a total-energy-conserving formulation with a radiative cooling function, which includes bremsstrahlung, synchrotron, and Compton effects. By comparison with other simulations we show that in optically thin advection-dominated accretion flows, radiative cooling can significantly affect the structure, without necessarily leading to an optically thick, geometrically thin accretion disk. We further compare the results of our radiatively-cooled simulation to the predictions of a previously developed analytic model for such flows. For the very low stress parameter and accretion rate found in our simulated disk, we closely match a state called the "transition" solution between an outer advection-dominated accretion flow and what would be a magnetically-dominated accretion flow (MDAF) in the interior. The qualitative and quantitative agreement between the numerical and analytic models is quite good, with only a few well-understood exceptions. According to the analytic model then, at significantly higher stress or accretion, we would expect a full MDAF to form. The collection of simulations in this work also provide important data for interpreting other numerical results in the literature, as they span the most common treatments of thermodynamics, including simulations evolving: 1) the internal energy only; 2) the internal energy plus an explicit cooling function; 3) the total energy without cooling; and 4) total energy including cooling. We find that the total energy formulation is a necessary prerequisite for proper treatment of radiative cooling in MRI accretion flows.

Using new, high-resolution interferometric observations of the CO and HCN molecules, we directly compare the molecular and ionized components of the interstellar medium in the center of the nearby spiral galaxy IC342, on spatial scales of ~ 10pc. The morphology of the tracers suggests that the molecular gas flow caused by a large-scale stellar bar has been strongly affected by the mechanical feedback from recent star formation activity within the central 100pc in the nucleus of the galaxy. Possibly, stellar winds and/or supernova shocks originating in the nuclear star cluster have compressed, and likely pushed outward, the infalling molecular gas, thus significantly reducing the gas supply to the central 10pc. Although our analysis currently lacks kinematic confirmation due to the face-on orientation of IC342, the described scenario is supported by the generally observed repetitive nature of star formation in the nuclear star clusters of late-type spiral galaxies.

We have carried out a multiwavelength observational campaign demonstrating some of the remarkable properties of the infrared-bright variable star BP Psc. Surrounded by a compact dusty, gaseous disk, this little-studied late-G (or early-K) type star emits about 75% of its detected energy flux at infrared wavelengths. Evidence for accretion of gas in conjunction with narrow bi-polar jets and Herbig-Haro objects is apparently consistent with classification of BP Psc as a pre-main sequence star, as postulated in most previous studies. If young, then BP Psc would be one of the nearest and oldest known classical T Tauri stars. However, such an evolutionary classification encounters various problems that are absent or much less severe if BP Psc is instead a luminosity class III post-main sequence star. In this case, it would be the first known example of a first ascent giant surrounded by a massive molecular disk with accompanying rapid gas accretion and prominent jets and HH objects. In this model, the genesis of the massive dusty gaseous disk could be a consequence of the envelopment of a low mass companion star. Properties in the disk may be conducive to the current formation of planets, a gigayear or more after the formation of BP Psc itself.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

(Abridged) We present resolution maps of CO, its isotopologues, and HCN from in the center of Maffei 2. The J=1-0 rotational lines of 12^CO, 13^CO, C18^O and HCN, and the J=2-1 lines of 13^CO and C18^O were observed with the OVRO and BIMA arrays. The 2-1/1-0 line ratios of the isotopologues constrain the bulk of the molecular gas to originate in low excitation, subthermal gas. From LVG modeling, we infer that the central GMCs have n(H_2) ~10^2.75 cm^-3 and T_k ~ 30 K. Continuum emission at 3.4 mm, 2.7 mm and 1.4 mm was mapped to determine the distribution and amount of HII regions and dust. Column densities derived from C18^O and 1.4 mm dust continuum fluxes indicate the CO conversion factor in the center of Maffei 2 is lower than Galactic by factors of ~2-4. Gas morphology and the clear ``parallelogram'' in the Position-Velocity diagram shows that molecular gas orbits within the potential of a nuclear (~220 pc) bar. The nuclear bar is distinct from the bar that governs the large scale morphology of Maffei 2. Giant molecular clouds in the nucleus are nonspherical and have large linewidths. Dense gas and star formation are concentrated at the sites of the x_1-x_2 orbit intersections of the nuclear bar, suggesting that the starburst is dynamically triggered.

ALMA will be able to detect a broad spectrum of molecular lines in galaxies. Current aperture synthesis observations indicate that the molecular line emission from galaxies is remarkably variable, even on kpc scales. Imaging spectroscopy at resolutions of an arcsecond or better will reduce the chemical complexity by allowing regions of physical conditions to be defined and classified.

We have measured the mid-infrared radiation from an orientation-unbiased sample of powerful 3C RR galaxies and quasars using the IRS and MIPS instruments aboard the Spitzer Space Telescope. We fit the Spitzer data as well as other measurements from the literature with synchrotron and dust components. At 15 microns, quasars are typically four times brighter than radio galaxies with the same isotropic radio power. Based on our fits, half of this difference can be attributed to the presence of non-thermal emission in the quasars but not the radio galaxies. The other half is consistent with dust absorption in the radio galaxies but not the quasars.

We have measured mid-infrared radiation from an orientation-unbiased sample of 3CRR galaxies and quasars at redshifts 0.4 < z < 1.2 with the IRS and MIPS instruments on the Spitzer Space Telescope. Powerful emission (L_24micron > 10^22.4 W/Hz/sr) was detected from all but one of the sources. We fit the Spitzer data as well as other measurements from the literature with synchrotron and dust components. The IRS data provide powerful constraints on the fits. At 15 microns, quasars are typically four times brighter than radio galaxies with the same isotropic radio power. Based on our fits, half of this difference can be attributed to the presence of non-thermal emission in the quasars but not the radio galaxies. The other half is consistent with dust absorption in the radio galaxies but not the quasars. Fitted optical depths are anti-correlated with core dominance, from which we infer an equatorial distribution of dust around the central engine. The median optical depth at 9.7 microns for objects with core-dominance factor R > 10^-2 is approximately 0.4; for objects with R < 10^-2, it is 1.1. We have thus addressed a long-standing question in the unification of FR II quasars and galaxies: quasars are more luminous in the mid-infrared than galaxies because of a combination of Doppler-boosted synchrotron emission in quasars and extinction in galaxies, both orientation-dependent effects.

We present high spectral resolution (v~12-16 km/s) Brackett line spectroscopy of the blue compact dwarf galaxy Henize 2-10 made with NIRSPEC on the Keck Telescope. The spatial resolution is seeing limited, at 1". We detect two distinct kinematic features separated by approximately 3", with heliocentric velocities of ~860 and ~890 km/s. In addition to a narrow core, the line profiles also display a broad, low intensity feature on the blue side of the centroid, which we attribute to an outflow. This may be a sign of aging in the clusters. We compare to archival high resolution Very Large Array (VLA) data at 1.3 cm, and find that the centimeter wavelength emission is resolved into six sources. These radio sources are organized into two larger groups, which we associate with the two kinematic peaks in the Brackett spectrum. We estimate a Lyman continuum rate of at least 7 x 10^52 s^-1, with a corresponding stellar mass of 6 x 10^6 M_sun is required to ionize the nebulae. We also estimate the size of the nebulae from the radio continuum brightness and find that the observed sources probably contain many HII regions in smaller, unresolved clumps. Brackett line profiles have supersonic line widths, but, aside from the blue wing, are comparable to line widths observed in Galactic ultracompact HII regions, which are excited by a single star, or a few stars.

We report subarcsecond-resolution VLA imaging of the centers of four nearby spiral galaxies: IC 342, Maffei II, NGC 2903, and NGC 6946. In each galaxy, 7 - 12 compact radio continuum sources were identified within the central 15" x 15". Slightly over half of the compact sources appear to be HII regions with flat or positive spectral indices (alpha >~ -0.1). The HII regions with rising spectra are optically thick at centimeter wavelengths, and thus dense (n_i ~ 10^4 cm^(-3)) and young. The largest of these HII regions require the excitation of 500 - 800 O stars, within regions of only few parsecs extent. These clusters approach the sizes expected for globular clusters.

We report subarcsecond-resolution VLA imaging of four nearby spiral galaxies: IC 342, Maffei II, NGC 2903, and NGC 6946. In each galaxy, compact radio continuum sources are identified in the central ~ 15" x 15" region. These compact sources are responsible for 20 - 30 % of the total emission from the central kpc of the host galaxies at 2 cm, but only ~ 5 - 10 % at 6 cm. More than half of the compact sources appear to be HII regions. The HII regions with rising spectra must be fairly dense (n_i ~ 10^4 cm^(-3)) and are presumably very young. The largest of these HII regions require the excitation of 500 - 800 O stars, within regions of only few parsecs extent. These clusters approach the sizes expected for globular clusters. Thermal free-free emission from compact sources contributes more significantly at 2 cm, while diffuse synchrotron emission dominates at 6 cm. The radio HII regions are found near the centers of giant molecular clouds in projection, and do not have obvious visual counterparts.

We present the results of a VLBI Space Observatory Programme (VSOP) observation of the subparsec structure in Centaurus A at 4.9 GHz. Owing to its proximity, our Centaurus A space-VLBI image is one of the highest spatial resolution images of an AGN ever made -- 0.01 pc per beam. The elongated core region is resolved into several components over 10 milli-arcseconds long (0.2 pc) including a compact component of brightness temperature 2.2x10^10K. We analyze the jet geometry in terms of collimation. Assuming the strongest component to be the core, the jet opening angle at ~ 5,000 r_s (Schwarzchild radii) from the core is estimated to be ~ 12 degree, with collimation of the jet to ~ 3 degree continuing out to ~ 20,000 r_s. This result is consistent with previous studies of the jet in M87, which favor MHD disk outflow models. Future space VLBI observations at higher frequencies will probably be able to image the collimation region, within 1,000 r_s of the center of Centaurus A, together with the accretion disk itself.

In this paper I propose that the inner part of a black hole accretion inflow ($< 100 r_{g}$) may enter a magnetically-dominated, magnetosphere-like phase in which the strong, well-ordered fields play a more important role than weak, turbulent fields. In the low/hard state this flow is interior to the standard ADAF usually invoked to explain the observed hot, optically thin emission. Preliminary solutions for these new MDAFs are presented. Time-dependent X-ray and radio observations give considerable insight into these processes, and a new interpretation of the X-ray power spectrum (as arising from many disk radii) may be in order. While an evaporative ADAF model explains the noise power above 0.01 Hz, an inner MDAF is needed to explain the high frequency cutoff near 1 Hz, the presence of a QPO, and the production of a jet. The MDAF scenario also is consistent with the data-based, phemonenological models presented at this meeting by several authors.

We have imaged emission from the millimeter lines of eight molecules--C2H, C34S, N2H+, CH3OH, HNCO, HNC, HC3N, and SO--in the central half kpc of the nearby spiral galaxy IC 342. The 5" (~50 pc) resolution images were made with OVRO. Using these maps we obtain a picture of the chemistry within the nuclear region on the sizescales of individual GMCs. Bright emission is detected from all but SO. There are marked differences in morphology for the different molecules. A principal component analysis is performed to quantify similarities and differences among the images. This analysis reveals that while all molecules are to zeroth order correlated, that is, they are all found in dense molecular clouds, there are three distinct groups of molecules distinguished by the location of their emission within the nuclear region. N2H+, C18O, HNC and HCN are widespread and bright, good overall tracers of dense molecular gas. C2H and C34S, tracers of PDR chemistry, originate exclusively from the central 50-100 pc region, where radiation fields are high. The third group of molecules, CH3OH and HNCO, correlates well with the expected locations of bar-induced orbital shocks. The good correlation of HNCO with the established shock tracer molecule CH3OH is evidence that this molecule, whose chemistry has been uncertain, is indeed produced by processing of grains. HC3N is observed to correlate tightly with 3mm continuum emission, demonstrating that the young starbursts are the sites of the warmest and densest molecular gas. We compare our HNC images with the HCN images of Downes et al. (1992) to produce the first high resolution, extragalactic HCN/HNC map: the HNC/HCN ratio is near unity across the nucleus and the correlation of both of these gas tracers with the star formation is excellent. (Abridged).

Non-relativistic three-dimensional magnetohydrodynamical (MHD) simulations of Poynting flux dominated (PFD) jets are presented. Our study focuses on the propagation of strongly magnetized hypersonic, but sub-Alfv\'enic ($C^{2}_{\rm s} \ll V^{2}_{\rm jet} < V^{2}_{\rm A}$) flow and on the subsequent development of a current-driven (CD) kink instability. This instability may be responsible for the ``wiggled'' structures seen in sub-parsec scale (VLBI) jets. In the present paper, we investigate the nonlinear behavior of PFD jets in a variety of external ambient magnetized gas distributions, including those with density, pressure, and temperature gradients. Our numerical results show that the jets can develop CD distortions in the trans-Alfv\'enic flow case, even when the flow itself is still strongly magnetically dominated. An internal non-axisymmetric body mode grows on time scales of order of the Alfv\'en crossing time and distorts the structure and magnetic configuration of the jet. The kink ($m=1$) mode of the CD instability, driven by the radial component of the Lorentz force, grows faster than other higher order modes ($m>1$). In the jet frame the mode grows locally and expands radially at each axial position where the jet is unstable: the instability, therefore, does not propagate as a wave along the jet length. A naturally-occurring, external helically magnetized wind, which is (quasi-) axially current-free, surrounds the well-collimated current-carrying jet and reduces velocity shear between the jet and external medium. This stabilizes the growth of MHD Kelvin-Helmholtz surface modes in the inner jet flow.

We present high resolution (~5'') maps of the J = 1 - 0 transitions of ^{13}CO and C^{18}O towards the nucleus of NGC 6946, made with the Owens Valley Millimeter Array. The images are compared with existing ^{12}CO(1-0) maps to investigate localized changes in gas properties across the nucleus. As compared to ^{12}CO, both ^{13}CO and C^{18}O are more confined to the central ring of molecular gas associated with the nuclear star formation; that is, ^{12}CO is stronger relative to ^{13}CO and C^{18}O away from the nucleus and along the spiral arms. The ^{12}CO(1-0)/^{13}CO(1-0) line ratio reaches very high values of >40. We attribute the relative ^{13}CO weakness to a rapid change in the interstellar medium from dense star forming cores in a central ring to diffuse, low density molecular gas in and behind the molecular arms. This change is abrupt, occurring in less than a beamsize (90 pc), about the size of a giant molecular cloud. Column densities determined from ^{13}CO(1-0), C^{18}O(1-0), and 1.4 mm dust continuum all indicate that the standard Galactic conversion factor, X_{CO}, overestimates the amount of molecular gas in NGC 6946 by factors of ~3-5 towards the central ring and potentially even more so in the diffuse gas away from the central starburst. We suggest that the nuclear bar acts to create coherent regions of molecular clouds with distinct and different physical conditions. The ^{12}CO(1-0)/^{13}CO(1-0) line ratio in galactic nuclei can be a signpost of a dynamically evolving ISM.

I review recent progress in the theory of relativistic jet production. The presently favored mechanism is an electrodynamic one, in which charged plasma is accelerated by electric fields that are generated by a rotating magnetic field. The most pressing issues of current interest are understanding what factors control the jet power, its speed, and its degree of collimation, and how these properties determine the type of jet observed and its effect on its environment. Recent observations of microquasars, pulsars, gamma-ray bursts (GRBs) and core-collapse supernovae indicate that jets play an important, and in some cases possibly dominant, role in each of these. The presence of a jet in all cases may provide an important clue to how these sources may be related. Based on these observations, and on recent theoretical investigations, I propose an evolutionary scheme that attempts to unify all of these relativistic galactic jet sources. I also discuss several important issues that must be resolved before this (or another scheme) can be adopted.

We present the basic elements and first results of an end-to-end simulation package whose purpose is to test the validity of the Space Interferometer Mission design. The fundamental simulation time step is one millisecond, with substructure at 1/8 ms, and the total duration of the simulation is five years. The end product of a given "wide-angle" astrometry run is an estimated grid star catalog over the entire sky with an accuracy of about 4 micro-arcseconds. SIMsim is divided into five separate modules that communicate via data pipes. The first generates the 'truth' data on the spacecraft structure and laser metrology. The second module generates uncorrupted fringes for the three SIM interferometers, based on the current spacecraft orientation, target stars' positions, etc. The third module reads out the CCD white light fringe data at specified times, corrupting that and the metrology data with appropriate errors. The data stream out of this module represents the basic data stream on the simulated spacecraft. The fourth module performs fringe-fitting tasks on this data, recovering the total path delay, and the fifth and final module inverts the entire metrology/delay dataset to ultimately determine the instantaneous path delay on a fiducial baseline fixed in space. (Pathlength feed forward is used every few milliseconds to re-position the interferometer to keep the fringes in the delay window.) The average of all such delays over an integration time (typically 30 s) is reported as one of several hundred thousand measured stellar delays over the five-year period, which are then inverted to produce the simulated catalog. Future plans include taking into account more sources of error from the SIM error budget and including narrow angle observations in the observing plan.

I review recent progress in the theory of relativistic jet production, with special emphasis on unifying black hole sources of stellar and supermassive size. Observations of both classes of objects, as well as theoretical considerations, indicate that such jets may be launched with a spine/sheath flow structure, having a much higher Lorentz factor ($\sim 50$) near the axis and a lower speed ($\Gamma \sim 10$ or so) away from the axis. It has become clear that one can no longer consider models of accretion flows without also considering the production of a jet by that flow. Furthermore, the rotation rate of the black hole also must be taken into account. It provides a third parameter that should break the mass/accretion rate degeneracy and perhaps explain why some sources are radio loud and some radio quiet. Slow jet acceleration and collimation is expected theoretically, and can explain some of the observed features of AGN jet sources. Finally, relativistic jets launched by MHD/ED processes are Poynting flux dominated by nature, and are potentially unstable if there is significant entrainment of thermal material.

Non-relativistic three-dimensional magnetohydrodynamic simulations of Poynting-flux-dominated (PFD) jets are presented. Our study focuses on the propagation of strongly magnetized hypersonic but sub-Alfv\'enic flow ($C_{\rm s}^2 << V_{\rm jet}^2 < V_{\rm A}^2$) and the development of a current-driven (CD) kink instability. This instability may be responsible for the "wiggled" structures seen in VLBI-scale AGN jets. In the present paper we investigate the nonlinear behavior of PFD jets in a variety of external ambient magnetized gas distributions, including those with density, pressure, and temperature gradients. Our numerical results show that PFD jets can develop kink distortions in the trans-Alfv\'enic flow case, even when the flow itself is still strongly magnetically dominated. In the nonlinear development of the instability, a non-axisymmetric mode grows on time scales of order the Alfv\'en crossing time (in the jet frame) and proceeds to disrupt the kinematic and magnetic structure of the jet. Because of a large scale poloidal magnetic field in the ambient medium, the growth of surface modes ({\it i.e.}, MHD Kelvin-Helmholtz instabilities) is suppressed. The CD kink mode ($m = 1$) grows faster than the other higher order modes ($m > 1$), driven in large part by the radial component of the Lorentz force.

Fully relativistic and causal equations for the flow of charge in curved spacetime are derived. It is believed that this is the first set of equations to be published that correctly describes the flow of charge, and evolution of the electromagnetic field, in highly dynamical relativistic environments on time scales much shorter than the collapse time ($GM/c^3$). The equations will be important for correctly investigating problems such as the dynamical collapse of magnetized stellar cores to black holes and the production of jets and gravitational waves. This system of equations, given the name of `charge dynamics', is analogous to those of hydrodynamics (which describe the flow of {\em mass} in spacetime rather than the flow of charge). The most important one in the system is the relativistic generalized Ohm's law, which is used to compute time-dependent four-current. Unlike previous equations for the current, this one ensures that charge drift velocities remain less than the speed of light, takes into account the finite current rise time, is expressed in a covariant form, and is suitable for general relativistic computations in an arbitrary metric. It includes the standard known effects (Lorentz force, Hall effect, pressure effect, and resistivity) and reduces to known forms of Ohm's law in the appropriate limits. In addition, the plasma particles are allowed to have highly relativistic drift velocities, resulting in an implicit equation for the `current beaming factor' $\gamma_q$.

A scheme is presented for accurately propagating the gravitational field constraints in finite difference implementations of numerical relativity. The method is based on similar techniques used in astrophysical magnetohydrodynamics and engineering electromagnetics, and has properties of a finite differential calculus on a four-dimensional manifold. It is motivated by the arguments that 1) an evolutionary scheme that naturally satisfies the Bianchi identities will propagate the constraints, and 2) methods in which temporal and spatial derivatives commute will satisfy the Bianchi identities implicitly. The proposed algorithm exactly propagates the constraints in a local Riemann normal coordinate system; {\it i.e.}, all terms in the Bianchi identities (which all vary as $\partial^3 g$) cancel to machine roundoff accuracy at each time step. In a general coordinate basis, these terms, and those that vary as $\partial g\partial^2 g$, also can be made to cancel, but differences of connection terms, proportional to $(\partial g)^3$, will remain, resulting in a net truncation error. Detailed and complex numerical experiments with four-dimensional staggered grids will be needed to completely examine the stability and convergence properties of this method. If such techniques are successful for finite difference implementations of numerical relativity, other implementations, such as finite element (and eventually pseudo-spectral) techniques, might benefit from schemes that use four-dimensional grids and that have temporal and spatial derivatives that commute.

We review recent progress in the theory of jet production, with particular emphasis on the possibility of 1) powerful jets being produced in the first few seconds after collapse of a supernova core and 2) those jets being responsible for the asymmetric explosion itself. The presently favored jet-production mechanism is an electrodynamic one, in which charged plasma is accelerated by electric fields that are generated by a rotating magnetic field anchored in the protopulsar. Recent observations of Galactic jet sources provide important clues to how {\em all} such sources may be related, both in the physical mechanism that drives the jet in the astrophysical mechanisms that create conditions conducive to jet formation. We propose a grand evolutionary scheme that attempts to unify these sources on this basis, with MHD supernovae providing the missing link. We also discuss several important issues that must be resolved before this (or another scheme) can be adopted.

Little is known about the origins of the giant star clusters known as globular clusters. How can hundreds of thousands of stars form simultaneously in a volume only a few light years across the distance of the sun to its nearest neighbor? Radiation pressure and winds from luminous young stars should disperse the star-forming gas and disrupt the formation of the cluster. Globular clusters in our Galaxy cannot provide answers; they are billions of years old. Here we report the measurement of infrared hydrogen recombination lines from a young, forming super star cluster in the dwarf galaxy, NGC 5253. The lines arise in gas heated by a cluster of an estimated million stars, so young that it is still enshrouded in gas and dust, hidden from optical view. We verify that the cluster contains 4000-6000 massive, hot "O" stars. Our discovery that the gases within the cluster are bound by gravity may explain why these windy and luminous O stars have not yet blown away the gases to allow the cluster to emerge from its birth cocoon. Young clusters in "starbursting" galaxies in the local and distant universe may be similarly gravitationally confined and cloaked from view.

We present high spatial resolution (1.2") mm-interferometric observations of the CO(2-1) line emission in the central 300pc of the late-type spiral galaxy IC342. The data, obtained with the Owens Valley Radio Observatory, allow first-time detection of a molecular gas disk that coincides with the luminous young stellar cluster in the nucleus of IC342. The nuclear CO disk has a diameter of ~30pc and a molecular gas mass of M(H2) ~ 2x10^5 M_sun. It connects via two faint CO bridges to the well-known, 100pc diameter circumnuclear gas ring. Analysis of the gas kinematics shows that the line-of-nodes in the inner 15pc is offset by about 45degree from the major kinematic axis, indicating non-circular motion of the gas within a few parsec of the dynamical center of IC342. Both the morphology and the kinematics of the CO gas indicate ongoing inflow of molecular gas into the central few parsec of IC342. We infer a gas inflow rate between 0.003 and 0.14 M_sun/yr, based on the observed surface density of the nuclear gas disk and estimates of the radial velocities of the surrounding gas. Inflow rates of this order can support repetitive star formation events in the nucleus of IC342 on timescales much smaller than a Hubble time.

We investigate the action of the magnetorotational instability (MRI) in the context of iron-core collapse. Exponential growth of the field on the rotation time scale by the MRI will dominate the linear growth process of field line "wrapping" with the same characteristic time. We examine a variety of initial rotation states, with solid body rotation or a gradient in rotational velocity, that correspond to models in the literature. A relatively modest value of the initial rotation, a period of ~ 10 s, will give a very rapidly rotating PNS and hence strong differential rotation with respect to the infalling matter. We assume conservation of angular momentum on spherical shells. Results are discussed for two examples of saturation fields, a fiducial field that corresponds to Alfven velocity = rotational velocity and a field that corresponds to the maximum growing mode of the MRI. Modest initial rotation velocities of the iron core result in sub-Keplerian rotation and a sub-equipartition magnetic field that nevertheless produce substantial MHD luminosity and hoop stresses: saturation fields of order 10^{15} - 10^{16} G develop within 300 msec after bounce with an associated MHD luminosity of about 10^{52} erg/s. Bi-polar flows driven by this MHD power can affect or even cause the explosions associated with core-collapse supernovae.

We report subarcsecond-resolution VLA and Keck mid-infrared imaging of the dwarf starburst galaxy II Zw 40. II Zw 40 contains a bright, compact thermal radio and infrared source with all the characteristics of a collection of dense HII regions ionized by at least 14,000 O stars. The supernebula is revealed to be multiple sources within an envelope of weaker emission. The radio emission is dominated by free-free emission at 2cm, and the spectrum of this emission appears to be rising. This suggests that the free-free emission is optically thick at 2cm, and that the individual HII regions are ~1pc in size. This complex of "supernebulae" dominates the total infrared luminosity of II Zw 40, although the radio source is less than ~150pc in diameter. Multiple super star clusters appear to be forming here, the much larger analogues of large Galactic HII region complexes.

A small sample of 7 southern Compact Steep Spectrum (CSS) radio sources has been selected as part of the study of a larger flux-limited complete sample of radio sources. High resolution images, using the VLBI network in the southern hemisphere and the high resolution MERLIN array, are presented for all sources in the CSS sample. The overall morphology of each source consists of well-defined double lobes but with substantial diffuse and extended components present. In the majority of cases only a fraction of the total flux density is detected on the VLBI baselines, indicating the presence of larger extended radio structures. However, all sources are unresolved at arcsecond scales and are of sub-galactic size, with linear size in the range 0.1-2 kpc. The radio properties of the sources agree well with CSS sources in other samples.

We report the detection of CO(2-1) and 3.1 mm and 1.3 mm continuum emission towards the extremely young starburst in NGC 5253, with data taken from the Owens Valley Millimeter Array. Faint CO emission originates in five molecular clouds distributed along the prominent dust lane seen in visual images. With the gas, the morphology of NGC 5253 looks much like a dwarf elliptical version of the ``dust-lane ellipticals'' or ``polar-ring'' class of galaxies. The molecular gas resides in GMCs well away from the radio-IR super-star cluster/supernebula seen in the radio and infrared. The millimeter continuum data confirm that the 2 cm flux from the supernebula is optically thick; the Lyman continuum rate derived from the 1.3 mm continuum is N$_{Lyc} \sim 6x10^{52} s^{-1}$ for the central ~20''. CO may underestimate the true molecular column density, as expected for a low metallicity system, although there are regions along the dust lane that appear to have near-Galactic conversion factors. We estimate a total molecular gas mass of $M_{H_{2}}\lsim 10^{7} M_{\odot}$. The molecular gas in the dust lane is falling into the galaxy, supporting an accretion hypothesis. The dust lane gas cannot therefore be causally associated with the current burst of star formation. A relatively small amount, $M_{H_{2}}\lsim 5x10^{5} M_{\odot}$, of molecular gas is associated with the current starburst. We estimate a star formation efficiency of at least 25 % and more likely ~75 %, consistent with the formation of a bound cluster. Despite the extreme youth of the starburst, the specific trigger of the starburst remains elusive, although the infall of gas in the dust lane suggests that there is more star formation to come in NGC 5253.

Strong toroidal magnetic fields generated in stellar collapse can generate magneto-centrifugal jets in analogy to those found in simulations of black hole accretion and explain why all core collapse supernovae are found to be substantially asymmetric and predominantly bi-polar. We describe two phases: the initial LeBlanc-Wilson jet and a subsequent protopulsar or toroidal jet that propagates at about the core escape velocity. The jets will produce bow shocks that tend to expel matter, including iron and silicon, into equatorial tori, accounting for observations of the element distribution in Cas A. A magnetic ``switch'' mechanism may apply in instances of low density and large magnetic field with subsequent increase in the speed and collimation of the toroidal jet, depositing relatively little momentum. The result could be enough infall to form a black hole with a third, highly relativistic jet that could catch up to the protopulsar jet after it has emerged from the star. The interaction of these two jets could generate internal shocks and explain the presence of iron lines in the afterglow. Recent estimates that typical gamma-ray burst energy is about 3x10^50 erg imply either a very low efficiency for conversion of rotation into jets, or a rather rapid turnoff of the jet process even though the black hole still rotates rapidly. Magnetars and ``hypernovae'' might arise in an intermediate parameter regime of energetic jets that yield larger magnetic fields and provide more energy than the routine case, but that are not so tightly collimated that they yield failed supernova. (slightly abridged)

We present high resolution (~2.5'') CO(1-0) and CO(2-1) images of the central kiloparsec of NGC 3077 made with the Owens Valley Millimeter Array. CO emission is distributed in three major complexes which resolve into at least seven GMCs. Two complexes are associated with the central starburst. A third, more distant complex is not associated with strong star formation. The GMCs are ~70 pc in size and contain ~10^{6} M_{\odot} of molecular gas. The Galactic conversion factor appears applicable to NGC 3077, consistent with its solar metallicity. Galactic rotation in NGC 3077 is detected for the first time in the molecular gas. The molecular gas counterrotates with respect to the large scale HI tidal bridge. The molecular clouds closest to the starburst have super-virial linewidths, possibly related to the turbulence generated by the starburst. The 2.6 mm radio continuum flux indicates that thermal Bremsstrahlung dominates the emission from the starburst region below 5 cm, and that the N_{Lyc} \simeq 3.7\times 10^{52} s^{-1} or ~3000 O7 stars, corresponds to a star formation rate of 0.4 M_{\odot} yr^{-1}. At this rate the amount of molecular gas can sustain star formation for only ~10 Myrs; thus NGC 3077 is a true starburst galaxy. The derived age of the starburst is consistent with the inferred ages of the superbubbles, which suggests that the burst is much younger than the age of the M 81-M 82-NGC 3077 interaction. We suggest that it is caused by gas that was pulled out of NGC 3077 during the interaction with M 81, raining back down onto the galaxy.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption (by ionized gas in the disk) of background emission from a radio counterjet. We report high angular resolution VLBI observations of the nearby (about 40 Mpc) radio galaxy NGC 4261 that confirm free-free absorption of radio emission from a counterjet by a geometrically thin, nearly edge-on disk at 1.6, 4.8, and 8.4 GHz. The angular width and depth of the absorption appears to increase with decreasing frequency, as expected. We derive an average electron density of ~10E4 per cc at a disk radius of about 0.2 pc, assuming that the inner disk inclination and opening angles are the same as at larger radii. Pressure balance between the thermal gas and the magnetic field in the disk implies an average field strength of 0.1 milligauss at a radius of 0.2 pc. These are the closest-in free-free absorption measurements to date of the conditions in an extragalactic accretion disk orbiting a black hole with a well-determined mass. If a standard advection-dominated accretion flow exists in the disk center, then the transition between thin and thick disk regions must occur at a radius less than 0.2 pc (4000 Schwarzschild radii).

We present high resolution maps ($\sim 2^{''}$) maps of the J=1-0 and J=2-1 transitions of C$^{18}$O in the central $\sim 150$ pc of the gas-rich nucleus of IC 342, made with the Owens Valley Millimeter Array. From the C$^{18}$O maps we are able to obtain the most accurate map of N$_{H_{2}}$ to date for IC 342. Due to their low opacities, the transitions of C$^{18}$O give a more reliable estimate of the true molecular gas column density distribution than the more common $^{12}$CO and $^{13}$CO isotopomers. The morphology of the C$^{18}$O emission in the nucleus is a mini-spiral similar to that of the main isotopomer, $^{12}$CO, except it is more symmetric, lacking the enhancements to north. We suggest that the asymmetries present in $^{12}$CO images may reflect the viewing perspective of the starburst region biased by the high optical depths of $^{12}$CO, rather than true asymmetries in the amount of molecular gas present. The giant molecular clouds seen in C$^{18}$O appear to be nonspherical, probably due to tidal arm shearing. Column densities determined from C$^{18}$O observations, 1.3 mm dust continuum, and the virial theorem indicate that the standard Galactic conversion factor, X$_{CO}$, overestimates the amount of molecular gas in the center of IC 342 by a factor of $\sim$2-3 at the molecular cloud peaks, and by $>$3 in the diffuse gas away from the starburst. Revised molecular masses based on this conversion factor imply that star formation efficiencies in the starburst region are very high. From the distribution of gas and star formation it appears that the sites of star formation are dynamically determined, rather than driven by density peaks. Near the central star-forming region, evidence is seen for chemical enrichment of C$^{18}$O due to massive stars.

Eight dwarf starburst galaxies have been observed with the Caltech Submillimeter Observatory (CSO) telescope in the CO J= 3 - 2 transition. The galaxies observed are He 2-10, NGC 5253, NGC 1569, NGC 3077, Haro 2, Haro 3, II Zw 40 and Mrk 86; all but the last two are detected. The central regions of He 2-10 and NGC 5253 were mapped and a CO(2-1) spectrum of NGC 5253 was obtained. The error weighted mean CO(3-2)/CO(1-0) ratio of the detected galaxies is 0.60$\pm$0.06, which is virtually identical to what is found for starbursts in the nuclei of nearby spirals, and suggests that the molecular gas is optically thick, warm (T$_{k}>$20 K), and moderately dense ($n_{H_{2}}\sim 10^{3-4} cm^{-3}$). The CO(3-2)/CO(1-0) ratio peaks at or close to the starburst in all cases. CO emission does not appear to be optically thin in these dwarfs, despite the low metallicity and intense radiation fields, which is probably because in order for CO to exist in detectable amounts it must be self-shielding and hence optically thick. Physical properties of the molecular clouds in these dwarf starbursts appear to be essentially the same as nearby spiral nuclei, with the possible exception that CO is more confined to the cloud cores.

A model is presented in which the strongest radio-emitting jet outflows are produced in black hole systems when the accretion is a geometrically thick ($H/R \sim 1$) inflow ({\it e.g.}, ADAF, CDAF) {\em and} if the black hole is rotating. For galactic black hole candidates, the model naturally accounts for the observed correlation of jet outflow with the black hole hard emission state and predicts an association of strong jets with rapid black hole rotation. When extended to the supermassive case, the model accounts for the highest radio galaxy and quasar jet powers and provides additional theoretical support for the ``spin paradigm'', which asserts that radio loud quasars are produced by Kerr holes and radio quiet ones by Schwarzschild holes. In some cases, the angular momentum and energy outflow in the jet may be large enough to significantly alter the structure of the accretion flow from that predicted by current models.

The physical conditions in the inner parsec of accretion disks believed to orbit the central black holes in active galactic nuclei can be probed by imaging the absorption of background radio emission by ionized gas in the disk. We report high angular resolution observations of the nearby galaxy NGC 4261 which show evidence for free-free absorption by a thin, nearly edge-on disk at several frequencies. The angular width, and probably the depth, of the absorption appears to increase with decreasing frequency, as expected. Because free-free absorption is much larger at lower frequencies, the longest possible baselines are needed to provide adequate angular resolution; observing at higher frequencies to improve resolution will not help.

The structure of AGN accretion disks on sub-parsec scales can be probed through free-free absorption of synchrotron emission from the base of symmetric radio jets. We present new VLBA observations of the nearby FR-I radio galaxy NGC 4261 at 22 and 43 GHz, and combine these with previous VLBA observations at 1.6 and 8.4 GHz to map absorption caused by an inner accretion disk. Assuming the disk is geometrically and optically thin and composed of a uniform 10^4 K plasma, the average electron density in the inner 0.1 pc is 10^3 - 10^8 cm^-3. Equating thermal gas pressure and magnetic field strength gives a disk magnetic field of 10^-4 - 10^-2 Gauss at 0.1 pc. The jet opening angle is between 0.3 and 20 degrees during the first 0.2 pc of the jet, and must be less than 5 degrees during the first 0.8 pc. We include an appendix containing expressions for a simple, optically thin, gas pressure dominated accretion disk model which may be applicable to other galaxies in addition to NGC 4261.

While attempts to unify certain classes of AGN using orientation and environmental effects have been successful, it is widely recognized that intrinsic properties of the accreting black hole system also must play a role in determining the appearance of such an object. In addition to mass and accretion rate, the angular momentum (or spin) of the black hole can play a crucial role in determining the power of a relativistic jet that is generated by magnetohydrodynamic acceleration near the hole. In this paper a scenario is presented, based on accretion theory and recent models of MHD jet production, in which the primary (although not only) parameter differentiating between radio loud and quiet objects is the black hole spin, and that determining quasar vs. radio galaxy is the accretion rate. A surprising number of desirable features result from these simple concepts and the accompanying equations. In addition, there are several testable predictions that can determine whether this grand unification scheme has further merit.

A model is presented in which both Fanaroff and Riley class I and II extragalactic jets are produced by magnetized accretion disk coronae in the ergospheres of rotating black holes. While the jets are produced in the accretion disk itself, the output power still is an increasing function of the black hole angular momentum. For high enough spin, the black hole triggers the magnetic switch, producing highly-relativistic, kinetic-energy-dominated jets instead of Poynting-flux-dominated ones for lower spin. The coronal mass densities needed to trigger the switch at the observed FR break power are quite small ($\sim 10^{-15} g cm^{-3}$), implying that the source of the jet material may be either a pair plasma or very tenuous electron-proton corona, not the main accretion disk itself. The model explains the differences in morphology and Mach number between FR I and II sources and the observed trend for massive galaxies to undergo the FR I/II transition at higher radio power. It also is consistent with the energy content of extended radio lobes and explains why, because of black hole spindown, the space density of FR II sources should evolve more rapidly than that of FR I sources. If the present model is correct, then the ensemble average speed of parsec-scale jets in sources distinguished by their FR I morphology (not luminosity) should be distinctly slower than that for sources with FR II morphology. The model also suggests the existence of a population of high-redshift, sub-mJy FR I and II radio sources associated with spiral or pre-spiral galaxies that flared once when their black holes were formed but were never again re-kindled by mergers.

Recent observations of Galactic Black Hole Candidates (BHCs) suggest that those that are superluminal jet sources have more rapid black hole spin rates than otherwise normal BHCs. This provides observational support for models of relativistic jet formation that extract rotational energy of the central black hole. To investigate this mechanism, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of the formation of a relativistic jet in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We assume that the initial velocity of the disk is zero. We find that the maximum velocity of the jet reaches 0.93c (Lorentz factor, 2.7) and the terminal velocity of the jet is 0.85c (Lorentz factor, 1.9). On the other hand, for a non-rotating (a=0) Schwarzschild black hole, the maximum outflow velocity is less than 0.6c for initial magnetospheric conditions similar to to those of the Kerr black hole case. These numerical results show the importance of the rapidly rotating black hole for the relativistic jet formation.

To investigate the formation mechanism of relativistic jets in active galactic nuclei and micro-quasars, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of jet formation in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We study cases in which the Keplerian accretion disk is both co-rotating and counter-rotating with respect to the black hole rotation. In the co-rotating disk case, our results are almost the same as those in Schwarzschild black hole cases: a gas pressure-driven jet is formed by a shock in the disk, and a weaker magnetically-driven jet is also generated outside the gas pressure-driven jet. On the other hand, in the counter-rotating disk case, a new powerful magnetically-driven jet is formed inside the gas pressure-driven jet. The newly found magnetically-driven jet in the latter case is accelerated by a strong magnetic field created by frame dragging in the ergosphere. Through this process, the magnetic field extracts the energy of the black hole rotation.

A new field of numerical astrophysics is introduced which addresses the solution of large, multidimensional structural or slowly-evolving problems (rotating stars, interacting binaries, thick advective accretion disks, four dimensional spacetimes, etc.). The technique employed is the Finite Element Method (FEM), commonly used to solve engineering structural problems. The approach developed herein has the following key features: 1. The computational mesh can extend into the time dimension, as well as space, perhaps only a few cells, or throughout spacetime. 2. Virtually all equations describing the astrophysics of continuous media, including the field equations, can be written in a compact form similar to that routinely solved by most engineering finite element codes. 3. The transformations that occur naturally in the four-dimensional FEM possess both coordinate and boost features, such that (a) although the computational mesh may have a complex, non-analytic, curvilinear structure, the physical equations still can be written in a simple coordinate system independent of the mesh geometry. (b) if the mesh has a complex flow velocity with respect to coordinate space, the transformations will form the proper arbitrary Lagrangian- Eulerian advective derivatives automatically. 4. The complex difference equations on the arbitrary curvilinear grid are generated automatically from encoded differential equations. This first paper concentrates on developing a robust and widely-applicable set of techniques using the nonlinear FEM and presents some examples.

We report the results of a 5-GHz southern-hemisphere snapshot VLBI observation of a sample of blazars. The observations were performed with the Southern Hemisphere VLBI Network plus the Shanghai station in 1993 May. Twenty-three flat-spectrum, radio-loud sources were imaged. These are the first VLBI images for 15 of the sources. Eight of the sources are EGRET (> 100 MeV) gamma-ray sources. The milliarcsecond morphology shows a core-jet structure for 12 sources, and a single compact core for the remaining 11. No compact doubles were seen. Compared with other radio images at different epochs and/or different frequencies, 3 core-jet blazars show evidence of bent jets, and there is some evidence for superluminal motion in the cases of 2 blazars. The detailed descriptions for individual blazars are given. This is the second part of a survey: the first part was reported by Shen et al. (AJ 114(1997)1999).

We report the results of a 5 GHz southern hemisphere VLBI survey of compact extragalactic radio sources. These observations were undertaken with the SHEVE array plus Shanghai station in November 1992. A sample of 22 sources was observed and images of 20 of them were obtained. Of the 20 sources imaged, 15 showed core-jet structure, one had a two-sided jet and 4 had only single compact cores. Eleven of the 16 core-jet (including one two-sided jet) sources show some evidence of bent jets. No compact doubles were found. A comparison with previous images and the temporal variability of the radio flux density showed evidence for superluminal motion in 4 of the sources. Five sources were high energy (>100 MeV) gamma-ray sources. Statistical analysis showed the dominance of highly polarized quasars among the detected gamma-ray sources, which emphasizes the importance of beaming effect in the gamma-ray emission.

A sub-parsec scale radio counterjet has been detected in the nucleus of the closest radio galaxy, Centaurus A (NGC 5128), with VLBI imaging at 2.3 and 8.4 GHz. This is one of the first detections of a VLBI counterjet and provides new constraints on the kinematics of the radio jets emerging from the nucleus of Cen A. A bright, compact core is seen at 8.4 GHz, along with a jet extending along P.A. 51 degrees. The core is completely absorbed at 2.3 GHz. Our images show a much wider gap between the base of the main jet and the counterjet at 2.3 GHz than at 8.4 GHz and also that the core has an extraordinarily inverted spectrum. These observations provide evidence that the innermost 0.4-0.8 pc of the source is seen through a disk or torus of ionized gas which is opaque at low frequencies due to free-free absorption.