We present the galaxy luminosity functions (LFs) of four Hickson Compact Groups using image data from the Subaru Hyper Suprime-Cam. A distinct dip appeared in the faint-ends of all the LFs at $M_g\sim-12$. A similar dip was observed in the LFs of the galaxy clusters Coma and Centaurus. However, LFs in the Virgo, Hydra, and the field had flatter slopes and no dips. As the relative velocities among galaxies are lower in compact groups than in clusters, the effect of galaxy-galaxy interactions would be more significant in compact groups. The $M_g\sim-12$ dip of compact groups may imply that frequent galaxy-galaxy interactions would affect the evolution of galaxies, and the dip in LF could become a boundary between different galaxy populations.

We present here the Calar Alto Schmidt-Lemaitre Telescope (CASTLE) concept, a technology demonstrator for curved detectors, that will be installed at the Calar Alto Observatory (Spain). This telescope has a wide field of view (2.36x1.56 deg^2) and a design, optimised to generate a Point Spread Function with very low level wings and reduced ghost features, which makes it considerably less susceptible to several systematic effects usually affecting similar systems. These characteristics are particularly suited to study the low surface brightness Universe. CASTLE will be able to reach surface brightness orders of magnitude fainter than the sky background level and observe the extremely extended and faint features around galaxies such as tidal features, stellar halos, intra-cluster light, etc. CASTLE will also be used to search and detect astrophysical transients such as gamma ray bursts (GRB), gravitational wave optical counterparts, neutrino counterparts, etc. This will increase the number of precisely localized GRBs from 20% to 60% (in the case of Fermi/GMB GRBs).

Malin 1 is the largest known low surface brightness (LSB) galaxy, the archetype of so-called giant LSBs. The structure and the origin of such galaxies are still poorly understood, especially due to the lack of high-resolution kinematics and spectroscopic data. We use emission lines from spectroscopic observations of Malin 1 aiming to bring new constraints on the internal dynamics and star formation history of Malin 1. We have extracted a total of 16 spectra from different regions of Malin 1 and calculated the rotational velocities of these regions from the wavelength shifts and star formation rates from the observed H$\alpha$ emission line fluxes. We compare our data with existing data and models for Malin 1. For the first time we present the inner rotation curve of Malin 1, characterized in the radial range r < 10 kpc by a steep rise in the rotational velocity up to at least 350 km/s (with a large dispersion), which had not been observed previously. We use these data to study a suite of new mass models for Malin 1. We show that in the inner regions dynamics may be dominated by the stars (although none of our models can explain the highest velocities measured) but that at large radii a massive dark matter halo remains necessary. The H$\alpha$ fluxes derived star formation rates are consistent with an early-type disk for the inner region, and with the level found in extended UV galaxies for the outer parts of the giant disk of Malin 1. We also find signs of high metallicity but low dust content for the inner regions.

We present spatial variations of the CO J=2-1/1-0 line ratio in M83 using Total Power array data from ALMA. While the intensities of these two lines correlate tightly, the ratio varies over the disk, with a disk average ratio of 0.69, and shows the galactic center and a two-arm spiral pattern. It is high (>0.7) in regions of high molecular gas surface density, but ranges from low to high ratios in regions of low surface density. The ratio correlates well with the spatial distributions and intensities of FUV and IR emissions, with FUV being the best correlated. It also correlates better with the ratio of IR intensities (70/350mic), a proxy for dust temperature, than with the IR intensities. Taken together, these results suggest either a direct or indirect link between the dust heating by the interstellar radiation field and the condition of GMCs, even though no efficient mechanism is known for a thermal coupling of dust and bulk gas in GMCs.

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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

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 observed five giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) in the 12CO J=1-0 line using the Atacama Large Millimeter/submillimeter Array (ALMA). The sample includes four GMCs with some signs of star formation -- either YSOs, H II regions, and/or young clusters -- and one quiescent GMC without any sign of massive star formation. The data from the ALMA 12 m, 7 m, and Total-Power arrays are jointly deconvolved to obtain high-fidelity images at high spatial resolution (3 arcsec = 0.7 pc). The four star-forming GMCs show very complex structures with clumps and filaments. The quiescent GMC shows a relatively diffuse, extended emission distribution without prominent clumps or filaments. This difference is similar to that between structured molecular gas in Milky Way spiral arms and unstructured gas in the inter-arm regions. We characterize the difference with the brightness distribution function and brightness distribution index. In conjunction with other ALMA studies of GMCs in the LMC, the five GMCs tentatively form an evolutionary trend: from less structured, quiescent GMCs to more structured, actively star-forming GMCs. A future ALMA study will be able to map molecular clouds over the LMC and reveal the evolutionary sequence of molecular clouds.

We present a [CII] 158um map of the entire M51 (including M51b) grand--design spiral galaxy observed with the FIFI-LS instrument on SOFIA. We compare the [CII] emission with the total far--infrared (TIR) intensity and star formation rate(SFR) surface density maps (derived using H_alpha and 24um emission) to study the relationship between [CII] and the star formation activity in a variety of environments within M51 on scales of 16" corresponding to ~660 pc. We find that [CII] and the SFR surface density are well correlated in the central, spiral arm, and inter-arm regions. The correlation is in good agreement with that found for a larger sample of nearby galaxies at kpc scales. We find that the SFR, and [CII] and TIR luminosities in M51 are dominated by the extended emission in M51's disk. The companion galaxy M51b, however, shows a deficit of [CII] emission compared with the TIR emission and SFR surface density, with [CII] emission detected only in the S-W part of this galaxy. The [CII] deficit is associated with an enhanced dust temperature in this galaxy. We interpret the faint [CII] emission in M51b to be a result of suppressed star formation in this galaxy, while the bright mid- and far-infrared emission, which drive the TIR and SFR values, are powered by other mechanisms. A similar but less pronounced effect is seen at the location of the black hole in M51's center. The observed [CII] deficit in M51b suggests that this galaxy is a valuable laboratory to study the origin of the apparent [CII] deficit observed in ultra-luminous galaxies.

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

The peculiar velocities of galaxies cause their redshift-space clustering to depend on the angle to the line-of-sight, providing a key test of gravitational physics on cosmological scales. These effects may be described using a multipole expansion of the clustering measurements. Focussing on Fourier-space statistics, we present a new analysis of the effect of the survey window function, and the variation of the line-of-sight across a survey, on the modelling of power spectrum multipoles. We determine the joint covariance of the Fourier-space multipoles in a Gaussian approximation, and indicate how these techniques may be extended to studies of overlapping galaxy populations via multipole cross-power spectra. We apply our methodology to one of the widest-area galaxy redshift surveys currently available, the 6-degree Field Galaxy Survey, deducing a normalized growth rate f*sigma_8(z=0.06) = 0.38 +/- 0.12 in the low-redshift Universe, in agreement with previous analyses of this dataset using different techniques. Our framework should be useful for processing future wide-angle galaxy redshift surveys.

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

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Low redshift measurements of Baryon Acoustic Oscillations (BAO) test the late time evolution of the Universe and are a vital probe of Dark Energy. Over the past decade both the 6-degree Field Galaxy Survey (6dFGS) and Sloan Digital Sky Survey (SDSS) have provided important distance constraints at $z < 0.3$. In this paper we re-evaluate the cosmological information from the BAO detection in 6dFGS making use of HOD populated COLA mocks for a robust covariance matrix and taking advantage of the now commonly implemented technique of density field reconstruction. For the 6dFGS data, we find consistency with the previous analysis, and obtain an isotropic volume averaged distance measurement of $D_{V}(z_{\mathrm{eff}}=0.097) = 372\pm17(r_{s}/r_{s}^{\mathrm{fid}})\,\mathrm{Mpc}$, which has a non-Gaussian likelihood outside the $1\sigma$ region. We combine our measurement from both the post-reconstruction clustering of 6dFGS and SDSS MGS offering the most robust constraint to date in this redshift regime, $D_{V}(z_{\mathrm{eff}}=0.122)=539\pm17(r_{s}/r^{\mathrm{fid}}_{s})\,\mathrm{Mpc}$. These measurements are consistent with standard $\Lambda\mathrm{CDM}$ and after fixing the standard ruler using a Planck prior on $\Omega_{m}h^{2}$, the joint analysis gives $H_{0}=64.0\pm3.5\,\mathrm{kms}^{-1}\mathrm{Mpc}^{-1}$. In the near future both the Taipan Galaxy Survey and the Dark Energy Spectroscopic Instrument (DESI) will improve this measurement to $1\%$ at low redshift.

We study the presence of tidal features associated with ultra diffuse galaxies (UDGs) in galaxy clusters. Specifically, we stack deep Subaru images of UDGs in the Coma cluster to determine whether they show position angle twists at large radii. Selecting galaxies with central surface brightness (g-band) >24 magarcsec^{-2} and projected half-light radius >1.5 kpc, we identify 287 UDGs in the Yagi et al. 2016 catalog of low surface brightness Coma objects. The UDGs have apparent spheroidal shapes with median Sersic index ~ 0.8 and median axis-ratio ~ 0.7. The images are processed by masking all background objects and rotating to align the major axis before stacking them in bins of properties such as axis ratio, angle of major axis with respect to the cluster center and separation from cluster center. Our image stacks reach further than 7 kpc (~4r_e). Analysis of the isophotes of the stacks reveal that the ellipticity remains constant up to the last measured point, which means that the individual galaxies have a non-varying position angle and axis ratio and show no evidence for tidal disruption out to ~ 4r_e. We demonstrate this explicitly by comparing our stacks with stacks of model UDGs with and without tidal features in their outskirts. We infer that the average tidal radius of the Coma UDGs is >7 kpc and estimate that the average dark matter fraction within the tidal radius of the UDGs inhabiting the innermost 0.5 Mpc of Coma is > 99%.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]

We present a new test of gravitational physics by comparing the growth rate of cosmic structure measured around voids with that measured around galaxies in the same large-scale structure dataset, the low-redshift 6-degree Field Galaxy Survey. By fitting a self-consistent Redshift Space Distortion model to the 2D galaxy-galaxy and void-galaxy correlation functions, we recover growth rate values f\sigma_8 = 0.42 \pm 0.06 and 0.39 \pm 0.11, respectively. The environmental-dependence of cosmological statistics can potentially discriminate between modified-gravity scenarios which modulate the growth rate as a function of scale or environment and test the underlying assumptions of homogeneity and isotropy.

We present the Palomar Transient Factory discoveries and the photometric and spectroscopic observations of PTF11kmb and PTF12bho. We show that both transients have properties consistent with the class of calcium-rich gap transients, specifically lower peak luminosities and rapid evolution compared to ordinary supernovae, and a nebular spectrum dominated by [Ca II] emission. A striking feature of both transients is their host environments: PTF12bho is an intra-cluster transient in the Coma Cluster, while PTF11kmb is located in a loose galaxy group, at a physical offset ~150 kpc from the most likely host galaxy. Deep Subaru imaging of PTF12bho rules out an underlying host system to a limit of $M_R > -8.0$ mag, while Hubble Space Telescope imaging of PTF11kmb reveals a marginal counterpart that, if real, could be either a background galaxy or a globular cluster. We show that the offset distribution of Ca-rich gap transients is significantly more extreme than that seen for Type Ia supernovae or even short-hard gamma-ray bursts (sGRBs). Thus, if the offsets are caused by a kick, they require larger kick velocities and/or longer merger times than sGRBs. We also show that almost all Ca-rich gap transients found to date are in group and cluster environments with elliptical host galaxies, indicating a very old progenitor population; the remote locations could partially be explained by these environments having the largest fraction of stars in the intra-group/intra-cluster light following galaxy-galaxy interactions.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

