We present an empirical relation between the cold gas surface density ($\Sigma_{\rm gas}$) and the optical extinction (${\rm A_V}$) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the $\Sigma_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$\alpha$/H$\beta$). We find that the best fit for this relation is $\Sigma_{\rm gas} ({\rm M_\odot pc^{-2}})\sim~26~\times~ {\rm A_V}({\rm mag})$, and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the $\Sigma_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(H\alpha)}$. We find that at larger values of ${\rm EW(H\alpha)}$ (i.e., actively star-forming regions) this ratio tend to converge to the value expected for dust-star mixed geometries ($\sim$ 30 $\mathrm{M_{\odot} \,pc^{-2}\,mag^{-1}}$). On radial scales, we do not find a significant relation between the $\Sigma_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $\Sigma_{\rm gas}$ using ${\rm A_V}$ with compilations in the literature of the gas fraction on global and radial scales as well as with well known scaling relations such as the radial star-formation law and the $\Sigma_{\rm gas}$-$\Sigma_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $\Sigma_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.

About 60 ADASS participants discussed the evolving roles of users, developers, and managers of astronomical data reduction software. Special emphasis was placed on the role of the user in the era of Big Data. Is it really inevitable that the increasingly ignorant user will just have to be satisfied with the standard products from data processing centres?

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We investigate the prevalence, properties, and kinematics of extraplanar diffuse ionized gas (eDIG) in a sample of 25 edge-on galaxies selected from the CALIFA survey. We measure ionized gas scale heights from ${\rm H\alpha}$ and find that 90% have measurable scale heights with a median of $0.8^{+0.7}_{-0.4}$ kpc. From the ${\rm H\alpha}$ kinematics, we find that 60% of galaxies show a decrease in the rotation velocity as a function of height above the midplane. This lag is characteristic of eDIG, and we measure a median lag of 21 km s$^{-1}$ kpc$^{-1}$ which is comparable to lags measured in the literature. We also investigate variations in the lag with radius. $\rm H{\small I}$ lags have been reported to systematically decrease with galactocentric radius. We find both increasing and decreasing ionized gas lags with radius, as well as a large number of galaxies consistent with no radial lag variation, and investigate these results in the context of internal and external origins for the lagging ionized gas. We confirm that the ${\rm [S{\small II}]}$/${\rm H\alpha}$ and ${\rm [N{\small II}]}$/${\rm H\alpha}$ line ratios increase with height above the midplane as is characteristic of eDIG. The ionization of the eDIG is dominated by star-forming complexes (leaky ${\rm H{\small II}}$ regions). We conclude that the lagging ionized gas is turbulent ejected gas likely resulting from star formation activity in the disk as opposed to gas in the stellar thick disk or bulge. This is further evidence for the eDIG being a product of stellar feedback and for the pervasiveness of this WIM-like phase in many local star-forming galaxies.

We use NH3 inversion transitions to trace the dense gas in the NGC 1333 region of the Perseus molecular cloud. NH3(1,1) and NH3(2,2) maps covering an area of 102 square arcminutes at an angular resolution of ~3.7" are produced by combining VLA interferometric observations with GBT single dish maps. The combined maps have a spectral resolution of 0.14 km/s and a sensitivity of 4 mJy/beam. We produce integrated intensity maps, peak intensity maps and dispersion maps of NH3(1,1) and NH3(2,2) and a line-of-sight velocity map of NH3(1,1). These are used to derive the optical depth for the NH3(1,1) main component, the excitation temperature of NH3(1,1), and the rotational temperature, kinetic temperature and column density of NH3 over the mapped area. We compare these observations with the CARMA J=1-0 observations of N2H+ and H13CO+ and conclude that they all trace the same material in these dense star forming regions. From the NH3(1,1) velocity map, we find that a velocity gradient ridge extends in an arc across the entire southern part of NGC 1333. We propose that a large scale turbulent cell is colliding with the cloud, which could result in the formation of a layer of compressed gas. This region along the velocity gradient ridge is dotted with Class 0/I YSOs, that could have formed from local overdensities in the compressed gas leading to gravitational instabilities. The NH3(1,1) velocity dispersion map also has relatively high values along this region, thereby substantiating the shock layer argument.

Astrophysics papers often rely on software which may or may not be available, and URLs are often used as proxy citations for software and data. We extracted all URLs from two journals' 2015 research articles, removed those from certain long-term reliable domains, and tested the remainder to determine what percentage of these URLs were accessible in October 2018.

We present a new package for joint deconvolution of ALMA 12m, 7m, and Total Power (TP) data, dubbed ``Total Power Map to Visibilities (TP2VIS)". It converts a TP (single-dish) map into visibilities on the CASA platform, which can be input into deconvolvers (e.g., CLEAN) along with 12m and 7m visibilities. A manual is presented in the Github repository (https://github.com/tp2vis/distribute). Combining data from the different ALMA arrays is a driver for a number of science topics, namely those that probe size scales of extended and compact structures simultaneously. We test TP2VIS using model images, one with a single Gaussian and another that mimics the internal structures of giant molecular clouds. The result shows that the better uv coverage with TP2VIS visibilities helps the deconvolution process and reproduces the model image within errors of only 5% over two orders of magnitude in flux.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

We examined software usage in a sample set of astrophysics research articles published in 2015 and searched for source code for the software mentioned in these research papers. We categorized the software to indicate whether source code is available for download and whether there are restrictions to accessing it, and if source code is not available, whether some other form of the software, such as a binary, is. We also extracted hyperlinks from one journal's 2015 research articles, as links in articles can serve as an acknowledgment of software use and lead to data used in the research, and tested them to determine which of these URLs are still accessible. For our sample of 715 software instances in the 166 articles we examined, we were able to categorize 418 records as to availability of source code and found that 285 unique codes were used, 58% of which offer source code available online for download. Of the 2,558 hyperlinks extracted from 1,669 research articles, at best, 90% of them were available over our testing period.

