Magnetic field plays an important role in star formation and galaxy evolution. Previous studies discussed about the origin of magnetic field and its effect to the environment. With the recent advancement of supercomputers, adding the magnetic field to a cosmological hydrodynamic simulations only become feasible. In this proceeding, we present the results of high-resolution magneto-hydrodynamic simulation with GIZMO and compare our simulation result with the previous literature and the observations.

We construct an X-ray spectral model for the clumpy torus in an active galactic nucleus (AGN), utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al.2016). The geometry of the torus is the same as that in Nenkova et al. (2008), which assumes a power law distribution of clumps in the radial direction and a normal distribution in the elevation direction. We apply our model to the broadband X-ray spectrum of the Circinus galaxy observed with XMM-Newton, Suzaku, and NuSTAR. Our model can well reproduce the observed X-ray spectrum, yielding a hydrogen column density along the line-of-sight ${N_{\rm{H}}^{\rm{LOS}}} = 4.86_{ - 0.04}^{ + 0.07} \times {10^{24}}$ cm−2 and a torus angular width ${sigma = 14.7_{ - 0.39}^{ + 0.44}}$ degree.

The derivation of accurate stellar populations of galaxies is a non-trivial task because of the well-known age-metallicity degeneracy. We aim to break this degeneracy by invoking a chemical evolution model (CEM) for isolated disk galaxy, where its metallicity enrichment history (MEH) is modelled to be tightly linked to its star formation history (SFH). Our CEM has been successfully tested on several local group dwarf galaxies whose SFHs and MEHs have been both independently measured from deep colour-magnitude diagrams of individual stars. By introducing the CEM into the stellar population fitting algorithm as a prior, we expect that the SFH of galaxies could be better constrained.

Early-type galaxies (ETGs) are of crucial importance to trace back the galaxy mass assembly across cosmic time, yet their formation and quenching remain remarkably elusive. The discoveries of massive, dead galaxies at ever-growing redshifts provided compelling evidence to push their formation up to redshift > 4–5 when the Universe was barely 1 Gyr old. In this talk I will present our results on the ages of a new sample of ETGs at z ∼ 3, built by exploiting HST WFC3/G141 rest-frame optical/near-UV grism spectroscopy to study the nature of 10 passive galaxy candidates at 2.5 < z < 3.5 in COSMOS.

This work is part of a PhD project aimed at quantifying the parent space density of distant genuine passive galaxies. I will also discuss the importance of multi-wavelength data in clarifying the degree of contamination by dusty star-forming galaxies affecting the color selection.

By considering a modified version of the evolutionary population synthesis (EPS) model for stellar populations (SPs) comprising binary stars, the retrieved galaxy and HII-region parameters/properties differ from the case of neglecting binary stars. The retrieved age, stellar metallicity and mass of galaxies increase (e.g. ∼ 0.2 dex when using spectral fitting algorithm), whilst the star formation rate decreases (∼0.2 dex). The radiation fields from intermediate-age SPs with binary stars could be potentially important ionizing sources in HII regions. Under this possibility, the theoretical division between star forming galaxy and AGN on the diagnostic diagrams would move towards the up-right corner and the retrieved gaseous metallicity would decrease.

Our prediction for the birth rate of binary neutron stars in SPs ranges from 10−9 to 10−6${\M {^\minus 1_\odot}} $ yr−1 when the kick velocity is from 0 to 190 km s−1.

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.

We highlight the challenges as well as lessons learnt in the derivation of the photometric redshifts for ∼4 million galaxies at 0 < z ≲ 6 contained in the Spitzer Extragalactic Representative Volume Survey (SERVS) and summarise the photometric redshift results recently published in Pforr et al. (2019). The inhomogeneous nature of the ancillary photometry for SERVS presents a similar situation to the one future, large, extragalactic surveys with e.g. LSST and JWST will face. We employ template SED-fitting to determine photometric redshifts. Our comparison of photometric redshifts to ∼75.000 public, spectroscopic redshifts results in an average σNMAD of ∼0.038 and outlier fraction of 3.7% for sources with the best photometric coverage. We find that photometric redshifts are determined most robustly when filter bands are numerous and cover a wide wavelength range. We highlight some possible improvements for the photometric redshifts in SERVS in the future.

The Herschel Extragalactic Legacy Project (HELP) focuses to publish an astronomical multiwavelength catalogue of millions of objects over 1300 deg2 of the Herschel Space Observatory survey fields. Millions of galaxies with ultraviolet–far infrared photometry make HELP a perfect sample for testing spectral energy distribution fitting models, and to prepare tools for next-generation data. In the frame of HELP collaboration we estimated the main physical properties of all galaxies from the HELP database and we checked a new procedure to select peculiar galaxies from large galaxy sample and we investigated the influence of used modules for stellar mass estimation.