We present measurements of the growth rate of cosmological structure from the modelling of the anisotropic galaxy clustering measured in the final data release of the VIPERS survey. The analysis is carried out in configuration space and based on measurements of the first two even multipole moments of the anisotropic galaxy auto-correlation function, in two redshift bins spanning the range $0.5 < z < 1.2$. We provide robust and cosmology-independent corrections for the VIPERS angular selection function, allowing recovery of the underlying clustering amplitude at the percent level down to the Mpc scale. We discuss several improvements on the non-linear modelling of redshift-space distortions (RSD) and perform detailed tests of a variety of approaches against a set of realistic VIPERS-like mock realisations. This includes using novel fitting functions to describe the velocity divergence and density power spectra $P_{\theta\theta}$ and $P_{\delta\theta}$ that appear in RSD models. These tests show that we are able to measure the growth rate with negligible bias down to separations of $5h^{-1}Mpc$. Interestingly, the application to real data shows a weaker sensitivity to the details of non-linear RSD corrections compared to mock results. We obtain consistent values for the growth rate times the matter power spectrum normalisation parameter of $f\sigma_8=0.55\pm 0.12$ and $0.40\pm0.11$ at effective redshifts of $z = 0.6$ and $z=0.86$ respectively. These results are in agreement with standard cosmology predictions assuming Einstein gravity in a $\Lambda \rm{CDM}$ background.

The outskirts of galaxies offer extreme environments where we can test our understanding of the formation, evolution, and destruction of molecules and their relationship with star formation and galaxy evolution. We review the basic equations that are used in normal environments to estimate physical parameters like the molecular gas mass from CO line emission and dust continuum emission. Then we discuss how those estimates may be affected when applied to the outskirts, where the average gas density, metallicity, stellar radiation field, and temperature may be lower. We focus on observations of molecular gas in the outskirts of the Milky Way, extragalactic disk galaxies, early-type galaxies, groups, and clusters. The scientific results show the versatility of molecular gas, as it has been used to trace Milky Way spiral arms out to a galactocentric radius of 15 kpc, to study star formation in extended ultraviolet disk galaxies, to probe galaxy interactions in polar ring S0 galaxies, and to investigate ram pressure stripping in clusters. We highlight the physical stimuli that accelerate the formation of molecular gas, including internal processes such as spiral arm compression and external processes such as interactions.