We are members of the Astrophysics Source Code Library's Advisory Committee and its editor-in-chief. The Astrophysics Source Code Library (ASCL, ascl.net) is a successful initiative that advocates for open research software and provides an infrastructure for registering, discovering, sharing, and citing this software. Started in 1999, the ASCL has been expanding in recent years, with an average of over 200 codes added each year, and now houses over 1,600 code entries.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We present H13CO+ (J=1-0) and HNC (J=1-0) maps of regions in Serpens South, Serpens Main and NGC 1333 containing filaments. We also observe the Serpens regions using H13CN (J=1-0). These dense gas tracer molecular line observations carried out with CARMA have an angular resolution of ~7", a spectral resolution of ~0.16 km/s and a sensitivity of 50-100 mJy/beam. Although the large scale structure compares well with the Herschel dust continuum maps, we resolve finer structure within the filaments identified by Herschel. The H13CO+ emission distribution agrees with the existing CARMA N2H+ (J=1-0) maps; so they trace the same morphology and kinematics of the filaments. The H13CO+ maps additionally reveal that many regions have multiple structures partially overlapping in the line-of-sight. In two regions, the velocity differences are as high as 1.4 m/s. We identify 8 filamentary structures having typical widths of 0.03-0.08 pc in these tracers. At least 50% of the filamentary structures have distinct velocity gradients perpendicular to their major axis with average values in the range 4-10 km/s/pc. These findings are in support of the theoretical models of filament formation by 2-D inflow in the shock layer created by colliding turbulent cells. We also find evidence of velocity gradients along the length of two filamentary structures; the gradients suggest that these filaments are inflowing towards the cloud core.

The Astrophysics Source Code Library (ASCL, ascl.net), established in 1999, is a citable online registry of source codes used in research that are available for download; the ASCL's main purpose is to improve the transparency, reproducibility, and falsifiability of research. In 2017, improvements to the resource included real-time data backup for submissions and newly-published entries, improved cross-matching of research papers with software entries in ADS, and expansion of preferred citation information for the software in the ASCL.

Broad efforts are underway to capture metadata about research software and retain it across services; notable in this regard is the CodeMeta project. What metadata are important to have about (research) software? What metadata are useful for searching for codes? What would you like to learn about astronomy software? This BoF sought to gather information on metadata most desired by researchers and users of astro software and others interested in registering, indexing, capturing, and doing research on this software. Information from this BoF could conceivably result in changes to the Astrophysics Source Code Library (ASCL) or other resources for the benefit of the community or provide input into other projects concerned with software metadata.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

The Astrophysics Source Code Library (ASCL) is a free online registry of research codes; it is indexed by ADS and Web of Science and has over 1300 code entries. Its entries are increasingly used to cite software; citations have been doubling each year since 2012 and every major astronomy journal accepts citations to the ASCL. Codes in the resource cover all aspects of astrophysics research and many programming languages are represented. In the past year, the ASCL added dashboards for users and administrators, started minting Digital Objective Identifiers (DOIs) for software it houses, and added metadata fields requested by users. This presentation covers the ASCL's growth in the past year and the opportunities afforded it as one of the few domain libraries for science research codes.

Improving software citation and credit continues to be a topic of interest across and within many disciplines, with numerous efforts underway. In this Birds of a Feather (BoF) session, we started with a list of actionable ideas from last year's BoF and other similar efforts and worked alone or in small groups to begin implementing them. Work was captured in a common Google document; the session organizers will disseminate or otherwise put this information to use in or for the community in collaboration with those who contributed.

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.

The past year has seen movement on several fronts for improving software citation, including the Center for Open Science's Transparency and Openness Promotion (TOP) Guidelines, the Software Publishing Special Interest Group that was started at January's AAS meeting in Seattle at the request of that organization's Working Group on Astronomical Software, a Sloan-sponsored meeting at GitHub in San Francisco to begin work on a cohesive research software citation-enabling platform, the work of Force11 to "transform and improve" research communication, and WSSSPE's ongoing efforts that include software publication, citation, credit, and sustainability. Brief reports on these efforts were shared at the BoF, after which participants discussed ideas for improving software citation, generating a list of recommendations to the community of software authors, journal publishers, ADS, and research authors. The discussion, recommendations, and feedback will help form recommendations for software citation to those publishers represented in the Software Publishing Special Interest Group and the broader community.

We describe a framework that allows a scientist-user to easily query for information across all Virtual Observatory (VO) repositories and pull it back for analysis. This framework hides the gory details of meta-data remediation and data formatting from the user, allowing them to get on with search, retrieval and analysis of VO data as if they were drawn from a single source using a science based terminology rather than a data-centric one.

The Astrophysics Source Code Library (ASCL; ascl.net) is a free online registry of codes used in astronomy research; it currently contains over 900 codes and is indexed by ADS. The ASCL has recently moved a new infrastructure into production. The new site provides a true database for the code entries and integrates the WordPress news and information pages and the discussion forum into one site. Previous capabilities are retained and permalinks to ascl.net continue to work. This improvement offers more functionality and flexibility than the previous site, is easier to maintain, and offers new possibilities for collaboration. This presentation covers these recent changes to the ASCL.

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 the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

How do we as a community encourage the reuse of software for telescope operations, data processing, and calibration? How can we support making codes used in research available for others to examine? Continuing the discussion from last year Bring out your codes! BoF session, participants separated into groups to brainstorm ideas to mitigate factors which inhibit code sharing and nurture those which encourage code sharing. The BoF concluded with the sharing of ideas that arose from the brainstorming sessions and a brief summary by the moderator.

