In this paper, I review several dust evolution studies based on the DustPedia nearby galaxy sample. I first present the dust spectral energy distribution model, implementing a hierarchical Bayesian method, that we have developed. I then discuss the dust evolution trends we have derived among (integrated) and within (resolved) galaxies. In particular, we show that the trend of dust-to-gas ratio with metallicity is clearly non-linear, indicating the need for grain growth in the interstellar medium. Our trend is closer to the one derived with damped Lyα systems than what was suggested by previous studies. We finally demonstrate the universal processing of small amorphous carbon grains by stellar photons.

In Panda et al.2018a, we constructed a refined sample from the original Shen et al.(2011) QSO catalog. Based on our hypothesis — the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature (TBBB) defined by the Big Blue Bump on the Spectral Energy Distribution (Panda et al.2017; Panda et al.2018b). We select the four extreme sources that have RFeII ⩾ 4.0 and use {CIGALE (Boquien et al.2018) to fit their multi—band photometric data. We also perform detailed spectral fitting including the Fe II pseudo—continuum (based on Śniegowska et al.2018)) to estimate and compare the value of RFEII. We show the dependence of FeII strength on changing metallicity.

Dust absorbs stellar emission and reradiates this energy in the far-infrared (FIR). FIR observations hence give us a direct view of the dust, and allow us to study its properties. Unfortunately, FIR observations are only available for a small subset of galaxies. In this work, we estimate the global FIR emission from global UV-NIR observations. We show that a machine learning method clearly outperforms a SED modelling approach. For each galaxy, we not only predict the FIR flux across the 6 Herschel bands, but also estimate individual uncertainties. We inspect the worst predictions, and investigate how the machine learning predictor generalizes on new data. Our predictor can be used as a virtual observatory, which is especially useful now that there is still no confirmed next-generation FIR telescope.

Recently huge amount of dust Mdust ≃ 106−7M⊙ in galaxies at z = 7–8 has been discovered by ALMA observations. The suggested timescale of the dust production was a few–several×108 yr, while the stellar mass was several × 109M⊙. This amount of dust cannot be easily explained only by a supply from supernovae if we consider the dust destruction by reverse shocks. We propose that these values can be consistently explained if we take into account the grain growth in the interstellar medium (ISM). This scenario successfully reproduces the evolution of the dust mass, as well as the SFR, and stellar mass simultaneously. We conclude that even at such an early epoch of the Universe, the dust grain growth in the ISM plays a significant role in galaxies.

Numerical simulations of disk galaxies with steady (long-lived) and dynamic (short-lived) spiral arms suggest that offsets between stellar and gas spiral arms depend on their nature or lifetime (Baba et al.2015). Based on this theoretical study, we investigated gas-star offsets in the nearby grand-design spiral galaxy M51, and found that its two spiral arms exhibit different offset dependences against radius. One arm is consistent with a steady arm, while the other is consistent with a dynamic arm. We deduce that this difference is likely due to a tidal interaction with the companion galaxy (Egusa et al.2017). For this study, a stellar mass distribution with a high accuracy at a high spatial resolution is essential, which has come to be available by applying recent SED fitting techniques to multi-wavelength images. We are now working to extend this study to other nearby spiral galaxies.

