Work in time-domain astronomy necessitates robust, automated data processing pipelines that operate in real time. We present the BANZAI pipeline which processes the thousands of science images produced across the Las Cumbres Observatory Global Telescope (LCOGT) network of robotic telescopes each night. BANZAI is designed to perform near real-time preview and end-of-night final processing for four types of optical CCD imagers on the three LCOGT telescope classes. It performs instrumental signature removal (bad pixel masking, bias and dark removal, flat-field correction), astrometric fitting and source catalog extraction. We discuss the design considerations for BANZAI, including testing, performance, and extensibility. BANZAI is integrated into the observatory infrastructure and fulfills two critical functions: (1) real-time data processing that delivers data to users quickly and (2) derive metrics from those data products to monitor the health of the telescope network. In the era of time-domain astronomy, to get from these observations to scientific results, we must be able to automatically reduce data with minimal human interaction, but still have insight into the data stream for quality control.

The Las Cumbres Observatory operates a fleet of robotically controlled telescopes currently two 2m, nine 1m, and ten 0.4m telescopes, distributed amongst six sites covering both hemispheres. Telescopes of an aperture class are equipped with an identical set of optical imagers, and those data are subsequently processed by a common pipeline (BANZAI). The telescopes operate without direct human supervision, and assessing the daily and long-term scientific productivity of the fleet of telescopes and instruments poses an operational challenge. One key operational metric of a telescope/instrument system is throughput. We present a method of long-term performance monitoring based on nightly science observations: For every image taken in matching filters and within the footprint of the PANSTARRS DR1 catalog we derive a photometric zeropoint, which is a good proxy for system throughput. This dataset of over $250000$ data points enables us to answer questions about general throughput degradation trends, and how individual telescopes perform at the various sites. This particular metric is useful to plan the effort level for on-site support and to prioritize the cleaning and re-aluminizing schedule of telescope optics and mirrors respectively.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

Phase fluctuations introduced by the atmosphere are the main limiting factor in attaining diffraction limited performance in extended interferometric arrays at millimeter and submillimeter wavelengths. We report the results of C-PACS, the Combined Array for Research in Millimeter-Wave Astronomy Paired Antenna Calibration System. We present a systematic study of several hundred test observations taken during the 2009-2010 winter observing season where we utilize CARMA's eight 3.5-m antennas to monitor an atmospheric calibrator while simultaneously acquiring science observations with 6.1-m and 10.4-m antennas on baselines ranging from a few hundred meters to ~2 km. We find that C-PACS is systematically successful at improving coherence on long baselines under a variety of atmospheric conditions. We find that the angular separation between the atmospheric calibrator and target source is the most important consideration, with consistently successful phase correction at CARMA requiring a suitable calibrator located $\lesssim$6$^\circ$ away from the science target. We show that cloud cover does not affect the success of C-PACS. We demonstrate C-PACS in typical use by applying it to the observations of the nearby very luminous infrared galaxy Arp 193 in $^{12}$CO(2-1) at a linear resolution of ~70 pc (0.12" x 0.18"), 3 times better than previously published molecular maps of this galaxy. We resolve the molecular disk rotation kinematics and the molecular gas distribution and measure the gas surface densities and masses on 90 pc scales. We find that molecular gas constitutes $\sim30\%$ of the dynamical mass in the inner 700 pc of this object with a surface density $\sim10^4 M_\odot$ pc$^{-2}$; we compare these properties to those of the starburst region of NGC 253.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

We present 350 micron polarization observations of four low-mass cores containing Class 0 protostars: L483, L1157, L1448-IRS2, and Serp-FIR1. This is the second paper in a larger survey aimed at testing magnetically regulated models for core-collapse. One key prediction of these models is that the mean magnetic field in a core should be aligned with the symmetry axis (minor axis) of the flattened YSO inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hour-glass shaped morphology as gravity drags the field inward towards the central protostar. We combine our results for the four cores with results for three similar cores that were published in the first paper from our survey. An analysis of the 350 micron polarization data for the seven cores yields evidence of a positive correlation between mean field direction and pseudodisk symmetry axis. Our rough estimate for the probability of obtaining by pure chance a correlation as strong as the one we found is about 5%. In addition, we combine together data for multiple cores to create a source-averaged magnetic field map having improved signal-to-noise ratio, and this map shows good agreement between mean field direction and pseudodisk axis (they are within 15 degrees). We also see hints of a magnetic pinch in the source-averaged map. We conclude that core-scale magnetic fields appear to be strong enough to guide gas infall, as predicted by the magnetically regulated models. Finally, we find evidence of a positive correlation between core magnetic field direction and bipolar outflow axis.

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

