We present the third phase of the largest high-frequency, high-resolution imaging survey of 231 nearby, hard X-ray selected active galactic nuclei (AGNs), with a very high 98% ± 1% detection fraction. This survey presents VLA 22 GHz radio observations with 1 ″ spatial resolution covering over 6 orders of magnitude in radio luminosity in nearby AGNs that span ∼4 orders of magnitude in black hole mass and X-ray luminosity. We identify three different radio morphologies: 44% ± 3% (102/231) are compact or unresolved, 46% ± 3% (106/231) show an extended structure (star formation, possible one-sided jets, etc.), and 8% ± 2% (19/231) have a biconical or two-sided jet-like morphology. The remaining 2% ± 1% (4/231) sources are nondetections. The radio-to-X-ray luminosity ratios of the Swift-BAT AGNs in our sample ( L _R / L _14−195keV ∼ 10 ^−5.5 and L _R / L _2−10keV ∼ 10 ^−5 with a scatter of ∼0.5 dex) are similar to that of coronally active stars ( L _R / L _X ∼ 10 ^−5 ). For most targets, extended emission in radio-quiet objects is broadly consistent with the expectation for star formation from previous far-infrared observations once the contribution from the radio core has been subtracted. Our sample represents nearby analogs of distant AGNs at the peak of black hole growth, and thus the high detection fraction in our work has important implications for future high frequency AGN radio surveys with the next generation Very Large Array or Square Kilometre Array, both of which should detect large fractions of more distant AGNs.

The EDGE-CALIFA survey provides spatially resolved optical integral-field unit and CO spectroscopy for 125 galaxies selected from the Calar Alto Legacy Integral Field Area Survey (CALIFA) Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different images. By joining these pixel tables to lower-dimensional tables that provide radial profiles, integrated spectra, or global properties, it is possible to investigate the dependence of local conditions on large-scale properties. The database is freely accessible and has been utilized in several publications. We illustrate the use of this database and highlight the effects of CO upper limits on the inferred slopes of the local scaling relations between the stellar mass, star formation rate (SFR), and H _2 surface densities. We find that the correlation between H _2 and SFR surface density is the tightest among the three relations.

We present a $^{12}$CO($J$=2-1) survey of 60 local galaxies using data from the Atacama Large Millimeter/submillimeter Compact Array as part of the Extragalactic Database for Galaxy Evolution: the ACA EDGE survey. These galaxies all have integral field spectroscopy from the CALIFA survey. Compared to other local galaxy surveys, ACA EDGE is designed to mitigate selection effects based on CO brightness and morphological type. Of the 60 galaxies in ACA EDGE, 36 are on the star-formation main sequence, 13 are on the red sequence, and 11 lie in the ``green valley" transition between these sequences. We test how star formation quenching processes affect the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}=$SFR/$M_{\rm mol}$, and related quantities in galaxies with stellar masses $10\leq$log[$M_\star/$M$_\odot$]$\leq11.5$ covering the full range of morphological types. We observe a systematic decrease of the molecular-to-stellar mass fraction ($R^{\rm mol}_{\star}$) with decreasing level of star formation activity, with green valley galaxies having also lower SFE$_{\rm mol}$ than galaxies on the main sequence. On average, we find that the spatially resolved SFE$_{\rm mol}$ within the bulge region of green valley galaxies is lower than in the bulges of main sequence galaxies if we adopt a constant CO-to-H$_2$ conversion factor, $\alpha_{\rm CO}$. While efficiencies in main sequence galaxies remain almost constant with galactocentric radius, in green valley galaxies we note a systematic increase of SFE$_{\rm mol}$, $R^{\rm mol}_{\star}$, and specific star formation rate, sSFR, with increasing radius. Our results suggest that although gas depletion (or removal) seems to be the most important driver of the star-formation quenching in galaxies transiting through the green valley, a reduction in star formation efficiency is also required during this stage.

We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately (∼1 kpc) spatially resolved measurements of CO (1–0) and CO (2–1) intensity, I CO, and mid-IR intensity, I MIR, at 8, 12, 22, and 24 μm. The I CO versus I MIR relationship is reasonably described by a power law with slopes 0.7–1.2 and normalization I CO ∼ 1 K km s−1 at I MIR ∼ 1 MJy sr−1. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO (1–0) and CO (2–1) allow us to infer that R21∝IMIR0.2 , in good agreement with other work. The 8 μm and 12 μm bands, with strong polycyclic aromatic hydrocarbon (PAH) features, show steeper CO versus mid-IR slopes than the 22 and 24 μm, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass (M ⋆) and anticorrelates with star formation rate/M ⋆. At ∼1 kpc resolution, the first four PHANGS–JWST targets show CO-to-mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities.