There is a linear relation between the mass of dense gas, traced by the HCN(1-0) luminosity, and the star formation rate (SFR), traced by the far-infrared luminosity. Recent observations of galactic disks have shown some systematic variations. In order to explore the SFR-dense gas link at high resolution ($\sim 4"$, $\sim 150$ pc) in the outer disk of an external galaxy, we have mapped a region about 5 kpc from the center along the northern spiral arm of M51 in the HCN(1-0), HCO$^+$(1-0) and HNC(1-0) emission lines using the Northern Extended Millimeter Array (NOEMA) interferometer. The HCN and HCO$^+$ lines were detected in 6 giant molecular associations (GMAs) while HNC emission was only detected in the two brightest GMAs. One of the GMAs hosts a powerful HII region and HCN is stronger than HCO$^+$ there. Comparing with observations of GMAs in the disks of M31 and M33 at similar angular resolution ($\sim 100$ pc), we find that GMAs in the outer disk of M51 are brighter in both HCN and HCO$^+$ lines by a factor of 3 on average. However, the $I_{HCN}/I_{CO}$ and $I_{HCO^+}/I_{CO}$ ratios are similar to the ratios in nearby galactic disks and the Galactic plane. Using the Herschel 70 $\mu$m data to trace the total IR luminosity at the resolution of the GMAs, we find that both the L$_{IR}$-L$_{HCN}$ and L$_{IR}$-L$_{HCO^+}$ relations in the outer disk GMAs are consistent with the proportionality between the L$_{IR}$ and the dense gas mass established globally in galaxies within the scatter. The IR/HCN and IR/HCO$^+$ ratios of the GMAs vary by a factor of 3, probably depending on whether massive stars are forming or not.

Theoretical studies on the response of interstellar gas to a gravitational potential disc with a quasi-stationary spiral arm pattern suggest that the gas experiences a sudden compression due to standing shock waves at spiral arms. This mechanism, called a galactic shock wave, predicts that gas spiral arms move from downstream to upstream of stellar arms with increasing radius inside a corotation radius. In order to investigate if this mechanism is at work in the grand-design spiral galaxy M51, we have measured azimuthal offsets between the peaks of stellar mass and gas mass distributions in its two spiral arms. The stellar mass distribution is created by the spatially resolved spectral energy distribution fitting to optical and near infrared images, while the gas mass distribution is obtained by high-resolution CO and HI data. For the inner region (r < 150"), we find that one arm is consistent with the galactic shock while the other is not. For the outer region, results are less certain due to the narrower range of offset values, the weakness of stellar arms, and the smaller number of successful offset measurements. The results suggest that the nature of two inner spiral arms are different, which is likely due to an interaction with the companion galaxy.

Low surface brightness galaxies (LSBGs) represent a significant percentage of local galaxies but their formation and evolution remain elusive. They may hold crucial information for our understanding of many key issues (i.e., census of baryonic and dark matter, star formation in the low density regime, mass function). The most massive examples - the so called giant LSBGs - can be as massive as the Milky Way, but with this mass being distributed in a much larger disk. Malin 1 is an iconic giant LSBG, perhaps the largest disk galaxy known. We attempt to bring new insights on its structure and evolution on the basis of new images covering a wide range in wavelength. We have computed surface brightness profiles (and average surface brightnesses in 16 regions of interest), in six photometric bands (FUV, NUV, u, g, i, z). We compared these data to various models, testing a variety of assumptions concerning the formation and evolution of Malin 1. We find that the surface brightness and color profiles can be reproduced by a long and quiet star-formation history due to the low surface density; no significant event, such as a collision, is necessary. Such quiet star formation across the giant disk is obtained in a disk model calibrated for the Milky Way, but with an angular momentum approximately 20 times larger. Signs of small variations of the star-formation history are indicated by the diversity of ages found when different regions within the galaxy are intercompared.For the first time, panchromatic images of Malin 1 are used to constrain the stellar populations and the history of this iconic example among giant LSBGs. Based on our model, the extreme disk of Malin 1 is found to have a long history of relatively low star formation (about 2 Msun/yr). Our model allows us to make predictions on its stellar mass and metallicity.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $\mu$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

We analyze radial and azimuthal variations of the phase balance between the molecular and atomic ISM in the Milky Way. In particular, the azimuthal variations -- between spiral arm and interarm regions -- are analyzed without any explicit definition of spiral arm locations. We show that the molecular gas mass fraction, i.e., fmol=H2/ (HI+H2) in mass, varies predominantly in the radial direction: starting from ~100% at the center, remaining ~>50% (~>60%) to R~6kpc, and decreasing to ~10-20% (~50%) at R=8.5 kpc when averaged over the whole disk thickness (in the mid plane). Azimuthal, arm-interarm variations are secondary: only ~20%, in the globally molecule-dominated inner MW, but becoming larger, ~40-50%, in the atom-dominated outskirts. This suggests that in the inner MW, the gas stays highly molecular (fmol>50%) as it goes from an interarm region, into a spiral arm, and back into the next interarm region. Stellar feedback does not dissociate molecules much, and the coagulation and fragmentation of molecular clouds dominate the evolution of the ISM at these radii. The trend differs in the outskirts, where the gas phase is globally atomic (fmol<50%). The HI and H2 phases cycle through spiral arm passage there. These different regimes of ISM evolution are also seen in external galaxies (e.g., LMC, M33, and M51). We explain the radial gradient of fmol by a simple flow continuity model. The effects of spiral arms on this analysis are illustrated in Appendix.

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

We report development of a simple and affordable radio interferometer suitable as an educational laboratory experiment. With the increasing importance of interferometry in astronomy, the lack of educational interferometers is an obstacle to training the future generation of astronomers. This interferometer provides the hands-on experience needed to fully understand the basic concepts of interferometry. The design of this interferometer is based on the Michelson & Pease stellar optical interferometer, but operates at a radio wavelength (~11 GHz; ~2.7cm); thus the requirement for optical accuracy is much less stringent. We utilize a commercial broadcast satellite dish and feedhorn. Two flat side mirrors slide on a ladder, providing baseline coverage. This interferometer resolves and measures the diameter of the Sun, a nice daytime experiment which can be carried out even in marginal weather (i.e., partial cloud cover). Commercial broadcast satellites provide convenient point sources for comparison to the Sun's extended disk. We describe the mathematical background of the adding interferometer, the design and development of the telescope and receiver system, and measurements of the Sun. We present results from a students' laboratory report.

We measure the bulk flow of the local Universe using the 6dF Galaxy Survey peculiar velocity sample (6dFGSv), the largest and most homogeneous peculiar velocity sample to date. 6dFGSv is a Fundamental Plane sample of $\sim10^4$ peculiar velocities covering the whole southern hemisphere for galactic latitude $|b| > 10^\circ$, out to redshift ${z=0.0537}$. We apply the `Minimum Variance' bulk flow weighting method, which allows us to make a robust measurement of the bulk flow on scales of $50$ and $70\,h^{-1}{\rm Mpc}$. We investigate and correct for potential bias due to the lognormal velocity uncertainties, and verify our method by constructing $\Lambda{\rm CDM}$ 6dFGSv mock catalogues incorporating the survey selection function. For a hemisphere of radius $50\,h^{-1}{\rm Mpc}$ we find a bulk flow amplitude of $U=248\pm58\,{\rm km}\,{\rm s}^{-1}$ in the direction $(l,b) = (318^\circ\pm20^\circ,40^\circ\pm13^\circ)$, and for $70\,h^{-1}{\rm Mpc}$ we find $U=243\pm58\,{\rm km}\,{\rm s}^{-1}$, in the same direction. Our measurement gives us a constraint on $\sigma_8$ of $1.01^{+1.07}_{-0.58}$. Our results are in agreement with other recent measurements of the direction of the bulk flow, and our measured amplitude is consistent with a $\Lambda{\rm CDM}$ prediction.

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langle\delta\delta v v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

In this work, we investigate the molecular gas and star formation properties in the barred spiral galaxy NGC 6946 using multiple molecular lines and star formation tracers. High-resolution image (100 pc) of $^{13}$CO (1-0) is created by single dish NRO45 and interferometer CARMA for the inner 2 kpc disk, which includes the central region (nuclear ring and bar) and the offset ridges of the primary bar. Single dish HCN (1-0) observations were also made to constrain the amount of dense gas. Physical properties of molecular gas are inferred by (1) the Large Velocity Gradient (LVG) calculations using our observations and archival $^{12}$CO (1-0), $^{12}$CO(2-1) data, (2) dense gas fraction suggested by HCN to $^{12}$CO (1-0) luminosity ratio, and (3) infrared color. The results show that the molecular gas in the central region is warmer and denser than that of the offset ridges. Dense gas fraction of the central region is similar with that of LIRGs/ULIRGs, while the offset ridges are close to the global average of normal galaxies. The coolest and least dense region is found in a spiral-like structure, which was misunderstood to be part of the southern primary bar in previous low-resolution observations. Star formation efficiency (SFE) changes by ~ 5 times in the inner disk. The variation of SFE agrees with the prediction in terms of star formation regulated by galactic bar. We find a consistency between star-forming region and the temperature inferred by the infrared color, suggesting that the distribution of sub-kpc scale temperature is driven by star formation.

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

We study the large-scale clustering of galaxies in the overlap region of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample and the WiggleZ Dark Energy Survey. We calculate the auto-correlation and cross-correlation functions in the overlap region of the two datasets and detect a Baryon Acoustic Oscillation (BAO) signal in each of them. The BAO measurement from the cross-correlation function represents the first such detection between two different galaxy surveys. After applying density-field reconstruction we report distance-scale measurements $D_V r_s^{\rm fid} / r_s = (1970 \pm 47, 2132 \pm 67, 2100 \pm 200)$ Mpc from CMASS, the cross-correlation and WiggleZ, respectively. We use correlated mock realizations to calculate the covariance between the three BAO constraints. The distance scales derived from the two datasets are consistent, and are also robust against switching the displacement fields used for reconstruction between the two surveys. This approach can be used to construct a correlation matrix, permitting for the first time a rigorous combination of WiggleZ and CMASS BAO measurements. Using a volume-scaling technique, our result can also be used to combine WiggleZ and future CMASS DR12 results. Finally, we use the cross-correlation function measurements to show that the relative velocity effect, a possible source of systematic uncertainty for the BAO technique, is consistent with zero for our samples.

We present an accurate and fast framework for generating mock catalogues including low-mass halos, based on an implementation of the COmoving Lagrangian Acceleration (COLA) technique. Multiple realisations of mock catalogues are crucial for analyses of large-scale structure, but conventional N-body simulations are too computationally expensive for the production of thousands of realisations. We show that COLA simulations can produce accurate mock catalogues with a moderate computation resource for low- to intermediate- mass galaxies in $10^{12} M_\odot$ haloes, both in real and redshift space. COLA simulations have accurate peculiar velocities, without systematic errors in the velocity power spectra for k < 0.15 h/Mpc, and with only 3-per-cent error for k < 0.2 h/Mpc. We use COLA with 10 time steps and a Halo Occupation Distribution to produce 600 mock galaxy catalogues of the WiggleZ Dark Energy Survey. Our parallelized code for efficient generation of accurate halo catalogues is publicly available at github.com/junkoda/cola_halo.

We report high spatial resolution observations of Giant Molecular Clouds (GMCs) in the nearby spiral galaxies NGC 6946, M101 and NGC 628 obtained with the CARMA telescope. We observed CO(1-0) over regions with active star formation, and higher resolution observations of CO(2-1) have allowed us to resolve some of the largest GMCs. Using a Bayesian fitting approach, we generate scaling relations for the sizes, line widths, and virial masses of the structures identified in this work. We do not find evidence for a tight power law relation between size and line width, although the limited dynamic range in cloud size remains a clear issue in our analysis. Additionally, we use a Bayesian approach to analyze the Kennicutt-Schmidt (K-S) relation for the identified structures. We find that the distribution of slopes are broadly distributed, mainly due to the limited dynamic range of our measured H2 mass surface density, and being most consistent with super-linear relations. On the other hand, when we use the Bayesian approach to analyze the K-S relation for a uniform grid, the distributions of slopes is consistent with sub-linear relations. On-arm regions tend to have higher star formation rates than inter-arm regions. As in NGC 6946, in M101 we find regions where the star formation efficiency shows marked peaks at specific galoctocentric radii. On the other hand, the distribution of SFE in NGC 628 is more contiguous. We hypothesize that differences in the distribution of SFE may be indicative of different processes driving the spiral structure.

The anisotropic galaxy 2-point correlation function (2PCF) allows measurement of the growth of large-scale structures from the effect of peculiar velocities on the clustering pattern. We present new measurements of the auto- and cross- correlation function multipoles of 69,180 WiggleZ and 46,380 BOSS-CMASS galaxies sharing an overlapping volume of ~0.2 (Gpc/h)^3. Analysing the redshift-space distortions (RSD) of galaxy 2-point statistics for these two galaxy tracers, we test for systematic errors in the modelling depending on galaxy type and investigate potential improvements in cosmological constraints. We build a large number of mock galaxy catalogues to examine the limits of different RSD models in terms of fitting scales and galaxy type, and to study the covariance of the measurements when performing joint fits. For the galaxy data, fitting the monopole and quadrupole of the WiggleZ 2PCF on scales 24<s<80 Mpc/h produces a measurement of the normalised growth rate $f\sigma_8$(z=0.54)=0.409$\pm$0.059, whereas for the CMASS galaxies we found a consistent constraint of $f\sigma_8$(z=0.54)=0.466$\pm$0.074. When combining the measurements, accounting for the correlation between the two surveys, we obtain $f\sigma_8$(z=0.54)=0.413$\pm$0.054, in agreement with the LCDM-GR model of structure growth and with other survey measurements.

We report the discovery of 854 ultra diffuse galaxies (UDGs) in the Coma cluster using deep R band images, with partial B, i, and Halpha band coverage, obtained with the Subaru telescope. Many of them (332) are Milky Way-sized with very large effective radii of r_e>1.5kpc. This study was motivated by the recent discovery of 47 UDGs by van-Dokkum et al. (2015); our discovery suggests >1,000 UDGs after accounting for the smaller Subaru field. The new UDGs show a distribution concentrated around the cluster center, strongly suggesting that the great majority are (likely longtime) cluster members. They are a passively evolving population, lying along the red sequence in the CM diagram with no Halpha signature. Star formation was, therefore, quenched in the past. They have exponential light profiles, effective radii re ~ 800 pc- 5 kpc, effective surface brightnesses mu_e(R)=25-28 mag arcsec-2, and stellar masses ~1x10^7 - 5x10^8Msun. There is also a population of nucleated UDGs. Some MW-sized UDGs appear closer to the cluster center than previously reported; their survival in the strong tidal field, despite their large sizes, possibly indicates a large dark matter fraction protecting the diffuse stellar component. The indicated baryon fraction ~<1% is less than the cosmic average, and thus the gas must have been removed from the possibly massive dark halo. The UDG population appears to be elevated in the Coma cluster compared to the field, indicating that the gas removal mechanism is related primarily to the cluster environment.

We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $\gamma = 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $\gamma = 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $\Lambda$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2\sigma$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.

We report the discovery of a new dwarf galaxy (NGC6503-d1) during the Subaru extended ultraviolet (XUV) disk survey. It is a likely companion of the spiral galaxy NGC6503. The resolved images, in B, V, R, i, and Halpha, show an irregular appearance due to bright stars with underlying, smooth and unresolved stellar emission. It is classified as the transition type (dIrr/dSph). Its structural properties are similar to those of the dwarfs in the Local Group, with a V absolute magnitude ~ -10.5, half-light radius ~400 pc, and central surface brightness ~25.2. Despite the low stellar surface brightness environment, one HII region was detected, though its Halpha luminosity is low, indicating an absence of any appreciable O-stars at the current epoch. The presence of multiple stellar populations is indicated by the color-magnitude diagram of ~300 bright resolved stars and the total colors of the dwarf, with the majority of its total stellar mass ~4x10^6 Msun in an old stellar population.

We report short-time variations in the plasma tail of C/2013 R1(Lovejoy). A series of short (two to three minutes) exposure images with the 8.2-m Subaru telescope shows faint details of filaments and their motions over 24 minutes observing duration. We identified rapid movements of two knots in the plasma tail near the nucleus (~ 3 x 10^5 km). Their speeds are 20 and 25 km/s along the tail and 3.8 and 2.2 km/s across it, respectively. These measurements set a constraint on an acceleration model of plasma tail and knots as they set the initial speed just after their formation. We also found a rapid narrowing of the tail. After correcting the motion along the tail, the narrowing speed is estimated to be ~ 8 km/s. These rapid motions suggest the need for high time-resolution studies of comet plasma tails with a large telescope.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

We present scale-dependent measurements of the normalised growth rate of structure $f\sigma_{8}(k, z=0)$ using only the peculiar motions of galaxies. We use data from the 6-degree Field Galaxy Survey velocity sample (6dFGSv) together with a newly-compiled sample of low-redshift $(z < 0.07)$ type Ia supernovae. We constrain the growth rate in a series of $\Delta k \sim 0.03 h{\rm Mpc^{-1}}$ bins to $\sim35\%$ precision, including a measurement on scales $>300 h^{-1}{\rm Mpc}$, which represents one of the largest-scale growth rate measurement to date. We find no evidence for a scale dependence in the growth rate, or any statistically significant variation from the growth rate as predicted by the {\it Planck} cosmology. Bringing all the scales together, we determine the normalised growth rate at $z=0$ to $\sim15\%$ in a manner {\it independent} of galaxy bias and in excellent agreement with the constraint from the measurements of redshift-space distortions from 6dFGS. We pay particular attention to systematic errors. We point out that the intrinsic scatter present in Fundamental-Plane and Tully-Fisher relations is only Gaussian in logarithmic distance units; wrongly assuming it is Gaussian in linear (velocity) units can bias cosmological constraints. We also analytically marginalise over zero-point errors in distance indicators, validate the accuracy of all our constraints using numerical simulations, and demonstrate how to combine different (correlated) velocity surveys using a matrix `hyper-parameter' analysis. Current and forthcoming peculiar velocity surveys will allow us to understand in detail the growth of structure in the low-redshift universe, providing strong constraints on the nature of dark energy.

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent \beta values (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.

We present significant improvements in cosmic distance measurements from the WiggleZ Dark Energy Survey, achieved by applying the reconstruction of the baryonic acoustic feature technique. We show using both data and simulations that the reconstruction technique can often be effective despite patchiness of the survey, significant edge effects and shot-noise. We investigate three redshift bins in the redshift range 0.2<$z$<1, and in all three find improvement after reconstruction in the detection of the baryonic acoustic feature and its usage as a standard ruler. We measure model independent distance measures $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ of 1716 $\pm$ 83 Mpc, 2221 $\pm$ 101 Mpc, 2516 $\pm$ 86 Mpc (68% CL) at effective redshifts z = 0.44, 0.6, 0.73, respectively, where $D_{\mathrm V}$ is the volume-average-distance, and $r_{\mathrm s}$ is the sound horizon at the end of the baryon drag epoch. These significantly improved 4.8, 4.5 and 3.4 percent accuracy measurements are equivalent to those expected from surveys with up to 2.5 times the volume of WiggleZ. These measurements are fully consistent with cosmologies allowed by the analyses of the Planck Collaboration and the Sloan Digital Sky Survey.We provide the $D_{\mathrm V}(r_{\mathrm s}^\mathrm{fid}/r_{\mathrm s})$ posterior probability distributions and their covariances. When combining these measurements with temperature fluctuations measurements of Planck, the polarization of WMAP9, and the 6dF Galaxy Survey baryonic acoustic feature, we do not detect deviations from a flat LCDM model. Assuming this model we constrain the current expansion rate to $H_0$ = 67.15 $\pm$ 0.98 kms$^{-1}$Mpc$^{-1}$. Allowing the equation of state of dark energy to vary we obtain $w_\mathrm{DE}$ = -1.080 $\pm$ 0.135. When assuming a curved LCDM model we obtain a curvature value of $\Omega_{\mathrm K}$ = -0.0043 $\pm$ 0.0047.

Peculiar velocity surveys, which measure galaxy peculiar velocities directly from standard candles in addition to redshifts, can provide strong constraints on the linear growth rate of cosmological large-scale structure at low redshift. The improvement originates from the physical relationship between galaxy density and peculiar velocity, which substantially reduces cosmic variance. We present the results of Fisher matrix forecasts of correlated fields of galaxy density and velocity. Peculiar velocity can improve the growth rate constraints by about a factor of two compared to density alone, if we can use all the information for wavenumber k < 0.2 h/Mpc. Future peculiar velocity surveys, TAIPAN, and the all-sky HI surveys, WALLABY and WNSHS, can measure the growth rate, f*sigma8 to 3 per cent at z ~ 0.025. Although the velocity subsample is about an order of magnitude smaller than the redshift sample from the same survey, it improves the constraint by 40 per cent compared to the same survey without velocity measurements. Peculiar velocity surveys can also measure the growth rate as a function of wavenumber with 15-30 per cent uncertainties in bins with widths 0.01 h/Mpc in the range 0.02 h/Mpc < k < 0.1 h/Mpc, which is a large improvement over galaxy density only. Such measurements on very large scales can detect signatures of modified gravity or non-Gaussianity through scale-dependent growth rate or galaxy bias. We use N-body simulations to improve the modelling of auto- and cross-power spectra of galaxy density and peculiar velocity by introducing a new redshift-space distortion term to the velocity, which has been neglected in previous studies. The velocity power spectrum has a damping in redshift space, which is larger than that naively expected from the similar effect in the galaxy power spectrum.

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

Submillimeter excess emission has been reported at 500 microns in a handful of local galaxies, and previous studies suggest that it could be correlated with metal abundance. We investigate the presence of an excess submillimeter emission at 500 microns for a sample of 20 galaxies from the Key Insights on Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) that span a range of morphologies and metallicities (12+log(O/H)=7.8-8.7). We probe the far-infrared (IR) emission using images from the Spitzer Space Telescope and Herschel Space Observatory in the wavelength range 24-500 microns. We model the far-IR peak of the dust emission with a two-temperature modified blackbody and measure excess of the 500 micron photometry relative to that predicted by our model. We compare the submillimeter excess, where present, with global galaxy metallicity and, where available, resolved metallicity measurements. We do not find any correlation between the 500 micron excess and metallicity. A few individual sources do show excess (10-20%) at 500 microns; conversely, for other sources, the model overpredicts the measured 500 micron flux density by as much as 20%, creating a 500 micron "deficit". None of our sources has an excess larger than the calculated 1-sigma uncertainty, leading us to conclude that there is no substantial excess at submillimeter wavelengths at or shorter than 500 microns in our sample. Our results differ from previous studies detecting 500 micron excess in KINGFISH galaxies largely due to new, improved photometry used in this study.

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.

We report a super-linear correlation for the star formation law based on new CO($J$=1-0) data from the CARMA and NOBEYAMA Nearby-galaxies (CANON) CO survey. The sample includes 10 nearby spiral galaxies, in which structures at sub-kpc scales are spatially resolved. Combined with the star formation rate surface density traced by H$\alpha$ and 24 $\mu$m images, CO($J$=1-0) data provide a super-linear slope of $N$ = 1.3. The slope becomes even steeper ($N$ = 1.8) when the diffuse stellar and dust background emission is subtracted from the H$\alpha$ and 24 $\mu$m images. In contrast to the recent results with CO($J$=2-1) that found a constant star formation efficiency (SFE) in many spiral galaxies, these results suggest that the SFE is not independent of environment, but increases with molecular gas surface density. We suggest that the excitation of CO($J$=2-1) is likely enhanced in the regions with higher star formation and does not linearly trace the molecular gas mass. In addition, the diffuse emission contaminates the SFE measurement most in regions where star formation rate is law. These two effects can flatten the power law correlation and produce the apparent linear slope. The super linear slope from the CO($J$=1-0) analysis indicates that star formation is enhanced by non-linear processes in regions of high gas density, e.g., gravitational collapse and cloud-cloud collisions.

We resolve 182 individual giant molecular clouds (GMCs) larger than 2.5 $\times$ 10$^{5}$ \Msun in the inner disks of five large nearby spiral galaxies (NGC 2403, NGC 3031, NGC 4736, NGC 4826, and NGC 6946) to create the largest such sample of extragalactic GMCs within galaxies analogous to the Milky Way. Using a conservatively chosen sample of GMCs most likely to adhere to the virial assumption, we measure cloud sizes, velocity dispersions, and $^{12}$CO (J=1-0) luminosities and calculate cloud virial masses. The average conversion factor from CO flux to H$_{2}$ mass (or \xcons) for each galaxy is 1-2 \xcounits, all within a factor of two of the Milky Way disk value ($\sim$2 \xcounits). We find GMCs to be generally consistent within our errors between the galaxies and with Milky Way disk GMCs; the intrinsic scatter between clouds is of order a factor of two. Consistent with previous studies in the Local Group, we find a linear relationship between cloud virial mass and CO luminosity, supporting the assumption that the clouds in this GMC sample are gravitationally bound. We do not detect a significant population of GMCs with elevated velocity dispersions for their sizes, as has been detected in the Galactic center. Though the range of metallicities probed in this study is narrow, the average conversion factors of these galaxies will serve to anchor the high metallicity end of metallicity-\xco trends measured using conversion factors in resolved clouds; this has been previously possible primarily with Milky Way measurements.

Studies of GMC evolution in galactic disks were limited to local, predominantly atom-rich small galaxies in the past, but have now been expanded to typical spiral galaxies with a rich molecular content. The evolution appears quite different between the two environments. GMCs exist almost exclusively along HI spiral arms and filaments in the disks of local small galaxies (LMC, M33), suggesting that GMCs form and end their short lives there. However, in a more molecular-rich environment (MW, M51), GMCs are present everywhere independent of HI structures. Indeed, the molecular gas fraction remains high and almost constant during arm passage into the next inter-arm region. The gas remains molecular, presumably in GMCs, for a long time. A transitional case has been found recently in the central region of the atom-rich galaxy M33 - GMCs do not coincide with HI there. Evolution of the physical conditions of molecular gas from spiral arms to inter-arm regions is also being revealed in molecule-rich galaxies. An increase of the CO J=2-1 and 1-0 line ratio in spiral arms in M51 suggests density and/or temperature increases by a factor of 2-3 in GMCs in the arms, compared to their counterparts in the inter-arm regions. An analysis of high-resolution Milky Way survey data revealed that the fraction of dense (or warm) clumps increases dramatically in the spiral arms.

[abridged] Aims. The aim of the present paper is to provide new and more detailed relations at the kpc scale between the gas surface density and the face-on optical depth directly calibrated on galaxies, in order to compute the attenuation not only for semi-analytic models but also observationally as new and upcoming radio observatories are able to trace gas ever farther in the Universe. Methods. We have selected a sample of 4 nearby resolved galaxies and a sample of 27 unresolved galaxies from the Herschel Reference Survey and the Very Nearby Galaxies Survey, for which we have a large set of multi-wavelength data from the FUV to the FIR including metallicity gradients for resolved galaxies, along with radio HI and CO observations. For each pixel in resolved galaxies and for each galaxy in the unresolved sample, we compute the face-on optical depth from the attenuation determined with the CIGALE SED fitting code and an assumed geometry. We determine the gas surface density from HI and CO observations with a metallicity-dependent XCO factor. Results. We provide new, simple to use, relations to determine the face-on optical depth from the gas surface density, taking the metallicity into account, which proves to be crucial for a proper estimate. The method used to determine the gas surface density or the face-on optical depth has little impact on the relations except for galaxies that have an inclination over 50d. Finally, we provide detailed instructions on how to compute the attenuation practically from the gas surface density taking into account possible information on the metallicity. Conclusions. Examination of the influence of these new relations on simulated FUV and IR luminosity functions shows a clear impact compared to older oft-used relations, which in turn could affect the conclusions drawn from studies based on large scale cosmological simulations.

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

We report measurements of geometric offsets between gas spiral arms and associated star forming regions in the grand-design spiral galaxy M51. These offsets are a suggested measure of the star formation timescale after the compression of gas at spiral arm entry. A surprising discrepancy, by an order of magnitude, has been reported in recent offset measurements in nearby spiral galaxies. Measurements using CO and H-alpha emission find large and ordered offsets in M51. On the contrary, small or non-ordered offsets have been found using the HI 21cm and 24 micron emissions, possible evidence against gas flow through spiral arms, and thus against the conventional density-wave theory with a stationary spiral pattern. The goal of this paper is to understand the cause of this discrepancy. We investigate potential causes by repeating those previous measurements using equivalent data, methods, and parameters. We find offsets consistent with the previous measurements and conclude that the difference of gas tracers, i.e., HI versus CO, is the primary cause. The HI emission is contaminated significantly by the gas photo-dissociated by recently-formed stars and does not necessarily trace the compressed gas, the precursor of star formation. The HI gas and star forming regions coincide spatially and tend to show small offsets. We find mostly positive offsets with substantial scatter between CO and H-alpha, suggesting gas flow through spiral arms (i.e., density-wave) though the spiral pattern may not necessarily be stationary.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

We report systematic variations in the CO(2-1)/CO(1-0) line ratio (R) in M51. The ratio shows clear evidence for the evolution of molecular gas from the upstream interarm regions, passage into the spiral arms and back into the downstream interarm regions. In the interarm regions, R is typically low <0.7 (and often 0.4-0.6); this is similar to the ratios observed in Galactic giant molecular clouds (GMCs) with low far-IR luminosities. However, the ratio rises to >0.7 (often 0.8-1.0) in the spiral arms, particularly at their leading (downstream) edge. R is also high, 0.8-1.0, in the central region. An LVG calculation provides insight into the changes in the gas physical conditions between the arm and interarm regions: cold and low density gas (~10K, ~300cm-3) is required for the interarm GMCs, but this gas must become warmer and/or denser in the more active star forming spiral arms. R is higher in areas of high 24micron brightness (an approximate tracer of star formation rate surface density) and high CO(1-0) integrated intensity (a well-calibrated tracer of total molecular gas surface density). The systematic enhancement of the CO(2-1) line relative to CO(1-0) in luminous star forming regions suggests that some caution is needed when using CO(2-1) as a tracer of bulk molecular gas mass.

We present the results of a Nobeyama 45-m water maser and ammonia survey of all 94 northern GLIMPSE Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based on their extended 4.5 micron emission. We observed the ammonia (1,1), (2,2), and (3,3) inversion lines, and detect emission towards 97%, 63%, and 46% of our sample, respectively (median rms ~50 mK). The water maser detection rate is 68% (median rms ~0.11 Jy). The derived water maser and clump-scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser associations. Water masers and warm dense gas, as indicated by emission in the higher-excitation ammonia transitions, are most frequently detected towards EGOs also associated with both Class I and II methanol masers. 95% (81%) of such EGOs are detected in water (ammonia(3,3)), compared to only 33% (7%) of EGOs without either methanol maser type. As populations, EGOs associated with Class I and/or II methanol masers have significantly higher ammonia linewidths, column densities, and kinetic temperatures than EGOs undetected in methanol maser surveys. However, we find no evidence for statistically significant differences in water maser properties (such as maser luminosity) among any EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, we find no correlation between isotropic water maser luminosity and clump number density. Water maser luminosity is weakly correlated with clump (gas) temperature and clump mass.

We explore the development of structures in molecular gas in the Milky Way by applying the analysis of the brightness distribution function (BDF) and the brightness distribution index (BDI) in the archival data from the Boston University-Five College Radio Astronomy Observatory 13CO J=1-0 Galactic Ring Survey. The BDI measures the fractional contribution of spatially confined bright molecular emission over faint emission extended over large areas. This relative quantity is largely independent of the amount of molecular gas and of any conventional, pre-conceived structures, such as cores, clumps, or giant molecular clouds. The structured molecular gas traced by higher BDI is located continuously along the spiral arms in the Milky Way in the longitude-velocity diagram. This clearly indicates that molecular gas changes its structure as it flows through the spiral arms. Although the high-BDI gas generally coincides with H II regions, there is also some high-BDI gas with no/little signature of ongoing star formation. These results support a possible evolutionary sequence in which unstructured, diffuse gas transforms itself into a structured state on encountering the spiral arms, followed by star formation and an eventual return to the unstructured state after the spiral arm passage.

We present high spatial resolution observations of Giant Molecular Clouds (GMCs) in the eastern part of the nearby spiral galaxy NGC 6946 obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We have observed 12CO(1-0), 12CO(2-1) and 13CO(1-0), achieving spatial resolutions of 5.4" x 5.0", 2.5" x 2.0" and 5.6" x 5.4" respectively over a region of 6 x 6 kpc. This region extends from 1.5 kpc to 8 kpc galactocentric radius, thus avoiding the intense star formation in the central kpc. We have recovered short-spacing u-v components by using single dish observations from the Nobeyama 45m and IRAM 30m telescopes. Using the automated CPROPS algorithm we identified 44 CO cloud complexes in the 12CO(1-0) map and 64 GMCs in the 12CO(2-1) maps. The sizes, line widths, and luminosities of the GMCs are similar to values found in other extragalactic studies. We have classified the clouds into on-arm and inter-arm clouds based on the stellar mass density traced by the 3.6 um map. On-arm clouds present in general higher star formation rates than clouds located on inter-arm regions. Although the star formation efficiency shows no systematic trend with galactocentric radius, some on-arm clouds -- which are more luminous and more massive compared to inter-arm GMCs -- are also forming stars more efficiently than the rest of the identified GMCs. These structures appear to be located on two specific regions in the spiral arms. One of them shows a strong gradient, suggesting that this region of high star formation efficiency may be the result of gas flow convergence.

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

Using infrared imaging from the Herschel Space Observatory, observed as part of the VNGS, we investigate the spatially resolved dust properties of the interacting Whirlpool galaxy system (NGC 5194 and NGC 5195), on physical scales of 1 kpc. Spectral energy distribution modelling of the new infrared images in combination with archival optical, near- through mid-infrared images confirms that both galaxies underwent a burst of star formation ~370-480 Myr ago and provides spatially resolved maps of the stellar and dust mass surface densities. The resulting average dust-to-stellar mass ratios are comparable to other spiral and spheroidal galaxies studied with Herschel, with NGC 5194 at log M(dust)/M(star)= -2.5+/-0.2 and NGC 5195 at log M(dust)/M(star)= -3.5+/-0.3. The dust-to-stellar mass ratio is constant across NGC 5194 suggesting the stellar and dust components are coupled. In contrast, the mass ratio increases with radius in NGC 5195 with decreasing stellar mass density. Archival mass surface density maps of the neutral and molecular hydrogen gas are also folded into our analysis. The gas-to-dust mass ratio, 94+/-17, is relatively constant across NGC 5194. Somewhat surprisingly, we find the dust in NGC 5195 is heated by a strong interstellar radiation field, over 20 times that of the ISRF in the Milky Way, resulting in relatively high characteristic dust temperatures (~30 K). This post-starburst galaxy contains a substantial amount of low-density molecular gas and displays a gas-to-dust ratio (73+/-35) similar to spiral galaxies. It is unclear why the dust in NGC 5195 is heated to such high temperatures as there is no star formation in the galaxy and its active galactic nucleus is 5-10 times less luminous than the one in NGC 5194, which exhibits only a modest enhancement in the amplitude of its ISRF.

We present our new deep optical imaging and long-slit spectroscopy for Arp 220 that is the archetypical ULIRG in the local universe. Our sensitive Ha imaging has newly revealed large-scale, Ha absorption, i.e., post-starburst regions in this merger; one is found in the eastern superbubble and the other is in the two tidal tails that are clearly reveled in our deep optical imaging. The size of Ha absorption region in the eastern bubble is 5 kpc x 7.5 kpc and the observed Ha equivalent widths are ~2 A +- 0.2 A. The sizes of the northern and southern Ha-absorption tidal tails are ~5 kpc x 10 kpc and ~6 kpc x 20 kpc, respectively. The observed Ha equivalent widths range from 4 A to 7 A. In order to explain the presence of the two post-starburst tails, we suggest a possible multiple-merger scenario for Arp 220 in which two post-starburst disk-like structures merged into one, and then caused the two tails. This favors that Arp 220 is a multiple merging system composed of four or more galaxies, arising from a compact group of galaxies. Taking our new results into account, we discuss a star formation history in the last 1 Gyr in Arp 220.

Recent observational results indicate that the functional shape of the spatially-resolved star formation-molecular gas density relation depends on the spatial scale considered. These results may indicate a fundamental role of sampling effects on scales that are typically only a few times larger than those of the largest molecular clouds. To investigate the impact of this effect, we construct simple models for the distribution of molecular clouds in a typical star-forming spiral galaxy, and, assuming a power-law relation between SFR and cloud mass, explore a range of input parameters. We confirm that the slope and the scatter of the simulated SFR-molecular gas surface density relation depend on the size of the sub-galactic region considered, due to stochastic sampling of the molecular cloud mass function, and the effect is larger for steeper relations between SFR and molecular gas. There is a general trend for all slope values to tend to ~unity for region sizes larger than 1-2 kpc, irrespective of the input SFR-cloud relation. The region size of 1-2 kpc corresponds to the area where the cloud mass function becomes fully sampled. We quantify the effects of selection biases in data tracing the SFR, either as thresholds (i.e., clouds smaller than a given mass value do not form stars) or backgrounds (e.g., diffuse emission unrelated to current star formation is counted towards the SFR). Apparently discordant observational results are brought into agreement via this simple model, and the comparison of our simulations with data for a few galaxies supports a steep (>1) power law index between SFR and molecular gas.

We have carried out survey observations toward the Galactic plane at l~38 deg in the 12CO and 13CO J=1-0 lines using the Nobeyama Radio Observatory 45-m telescope. A wide area (0.8 x 0.8 deg) was mapped with high spatial resolution (17"). The line of sight samples the gas in both the Sagittarius arm and the inter-arm regions. The present observations reveal how the structure and physical conditions vary across a spiral arm. We classify the molecular gas in the line of sight into two distinct components based on its appearance: the bright and compact B component and the fainter and diffuse (i.e., more extended) D component. The B component is predominantly seen at the spiral arm velocities, while the D component dominates at the inter-arm velocities and is also found at the spiral arm velocities. We introduce the brightness distribution function and the brightness distribution index (BDI, which indicates the dominance of the B component) in order to quantify the map's appearance. The radial velocities of BDI peaks coincide with those of high 12CO J=3-2/12CO J=1-0 intensity ratio (i.e., warm gas) and H II regions, and tend to be offset from the line brightness peaks at lower velocities (i.e., presumably downstream side of the arm). Our observations reveal that the gas structure at small scales changes across a spiral arm: bright and spatially confined structures develop in a spiral arm, leading to star formation at downstream side, while extended emission dominates in the inter-arm region.

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{\mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

We report deep Subaru Halpha observations of the XUV disk of M83. These new observations enable the first complete census of very young stellar clusters over the entire XUV disk. Combining Subaru and GALEX data with a stellar population synthesis model, we find that (1) the standard, but stochastically-sampled, initial mass function (IMF) is preferred over the truncated IMF, because there are low mass stellar clusters (10^{2-3}Msun) that host massive O-type stars; that (2) the standard Salpeter IMF and a simple aging effect explain the counts of FUV-bright and Halpha-bright clusters with masses >10^3Msun; and that (3) the Halpha to FUV flux ratio over the XUV disk supports the standard IMF. The Subaru Prime Focus Camera (Suprime-Cam) covers a large area even outside the XUV disk -- far beyond the detection limit of the HI gas. This enables us to statistically separate the stellar clusters in the disk from background contamination. The new data, model, and previous spectroscopic studies provide overall consistent results with respect to the internal dust extinction (Av~0.1 mag) and low metallicity (~0.2Zsun) using the dust extinction curve of SMC.

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $\mu$m, [O \one] 63 $\mu$m, and [N \ii] 122 $\mu$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\lesssim$ 0.95. In the warmest regions, where $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $\mu$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

We present the largest sample to date of giant molecular clouds (GMCs) in a substantial spiral galaxy other than the Milky Way. We map the distribution of molecular gas with high resolution and image fidelity within the central 5 kpc of the spiral galaxy NGC 6946 in the 12CO (J=1-0) transition. By combining observations from the Nobeyama Radio Observatory 45-meter single dish telescope and the Combined Array for Research in Millimeter Astronomy (CARMA) interferometer, we are able to obtain high image fidelity and accurate measurements of LCO compared with previous purely interferometric studies. We resolve individual giant molecular clouds (GMCs), measure their luminosities and virial masses, and derive Xco - the conversion factor from CO measurements to H2 masses - within individual clouds. On average, we find that Xco = 1.2 \times 10^20 cm-2 / (K km s-1), which is consistent within our uncertainties with previously derived Galactic values as well as the value we derive for Galactic GMCs above our mass sensitivity limit. The properties of our GMCs are largely consistent with the trends observed for molecular clouds detected in the Milky Way disk, with the exception of six clouds detected within \sim400 pc of the center of NGC 6946, which exhibit larger velocity dispersions for a given size and luminosity, as has also been observed at the Galactic center.

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

Dynamical interactions between planets and debris disks may sculpt the disk structure and impact planetary orbits, but only a few systems with both imaged planets and spatially resolved debris disks are known. With the Caltech Submm Observatory (CSO), we have observed the HR 8799 debris disk at 350{\mu}m. The 350{\mu}m map is the first spatially resolved measurement of the debris disk encircling the HR 8799 planetary system at this wavelength. Both the flux and size of the emission are consistent with a Kuiper belt of dust extending from ~100-300 AU. Although the resolution of the current map is limited, the map shows an indication of offset asymmetric emission, and several scenarios for this possibility are explored with radiative transfer calculations of a star-disk system and N-body numerical simulations of planet-disk interactions with parameters representative of the HR 8799 system.

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of H\alpha+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, H\alpha+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

We have conducted interferometric observations with CARMA and an OTF mapping with the 45-m telescope at NRO in the CO (J=1-0) emission line of NGC 3521. Combining these new data, together with CARMA+NRO45 data for M51a and archival SINGS H$\alpha$, 24$\mu$m, THINGS H I and GALEX FUV data for both galaxies, we investigate the empirical scaling law that connects the surface density of SFR and cold gas (the Schmidt-Kennicutt law) on a spatially-resolved basis. We argue that plausibly deriving SFR maps of nearby galaxies requires the diffuse stellar/dust background emission to be carefully subtracted. An approach to complete this task is presented and applied in our pixel-by-pixel analysis on both galaxies, showing that the controversial results whether the molecular S-K law is super-linear or basically linear is a result of removing or preserving the local background. In both galaxies, the power index of the molecular S-K law is 1.5-1.9 at the highest available resolution (230 pc), and decreases monotonically for decreasing resolution; while the scatter (mainly intrinsic) increases as the resolution becomes higher, indicating a trend for which the S-K law breaks down below some scale. Both quantities are systematically larger in M51a than in NGC 3521, but when plotted against the de-projected scale (\delta_{dp}), they become highly consistent between the two galaxies, tentatively suggesting that the sub-kpc molecular S-K law in spiral galaxies depends only on the considered scale, without varying amongst spiral galaxies. A logarithmic function \gamma_{H_2}=-1.1 log[\delta_{dp}/kpc]+1.4 and a linear relation \sigma_{H_2}=-0.2 [\delta_{dp}/kpc]+0.7 are obtained through fitting to the M51a data, which describes both galaxies impressively well on sub-kpc scales. A larger sample of galaxies with better quality data are required to test the general applicability of these relations.

Self-Interacting Dark Matter (SIDM) is a collisional form of cold dark matter (CDM), originally proposed to solve problems that arose when the collisionless CDM theory of structure formation was compared with observations of galaxies on small scales. The quantitative impact of the proposed elastic collisions on structure formation has been estimated previously by Monte Carlo N-body simulations and by a conducting fluid model, with apparently diverging results. To improve this situation, we make direct comparisons between new Monte Carlo N-body simulations and solutions of the conducting fluid model, for isolated SIDM haloes of fixed mass. This allows us to separate cleanly the effects of gravothermal relaxation from those of continuous mass accretion in an expanding background universe. When these two methods were previously applied to halo formation with cosmological boundary conditions, they disagreed by an order of magnitude about the size of the scattering cross section required to solve the so-called 'cusp-core problem.' We show here, however, that the methods agree with each other within 20 per cent for isolated haloes. This suggests that the two methods are consistent, and that their disagreement for cosmological haloes is not caused by a breakdown of their validity. The isolated haloes studied here undergo gravothermal collapse. We compare the solutions calculated by these two methods for gravothermal collapse starting from several initial conditions. This allows us to calibrate the heat conduction which accounts for the effect of elastic hard-sphere scattering in the fluid model. The amount of tuning of the thermal conductivity parameters required to bring the two methods into close agreement for isolated haloes, however, is too small to explain the discrepancy found previously in the cosmological context.

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 the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

We present new $^{12}$CO(J=1-0) observations of the barred galaxy NGC 4303 using the Nobeyama 45m telescope (NRO45) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The H$\alpha$ images of barred spiral galaxies often show active star formation in spiral arms, but less so in bars. We quantify the difference by measuring star formation rate and efficiency at a scale where local star formation is spatially resolved. Our CO map covers the central 2$\farcm$3 region of the galaxy; the combination of NRO45 and CARMA provides a high fidelity image, enabling accurate measurements of molecular gas surface density. We find that star formation rate and efficiency are twice as high in the spiral arms as in the bar. We discuss this difference in the context of the Kennicutt-Schimidt (KS) law, which indicates a constant star formation rate at a given gas surface density. The KS law breaks down at our native resolution ($\sim$ 250 pc), and substantial smoothing (to 500 pc) is necessary to reproduce the KS law, although with greater scatter.

We present sensitive and high angular resolution CO(1-0) data obtained by the Combined Array for Research in Millimeter-wave Astronomy (CARMA) observations toward the nearby grand-design spiral galaxy M 51. The angular resolution of 0.7" corresponds to 30 pc, which is similar to the typical size of Giant Molecular Clouds (GMCs), and the sensitivity is also high enough to detect typical GMCs. Within the 1' field of view centered on a spiral arm, a number of GMC-scale structures are detected as clumps. However, only a few clumps are found to be associated with each Giant Molecular Association (GMA), and more than 90% of the total flux is resolved out in our data. Considering the high sensitivity and resolution of our data, these results indicate that GMAs are not mere confusion of GMCs but plausibly smooth structures. In addition, we have found that the most massive clumps are located downstream of the spiral arm, which suggests that they are at a later stage of molecular cloud evolution across the arm and plausibly are cores of GMAs. By comparing with H-alpha and Pa-alpha images, most of these cores are found to have nearby star forming regions. We thus propose an evolutionary scenario for the interstellar medium, in which smaller molecular clouds collide to form smooth GMAs at spiral arm regions and then star formation is triggered in the GMA cores. Our new CO data have revealed the internal structure of GMAs at GMC scales, finding the most massive substructures on the downstream side of the arm in close association with the brightest H II regions.

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

We present some preliminary results from a series of extremely large, high-resolution N-body simulations of the formation of early nonlinear structures. We find that the high-z halo mass function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the abundance of rare halos. This discrepancy is in rough agreement with previous results based on smaller simulations. We also show that the number density of minihaloes is correlated with local matter density, albeit with a significant scatter that increases with redshift, as minihaloes become increasingly rare. The average correlation is in rough agreement with a simple analytical extended Press-Schechter model, but can differ by up to factor of 2 in some regimes.

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $\mu$m, [OIII] 88 $\mu$m, [NII] 122 $\mu$m, [CII] 158 $\mu$m and [NII] 205 $\mu$m. We observe that in the [OI] 63 $\mu$m, [OIII] 88 $\mu$m, and [NII] 122 $\mu$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $\mu$m/[CII] 158 $\mu$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $\mu$m and [OIII] 88 $\mu$m maps. Variations of the [NII] 122 $\mu$m/[NII] 205 $\mu$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

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.

We study the formation of disk-dominated galaxies in a LCDM universe. Their existence is considered to be a challenge for the LCDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gasrich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in LCDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion, and show that the dependence of the disk-dominated fraction on galaxy mass provides a further test of the merger hypothesis. [For additional details, see, Koda et al. (2009).]

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.

Spiral arm spurs are prominent features that have been observed in extinction and 8$\mu$m emission in nearby galaxies. In order to understand their molecular gas properties, we used the Owens Valley Radio Observatory to map the CO(J=1--0) emission in three spurs emanating from the inner northwestern spiral arm of M51. We report CO detections from all three spurs. The molecular gas mass and surface density are M$_{H2} \sim3\times10^6$ M$_{\sun}$ and $\Sigma_{H2} \sim$50 M$_{\sun}$ pc$^{-2}$. Thus, relative to the spiral arms, the spurs are extremely weak features. However, since the spurs are extended perpendicular to the spiral arms for $\sim$500 pc and contain adequate fuel for star formation, they may be the birthplace for observed inter-arm HII regions. This reduces the requirement for the significant time delay that would be otherwise needed if the inter-arm star formation was initiated in the spiral arms. Larger maps of galaxies at similar depth are required to further understand the formation and evolution of these spurs and their role in star formation - such data should be forthcoming with the new CARMA and future ALMA telescopes and can be compared to several recent numerical simulations that have been examining the evolution of spiral arm spurs.

We report the detection of CO(2-1) emission from the z=1.2 ultraluminous infrared galaxy (ULIRG) GOODS J123634.53+621241.3 (also known as the sub-millimeter galaxy GN26). These observations represent the first discovery of high-redshift CO emission using the new Combined Array for Research in Millimeter-Wave Astronomy (CARMA). Of all high-redshift (z>1) galaxies within the GOODS-North field, this source has the largest far-infrared (FIR) flux observed in the Spitzer 70um and 160um bands. The CO redshift confirms the optical identification of the source, and the bright CO(2-1) line suggests the presence of a large molecular gas reservoir of about 7x10^10 M(sun). The infrared-to-CO luminosity ratio of L(IR)/L'(CO) = 80+/-30 L(sun) (K Km/s pc^2)^-1 is slightly smaller than the average ratio found in local ULIRGs and high-redshift sub-millimeter galaxies. The short star-formation time scale of about 70 Myr is consistent with a starburst associated with the merger event and is much shorter than the time scales for spiral galaxies and estimates made for high-redshift galaxies selected on the basis of their B-z and z-K colors.

The merging galaxy cluster 1E 0657-56, known as the "bullet cluster," is one of the hottest clusters known. The X-ray emitting plasma exhibits bow-shock-like temperature and density jumps. The segregation of this plasma from the peaks of the mass distribution determined by gravitational lensing has been interpreted as a direct proof of collisionless dark matter. If the high shock speed inferred from the shock jump conditions equals the relative speed of the merging CDM halos, however, this merger is predicted to be such a rare event in a LCDM universe that observing it presents a possible conflict with the LCDM model. We examined this question using high resolution, 2D simulations of gas dynamics in cluster collisions to analyze the relative motion of the clusters, the bow shock, and the contact discontinuity, and relate these to the X-ray data for the bullet cluster. We find that the velocity of the fluid shock need not equal the relative velocity of the CDM components. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profiles of the colliding clusters, such a decrease of the inferred halo relative velocity would significantly increase the likelihood of finding 1E 0657-56 in a LCDM universe. (Conference proceedings based on a poster at Bash Symposium 2007)

We study the formation of disk-dominated galaxies in a Lambda CDM universe. Their existence is considered to be a challenge for the Lambda CDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gas-rich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in Lambda CDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruction of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact discontinuity and relate these to X-ray data. We focus on the "bullet cluster," 1E 0657-56, a near head-on collision of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock imply a high shock velocity 4,700 km/s. The velocity of the fluid shock has been widely interpreted as the relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the inferred halo relative velocity would increase the likelihood of finding 1E 0657-56 in a LambdaCDM universe.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

We present results from a narrow-band optical survey of a contiguous area of 1.95 deg^2, covered by the Cosmic Evolution Survey (COSMOS). Both optical narrow-band (lambda_c = 8150 AA and Delta_lambda = 120 AA) and broad-band (B, V, g', r', i', and z') imaging observations were performed with the Subaru prime-focus camera, Suprime-Cam on the Subaru Telescope. We provide the largest contiguous narrow-band survey, targetting Ly alpha emitters (LAEs) at z~5.7. We find a total of 119 LAE candidates at z~5.7. Over the wide-area covered by this survey, we find no strong evidence for large scale clustering of LAEs. We estimate a star formation rate (SFR) density of ~7*10^-4 M_sun yr^-1 Mpc^-3 for LAEs at z~5.7, and compare it with previous measurements.

We report on the discovery of a new potential galaxy threshing system in the COSMOS 2 square degree field using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru Telescope. This system consists of a giant elliptical galaxy with $M_V \approx -21.6$ and a tidally disrupted satellite galaxy with $M_V \approx -17.7$ at a photometric redshift of $z \approx 0.08$. This redshift is consistent with the spectroscopic redshift of 0.079 for the giant elliptical galaxy obtained from the Sloan Digital Sky Survey (SDSS) archive. The luminosity masses of the two galaxies are $3.7 \times 10^{12} \cal{M}_{\odot}$ and $3.1 \times 10^{9} \cal{M}_{\odot}$, respectively. The distance between the two galaxies is greater than 100 kpc. The two tidal tails emanating from the satellite galaxy extend over 150 kpc. This system would be the second well-defined galaxy threshing system found so far.

We present deep optical imaging observations of 2 square degree area, covered by the Cosmic Evolution Survey (COSMOS), made by the prime-focus Camera (Supreme-Cam) on the 8.2m Subaru Telescope. Observations were done in six broad-band [B (4459.7 AA), g' (4723.1 AA), V (5483.8 AA), r' (6213.0 AA), i' (7640.8 AA), z' (8855.0 AA)], and one narrow-band (NB816) filters. A total of 10^6 galaxies were detected to i'~26.5 mag. These data, combined with observations at u* and K-band are used to construct the photometric catalogs for the COSMOS and to measure their photometric redshifts, multi-band spectral energy distributions, stellar masses and identification of high redshift candidates. This catalog provides multi-waveband data for scientific analysis of the COSMOS survey.

We present CO interferometer observations of the barred galaxy NGC 4303 (M61). This galaxy has a gas concentration at the central region and offset ridges in the bar. Sharp velocity gradients are apparent across the ridges. Analyses of the CO data and the newborn stellar clusters revealed in HST images indicate the existence of unresolved molecular clouds with masses of 10^4-6Msun. The observed shear velocity gradient across the ridges is too small to break up giant molecular clouds. Therefore, the clouds are likely to survive passage through the ridges. We discuss a cloud orbit model in a bar potential. The model reproduces the narrow offset ridges and sharp velocity gradients across the ridges in NGC 4303. We discuss cloud-cloud collisions (and close interactions) as a possible triggering mechanism for star formation. The newborn stellar clusters in NGC 4303 are located predominantly at the leading sides of the offset ridges, where cloud orbits are densely populated and suggest a high collisional frequency and possibly a high rate of triggered star formation. Cloud-based dynamics is less dissipative than smooth hydrodynamic models, possibly extending the timescales of gas dynamical evolution and gas fueling to central regions in barred galaxies.

We present high-resolution ($1\farcs8$ -- $4\farcs5$) CO data of the Virgo spiral galaxy NGC 4569, obtained using the Nobeyama Millimeter Array. We found that the molecular gas is highly concentrated in the circumnuclear region with two off-center peaks. A CO image with the highest angular resolution of $2\farcs0 \times 1\farcs8$ shows that six blobs likely form a circumnuclear elliptical ring (CER) with a semi-major axis radius of 720 pc. The CER shows a strongly twisted velocity field, and the position--velocity (PV) diagram shows significant forbidden velocity components. These kinetic features are understood as being non-circular motion due to a bar-potential. We found that the CER coincides with the H$\alpha$ bright central core and that the mass ratio of the molecular gas to the dynamical mass is about 18%. These results support a gaseous inflow scenario induced by a weak bar potential and self-gravity of the gas.

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations assuming the standard cold dark matter model ($\Omega_0 = 1$). Our analysis on formation and evolution of the galaxy and its substructures archives until $z = 2.0$. In such a high-redshift universe, the Einstein-de Sitter universe is still a good approximation for a recently favored $\Lambda$-dominated universe, and then our results does not depend on the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider radiative cooling below $T = 10^4 K$, we adopt a fine mass resolution ($m_{\rm SPH} = 1.12 \times 10^3 \Msun$). Because of the cooling, each clump immediately forms a `core-halo' structure which consists of a baryonic core and a dark matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically decreases each pericenter passage. At $z = 2$, more than half of the clumps become baryon dominated systems (baryon mass/total mass $> 0.5$). Our results support the tidal evolution scenario of the formation of globular clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe.

New 13CO data from the BU-FCRAO Milky Way Galactic Ring Survey (GRS) are analyzed to understand the shape and internal motions of molecular clouds. For a sample of more than five hundred molecular clouds, we find that they are preferentially elongated along the Galactic plane. On the other hand, their spin axes are randomly oriented. We therefore conclude that the elongation is not supported by internal spin but by internal velocity anisotropy. It has been known that some driving mechanisms are necessary to sustain the supersonic velocity dispersion within molecular clouds. The mechanism for generating the velocity dispersion must also account for the preferred elongation. This excludes some driving mechanisms, such as stellar winds and supernovae, because they do not produce the systemic elongation along the Galactic plane. Driving energy is more likely to come from large scale motions, such as the Galactic rotation.

The Cosmic Evolution Survey (COSMOS) is a Hubble Space Telescope (HST) treasury project.The COSMOS aims to perform a 2 square degree imaging survey of an equatorial field in $I$(F814W) band, using the Advanced Camera for Surveys (ACS). Such a wide field survey, combined with ground-based photometric and spectroscopic data, is essential to understand the interplay between large scale structure, evolution and formation of galaxies and dark matter. In 2004, we have obtained high-quality, broad band images of the COSMOS field ($B, V, r^\prime, i^\prime,$ and $ z^\prime$) using Suprime-Cam on the Subaru Telescope, and we have started our new optical multi-band program, COSMOS-21 in 2005. Here, we present a brief summary of the current status of the COSMOS project together with contributions from the Subaru Telescope. Our future Subaru program, COSMOS-21, is also discussed briefly.

We have initiated the first CO interferometer survey of early-type spiral galaxies (S0-Sab). We observed five early-type spiral galaxies with HII nuclei (indicating circumnuclear starburst activities). These observations indicate gas masses for the central kiloparsec of \sim 1-5% of the dynamical masses. Such low gas mass fractions suggest that large-scale gravitational instability in the gas is unlikely to be the driving cause for the starburst activities. The Toomre Q values were >1 (mostly >3) within the central kiloparsec, indicating that the gas disks are globally gravitationally stable. The area filling factor of the gas disks is estimated to be about 0.05. This small value indicates the existence of lumpy structure, i.e. molecular clouds, in the globally-gravitationally stable disks. The typical surface density of the molecular clouds is as high as \sim 3000 Msun pc^{-2}. We reconsider the nature of the Toomre Q criterion, and conclude that the Q derived from CO observations indicates neither star formation nor molecular cloud formation. This argument should be valid not only for the circumnuclear disks but also for any region in galactic disks. We tentatively explore an alternative model, i.e. cloud-cloud collisions, as an initiating mechanism of star formation.

CO emission has been imaged at 2" resolution in the central 10 kpc of the barred spiral galaxy NGC 3367. This galaxy has bipolar synchrotron lobes out to a radii of 6 kpc, straddling the compact nucleus. The peak molecular emission is in a source of radius 2" (425 pc) centered on the galaxy nucleus. The molecular mass is about 3 times 10**8 Msun in this peak an about 5.9 times 10**8 Msun within a radius of 4".5 (950 pc). The very large gas masses in the central source imply extinctions sufficiently high to completely obscure optical emission lines (e.g.BLR) associated with the nuclear radio source. The observed Balmer lines probably originate in the narrow line region few hundred pc from the nucleus. The CO emission in the central region is elongated NE-SW, very similar to the PA of the large scale synchrotron lobes. This elongation is likely due to non-axisymmetric gravitational potential of the stellar bar. We infer that the NE radio continuum lobe is on the far side of the galaxy and the SW lobe is on the near side. The central mass of molecular gas is sufficient mass to power the AGN accretion luminosity for over 10**8 yrs at 3 Msun per yr.

We report on high-resolution interferometer observations of the 12CO(J=1-0) emission in the central 5 kpc region of the Seyfert 2 galaxy NGC 4501. The observations were made using the Nobeyama Millimeter Array during a long-term CO line survey of Virgo spirals. The major features are: (1) a nuclear concentration with a radius of r~5'' (390 pc), which is resolved into double peaks, and (2) spiral arms which extend out from the nuclear region. The nuclear component has a mass of 1.3x10^8 Mo, which corresponds to ~3.5% of the dynamical mass, and shows a slight non-circular motion. The double peaks are separated by ~4''.7. (370 pc), and located on the root of optical spiral arms in a HST image. The gas arms are associated with the spiral dust lanes, and are linked to the central double peaks. The non-circular motions along the molecular arms indicate the fact that the gas is driven by the density wave, rather than the stochastic processes. We calculated the gas cloud orbits in a stellar spiral potential, and explained the observed CO spiral arms and non-circular motions. We suggest that the central gas condensation arises from spiral-driven gas transfer. We estimated and compared the effect of two possible mechanisms of angular-momentum transfer: galactic shock, and gravitational torques. We discuss that the galactic shock is dominant.

We found that loosely wound spiral shocks in an isothermal gas disk caused by a non-axisymmetric potential are hydrodynamically unstable, if the shocks are strong enough. High resolution, global hydrodynamical simulations using three different numerical schemes, i.e. AUSM, CIP, and SPH, show similarly that trailing spiral shocks with the pitch angle of larger than ~10 deg wiggle, and clumps are developed in the shock-compressed layer. The numerical simulations also show clear wave crests that are associated with ripples of the spiral shocks. The spiral shocks tend to be more unstable in a rigidly rotating disk than in a flat rotation. This instability could be an origin of the secondary structures of spiral arms, i.e. the spurs/fins, observed in spiral galaxies. In spite of this local instability, the global spiral morphology of the gas is maintained over many rotational periods. The Kelvin-Helmholtz (K-H) instability in a shear layer behind the shock is a possible mechanism for the wiggle instability. The Richardson criterion for the K-H stability is expressed as a function of the Mach number, the pitch angle, and strength of the background spiral potential. The criterion suggests that spiral shocks with smaller pitch angles and smaller Mach numbers would be more stable, and this is consistent with the numerical results.

We describe a new method to accelerate neighbor searches on GRAPE, i.e. a special purpose hardware that efficiently calculates gravitational forces and potentials in $N$-body simulations. In addition to the gravitational calculations, GRAPE simultaneously constructs the lists of neighbor particles that are necessary for Smoothed Particle Hydrodynamics (SPH). However, data transfer of the neighbor lists from GRAPE to the host computer is time consuming, and can be a bottleneck. In fact, the data transfer can take about the same time as the calculations of forces themselves. Making use of GRAPE's special treatment of neighbor lists, we can reduce the amount of data transfer if we search neighbors in the order that the neighbor lists, constructed in a single GRAPE run, overlap each other. We find that the Morton-ordering requires very low additional calculation and programming costs, and results in successful speed-up on data transfer. We show some benchmark results in the case of GRAPE-5. Typical reduction in transferred data becomes as much as 90%. This method is suitable not only for GRAPE-5, but also GRAPE-3 and the other versions of GRAPE.

Based on a high-resolution CO survey of Virgo spirals with the Nobeyama Millimeter-wave Array, we determined the dynamical centers using velocity fields, and derived position-velocity diagrams (PVDs) along the major axes of the galaxies across their dynamical centers. We applied a new iteration method to derive rotation curves (RCs), which reproduce the observed PVDs. The obtained high-accuracy RCs generally show steep rise in the central 100 to 200 pc regions, followed by flat rotation in the disk. We applied a deconvolution method to calculate the surface-mass density (SMD) using the RCs based on two extreme assumptions that the mass distribution is either spherical or thin-disk shaped. Both assumptions give nearly identical results, agreeing with each other within a factor of two at any radii. The SMD distributions revealed central massive cores with peak SMD of 10^4 - 10^5 Msun pc^-2 and total mass within 200 pc radius of the order of about 10^9 Msun Correlation analysis among the derived parameters show that the central CO-line intensity is positively correlated with the central SMD, which suggests that the deeper is the gravitational potential, the higher is the molecular gas concentration in the nuclei regardless morphological types.

We present high-angular-resolution (1".5 - 5") interferometer observations of the \co emission in the central region of the SA(s)c galaxy NGC 4254. The observations were obtained using the Nobeyama Millimeter-wave Array (NMA) during the course of a long-term CO line survey of Virgo spirals. We present the spectra, channel maps, integrated intensity distributions, velocity fields, position--velocity diagrams, and compare the data with various optical images. The rotation velocity is already finite at the nucleus, or at least it rises steeply to 80 km/s within the central 1", indicating the existence of a massive core of 10^8 Msun within 1" (80 pc) radius. The CO intensity maps show that the inner disk has well-developed multiple spiral arms, winding out from a bar-shaped elongated molecular complex. In addition to the bisymmetric spiral arms, an asymmetric tightly wound arm with high molecular gas density is found to wind out from the molecular bar. The molecular spiral arms, particularly the tightly wound arm, well traces optical dark lanes, and are associated with H$\alpha$ arms having many H {\sc ii} regions. The inner asymmetric spiral structures can be explained by ram-pressure distortion of inter-arm low density regions of the inner disk by the intra-cluster gas wind, and is indeed reproduced by a hydrodynamical simulation.

We present the results of the Virgo high-resolution CO survey (ViCS) obtained with the Nobeyama Millimeter-wave Array (NMA). This survey was made in the course of a long-term project at Nobeyama from 1999 December through 2002 April. The objects were selected from Virgo cluster members, considering CO richness from single dish flux, mild inclination, and lack of strong tidal perturbations. The central 1 arc min regions (4.7 kpc) of 15 spiral galaxies were observed with resolutions of 2 to 5 arcsec and 10 to 20 km/s, and sensitivities of 20 mJy/ beam for a 10 km/s channel. The objects lie at the same distance of the Virgo cluster (16.1 Mpc), which is advantageous for comparisons among individual galaxies. We describe the details of observations and data reduction, and present an atlas of integrated CO intensity maps, velocity fields and position-velocity diagrams along the major axes. The molecular gas morphology in the Virgo galaxies shows a wealth of variety, not specifically depending on the Hubble types.

within 24 lines with 80 characters) We present ^12CO (1-0) observations in the central 4.5 kpc (1 arcmin) of the Halpha/Radio lobe galaxy NGC 3079 with the Nobeyama Millimeter Array. The molecular gas shows four components: a main disk, spiral arms, a nuclear disk, and a nuclear core. The main disk extends beyond our spatial coverage. The spiral arms are superimposed on the main disk. The nuclear disk with about 600 pc radius has an intense concentration of molecular gas. Its appearance on PV diagrams is indicative of oval motions of the gas, rather than circular. The nuclear core is more compact than our resolution. Though it is unresolved, the nuclear core shows a very high velocity about 200 km/s even at the radius of about 100 pc on the PV diagram. We propose a model that NGC 3079 contains a weak bar. This model explains the observed features: the main disk and spiral arms result from gaseous x1-orbits and associated crowding respectively. The nuclear disk arises from gaseous x2-orbits. From the appearance of the spiral arms on the PV diagram, the pattern speed of the bar is estimated to be 55+-10 km/s/kpc. The high velocity of the nuclear core cannot be explained by our model for a bar. Thus we attribute it to a central massive core with a dynamical mass of 10^9 Msun within the central 100 pc. This mass is three orders of magnitude more massive than that of a central black hole in this galaxy.

Rotation curves are often used to estimate the mass distribution of spiral galaxies, assuming circular rotation of disks. However, non-circular motions caused by a non-axisymmetric gravitational potential, such as a stellar bar, may disturb the velocity field, resulting in errors in mass estimation, especially in the central regions of galaxies, because the line-of-sight velocity depends on the viewing angles in a non-axisymmetric flow. Observing rotation curves of edge-on galaxies in time-dependent numerical simulations from different viewing angles, we obtain errors in the estimation of galactic mass from the rotation curves. In the most extreme case, the ellipticity of gas orbits is as high as 0.8 in the central regions, even if the bar potential is weak. When rotation curves are defined as the highest velocity envelope of position-velocity diagrams, the mass estimated from the rotation curves is larger than the true mass by a factor of five for 15% of the viewing angles, and the ratio between the apparent mass and true mass is less than six for any viewing angle. The overestimation in mass occurs more frequently than the underestimation.

The structure of the interstellar medium in the central kpc region of a galaxy with a weak bar-like potential is investigated taking into account realistic cooling and heating processes and the self-gravity of the gas. Using high resolution hydrodynamical simulations, it is revealed that the resonant structures (e.g. smooth spiral shocks and a nuclear ring) are very different from those seen in past numerical models where simple models of the ISM, i.e. non-self-gravitating, isothermal gas were assumed. We find that the pc-scale filaments and clumps form large scale spirals, which resemble those seen in real galaxies. The fine structures are different between the arms and in the nuclear region. The next generation millimeter interferometer (ALMA) may reveal the fine structures of the cold gas in nearby galaxies. We also find a large scale anisotropy in the gas temperature, which is caused due to non-circular velocity field of the gas.The damped orbit model based on the epicyclic approximation explains the distribution of the hot (> 10^4 K) and cold (< 100 K) gases appearing alternately around the galactic center. Because of the temperature anisotropy, cold gases observed by molecular lines do not necessarily represent the real gas distribution in galaxies. Position-Velocity diagrams depend strongly on the viewing angles. As a result, the rotational velocity inferred from the PV maps could be two times larger or smaller than the true circular velocity.

We have performed high-resolution synthesis observations of the 12CO (J=1-0) line emission from the radio lobe edge-on spiral NGC 3079 using a 7-element mm-wave interferometer at Nobeyama, which consisted of the 45-m telescope and 6-element Array. The molecular nuclear disk (NMD) of 750 pc radius is found to be inclined by 20 degrees from the optical disk, and the NMD has spiral arms. An ultra-high-density molecular core (UHC) was found at the nucleus. The gaseous mass of the UHC within 125 pc radius is as large as about 3 x 10^8 Msun, an order of magnitude more massive than that in the same area of the Galactic Center, and the mean density is as high as approximately 3 x 10^3 H2 cm^-3. A position-velocity diagram along the major axis indicates that the rotation curve starts at a finite velocity exceeding 300 \kms already from the nucleus. The surface mass density in the central region is estimated to be as high as about 10^5 Msun pc^-2, producing a very deep gravitational potential. We argue that the very large differential rotation in such a deep potential will keep the UHC gravitationally stable from current star formation.

We investigate the origin of a unified scaling relation in spiral galaxies. Observed spiral galaxies are spread on a plane in the three-dimensionallogarithmic space of luminosity L, radius R and rotation velocity V. The plane is expressed as $L \propto (V R)^{\alpha}$ in I-passband, where $\alpha$ is a constant. On the plane, observed galaxies are distributed in an elongated region which looks like the shape of a surfboard. The well-known scaling relations, L-V (Tully-Fisher relation), V-R (also the Tully-Fisher relation) and R-L (Freeman's law), can be understood as oblique projections of the surfboard-like plane into 2-D spaces. This unified interpretation of the known scaling relations should be a clue to understand the physical origin of all the relations consistently. Furthermore, this interpretation can also explain why previous studies could not find any correlation between TF residuals and radius. In order to clarify the origin of this plane, we simulate formation and evolution of spiral galaxies with the N-body/SPH method, including cooling, star formation and stellar feedback. Initial conditions are set to isolated 14 spheres with two free parameters, such as mass and angular momentum. The CDM (h=0.5, $\Omega_0=1$) cosmology is considered as a test case. The simulations provide the following two conclusions: (a) The slope of the plane is well reproduced but the zero-point is not. This zero-point discrepancy could be solved in a low density ($\Omega_0<1$) and high expansion (h>0.5) cosmology. (b) The surfboard-shaped plane can be explained by the control of galactic mass and angular momentum.

We discuss the origin of the Tully-Fisher (TF) relation using the $N$-body/SPH method, which includes cooling, star formation and stellar feedback of energy, mass and metals. We consider initially rotating overdense spheres, and trace formation processes of disk galaxies from $z=25$ to $z=0$ in the Cold Dark Matter (CDM) cosmology. To clarify the origin of the TF relation, we simulate formation of 14 galaxies with different masses and spin parameters, and compute observable values, such as the total magnitude and the line-width. We find that the simulated galaxies reproduce the slope and scatter of the TF relation: the slope is originated in the difference of total galactic masses, and the scatter is produced by the difference of initial spin parameters. As well as the TF relation, observed features of spiral galaxies, such as the exponential light-profile and the flat rotation curve, are reproduced in our simulations, which were assumed {\it a priori} in past semi-analytical approaches.

We present high-resolution central-to-outer rotation curves for Sb, SBb, Sc, and SBc galaxies. We discuss their general characteristics, particularly their central behavior, as well as dependencies on morphological types, activity, and peculiarity. The rotation curves generally show a steep nuclear rise, and high-velocity central rotation, followed by a broad maximum in the disk and then by a flat rotation due to the massive halo. Since the central high velocity and steep rise are common to all massive galaxies, it cannot be due to non-circular motions. Disk rotation curves of barred galaxies show larger dispersion than those of normal galaxies which is probably due to non-circular motions. Interacting galaxies show often perturbed outer rotation curves, while their central rotation shows no particular peculiarity. Also, central activities, such as starbursts and AGN, appear to show no particular correlation with the property of rotation curves. This would suggest that the central activities are triggered by a more local effect than the global dynamical property.