The spectral energy distribution (SED) model should treat the evolution of a galaxy from its birth. Dust in galaxies affects the formation and evolution of galaxies in various ways. For example, dust grains scatter and absorb stellar emitted ultraviolet (UV) photons and re-emit the radiation at infrared (IR) wavelengths. In this work, we construct a galaxy SED model based on our dust evolution model (Asano et al. 2013a,b, 2014) with a rigorous treatment of the chemical evolution. To reduce the computational cost, we adopt mega-grain approximation (MGA; (MGA; Inoue, 2005). MGA regards a high density dusty region as a huge size (10 pc) dust grain for calculating dust scattering. In this approximation, we can solve the radiative transfer easily and provide SEDs and attenuation curves of galaxies. This model can be used to fit any galaxy in the wavelength range of 10 nm-3 mm.

I will review the tips learned from panchromatic modeling of active galactic nuclei (AGNs), based on our recent work to study the relationship between AGN and star formation (SF). Several AGN SED models are compared, and a significant AGN contribution is found in the IR luminosities and corresponding star formation rates (SFRs). I will review the AGN-SF relation and how different parameters and sample selections affect the observed correlation. I will then report on the constant ratio discovered between the SFR and the black hole mass accretion rate (BHAR), and their implications on the gas supply and galaxy formation history of these systems. Caveats and important questions to answer are summarized at the end.

Recently, spectroscopic detections of O[III] 88 μm and Ly-α emission lines from the z ≍ 9.1 galaxy MACS1149-JD1 have been presented, and with these, some interesting properties of this galaxy were uncovered. One such property is that MACS1149-JD1 exhibits a significant Balmer break at around rest-frame 4000 Å, which may indicate that the galaxy has experienced large variations in star formation rate prior to z ∼ 9, with a rather long period of low star formation activity. While some simulations predict large variations in star formation activity in high-redshift galaxies, it is unclear whether the simulations can reproduce the kind of variations seen in MACS1149-JD1. Here, we utilize synthetic spectra of simulated galaxies from two simulation suites in order to study to what extent these can accurately reproduce the spectral features (specifically the Balmer break) observed in MACS1149-JD1. We show that while the simulations used in this study produce galaxies with varying star formation histories, galaxies such as MACS1149-JD1 would be very rare in the simulations. In principle, future observations with the James Webb Space Telescope may tell us if MACS1149-JD1 represents something rare, or if such galaxies are more common than predicted by current simulations.

Despite decreasing cosmic star formation rate density over the last 10 Gyr, the stellar mass (M*) buildups in galaxies were still progressing during this epoch. About 50% of the current M* density in the universe was built over the last ∼8.7 Gyr. In this research, we investigated the stellar mass buildup and quenching of spatially resolved regions within massive disk galaxies over the last 10 Gyr. We apply the spectral energy distribution (SED) fitting method to SEDs of sub-galactic regions in galaxies to derive the spatially resolved distributions of SFR and M* in the galaxies. This namely pixel-to-pixel SED fitting method is applied to massive disk galaxies at 0.01 < z < 0.02 and 0.8 < z < 1.8. We found that massive disk galaxies tend to build their M* and quench their star formation progressively from the central region to the outskirts, i.e. inside-out stellar mass buildup and quenching.

Stellar masses are crucial ingredients for putting galaxies in the context of galaxy evolution and are commonly evaluated via Spectral Energy Distribution (SED)-fitting analyses which are hampered by dust attenuation. Observational constraints of attenuation in various galaxy classes provide key inputs for fitting a SED. I will present recent results about the attenuation properties of a sample of Herschel-selected galaxies at 0.7 ⩽ z ⩽ 1.6 widely spanning the star-forming Main Sequence (MS). I will show that far-IR selected galaxies on the MS are well described with local attenuation recipes. Conversely, common recipes cannot recover the SFR of far-IR selected starburst galaxies well above the MS. The SFR of these outliers appears to be hidden by the ∼90% in optically thick cores. These findings pose challenges for SED-fitting codes based on energy balance assumptions that might break in these peculiar sources.

We have developed a new SED fitting tool specialized for frontier redshift galaxies. It is a common case for high-z galaxies that the available data are restricted to rich optical to near-infrared photometry and few far-infrared (FIR) data deep enough to detect the faint object (e.g., HST/WFC3 + Spitzer/IRAC + ALMA). In such situation, one cannot perform a complicated modeling of dust emission in FIR regime. We then adopt simple treatment for the dust emission using empirical LIRG templates. Instead, we adopt a sophisticated and physically motivated modeling for stellar and nebular emission parts in rest-frame UV-to-optical regime. Our new code fits not only broad band photometry but also spectral emission line flux. There is an option to fit observed SED with two templates with different physical properties. Our new code, PANHIT, is now in public, and was already applied to some high-z frontier galaxies.

NGC 300 is a near-optical twin of the Local Group galaxy M33, which are benchmarks for understanding late-type spiral galaxies. They are two bulgeless and low-mass spiral galaxies in different environments. In order to explore the common properties and differences between the two nearby low-mass systems, we first use the simple chemical evolution model to explore the star formation history (SFH) of NGC300 and M33, and then compare the feasible model predicted SFH of NGC 300 with that of M33. Through comparing the SFHs between them, it can be found that the mean stellar age of NGC 300 is older than that of M33, there is a recent lack of primordial gas infall onto the disk of NGC 300, recent star formation along the disk of NGC 300 is less active than that of M33, and the local environment may play a key role in the secular evolution of a galaxy.

Modern cosmological simulations suggest that the hierarchical assembly of dark matter halos provided the gravitational wells that allowed the primordial gases to form stars and galaxies inside them. The first galaxies comprised of the first systems of stars gravitationally bound in dark matter halos are naturally recognized as the building blocks of early Universe. To understand the formation of the first galaxies, we use an adaptive mesh refinement (AMR) cosmological code, Enzo to simulate the formation and evolution of the first galaxies. We first model an isolated galaxy by considering much microphysics such as star formation, stellar feedback, and primordial gas cooling. To examine the effect of Pop III stellar feedback to the first galaxy formation, we adjust the initial temperature, density distribution and metallicity distributions by assuming different IMFs of the first stars. Our results suggest that star formation in the first galaxies is sensitive to the initial conditions of Pop III supernovae and their remnants. Our study can help to correlate the populations of the first stars and supernovae to star formation inside these first galaxies which may be soon observed by the (James Webb Space Telescope JWST).

Using the FirstLight database of 300 zoom-in cosmological simulations we provide rest-frame UV-optical spectral energy distributions of galaxies with complex star-formation histories that are coupled to the non-uniform gas accretion history of galactic halos during cosmic dawn. The population at any redshift is very diverse ranging from starbursts to quiescent galaxies even at a fixed stellar mass. The FirstLight simulations make predictions on the rest-frame UV-optical absolute magnitudes, colors and optical emission lines of galaxies at z = 6–12 that will be observed for the first time with JWST and the next generation of telescopes in the coming decade.

We present a velocity of galactic outflows in star-forming galaxies at the highest redshift, z ∼ 6, so far studied with metal absorption lines. Absorption-line studies of galactic outflows need well-determined redshifts, but there are few strong emission lines in the observed-frame optical spectra of galaxies at high redshifts. In this work, we use the systemic redshifts determined by the ALMA [CII]158 μm emission lines. The sample consists of seven Lyman break galaxies at 5.1 < z < 5.7 whose Keck/DEIMOS and ALMA data are available in the archive. The outflow maximum velocity (νmax) is estimated by a fitting of line profiles to metal absorption lines in a composite spectrum. We find that νmax monotonically increases from z ∼ 0 to 6 and that νmax tightly correlates with the halo circular velocity estimated from the stellar mass.

Panchromatic modeling is one of the most powerful tools at our disposal to measure reliably the physical properties of galaxies across cosmic times. We present here an entirely new implementation in python of one such tool: CIGALE. Developed along three main design principles: simplicity, modularity, and efficiency, it has proven to be a versatile code that in addition to estimating the physical properties of galaxies (or regions within galaxies), can generate arbitrary sets of theoretical models or be used as a library to build other tools. Among its defining features, it is a truly panchromatic code ranging from the far-ultraviolet to the radio that takes into account numerous physical components (including active nuclei or synchrotron emission), that can fit non-photometric data, handle upper limits, determine photometric redshifts, and even build mock catalogs.

Studying the density profiles of galaxy groups offers an important insight on how large-scale structure in the Universe formed and evolved, since galaxy groups bridge the gap between individual galaxies and galaxy clusters. We aim to probe the total density profile of the galaxy group that is gravitational lensing HELMS18, a submillimeter galaxy at z = 2.39 from the Herschel’s HerMES Large Mode Survey (HELMS), by combining strong gravitational lensing with kinematics of the centrally-located galaxies and kinematics of the group members. We will use high-resolution data of HELMS18 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and multi-object spectroscopic data of the group members from Gemini-GMOS. Our final goal is to combine these observations to probe the stellar and dark matter density profiles and to build a complete description of this galaxy group.

Radio continuum emission from galaxies is powered by a combination of distinct physical processes, each providing unique diagnostic information. Over frequencies spanning ∼ 1–120 GHz, radio spectra of star-forming galaxies are primarily comprised of: (1) non-thermal synchrotron emission powered by accelerated cosmic-ray electrons/positrons; (2) free-free emission from young massive star-forming (H ii) regions; (3) anomalous microwave emission, which is a dominant, but completely unconstrained, foreground in cosmic microwave background experiments; and (4) cold, thermal dust emission that accounts for most of the dust and total mass content in the interstellar medium in galaxies. In this proceeding, we discuss these key energetic processes that contribute to the radio emission from star-forming galaxies, with an emphasis on frequencies ≳30 GHz, where current investigations of star formation within nearby galaxies show that the free-free emission begins to dominate over non-thermal synchrotron emission. We also discuss how planned radio facilities that will access these frequencies, such as a next-generation Very Large Array (ngVLA), will be transformative to our understanding of the star formation process in galaxies.

The star formation history (SFH) of galaxies allow us to investigate when galaxies formed their stars and assembled their mass. We can constrain the SFH with high level of precision from galaxies with resolved stellar populations, since we are able to discriminate between stars of different ages from the spectrum they emit. However, the relative importance of secular evolution (nature) over nurture is not yet clear, and separating the effects of interaction-driven evolution in the observed galaxy properties is not trivial. The aim of this study is to use MaNGA (Mapping Nearby Galaxies at APO) Integral Field Unit (IFU) data, in combination with multi-wavelength data, to constrain the SFH of nearby isolated galaxies. We present here the new techniques we are developing to constrain the SFH with high level of precision from Spectral Energy Distribution (SED) fitting. This study is part of a China-Chile collaboration program where we are applying these new techniques to investigate how galaxies formed and evolve in different environments.

We build a theoretical picture of how the light from galaxies evolves across cosmic time. In particular, we predict the evolution of the galaxy spectral energy distribution (SED) by carefully integrating the star formation and metal enrichment histories of semi-analytic model (SAM) galaxies and combining these with stellar population synthesis models which we call mentari. Our SAM combines prescriptions to model the interplay between gas accretion, star formation, feedback process, and chemical enrichment in galaxy evolution. From this, the SED of any simulated galaxy at any point in its history can be constructed and compared with telescope data to reverse engineer the various physical processes that may have led to a particular set of observations. The synthetic SEDs of millions of simulated galaxies from mentari can cover wavelengths from the far UV to infrared, and thus can tell a near complete story of the history of galaxy evolution.