The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

Producing images from interferometer data requires accurate modeling of the sources in the field of view, which is typically done using the CLEAN algorithm. Given the large number of degrees of freedom in interferometeric images, one constrains the possible model solutions for CLEAN by masking regions that contain emission. Traditionally this process has largely been done by hand. This approach is not possible with today's large data volumes which require automated imaging pipelines. This paper describes an automated masking algorithm that operates within CLEAN called AUTO-MULTITHRESH. This algorithm was developed and validated using a set of ~1000 ALMA images chosen to span a range of intrinsic morphology and data characteristics. It takes a top-down approach to producing masks: it uses the residual images to identify significant peaks and then expands the mask to include emission associated with these peaks down to lower signal-to-noise noise. The AUTO-MULTITHRESH algorithm has been implemented in CASA and has been used in production as part of the ALMA Imaging Pipeline starting with Cycle 5. It has been shown to be able to mask a wide range of emission ranging from simple point sources to complex extended emission with minimal tuning of the parameters based on the point spread function of the data. Although the algorithm was developed for ALMA, it is general enough to have been used successfully with data from other interferometers with appropriate parameter tuning. Integrating the algorithm more deeply within the minor cycle could lead to future performance improvements.

We present detailed submillimeter- through centimeter-wave observations of the extraordinary extragalactic transient AT2018cow. The apparent characteristics -- the high radio luminosity, the long-lived emission plateau at millimeter bands, and the sub-relativistic velocity -- have no precedent. A basic interpretation of the data suggests $E_k \gtrsim 10^{48}$ erg coupled to a fast but sub-relativistic ($v \approx 0.13c$) shock in a dense ($n_e \approx 3 \times 10^{5}\,$cm$^{-3}$) medium. We find that the X-ray emission is not naturally explained by an extension of the radio-submm synchrotron spectrum, nor by inverse Compton scattering of the dominant blackbody UVOIR photons by energetic electrons within the forward shock. By $\Delta t \approx20\,$days, the X-ray emission shows spectral softening and erratic inter-day variability. Taken together, we are led to invoke an additional source of X-ray emission: the central engine of the event. Regardless of the nature of this central engine, this source heralds a new class of energetic transients shocking a dense medium, which at early times are most readily observed at millimeter wavelengths.

We present results from a joint ALMA-HST study of the nearby spiral galaxy NGC 628. We combine the HST LEGUS database of over 1000 stellar clusters in NGC 628 with ALMA Cycle 4 millimeter/submillimeter observations of the cold dust continuum that span ~15 square kpc including the nuclear region and western portions of the galaxy's disk. The resolution -- 1.1" or approximately 50 pc at the distance of NGC 628 -- allows us to constrain the spatial variations in the slope of the millimeter dust continuum as a function of the ages and masses of the nearby stellar clusters. Our results indicate an excess of dust emission in the millimeter assuming a typical cold dust model for a normal star-forming galaxy, but little correlation of the dust continuum slope with stellar cluster age or mass. For the depth and spatial coverage of these observations, we cannot substantiate the millimeter/submillimeter excess arising from the processing of dust grains by the local interstellar radiation field. We detect a bright unknown source in NGC 628 in ALMA bands 4 and 7 with no counterparts at other wavelengths from ancillary data. We speculate this is possibly a dust obscured supernova.

We report the results of an investigation to determine the nature of the offset active galactic nucleus (AGN) found in the source CXO J101527.2+625911. Hubble Space Telescope and Chandra X-ray observatory data had suggested that the offset AGN, which has an angular separation of only 0\farcs26 from the center of the host galaxy, is a recoiled Super Massive Black Hole (rSMBH). We carried out high angular resolution observations with both the VLBA (1.54 GHz) and the VLA (10.0 GHz \& 33.0 GHz) and detected a single compact radio source in the center of the host galaxy, with no radio continuum emission associated with the offset AGN. The detected radio source has a high brightness temperature value of $T_b=7.2\times10^{7}$ K, indicating that the radio emission is associated with an AGN. Furthermore, we present the decomposition of high-resolution KECK spectra of the [O III]5007\AA line into two narrow emission line components, which is a characteristic sign of a dual black hole system. These new radio and optical wavelength results suggest that CXO J101527.2+625911 is the host of a dual black hole system rather than a rSMBH.

We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$\sigma$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langle\alpha\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems.

