Feedback and outflows associated with a quasar phase are expected to be critical in quenching the most massive galaxies at high-z. Observations targeting the cool molecular and atomic phases, which dominate the mass and momentum budget of massive galaxy outflows and remove the direct fuel for star formation are, however, severely limited in high-z QSO hosts. We discuss two recent ALMA programs: one targeting molecular outflows in 3 z ∼ 6 QSO hosts using the OH 119 μm absorption line and another targeting the diffuse, predominantly atomic gas in the halos surrounding 5 QSO host between z ∼ 2 – 4 using the OH+(11 – 10) absorption line. Outflows are successfully detected in both samples and compared with outflows driven by high-z star-forming galaxies observed in the same lines. Both studies indicate that observing QSOs during the blow-out phase is crucial for studying the impact of the active nucleus on the ejection of gas from the host galaxy.

Using the SDSS spectroscopy, we have carried out fine optical spectral classification for activity types for 710 AGN candidates. These objects come from a larger sample of some 2,500 candidate AGN using pre-selection by various samples; bright objects of the Catalog of Quasars and Active Galactic Nuclei, AGN candidates among X-ray sources, optically variable radio sources, IRAS extragalactic objects, etc. A number of papers have been published with the results of this spectral classification. More than 800 QSOs have been identified and classified, including 710 QSOs, Seyferts and Composites. The fine classification shows that many QSOs show the same features as Seyferts, i.e., subtypes between S1 and S2 (S1.2, S1.5, S1.8 and S1.9). We have introduced subtypes for the QSOs: QSO1.2, QSO1.5, QSO1.8, QSO1.9, though the last subtype does not appear in SDSS wavelength range due to mostly highly redshifted Hα (the main line for identification of the 1.9 subtype). Thus, independent of the luminosity (which serves as a separator between QSOs and Seyferts), AGN show the same features. We also have classified many objects as Composites, spectra having composite characteristics between Sy and LINERs, Sy and HII or LINERs and HII; in some cases all three characteristics appear together resulting as Sy/LINER/HII subtype. The QSOs subtypes together with Seyfert ones allow to follow AGN properties along larger redshift range expanding our knowledge on the evolution of AGN to more distant Universe represented by QSOs.

Nebular Heii emission implies the presence of energetic photons (E≽54 eV). Despite the great deal of effort dedicated to understanding Heii ionization, its origin has remained mysterious, particularly in metal-deficient star-forming galaxies. Unfolding Heii-emitting, metal-poor starbursts at z∼0 can yield insight into the powerful ionization processes occurring in the primordial universe. Here we present a study on the origin of the extended nebular Heii emission in SBS 0335-052E, one of the most metal-poor (Z ∼ 3% Z⊙ Heii-emitter starbursts known locally. Based on optical VLT/MUSE spectroscopic and Chandra X-ray observations, and current stellar models we found that the Heii-ionization budget of SBS 0335-052E can only be produced by peculiar, nearly metal-free ionizing stars (called here “PopIII-like” stars) with a top-heavy initial mass function. This result is in line with recent simulations for PopIII star formation down to z=0.

We analyzed archival HST and IUE ultraviolet spectra of 29 nearby star-forming galaxies. The range of aperture sizes permits studies of the galaxy properties over pc to kpc scales. We measured line strengths and spectral energy distributions over the 1200 – 1300 Åwavelength range and established trends with galaxy properties. Updated oxygen abundances were measured from ancillary optical data. Star-formation rates and internal dust attenuations were derived from comparison with synthesis models. The interstellar absorption lines are heavily saturated, yet scale with oxygen abundance. We interpret this as due to macroscopic velocities arising in a turbulent ISM and large-scale outflows. The stellar-wind lines also scale with oxygen abundance. As these lines are shaped by mass loss, which is driven by the Fe abundance, we can study the α-element/Fe ratio in these galaxies.

Super star cluster (SSC) A1 in starburst galaxy NGC 3125 has the strongest broad He II λ1640 emission line ever observed in the nearby Universe and constitutes an important template for interpreting observations of galaxies that are located out to a redshift of z∼3. We use observations of SSC A1 obtained with the Cosmic Origins Spectrograph (COS) on board of the Hubble Space Telescope (HST) in order to check if there is a contribution of nebular emission to the He II line. In addition, we compare the COS G130M + G160M observations of A1 (1150 – 1750∘A) to the latest Charlot & Bruzual population synthesis models, which account for Very Massive Stars (VMS) of up to 300 Mȯ. A model with Z = 0.008 and age = 2.4 Myr provides a very reasonable fit to the C III λ1175, N V λ1240, C IV λ1550, He II λ1640, and N IV λ1718 stellar-wind features, although the O V λ1371 line is not well reproduced. Overall, our results show the great improvement of stellar evolution and population synthesis models over the past decade, and in particular, the improved formulation of stellar mass loss rates.

The GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) is a radio continuum survey at 76–227 MHz of the entire southern sky (Declination $<\!{+}30^{\circ}$ ) with an angular resolution of ${\approx}2$ arcmin. In this paper, we combine GLEAM data with optical spectroscopy from the 6dF Galaxy Survey to construct a sample of 1 590 local (median $z \approx 0.064$ ) radio sources with $S_{200\,\mathrm{MHz}} > 55$ mJy across an area of ${\approx}16\,700\,\mathrm{deg}^{2}$ . From the optical spectra, we identify the dominant physical process responsible for the radio emission from each galaxy: 73% are fuelled by an active galactic nucleus (AGN) and 27% by star formation. We present the local radio luminosity function for AGN and star-forming (SF) galaxies at 200 MHz and characterise the typical radio spectra of these two populations between 76 MHz and ${\sim}1$ GHz. For the AGN, the median spectral index between 200 MHz and ${\sim}1$ GHz, $\alpha_{\mathrm{high}}$ , is $-0.600 \pm 0.010$ (where $S \propto \nu^{\alpha}$ ) and the median spectral index within the GLEAM band, $\alpha_{\mathrm{low}}$ , is $-0.704 \pm 0.011$ . For the SF galaxies, the median value of $\alpha_{\mathrm{high}}$ is $-0.650 \pm 0.010$ and the median value of $\alpha_{\mathrm{low}}$ is $-0.596 \pm 0.015$ . Among the AGN population, flat-spectrum sources are more common at lower radio luminosity, suggesting the existence of a significant population of weak radio AGN that remain core-dominated even at low frequencies. However, around 4% of local radio AGN have ultra-steep radio spectra at low frequencies ( $\alpha_{\mathrm{low}} < -1.2$ ). These ultra-steep-spectrum sources span a wide range in radio luminosity, and further work is needed to clarify their nature.