We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio R12/13 and the properties of the stars and ionized gas. Higher R12/13 values are found in interacting galaxies than in non-interacting galaxies. The global R12/13 slightly increases with infrared color F60/F100, but appears insensitive to other host galaxy properties such as morphology, stellar mass, or galaxy size. We also present annulus-averaged R12/13 profiles for our sample up to a galactocentric radius of 0.4r25 ( 6 kpc), taking into account the 13CO(1-0) non-detections by spectral stacking. The radial profiles of R12/13 are quite flat across our sample. Within galactocentric distances of 0.2r25, azimuthally-averaged R12/13 increases with star formation rate. However, the Spearman rank correlation tests show the azimuthally-averaged R12/13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25 from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, R12/13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on R12/13, which further complicate the interpretations of R12/13 variations.

Shocks and torques produced by non-axisymmetric structures such as spiral arms and bars may transport gas to galaxy central regions. We test this hypothesis by studying the dependence of the concentration of CO luminosity (CCO) and molecular gas (Cmol) and the star formation rate (CSFR) in the central ∼2 kpc on the strength of non-axisymmetric disk structure using a sample of 57 disk galaxies selected from the EDGE-CALIFA survey. TheCmolis calculated using a CO-to-H2conversion factor that decreases with higher metallicity and higher stellar surface density. We find thatCmolis systematically 0.22 dex lower thanCCO. We confirm that highCmoland strong non-axisymmetric disk structure are more common in barred galaxies than in unbarred galaxies. However, we find that spiral arms also increaseCmol. We show that there is a good correlation betweenCmoland the strength of non-axisymmetric structure (which can be due to a bar, spiral arms, or both). This suggests that the stronger the bars and spirals, the more efficient the galaxy is at transporting cold gas to its center. Despite the small subsample size, theCmolof the four Seyferts are not significantly reduced compared to inactive galaxies of similar disk structure, implying that the active galactic nucleus feedback in Seyferts may not notably affect the molecular gas distribution in the central ∼2 kpc. We find thatCSFRtightly correlates withCmolin both unbarred and barred galaxies. Likewise, elevatedCSFRis found in galaxies with strong disk structure. Our results suggest that the disk structure, either spirals or bars, can transport gas to the central regions, with higher inflow rates corresponding to stronger structure, and consequently boost central star formation. Both spirals and bars play, therefore, an essential role in the secular evolution of disk galaxies.

