We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $\mu$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy\ radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $\alpha=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 \sigma$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

The mass-luminosity relation for late-type stars has long been a critical tool for estimating stellar masses. However, there is growing need for both a higher-precision relation and a better understanding of systematic effects (e.g., metallicity). Here we present an empirical relationship between Mks and mass spanning $0.075M_\odot<M<0.70M_\odot$. The relation is derived from 62 nearby binaries, whose orbits we determine using a combination of Keck/NIRC2 imaging, archival adaptive optics data, and literature astrometry. From their orbital parameters, we determine the total mass of each system, with a precision better than 1% in the best cases. We use these total masses, in combination with resolved Ks magnitudes and system parallaxes, to calibrate the mass-Mks relation. The result can be used to determine masses of single stars with a precision of 2-3%, which we confirm by a comparison to dynamical masses from the literature. The precision is limited by scatter around the best-fit relation beyond mass uncertainties, perhaps driven by intrinsic variation in the mass-Mks relation or underestimated measurement errors. We find the effect of [Fe/H] on the mass-Mks relation is likely negligible for metallicities in the Solar neighborhood (0.0+/-2.2% change in mass per dex change in [Fe/H]). This weak effect is consistent with predictions from the Dartmouth Stellar Evolution Database, but inconsistent with those from MESA Isochrones and Stellar Tracks. A sample of binaries with a wider range of abundances will be required to discern the importance of metallicity in extreme populations (e.g., in the Galactic Halo or thick disk).

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.

Observations of the baryon to dark matter fraction in galaxies through cosmic time are a fundamental test for galaxy formation models. Recent observational studies have suggested that some disk galaxies at $z>1$ host declining rotation curves, in contrast with observations of low redshift disk galaxies where stellar or HI rotation curves flatten at large radii. We present an observational counterexample, a galaxy named DSFG850.95 at $z=1.555$ (4.1 Gyr after the big bang) that hosts a flat rotation curve between radii of $\sim$6--14 kpc (1.2--2.8 disk scale lengths) and has a dark matter fraction of $0.44\pm0.08$ at the H-band half light radius, similar to the Milky Way. We create position-velocity and position-dispersion diagrams using Keck/MOSFIRE spectroscopic observations of H$\alpha$ and [NII] emission features, which reveal a flat rotation velocity of $V_{\rm flat}=285\pm12$ km/s and an ionized gas velocity dispersion of $\sigma_{0}=48\pm4$ km/s. This galaxy has a rotation-dominated velocity field with $V_{\rm flat}/\sigma_{0}\sim6$. Ground-based H-band imaging reveals a disk with S\'ersic index of $1.29\pm0.03$, an edge-on inclination angle of $87\pm2^{\circ}$, and an H-band half light radius of $8.4\pm0.1$ kpc. Our results point to DSFG850.95 being a massive, rotationally-supported disk galaxy with a high dark-matter-to-baryon fraction in the outer galaxy, similar to disk galaxies at low redshift.

The next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of $\sim10^{5}\,M_\odot$. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations of a $z\sim4.5$ gas rich star-forming galaxy. Using the DESPOTIC radiative transfer code that encompasses simultaneous thermal, chemical, and statistical equilibrium in calculating the molecular and atomic level transitions of CO from ALMA for comparison. We find that observations of CO(1-0) are especially important for tracing the systemic redshift of the galaxy and the total mass of the well-shielded molecular gas reservoir, while even CO(2-1) can predominantly trace denser gas regions distinct from CO(1-0). The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies of CO(1-0) routine at $z\sim2-5$.

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe.