The Astrophysics Source Code Library (ASCL, http://ascl.net/) is an online registry of over 700 source codes that are of interest to astrophysicists, with more being added regularly. The ASCL actively seeks out codes as well as accepting submissions from the code authors, and all entries are citable and indexed by ADS. All codes have been used to generate results published in or submitted to a refereed journal and are available either via a download site or froman identified source. In addition to being the largest directory of scientist-written astrophysics programs available, the ASCL is also an active participant in the reproducible research movement with presentations at various conferences, numerous blog posts and a journal article. This poster provides a description of the ASCL and the changes that we are starting to see in the astrophysics community as a result of the work we are doing.

The Astrophysics Source Code Library, started in 1999, has in the past three years grown from a repository for 40 codes to a registry of over 700 codes that are now indexed by ADS. What comes next? We examine the future of the ASCL, the challenges facing it, the rationale behind its practices, and the need to balance what we might do with what we have the resources to accomplish.

We examine the formation of clusters of galaxies in numerical simulations of a QUMOND cosmogony with massive sterile neutrinos. Clusters formed in these exploratory simulations develop higher velocities than those found in {\Lambda}CDM simulations. The bulk motions of clusters attain about 1000 km/s by low redshift, comparable to observations whereas {\Lambda}CDM simulated clusters tend to fall short. Similarly, high pairwise velocities are common in cluster-cluster collisions like the Bullet cluster. There is also a propensity for the most massive clusters to be larger in QUMOND and to appear earlier than in {\Lambda}CDM, potentially providing an explanation for 'pink elephants' like El Gordo. However, it is not obvious that the cluster mass function can be recovered.

While software and algorithms have become increasingly important in astronomy, the majority of authors who publish computational astronomy research do not share the source code they develop, making it difficult to replicate and reuse the work. In this paper we discuss the importance of sharing scientific source code with the entire astrophysics community, and propose that journals require authors to make their code publicly available when a paper is published. That is, we suggest that a paper that involves a computer program not be accepted for publication unless the source code becomes publicly available. The adoption of such a policy by editors, editorial boards, and reviewers will improve the ability to replicate scientific results, and will also make the computational astronomy methods more available to other researchers who wish to apply them to their data.

We present the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey, an observational study of molecular gas in 31 star-forming galaxies from z=0.05 to z=0.5, with stellar masses of (4-30)x10^10 M_Sun and star formation rates of 4-100 M_Sun yr^-1. This survey probes a relatively un-observed redshift range in which the molecular gas content of galaxies is expected to have evolved significantly. To trace the molecular gas in the EGNoG galaxies, we observe the CO(1-0) and CO(3-2) rotational lines using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We detect 24 of 31 galaxies and present resolved maps of 10 galaxies in the lower redshift portion of the survey. We use a bimodal prescription for the CO to molecular gas conversion factor, based on specific star formation rate, and compare the EGNoG galaxies to a large sample of galaxies assembled from the literature. We find an average molecular gas depletion time of 0.76 \pm 0.54 Gyr for normal galaxies and 0.06 \pm 0.04 Gyr for starburst galaxies. We calculate an average molecular gas fraction of 7-20% at the intermediate redshifts probed by the EGNoG survey. By expressing the molecular gas fraction in terms of the specific star formation rate and molecular gas depletion time (using typical values), we also calculate the expected evolution of the molecular gas fraction with redshift. The predicted behavior agrees well with the significant evolution observed from z~2.5 to today.

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.

As observations of molecular gas in galaxies are pushed to lower star formation rate galaxies at higher redshifts, it is becoming increasingly important to understand the conditions of the gas in these systems to properly infer their molecular gas content. The rotational transitions of the carbon monoxide (CO) molecule provide an excellent probe of the gas excitation conditions in these galaxies. In this paper we present the results from the gas excitation sample of the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey at the Combined Array for Research in Millimeter-wave Astronomy (CARMA). This subset of the full EGNoG sample consists of four galaxies at z~0.3 with star formation rates of 40-65 M_Sun yr^-1 and stellar masses of ~2x10^11 M_Sun. Using the 3 mm and 1 mm bands at CARMA, we observe both the CO(1-0) and CO(3-2) transitions in these four galaxies in order to probe the excitation of the molecular gas. We report robust detections of both lines in three galaxies (and an upper limit on the fourth), with an average line ratio, r_31 = L'_CO(3-2) / L'_CO(1-0), of 0.46 \pm 0.07 (with systematic errors \lesssim 40%), which implies sub-thermal excitation of the CO(3-2) line. We conclude that the excitation of the gas in these massive, highly star-forming galaxies is consistent with normal star-forming galaxies such as local spirals, not starbursting systems like local ultra-luminous infrared galaxies. Since the EGNoG gas excitation sample galaxies are selected from the main sequence of star-forming galaxies, we suggest that this result is applicable to studies of main sequence galaxies at intermediate and high redshifts, supporting the assumptions made in studies that find molecular gas fractions in star forming galaxies at z~1-2 to be an order of magnitude larger than what is observed locally.

The Astrophysics Source Code Library (ASCL), founded in 1999, is a free on-line registry for source codes of interest to astronomers and astrophysicists. The library is housed on the discussion forum for Astronomy Picture of the Day (APOD) and can be accessed at http://ascl.net. The ASCL has a comprehensive listing that covers a significant number of the astrophysics source codes used to generate results published in or submitted to refereed journals and continues to grow. The ASCL currently has entries for over 500 codes; its records are citable and are indexed by ADS. The editors of the ASCL and members of its Advisory Committee were on hand at a demonstration table in the ADASS poster room to present the ASCL, accept code submissions, show how the ASCL is starting to be used by the astrophysics community, and take questions on and suggestions for improving the resource.

Progress is being made in code discoverability and preservation, but as discussed at ADASS XXI, many codes still remain hidden from public view. With the Astrophysics Source Code Library (ASCL) now indexed by the SAO/NASA Astrophysics Data System (ADS), the introduction of a new journal, Astronomy & Computing, focused on astrophysics software, and the increasing success of education efforts such as Software Carpentry and SciCoder, the community has the opportunity to set a higher standard for its science by encouraging the release of software for examination and possible reuse. We assembled representatives of the community to present issues inhibiting code release and sought suggestions for tackling these factors. The session began with brief statements by panelists; the floor was then opened for discussion and ideas. Comments covered a diverse range of related topics and points of view, with apparent support for the propositions that algorithms should be readily available, code used to produce published scientific results should be made available, and there should be discovery mechanisms to allow these to be found easily. With increased use of resources such as GitHub (for code availability), ASCL (for code discovery), and a stated strong preference from the new journal Astronomy & Computing for code release, we expect to see additional progress over the next few years.

The isolation, simple apparent structure, and low luminosity of the nearby spiral galaxy NGC 6503 make it an ideal candidate for an in-depth kinematic and photometric study. We introduce a new publicly available code, DiskFit, that implements procedures for fitting non-axisymmetries in either kinematic or photometric data. We use DiskFit to analyze new Halpha and CO velocity field data as well as HI kinematics from Greisen et al. to search for non-circular motions in the disc of NGC 6503. We find NGC 6503 to have remarkably regular gas kinematics that are well-described by rotation only. We also use DiskFit and a new Ks-band image of NGC 6503 to constrain photometric models of the disc, bar and bulge. We find the galaxy to be photometrically dominated by the disc. We find NGC 6503 to contain a faint bar and an exponential bulge which together contribute only ~5% of the total galaxy light. The combination of our kinematic and photometric DiskFit models suggest that NGC 6503 contains a weak, end-on bar that may have produced its Type II surface brightness profile but is unlikely to be responsible for its strong sigma-drop.

Based on Bremer et al. (2011) and Eckart et al. (2012) we report on simultaneous observations and modeling of the millimeter, near-infrared, and X-ray flare emission of the source Sagittarius A* (SgrA*) associated with the super-massive black hole at the Galactic Center. We study physical processes giving rise to the variable emission of SgrA* from the radio to the X-ray domain. To explain the statistics of the observed variability of the (sub-)mm spectrum of SgrA*, we use a sample of simultaneous NIR/X-ray flare peaks and model the flares using a synchrotron and SSC mechanism. The observations reveal flaring activity in all wavelength bands that can be modeled as the signal from adiabatically expanding synchrotron self-Compton (SSC) components. The model parameters suggest that either the adiabatically expanding source components have a bulk motion larger than v_exp or the expanding material contributes to a corona or disk, confined to the immediate surroundings of SgrA*. For the bulk of the synchrotron and SSC models, we find synchrotron turnover frequencies in the range 300-400 GHz. For the pure synchrotron models this results in densities of relativistic particles of the order of 10^6.5 cm^-3 and for the SSC models, the median densities are about one order of magnitude higher. However, to obtain a realistic description of the frequency-dependent variability amplitude of SgrA*, models with higher turnover frequencies and even higher densities are required. We discuss the results in the framework of possible deviations from equilibrium between particle and magnetic field energy. We also summarize alternative models to explain the broad-band variability of SgrA*.

Much of scientific progress now hinges on the reliability, falsifiability and reproducibility of computer source codes. Astrophysics in particular is a discipline that today leads other sciences in making useful scientific components freely available online, including data, abstracts, preprints, and fully published papers, yet even today many astrophysics source codes remain hidden from public view. We review the importance and history of source codes in astrophysics and previous efforts to develop ways in which information about astrophysics codes can be shared. We also discuss why some scientist coders resist sharing or publishing their codes, the reasons for and importance of overcoming this resistance, and alert the community to a reworking of one of the first attempts for sharing codes, the Astrophysics Source Code Library (ASCL). We discuss the implementation of the ASCL in an accompanying poster paper. We suggest that code could be given a similar level of referencing as data gets in repositories such as ADS.

Building on the "drPACS" contribution at ADASS XX of a simple Unix pipeline infrastructure, we implemented a pipeline toolkit using the package MIRIAD to combine Interferometric and Single Dish data (MIS). This was prompted by our observations made with the Combined Array For Research in Millimeter-wave Astronomy (CARMA) interferometer of the star-forming region NGC 1333, a large survey highlighting the new 23-element and singledish observing modes. The project consists of 20 CARMA datasets each containing interferometric as well as simultaneously obtained single dish data, for 3 molecular spectral lines and continuum, in 527 different pointings, covering an area of about 8 by 11 arcminutes. A small group of collaborators then shared this toolkit and their parameters via CVS, and scripts were developed to ensure uniform data reduction across the group. The pipeline was run end-to-end each night as new observations were obtained, producing maps that contained all the data to date. We will show examples of the scripts and data products. This approach could serve as a model for repeated calibration and mapping of large mixed-mode correlation datasets from ALMA.

Here we describe the Astrophysics Source Code Library (ASCL), which takes an active approach to sharing astrophysical source code. ASCL's editor seeks out both new and old peer-reviewed papers that describe methods or experiments that involve the development or use of source code, and adds entries for the found codes to the library. This approach ensures that source codes are added without requiring authors to actively submit them, resulting in a comprehensive listing that covers a significant number of the astrophysics source codes used in peer-reviewed studies. The ASCL now has over 340 codes in it and continues to grow. In 2011, the ASCL (http://ascl.net) has on average added 19 new codes per month. An advisory committee has been established to provide input and guide the development and expansion of the new site, and a marketing plan has been developed and is being executed. All ASCL source codes have been used to generate results published in or submitted to a refereed journal and are freely available either via a download site or from an identified source. This paper provides the history and description of the ASCL. It lists the requirements for including codes, examines the benefits of the ASCL, and outlines some of its future plans.

The mini-spiral is a feature of the interstellar medium in the central ~2 pc of the Galactic center. It is composed of several streamers of dust and ionised and atomic gas with temperatures between a few 100 K to 10^4 K. There is evidence that these streamers are related to the so-called circumnuclear disk of molecular gas and are ionized by photons from massive, hot stars in the central parsec. We attempt to constrain the emission mechanisms and physical properties of the ionized gas and dust of the mini-spiral region with the help of our multiwavelength data sets. Our observations were carried out at 1.3 mm and 3 mm with the mm interferometric array CARMA in California in March and April 2009, with the MIR instrument VISIR at ESO's VLT in June 2006, and the NIR Br-gamma with VLT NACO in August 2009. We present high resolution maps of the mini-spiral, and obtain a spectral index of 0.5 for Sgr A*, indicating an inverted synchrotron spectrum. We find electron densities within the range 0.8-1.5x10^4 cm-3 for the mini-spiral from the radio continuum maps, along with a dust mass contribution of ~0.25 solar masses from the MIR dust continuum, and extinctions ranging from 1.8-3 at 2.16 micron in the Br-gamma line. We observe a mixture of negative and positive spectral indices in our 1.3 mm and 3 mm observations of the extended emission of the mini-spiral, which we interpret as evidence that there are a range of contributions to the thermal free-free emission by the ionized gas emission and by dust at 1.3 mm.

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

We present high-resolution CARMA 230 GHz continuum imaging of nine deeply embedded protostars in the Serpens Molecular Cloud, including six of the nine known Class 0 protostars in Serpens. This work is part of a program to characterize disk and envelope properties for a complete sample of Class 0 protostars in nearby low-mass star forming regions. Here we present CARMA maps and visibility amplitudes as a function of uv-distance for the Serpens sample. Observations are made in the B, C, D, and E antenna configurations, with B configuration observations utilizing the CARMA Paired Antenna Calibration System. Combining data from multiple configurations provides excellent uv-coverage (4-500 klam), allowing us to trace spatial scales from 1e2 to 1e4 AU. We find evidence for compact disk components in all of the observed Class 0 protostars, suggesting that disks form at very early times (t<0.2 Myr) in Serpens. We make a first estimate of disk masses using the flux at 50 klam, where the contribution from the envelope should be negligible, assuming an unresolved disk. The resulting disk masses range from 0.04 Msun to 1.7 Msun, with a mean of approximately 0.2 Msun. Our high resolution maps are also sensitive to binary or multiple sources with separations > 250 AU, but significant evidence of multiplicity on scales <2000 AU is seen in only one source.

We report the CO(J=1-0) observations of M51 using both the Combined Array for Research in Millimeter Astronomy (CARMA) and the Nobeyama 45m telescope (NRO45). We describe a procedure for the combination of interferometer and single-dish data. In particular, we discuss (1) the joint imaging and deconvolution of heterogeneous data, (2) the weighting scheme based on the root-mean-square (RMS) noise in the maps, (3) the sensitivity and uv-coverage requirements, and (4) the flux recovery of a combined map. We generate visibilities from the single-dish map and calculate the noise of each visibility based on the RMS noise. Our weighting scheme, though it is applied to discrete visibilities in this paper, is applicable to grids in uv-space, and this scheme may advance in future software development. For a realistic amount of observing time, the sensitivities of the NRO45 and CARMA visibility data sets are best matched by using the single dish baselines only up to 4-6 kilo-lambda (about 1/4-1/3 of the dish diameter). The synthesized beam size is determined to conserve the flux between synthesized beam and convolution beam. The superior uv-coverage provided by the combination of CARMA long baseline data with 15 antennas and NRO45 short spacing data results in the high image fidelity, which is evidenced by the excellent overlap between even the faint CO emission and dust lanes in an optical HST image and PAH emission in an Spitzer 8 micron image.

We present 0.15" resolution observations of the 227 GHz continuum emission from the circumstellar disk around the FU-Orionis star PP 13S*. The data were obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) Paired Antenna Calibration System (C-PACS), which measures and corrects the atmospheric delay fluctuations on the longest baselines of the array in order to improve the sensitivity and angular resolution of the observations. A description of the C-PACS technique and the data reduction procedures are presented. C-PACS was applied to CARMA observations of PP 13S*, which led to a factor of 1.6 increase in the observed peak flux of the source, a 36% reduction in the noise of the image, and a 52% decrease in the measured size of the source major axis. The calibrated complex visibilities were fitted with a theoretical disk model to constrain the disk surface density. The total disk mass from the best fit model corresponds to 0.06 \msun, which is larger than the median mass of a disk around a classical T Tauri star. The disk is optically thick at a wavelength of 1.3 mm for orbital radii less than 48 AU. At larger radii, the inferred surface density of the PP 13S* disk is an order of magnitude lower than that needed to develop a gravitational instability.

We report on a successful, simultaneous observation and modelling of the millimeter (mm) to near-infrared (NIR) flare emission of the Sgr A* counterpart associated with the supermassive black hole at the Galactic centre (GC). We present a mm/sub-mm light curve of Sgr A* with one of the highest quality continuous time coverages and study and model the physical processes giving rise to the variable emission of Sgr A*.

We report the highest-fidelity observations of the spiral galaxy M51 in CO emission, revealing the evolution of giant molecular clouds (GMCs) vis-a-vis the large-scale galactic structure and dynamics. The most massive GMCs (so-called GMAs) are first assembled and then broken up as the gas flow through the spiral arms. The GMAs and their H2 molecules are not fully dissociated into atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large-scale galactic dynamics --their coagulation into GMAs is due to spiral arm streaming motions upon entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream side. In M51, the majority of the gas remains molecular from arm entry through the inter-arm region and into the next spiral arm passage.

Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicly and to embedding it in the scientific literature, which can lead to significant duplication of effort and an incomplete scientific record. We identify some of these structural disincentives and suggest a variety of approaches to address them, with the goals of raising the quality of astronomical software, improving the lot of scientist-authors, and providing benefits to the entire community, analogous to the benefits provided by open access to large survey and simulation datasets. Our aim is to open a conversation on how to move forward. We advocate that: (1) the astronomical community consider software as an integral and fundable part of facility construction and science programs; (2) that software release be considered as integral to the open and reproducible scientific process as are publication and data release; (3) that we adopt technologies and repositories for releasing and collaboration on software that have worked for open-source software; (4) that we seek structural incentives to make the release of software and related publications easier for scientist-authors; (5) that we consider new ways of funding the development of grass-roots software; (6) and that we rethink our values to acknowledge that astronomical software development is not just a technical endeavor, but a fundamental part of our scientific practice.

We present resolved images of the dust continuum emission from the debris disk around the young (80-200 Myr) solar-type star HD 107146 with CARMA at $\lambda$1.3 mm and the CSO at $\lambda$350 $\mu$m. Both images show that the dust emission extends over an $\sim$10\arcsec diameter region. The high resolution (3\arcsec) CARMA image further reveals that the dust is distributed in a partial ring with significant decrease in flux inward of 97 AU. Two prominent emission peaks appear within the ring separated by $\sim$140 degrees in position angle. The morphology of the dust emission is suggestive of dust captured into a mean motion resonance, which would imply the presence of a planet at an orbital radius of $\sim$45-75 AU.

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs.

We use CO and H alpha velocity fields to study the gas kinematics in the spiral arms and interarms of M51 (NGC 5194), and fit the 2D velocity field to estimate the radial and tangential velocity components as a function of spiral phase (arm distance). We find large radial and tangential streaming velocities, which are qualitatively consistent with the predictions of density wave theory and support the existence of shocks. The streaming motions are complex, varying significantly across the galaxy as well as along and between arms. Aberrations in the velocity field indicate that the disk is not coplanar, perhaps as far in as 20\arcsec\ (800 pc) from the center. Velocity profile fits from CO and H alpha are typically similar, suggesting that most of the H alpha emission originates from regions of recent star formation. We also explore vortensity and mass conservation conditions. Vortensity conservation, which does not require a steady state, is empirically verified. The velocity and density profiles show large and varying mass fluxes, which are inconsistent with a steady flow for a single dominant global spiral mode. We thus conclude that the spiral arms cannot be in a quasi-steady state in any rotating frame, and/or that out of plane motions may be significant.

Miriad is a radio interferometry data-reduction package, designed for taking raw data through to the image analysis stage. The Miriad project, begun in 1988, is now middle-aged. With the wisdom of hindsight, we review design decisions and some of Miriad's characteristics.

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 microns with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete dataset ranging from 3.6 to 1300 microns and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_int = 0.08 +/- 0.04 L_sun, placing it among the lowest luminosity protostars known. Though it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of Very Low Luminosity Objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_int > 1 L_sun) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation, and suggest a possible explanation for the very low luminosity of this source.