We measure the molecular-to-atomic gas ratio, R mol , and the star formation rate (SFR) per unit molecular gas mass, SFE mol , in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We stack ALMA 12 CO ( J = 2−1) spectra coherently using H i velocities from the VIVA survey to detect faint CO emission out to galactocentric radii r gal ∼ 1.2 r 25 . We determine the scale lengths for the molecular and stellar components, finding a ∼3:5 relation compared to ∼1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute R mol as a function of different physical quantities. While the spatially resolved R mol on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius R e , R mol ( r < R e ), shows a systematic increase with the level of H i , truncation and/or asymmetry (H I perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within R e , R ⋆ mol ( r < R e ) and R ⋆ atom ( r < R e ) , shows that VERTICO galaxies have increasingly lower R ⋆ atom ( r < R e ) for larger levels of H I perturbation (compared to field galaxies matched in stellar mass), but no significant change in R ⋆ m o l ( r < R e ) . We also measure a clear systematic decrease of the SFE mol within R e , SFE mol ( r < R e ), with increasingly perturbed H i . Therefore, compared to field galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their R mol and decreases the efficiency with which their molecular gas forms stars.

Feedback from an active galactic nucleus (AGN) is often implicated as a mechanism that leads to the quenching of galactic star formation. However, AGN-driven quenching is challenging to reconcile with observations that AGN hosts tend to harbour equal (or even excess) amounts of gas compared with inactive galaxies of similar stellar mass. In this paper, we investigate whether AGN feedback happens on sub-galactic (kpc) scales, an effect that might be difficult to detect with global gas measurements. Using kpc-scale measurements of molecular gas (Sigma_H2) and stellar mass (Sigma_*) surface densities taken from the EDGE-CALIFA survey, we show that the gas fractions of central AGN regions are typically a factor of ~2 lower than in star-forming regions. Based on four galaxies with the best spaxel statistics, the difference between AGN and star-forming gas fractions is seen even within a given galaxy, indicating that AGN feedback is able to deplete the molecular gas reservoir in the central few kpc.

We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map 12CO (2–1), 13CO (2–1), and C18O (2–1) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array. The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, 12CO (2–1) observations, and data reduction procedures. We characterize global 12CO (2–1) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two nondetections, and produce resolved 12CO (2–1) data products (median resolution = 8″ ≈ 640 pc). Azimuthally averaged 12CO (2–1) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size–mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size–mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy data set for studying galaxy evolution in our closest cluster.

It is not clear whether gamma-ray bursts (GRBs) are unbiased tracers of cosmic star formation at z > 2. Since dusty starburst galaxies are significant contributors to the cosmic star formation at these redshifts, they should form a major part of the GRB host population. However, recent studies at z ≤ 2 have shown that the majority of the star formation activity in GRB hosts is not obscured by dust. Here, we investigate the galaxy-scale dust obscuration in z ∼ 2–3.5 GRB hosts preselected to have high-resolution, high signal-to-noise afterglow spectra in the rest-frame ultraviolet (UV) and thus relatively low line-of-sight dust obscuration. We present new deep Very Large Array observations of four GRB hosts, and compare the radio-based (upper limits on the) total star formation rates (SFRs) to the unobscured SFRs derived from fits to the optical-UV spectral energy distribution. The fraction of the total SFR that is obscured by dust in these galaxies is found to be <90% in general, and ≲50% for GRB 021004 in particular. These observations suggest that z ∼ 2–3.5 GRBs with UV-unobscured sightlines originate in star-forming galaxies with low overall dust obscuration, unlike the dusty starburst population.

We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.

We present multiwavelength modeling of the afterglow from the long γ -ray burst (GRB) 160625B using Markov Chain Monte Carlo techniques of the afterglowpy Python package. GRB 160625B is an extremely bright burst with a rich set of observations spanning from radio to γ -ray frequencies. These observations range from ∼0.1 days to >1000 days, thus making this event extremely well suited to such modeling. In this work we compare top-hat and Gaussian jet structure types in order to find best-fit values for the GRB jet collimation angle, viewing angle, and other physical parameters. We find that a Gaussian-shaped jet is preferred (2.7 σ –5.3 σ ) over the traditional top-hat model. Our estimate for the opening angle of the burst ranges from 1.°26 to 3.°90, depending on jet-shape model. We also discuss the implications that assumptions on jet shape, viewing angle, and particularly the participation a fraction of electrons have on the final estimation of GRB intrinsic energy release and the resulting energy budget of the relativistic outflow. Most notably, allowing the participation fraction to vary results in an estimated total relativistic energy of ∼10 53 erg. This is two orders of magnitude higher than when the total fraction is assumed to be unity; thus, this parameter has strong relevance for placing constraints on long GRB central engines, details of the circumburst media, and host environment.

We present post-jet-break Hubble Space Telescope, Very Large Array, and Chandra observations of the afterglow of the long gamma-ray bursts GRB 160625B (between 69 and 209 days) and GRB 160509A (between 35 and 80 days). We calculate the post-jet-break decline rates of the light curves and find the afterglow of GRB 160625B is inconsistent with a simple t(-3/4) steepening over the break, expected from the geometric effect of the jet edge entering our line of sight. However, the favored optical post-break decline is also inconsistent with the f(nu) proportional to t(-p) decline (where p 2.3 from the pre-break light curve), which is expected from exponential lateral expansion of the jet; perhaps suggesting lateral expansion that only affects a fraction of the jet. The post-break decline of GRB 160509A is consistent with both the t(-3/4) steepening and with f(nu) proportional to t(-p). We also use boxfit to fit afterglow models to both light curves and find both to be energetically consistent with a millisecond magnetar central engine, but the magnetar parameters need to be extreme (i.e., E similar to 3 x 10(52) erg). Finally, the late-time radio light curves of both afterglows are not reproduced well by boxfit and are inconsistent with predictions from the standard jet model; instead, both are well represented by a single power-law decline (roughly f(nu) proportional to t(-1)) with no breaks. This requires a highly chromatic jet break and possibly a two-component jet for both bursts.

Understanding how galaxies cease to form stars represents an outstanding challenge for galaxy evolution theories. This process of “star formation quenching” has been related to various causes, including active galactic nuclei activity, the influence of large-scale dynamics, and the environment in which galaxies live. In this paper, we present the first results from a follow-up of CALIFA survey galaxies with observations of molecular gas obtained with the APEX telescope. Together with the EDGE-CARMA observations, we collected 12CO observations that cover approximately one effective radius in 472 CALIFA galaxies. We observe that the deficit of galaxy star formation with respect to the star formation main sequence (SFMS) increases with the absence of molecular gas and with a reduced efficiency of conversion of molecular gas into stars, which is in line with the results of other integrated studies. However, by dividing the sample into galaxies dominated by star formation and galaxies quenched in their centres (as indicated by the average value of the Hα equivalent width), we find that this deficit increases sharply once a certain level of gas consumption is reached, indicating that different mechanisms drive separation from the SFMS in star-forming and quenched galaxies. Our results indicate that differences in the amount of molecular gas at a fixed stellar mass are the primary drivers for the dispersion in the SFMS, and the most likely explanation for the start of star formation quenching. However, once a galaxy is quenched, changes in star formation efficiency drive how much a retired galaxy differs in its star formation rate from star-forming ones of similar masses. In other words, once a paucity of molecular gas has significantly reduced star formation, changes in the star formation efficiency are what drives a galaxy deeper into the red cloud, hence retiring it.

We present a multilinear analysis to determine the significant predictors of star formation in galaxies using the combined EDGE-CALIFA sample of galaxies. We analyze 1845 kpc-scale lines of sight across 39 galaxies with molecular line emission measurements from EDGE combined with optical IFU data drawn from CALIFA. We use the Least Absolute Shrinkage and Selection Operator (LASSO) to identify significant factors in predicting star formation rates. We find that the local star formation rate surface density is increased by higher molecular gas surface densities and stellar surface densities. In contrast, we see lower star formation rates in systems with older stellar populations, higher gas- and stellar-phase metallicities and larger galaxy masses. We also find a significant increase in star formation rate with galactocentric radius normalized by the disk scale length, which suggests additional parameters regulating star formation rate not explored in this study.

The purpose of this white paper is to summarize the current landscape of astrophotonic devices and their scientific impact, highlight the key issues, and outline specific technological and organizational approaches to address these issues in the coming decade and thereby enable new discoveries as we embark on the era of extremely large telescopes.

Pradip Gatkine1*, Sylvain Veilleux1,2,3, John Mather4, Christopher Betters5, Jonathan Bland-Hawthorn5, Julia Bryant5, S. Bradley Cenko4,2, Mario Dagenais6, Drake Deming1, Simon Ellis7, Matthew Greenhouse4, Andrew Harris1, Nemanja Jovanovic8, Steve Kuhlmann9, Alexander Kutyrev4, Sergio Leon-Saval5, Kalaga Madhav10, Samuel Moseley4, Barnaby Norris5, Bernard Rauscher4, Martin Roth10, and Stuart Vogel1

Astrophotonics is the application of versatile photonic technologies to channel, manipulate, and disperse guided light from one or more telescopes to achieve scientific objectives in astronomy in an efficient and cost-effective way. The developments and demands from the telecommunication industry have driven a major boost in photonic technology and vice versa in the last 40 years. The photonic platform of guided light in fibers and waveguides has opened the doors to next-generation instrumentation for both ground- and space-based telescopes in optical and near/mid-IR bands, particularly for the upcoming extremely large telescopes (ELTs). The large telescopes are pushing the limits of adaptive optics to reach close to a near-diffraction-limited performance. The photonic devices are ideally suited for capturing this AO-corrected light and enabling new and exciting science such as characterizing exoplanet atmospheres. The purpose of this white paper is to summarize the current landscape of astrophotonic devices and their scientific impact, highlight the key issues, and outline specific technological and organizational approaches to address these issues in the coming decade and thereby enable new discoveries as we embark on the era of extremely large telescopes.

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