One important result from recent large integral field spectrograph (IFS) surveys is that the intrinsic velocity dispersion of galaxies traced by star-forming gas increases with redshift. Massive, rotation-dominated discs are already in place at z~2, but they are dynamically hotter than spiral galaxies in the local Universe. Although several plausible mechanisms for this elevated velocity dispersion (e.g. star formation feedback, elevated gas supply, or more frequent galaxy interactions) have been proposed, the fundamental driver of the velocity dispersion enhancement at high redshift remains unclear. We investigate the origin of this kinematic evolution using a suite of cosmological simulations from the FIRE (Feedback In Realistic Environments) project. Although IFS surveys generally cover a wider range of stellar masses than in these simulations, the simulated galaxies show trends between intrinsic velocity dispersion, SFR, and redshift in agreement with observations. In both the observed and simulated galaxies, intrinsic velocity dispersion is positively correlated with SFR. Intrinsic velocity dispersion increases with redshift out to z~1 and then flattens beyond that. In the FIRE simulations, intrinsic velocity dispersion can vary significantly on timescales of <100 Myr. These variations closely mirror the time evolution of the SFR and gas inflow rate. By cross-correlating pairs of intrinsic velocity dispersion, gas inflow rate, and SFR, we show that increased gas inflow leads to subsequent enhanced star formation, and enhancements in intrinsic velocity dispersion tend to temporally coincide with increases in gas inflow rate and SFR.

