Search

Using a combination of HST, JWST, and ALMA data, we perform spatially resolved spectral energy distributions (SED) fitting of fourteen 4

Most star formation in the local Universe occurs in spiral galaxies, but their origin remains an unanswered question. Various theories have been proposed to explain the development of spiral arms, each predicting different spatial distributions of the interstellar medium. This study maps the star formation rate (SFR) and gas-phase metallicity of nine spiral galaxies with the TYPHOON survey to test two dominating theories: density wave theory and dynamic spiral theory. We discuss the environmental effects on our galaxies, considering reported environments and merging events. Taking advantage of the large field of view covering the entire optical disk, we quantify the fluctuation of SFR and metallicity relative to the azimuthal distance from the spiral arms. We find higher SFR and metallicity in the trailing edge of NGC~1365 (by 0.117~dex and 0.068~dex, respectively) and NGC~1566 (by 0.119~dex and 0.037~dex, respectively), which is in line with density wave theory. NGC~2442 shows a different result with higher metallicity (0.093~dex) in the leading edge, possibly attributed to an ongoing merging. The other six spiral galaxies show no statistically significant offset in SFR or metallicity, consistent with dynamic spiral theory. We also compare the behaviour of metallicity inside and outside the co-rotation radius (CR) of NGC~1365 and NGC~1566. We find comparable metallicity fluctuations near and beyond the CR of NGC~1365, indicating gravitational perturbation. NGC~1566 shows the greatest fluctuation near the CR, in line with the analytic spiral arms. Our work highlights that a combination of mechanisms explains the origin of spiral features in the local Universe., Comment: 22 pages, 10 figures, 4 tables, accepted for publication in MNRAS

We measure the ionised gas velocity dispersions of star-forming galaxies in the MAGPI survey ($z\sim0.3$) and compare them with galaxies in the SAMI ($z\sim0.05$) and KROSS ($z\sim1$) surveys to investigate how the ionised gas velocity dispersion evolves. For the first time, we use a consistent method that forward models galaxy kinematics from $z=0$ to $z=1$. This method accounts for spatial substructure in emission line flux and beam smearing. We investigate the correlation between gas velocity dispersion and galaxy properties to understand the mechanisms that drive gas turbulence. We find that in both MAGPI and SAMI galaxies, the gas velocity dispersion more strongly correlates with the star-formation rate surface density ($\Sigma_{\rm SFR}$) than with a variety of other physical properties, and the average gas velocity dispersion is similar, at the same $\Sigma_{\rm SFR}$, for SAMI, MAGPI and KROSS galaxies. The results indicate that mechanisms related to $\Sigma_{\rm SFR}$ could be the dominant driver of gas turbulence from $z\sim1$ to $z\sim0$, for example, stellar feedback and/or gravitational instability. The gas velocity dispersion of MAGPI galaxies is also correlated with the non-rotational motion of the gas, illustrating that in addition to star-formation feedback, gas transportation and accretion may also contribute to the gas velocity dispersion for galaxies at $z\sim 0.3$. KROSS galaxies only have a moderate correlation between gas velocity dispersion and $\Sigma_{\rm SFR}$ and a higher scatter of gas velocity dispersion with respect to $\Sigma_{\rm SFR}$, in agreement with the suggestion that other mechanisms, such as gas transportation and accretion, are relatively more important at higher redshift galaxies., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