We present an analysis of the molecular gas distributions in the 29 barred and 15 unbarred spirals in BIMA SONG. For CO-bright galaxies, we confirm the conclusion by Sakamoto et al. (1999b) that barred spirals have higher molecular gas concentrations in the central kiloparsec. The SONG sample also includes 27 galaxies below the CO brightness limit used by Sakamoto et al. Even in these CO-faint galaxies we show that high central gas concentrations are more common in barred galaxies, consistent with radial inflow driven by the bar. However, there is a significant population of early-type (Sa--Sbc) barred spirals (6 of 19) that have little or no molecular gas detected in the nuclear region and out to the bar co-rotation radius. In these galaxies, the bar has already driven most of the gas within the bar to the nuclear region, where it has been consumed by star formation. The median nuclear gas mass is over four times higher in early type bars; since the gas consumption rate is an order of magnitude higher in early type bars, early types must have significantly higher bar-driven inflows. The lower inflow rates in late type bars can be attributed to differences in bar structure between early and late types. Despite bar-driven inflows, the data indicate that it is highly unlikely for a late type galaxy to evolve into an early type via bar-induced gas inflow. Nonetheless, secular evolutionary processes are undoubtedly present, and pseudo-bulges are inevitable; evidence for pseudo-bulges is likely to be clearest in early-type galaxies because of their high gas inflow rates and higher star formation activity (abridged).

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics.