Mass-loss and radiation feedback from evolving massive stars produce galactic-scale superwinds, sometimes surrounded by pressure-driven bubbles. Using the time-dependent stellar population typically seen in star-forming regions, we conduct hydrodynamic simulations of a starburst-driven superwind model coupled with radiative efficiency rates to investigate the formation of radiative cooling superwinds and bubbles. Our numerical simulations depict the parameter space where radiative cooling superwinds with or without bubbles occur. Moreover, we employ the physical properties and time-dependent ionization states to predict emission line profiles under the assumption of collisional ionization and non-equilibrium ionization caused by wind thermal feedback in addition to photoionization created by the radiation background. We see the dependence of non-equilibrium ionization structures on the time-evolving ionizing source, leading to a deviation from collisional ionization in radiative cooling wind regions over time.

Ultracompact Hii regions (UC-HII) are the young, very dense cores of massive star-forming regions in dwarf galaxies, where newly formed massive OB stars are surrounded by natal molecular clouds. Thermal energy deposited by mechanical feedback from a cluster of massive OB stars can form a superwind, which may lead to a wind-blown bubble as well as radiative cooling. We investigate the formation of radiatively cooling superwinds in UC-HII using a radiative cooling module in the hydrodynamics program. We built a grid of hydrodynamic simulations to determine the dependence of radiative cooling on the cluster radius, mass-deposition rate, wind velocity, and ambient medium in UC-HII. Our findings could help to better understand star formation in massive star-forming regions, where cool superwinds could trigger the formation of molecular clumpy regions.

IC 10 as a starburst dwarf galaxy in the Local Group (LG) has a large population of newly formed stars that are massive and intrinsically very bright in comparison with other LG galaxies. Using the Isaac Newton Telescope (INT) with the Wide Field Camera (WFC) in the i-band and V-band, we performed an optical monitoring survey to identify the most evolved asymptotic giant branch stars (AGBs) and red supergiant stars (RSGs) in this star-forming galaxy, which can be used to determine the star formation history (SFH). The E(B - V) as an effective factor for obtaining the precise magnitude of stars is measured for each star using a 2D dust map (SFD98) to obtain a total extinction for each star in both the i-band and V-band. We obtained the photometric catalog for 53579 stars within the area of 0.07 deg2 (13.5 kpc2), of which 762 stars are classified as variable candidates after removing the foreground stars and saturated ones from our catalog. To reconstruct the SFH for IC 10, we first identified 424 long-period variable (LPV) candidates within the area of two half-light radii (2rh) from the center of the galaxy. We estimated the recent star formation rate (SFR) at ∼ 0.32 yr-1 for a constant metallicity Z = 0.0008, showing the galaxy is currently undergoing high levels of star formation. Also, a total stellar mass of is obtained within 2rh for that metallicity.