The estimation of interstellar dust masses is an important pursuit in our understanding of both local and early Universe – see e.g. the “dust budget crisis”. One of the most used methods of estimating dust masses – dust emission fitting – requires an estimate of the dust opacity at far-infrared and submillimeter wavelengths, but in most models this quantity is based on extrapolation rather than on actual measurements. It is becoming more and more evident that the opacity in typical dust models differs from that of dust analogs measured in the lab, meaning that astronomical dust mass estimations may need to be revised. To estimate the systematic errors introduced by this mismatch, we calculated dust emission for a model where dust far-infrared opacity is the same as that measured in lab samples, then we fit the synthetic emission with a typical (modified blackbody) dust model. Our results show that, if interstellar dust is indeed similar to the lab dust analogs, most fits may overestimate dust masses by as much as an order of magnitude.

We analyse the dust-to-gas mass ratio (DGR) in nearby galaxies on kiloparsec scales. We focus on their dependence on metallicity and the CO-to-H2 conversion factor, αco. We use a sample of 25 nearby galaxies from SINGS and combine our data with CO (2-1) and H I observations from the HERACLES and THINGS surveys. We implement a Hierarchical Bayesian method to derive the dust mass via fitting the infrared data from 100 to 500 μm with a single modified blackbody. We find that the DGR-metallicity relation follows a power law and we study its strong dependency on the conversion factor αco. Our results indicate a strong connection between interstellar dust and gas. The resolved DGR-metallicity relation cannot be represented with a single power law. The scatter in this relation shows the strong impact of several processes that take place in every galaxy.

The infrared (IR) galaxies detected at Herschel/SPIRE 250 μm band over the AKARI’s NEP-Wide field are various types of dusty star-forming (SF) galaxies ranging from quiescent to starbursts, having mid-IR polycyclic aromatic hydrocarbon (PAH) features near 8 μm. The measurements of the 8 μm luminosity (L8μm) along with the total infrared luminosity (LIR) based on the physical modeling of SEDs a take unique advantage of the continuous near- to mid-IR coverage, far-IR data points, and spectroscopically determined accurate redshifts. Our sample shows shortage of 8 μm luminosity compared to the total IR luminosity. This deficit gets severe in more luminous IR galaxies, suggesting PAH molecules in these galaxies are destroyed by a strong radiation field from SF regions, or the existence of a unexpectedly large amount of cold dust in the ISM that contributes to LIR.

We calculate the spectral energy distribution of the first galaxies which contain pre-main-sequence stars by using the stellar evolution code Modules for Experiments in Stellar Astrophysics, the spectra model BT-Settl, and the stellar population synthesis code PEGASE. We calculate the galaxy spectral energy distribution for Salpeter Initial Mass Function. We find that very young first galaxies are bright also in mid-infrared, and the contribution of pre-main-sequence stars can be significant over 0.1 Myr after a star-formation episode.

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.

Interstellar dust is traced by not only thermal emission but also scattered light. The scattered light spectrum observed from ultraviolet (UV) to near-infrared (IR) is useful to constrain some dust properties, such as size distribution, albedo, and composition. Milky Way Galaxy is a unique environment to observe the diffuse scattered light because we can extract it by removing the contribution of starlight. We have observed the UV to near-IR scattered light with space instruments, including Diffuse Infrared Background Experiment (DIRBE), Hubble Space Telescope (HST), and Multi-purpose Infra-Red Imaging System (MIRIS). The scattered light spectrum is marginally consistent with prediction from a recent dust model including carbonaceous and silicate grains with polycyclic aromatic hydrocarbon (PAH). Based on the MIRIS observation of a diffuse cloud, we compare the scattered light color with the dust model with or without grains larger than 1 micrometer. The result shows that the color is consistent with the model without the large grains, which is consistent with recent simulations of dust growth in low-density regions. However, some observations have shown the spectral excess at ∼ 0.6 micrometer wavelength, suggesting the presence of extended red emission (ERE) which cannot be explained by the conventional dust model.

We present a 4.7σ detection of the [OIII] 88 μm line in a gravitationally-lensed Lyman break galaxy, RXC J2248-ID3, using the Atacama Large Millimeter/submillimeter Array (ALMA). We did not detect [CII] 158 μm and rest-frame 90 μm dust continuum emission, suggesting that the bulk of the interstellar medium (ISM) is ionized. Our two-component SED model combining the previous Hubble Space Telescope (HST) data and new photometry obtained from Very Large Telescope (VLT), Spitzer and ALMA suggests the presence of young (∼2 Myr) and mature (∼600 Myr) stellar components with the metallicity of Z = 0.2Z⊙. Our findings are in contrast with previous results claiming a very young, metal-poor stellar component.

Dust attenuation shapes the spectral energy distributions of galaxies and any modelling and fitting procedure of their spectral energy distributions must account for this process. We present results of two recent works dedicated at measuring the dust attenuation curves in star forming galaxies at redshift from 0.5 to 3, by fitting continuum (photometric) and line (spectroscopic) measurements simultaneously with CIGALE using variable attenuation laws based on flexible recipes. Both studies conclude to a large variety of effective attenuation laws with an attenuation law flattening when the obscuration increases. An extra attenuation is found for nebular lines. The comparison with radiative transfer models implies a flattening of the attenuation law up to near infrared wavelengths, which is well reproduced with a power-laws recipe inspired by the Charlot and Fall recipe. Here we propose a global modification of the Calzetti attenuation law to better reproduce the results of radiative transfer models.