Although nuclear rings of gas and star formation are common in barred spiral galaxies, current theories of why and how they form do not provide the level of detail needed to quantify the effect that these rings can have on the fueling of active galactic nuclei and on the evolution of their host galaxy. In this paper we use detailed modeling to show that existence of nuclear rings is directly related to the existence of the orbit family whose major axis is perpendicular to the major axis of the bar (x_2). We explore a large range of barred galaxy potentials and for each potential we use a two-dimensional hydrodynamic simulation to determine whether and at what radius a nuclear ring forms. We compare the results of the hydrodynamic simulations to numerical integrations of periodic orbits in a barred potential and show that the rings only form when a minimum amount of x_2 orbits exists. Because the rings migrate inwards with time as they accumulate gas, the radius at which a nuclear ring is seen does not give direct information on the shape of the rotation curve. We also show that the common assumption that nuclear rings are related to an inner Lindblad resonance is incorrect. In fact, we show that there is no RESONANCE at the inner Lindblad resonance in barred galaxies. We also compare the predictions of this theory to HST observations and show that it correctly predicts the observed gas and star formation morphology of nuclear rings.

We use data from the BIMA Survey of Nearby Galaxies (SONG) to investigate the relationship between ellipticity and central mass concentration in barred spirals. Existing simulations predict that bar ellipticity decreases as inflowing mass driven by the bar accumulates in the central regions, ultimately destroying the bar. Using the ratio of the bulge mass to the mass within the bar radius as an estimate of the central mass concentration, we obtain dynamical mass estimates from SONG CO 1-0 rotation curve data. We find an inverse correlation between bar ellipticity and central mass concentration, consistent with simulations of bar dissolution.

