Search

This paper presents a 610 MHz radio survey covering 1.94 square degrees around the North Ecliptic Pole (NEP), which includes parts of the AKARI (ADF-N) and Euclid, Deep Fields North. The median 5-sigma sensitivity is 28 microJy beam per beam, reaching as low as 19 microJy per beam, with a synthesised beam of 3.6 x 4.1 arcsec. The catalogue contains 1675 radio components, with 339 grouped into multi-component sources and 284 isolated components likely part of double radio sources. Imaging, cataloguing, and source identification are presented, along with preliminary scientific results. From a non-statistical sub-set of 169 objects with multi-wavelength AKARI and other detections, luminous infrared galaxies (LIRGs) represent 66 percent of the sample, ultra-luminous infrared galaxies (ULIRGs) 4 percent, and sources with L_IR < 1011 L_sun 30 percent. In total, 56 percent of sources show some AGN presence, though only seven are AGN-dominated. ULIRGs require three times higher AGN contribution to produce high-quality SED fits compared to lower luminosity galaxies, and AGN presence increases with AGN fraction. The PAH mass fraction is insignificant, although ULIRGs have about half the PAH strength of lower IR-luminosity galaxies. Higher luminosity galaxies show gas and stellar masses an order of magnitude larger, suggesting higher star formation rates. For LIRGs, AGN presence increases with redshift, indicating that part of the total luminosity could be contributed by AGN activity rather than star formation. Simple cross-matching revealed 13 ROSAT QSOs, 45 X-ray sources, and 61 sub-mm galaxies coincident with GMRT radio sources.

We present a catalog of optical and infrared identifications (ID) of X-ray sources in the AKARI North Ecliptic Pole (NEP) Deep field detected with Chandra covering $\sim 0.34\,{\rm deg^{2}}$ with 0.5-2 keV flux limits ranging $\sim 2 \mathrm{-} 20\times 10^{-16}\,{\rm erg\,s^{-1}\,cm^{-2}}$. The optical/near-infrared counterparts of the X-ray sources are taken from our Hyper Suprime Cam (HSC)/Subaru and Wide-Field InfraRed Camera (WIRCam)/Canada-France-Hawaii Telescope (CFHT) data because these have much more accurate source positions due to their spatial resolution than that of {Chandra} and longer wavelength infrared data. We concentrate our identifications in the HSC $g$ band and WIRCam $K_{\rm s}$ band-based catalogs. To select the best counterpart, we utilize a novel extension of the likelihood-ratio (LR) analysis, where we use the X-ray flux as well as $g - K_{\rm s}$ colors to calculate the likelihood ratio. Spectroscopic and photometric redshifts of the counterparts are summarized. Also, simple X-ray spectroscopy is made on the sources with sufficient source counts. We present the resulting catalog in an electronic form. The main ID catalog contains 403 X-ray sources and includes X-ray fluxes, luminosities, $g$ and $K_{\rm s}$ band magnitudes, redshifts, and their sources, optical spectroscopic properties, as well as intrinsic absorption column densities and power-law indices from simple X-ray spectroscopy. The identified X-ray sources include 27 Milky-Way objects, 57 type I AGNs, 131 other AGNs, and 15 galaxies. The catalog serves as a basis for further investigations of the properties of the X-ray and near-infrared sources in this field. (Abridged), Comment: 16 pages, 9 figures, Three electronic (fits) tables are included in src. Accepted to Astronomy and Astrophysics

We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and sensitivity in the infrared, now reveal that there are more than ten discrete regions with outflows, nearly doubling the number of such known regions. Outflows exist near the radio lobes, at the nucleus, in a biconical structure perpendicular to the jet, and in a bubble moving against the disk. In some of them, velocities above escape velocity are observed. Stratification is also observed, with higher ionization or excitation gas attaining higher velocities. More outflows and bow shocks, found further away from the nucleus than the radio lobes, in regions without significant radio emission, reveal the existence of past or weak radio jets that interacted with the interstellar medium. The coincidence of the bow shocks with the optical extended emission line region (EELR) suggests that the jets also contributed to the gas ionization. Maps of the H2 gas excitation temperature, T_ex, indicate that the molecular gas is most excited in regions with radio emission. There, T_ex is more than 100K higher than in the EELR interior. We argue that a combination of jet-related shocks and cosmic rays is likely responsible for this excess molecular gas excitation., Comment: Submitted to ApJ