We present JWST NIRCam imaging targeting 13 $z\sim3$ infrared-luminous ($L_{\rm IR}\sim5\times10^{12}L_{\odot}$) galaxies from the ALESS survey with uniquely deep, high-resolution (0.08$''$$-$0.16$''$) ALMA 870$\mu$m imaging. The 2.0$-$4.4$\mu$m (observed frame) NIRCam imaging reveals the rest-frame near-infrared stellar emission in these submillimeter-selected galaxies (SMGs) at the same (sub-)kpc resolution as the 870$\mu$m dust continuum. The newly revealed stellar morphologies show striking similarities with the dust continuum morphologies at 870$\mu$m, with the centers and position angles agreeing for most sources, clearly illustrating that the spatial offsets reported previously between the 870$\mu$m and HST morphologies were due to strong differential dust obscuration. The F444W sizes are 78$\pm$21% larger than those measured at 870$\mu$m, in contrast to recent results from hydrodynamical simulations that predict larger 870$\mu$m sizes. We report evidence for significant dust obscuration in F444W for the highest-redshift sources, emphasizing the importance of longer-wavelength MIRI imaging. The majority of the sources show evidence that they are undergoing mergers/interactions, including tidal tails/plumes -- some of which are also detected at 870$\mu$m. We find a clear correlation between NIRCam colors and 870$\mu$m surface brightness on $\sim$1 kpc scales, indicating that the galaxies are primarily red due to dust -- not stellar age -- and we show that the dust structure on $\sim$kpc-scales is broadly similar to that in nearby galaxies. Finally, we find no strong stellar bars in the rest-frame near-infrared, suggesting the extended bar-like features seen at 870$\mu$m are highly obscured and/or gas-dominated structures that are likely early precursors to significant bulge growth., Comment: 26 pages, 10 figures. Submitted to ApJ. Comments welcome

The 2175 {\AA} bump shows considerable variations in its strength, width, and central wavelength when observed along different sightlines in the Milky Way and other galaxies. These variations offer valuable insights into the composition, size distribution, and processing of interstellar dust grains along different sightlines. This paper introduces a mission concept called UVESS (Ultraviolet Extinction Sky Survey) aimed at exploring the composition of the interstellar medium (ISM) within both the Milky Way and nearby Local Group Galaxies by mapping the variation of UV extinction curve slopes and the 2175 {\AA} feature across a majority of the sky to gain insights into the makeup of the ISM. Recent advancements in UV instrumentation and technologies pave the way for the development of high-throughput instruments in compact form factors. In this paper, we outline mission science goals and instrument concept tailored for a small satellite-based platform dedicated to the study of UV extinction.

We test the influence of spiral arms on the star formation activity of disk galaxies by constructing and fitting multi-wavelength SEDs for the two nearby spiral galaxies NGC 628 and NGC 4321, at a spatial scale of 1-1.5kpc scale. Recent results in the literature support the 'gatherers' picture, i.e., that spiral arms gather material but do not trigger star formation. However, ambiguities in the diagnostics used to measure star formation rates (SFRs) and other quantities have hampered attempts at reaching definite conclusions. We approach this problem by utilizing the physical parameters output of the MAGPHYS fitting code, which we apply to the Ultraviolet-to-Far Infrared (UV-to-FIR) photometry, in $\geq$20 bands, of spatially-resolved regions in the two galaxies. We separate the regions in arm and interarm, and study the distributions of the specific SFRs (sSFRs=SFR/M$_{star}$), stellar ages and star formation efficiency (SFE=SFR/M$_{gas}$). We find that the distributions of these parameters in the arm regions are almost indistinguishable from those in the interarm regions, with typical differences of a factor 2 or less in the medians. These results support the 'gatherer' scenario of spiral arms, which we plan to test with a larger sample in the near future., Comment: 28 pages, 23 figures, 3 tables; accepted for publication in ApJ

We use the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3-4 Gyr, as probed by the observational proxy for spin, $\lambda_{R}$. We use 2D stellar kinematics to measure $\lambda_{R}$ along with detailed photometric models to estimate galaxy ellipticity. The combination of these measurements quantifies the kinematic classes of `fast rotators' and the rarer `slow rotators', which show no regular rotation in their line-of-sight velocity fields. We compare 51 MAGPI galaxies with $\log_{10} (M_{\star}/\mathrm{M}_\odot) > 10$ to carefully drawn samples of MaNGA galaxies in the local Universe, selected to represent possible descendants of the MAGPI progenitors. The EAGLE simulations are used to identify possible evolutionary pathways between the two samples, explicitly accounting for progenitor bias in our results and the varied evolutionary pathways a galaxy might take between the two epochs. We find that the occurrence of slow rotating galaxies is unchanged between the MAGPI ($z \sim 0.3$) and MaNGA ($z \sim 0$) samples, suggesting the massive slow rotator population was already in place $\sim 4$ Gyr ago and has not accumulated since. There is a hint of the MAGPI sample having an excess of high $\lambda_{R}$ galaxies compared to the MaNGA sample, corresponding to more ordered rotation, but statistically the samples are not significantly different. The large-scale stellar kinematics, as quantified through the $\lambda_{R}$ parameter, of galaxies at $z \sim 0.3$ have already evolved into the diversity of structures seen today in the local Universe., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