We compare the distribution of molecular gas and star formation activity in the bar region of six spirals (NGC 2903, NGC 3627, NGC 4321, NGC 5457, NGC 6946, & IC 342) from the BIMA Survey of Nearby Galaxies (SONG). The molecular gas, traced using the CO (J=1--0) emission line, is brightest along the leading edge of the stellar bar in the bar dust lanes. The star formation activity, traced using the Halpha emission line, is offset towards the leading side of the CO emission. A cross-correlation analysis shows that a) the HII regions are offset 0--800 pc on the leading side of the CO emission, b) the largest offsets are found in the strongest bars, and c) there is a wide range in offsets in a single bar with no systematic pattern as a function of the galacto-centric radius. The CO-Halpha offset constrains how stars may form depending on the gas flow. We examine possible star formation scenarios in context of the two main classes of bar gas flow simulations, the N-body/sticky particle and hydrodynamic models. Though both model gas flows are generally consistent with the observed offsets, we suggest the inclusion of a two- or multi-phase medium to improve the agreement between models and observation.

Like hardware, evolution of software has had a major impact on the field of particle simulations. This paper illustrates how simulation software has evolved, and where it can go. In addition, with the various ongoing Virtual Observatory efforts, producers of theory data should think more about sharing their data! Some examples are given of what we can do with our data and how to share it with our colleagues and observers. In the Appendix we summarize the findings of an informal data and software usage survey that we took during this conference.