Though half of cosmic starlight is absorbed by dust and reradiated at long wavelengths (3$\mu$m-3mm), constraints on the infrared through millimeter galaxy luminosity function (the `IRLF') are poor in comparison to the rest-frame ultraviolet and optical galaxy luminosity function, particularly at z>2.5. Here we present a backward evolution model for interpreting number counts, redshift distributions, and cross-band flux density correlations in the infrared and submillimeter sky, from 70$\mu$m-2mm, using a model for the IRLF out to the epoch of reionization. Mock submillimeter maps are generated by injecting sources according to the prescribed IRLF and flux densities drawn from model spectral energy distributions that mirror the distribution of SEDs observed in $0<z<5$ dusty star-forming galaxies (DSFGs). We explore two extreme hypothetical case-studies: a dust-poor early Universe model, where DSFGs contribute negligibly ($<$10%) to the integrated star-formation rate density at $z>4$, and an alternate dust-rich early Universe model, where DSFGs dominate $\sim$90% of $z>4$ star-formation. We find that current submm/mm datasets do not clearly rule out either of these extreme models. We suggest that future surveys at 2mm will be crucial to measuring the IRLF beyond $z\sim4$. The model framework developed in this paper serves as a unique tool for the interpretation of multiwavelength IR/submm extragalactic datasets and will enable more refined constraints on the IRLF than can be made from direct measurements of individual galaxies' integrated dust emission.

We present a deep survey of the SuperCLASS super-cluster - a region of sky known to contain five Abell clusters at redshift $z\sim0.2$ - performed using the Arcminute Microkelvin Imager (AMI) Large Array (LA) at 15.5$~$GHz. Our survey covers an area of approximately 0.9 square degrees. We achieve a nominal sensitivity of $32.0~\mu$Jy beam$^{-1}$ toward the field centre, finding 80 sources above a $5\sigma$ threshold. We derive the radio colour-colour distribution for sources common to three surveys that cover the field and identify three sources with strongly curved spectra - a high-frequency-peaked source and two GHz-peaked-spectrum sources. The differential source count (i) agrees well with previous deep radio source count, (ii) exhibits no evidence of an emerging population of star-forming galaxies, down to a limit of 0.24$~$mJy, and (iii) disagrees with some models of the 15$~$GHz source population. However, our source count is in agreement with recent work that provides an analytical correction to the source count from the SKADS Simulated Sky, supporting the suggestion that this discrepancy is caused by an abundance of flat-spectrum galaxy cores as-yet not included in source population models.

We present in this work the development of a solar data assimilation method based on an axisymmetric mean field dynamo model and magnetic surface data, our mid-term goal is to predict the solar quasi cyclic activity. Here we focus on the ability of our algorithm to constrain the deep meridional circulation of the Sun based on solar magnetic observations. To that end, we develop a variational data assimilation technique. Within a given assimilation window, the assimilation procedure minimizes the differences between data and the forecast from the model, by finding an optimal meridional circulation in the convection zone, and an optimal initial magnetic field, via a quasi-Newton algorithm. We demonstrate the capability of the technique to estimate the meridional flow by a closed-loop experiment involving 40 years of synthetic, solar-like data. By assimilating the synthetic magnetic proxies annually, we are able to reconstruct a (stochastic) time-varying meridional circulation which is also slightly equatorially asymmetric. We show that the method is robust in estimating a flow whose level of fluctuation can reach 30% about the average, and that the horizon of predictive capability of the method is of the order of 1 cycle length.

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow which best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires the knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique, for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations, by varying their spatial density, temporal density, noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step towards using this technique to i) better constrain the physical processes occurring inside the Sun and ii) better predict solar activity on decadal time scales.

We present near-infrared and optical spectroscopic observations of a sample of 450$\mu$m and 850$\mu$m-selected dusty star-forming galaxies (DSFGs) identified in a 400 arcmin$^2$ area in the COSMOS field. Thirty-one sources of the 102 targets were spectroscopically confirmed at $0.2<z<4$, identified primarily in the near-infrared with Keck MOSFIRE and some in the optical with Keck LRIS and DEIMOS. The low rate of confirmation is attributable both to high rest-frame optical obscuration in our targets and limited sensitivity to certain redshift ranges. The high-quality photometric redshifts available in the COSMOS field allow us to test the robustness of photometric redshifts for DSFGs. We find a subset (11/31$\approx35$%) of DSFGs with inaccurate ($\Delta z/(1+z)>0.2$) or non-existent photometric redshifts; these have very distinct spectral energy distributions from the remaining DSFGs, suggesting a decoupling of highly obscured and unobscured components. We present a composite rest-frame 4300--7300\AA\ spectrum for DSFGs, and find evidence of 200$\pm$30 km s$^{-1}$ gas outflows. Nebular line emission for a sub-sample of our detections indicate that hard ionizing radiation fields are ubiquitous in high-z DSFGs, even more so than typical mass or UV-selected high-z galaxies. We also confirm the extreme level of dust obscuration in DSFGs, measuring very high Balmer decrements, and very high ratios of IR to UV and IR to H$\alpha$ luminosities. This work demonstrates the need to broaden the use of wide bandwidth technology in the millimeter to the spectroscopic confirmations of large samples of high-z DSFGs, as the difficulty in confirming such sources at optical/near-infrared wavelengths is exceedingly challenging given their obscuration.

We present the results of 325 MHz GMRT observations of a super-cluster field, known to contain five Abell clusters at redshift $z \sim 0.2$. We achieve a nominal sensitivity of $34\,\mu$Jy beam$^{-1}$ toward the phase centre. We compile a catalogue of 3257 sources with flux densities in the range $183\,\mu\rm{Jy}\,-\,1.5\,\rm{Jy}$ within the entire $\sim 6.5$ square degree field of view. Subsequently, we use available survey data at other frequencies to derive the spectral index distribution for a sub-sample of these sources, recovering two distinct populations -- a dominant population which exhibit spectral index trends typical of steep-spectrum synchrotron emission, and a smaller population of sources with typically flat or rising spectra. We identify a number of sources with ultra-steep spectra or rising spectra for further analysis, finding two candidate high-redshift radio galaxies and three gigahertz-peaked-spectrum radio sources. Finally, we derive the Euclidean-normalised differential source counts using the catalogue compiled in this work, for sources with flux densities in excess of $223 \, \mu$Jy. Our differential source counts are consistent with both previous observations at this frequency and models of the low-frequency source population. These represent the deepest source counts yet derived at 325 MHz. Our source counts exhibit the well-known flattening at mJy flux densities, consistent with an emerging population of star-forming galaxies; we also find marginal evidence of a downturn at flux densities below $308 \, \mu$Jy, a feature so far only seen at 1.4 GHz.

The most prodigious starburst galaxies are absent in massive galaxy clusters today, but their connection with large scale environments is less clear at $z\gtrsim2$. We present a search of large scale structure around a galaxy cluster core at $z=2.095$ using a set of spectroscopically confirmed galaxies. We find that both color-selected star-forming galaxies (SFGs) and dusty star-forming galaxies (DSFGs) show significant overdensities around the $z=2.095$ cluster. A total of 8 DSFGs (including 3 X-ray luminous active galactic nuclei, AGNs) and 34 SFGs are found within a 10 arcmin radius (corresponds to $\sim$15 cMpc at $z\sim2.1$) from the cluster center and within a redshift range of $\Delta z=0.02$, which leads to galaxy overdensities of $\delta_{\rm DSFG}\sim12.3$ and $\delta_{\rm SFG}\sim2.8$. The cluster core and the extended DSFG- and SFG-rich structure together demonstrate an active cluster formation phase, in which the cluster is accreting a significant amount of material from large scale structure while the more mature core may begin to virialize. Our finding of this DSFG-rich structure, along with a number of other protoclusters with excess DSFGs and AGNs found to date, suggest that the overdensities of these rare sources indeed trace significant mass overdensities. However, it remains puzzling how these intense star formers are triggered concurrently. Although an increased probability of galaxy interactions and/or enhanced gas supply can trigger the excess of DSFGs, our stacking analysis based on 850 $\mu$m images and morphological analysis based on rest-frame optical imaging do not show such enhancements of merger fraction and gas content in this structure.

(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z.

We present a survey toward 16 Low Luminosity Objects (LLOs with an internal luminosity, Lint, lower than 0.2 Lsun) with N2H+ (1-0), N2H+ (3-2), N2D+ (3-2), HCO+ (3-2) and HCN (3-2) using the Arizona Radio Observatory Kitt Peak 12m Telescope and Submillimeter Telescope. Our goal is to probe the nature of these faint protostars which are believed to be either very low mass or extremely young protostars. We find that the N2D+/N2H+ column density ratios of LLOs are similar to those of typical starless cores and Class 0 objects. The N2D+/N2H+ column density ratios are relatively high (> 0.05) for LLOs with kinetic temperatures less than 10 K in our sample. The distribution of N2H+ (1-0) line widths spreads between that of starless cores and young Class 0 objects. If we use the line width as a dynamic evolutionary indicator, LLOs are likely young Class 0 protostellar sources. We further use the optically thick tracers, HCO+ (3-2) and HCN (3-2), to probe the infall signatures of our targets. We derive the asymmetry parameters from both lines and estimate the infall velocities by fitting the HCO+ (3-2) spectra with two-layer models. As a result, we identify eight infall candidates based on the infall velocities and seven candidates have infall signatures supported by asymmetry parameters from at least one of HCO+ (3-2) and HCN (3-2).

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $\mu$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

While performing a survey to detect recoiling supermassive black holes, we have identified an unusual source having a projected offset of 800 pc from a nearby dwarf galaxy. The object, SDSS J113323.97+550415.8, exhibits broad emission lines and strong variability. While originally classified as a supernova (SN) because of its nondetection in 2005, we detect it in recent and past observations over 63 yr and find over a magnitude of rebrightening in the last 2 years. Using high-resolution adaptive optics observations, we constrain the source emission region to be <12 pc and find a disturbed host-galaxy morphology indicative of recent merger activity. Observations taken over more than a decade show narrow [O III] lines, constant ultraviolet emission, broad Balmer lines, a constant putative black hole mass over a decade of observations despite changes in the continuum, and optical emission-line diagnostics consistent with an active galactic nucleus (AGN). However, the optical spectra exhibit blueshifted absorption, and eventually narrow Fe II and [Ca II] emission, each of which is rarely found in AGN spectra. While this peculiar source displays many of the observational properties expected of a potential black hole recoil candidate, some of the properties could also be explained by a luminous blue variable star (LBV) erupting for decades since 1950, followed by a Type IIn SN in 2001. Interpreted as an LBV followed by a SN analogous to SN 2009ip, the multi-decade LBV eruptions would be the longest ever observed, and the broad Halpha emission would be the most luminous ever observed at late times (>10 yr), larger than that of unusually luminous supernovae such as SN 1988Z, suggesting one of the most extreme episodes of pre-SN mass loss ever discovered.

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

Sakharov oscillations, conventionally discussed in the context of early universe evolution and the anisotropy of cosmic microwave background radiation, is the manifestation of interfering acoustic waves synchronously generated in an ideal fluid. Here we report the laboratory demonstration of Sakharov oscillations in a quenched atomic superfluid. We quench the sample by Feshbach tuning and monitor the subsequent density fluctuations at different time and length scales by in situ imaging. Sakharov oscillations are identified as the multi-peak structure in the atomic density power spectrum, resembling that of the cosmic microwave background. We also observe Sakharov oscillations in the time domain, from which we extract the energy dispersion of the superfluid, and determine the sonic horizon of the excitations.

We present the results of a search for galaxy alignments in 12 galaxy clusters at z>0.5, a statistically complete subset of the very X-ray luminous clusters from the MAssive Cluster Survey (MACS). Using high-quality images taken with the Hubble Space Telescope (HST) that render measurement errors negligible, we find no radial galaxy alignments within 500 kpc of the cluster centres for a sample of 545 spectroscopically confirmed cluster members. A mild, but statistically insignificant trend favouring radial alignments is observed within a radius of 200 kpc and traced to galaxies on the cluster red sequence. Our results for massive clusters at z>0.5 are in stark contrast to the findings of previous studies which find highly significant radial alignments of galaxies in nearby clusters at z~0.1 out to at least half the virial radius using imaging data from the SDSS. The discrepancy becomes even more startling if radial alignment becomes more prevalent at decreasing clustercentric distance, as suggested by both our and previous work. We investigate and discuss potential causes for the disparity between our findings based on HST images of clusters at z>0.5 and those obtained using groundbased images of systems at z~0.1. We conclude that the most likely explanation is either dramatic evolution with redshift (in the sense that radial alignments are less pronounced in dynamically younger systems) or the presence of systematic biases in the analysis of SDSS imaging data that cause at least partly spurious alignment signals.

We present a near-infrared extinction study of nine dense cores at evolutionary stages between starless to Class I. Our results show that the density structure of all but one observed cores can be modeled with a single power law rho \propto r^p between ~ 0.2R-R of the cores. The starless cores in our sample show two different types of density structures, one follows p ~ -1.0 and the other follows p ~ -2.5, while the protostellar cores all have p ~ -2.5. The similarity between the prestellar cores with p ~ -2.5 and protostellar cores implies that those prestellar cores could be evolving towards the protostellar stage. The slope of p ~ -2.5 is steeper than that of an singular isothermal sphere, which may be interpreted with the evolutionary model of cores with finite mass.

We examine the possibility that Ultraluminous X-ray sources (ULXs) represent the extreme end of the black hole X-ray binary (XRB) population. Based on their X-ray properties, we suggest that ULXs are persistently in a high/hard spectral state and we propose a new disk-jet model that can accomodate both a high accretion rate and a hard X-ray spectrum. Our model predicts that the modified disk emission can be substantially softer than that predicted by a standard disk as a result of jet cooling and this may explain the unusually soft components that are sometimes present in the spectra of bright ULXs. We also show that relativistic beaming of jet emission can indeed account for the high X-ray luminosities of ULXs, but strong beaming produces hard X-ray spectra that are inconsistent with observations. We predict the beamed synchrotron radio emission should have a flat spectrum with a flux density < 0.01 mJy.