The nature of the first galaxies that reionized the universe during the Epoch of Reionization (EoR) remains unclear. Attempts to directly determine spectral properties of these early galaxies are affected by both limited photometric constraints across the spectrum and by the opacity of the intergalactic medium (IGM) to the Lyman Continuum (LyC) at high redshift. We approach this by analysing properties of analogous extreme emission line galaxies (EELGs, [OIII]+Hbeta EW $>400$ Angstrom) at $2.55.5$ observed with the James Webb Space Telesope (JWST) with self-consistent spectral energy distribution fitting methodology. This work focuses on the comparison of their UV slopes ($\beta_P$), ionizing photon production efficiencies $\xi_{ion}$, star formation rates and dust properties to determine the effectiveness of this analogue selection technique. We report the median ionizing photon production efficiencies as log$_{10}(\xi_{ion}/(Hz\ {\rm erg}^{-1}))=$$25.14^{+0.06}_{-0.04}$,$25.16^{+0.06}_{-0.05}$,$25.16^{+0.04}_{-0.05}$,$25.18^{+0.06}_{-0.07}$ for our ZFOURGE control, ZFOURGE EELG, JADES and CEERS samples respectively. ZFOURGE EELGs are 0.57 dex lower in stellar mass and have half the dust extinction, compared to their ZFOURGE control counterparts. They also have a similar specific star formation rates and $\beta_P$ to the $z>5.5$ samples. We find that EELGs at low redshift ($2.5

We present the molecular gas content and ISM conditions of MACSJ0717 Az9, a strong gravitationally lensed $z=4.273$, $M_{*} \simeq 2\times10^9M_{\odot}$ star-forming galaxy with an unusually high ($\sim 80\%$) obscured star formation fraction. We detect CO(4-3) in two independent lensed images, as well as [N II]205$\mu$m, with ALMA. We derive a molecular gas mass of log$_{10}[M_{H_{2}} (M_{\odot})] = 9.77$ making it moderately deficient in molecular gas compared to the lower redshift gas fraction scaling relation. Leveraging photodissociation region (PDR) models, we combine our CO(4-3) measurements with existing measurements of the [C II] 158$\mu$m line and total infrared luminosity to model the PDR conditions. We find PDR conditions similar to local star-forming galaxies, with a mean hydrogen density log$_{10}$[$n_H$ $cm^{-3}$] = $4.80\pm0.39$ and a mean radiation field strength log$_{10}$[G$_0$ Habing] = $2.83\pm0.26$. Based on Band 3 continuum data, we derive an upper limit on the intrinsic dust mass of log$_{10}[M_{\rm dust} (M_{\odot})] < 7.73$, consistent with existing estimates. We use the 3D tilted-ring model fitting code 3D-Barolo to determine the kinematic properties of the CO(4-3) emitting gas. We find that it is rotationally dominated, with a $V/\sigma=4.6 \pm 1.7$, consistent with the kinematics of the [C II]. With PDR conditions remarkably similar to normal dusty star-forming galaxies at z ~ 0.2 and a stable molecular disk, our observations of Az9 suggest that the dust-obscured phase for a low-mass galaxy at z$\sim$4 is relatively long. Thus, Az9 may be representative of a more widespread population that has been missed due to insufficiently deep existing millimeter surveys., Comment: 19 pages, 8 figures

Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial trends in the star formation (SF) activity and luminosity-weighted stellar ages as a function of offset from the star-forming main sequence (SFMS) for a total of 294 galaxies. Using both H$\alpha$ emission and the 4000 Angstrom break (i.e., D4000) as star formation rate (SFR) tracers, we find overall flat radial profiles for galaxies lying on and above the SFMS, suggestive of physical processes that enhance/regulate SF throughout the entire galaxy disc. However, for galaxies lying below the SFMS, we find positive gradients in SF suggestive of inside-out quenching. Placing our results in context with results from other redshift regimes suggests an evolution in radial trends at $z$ ~ 0.3 for SF galaxies above the SFMS, from uniformly enhanced SF at $z$ ~ 1 and $z$ ~ 0.3 to centrally enhanced SF at $z$ ~ 0 (when averaged over a wide range of mass). We also capture higher local SFRs for galaxies below the SFMS compared to that of $z$ ~ 0, which can be explained by a larger population of quenched satellites in the local Universe and/or different treatments of limitations set by the D4000-sSFR relation., Comment: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

Stellar feedback-driven outflows are important regulators of the gas-star formation cycle. However, resolving outflow physics requires high resolution observations that can only be achieved in very nearby galaxies, making suitable targets rare. We present the first results from the new VLT/MUSE large program MAUVE (MUSE and ALMA Unveiling the Virgo Environment), which aims to understand the gas-star formation cycle within the context of the Virgo cluster environment. Outflows are a key part of this cycle, and we focus on the peculiar galaxy NGC 4383, which hosts a $\sim6\,$kpc bipolar outflow fuelled by one of Virgo's most HI-rich discs. The spectacular MUSE data reveal the clumpy structure and complex kinematics of the ionised gas in this M82-like outflow at 100 pc resolution. Using the ionised gas geometry and kinematics we constrain the opening half-angle to $\theta=25-35^\circ$, while the average outflow velocity is $\sim210$ kms$^{-1}$. The emission line ratios reveal an ionisation structure where photoionisation is the dominant excitation process. The outflowing gas shows a marginally elevated gas-phase oxygen abundance compared to the disc but is lower than the central starburst, highlighting the contribution of mixing between the ejected and entrained gas. Making some assumptions about the outflow geometry, we estimate an integrated mass outflow-rate of $\sim1.8~$M$_\odot$yr$^{-1}$ and a corresponding mass-loading factor in the range 1.7-2.3. NGC 4383 is a useful addition to the few nearby examples of well-resolved outflows, and will provide a useful baseline for quantifying the role of outflows within the Virgo cluster., Comment: 16 pages, 7 figures, 1 appendix. Accepted to MNRAS

Geostatistical methods are powerful tools for understanding the spatial structure of the metallicity distribution of galaxies, and enable construction of accurate predictive models of the 2D metallicity distribution. However, so far these methods have only been applied to very high spatial resolution metallicity maps, leaving it uncertain if they will work on lower quality data. In this study, we apply geostatistical techniques to high-resolution spectroscopic maps of three local galaxies convolved to eight different spatial resolutions ranging from ~40pc to ~1 kpc per pixel. We fit a geostatistical model to the data at all resolutions, and find that for metallicity maps where small scale structure is visible by eye (with > ~10 resolution elements per Re), all parameters, including the metallicity correlation scale, can be recovered accurately. At all resolutions tested, we find that point metallicity predictions from such a geostatistical model outperform a circularly symmetric metallicity gradient model. We also explore dependence on the number of data points, and find that N > ~100 spatially resolved metallicity values are sufficient to train a geostatistical model that yields more accurate metallicity predictions than a radial gradient model. Finally, we investigate the potential detrimental effects of having spaxels smaller than an individual Hii region by repeating our analysis with metallicities integrated over Hii regions. We see that spaxel-based measurements have more noise, as expected, but the underlying spatial metallicity distribution can be recovered regardless of whether spaxels or integrated regions are used., Comment: 27 pages, including three appendices. 17 figures in main text, and 7 supplementary figures in appendices. Accepted for publication in MNRAS

We analyze TYPHOON long slit absorption line spectra of the starburst barred spiral galaxy NGC 1365 obtained with the Progressive Integral Step Method covering an area of 15 square kpc. Applying a population synthesis technique, we determine the spatial distribution of ages and metallicity of the young and old stellar population together with star formation rates, reddening, extinction and the ratio R$_V$ of extinction to reddening. We detect a clear indication of inside-out growth of the stellar disk beyond 3 kpc characterized by an outward increasing luminosity fraction of the young stellar population, a decreasing average age and a history of mass growth, which was finished 2 Gyrs later in the outermost disk. The metallicity of the young stellar population is clearly super solar but decreases towards larger galactocentric radii with a gradient of -0.02 dex/kpc. On the other hand, the metal content of the old population does not show a gradient and stays constant at a level roughly 0.4 dex lower than that of the young population. In the center of NGC 1365 we find a confined region where the metallicity of the young population drops dramatically and becomes lower than that of the old population. We attribute this to infall of metal poor gas and, additionally, to interrupted chemical evolution where star formation is stopped by AGN and supernova feedback and then after several Gyrs resumes with gas ejected by stellar winds from earlier generations of stars. We provide a simple model calculation as support for the latter.

Spiral structures are important drivers of the secular evolution of disc galaxies, however, the origin of spiral arms and their effects on the development of galaxies remain mysterious. In this work, we present two three-armed spiral galaxies at z~0.3 in the Middle Age Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey. Taking advantage of the high spatial resolution (~0.6'') of the Multi-Unit Spectroscopic Unit (MUSE), we investigate the two-dimensional distributions of different spectral parameters: Halpha, gas-phase metallicity, and D4000. We notice significant offsets in Halpha (~0.2 dex) as well as gas-phase metallicities (~0.05 dex) among the spiral arms, downstream and upstream of MAGPI1202197197 (SG1202). This observational signature suggests the spiral structure in SG1202 is consistent with arising from density wave theory. No azimuthal variation in Halpha or gas-phase metallicities is observed in MAGPI1204198199 (SG1204), which can be attributed to the tighter spiral arms in SG1204 than SG1202, coming with stronger mixing effects in the disc. The absence of azimuthal D4000 variation in both galaxies suggests the stars at different ages are well-mixed between the spiral arms and distributed around the disc regions. The different azimuthal distributions in Halpha and D4000 highlight the importance of time scales traced by various spectral parameters when studying 2D distributions in spiral galaxies. This work demonstrates the feasibility of constraining spiral structures by tracing interstellar medium (ISM) and stellar population at z~0.3, with a plan to expand the study to the full MAGPI survey., Comment: 15 pages, 11 figures, 2 tables, accepted for publication in MNRAS

We report ALMA detections of [CII] and dust continuum in Az9, a multiply-imaged galaxy behind the Frontier Field cluster MACSJ0717.5+3745. The bright [CII] emission line provides a spectroscopic redshift of z = 4.274. This strongly lensed (mu = 7 +/- 1) galaxy has an intrinsic stellar mass of only 2e9 Msun and a total star formation rate of 26 Msun/yr (~80% of which is dust obscured). Using public magnification maps, we reconstruct the [CII] emission in the source plane to reveal a stable, rotation-dominated disk with V/sigma = 5.3, which is > 2x higher than predicted from simulations for similarly high-redshift, low-mass galaxies. In the source plane, the [CII] disk has a half-light radius of 1.8 kpc and, along with the dust, is spatially offset from the peak of the stellar light by 1.4 kpc. Az9 is not deficient in [CII]; L[CII]/LIR = 0.0027 consistent with local and high redshift normal star forming galaxies. While dust-obscured star formation is expected to dominate in higher mass galaxies, such a large reservoir of dust and gas in a lower mass disk galaxy 1.4 Gyr after the Big Bang challenges our picture of early galaxy evolution. Furthermore, the prevalence of such low-mass dusty galaxies has important implications for the selection of the highest redshift dropout galaxies with JWST. As one of the lowest stellar mass galaxies at z > 4 to be detected in dust continuum and [CII], Az9 is an excellent laboratory in which to study early dust enrichment in the interstellar medium., Comment: Accepted for publication in ApJ Letters

We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameterised by the combined mass density slope $\gamma$ (baryons $+$ dark matter), which is the logarithmic change of density with radius. We use a MAGPI sample of 28 galaxies from low-to-mid density environments and compare to density slopes derived from galaxies in the high density Frontier Fields clusters in the redshift range $0.29

We report spatially resolved dust properties of the quasar host galaxy BRI 1335-0417 at redshift $z = 4.4$ constrained by the ALMA observations. The dust temperature map, derived from a greybody fit to rest-frame 90 and 161~$\mu$m continuum images, shows a steep increase towards the centre, reaching $57.1 \pm 0.3$ K and a flat median profile at the outer regions of $\sim$38 K. Image decomposition analysis reveals the presence of a point source in both dust continuum images spatially coincident with the highest temperature peak and the optical quasar position, which we attribute to warm dust heated by an active galactic nucleus (AGN). We show that a model including this warm component along with cooler dust heated by star formation describes the global SED better than a single component model, with dust temperatures of 87.1$^{+34.1}_{-18.3}$ K (warm component) and 52.6$^{+10.3}_{-11.0}$ K (cold component). The star formation rate (SFR) estimated from the cold dust component is $1700_{-400}^{+500} M_\odot$ yr$^{-1}$, a factor of three smaller than previous estimates due to a large AGN contribution ($53^{+14}_{-15}$\%). The unresolved warm dust component also explains the steep temperature gradient, as the temperature profile derived after the point source subtraction is flat. The point source subtraction also reduces the estimated central SFR surface density $\Sigma_{\mathrm{SFR}}$ by over a factor of three. With this correction, spatially resolved measurements of $\Sigma_{\mathrm{SFR}}$ and the surface gas mass density $\Sigma_{\mathrm{gas}}$ form a roughly linear sequence in the Kennicutt-Schmidt diagram with a constant gas depletion time of 50-200 Myr. The demonstrated AGN-host galaxy decomposition reveals the importance of spatially resolved data for accurate measurements of quasar host galaxy properties, including dust temperature, star-formation rates, and size., Comment: 26pages, 27figures, 4tables, Published in MNRAS; comments are warmly welcomed

Previous studies have shown that the normalization and scatter of the galaxy 'main sequence' (MS), the relation between star formation rate (SFR) and stellar mass ($M_*$), evolves over cosmic time. However, such studies often rely on photometric redshifts and/or only rest-frame UV to near-IR data, which may underestimate the SFR and $M_*$ uncertainties. We use MAGPHYS+photo-z to fit the UV to radio spectral energy distributions of 12,380 galaxies in the COSMOS field at $0.5

Bars play an important role in mixing material in the inner regions of galaxies and stimulating radial migration. Previous observations have found evidence for the impact of a bar on metallicity gradients but the effect is still inconclusive. We use the TYPHOON/PrISM survey to investigate the metallicity gradients along and beyond the bar region across the entire star-forming disk of five nearby galaxies. Using emission line diagrams to identify star-forming spaxels, we recover the global metallicity gradients ranging from -0.0162 to -0.073 dex/kpc with evidence that the galactic bars act as an agent in affecting in-situ star formation as well as the motions of gas and stars. We observe cases with a `shallow-steep' metallicity radial profile, with evidence of the bar flattening the metallicity gradients inside the bar region (NGC~5068 and NGC~1566) and also note instances where the bar appears to drive a steeper metallicity gradient producing `steep-shallow' metallicity profiles (NGC~1365 and NGC~1744). For NGC~2835, a `steep-shallow' metallicity gradient break occurs at a distance $\sim$ 4 times the bar radius, which is more likely driven by gas accretion to the outskirt of the galaxy instead of the bar. The variation of metallicity gradients around the bar region traces the fluctuations of star formation rate surface density in NGC~1365, NGC~1566 and NGC~1744. A larger sample combined with hydrodynamical simulations is required to further explore the diversity and the relative importance of different ISM mixing mechanisms on the gas-phase metallicity gradients in local galaxies., Comment: 17 pages, 6 figures, 2 tables, published on MNRAS

The metallicity of diffuse ionised gas (DIG) cannot be determined using strong emission line diagnostics, which are calibrated to calculate the metallicity of Hii regions. Because of this, resolved metallicity maps from integral field spectroscopy (IFS) data remain largely incomplete. In this paper (the second of a series), we introduce the geostatistical technique of universal kriging, which allows the complete 2D metallicity distribution of a galaxy to be reconstructed from metallicities measured at Hii regions, accounting for spatial correlations between nearby data points. We apply this method to construct high-fidelity metallicity maps of the local spiral galaxy NGC 5236 using data from the TYPHOON/PrISM survey. We find significant correlation in the metallicity of Hii regions separated by up to 0.4-1.2 kpc. Predictions constructed using this method were tested using cross-validation in Hii regions, and we show that they outperform significantly interpolation based on metallicity gradients. Furthermore, we apply kriging to predict the metallicities in regions dominated by DIG emission, considering seven additional spiral galaxies with high resolution (<100pc) metallicity maps. We compare kriging maps to DIG metallicities computed with novel ionisation corrections, and find that such corrections introduce a systematic offset of up to $\pm0.1$ dex for any individual galaxy, with a scatter of 0.02-0.07 dex for the sample. Overall we recommend universal kriging, together with a calibrated geostatistical model, as the superior method for inferring the metallicities of DIG-dominated regions in local spiral galaxies, demonstrating further the potential of applying geostatistical methods to spatially resolved galaxy observations., Comment: 26 pages, 11 figures, including 2 appendices. Accepted for publication in MNRAS

The UV bump is a broad absorption feature centered at 2175{\AA} that is seen in the attenuation/extinction curve of some galaxies, but its origin is not well known. Here, we use a sample of 86 star-forming galaxies at z=1.7-2.7 with deep rest-frame UV spectroscopy from the MUSE HUDF Survey to study the connection between the strength of the observed UV 2175{\AA} bump and the Spitzer/MIPS 24 micron photometry, which at the redshift range of our sample probes mid-IR polycyclic aromatic hydrocarbon (PAH) emission at ~6-8 micron. The sample has robust spectroscopic redshifts and consists of typical main-sequence galaxies with a wide range in stellar mass (log(Mstar/Msun) ~ 8.5-10.7) and star formation rates (SFRs; SFR ~ 1-100 Msun/yr). Galaxies with MIPS detections have strong UV bumps, except for those with mass-weighted ages younger than ~150 Myr. We find that the UV bump amplitude does not change with SFR at fixed stellar mass but increases with mass at fixed SFR. The UV bump amplitude and the PAH strength (defined as mid-IR emission normalized by SFR) are highly correlated and both also correlate strongly with stellar mass. We interpret these correlations as the result of the mass-metallicity relationship, such that at low metallicities PAH emission is weak due to a lower abundance of PAH molecules. The weak or complete absence of the 2175{\AA} bump feature on top of the underlying smooth attenuation curve at low mass/metallicities is then expected if the PAH carriers are the main source of the additional UV absorption., Comment: 9 pages, 7 figures, accepted for publication in mnras

We identify a sample of spectroscopically measured emission line galaxy (ELG) pairs up to z=1.6 from the WFC3 Infrared Spectroscopic Parallels (WISP) survey. WISP obtained slitless, near-infrared grism spectroscopy along with direct imaging in the J and H bands by observing in the pure-parallel mode with the Wide Field Camera Three (WFC3) on the Hubble Space Telescope (HST). From our search of 419 WISP fields covering an area of ~0.5 deg$^{2}$, we find 413 ELG pair systems, mostly Halpha emitters. We then derive reliable star formation rates (SFRs) based on the attenuation-corrected Halpha fluxes. Compared to isolated galaxies, we find an average SFR enhancement of 40%-65%, which is stronger for major pairs and pairs with smaller velocity separations (Delta_v < 300 km/s). Based on the stacked spectra from various subsamples, we study the trends of emission line ratios in pairs, and find a general consistency with enhanced lower-ionization lines. We study the pair fraction among ELGs, and find a marginally significant increase with redshift $f \propto (1+z)^\alpha$, where the power-law index \alpha=0.58$\pm$0.17 from $z\sim$0.2 to $z\sim$1.6. The fraction of Active galactic Nuclei (AGNs), is found to be the same in the ELG pairs as compared to isolated ELGs., Comment: 22 pages, 9 figures, 6 tables, accepted for publication in ApJ

Multi-wavelength images from the farUV (~0.15 micron) to the sub-millimeter of the central region of the galaxy NGC 3351 are analyzed to constrain its stellar populations and dust attenuation. Despite hosting a ~1 kpc circumnuclear starburst ring, NGC 3351 deviates from the IRX-beta relation, the relation between the infrared-to-UV luminosity ratio and the UV continuum slope (beta) that other starburst galaxies follow. To understand the reason for the deviation, we leverage the high angular resolution of archival nearUV-to-nearIR HST images to divide the ring into ~60-180 pc size regions and model each individually. We find that the UV slope of the combined intrinsic (dust-free) stellar populations in the central region is redder than what is expected for a young model population. This is due to the region's complex star formation history, which boosts the nearUV emission relative to the farUV. The resulting net attenuation curve has a UV slope that lies between those of the starburst attenuation curve (Calzetti et al. 2000) and the Small Magellanic Cloud extinction curve; the total-to-selective attenuation value, R'(V)=4.93, is larger than both. As found for other star-forming galaxies, the stellar continuum of NGC 3351 is less attenuated than the ionized gas, with E(B-V)_{star}=0.40 E(B-V)_{gas}. The combination of the `red' intrinsic stellar population and the new attenuation curve fully accounts for the location of the central region of NGC 3351 on the IRX-beta diagram. Thus, the observed characteristics result from the complex mixture of stellar populations and dust column densities in the circumnuclear region. Despite being a sample of one, these findings highlight the difficulty of defining attenuation curves of general applicability outside the regime of centrally-concentrated starbursts., Comment: 45 page, 17 figures; accepted for publication on the Astrophysical Journal, March 21st, 2021

We provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) - stellar mass ($M_*$) relation, called the main sequence of star-forming galaxies (MS), for galaxies out to $z\sim5$. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for $\sim$200,000 mass-selected galaxies at $z>0.3$ in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log($M_*$/M$_\odot$)$<$10) and a flattening at high stellar mass that becomes more prominent at low redshift ($z<1.5$). We find that the SFR density peaks at $1.5

We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies., Comment: 20 pages and 9 figures and 1 appendix, accepted for publication in ApJ

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation., Comment: 33 pages, 10 figures; accepted for publication in ApJS. The catalog tables will be available on ApJS and VizieR; authors' version is available now upon request

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe., Comment: 12 pages, 6 figures, accepted for publication in ApJ

The mass fraction of dense gas within giant molecular clouds (GMCs) of the Milky Way is investigated using 13CO data from the FCRAO Galactic Plane Surveys and the Bolocam Galactic Plane Survey (BGPS) of 1.1 mm dust continuum emission. A sample of 860 compact dust sources are selected from the BGPS catalog and kinematically linked to 344 clouds of extended (>3') 13CO J=1-0 emission. Gas masses are tabulated for the full dust source and subregions within the dust sources with mass surface densities greater than 200 M_solar/pc^2, which are assumed to be regions of enhanced volume density. Masses of the parent GMCs are calculated assuming optically thin 13CO J=1-0 emission and LTE conditions. The mean fractional mass of dust sources to host GMC mass is 0.11 (-0.06, +0.12). The high column density subregions comprise 0.07 (-0.05, +0.13) of the mass of the cloud. Owing to our assumptions, these values are upper limits to the true mass fractions. The fractional mass of dense gas is independent of GMC mass and gas surface density. The low dense gas mass fraction suggests that the formation of dense structures within GMCs is the primary bottleneck for star formation. The distribution of velocity differences between the dense gas and the low density material along the line of sight is also examined. We find a strong, centrally peaked distribution centered on zero velocity displacement. This distribution of velocity differences is modeled with radially converging flows towards the dense gas position that are randomly oriented with respect to the observed line of sight. These models constrain the infall velocities to be 2-4 km/s for various flow configurations., Comment: 12 pages, 6 figures, 1 table. Accepted for publication in ApJ

We introduce the Stony Brook / SMARTS Atlas of (mostly) Southern Novae. This atlas contains both spectra and photometry obtained since 2003. The data archived in this atlas will facilitate systematic studies of the nova phenomenon and correlative studies with other comprehensive data sets. It will also enable detailed investigations of individual objects. In making the data public we hope to engender more interest on the part of the community in the physics of novae. The atlas is on-line at \url{http://www.astro.sunysb.edu/fwalter/SMARTS/NovaAtlas/} ., Comment: 11 figures; 5 tables