Two dimensional velocity fields have been an important tool for nearly 30 years and are instrumental in understanding galactic mass distributions and deviations from an ideal galactic disk. Recently a number of new instruments have started to produce more detailed velocity fields of the disks and nuclear regions of galaxies. This paper summarizes some of the underlying techniques for constructing velocity fields and deriving rotation curves. It also urges to simulate observations from the data-cube stage to reject subtle biases in derived quantities such as rotation curves.

Although inner star-forming rings are common in optical images of barred spiral galaxies, observational evidence for the accompanying molecular gas has been scarce. In this paper we present images of molecular inner rings, traced using the CO (1-0) emission line, from the Berkeley-Illinois-Maryland-Association Survey of Nearby Galaxies (BIMA SONG). We detect inner ring CO emission from all five SONG barred galaxies classified as inner ring (type (r)). We also examine the seven SONG barred galaxies classified as inner spiral (type (s)); in one of these, NGC 3627, we find morphological and kinematic evidence for a molecular inner ring. Inner ring galaxies have been classified as such based on optical images, which emphasize recent star formation. We consider the possibility that there may exist inner rings in which star formation efficiency is not enhanced. However, we find that in NGC 3627 the inner ring star formation efficiency is enhanced relative to most other regions in that galaxy. We note that the SONG (r) galaxies have a paucity of CO and H alpha emission interior to the inner ring (except near the nucleus), while NGC 3627 has relatively bright bar CO and H alpha emission; we suggest that galaxies with inner rings such as NGC 3627 may be misclassified if there are significant amounts of gas and star formation in the bar.

During the last couple of years, observers have started to make plans for a Virtual Observatory, as a federation of existing data bases, connected through levels of software that enable rapid searches, correlations, and various forms of data mining. We propose to extend the notion of a Virtual Observatory by adding archives of simulations, together with interactive query and visualization capabilities, as well as ways to simulate observations of simulations in order to compare them with observations. For this purpose, we have already organized two small workshops, earlier in 2001, in Tucson and Aspen. We have also provided concrete examples of theory data, designed to be federated with a Virtual Observatory. These data stem from a project to construct an archive for our large-scale simulations using the GRAPE-6 (a 32-Teraflops special purpose computer for stellar dynamics). We are constructing interfaces by which remote observers can observe these simulations. In addition, these data will enable detailed comparisons between different simulations.

We report here results of high-resolution hydrodynamical simulations of gas flows in barred galaxies, with a focus on gas dynamics in the central kiloparsec. In a single bar with an Inner Lindblad Resonance, we find either near-circular motion of gas in the nuclear ring, or a spiral shock extending towards the galaxy center, depending on the sound speed in the gas. From a simple model of a dynamically-possible doubly barred galaxy with resonant coupling, we infer that the secondary bar is likely to end well inside its corotation. Such a bar cannot create shocks in the gas flow, and therefore will not reveal itself in color maps through straight dust lanes: the gas flows induced by it are different from those caused by the rapidly rotating main bars. In particular, we find that secondary stellar bars are unlikely to increase the mass inflow rate into the galactic nucleus.