The 229 GHz (lambda 1.3mm) radio emission from Orion-KL was mapped with up to 0.14" angular resolution with CARMA, allowing measurements of the flux densities of Source I ('SrcI') and the Becklin-Neugebauer Object (BN), the 2 most massive stars in this region. We find integrated flux densities of 310 +/- 45 mJy for SrcI and 240 +/- 35 mJy for BN. SrcI is optically thick even at 229 GHz. No trace of the H30alpha recombination line is seen in its spectrum, although the v_2=1, 5(5,0)-6(4,3) transition of H2O, 3450 K above the ground state, is prominent. SrcI is elongated at position angle 140 degrees, as in 43 GHz images. These results are most easily reconciled with models in which the radio emission from SrcI arises via the H- free-free opacity in a T < 4500 K disk, as considered by Reid et al. (2007). By contrast, the radio spectrum of BN is consistent with p+/e- free-free emission from a dense (n_e ~ 5x10^7 cm^{-3}), but otherwise conventional, hypercompact HII region. The source is becoming optically thin at 229 GHz, and the H30alpha recombination line, at VLSR = 23.2 +/- 0.5 km/sec, is prominent in its spectrum. A Lyman continuum flux of 5x10^{45} photons/sec, consistent with that expected from a B star, is required to maintain the ionization. Supplementary 90 GHz observations were made to measure the H41alpha and H42alpha recombination lines toward BN. Published 43 and 86 GHz data suggest that SrcI brightened with respect to BN over the 15 year period from 1994 to 2009.

We report the results of radio interferometric observations of the 21-micron source IRAS 22272+5435 in the CO J=2-1 line. 21-micron sources are carbon-rich objects in the post-AGB phase of evolution which show an unidentified emission feature at 21 micron. Since 21-micron sources usually also have circumstellar molecular envelopes, the mapping of CO emission from the envelope will be useful in tracing the nebular structure. From observations made with the Combined Array for Research in Millimeter-wave Astronomy (CARMA), we find that a torus and spherical wind model can explain only part of the CO structure. An additional axisymmetric region created by the interaction between an invisible jet and ambient material is suggested.

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift {\gamma}-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to {\gamma}-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.

We report new 350 micron polarization observations of the thermal dust emission from the cores surrounding the low-mass, Class 0 YSOs L1527, IC348-SMM2 and B335. We have inferred magnetic field directions from these observations, and have used them together with results in the literature to determine whether magnetically regulated core-collapse and star-formation models are consistent with the observations. These models predict a pseudo-disk with its symmetry axis aligned with the core magnetic field. The models also predict a magnetic field pinch structure on a scale less than or comparable to the infall radii for these sources. In addition, if the core magnetic field aligns (or nearly aligns) the core rotation axis with the magnetic field before core collapse, then the models predict the alignment (or near alignment) of the overall pinch field structure with the bipolar outflows in these sources. We show that if one includes the distorting effects of bipolar outflows on magnetic fields, then in general the observational results for L1527 and IC348-SMM2 are consistent with these magnetically regulated models. We can say the same for B335 only if we assume the distorting effects of the bipolar outflow on the magnetic fields within the B335 core are much greater than for L1527 and IC348-SMM2. We show that the energy densities of the outflows in all three sources are large enough to distort the magnetic fields predicted by magnetically regulated models.

NGC 2024, a sites of massive star formation, have complex internal structures caused by cal heating by young stars, outflows, and stellar winds. These complex cloud structures lead to intricate emission line shapes. The goal of this paper is to show that the complex line shapes of 12 CO lines in NGC 2024 can be explained consistently with a model, whose temperature and velocity structure are based on the well-established scenario of a PDR and the Blister model. We present velocity-resolved spectra of seven CO lines ranging from J=3 to J=13, and we combined these data with CO high-frequency data from the ISO satellite. We find that the bulk of the molecular cloud associated with NGC 2024 consists of warm (75 K) and dense (9e5 cm-3) gas. An additional hot (~ 300 K) component, located at the interface of the HII region and the molecular cloud, is needed to explain the emission of the high-J CO lines. Deep absorption notches indicate that very cold material (20 K) exists in front of the warm material, too. A temperature and column density structure consistent with those predicted by PDR models, combined with the velocity structure of a Blister model, appropriately describes the observed emission line profiles of this massive star forming region. This case study of NGC 2024 shows that, with physical insights into these complex regions and careful modeling, multi-line observations of CO can be used to derive detailed physical conditions in massive star forming regions.

Deuterated ions are abundant in cold (T=10 K), dense (n=10^5 cm^-3) regions, in which CO is frozen out onto dust grains. In such environments, the deuterium fractionation of such ions can exceed the elemental abundance ratio of D/H by a factor of 10^4. In this paper we use the deuterium fractionation to investigate the evolutionary state of Class 0 protostars. In a sample of 20 protostellar objects, we found a clear correlation between the N2D+/N2H+ ratio and evolutionary tracers. As expected, the coolest, i.e. the youngest, objects show the largest deuterium fractionation. Furthermore, we find that sources with a high N2D+/N2H+ ratio show clear indication for infall.

Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud (LMC) is ideally suited to study star formation in such an environment. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR spectral lines of the main carbon carrying species, CO, CI and CII, which originate in the surface layers of molecular clouds illuminated by the UV radiation of the newly formed, young stars. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 telescope of the 12CO J = 4-3, J = 7-6, and 13CO J = 4-3 rotational and [CI] 3P1-3P0 and 3P2-3P1 fine-structure transitions. The 13CO J =4-3 and [CI] 3P2-3P1 transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W region has temperatures of about 80 K and densities of about 10^4 cm^-3. The estimated C to CO abundance ratio close to unity is substantially higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [CII] line intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about \chi ~220 and an average density of the clump ensemble of about 10^5 cm^-3, thus confirming the presence of high density material in the LMC N159W region.