We processed the catalogue data for all snapshots of the Illustris TNG100 cosmological simulation and collected every calculated property of the galaxies formed at different redshifts. With this dataset we can statistically analyze parameters for galaxy samples at given redshifts, as well as trace sample parameters over the entire time range of the simulation. Focusing first on star formation rate (SFR) and metallicity, we see the cosmic star formation history with the mean maximum at around z ≈ 1.6 and the reionization bump at around z ≈ 5, while metallicity increases. For a sample of strongly star-forming galaxies with SFR > 10 M⊙ yr−1 we found different characteristics compared to the whole sample. The mean metallicity of highly star-forming galaxies is higher and changes less, and the mean SFR has its maximum at around the reionization bump.

Ultra-luminous infrared galaxies (ULIRGs) are extreme in many ways. The major mergers trigger star formation at very high rates that cause the ISM to be dominated by infrared (IR) photons. We show the ammonia spectra toward the two cores of Arp 220, the nearest ULIRG, in three Very Large Array (VLA) bands (Ku, K, Ka). Typical decay times of the non-metastable transitions ∼ 100 s and are therefore usually difficult to observe. The FIR excitation of Arp 220, however, shows that non-metastable states are widely populated up to a limiting energy of ∼ 1500 K. We assume that this atypical ammonia spectrum is due to the strong FIR field that re-excites the ammonia molecule on timescales much shorter than the already short decay times. The resulting level population causes a break-down of the typical assumptions made for the use of ammonia as a molecular thermometer.

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. Due to the rarity of mergers in the local Universe, samples of nearby merging galaxies suitable for studies of individual GMCs are limited. Idealized simulations provide us with a new window to study GMC evolution during a merger, and assist in interpreting observations. We conduct a pixel-by-pixel analysis of the simulated molecular gas properties in both undisturbed control galaxies and galaxy mergers. The simulated GMC-pixels follow a similar trend in a diagram of velocity dispersion (σv) versus gas surface density (Σmol) as observed in normal spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For simulated mergers, we see a significant increase in both the Σmol and σv for GMC-pixels by a factor of 5 – 10, which put these pixels to be above the trend of PHANGS galaxies in the σv vs Σmol diagram. This deviation indicates that GMCs in the simulated merger are more gravitationally unbound and have higher virial parameter (αvir) of 10 – 100, which is much larger than that of simulated control galaxies. Furthermore, we find that the increase in αvir generally happens at the same time as the increase in global star formation rate (SFR), which suggests feedback is playing a role in dispersing the gas. The correspondence between high αvir and SFR also suggests some other physical mechanisms besides self-gravity are helping the GMCs in starburst mergers to collapse and form stars.

Existence of the cold-mode gas accretion along with the hot-mode accretion can explain the diversity in the galactic star formation history across galaxy mass. We examine the role of various physical processes in producing the observed diversity.

Thermal energies deposited by OB stellar clusters in starburst galaxies lead to the formation of galactic superwinds. Multi-wavelength observations of starburst-driven superwinds pointed at complex thermal and ionization structures which cannot adequately be explained by simple adiabatic assumptions. In this study, we perform hydrodynamic simulations of a fluid model coupled to radiative cooling functions, and generate time-dependent non-equilibrium photoionization models to predict physical conditions and ionization structures of superwinds using the maihem atomic and cooling package built on the program flash. Time-dependent ionization states and physical conditions produced by our simulations are used to calculate the emission lines of superwinds for various parameters, which allow us to explore implications of non-equilibrium ionization for starburst regions with potential radiative cooling.