We analyze Chandra X-ray observatory data for a sample of 63 luminous infrared galaxies (LIRGs), sampling the lower-infrared luminosity range of the Great Observatories All-Sky LIRG survey (GOALS), which includes the most luminous infrared selected galaxies in the local universe. X-rays are detected for 84 individual galaxies within the 63 systems, for which arcsecond resolution X-ray images, fluxes, infrared and X-ray luminosities, spectra and radial profiles are presented. Using X-ray and MIR selection criteria, we find AGN in (31$\pm$5)% of the galaxy sample, compared to the (38$\pm$6)% previously found for GOALS galaxies with higher infrared luminosities (C-GOALS I). Using mid-infrared data, we find that (59$\pm$9)% of the X-ray selected AGN in the full C-GOALS sample do not contribute significantly to the bolometric luminosity of the host galaxy. Dual AGN are detected in two systems, implying a dual AGN fraction in systems that contain at least one AGN of (29$\pm$14)%, compared to the (11$\pm$10)% found for the C-GOALS I sample. Through analysis of radial profiles, we derive that most sources, and almost all AGN, in the sample are compact, with half of the soft X-ray emission generated within the inner $\sim 1$ kpc. For most galaxies, the soft X-ray sizes of the sources are comparable to those of the MIR emission. We also find that the hard X-ray faintness previously reported for the bright C-GOALS I sources is also observed in the brightest LIRGs within the sample, with $L_{\rm FIR}>8\times10^{10}$ L$_{\odot}$.

Post-starburst or "E+A" galaxies are rapidly transitioning from star-forming to quiescence. While the current star formation rate of post-starbursts is already at the level of early type galaxies, we recently discovered that many have large CO-traced molecular gas reservoirs consistent with normal star forming galaxies. These observations raise the question of why these galaxies have such low star formation rates. Here we present an ALMA search for the denser gas traced by HCN (1--0) and HCO+ (1--0) in two CO-luminous, quiescent post-starburst galaxies. Intriguingly, we fail to detect either molecule. The upper limits are consistent with the low star formation rates and with early-type galaxies. The HCN/CO luminosity ratio upper limits are low compared to star-forming and even many early type galaxies. This implied low dense gas mass fraction explains the low star formation rates relative to the CO-traced molecular gas and suggests the state of the gas in post-starburst galaxies is unusual, with some mechanism inhibiting its collapse to denser states. We conclude that post-starbursts galaxies are now quiescent because little dense gas is available, in contrast to the significant CO-traced lower density gas reservoirs that still remain.

The spectral type of Mrk 1018 changed from Type 1.9 to 1 and returned back to 1.9 over a period of 40 years. We have investigated physical mechanisms responsible for the spectral change in Mrk 1018 by analyzing archival spectral and imaging data. Two kinematically distinct broad-line components, blueshifted and redshifted components, are found from spectral decomposition. The velocity offset curve of the broad-line as a function of time shows a characteristic pattern. An oscillating recoiled supermassive black hole (rSMBH) scenario is proposed to explain the observed velocity offset in broad emission lines. A Bayesian Markov-Chain Monte Carlo simulation is performed to derive the best fit orbital parameters, we find that the rSMBH has a highly eccentric orbit with a period of $\sim$29 years. The AGN activity traced by variation of broad H$\beta$ emission line is found to increase and decrease rapidly at the start and end of the cycle, and peaked twice at the start and near the end of the cycle. Extinction at the start and end of the cycle (when its spectral type is Type 1.9) is found to increase due to increased covering factor. Perturbations of the accretion disk caused by pericentric passage can reasonably explain the AGN activity and spectral change in Mrk 1018. Currently, the spectral type of Mrk 1018 is Type 1.9, and we do not know if it will repeat a similar pattern of spectral change in the future, but if it does, then spectral type will turn to Type 1 around mid 2020's.

We combine information from the clustering of HI galaxies in the 100% data release of the Arecibo Legacy Fast ALFA survey (ALFALFA), and from the HI content of optically-selected galaxy groups found in the Sloan Digital Sky Survey (SDSS) to constrain the relation between halo mass $M_h$ and its average total HI mass content $M_{\rm HI}$. We model the abundance and clustering of neutral hydrogen through a halo-model-based approach, parametrizing the $M_{\rm HI}(M_h)$ relation as a power law with an exponential mass cutoff. To break the degeneracy between the amplitude and low-mass cutoff of the $M_{\rm HI}(M_h)$ relation, we also include a recent measurement of the cosmic HI abundance from the $\alpha$.100 sample. We find that all datasets are consistent with a power-law index $\alpha=0.44\pm 0.08$ and a cutoff halo mass $\log_{10}M_{\rm min}/(h^{-1}M_\odot)=11.27^{+0.24}_{-0.30}$. We compare these results with predictions from state-of-the-art magneto-hydrodynamical simulations, and find both to be in good qualitative agreement, although the data favours a significantly larger cutoff mass that is consistent with the higher cosmic HI abundance found in simulations. Both data and simulations seem to predict a similar value for the HI bias ($b_{\rm HI}=0.875\pm0.022$) and shot-noise power ($P_{\rm SN}=92^{+20}_{-18}\,[h^{-1}{\rm Mpc}]^3$) at redshift $z=0$.