Heiiλ1640 emission in the absence of other metal lines is the most sought-after emission line to detect and characterize metal free stellar populations. However, even recent stellar population models with sophisticated treatment of stellar evolution also lack sufficient He+ ionising photons to reproduce observed He 0.1em ii fluxes. We use VLT/MUSE GTO observations to compile a catalogue of 15 z ∼ 2–4 He ii λ1640 emitters from ∼10–30 hour pointings. We show that both He ii λ1640 detections and non-detections occupy similar distribution in UV absolute magnitudes. Rest-UV emission line analysis of our sample shows that the emission lines of our He ii λ1640 emitters are driven by star-formation in solar to moderately sub-solar (∼1/20th) metallicity conditions. However, we find that even after considering effects from binary stars, we are unable to reproduce the He ii λ1640 equivalent widths. Alternative mechanisms are necessary to compensate for the missing He+ ionising photons.

We used the spectroscopic and astrometric data provided from the GALactic Archaeology with HERMES (GALAH) Data Release (DR2) and Gaia DR2, respectively, for a large sample of stars to investigate the behaviour of the [ $\alpha$ /Fe] abundances via two procedures, that is, kinematically and spectroscopically. With the kinematical procedure, we investigated the distribution of the [ $\alpha$ /Fe] abundances into the high-/low-probability thin disc, and high-/low-probability thick-disc populations in terms of total space velocity, [Fe/H] abundance, and age. The high-probability thin-disc stars dominate in all sub-intervals of [ $\alpha$ /Fe], including the rich ones: [ $\alpha$ /Fe] $\,>\,0.3$ dex, where the high-probability thick-disc stars are expected to dominate. This result can be explained by the limiting apparent magnitude of the GALAH DR2 ( $V \lt 14$ mag) and intermediate galactic latitude of the star sample. Stars in the four populations share equivalent [ $\alpha$ /Fe] and [Fe/H] abundances, total space velocities, and ages. Hence, none of these parameters can be used alone for separation of a sample of stars into different populations. High-probability thin-disc stars with abundance $-1.3 \lt {\rm[Fe/H]}\leq -0.5$ dex and age $9 \lt \tau\leq13$ Gyr are assumed to have different birth places relative to the metal-rich and younger ones. With the spectroscopic procedure, we separated the sample stars into $\alpha$ -rich and $\alpha$ -poor categories by means of their ages as well as their [ $\alpha$ /Fe] and [Fe/H] abundances. Stars older than 8 Gyr are richer in [ $\alpha$ /Fe] than the younger ones. We could estimate the abundance [ $\alpha$ /Fe] = 0.14 dex as the boundary separating the $\alpha$ -rich and $\alpha$ -poor sub-samples in the [ $\alpha$ /Fe] $\,\times\,$ [Fe/H] plane.

Although AGN do not typically dominate the bolometric emission of dusty star forming galaxies, large AGN fractions (sometimes > 40%) have been observed in various sub-millimeter surveys. These diagnostics have been however mostly based on X-ray counterpart selections and a complete multiwavength census of the fraction of AGN hosts is needed. I will present new advances in the modelling of panchromatic spectral energy distributions (SEDs) of active galactic nuclei (AGN), based on our publicly available code AGNfitter (Calistro-Rivera et al.2016). AGNfitter implements a fully Bayesian Markov Chain Monte Carlo method to fit the spectral energy distributions of AGNs pushing the wavelengths frontiers from the radio to the X-rays. I will present a recent application of AGNfitter on dusty star forming galaxies in the ALESS submillimeter survey to obtain an unbiased multiwavelength characterisation of the nuclear activity buried in dusty star formation. Our method reveals a significantly larger contribution of AGN activity to the emission in these galaxies than previously observed based on X-rays diagnostics. Our method represents a unique tool to potentially characterise an unbiased accretion history of the Universe when applied to larger populations of star-forming galaxies.

The cornerstone mission of the European Space Agency, Gaia, has revealed properties of 260 000 white dwarfs in the Galaxy, allowing us for the first time to constrain the evolution of white dwarfs with a large sample. Complementary surveys (CoRoT, Kepler, K2, APOGEE and Gaia-ESO), will revolutionize our understanding of the formation and history of our Galaxy, providing accurate stellar masses, radii, ages, distances, and chemical properties for very large samples of stars across different Galactic stellar populations. To exploit the potential of the combination of spectroscopic, seismic and astrometric observations, the population synthesis approach is a very crucial and efficient tool. We develop the Besançon Galaxy model (BGM, Lagarde et al.2017) for which stellar evolution predictions are included, providing the global asteroseismic properties and the surface chemical abundances along the evolution of low- and intermediate-mass stars. For the first time, the BGM can explore the effects of an extra-mixing occurring in red-giant stars Lagarde et al.2019, changing their stellar properties. The next step is to model a consistent treatment of giant stars and their remnants (e.g., white dwarfs). This kind of improvement would help us to constrain stellar and Galactic evolutions.