In dense clusters a bewildering variety of interactions between stars can be observed, ranging from simple encounters to collisions and other mass-transfer encounters. With faster and special-purpose computers like GRAPE, the amount of data per simulation is now exceeding 1TB. Visualization of such data has now become a complex 4D data-mining problem, combining space and time, and finding interesting events in these large datasets. We have recently starting using the virtual reality simulator, installed in the Hayden Planetarium in the American Museum for Natural History, to tackle some of these problem. This work (http://www.astro.umd.edu/nemo/amnh/) reports on our first ``observations'', modifications needed for our specific experiments, and perhaps field ideas for other fields in science which can benefit from such immersion. We also discuss how our normal analysis programs can be interfaced with this kind of visualization.

We present multi-wavelength (interferometer and single-dish CO J=1-0, Halpha, broadband optical and near-infrared) observations of the classic barred spiral NGC 5383. We compare the observed central gas and dust morphology to the predictions of recent hydrodynamic simulations. In the nuclear region, our observations reveal three peaks lying along a S-shaped gas and dust distribution. In contrast, the model predicts a circumnuclear ring, not the observed S-shaped distribution; moreover, the predicted surface density contrast between the central gas accumulation and the bar dust lanes is an order of magnitude larger than observed. The discrepancies are not due to unexplored model parameter space or a nuclear bar but are probably due to the vigorous (7 solar masses per year) star formation activity in the center. As is common in similar bars, the star formation rate in the bar between the bar ends and the central region is low (~0.5 solar masses per yr), despite the high gas column density in the bar dust lanes; this is generally attributed to shear and shocks. We note a tendency for the HII regions to be associated with the spurs feeding the main bar dust lanes, but these are located on the leading side of the bar. We propose that stars form in the spurs, which provide a high column density but low shear environment. HII regions can therefore be found even on the leading side of the bar because the ionizing stars pass ballistically through the dust lane

Hierarchical structure formation models predict the existence of large numbers of low velocity dispersion dark halos. Galaxy surveys find far fewer galaxies than predicted by analytical estimates and numerical simulations. In this paper, we suggest that these dark halos are not missing, but have been merely misplaced in the galactic astronomy section of the journals: they are the High Velocity Clouds (HVCs). We review the predictions of our model for the Local Group origin of the HVCs and its implications for the formation and the evolution of our Galaxy. We describe recent observations that confirm many of earlier predictions and discuss future tests of the model.

We suggest that the high--velocity clouds (HVCs) are large clouds, with typical diameters of 25 kpc and containing 5e7 solar masses of neutral gas and 3e8 solar masses of dark matter, falling onto the Local Group; altogether the HVCs contain 10$^{11}$ solar masses of neutral gas. Our reexamination of the Local-Group hypothesis for the HVCs connects their properties to the hierarchical structure formation scenario and to the gas seen in absorbtion towards quasars. We begin by showing that at least one HVC complex (besides the Magellanic Stream) must be extragalactic at a distance > 40 kpc from the Galactic center, with a diameter > 20 kpc and a mass > 1e8 solar masses. We interpret the more distant HVCs as dark matter ``mini--halos'' moving along filaments towards the Local Group. Most poor galaxy groups should contain HI structures to large distances bound to the group. The HVCs are local analogues of the Lyman--limit clouds We argue that there is a Galactic fountain in the Milky Way, but that the fountain does not explain the origin of the HVCs. Our analysis of the HI data leads to the detection of a vertical infall of low-velocity gas towards the plane. This implies that the chemical evolution of the Galactic disk is governed by episodic infall of metal-poor HVC gas that only slowly mixes with the rest of the interstellar medium. The Local--Group infall hypothesis makes a number of testable predictions. The HVCs should have sub-solar metallicities. Their H$\alpha $ emission should be less than that seen from the Magellanic Stream. The clouds should not be seen in absorption to nearby stars. The clouds should be detectable in both emission and absorption around other groups.

We report results from a multiwavelength campaign to measure the simultaneous spectrum of the super-massive black hole candidate Sgr A* in the Galactic Center from cm to mm-wavelengths using the VLA, BIMA, the Nobeyama 45m, and the IRAM 30m telescopes. The observations confirm that the previously detected mm-excess is an intrinsic feature of the spectrum of Sgr A*. The excess can be interpreted as due to the presence of an ultra-compact component of relativistic plasma with a size of a few Schwarzschild radii near the black hole. If so, Sgr A* might offer the unique possibility to image the putative black hole against the background of this component with future mm-VLBI experiments.

We report first results from a multiwavelength campaign to measure the simultaneous spectrum of Sgr A* from cm to mm wavelengths. The observations confirm that the previously detected submm-excess is not due to variability; the presence of an ultracompact component with a size of a few Schwarzschild radii is inferred. In a VLA survey of LINER galaxies, we found Sgr A*-like nuclei in one quarter of the galaxies searched, suggesting a link between those low-power AGN and the Galactic Center.

Mass inflow in barred galaxies has been invoked to account for a wide variety of phenomena, but until now direct evidence for inflow has been lacking. We present Fabry-Perot H-alpha observations of the barred spiral galaxy NGC 1530 from which we determine velocities of the ionized gas for the entire region swept by the bar. We compare the velocity field to models of gas flow in barred spirals and show that it is well reproduced by ideal gas hydrodynamic models. Inspection of the models and observations reveals that gas entering the bar dust lanes streams directly down the dust lanes toward the 2 kpc radius nuclear ring. The models predict that approximately 20% of the gas flowing down the dust lane enters the nuclear ring; the remaining gas sprays around the ring to the other bar dust lane. The fraction of the gas entering the ring is relatively insensitive to the shape or size of the bar. Our observations of the velocity field and dust optical depth yield a mass inflow rate into the nuclear ring of 1 solar mass per year.