The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data.

We have carried out a systematic search for recoiling supermassive black holes (rSMBH) using the Chandra Source and SDSS Cross Matched Catalog. From the survey, we have detected a potential rSMBH, 'CXO J101527.2+625911' at z=0.3504. The CXO J101527.2+625911 has a spatially offset (1.26$\pm$0.05 kpc) active SMBH and kinematically offset broad emission lines (175$\pm$25 km s$^{\rm -1}$ relative to systemic velocity). The observed spatial and velocity offsets suggest this galaxy could be a rSMBH, but we also have considered a possibility of dual SMBH scenario. The column density towards the galaxy center was found to be Compton thin, but no X-ray source was detected. The non-detection of the X-ray source in the nucleus suggests either there is no obscured actively accreting SMBH, or there exists an SMBH but has a low accretion rate (i.e. low-luminosity AGN (LLAGN)). The possibility of the LLAGN was investigated and found to be unlikely based on the H$\alpha$ luminosity, radio power, and kinematic arguments. This, along with the null detection of X-ray source in the nucleus supports our hypothesis that the CXO J101527.2+625911 is a rSMBH. Our GALFIT analysis shows the host galaxy to be a bulge-dominated elliptical. The weak morphological disturbance and small spatial and velocity offsets suggest that CXO J101527.2+625911 could be in the final stage of merging process and about to turn into a normal elliptical galaxy.

We study the limits of the spatial and velocity resolution of radio interferometry to infer the mass of supermassive black holes (SMBHs) in galactic centres using the kinematics of circum-nuclear molecular gas, by considering the shapes of the galaxy surface brightness profile, signal-to-noise ratios (S/Ns) of the position-velocity diagram (PVD) and systematic errors due to the spatial and velocity structure of the molecular gas. We argue that for fixed galaxy stellar mass and SMBH mass, the spatial and velocity scale that need to be resolved increase and decrease, respectively, with decreasing \sersic\ index of the galaxy surface brightness profile. We validate our arguments using simulated PVDs for varying beam size and velocity channel width. Furthermore, we consider the systematic effects on the inference of the SMBH mass by simulating PVDs including the spatial and velocity structure of the molecular gas, which demonstrates that their impacts are not significant for a PVD with good S/N unless the spatial and velocity scale associated with the systematic effects are comparable to or larger than the angular resolution and velocity channel width of the PVD from pure circular motion. Also, we caution that a bias in a galaxy surface brightness profile owing to the poor resolution of a galaxy photometric image can largely bias the SMBH mass by an order of magnitude. This study shows the promise and the limit of ALMA observations for measuring SMBH mass using molecular gas kinematics and provides a useful technical justification for an ALMA proposal with the science goal of measuring SMBH mass.

We combine data from the Sloan Digital Sky Survey (SDSS) and the Arecibo Legacy Fast ALFA Survey (ALFALFA) to study the cold atomic gas content of galaxies in groups and clusters in local universe. A careful cross-matching of galaxies in the SDSS, ALFALFA and SDSS group catalogs provides a sample of group galaxies with stellar masses $10^{8.4} M_{\odot} \le M_{*} \le 10^{10.6} M_{\odot}$ and group halo masses $10^{12.5} h^{-1} M_{\odot} \le M_h \le 10^{15.0} h^{-1} M_{\odot}$. Controlling our sample in stellar mass and redshift, we find no significant radial variation in the galaxy \hi\ gas-to-stellar mass ratio for the halo mass range in our sample. However, the fraction of galaxies detected in ALFALFA declines steadily towards the centers of groups with the effect being most prominent in the most massive halos. In the outskirts of massive halos a hint of a depressed detection fraction for low mass galaxies suggests pre-processing that decreases the \hi\ in these galaxies before they fall into massive clusters. We interpret the decline in the ALFALFA detection of galaxies in the context of a threshold halo mass for ram pressure stripping for a given galaxy stellar mass. The lack of an observable decrease in the galaxy \hi\ gas-to-stellar mass ratio with the position of galaxies within groups and clusters highlights the difficulty of detecting the impact of environment on the galaxy \hi\ content in a shallow \hi\ survey.