Studying molecular gas in the central regions of the star burst galaxies NGC4945 and Circinus enables us to characterize the physical conditions and compare them to previous local and high-z studies. We estimate temperature, molecular density and column densities of CO and atomic carbon. Using model predictions we give a range of estimated CO/C abundance ratios. Using the new NANTEN2 4m sub-millimeter telescope in Pampa La Bola, Chile, we observed for the first time CO 4-3 and [CI] 3P1-3 P0 at the centers of both galaxies at linear scale of 682 pc and 732 pc respectively. We compute the cooling curves of 12CO and 13CO using radiative transfer models and estimate the physical conditions of CO and [CI]. The centers of NGC4945 and Circinus are very [CI] bright objects, exhibiting [CI] 3P1 - 3 P0 luminosities of 91 and 67Kkms-1kpc2, respectively. The [CI] 3P1-3 P0/CO 4-3 ratio of integrated intensities are large at 1.2 in NGC4945 and 2.8 in Circinus. Combining previous CO J= 1-0, 2-1 and 3-2 and 13CO J= 1-0, 2-1 studies with our new observations, the radiative transfer calculations give a range of densities, n(H2) = 10^3-3*104^cm-3, and a wide range of kinetic temperatures, Tkin = 20 - 100K, depending on the density. Future CO J= 7-6 and [CI] 2-1 observations will be important to resolve the ambiguity in the physical conditions and confirm the model predictions.

The Carina region is an excellent astrophysical laboratory for studying the feedback mechanisms of newly born, very massive stars within their natal giant molecular clouds (GMCs) at only 2.35 kpc distance. We use a clumpy PDR model to analyse the observed intensities of atomic carbon and CO and to derive the excitation conditions of the gas. The NANTEN2-4m submillimeter telescope was used to map the [CI] 3P1-3P0, 3P2-3P1 and CO 4-3, 7-6 lines in two 4'x4' regions of Carina where molecular material interfaces with radiation from the massive star clusters. One region is the northern molecular cloud near the compact OB cluster Tr14, and the second region is in the molecular cloud south of etaCar and Tr16. These data were combined with 13CO SEST spectra, HIRES/IRAS 60um and 100um maps of the FIR continuum, and maps of 8um IRAC/Spitzer and MSX emission. We used the HIRES far-infrared dust data to create a map of the FUV field heating the gas. The northern region shows an FUV field of a few 1000 in Draine units while the field of the southern region is about a factor 10 weaker. We constructed models consisting of an ensemble of small spherically symmetric PDR clumps within the 38" beam (0.43pc), which follow canonical power-law mass and mass-size distributions. We find that an average local clump density of 2x10**5 cm-3 is needed to reproduce the observed line emission at two selected interface positions. Stationary, clumpy PDR models reproduce the observed cooling lines of atomic carbon and CO at two positions in the Carina Nebula.

The THz atmospheric windows centered at roughly 1.3 and 1.5~THz, contain numerous spectral lines of astronomical importance, including three high-J CO lines, the N+ line at 205 microns, and the ground transition of para-H2D+. The CO lines are tracers of hot (several 100K), dense gas; N+ is a cooling line of diffuse, ionized gas; the H2D+ line is a non-depleting tracer of cold (~20K), dense gas. As the THz lines benefit the study of diverse phenomena (from high-mass star-forming regions to the WIM to cold prestellar cores), we have built the CO N+ Deuterium Observations Receiver (CONDOR) to further explore the THz windows by ground-based observations. CONDOR was designed to be used at the Atacama Pathfinder EXperiment (APEX) and Stratospheric Observatory For Infrared Astronomy (SOFIA). CONDOR was installed at the APEX telescope and test observations were made to characterize the instrument. The combination of CONDOR on APEX successfully detected THz radiation from astronomical sources. CONDOR operated with typical Trec=1600K and spectral Allan variance times of 30s. CONDOR's first light observations of CO 13-12 emission from the hot core Orion FIR4 (= OMC1 South) revealed a narrow line with T(MB) = 210K and delta(V)=5.4km/s. A search for N+ emission from the ionization front of the Orion Bar resulted in a non-detection. The successful deployment of CONDOR at APEX demonstrates the potential for making observations at THz frequencies from ground-based facilities.

The presence of H2D+ in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of many molecular species. We determine the H2D+ abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and dynamical structure. We observed the ortho-H2D+ ground-state line 1_10-1_11, the N2H+ J=4-3 line, and the H13CO+ 4-3 line with the APEX telescope. We report the probable detection of the o-H2D+ line at an intensity Tmb=0.22 +- 0.08 K and exclusively thermal line width, and find only upper limits to the N2H+ 4-3 and H13CO+ 4-3 intensities. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and excitation simulations reproduce the observed intensities and that of o-H2D+ in particular.