This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$\alpha$ forest and Milky Way satellite galaxies.

We follow up emission line galaxies identified through the near-infrared slitless HST/WFC3 WISP survey with VLT/FORS2 optical spectroscopy. Over 4 WISP fields, we targetted 85 of 138 line emission objects at $0.4

We present reduced images and catalogues of photometric and emission line data ($\sim$230,000 and $\sim$8,000 sources, respectively) for the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. These data are made publicly available on the Mikulski Archive for Space Telescopes (MAST) and include reduced images from various facilities: ground-based $ugri$, HST WFC3, and Spitzer IRAC (Infrared Array Camera). Coverage in at least one additional filter beyond the WFC3/IR data are available for roughly half of the fields (227 out of 483), with $\sim$20% (86) having coverage in six or more filters from $u$-band to IRAC 3.6$\mu$m (0.35-3.6$\mu$m). For the lower spatial resolution (and shallower) ground-based and IRAC data, we perform PSF-matched, prior-based, deconfusion photometry (i.e., forced-photometry) using the TPHOT software to optimally extract measurements or upper limits. We present the methodology and software used for the WISP emission line detection and visual inspection. The former adopts a continuous wavelet transformation that significantly reduces the number of spurious sources as candidates before the visual inspection stage. We combine both WISP catalogues and perform SED fitting on galaxies with reliable spectroscopic redshifts and multi-band photometry to measure their stellar masses. We stack WISP spectra as functions of stellar mass and redshift and measure average emission line fluxes and ratios. We find that WISP emission line sources are typically `normal' star-forming galaxies based on the Mass-Excitation diagram ([OIII]/H$\beta$ vs. $M_\star$; $0.74

The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$\mu$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models.

Including millimeter-wave (mm-wave) data in multi-wavelength studies of the variability of active galactic nuclei (AGN) can provide insights into AGN physics that are not easily accessible at other wavelengths. We demonstrate in this work the potential of cosmic microwave background (CMB) telescopes to provide long-term, high-cadence mm-wave AGN monitoring over large fractions of sky. We report on a pilot study using data from the SPTpol instrument on the South Pole Telescope (SPT), which was designed to observe the CMB at arcminute and larger angular scales. Between 2013 and 2016, SPTpol was used primarily to observe a single 500 deg^2 field, covering the entire field several times per day with detectors sensitive to radiation in bands centered at 95 and 150 GHz. We use SPT 150 GHz observations to create AGN light curves, and we compare these mm-wave light curves to those at other wavelengths, in particular gamma-ray and optical. In this Letter, we focus on a single source, PKS 2326-502, which has extensive, day-timescale monitoring data in gamma-ray, optical, and now mm-wave between 2013 and 2016. We find PKS 2326-502 to be in a flaring state in the first two years of this monitoring, and we present a search for evidence of correlated variability between mm-wave, optical R band, and gamma-ray observations. This pilot study is paving the way for AGN monitoring with current and upcoming CMB experiments such as SPT-3G, Simons Observatory, and CMB-S4, including multi-wavelength studies with facilities such as VRO-LSST., Comment: 9 pages, 3 figures, accepted to Astrophysical Journal Letters

We present results from the James Webb Space Telescope (JWST) Director's Discretionary Time Early Release Science (ERS) program 1328 targeting the nearby, Luminous Infrared Galaxy (LIRG), VV 114. We use the MIRI and NIRSpec instruments to obtain integral-field spectroscopy of the heavily obscured Eastern nucleus (V114E) and surrounding regions. The spatially resolved, high-resolution, spectra reveal the physical conditions in the gas and dust over a projected area of 2-3 kpc that includes the two brightest IR sources, the NE and SW cores. Our observations show for the first time spectroscopic evidence that the SW core hosts an AGN as evidenced by its very low 6.2 {\mu}m and 3.3 {\mu}m PAH equivalent widths (0.12 and 0.017 {\mu}m respectively) and mid and near-IR colors. Our observations of the NE core show signs of deeply embedded star formation including absorption features due to aliphatic hydrocarbons, large quantities of amorphous silicates, as well as HCN due to cool gas along the line of sight. We detect elevated [Fe II]/Pf{\alpha} consistent with extended shocks coincident with enhanced emission from warm H$_{2}$, far from the IR-bright cores and clumps. We also identify broadening and multiple kinematic components in both H$_{2}$ and fine structure lines caused by outflows and previously identified tidal features., Comment: 12 pages, 4 figures, accepted for publication in ApJL

We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope SZ survey (SPT-SZ) and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array (ALMA). This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans $1.9$$<$z$<$$6.9$ and covers rest-frame frequencies of 240$-$800 GHz. Combining this data with low-J CO observations from the Australia Telescope Compact Array (ATCA), we detect multiple bright line features from $^{12}$CO, $[$CI$]$, and H$_2$O, as well as fainter molecular transitions from $^{13}$CO, HCN, HCO$^+$, HNC, CN, H$_2$O$^+$, and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright $^{12}$CO emission lines, and contain either warmer and more excited dense gas tracers, or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L$_{\mathrm{IR}}$$>$$10^{12}$L$_\odot$), and will inform studies of chemical evolution before the peak epoch of star formation at $z=2-3$., Comment: 26 pages, 10 figures, accepted for publication in ApJ

We present mid-infrared spectroscopic observations of the nucleus of the nearby Seyfert galaxy NGC 7469 taken with the MIRI instrument on the James Webb Space Telescope (JWST) as part of Directors Discretionary Time Early Release Science (ERS) program 1328. The high resolution nuclear spectrum contains 19 emission lines covering a wide range of ionization. The high ionization lines show broad, blueshifted emission reaching velocities up to 1700 km s$^{-1}$ and FWHM ranging from $\sim500 - 1100$ km s$^{-1}$. The width of the broad emission and the broad to narrow line flux ratios correlate with ionization potential. The results suggest a decelerating, stratified, AGN driven outflow emerging from the nucleus. The estimated mass outflow rate is one to two orders of magnitude larger than the current black hole accretion rate needed to power the AGN. Eight pure rotational H$_{2}$ emission lines are detected with intrinsic widths ranging from FWHM $\sim 125-330$ km s$^{-1}$. We estimate a total mass of warm H$_{2}$ gas of $\sim1.2\times10^{7}$M$_{\odot}$ in the central 100 pc. The PAH features are extremely weak in the nuclear spectrum, but a $6.2\mu$m PAH feature with an equivalent width $\sim0.07\mu$m and a flux of $2.7\times10^{-17}$ W m$^{-2}$ is detected. The spectrum is steeply rising in the mid-infrared, with a silicate strength $\sim0.02$, significantly smaller than seen in most PG QSOs, but comparable to other Seyfert 1's. These early MIRI mid-infrared IFU data highlight the power of JWST to probe the multi-phase interstellar media surrounding actively accreting supermassive black holes.

Multiwavelength variability studies of active galactic nuclei (AGN) can be used to probe their inner regions which are not directly resolvable. Dust reverberation mapping (DRM) estimates the size of the dust emitting region by measuring the delays between the infrared (IR) response to variability in the optical light curves. We measure DRM lags of Zw229-015 between optical ground-based and Kepler light curves and concurrent IR Spitzer 3.6 and 4.5 $\mu$m light curves from 2010-2015, finding an overall mean rest-frame lag of 18.3 $\pm$ 4.5 days. Each combination of optical and IR light curve returns lags that are consistent with each other within 1$\sigma$, which implies that the different wavelengths are dominated by the same hot dust emission. The lags measured for Zw229-015 are found to be consistently smaller than predictions using the lag-luminosity relationship. Also, the overall IR response to the optical emission actually depends on the geometry and structure of the dust emitting region as well, so we use Markov chain Monte Carlo (MCMC) modelling to simulate the dust distribution to further estimate these structural and geometrical properties. We find that a large increase in flux between the 2011-2012 observation seasons, which is more dramatic in the IR light curve, is not well simulated by a single dust component. When excluding this increase in flux, the modelling consistently suggests that the dust is distributed in an extended flat disk, and finds a mean inclination angle of 49$^{+3}_{-13}$ degrees., Comment: 32 pages, 32 Figures, 7 Tables; Accepted for publication in MNRAS

The GLASS James Webb Space Telescope Early Release Science (hereafter GLASS-JWST-ERS) Program will obtain and make publicly available the deepest extragalactic data of the ERS campaign. It is primarily designed to address two key science questions, namely, "what sources ionized the universe and when?" and "how do baryons cycle through galaxies?", while also enabling a broad variety of first look scientific investigations. In primary mode, it will obtain NIRISS and NIRSpec spectroscopy of galaxies lensed by the foreground Hubble Frontier Field cluster, Abell 2744. In parallel, it will use NIRCam to observe two fields that are offset from the cluster center, where lensing magnification is negligible, and which can thus be effectively considered blank fields. In order to prepare the community for access to this unprecedented data, we describe the scientific rationale, the survey design (including target selection and observational setups), and present pre-commissioning estimates of the expected sensitivity. In addition, we describe the planned public releases of high-level data products, for use by the wider astronomical community., Comment: 23 pages, 10 figures; ApJ in press

We present the multi-wavelength counterparts of 850-$\mu$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$\mu$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$\mu$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$., Comment: 24 pages, 16 figures. Accepted for publication in MNRAS

We present the first characterisation of the average dust attenuation curve at $z\sim1.3$ by combining rest-frame ultraviolet through near-IR photometry with Balmer decrement ($\mathrm{H}\alpha$/$\mathrm{H}\beta$) constraints for $\sim$900 galaxies with $8\lesssim\log (M_\star /M_\odot)<10.2$ at $0.75

The identification of an emission line is unambiguous when multiple spectral features are clearly visible in the same spectrum. However, in many cases, only one line is detected, making it difficult to correctly determine the redshift. We developed a freely available unsupervised machine-learning algorithm based on unbiased topology (UMLAUT) that can be used in a very wide variety of contexts, including the identification of single emission lines. To this purpose, the algorithm combines different sources of information, such as the apparent magnitude, size and color of the emitting source, and the equivalent width and wavelength of the detected line. In each specific case, the algorithm automatically identifies the most relevant ones (i.e., those able to minimize the dispersion associated with the output parameter). The outputs can be easily integrated into different algorithms, allowing us to combine supervised and unsupervised techniques and increasing the overall accuracy. We tested our software on WISP (WFC3 IR Spectroscopic Parallel) survey data. WISP represents one of the closest existing analogs to the near-IR spectroscopic surveys that are going to be performed by the future Euclid and Roman missions. These missions will investigate the large-scale structure of the universe by surveying a large portion of the extragalactic sky in near-IR slitless spectroscopy, detecting a relevant fraction of single emission lines. In our tests, UMLAUT correctly identifies real lines in 83.2% of the cases. The accuracy is slightly higher (84.4%) when combining our unsupervised approach with a supervised approach we previously developed.

High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_{\odot}$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar media. In this paper, we present high resolution ($\sim$0.4$^{\prime\prime}$) observations of CO(7-6), [CI](2-1) and dust continuum of 3 lensed galaxies from the SPT-SMG sample at $z\sim$ 3 with the Atacama Large Millimeter/submillimeter Array. Our sources have high intrinsic star-formation rates ($>$850 M$_{\odot}$yr$^{-1}$) and rather short depletion timescales ($<$100 Myr). Based on the L$_{[\rm CI](2-1)}$/L$_{\rm CO(7-6)}$ and L$_{[\rm CI](2-1)}$/L$_{\rm IR}$ ratios, our galaxy sample has similar radiation field intensities and gas densities compared to other submillimetre galaxies. We perform visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H$_2$ conversion factor $\alpha_{\rm CO}$ but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103-45 and SPT2147-50 are likely rotating disks while SPT2357-51 is possibly a major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-$z$ dusty star-forming galaxies., Comment: 29 pages, 26 figures, Accepted for publication in A&A. Figure 16 and the related conclusions have been revised

Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We used galaxies detected in 5 optical bands ($g$, $r$, $i$, $z$, and $Y$) by the Subaru Hyper Suprime-Cam (HSC), assisted with $u$-band from Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam, and IRAC1 and IRAC2 bands from the $Spitzer$ space telescope for photo-$z$ estimation. We calculated the local density around every galaxy using the 10$^{th}$-nearest neighbourhood. Cluster candidates were determined by applying the friends-of-friends algorithm to over-densities. 88 cluster candidates containing 4390 member galaxies below redshift 1.1 in 5.4 deg$^2$ have been detected. The reliability of our method was examined through false detection tests, redshift uncertainty tests, and applications on the COSMOS data, giving false detection rates of 0.01 to 0.05 and recovery rate of 0.9 at high richness. 3 X-ray clusters previously observed by $ROSAT$ and $Chandra$ were recovered. The cluster galaxies show higher stellar mass and lower star formation rate (SFR) compared to the field galaxies in two-sample Z-tests. These cluster candidates are useful for environmental studies of galaxy evolution and future astronomical surveys in the NEP, where $AKARI$ has performed unique 9-band mid-infrared photometry for tens of thousands of galaxies., Comment: 19 pages, 16 figures, has been accepted for publication in MNRAS

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

The $AKARI$ infrared (IR) space telescope conducted two surveys (Deep and Wide) in the North Ecliptic Pole (NEP) field to find more than 100,000 IR sources using its Infrared Camera (IRC). IRC's 9 filters, which cover wavebands from 2 to 24 $\mu$m continuously, make $AKARI$ unique in comparison with other IR observatories such as $Spitzer$ or $WISE$. However, studies of the $AKARI$ NEP-Wide field sources had been limited due to the lack of follow-up observations in the ultraviolet (UV) and optical. In this work, we present the Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam $u$-band source catalogue of the $AKARI$ NEP-Wide field. The observations were taken in 7 nights in 2015 and 2016, resulting in 82 observed frames covering 3.6 deg$^2$. The data reduction, image processing and source extraction were performed in a standard procedure using the \textsc{Elixir} pipeline and the \textsc{AstrOmatic} software, and eventually 351,635 sources have been extracted. The data quality is discussed in two regions (shallow and deep) separately, due to the difference in the total integration time (4,520 and 13,910 seconds). The 5$\sigma$ limiting magnitude, seeing FWHM, and the magnitude at 50 per cent completeness are 25.38 mag (25.79 mag in the deep region), 0.82 arcsec (0.94 arcsec) and 25.06 mag (25.45 mag), respectively. The u-band data provide us with critical improvements to photometric redshifts and UV estimates of the precious infrared sources from the $AKARI$ space telescope., Comment: 13 pages, 14 figures, accepted for publication in MNRAS

Dusty high-z galaxies are extreme objects with high star formation rates (SFRs) and luminosities. Characterising the properties of this population and analysing their evolution over cosmic time is key to understanding galaxy evolution in the early Universe. We select a sample of high-z dusty star-forming galaxies (DSFGs) and evaluate their position on the main sequence (MS) of star-forming galaxies, the well-known correlation between stellar mass and SFR. We aim to understand the causes of their high star formation and quantify the percentage of DSFGs that lie above the MS. We adopted a multi-wavelength approach with data from optical to submillimetre wavelengths from surveys at the North Ecliptic Pole (NEP) to study a submillimetre sample of high-redshift galaxies. Two submillimetre selection methods were used, including: sources selected at 850$\mathrm{\, \mu m}$ with the Sub-millimetre Common-User Bolometer Array 2) SCUBA-2 instrument and {\it Herschel}-Spectral and Photometric Imaging Receiver (SPIRE) selected sources (colour-colour diagrams and 500$\mathrm{\, \mu m}$ risers), finding that 185 have good multi-wavelength coverage. The resulting sample of 185 high-z candidates was further studied by spectral energy distribution (SED) fitting with the CIGALE fitting code. We derived photometric redshifts, stellar masses, SFRs, and additional physical parameters, such as the infrared luminosity and active galactic nuclei (AGN) contribution. We find that the different results in the literature are, only in part, due to selection effects. The difference in measured SFRs affects the position of DSFGs on the MS of galaxies; most of the DSFGs lie on the MS (60\%). Finally, we find that the star formation efficiency (SFE) depends on the epoch and intensity of the star formation burst in the galaxy; the later the burst, the more intense the star formation.

The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870\, \mathrm{\mu m}} > 25\, \mathrm{mJy}$) sample, initially selected at $1.4$ mm. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the $3$ mm spectral window, targeting carbon monoxide line emission. By combining these measurements with ancillary data, the SPT sample is now spectroscopically complete, with redshifts spanning $1.9$$<$$z$$<$$6.9$ and a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3\, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enabling the "high-redshift tail" of extremely luminous DSFGs to be measured. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts., Comment: 56 pages, 25 figures, accepted for publication in ApJ

Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., Comment: 19 pages, 9 figures, ApJ in press

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 \mu m}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy., Comment: 12 pages, 8 figures. Accepted for publication in A&A

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

We present Atacama Compact Array and Atacama Pathfinder Experiment observations of the [N II] 205 $\mu$m fine-structure line in 40 sub-millimetre galaxies lying at redshifts z = 3 to 6, drawn from the 2500 deg$^2$ South Pole Telescope survey. This represents the largest uniformly selected sample of high-redshift [N II] 205 $\mu$m measurements to date. 29 sources also have [C II] 158 $\mu$m line observations allowing a characterization of the distribution of the [C II] to [N II] luminosity ratio for the first time at high-redshift. The sample exhibits a median L$_{[C II]}$ /L$_{[N II]}$ $\approx$ 11 and interquartile range of 5.0 to 24.7. These ratios are similar to those observed in local (U)LIRGs, possibly indicating similarities in their interstellar medium. At the extremes, we find individual sub-millimetre galaxies with L$_{[C II]}$ /L$_{[N II]}$ low enough to suggest a smaller contribution from neutral gas than ionized gas to the [C II] flux and high enough to suggest strongly photon or X-ray region dominated flux. These results highlight a large range in this line luminosity ratio for sub-millimetre galaxies, which may be caused by variations in gas density, the relative abundances of carbon and nitrogen, ionization parameter, metallicity, and a variation in the fractional abundance of ionized and neutral interstellar medium., Comment: 8 pages, 4 figures, accepted to MNRAS

Galaxies grow inefficiently, with only a few percent of the available gas converted into stars each free-fall time. Feedback processes, such as outflowing winds driven by radiation pressure, supernovae or supermassive black hole accretion, can act to halt star formation if they heat or expel the gas supply. We report a molecular outflow launched from a dust-rich star-forming galaxy at redshift 5.3, one billion years after the Big Bang. The outflow reaches velocities up to 800 km/s relative to the galaxy, is resolved into multiple clumps, and carries mass at a rate within a factor of two of the star formation rate. Our results show that molecular outflows can remove a large fraction of the gas available for star formation from galaxies at high redshift., Comment: Author's version of paper to appear in Science on 7 September 2018

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis., Comment: 61 pages, 24 figures, accepted for publication in Astronomy and Astrophysics

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today., Comment: To appear in April 26 issue of Nature

For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups., Comment: 31 pages, 23 figures

The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex., Comment: 21 pages, 8 figures, 2 tables. Accepted for publication in ApJ

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch., Comment: Nature, published online on 06 December 2017

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light., Comment: 10 pages, 11 figures and 3 tables. Accepted to the ApJS

The recent updates of the North Ecliptic Pole deep (0.5~deg$^2$, NEP-Deep) multi-wavelength survey covering from X-ray to radio-wave is presented. The NEP-Deep provides us with several thousands of 15~$\mu$m or 18~$\mu$m selected sample of galaxies, which is the largest sample ever made at this wavelengths. A continuous filter coverage in the mid-infrared wavelength (7, 9, 11, 15, 18, and 24~$\mu$m) is unique and vital to diagnose the contributions from starbursts and AGNs in the galaxies out to $z$=2.The new goal of the project is to resolve the nature of the cosmic star formation history at the violent epoch (e.g. $z$=1--2), and to find a clue to understand its decline from $z$=1 to present universe by utilizing the unique power of the multi-wavelength survey. The progress in this context is briefly mentioned., Comment: 4 pages, 2 figures, AKARI conference proceedings; accepted to PKAS

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths., Comment: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features., Comment: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed., Comment: In press. This paper belongs to the "SPICA Special Issue" on PASA

We built a 8um selected sample of galaxies in the NEP-AKARI field by defining 4 redshift bins with the four AKARI bands at 11, 15, 18 and 24 microns (0.15

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization., Comment: Accepted by ApJL, Fixed compilation problem in v.2