Weinberg (2012) described a constructive algorithm for computing the marginal likelihood, Z, from a Markov chain simulation of the posterior distribution. Its key point is: the choice of an integration subdomain that eliminates subvolumes with poor sampling owing to low tail-values of posterior probability. Conversely, this same idea may be used to choose the subdomain that optimizes the accuracy of Z. Here, we explore using the simulated distribution to define a small region of high posterior probability, followed by a numerical integration of the sample in the selected region using the volume tessellation algorithm described in Weinberg (2012). Even more promising is the resampling of this small region followed by a naive Monte Carlo integration. The new enhanced algorithm is computationally trivial and leads to a dramatic improvement in accuracy. For example, this application of the new algorithm to a four-component mixture with random locations in 16 dimensions yields accurate evaluation of Z with 5% errors. This enables Bayes-factor model selection for real-world problems that have been infeasible with previous methods.

Isothermal and self-gravitating systems bound by non-conducting and conducting walls are known to be unstable if the density contrast between the center and the boundary exceeds critical values. We investigate the equilibrium and dynamical evolution of isothermal and self-gravitating system embedded in potential well, which can be the situation of many astrophysical objects such as the central parts of the galaxies, or clusters of galaxies with potential dominated by dark matter, but is still limited to the case where the potential well is fixed during the evolution. As the ratio between the depth of surrounding potential well and potential of embedded system becomes large, the potential well becomes effectively the same boundary condition as conducting wall, which behaves like a thermal heat bath. We also use the direct N-body simulation code, NBODY6 to simulate the dynamical evolution of stellar system embedded in potential wells and propose the equilibrium models for this system. In deep potential well, which is analogous to the heat bath with high temperature, the embedded self-gravitating system is dynamically hot, and loosely bound or can be unbound since the kinetic energy increases due to the heating by the potential well. On the other hand, the system undergoes core collapse by self-gravity when potential well is shallow. Binary heating can stop the collapse and leads to the expansion, but the evolution is very slow because the potential as a heat bath can absorb the energy generated by the binaries. The system can be regarded as quasi-static. Density and velocity dispersion profiles from the N-body simulations in the final quasi-equilibrium state are similar to our equilibrium models assumed to be in thermal equilibrium with the potential well.

We introduce a new galaxy image decomposition tool, GALPHAT (GALaxy PHotometric ATtributes), to provide full posterior probability distributions and reliable confidence intervals for all model parameters. GALPHAT is designed to yield a high speed and accurate likelihood computation, using grid interpolation and Fourier rotation. We benchmark this approach using an ensemble of simulated Sersic model galaxies over a wide range of observational conditions: the signal-to-noise ratio S/N, the ratio of galaxy size to the PSF and the image size, and errors in the assumed PSF; and a range of structural parameters: the half-light radius $r_e$ and the Sersic index $n$. We characterise the strength of parameter covariance in Sersic model, which increases with S/N and $n$, and the results strongly motivate the need for the full posterior probability distribution in galaxy morphology analyses and later inferences. The test results for simulated galaxies successfully demonstrate that, with a careful choice of Markov chain Monte Carlo algorithms and fast model image generation, GALPHAT is a powerful analysis tool for reliably inferring morphological parameters from a large ensemble of galaxies over a wide range of different observational conditions. (abridged)

We introduce a novel image decomposition package, GALPHAT, that provides robust estimates of galaxy surface brightness profiles using Bayesian Markov Chain Monte Carlo. The GALPHAT-determined posterior distribution of parameters enables us to assign rigorous statistical confidence intervals to maximum a posteriori estimates and to test complex galaxy formation and evolution hypotheses. We describe the GALPHAT algorithm, assess its performance using test image data, and demonstrate that it has sufficient speed for production analysis of a large galaxy sample. Finally we briefly introduce our ongoing science program to study the distribution of galaxy structural properties in the local universe using GALPHAT.

We have carried out N-body simulations for rotating star clusters with equal mass and compared the results with Fokker-Planck models. These two different approaches are found to produce fairly similar results, although there are some differences with regard to the detailed aspects. We confirmed the acceleration of the core collapse of a cluster due to an initial non-zero angular momentum and found a similar evolutionary trend in the central density and velocity dispersion in both simulations. The degree of acceleration depends on the initial angular momentum. Angular momentum is being lost from the cluster due to the evaporation of stars with a large angular momentum on a relaxation time scale.