We report the first detection of [NII] 122 μm line toward a QSO-SMG pair, BRI 1202-0725, at z = 4.7 using ALMA. Combining with [NII] 205 μm line detection and taking the line ratio of [NII]122/[NII]205, we constrain electron densities of both galaxies. The derived electron densities are ${26_{ - 11}^{ + 12}$ and ${134_{ - 39}^{ + 50}$ cm−3 for the SMG and the QSO, respectively, which are the first measurements for galaxies at z > 4. The electron density of the SMG is comparable to the Galactic plane and the average of local spiral galaxies, while the value for the QSO is comparable to local starbursts and optical-line based measurements for star-forming galaxies at z ∼ 2–3. Considering the similar star-formation rates (SFRs) of ≍ 1000 M⊙ yr−1 for both galaxies, our results suggest a large scatter of electron densities at fixed SFR and caution against using optical lines for dusty starbursts. The details of this report are presented in Lee et al. 2019 (submitted).

We investigate the stellar population of star-forming galaxies at z ∼ 4 by focusing on their slope of rest-frame ultraviolet continuum called UV spectral slope β. We analyze the sample of bright Lyman Break Galaxies (LBGs) with Subaru/i′≤26.0in the Subaru/XMM-Newton Deep Survey field. Our detailed SED fitting analysis indicates that the LBGs with observed UV slope > −1.7, , Av > 1.0, and intrinsic UV slope < −2.5 are the intrinsically active star-forming galaxies with star formation rates larger than a few × 102 M⊙yr−1. A significant fraction of the UV-selected LBGs at z ∼ 4 is on-going active and dust obscured star-forming galaxies.

Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim.

Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objects between redshifts of 0.005 and 0.1. The aim of this classification is to split the galaxy population into merger and non-merger systems and we are currently achieving an accuracy of 92.5%. Stellar masses and star formation rates are also estimated using panchromatic data for the entire galaxy population. With these preliminary data, the mergers are placed onto the full galaxy main sequence, where we find that merging systems lie across the entire star formation rate - stellar mass plane.

Optical properties of infrared-bright (IR-bright) dust-obscured galaxies (DOGs) are reported. DOGs are faint in optical but very bright in mid-IR, which are powered by active star formation (SF) or active galactic nucleus (AGN), or both. The DOGs is a candidate population that are evolving from a gas-rich merger to a quasar. By combining three catalogs of optical (Subaru Hyper Suprime-Cam), near-IR (VIKING), and mid-IR (ALLWISE), we have discovered 571 IR-bright DOGs. Using their spectral energy distributions, we classified the selected DOGs into the SF-dominated DOGs and the AGN-dominated DOGs. We found that the SF-dominated DOGs show a redder optical color than the AGN-dominated DOGs. Interestingly, some DOGs shows extremely blue color in optical (blue-excess DOGs: bluDOGs). A possible origin for this blue excess is either the leaked AGN light or stellar UV light from nuclear starbursts. The BluDOGs may be in the transition phase from obscured AGNs to unobscured AGNs.

To study the dust obscured phase of the galaxy evolution during the peak of the Star Formation Rate (SFR) and the Black Hole Accretion Rate (BHAR) density functions (z = 1–4), rest frame mid-to-far infrared (IR) spectroscopy is needed. At these frequencies, dust extinction is at its minimum and a variety of atomic and molecular transitions, tracing most astrophysical domains, occur. The future IR space telescope mission, SPICA, fully redesigned with its 2.5m mirror cooled down to T < 8K, will be able to perform such observations. With SPICA, we will: 1) obtain a direct spectroscopic measurement of the SFR and of the BHAR histories, 2) measure the evolution of metals and dust to establish the matter cycle in galaxies, 3) uncover the feedback and feeding mechanisms in large samples of distant galaxies, either AGN- or starburst-dominated, reaching lookback times of nearly 12 Gyr. SPICA large-area deep surveys will provide low-resolution, mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, and even the potential to uncover the youngest, most luminous galaxies in the first few hundred million years. In this paper a brief review of the scientific preparatory work that has been done in extragalactic astronomy by the SPICA Consortium will be given.

Spatially resolved studies of galaxies in the high-redshift Universe have traditionally been reliant on data at rest-frame optical and UV wavelengths, which can be biased towards the least dust-obscured galaxies. For several years now, we have been able to resolve and probe the morphology of longer-wavelength emission from distant galaxies with ALMA, and a number of recent ALMA studies were presented at the IAU Symposium No. 352. These included our study of the resolved multi-wavelength emission of galaxies at z ∼ 2. As part of the SHiZELS collaboration, we are mapping the Hα emission line (from SINFONI/VLT), UV continuum (from HST), and the far-infrared (from ALMA) emission from a small sample of Hα-selected galaxies. In this proceedings paper, we showcase the high quality of our data, and the spectacular structures displayed by one of our most dusty sources. We also provide an overview of some highly complementary simulation-based work, using galaxies drawn from the FIRE-2 zoom-in cosmological hydrodynamical simulations. Using sophisticated radiative transfer techniques, we have derived predictions for the spatially-resolved emission of a sample of star-forming galaxies, from rest-frame far-ultraviolet to the far-infrared. For both observed and simulated galaxies, emission maps show striking differences with wavelength, with the same galaxy appearing clumpy and extended in the far-ultraviolet yet compact at far-infrared wavelengths.