Synthetic Galaxy Images and Spectra from the Illustris Simulation
We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.
keywords:methods: numerical – cosmology: theory – cosmology: galaxy formation
The formation of dark matter haloes has been studied extensively using numerical dark matter only simulations (e.g., Springel et al., 2005b; Boylan-Kolchin et al., 2009; Fosalba et al., 2008; Teyssier et al., 2009; Klypin et al., 2011). Extending the insight from dark matter only simulations to include a theory of galaxy formation requires a method to link the formation of dark matter haloes to observable galaxy properties. The most direct method available is hydrodynamical simulations, which model the co-evolution of dark matter and baryons (e.g. Katz et al., 1992, 1996; Weinberg et al., 1997; Murali et al., 2002; Springel & Hernquist, 2003b; Kereš et al., 2005; Ocvirk et al., 2008; Crain et al., 2009; Croft et al., 2009; Schaye et al., 2010; Oppenheimer et al., 2010; Vogelsberger et al., 2012). By directly including hydrodynamics in structure formation simulations one can probe the thermal state and column density distribution of the intergalactic medium (IGM) via the Lyman- forest (Cen et al., 1994; Zhang et al., 1995; Hernquist et al., 1996; Theuns et al., 1998), the phase structure and heavy element composition of the circumgalactic medium (CGM) and IGM (Aguirre et al., 2001; Cen & Fang, 2006; van de Voort & Schaye, 2012; Shen et al., 2013), and the nature and rates of accretion of dark matter and baryons into galaxies (Kereš et al., 2005; van de Voort et al., 2011; Nelson et al., 2013). Since these simulations can follow the dynamics of both the dark matter and baryons down to small spatial scales, predictions can be made about their internal structure including the distribution of gas (Kereš et al., 2012; Torrey et al., 2012) and the formation of stellar disks and bulges (Abadi et al., 2003; Governato et al., 2004; Agertz et al., 2011; Sales et al., 2012; Marinacci et al., 2014). As a result of this detailed information, simulations are a valuable tool for interpreting observational data and placing observed galaxies into a more complete evolution based cosmological context.
Many important galaxy properties – such as stellar mass – have substantial uncertainties associated with their measurement (e.g., Conroy & Wechsler, 2009; Conroy & Gunn, 2010, and references therein). These uncertainties can be physical in origin (e.g., the slope of the IMF, treatment of the thermally pulsating asymptotic giant branch phase, etc.) or due to assumptions made during the spectral energy distribution (SED) fitting procedures (Papovich et al., 2001; Wuyts et al., 2007; Gallazzi & Bell, 2009; Michałowski et al., 2012; Banerji et al., 2013). Accurately estimating physical galaxy properties based on broadband photometric data points (or even full SEDs) is difficult because observed galaxy SEDs contain contributions from stars with a complex distribution of stellar ages and metallicities. Simplifying assumptions regarding the functional form of the star formation history, or a uniform stellar metallicity, are often assumed. These assumptions can lead to systematic errors on the derived galaxy stellar mass, independent of physical uncertainties (e.g., Mitchell et al., 2013). Similarly, derived galaxy structural properties – such as bulge-to-disk decompositions or galaxy sizes – are often determined inconsistently in simulations and observations (Scannapieco et al., 2010).
A crucial difficulty when relating hydrodynamical simulations and observations is the translation of results between the two, a necessary step for accurate comparisons (e.g., Conroy et al., 2010; Hayward et al., 2012, 2013a, 2013b). There are two fundamental ways to make this translation, by either converting observed quantities into the space of physical parameters, or by generating ‘mock’ (or synthetic) observations of the simulation output. Only the second approach provides a unique mapping, avoiding the uncertainties involved in a potentially non-unique inversion. It also allows observational effects to be accounted for (e.g. noise, survey geometry, …), and enables observational tools and data analysis techniques to be run on the simulated observations as is. For comparison to observations of stellar light, this can be achieved in simulations by adopting stellar population synthesis (SPS) models to assign light to all star particles within a galaxy based on their age and metallicity distributions (Tinsley, 1972; Bruzual A., 1983; Buzzoni, 1989; Bruzual A. & Charlot, 1993; Worthey, 1994; Maraston, 1998; Leitherer et al., 1999; Bruzual & Charlot, 2003; Thomas et al., 2003; Maraston, 2005).
These mock observations can be used to address potential issues in the way that SED fitting or photometric analysis is performed and even to directly fit galaxy SEDs rather than relying on SED modelling (Lanz et al., 2014). A large body of literature already exists within the idealized merger simulation community in this direction: Lotz et al. (2008) used mock broadband images of idealized merger simulations (Cox et al., 2006) to determine the conversion from observed galaxy pair counts into merger rates. Wild et al. (2009) and Snyder et al. (2011) used mock galaxy SEDs of idealized mergers to determine the lifetime and redshift dependent abundance of K+A galaxies, while Snyder et al. (2013) used synthetic galaxy spectra to propose a new Mid-IR diagnostic of AGN. Wuyts et al. (2010) investigated how well the true mass profiles of simulated massive compact galaxies could be recovered from their light profiles in different bands. Hayward et al. (2011) elucidated the physical meaning of the submillimetre galaxy selection, and Hayward et al. (2014b) tested how well the integrated IR luminosity traces the star formation rate.
Similar work exists within the semi analytic community including the production of mock galaxy luminosity functions (e.g., Henriques et al., 2011; Somerville et al., 2012) and mock light-cone data (Kitzbichler & White, 2007; Henriques et al., 2012; Overzier et al., 2013; Bernyk et al., 2014) based on dark matter only simulations (Springel et al., 2005b). The Millennium Run Observatory (MRObs, Overzier et al., 2013) built a theoretical virtual observatory capable of producing synthetic data products that could then be reduced and analyzed with observational software. The MRObs synthetic observations covered forty filters where the light was assigned to each galaxy by combining information about the star formation history with stellar population synthesis models. The resulting data products were formatted in a way that they could be analyzed using standard observational software packages (Overzier et al., 2013) and distributed using a web based relational SQL database (Lemson & Springel, 2006).
In this paper, we describe the production of a mock galaxy image catalog based on the recent Illustris simulation (Vogelsberger et al., 2014b, a; Genel et al., 2014). The Illustris project has two basic principles: (i) build a comprehensive physical framework that allows galaxy populations to regulate their stellar mass growth appropriately and (ii) apply that physical framework in hydrodynamical simulations of large cosmological volumes. A detailed description of our adopted feedback module, which includes a model for star formation driven winds as well as AGN driven outflows, was presented in Vogelsberger et al. (2013) and demonstrated to appropriately regulate galaxy populations across cosmic time in Torrey et al. (2014). Adopting this feedback model ensures that the galaxies within these simulations match a wide range of galaxy observable properties such as the redshift galaxy stellar mass function. The main Illustris simulation was run in a periodic box of size Mpc with sufficiently high mass and spatial resolution (i.e. and kpc) to resolve a wide range of galaxy structures (Vogelsberger et al., 2014b, a; Genel et al., 2014).
Using results from the Illustris simulation, we employ stellar population synthesis (SPS) templates to produce synthetic broadband galaxy images and integrated SEDs. These mock images encode information about the spatial distribution and formation history of galactic stellar components. The synthetic galaxy SEDs and images produced from the Illustris simulation are unique compared to previous synthetic data catalogs. Compared against idealized isolated and merging galaxy simulations (e.g., Lotz et al., 2008; Snyder et al., 2011), the Illustris simulation galaxies are embedded in a proper cosmological context. The galaxies in the Illustris simulation accrete gas from the surrounding intergalactic medium, undergo major and minor mergers, and experience feedback as described in Vogelsberger et al. (2013). All of these processes are reflected in the stellar mass distribution and star formation history, which define the resulting mock images. Previous generations of synthetic images and spectra have been generated from “zoom-in” style simulations of individual (or a small number of) objects (e.g., Guedes et al., 2011; Stinson et al., 2013; Snyder et al., 2014). Like the Illustris simulation, these simulations contain a proper cosmological context for galaxy growth, and they have there been used to study the growth of galactic disks and bulges. Our large-volume simulation approach is distinct from the “zoom-in” style simulations in that it allows us to analyze the properties of a large population of galaxies with varied formation histories and environments under uniform assumptions about the input physics (e.g., Pedrosa et al., 2014). In contrast to semi-analytic based mock data products (e.g., Overzier et al., 2013), our hydrodynamic approach allows us to self-consistently treat the internal galactic dynamics (e.g., bar formation, merger events, etc.) that drive galaxy morphology evolution. The resulting synthetic image catalog has many applications including building mass dependent average galaxy SED templates, isolating the origin of shelled galaxies, and understanding the viewing angle dependence of barred galaxies.
In this paper, our main goal is to detail the methods for producing these images, explain the data format, and explain how the image catalog will be accessible. We additionally demonstrate some initial science applications of the resulting synthetic images and spectra. This paper is outlined as follows. In section 2 we describe the simulation that has been used to produce our mock image catalog as well as our methods for producing the mock galaxy image catalog based on its output. In section 3 we describe the format for the mock galaxy images, present properties of the galaxy SEDs, and demonstrate multi-band image generation. In section 4 we provide a detailed application of this data, by determining stellar mass measurements using a popular SED fitting code to compare the SED derived stellar masses against the simulation based stellar mass determination. In section 5 and section 6 we summarize and conclude.
Preliminary JWST/NIRCAM filter curves were obtained from www.stsci.edu/jwst/instruments/nircam/instrumentdesign/filters
In this section we describe our procedure for producing mock galaxy SEDs and photometric images. We provide details on the simulation that forms the basis of our mock image catalog, explain how light is assigned to each galaxy to produce mock images, and outline the resulting data product format.
2.1 The Illustris Simulation
All mock images are produced based on galaxies formed in the Illustris simulation. A detailed description of the Illustris simulation can be found in Vogelsberger et al. (2014b), Vogelsberger et al. (2014a), and Genel et al. (2014). We briefly summarize the simulation properties below.
Illustris is a large-volume cosmological hydrodynamical simulation that was run in a periodic box of side length Mpc with dark matter particles () and baryon resolution elements (). Cosmological parameters consistent with the latest Wilkinson Microwave Anisotropy Probe (WMAP)-9 measurements were adopted (, , , , , and km s Mpc with ). The Illustris simulation was run using the moving mesh code AREPO which, in addition to gravity and hydrodynamics (Springel, 2010), includes comprehensive physics and feedback modules allowing for radiative gas cooling (Katz et al., 1996; Wiersma et al., 2009a), heating and ionisation by a UV background (Faucher-Giguère et al., 2008, 2009; McQuinn et al., 2009), star formation with associated feedback (Springel & Hernquist, 2003a), mass and metal return to the interstellar medium from aging stellar populations (Wiersma et al., 2009b), and active galactic nuclei (AGN) feedback (Springel et al., 2005a; Sijacki & Springel, 2006; Sijacki et al., 2007; Hayward et al., 2014a) as described in detail in Vogelsberger et al. (2013). In practice, the feedback included in the Illustris simulation was tuned to match the redshift galaxy stellar mass function and evolving cosmic star formation rate density (Vogelsberger et al., 2013), and shown to produce galaxy populations that evolve consistently with observations (Torrey et al., 2014).
2.2 Galaxy Image Production from Illustris Data
The main function of the image pipeline is to produce synthetic images and spectra for galaxies in the Illustris simulation. Haloes are defined in the Illustris simulation via a Friends-of-Friends (FoF) algorithm while galaxies are found via the SUBFIND halo finder (Springel et al., 2001). For each halo and galaxy we identify all associated cells and particles (gas, stars, dark matter, and black holes) and produce a ‘mini-snapshot’ file that contains all relevant particle information. These mini-snapshots files are in principle a redundant replication of the simulation data, however they are produced to make the data for the galaxies more easily accessible (each full Illustris snapshot is 1̃.4TB, and therefore it can be rather inefficient to extract individual galaxy data). We produce two types of mini-snapshot files including either the full FoF groups or only individual subhalos. Images and spectra made from the FoF mini-snapshots will contain contributions from neighboring and satellite galaxies. This can be an advantage when, e.g., visually identifying merger events, but is not favorable when, e.g., trying to derive a galaxy’s mass based on its SED (if it contains ”extra light” contributions from neighboring galaxies). These mini-snapshots can be opened and processed very efficiently and are in the natural format for simulation post-processing analysis tools that have been designed to run on individual galaxies.
assign full spectra to each star particle based on its age and metallicity (see Section 2.2.1);
generate images of arbitrary field of view and pixel size for several camera positions (see Section 2.2.4);
convolve spectra with pre-tabulated broadband filter transmission functions (see Section 2.2.4 and Table 1).
We perform these tasks without applying the dust absorption, scattering, or emission functionalities of the SUNRISE code. While these are important considerations for SED modeling, we do not expect the SUNRISE radiative transfer results to be fully converged for the mass and spatial resolutions employed by the Illustris simulation (Jonsson et al., 2010). Instead, we use only a simple empirically motivated dust obscuration correction (described below) and leave a more complete examination of dust and its consequence for mock images to future studies. Using SUNRISE allows us to rely upon an existing, tested, publicly available code, and to prepare our pipeline software for use on future hydro simulations where full dust modeling may be appropriate. Although AGN are present in the simulation, we do not include their light contributions.
We post-process all galaxies with more than 500 stellar particles – roughly corresponding to stellar masses of , and run SUNRISE in two modes depending on the particle count of the system. For systems with stellar masses ( stellar particles) we calculate both integrated SEDs and spatially resolved photometric maps in 36 broadband filters (see Table 1 and Section 2.2.4). The simulation volume contains 7,000 galaxies above this limit at redshift . Systems with stellar masses below this limit are not very well spatially resolved – making imaging less useful – so we store only integrated SEDs. The integrated SEDs do not allow us to study the internal structure of galaxies, but still allow us to identify the impact of metallicity and star formation history on galaxy SEDs which we use in Section 4.1.
When a star particle is formed in the simulation it inherits the metallicity of the local ISM gas and records its birth time and initial formation mass. Star particles then slowly return mass to the ISM to account for mass loss from aging stellar populations. We track both the formation mass (a fixed value), as well as the current mass (a time dependent value) for each star particle in the simulation. For the galaxies, light is then assigned to all constituent star particles based on their mass, age, and metallicity values, using single-age stellar population SED templates. We adopt the STARBURST99 (SB99) single-age stellar population models which includes an SED grid resolved with 240 wavelength bins at 308 age and 25 metallicity values (Leitherer et al., 1999; Vázquez & Leitherer, 2005; Leitherer et al., 2010), from which a single SED is assigned to each star particle, scaled according to its formation mass. The fiducial SB99 input SEDs also contain nebular emission and absorption lines according to the model outlined by Leitherer et al. (1999). Figure 1 shows the age dependent SEDs for three representative stellar metallicity values in the SB99 templates.
In our current simulations, the star particles have an initial stellar mass of and represent an unresolved full stellar population. As a result, it is not immediately obvious how the associated stellar light should be spatially distributed. Adopting an appropriate scheme for spreading stellar light is unimportant for integrated quantities (e.g., the integrated galaxy SED) when radiative transfer is not performed but can impact visual interpretation of galactic structure and bias quantitative measures. Over smoothing of stellar light will lead to an undesirable modification or deformation of, e.g., galaxy light profiles, whereas under smoothing light can have adverse consequences for the resulting galaxy images. For example, assigning all stellar light to a single pixel (i.e. assuming each simulation star particle is best represented by a point source) will result in a population of easily identified bright pixels in the galactic outskirts. This can bias our mock galaxy images toward easy detection of diffuse light components which may in fact have very low surface brightness. When using stellar light point sources the peak surface brightness and distribution of surface brightnesses for pixels in an image will be unphysically set by the adopted pixel scale and/or number. Applying a finite spatial extent to the light assigned from the star particles ensures that the surface brightness will not be dependent on the image resolution as long as the pixel scale is smaller than the smoothing scale for the star particles.
We consider six choices for spreading the stellar light as described below, and show an example of each in Figure 2. The top left panel of Figure 2 assigns the light from the star particles as a point source, a method which requires no additional assumptions, but leads to pixelated images where all light from any star particle is always contained within a single pixel. Image pixelation becomes particularly noticeable in the outskirts of the galaxy, and in general is not realistic since each star particle represents a population of stars, with finite spatial extent. Artificial pixelation can be reduced by spreading the light from the star particles following some radial distribution, e.g. a Gaussian kernel. This then requires a smoothing length, for which one natural choice is the gravitational softening length – the same example image resulting from this choice is shown in the middle panel of Figure 2. Employing this fixed smoothing length results in structure in the center of the galaxy being under-resolved, while the galaxy outskirts are represented by isolated, circular stellar light distributions. Increasing the fixed smoothing length further, as shown in the top right panel of Figure 2, leads to a smooth distribution of light in the galaxy outskirts, although with ever greater loss of visible galaxy structure, as one might expect.
As an alternative, we find that adaptive light spreading based on the local density of star particles yields a smooth and well-resolved distribution of light from the inner to outer parts of the galaxy. In practice we set the Gaussian FWHM for each star equal to the distance enclosing its nearest neighbor (in 3D). The level of smoothing is then set based on the neighbor number, , that we adopt. The bottom three panels of Figure 2 show examples of this adaptive light smoothing for three choices of the neighbor number (). We have modified SUNRISE to include the same tree build and neighbor finding steps included in GADGET-2 (Springel, 2005) to facilitate efficient neighbor finding for all systems.
None of the light distribution options are formally any more or less valid than the others, and we cannot rule out that a more physical light assignment procedure may exist. In all subsequent analysis in this paper we adopt as our fiducial choice for the adaptive light smoothing.
We use photon packets to perform the SUNRISE Monte Carlo photon propagation scheme. This results in some residual Monte Carlo noise in low surface brightness features that manifests as fluctuations in the pixel-to-pixel intensity. These low surface brightness features are in general poorly resolved and the Monte Carlo noise should not be mistaken for detailed information about the stellar mass/light distribution. We have tested and confirmed that photon packets is sufficient such that all of the results shown in this paper are well-converged and would not change significantly with increased photon packet count. However, caution should still be exercised when examining the detailed structure of low surface brightness features.
A proper treatment of dust and its modification of simulated galaxy spectra requires knowledge of the distribution of gas and dust on small ( 1 kpc) spatial scales. To most accurately determine the effects of dust attenuation, full radiative transfer should be performed on the simulated galaxies, but the spatial resolution of state-of-the-art cosmological simulations is not sufficient for such calculations (Jonsson et al., 2010). In principle, one can estimate the impact of dust obscuration by calculating the column density of dust (via the gas density, temperature, and metallicity spatial distribution), and thus the optical depth, along lines of sight to each star in the galaxies (e.g., Hopkins et al., 2005; Robertson et al., 2007; Wuyts et al., 2009a, b). However, in the Illustris simulation we are limited to kpc scale resolution, with a pressurized, effective equation-of-state treatment of the complex structure of the ISM (Springel & Hernquist, 2003a). As a result, the gas disks in our simulated galaxies are thicker (in terms of disk scale height) and more homogeneous than the ISM in real galaxies. A naive calculation of the optical depths to stars within galaxies based on the simulated gas distribution will therefore not necessarily lead to accurate dust obscuration measurements. In principle, a more robust dust correction can be constructed via a sub-grid dust-obscuration model that accounts for the unresolved properties of the gas disks. Such a model is beyond the scope of this paper, but is being considered in Snyder et al., (in prep).
Here, we instead consider the simple empirical dust model of Charlot & Fall (2000, hereafter CF00). The CF00 model accounts for the effects of dust by applying an effective absorption of to the standard SB99 SED. The normalisation of the absorption curve is lowered by a factor of three for stars older than yrs, which are assumed to have left their birth clouds. The only parameter in this model is stellar age; no knowledge of the local ISM conditions (including the local gas density or metallicity) is taken into account. Since most galaxies have only a very low fraction of star particles with ages yrs, applying this CF00 model typically results in a more-or-less uniform screen being placed in front of the fiducial galaxy image. As a result, the CF00 models have a spatial distribution of light which is similar to the fiducial non-dusty images (albeit, with an intensity offset determined by the CF00 attenuation curve).
To appropriately include the impact of the unresolved birth cloud on the emergent SED from young star particles, we adopt a model that takes into account nebular emission around young stars. This is potentially important because line emission from gas in HII regions (e.g., H) can contribute substantially to the flux in certain broadband filters. Line contamination in broadband filters can modify inferred ages and stellar masses, especially for high redshift galaxies where the average star formation rates and nebular emission line emission are much higher than in the local universe (Zackrisson et al., 2008; Schaerer & de Barros, 2009; Ono et al., 2010; Stark et al., 2013). We calculate SEDs for the galaxy population including a model (for a complete description see Jonsson et al., 2010) that includes the environment around young stellar populations. To account for this in our mock galaxy SEDs, SUNRISE reassigns the intrinsic SEDs of young star particles assuming a partially obscured young stellar spectra with added contributions from HII regions (Dopita et al., 2005, 2006a, 2006b) as calculated in Groves et al. (2008) via the MAPPINGS-III 1-D photoionisation and radiative transfer code (Binette et al., 1985; Sutherland & Dopita, 1993; Groves et al., 2004a, b). We assume a photodissociation region covering fraction of and cluster masses of . The MAPPINGS model replaces SB99 emission for all young ( yrs) star particles. As in the case of the CF00 model, the spatial distribution of light for our simulated galaxies with and without the MAPPINGS model is similar and so we store only integrated quantities for this model. In summary, our treatment models the partial obscuration of young stars by their birth clouds, where some of the obscured light is re-emitted in nebular emission lines.
Image Production and Broadband Definitions
Images are produced by projecting the stellar light using pinhole cameras with a resolution of pixels. We use a fiducial value of for the full image catalog. However, we have used higher pixel resolution () for many of the individual images included in this paper. Each galaxy is imaged with four cameras, which see the galaxy from four different viewing angles. The cameras sit at the vertices of a tetrahedron and are pointed toward the center of the tetrahedron. We move the galaxy’s position such that the galaxy’s potential minimum sits at the center of the tetrahedron, which results in the four cameras being randomly aligned with respect to the rotation axis of each galaxy. The first camera (CAMERA0, in our later notation) corresponds to the positive direction. Each camera is placed 50 Mpc away from the galaxy’s center, regardless of the galaxy’s intrinsic redshift, and the image field-of-view is set to be 10 times the stellar half mass radius for the galaxy.
The first imaging step is the production of a mock integral field unit (IFU) data cube – with full spectra being associated with each pixel. Integrated global SEDs are produced by adding the light contributions from all IFU pixels. As such, we note that any light contributions that fall outside of the initially selected field of view (10 times the stellar half mass radius) will not be included. We note that we have examined global SEDs computed using larger fields of view, and found them to be nearly indistinguishable from those produced using our fiducial field of view choice. We convolve the IFU mock data cubes with 36 broadband filter transmission functions – as listed in Table 1 – to produce an equal number of broadband images for the galaxy. The resulting broadband images can be used to generate RGB images based on various filter combinations, and to facilitate quantitative comparisons with observations made with particular filters.
The images generated from this process are idealized images, in the sense that they contain no sky noise, no camera point-spread-function (PSF) blurring, and lack galactic and extragalactic image contamination (from the foreground/background). We do not include these “image realism” contributions because (i) they can be added to the idealized images in “post processing” and (ii) the best choice for the, e.g., camera PSF blurring will depend on the specific instrument and/or application. In future work we will include such non-intrinsic contributions before analyzing the photometrically inferred galaxy properties.
The post processing of every galaxy takes roughly 30 minutes running on 4 processors – leading to a total compute time of roughly 80,000 CPU hours for the redshift population.
3 Results: Galaxy Images and SEDs
The redshift galaxy population has been run through the SUNRISE pipeline producing a catalog of 40,000 galaxy spectra and 7,000 galaxy images with associated photometric data.
This data is available through a web-based query and download service.
Photometric data for a galaxy is contained in a single, multi-extension FITS file containing four image extensions corresponding to the four camera views. Each extension contains an array of dimensionality . The ordering of the broadband filters can be found in Table 1 as well as in the first field in the FILTERS extension. Image units, pixel count, and other relevant details are stored in the camera extension header as outlined in Table 2. This initial release closely follows the SUNRISE output format, for which additional information and documentation was presented by Jonsson (2006); Jonsson et al. (2010).
|cameradist||Distance from camera to galaxy [kpc]|
|linear_fov||Image FOV [kpc]|
|CAMPOS[X,Y,Z]||Camera Position [kpc]|
|CAMDIR[X,Y,Z]||Camera Viewing Direction|
|CAMUP[X,Y,Z]||Camera Up Direction|
|TTYPE[1-12]||Names of included data fields|
|TUNIT[1-12]||Units for included data fields|
|L_bol_grid||Bolometric Luminosity [W]|
|lambda||SED array [m]|
|L_lambda||SED array [W/m]|
|TTYPE[1-20]||Names of included data fields|
|TUNIT[1-20]||Units for included data fields|
|Filter Name List||String array giving the broadband filter names|
|Effective Wavelengths||for each filter [m]|
|AB_mag_nonscatter[0,1,2,3]||Absolute magnitudes in each band|
|IMUNIT||Units for included images|
|data||36 x N x N array; Contains all broadband images|
3.1 Integrated SEDs
We store integrated global SEDs for each galaxy under the IntegratedQuantities field. Example SEDs are shown in Figure 3 for one star-forming blue galaxy (blue line, ID=350781), and one quiescent red galaxy (red line, ID=138415). Both galaxies have stellar masses of . The galaxies have star formation rates of yr (blue) and yr (red). While the blue galaxy still has a reservoir of star-forming gas, the red galaxy has quenched its star formation via AGN feedback. The two galaxies have noticeably different SED slopes with the star-forming galaxy having substantially more UV emission. Moreover, since a somewhat significant amount of light in the near IR comes from young stellar populations, an offset in the near IR SED can also be seen.
The fiducial SEDs do not take into account dust, which implies that the UV luminosities are overestimated due to missing SED attenuation. This is demonstrated in Figure 3, where the cyan and magenta lines show the star-forming galaxy’s SED when applying the MAPPINGS HII region model and the CF00 empirical dust correction, respectively. In both cases, these models reduce the UV luminosity by approximating the obscuration of young stellar spectra. At optical wavelengths, the CF00 model depresses the SED by , while the MAPPINGS model is nearly coincident with the pure SB99 spectra. However, since both galaxies are fairly massive and dominated by older stellar populations which contribute to the near IR light, the offset in the near IR portion of the SED is significantly less noticeable than in the UV.
To present a more complete overview of the synthetic spectra properties, we show all 40,000 Illustris galaxies by binning in narrow (0.01 dex) stellar mass bins and calculate the median galaxy spectra at each wavelength as a function of galaxy stellar mass in Figure 4. This figure shows the mean galaxy SED template that applies to our redshift synthetic galaxy spectra as a function of stellar mass. The intensity of the median NIR SED increases smoothly with increasing galaxy stellar mass. Little UV emission from both low- and high-mass galaxies can be identified. The least massive galaxies produce little UV emission owing to their low star formation rates derived from their position on the star formation main sequence. The high mass galaxies have low star formation rates owning to quenching from AGN feedback in our model. This view of the data also allows us to pick out several features of the mock spectra easily. Examples include the Lyman limit at 912 Å, Balmer Break at 4000 Å, and H emission line at 6563 Å– which are all included in the input SB99 template SEDs. Departures from the median trends identified in Figure 4 can be substantial in the UV part of the spectrum owing to fluctuations in the star formation rates of galaxies – even at similar stellar masses. This variability is pronounced for high-mass systems (i.e. ) where quenched and non-quenched galaxies sit in the same stellar mass bins. On the other hand, the variability of the median spectra are low for the most massive systems in the NIR region owing to their NIR luminosities being dominated by old stellar populations.
The median galaxy SED templates for the Illustris galaxy population presented in Figure 4 do somewhat depend on the adopted choice for the stellar population synthesis model. We have compared the findings of Figure 4 with the BC03-based outputs and found that qualitatively the mass binned median galaxy spectra generated using the BC03 stellar population synthesis models look very similar to those generated using the SB99 procedure. The main distinction between these models can be seen in their UV spectra. Whereas the SB99 models have little or no UV emission from old stellar populations, the BC03 catalogs contain non-negligible UV emission from old stars. This causes the UV brightness of massive (quenched) galaxies to be substantially brighter for the BC03 model. In principle, this difference also contributes to the UV emission of intermediate mass galaxies (i.e. ). However, since the intermediate mass galaxies have substantial ongoing star formation dominating their UV spectra, any contribution to the UV spectra from old stellar populations is not easily identified.
3.2 Galaxy Images
Spatially resolved galaxy images are stored in 36 bands, from which surface brightness maps and color images can be directly produced. Color galaxy images produced from our image pipeline were first presented in Figure 13 of Vogelsberger et al. (2014a) which shows 42 blue star forming galaxies and 42 red, passive galaxies as viewed through the SDSS-g, -r, -i bands using the Lupton et al. (2004) asinh scaling. Here we present additional applications of the broadband images by demonstrating: (1) the viewing angle dependence of a select set of galaxy types, (2) the visual appearance of a simulated galaxy in different filters, and (3) the time evolution of several elliptical galaxies exhibiting stellar shells.
Multiple Galaxy Viewing Angles
Figure 5 shows RGB color images of two disk galaxies made using the SDSS-g, -r, -i bands as seen from three different directions to demonstrate the impact of viewing angle on perceived galactic morphology. The left (right) most image is aligned (anti-aligned) with the angular momentum vector for the stars in the galaxy. The central images uniformly transition between these two polar views. We note that the full galaxy image catalog contains only four viewing angles for a galaxy, without a guarantee of any of the images being exactly (anti-)aligned with the rotation axes. The images presented here have been constructed after manually rotating the orientation of each galaxy to highlight the viewing angle dependence. While the top galaxy (ID=283832) is relatively axisymmetric, the bottom galaxy has a prominent bar (ID=261085). The observable strength of the bar feature is a function of viewing angle, with the bar appearing strongest when viewed face-on and being completely unidentifiable when the galaxy is seen edge-on. These mock images could then be used to derive quantitative correction factors for visually classified bar strengths, or to determine what fraction of strongly barred systems would have been classified as having weak bars, or no bars altogether, based on viewing angle.
Figure 6 shows an example galaxy as seen through 12 different broadband filters, as labeled within the figure. Each image has been scaled independently using an “asinh” scaling, with the non-linear transition set to the mean pixel intensity for that band. Bands increase in wavelength from top left to lower right (effective wavelengths associated with the bands can be found in Table 1). The UV bands reveal the emission from young star particles – tracing the location of ongoing star formation in this system. Longer-wavelength bands trace stellar mass more efficiently and exhibit less band-to-band variation. We note that this particular galaxy (ID=283832) was selected because it possesses a large stellar disk (), with moderate levels of ongoing star formation, enabling us to make a clear demonstration of the wavelength-dependent image variation. This is the largest and most massive disk galaxy in the simulation volume; only a small fraction of the galaxies in the simulation are as well-resolved as the one presented here.
Figure 7 shows a time series of RGB images from the rest frame SDSS-g, -r, -i bands for three particular galaxies that we visually identify to possess stellar shells at redshift . We trace these systems back in time using merger trees, identifying the most massive progenitor galaxy at each previous snapshot. We generate images for a system at every snapshot using a fixed field of view of 120 physical kpc. We find that merger events occur for all three systems within the past Gyr, which are then responsible for forming the observed present-day stellar shells. This is consistent with the established theory of stellar shell formation via major (Hernquist & Spergel, 1992) and minor (Quinn, 1984; Dupraz & Combes, 1986; Hernquist & Quinn, 1988; Kojima & Noguchi, 1997) mergers.
To identify the shell locations, we make monochromatic images of these systems in the Johnson-K band (as a good representation of the stellar mass) and apply an unsharp mask procedure. In detail, we convolve the original image with a two dimensional gaussian with a standard deviation of 5 pixels. We then subtract the gaussian blurred image from the original image, and discard all negative pixels to specifically highlight stellar light excesses relative to the blurred image. The resulting masked image is shown in Figure 8 using a logarithmic stretch. All high frequency image contributions – including noise from the discrete particle representation of stellar mass – are amplified in the masked image. As a result, the shells are more easily distinguished from the smoothly varying light component. Using the masked images, we identify the stellar shells down to fairly low galactocentric radii with 3-5 shells being visually apparent in the systems at redshift z=0. Stellar shells have been observed at low galactocentric radius using similar image masking techniques (see Canalizo et al., 2007, for a detailed description). We have confirmed the presence of low radii stellar shells by examination of a radial velocity versus radius phase diagram.
The presence of multiple systems with layered stellar shells in our cosmological simulation is of interest since the majority of previous theoretical work studying the evolution and characteristics of stellar shells formed via mergers has been done either with idealized merger simulations (Quinn, 1984; Hernquist & Spergel, 1992; Dupraz & Combes, 1986; Hernquist & Quinn, 1988; Kojima & Noguchi, 1997; Canalizo et al., 2007) or with dark matter only cosmological simulations (e.g., Cooper et al., 2011). Since these systems are realized here within the full cosmological context, we can address some previously inaccessible questions, including: (i) the origin of very tightly bound shells, (ii) the observability lifetime of shells and other tidal features following merger events, and (iii) the predicted frequency of visible shells in the galaxy population. A deeper exploration into these topics using the Illustris simulation and images presented here is deferred to a separate, targeted study.
4 Example Application: Derived Stellar Masses
An examination of the distribution of galaxies of various morphological types as a function of mass, the morphology-environment relation, and the role of mergers in driving morphological evolution will be considered elsewhere (Snyder et al., in prep.). Here we investigate a related problem which can be attacked with the image pipeline catalog: the accuracy of stellar mass estimates based on broadband photometric fitting. The simple question we ask is: How well do broadband fitting routines recover the stellar mass from the simulations? Similar topics have been addressed using semi-analytic models (e.g., Mitchell et al., 2013) and idealized hydrodynamical simulations (Wuyts et al., 2009a; Michałowski et al., 2014; Hayward & Smith, 2014).
For observational data, stellar masses are often determined by applying one of a few commonly adopted functional forms for the star formation history and finding the best-fitting stellar mass, galaxy age, star formation timescale, dust attenuation parameters, and metallicity to the observed SED. One common assumption is an exponentially declining star formation history which is characterized by two free parameters: the galaxy formation time and the star formation timescale. The simplest method for constraining these parameters is to perform a chi-squared minimization of the observed SED or broadband photometry relative to a grid of mock SEDs which uniformly cover the parameter space. This is the general approach that is taken in, e.g., FAST (Kriek et al., 2009).
To perform this exercise we assign light to a galaxy based on a stellar population synthesis template (e.g., the SB99 model), and then use the same (or a similar) stellar population synthesis model to find a best fit to the SED. One might therefore worry that this is a circular argument, which will not shed any light on the validity of “true” galaxy mass estimates. However, the star formation history and heavy element composition for our simulated galaxies is known and properly reflected in the integrated galaxy SED, while the observational models must assume a functional form for the star formation history and metallicity. The discrepancy leads, in fact, to the two recommendations for improved stellar mass determination, described below.
To make the comparison as clear and idealized as possible, we neglect dust in the simulated galaxy SEDs, and force FAST to neglect dust obscuration as well (i.e. set ). Neglecting dust removes a large age/dust degeneracy which is present in real extragalactic data and mass measurements. Furthermore, we ignore dust re-radiation at long wavelengths – which would add uncertainty into both our models and model SED fitting. Finally, we have not included any photometric and/or distance related errors in the galaxy’s positioning. We assume that we know a galaxy’s distance exactly (without error), and assign apparent magnitudes in each band appropriately. Owing to these simplifying assumptions, the comparisons that we perform below are highly idealized. We expect that introducing any additional SED complications (dust, distance errors, etc.) would bring the derived stellar masses into worse agreement than what we present below because there would be additional characteristics of the intrinsic galaxy SEDs that are not fully accounted for by FAST. Because of these simplifying assumptions, we will show below that we can obtain relatively good mass estimates for our mock galaxies based on their mock SEDs using only 5 bands. However, we find that restricting the allowable age range and metallicity values that FAST is permitted to sample can improve the derived stellar mass errors.
4.1 Comparison: Fast vs. Subfind stellar masses
We employ the code FAST (Kriek et al., 2009) to fit stellar population synthesis models to the artificially generated image pipeline broadband photometry. We assume an exponentially declining star formation history () to determine stellar masses based on the stellar-light-only galaxy photometry from the image pipeline. For the initial comparison, we allow FAST to select from the full range of input SED possibilities. Specifically, the galaxy age is taken in the range years with logarithmic steps of dex, the star formation timescale is chosen in the range years with logarithmic steps of dex, and the metallicity is selected from .
The mock galaxy photometry uses the SB99 templates while FAST uses the BC03 templates.
FAST input files are generated by calculating rest-frame broadband fluxes for the five SDSS bands (u, g, r, i, z) converted to units of flux density assuming the galaxies are at redshift using the same cosmology employed in the Illustris simulation.
The left panel of Figure 9 shows a two dimensional histogram indicating the relationship between FAST photometrically derived stellar mass and the SUBFIND tabulated stellar mass for all 41,517 galaxies with more than 500 star particles at redshift . The black line indicates a perfect 1:1 relationship (i.e. the FAST masses exactly equal the SUBFIND masses). There is a tight relationship between these two quantities that nearly follows the ideal 1:1 scaling, with a slight offset and scatter. Out of the 41,517 galaxies for which we perform this fitting, 1755 (4%) disagree by more than a factor of two and 10 (0.02%) disagree by more than a factor of five. We calculated the average (median) star formation rate histories for galaxies in a mass bin () for several Log(M/ M) values. The resulting plot is shown in Figure 10 and can be used to partially determine what drives the FAST mass estimate. The median star formation histories evolve with this “mass error” parameter. Galaxies with the lowest mass error parameter, which correspond to galaxies whose masses have been most severely underestimated, have the oldest stellar populations and earliest peak in their SFR histories. Galaxies with old stellar populations will have the largest mass-to-light ratios, which drives their systematically underestimated stellar masses if FAST does not choose a suitably large galaxy age. Moving toward larger (and positive) mass error parameters we find that the average SFR history becomes biased toward later times and the present day star formation rate monotonically increases. For the most extreme positive mass error parameters, we find that they are heavily peaked toward very recent star formation. The most extreme objects have either very early or late formation histories which are not captured in the FAST best fit parameters.
The left panel of Figure 11 shows the stellar half mass assembly time. This can be compared against the center panel of Figure 11 which shows the stellar half mass assembly time based on the assumed FAST star formation history. In general the two distributions look very different. The FAST age values span a larger range than found in SUBFIND. We plot the ratio of the SUBFIND age to the FAST age in the right panel of Figure 11 as a function of galaxy mass. There is significant scatter in this relation that is fairly mass dependent, with massive galaxies being assigned more accurate age estimates than low mass systems. This mis-estimation of galaxy formation times can cause problems for the inferred mass-to-light ratios.
The right panel of Figure 9 shows the ratio of the FAST masses to the SUBFIND masses as a function of the SUBFIND mass, which helps to identify some residual trends in the derived FAST masses. There is a mass dependence in the average FAST mass offset and scatter which can be ascertained from the distribution of galaxies outlined by the contour lines. We find that a principle origin of the mass dependent offset and scatter is the discrete and sparsely sampled metallicity values in the FAST BC03 stellar population synthesis templates. The standard BC03 templates used in FAST draw from discrete metallicities values of without interpolation, and the result is a discrete jump in the derived stellar masses depending on which template metallicity is applied. The left panel of Figure 12 shows the average (mean) template metallicity that FAST applied to galaxies in each pixel. Since FAST is forced to pick one of these sparsely populated metallicity values, the derived galaxy mass offsets correlate directly with different FAST metallicity values. The right panel of Figure 12 shows the mass weighted average stellar metallicity for these same galaxies. We can gather from this plot that, even in this highly idealized experiment, FAST is not able to accurately determine a galaxy’s metallicity when using only 5 broadband filters. Moreover, allowing for FAST to assign metallicity values to each galaxy independently results in galaxies of similar intrinsic mass being assigned factor of different stellar mass estimates. We have further checked – and we discuss in section 4.3 – that if we restrict FAST to using only a single metallicity () template, all derived galaxy masses fall along a single continuous relation, albeit with an offset from the desired relation.
4.2 Alternative Assumed Star Formation Histories
Here we re-compute FAST masses by employing delayed or single-burst star formation histories. Both the delayed and single burst star formation histories are described by two parameters: a galaxy formation time and star formation timescale, . The galaxy formation time gives the lookback time when the galaxy started forming stars (i.e. for all times before this, the star formation rate is assumed to be zero). In the case of the delayed star formation history (), the star formation timescale modulates the decay timescale for the late time star formation rate decline. In the case of the single burst star formation history, the star formation rate is assumed to be constant after the galaxy formation time for one star formation timescale, after which it drops back to zero ( for ).
Figure 13 shows the ratio of the FAST mass to the SUBFIND mass as a function of the SUBFIND mass for the delayed star formation history (left) and single burst star formation history (right). The distribution of galaxies in the delayed star formation history space looks similar to what was found for the previously assumed exponential star formation history. By examining a number of “best fit” star formation histories from the exponential and delayed models, we find that FAST prefers very similar (and young) star formation history in many cases, explaining why the distribution of derived galaxy masses is so similar.
The distribution of derived galaxy masses for the single burst star formation rate history is more complex. Figure 14 show the average metallicity (left) and age (right). There are three metallicities assigned to galaxies, causing mass estimate continuity issues as discussed above. The metallicity alone cannot explain the large spread in the galaxy population at . The right panel of Figure 14 shows that there is a large age gradient, where the location along this gradient determines a galaxy’s FAST mass relative to the intrinsic SUBFIND mass. Galaxies that are assigned old ages (and long star formation timescales) yield larger derived mass estimates (owing to their larger mass-to-light ratios). The origin of the large age spread for the low-mass galaxy population likely arises from the diversity of recent star formation histories of these low mass systems, which can substantially impact the best fit SED that FAST selects. It has been found, for instance, that when bursty systems are fit with single-component SFHs, their stellar masses can be significantly underestimated (e.g., Michałowski et al., 2012, 2014). Regardless of the underlying cause, we conclude that when FAST assigns a very young age to a galaxy it is likely that the galaxy’s mass will be underestimated.
4.3 Holding the Metallicity Fixed
In the previous subsections, we have found that allowing FAST to pick a galaxy’s metallicity during the SED fitting process can lead to mass-dependent errors in the derived galaxy masses. Here we re-compute FAST masses while holding the assumed metallicity fixed at (i.e. solar metallicity). This is a commonly adopted assumption when deriving stellar masses for real observational data (e.g., Wuyts et al., 2012). The ratio of the FAST mass to the SUBFIND mass as a function of the SUBFIND mass for a fixed metallicity with an assumed exponential star formation rate history is shown in the left panel of Figure 15. There is an offset between the derived mass and the intrinsic mass that evolves with intrinsic mass. This trend is stronger than what was found when the metallicity was allowed to vary (Figure 9). At the low mass end, the systematic error in the estimated mass is dex, while the scatter about this value is roughly dex.
Employing a fixed metallicity provides a more continuous and less broad distribution of mass estimates with respect to the intrinsic mass. In particular, holding the metallicity fixed in the FAST calculation reduces some of the variation in the mass ratios previously found in Section 4.1. However, the estimated mass error shows a clear mass dependent offset from the intrinsic mass and there is still substantial scatter in the derived masses. The right hand panel of Figure 15 shows the average age that FAST has assigned to each galaxy. For all exercises up until this point, we have allowed the FAST age selection to be selected from years. Examining FAST’s selection for the derived galaxy age value, there is a very clear residual trend where FAST underestimates the mass of a galaxy when it assigns a young age (i.e. years) to that galaxy. This is consistent with recent findings in the literature that young galaxy age assignment values can lead to underestimates in the stellar mass (e.g., Wuyts et al., 2012; Maraston et al., 2013; Price et al., 2014).
4.4 Restricting the Galaxy Age
In an attempt to improve the accuracy of the estimated masses, we restrict the galaxy age range from which FAST can select to years. The resulting ratio of the FAST mass to the SUBFIND mass as a function of the SUBFIND mass (still for a fixed metallicity with an assumed exponential star formation rate history) is shown in Figure 16. This case has the lowest systematic offset, mass dependence, and scatter in the estimated mass error of all of the cases presented in this paper. The origin of this reduced bias is fairly clear: a relatively low mass of recently formed stars (say, of a galaxy’s total stellar mass) can contribute substantially to the integrated SED – especially at short wavelengths – forcing FAST to assume a relatively young age for galaxy as a whole when only a single exponential star formation rate history is assumed. By restricting the minimum galaxy age to be years, we effectively prevent low galaxy ages (with correspondingly low mass-to-light ratios) from being assigned. We would therefore encourage caution when allowing a single galaxy age parameter to vary in an unconstrained manner, when determining the stellar masses of real systems. For redshift we find that a lower age bound of years works well for our simulated galaxies. However, this value will need to be lowered for high redshift galaxies. It is also possible that this issue could be solved by adopting multi-component star formation histories (e.g., an exponential with superimposed bursts Michałowski et al., 2012, 2014).
4.5 Spatially Resolved Derived Galaxy Properties
In this subsection we apply the FAST broadband fitting procedure to spatially resolved galaxy images. By performing the fitting procedure in each pixel we estimate the resolved surface density distribution which contains information about the dynamical state of the system as predicted by the simulation. We first perform this analysis on the idealized image of a single redshift galaxy without applying any restrictions on the FAST best fit parameters. The central panel of Figure 17 shows the stellar mass surface density as derived from applying FAST to each pixel of the image. For reference, the projected stellar surface density taken directly from the simulation for this galaxy is shown in the left panel of Figure 17. This galaxy (ID=283832) possesses a radially declining surface density profile as seen in the left panel of Figure 17. The stellar surface density profile is recovered in the FAST analysis, with some detailed differences. We find “holes” in the FAST stellar mass surface density profiles which are present in pixels which have significant contributions from young star particles. These pixels are assigned young ages owing to the influence of the luminous young stellar population despite the fact that the majority of the mass in these areas is characterized by an older average age. The low mass-to-light ratios associated with the young stellar ages result in an underestimation of the stellar mass surface density in these regions.
We next restrict the FAST best fit parameters based on the minimum age and fixed metallicity criteria that we prescribe in the previous section. The right panel of Figure 17 shows the stellar mass surface density map derived from FAST when we restrict the minimum stellar age to years and fix the stellar metallicity to solar. The derived stellar mass surface density map is qualitatively similar to that shown in the central panel of Figure 17, with the surface density holes being somewhat less pronounced. Because there are no pixels in this particular galaxy image which are mass dominated by young stellar populations we find that applying the minimum stellar age objectively improves the derived stellar mass surface density map. Applying a fixed metallicity impacts the derived stellar mass surface density maps at the 0.1 dex level, which is subdominant compared to the adopted age.
The Illustris Simulation Observatory presented in this paper is a large mock galaxy image catalog constructed directly from a hydrodynamical simulation. In the past, large mock galaxy catalogs have been produced by applying semi-analytic models of galaxy formation to the output of dark matter only simulations (e.g., Kitzbichler & White, 2007; Henriques et al., 2012; Overzier et al., 2013). The approach used in the MRObs is to assign light to a bulge and disk component of a galaxy separately based on assumptions about their radial surface brightness profiles by using the bulge-to-total ratio, bulge half light radius, and disk scale height (Overzier et al., 2013). These galactic light components are then included in images by assuming a sky position angle and inclination. After projecting all galaxies from a conic observational volume onto an image, point spread functions taken from the appropriate instruments are applied to yield mock light cones that compare against real observations favorably (Overzier et al., 2013). While these models have the advantage that they can be run on very large simulated volumes – allowing them to produce mock light-cones and mock survey data – they can only indirectly take into account the role that galactic dynamics (e.g., galaxy bar formation, major/minor mergers, disk instabilities, etc.) play in morphological evolution.
In contrast, hydrodynamic simulations model gas processes and internal galactic dynamics self-consistently down to the gravitational softening limit. The resulting sample of simulated galaxies can therefore be directly converted into synthetic images which encode information about the stellar age distribution, stellar metallicity distribution, and stellar spatial distribution. No further assumptions need to be made in order to generate a galaxy image. The resulting light distribution from galaxies may therefore end up as disk-dominated, bulge-dominated, irregular, or highly asymmetric as dictated by the formation history and internal galactic dynamics of each system. In this paper we have limited our focused to the construction of idealized images. In a companion paper, Snyder et al., (in prep) outlines our steps for adding image realism to the idealized galaxy images presented in this paper. Specifically, after creating the idealized images one can in post-processing rebin the images to an appropriate pixel scale, convolve with a point spread function, add background images or noise at a level consistent with a given instrument and exposure time. Analysis of the idealized images presented in this paper after being put through such post-processing steps will be soon presented in a companion paper (Snyder et al., in prep).
A series of important limitations need to be considered in order to properly exploit the Illustris output. An incomplete list of these is as follows:
The redshift gravitational softening length in the simulation is kpc, making small systems and small-scale features within large systems poorly resolved. While we do observe barred galaxies formed during merger/flyby events, force softening can suppress disk instabilities and bar formation in non-interacting galaxies (Kaufmann et al., 2007). This can lead to a low fraction of barred systems, and potentially modulate the efficiency of stellar bulge formation (to the extent that bars help drive bulge formation).
The mass resolution of our models may impact results, particularly for low mass galaxies. Given our gas/star particle mass resolution of , young stellar populations will be stochastically sampled with a small number of particles. This has been shown in the past to potentially bias non-parametric fitting parameters, such as the Gini coefficient (Lotz et al., 2008; Jonsson et al., 2010).
Our current models make use of a pressurized equation of state to handle the multi-phase ISM in a sub-grid fashion (Springel & Hernquist, 2003a). Features that are expected to exist in a true multiphase ISM, such as spiral arm structure driven by pressure waves, will not be properly captured. Some simulated galaxies – such as the disk galaxies presented in Figure 5 – do show signs of spiral features. However, this spiral arm structure is driven by swing amplification (Goldreich & Lynden-Bell, 1965; Julian & Toomre, 1966) of perturbations to the halo potential caused by the representation of dark matter by discrete simulation particles (D’Onghia et al., 2013).
Our assumed physics (and specifically our feedback implementation) can influence galaxy morphology. We find, for instance, a number of galaxies with enhanced rings of star formation in our simulation, which are likely caused by the interaction of star formation driven winds with the dense galactic ISM. The presence of these rings and/or the specific feedback formulation used in our model may influence visual and automated galaxy morphology classification or galaxy size growth.
Our mock images and SEDs do not contain consider detailed radiation transfer. The kilo parsec scale spatial resolution of the Illustris simulation is slightly too coarse to carry out self-consistent dust modeling. We have applied simple empirical corrections based on the CF00 model to the galaxy SEDs. Analysis of any systems spatial light distributions and/or SED properties where dust is expected to play a major role should therefore use caution.
Keeping in mind these caveats, we still expect the Illustris Simulation Observatory to be a useful and reliable resource in many situations, including: (i) addressing issues that cannot be accessed with real data (e.g. the viewing angle dependence of galaxy morphology and time evolution of galactic features), (ii) studying situations where the internal galaxy dynamics governing the spatial distribution of stars is of the essence (e.g. in the formation and evolution of stellar shells or the viewing angle dependence of bar classification, as explored briefly in this paper), and (iii) where the full hydro simulation approach of the Illustris simulation is likely robust (e.g. for the formation of bars via mergers/interactions, the observability timescales for galaxy pairs and tidal features, the puffing up of massive galaxies via dry mergers, the formation of disk galaxies via in-situ star formation, etc.). Moreover, using an even-handed determination of galaxy morphology based on synthetic images will help us further identify the environments and mass ranges where our simulated galaxy populations are successfully or poorly replicating observations (Snyder et al., in prep).
We expect that the Illustris Simulation Observatory will be an evolving resource, not only for the current Illustris simulation, but also for next generation simulations with higher resolution, larger volume, and more sophisticated physics implementations.
The main goal of this paper is the presentation of a large synthetic galaxy image catalog based on a high-resolution full-volume cosmological hydrodynamical simulation. We described our method for building this image catalog, in particular, the assignment of stellar light to each star particle within a galaxy. We have made 7,000 galaxy images in 36 broadbands plus integrated SEDs for 40,000 systems. We have demonstrated that these broadband images can be used to generate high quality galaxy images. Our intention when constructing this mock galaxy image database is that it may serve as a useful tool for observers in the interpretation of extragalactic data. In principle, the FITS file format of our mock galaxy images will make it easy to process these images with existing image analysis tools, similar to what was done in Overzier et al. (2013). This should facilitate galaxy selection from our catalog using selection criteria matched to observational targets. It also implies that spatially resolved galaxy properties (e.g., half-light-radii, bulge-to-disk ratios, etc.) can be derived using un-modified photometric image decomposition or fitting techniques.
As a first application, we showed that the SED fitting code FAST is able to derive the “correct” galaxy mass within a factor of 2 for 95% of galaxies in our idealized sample. We also uncovered some residual trends in the stellar mass estimates and argued that sparsely populated stellar population synthesis template metallicity values and low galaxy age assignments are responsible for the majority of the error in the photometrically derived galaxy mass values. Although our idealized comparison cannot address systematic uncertainties such as the dependence of the mass-to-light ratio for a galaxy population on the shape of the IMF, we have used our mock image pipeline to show that, overall, FAST is able to derive reasonably accurate galaxy stellar masses for complex star formation histories using a relatively simple assumption for those histories. We concluded that, for redshift , employing a lower age bound of years with a fixed solar metallicity produces the smallest errors in the stellar mass determination.
Several additional applications of the mock image pipeline are currently under way. We are undertaking automated non-parametric measurements of galaxy morphology for direct comparison with survey data (Snyder et al., in prep). Among other goals, this will allow us to quantitatively assess the extent to which the morphological mix of the simulated galaxies match reality. Visual classification of the galaxy population for galaxies with stellar particles is also planned. This selection criterion results in a sample of 1000 galaxies with stellar masses , which we have verified contains a mix of both disks and elliptical systems. Visual classification is an alternative to the automated non-parametric classifications, and enables us to directly address whether, for example, the disk-to-elliptical fraction for simulated galaxies is comparable to the values measured from SDSS in GalaxyZoo (Lintott et al., 2008).
By making a fair comparison between cosmological simulations and large galaxy surveys, we can begin to connect theoretical ideas on the origins of galaxy structure with actual observations in ever more accurate ways. The role that simulations play in shaping our understanding of galaxy formation will continue to increase as the models become increasingly more detailed and physically accurate. Comparing the detailed structure and properties of individual galaxies – within the context of their larger populations – is an important point of overlap between theory and observation, and mock catalogs based on realistic, synthetic observations of simulated systems are a fundamental requisite for maximizing the insight gained from that overlap.
We thank Dave Sanders and Arjun Dey for helpful suggestions on this work. We thank Patrik Jonsson for his significant efforts to write, develop, and support the SUNRISE code. GS acknowledges support from the HST grants program, number HST- AR-12856.01-A. Support for program #12856 (PI Lotz) was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. CCH is grateful to the Klaus Tschira Foundation for financial support and acknowledges the hospitality of the Aspen Center for Physics, which is supported by the National Science Foundation Grant No. PHY-1066293. VS acknowledges support by the DFG Research Centre SFB-881 The Milky Way System through project A1, and by the European Research Council under ERC-StG EXAGAL-308037. LH acknowledges support from NASA grant NNX12AC67G and NSF grant AST-1312095.
- E-mail: firstname.lastname@example.org
- Moore Prize Postdoctoral Scholar in Theoretical Astrophysics
- footnotetext: The synthetic images and spectra presented in this paper are available at http://www.illustris-project.org/galaxy_obs/
- Sunrise is freely available and documented at http://code.google.com/p/sunrise/.
- The SB99 and BC03 model templates agree well enough in the optical bands that this does not introduce a major source of error into our analysis. However, since the old stellar populations have different UV emission in the two models, very poor fits can be identified when including, e.g., the GALEX bands.
- We note that it is slightly inconsistent to use rest-frame broadband fluxes for input into FAST (which assumes redshifted fluxes), but we have checked that the neglected k-correction error is minimal for such a low redshift.
- Abadi, M. G., Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, ApJ, 597, 21
- Agertz, O., Teyssier, R., & Moore, B. 2011, MNRAS, 410, 1391
- Aguirre, A., Hernquist, L., Schaye, J., et al. 2001, ApJ, 561, 521
- Banerji, M., Glazebrook, K., Blake, C., et al. 2013, MNRAS, 431, 2209
- Bernyk, M., Croton, D. J., Tonini, C., et al. 2014, ArXiv e-prints
- Binette, L., Dopita, M. A., & Tuohy, I. R. 1985, ApJ, 297, 476
- Boylan-Kolchin, M., Springel, V., White, S. D. M., Jenkins, A., & Lemson, G. 2009, MNRAS, 398, 1150
- Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000
- Bruzual A., G. 1983, ApJ, 273, 105
- Bruzual A., G., & Charlot, S. 1993, ApJ, 405, 538
- Buzzoni, A. 1989, ApJS, 71, 817
- Canalizo, G., Bennert, N., Jungwiert, B., et al. 2007, ApJ, 669, 801
- Carter, D. 1979, MNRAS, 186, 897
- Cen, R., & Fang, T. 2006, ApJ, 650, 573
- Cen, R., Miralda-Escudé, J., Ostriker, J. P., & Rauch, M. 1994, ApJ, 437, L9
- Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718
- Conroy, C., & Gunn, J. E. 2010, ApJ, 712, 833
- Conroy, C., & Wechsler, R. H. 2009, ApJ, 696, 620
- Conroy, C., White, M., & Gunn, J. E. 2010, ApJ, 708, 58
- Cooper, A. P., Martínez-Delgado, D., Helly, J., et al. 2011, ApJ, 743, L21
- Cox, T. J., Jonsson, P., Primack, J. R., & Somerville, R. S. 2006, MNRAS, 373, 1013
- Crain, R. A., Theuns, T., Dalla Vecchia, C., et al. 2009, MNRAS, 399, 1773
- Croft, R. A. C., Di Matteo, T., Springel, V., & Hernquist, L. 2009, MNRAS, 400, 43
- D’Onghia, E., Vogelsberger, M., & Hernquist, L. 2013, ApJ, 766, 34
- Dopita, M. A., Groves, B. A., Fischera, J., et al. 2005, ApJ, 619, 755
- Dopita, M. A., Fischera, J., Sutherland, R. S., et al. 2006a, ApJ, 647, 244
- —. 2006b, ApJS, 167, 177
- Dupraz, C., & Combes, F. 1986, A&A, 166, 53
- Faucher-Giguère, C.-A., Lidz, A., Zaldarriaga, M., & Hernquist, L. 2009, ApJ, 703, 1416
- Faucher-Giguère, C.-A., Prochaska, J. X., Lidz, A., Hernquist, L., & Zaldarriaga, M. 2008, ApJ, 681, 831
- Fosalba, P., Gaztañaga, E., Castander, F. J., & Manera, M. 2008, MNRAS, 391, 435
- Gallazzi, A., & Bell, E. F. 2009, ApJS, 185, 253
- Genel, S., Vogelsberger, M., Springel, V., et al. 2014, ArXiv e-prints 1405.3749
- Goldreich, P., & Lynden-Bell, D. 1965, MNRAS, 130, 125
- Governato, F., Mayer, L., Wadsley, J., et al. 2004, ApJ, 607, 688
- Groves, B., Dopita, M. A., Sutherland, R. S., et al. 2008, ApJS, 176, 438
- Groves, B. A., Dopita, M. A., & Sutherland, R. S. 2004a, ApJS, 153, 9
- —. 2004b, ApJS, 153, 75
- Guedes, J., Callegari, S., Madau, P., & Mayer, L. 2011, ApJ, 742, 76
- Hayward, C. C., Behroozi, P. S., Somerville, R. S., et al. 2013a, MNRAS
- Hayward, C. C., Jonsson, P., Kereš, D., et al. 2012, MNRAS, 424, 951
- Hayward, C. C., Kereš, D., Jonsson, P., et al. 2011, ApJ, 743, 159
- Hayward, C. C., Narayanan, D., Kereš, D., et al. 2013b, MNRAS, 428, 2529
- Hayward, C. C., & Smith, D. J. B. 2014, arXiv:1409.6332
- Hayward, C. C., Torrey, P., Springel, V., Hernquist, L., & Vogelsberger, M. 2014a, MNRAS, 442, 1992
- Hayward, C. C., Lanz, L., Ashby, M. L. N., et al. 2014b, MNRAS, 445, 1598
- Henriques, B., Maraston, C., Monaco, P., et al. 2011, MNRAS, 415, 3571
- Henriques, B. M. B., White, S. D. M., Lemson, G., et al. 2012, MNRAS, 421, 2904
- Hernquist, L., Katz, N., Weinberg, D. H., & Miralda-Escudé, J. 1996, ApJ, 457, L51
- Hernquist, L., & Quinn, P. J. 1988, ApJ, 331, 682
- Hernquist, L., & Spergel, D. N. 1992, ApJ, 399, L117
- Hopkins, P. F., Hernquist, L., Martini, P., et al. 2005, ApJ, 625, L71
- Jonsson, P. 2006, MNRAS, 372, 2
- Jonsson, P., Groves, B. A., & Cox, T. J. 2010, MNRAS, 403, 17
- Julian, W. H., & Toomre, A. 1966, ApJ, 146, 810
- Katz, N., Hernquist, L., & Weinberg, D. H. 1992, ApJ, 399, L109
- Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJS, 105, 19
- Kaufmann, T., Mayer, L., Wadsley, J., Stadel, J., & Moore, B. 2007, MNRAS, 375, 53
- Kereš, D., Katz, N., Weinberg, D. H., & Davé, R. 2005, MNRAS, 363, 2
- Kereš, D., Vogelsberger, M., Sijacki, D., Springel, V., & Hernquist, L. 2012, MNRAS, 425, 2027
- Kitzbichler, M. G., & White, S. D. M. 2007, MNRAS, 376, 2
- Klypin, A. A., Trujillo-Gomez, S., & Primack, J. 2011, ApJ, 740, 102
- Kojima, M., & Noguchi, M. 1997, ApJ, 481, 132
- Kriek, M., van Dokkum, P. G., Labbé, I., et al. 2009, ApJ, 700, 221
- Lanz, L., Hayward, C. C., Zezas, A., et al. 2014, ApJ, 785, 39
- Leitherer, C., Ortiz Otálvaro, P. A., Bresolin, F., et al. 2010, ApJS, 189, 309
- Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3
- Lemson, G., & Springel, V. 2006, in Astronomical Society of the Pacific Conference Series, Vol. 351, Astronomical Data Analysis Software and Systems XV, ed. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique, 212
- Lintott, C. J., Schawinski, K., Slosar, A., et al. 2008, MNRAS, 389, 1179
- Lotz, J. M., Jonsson, P., Cox, T. J., & Primack, J. R. 2008, MNRAS, 391, 1137
- Lupton, R., Blanton, M. R., Fekete, G., et al. 2004, PASP, 116, 133
- Maraston, C. 1998, MNRAS, 300, 872
- —. 2005, MNRAS, 362, 799
- Maraston, C., Pforr, J., Henriques, B. M., et al. 2013, MNRAS, 435, 2764
- Marinacci, F., Pakmor, R., & Springel, V. 2014, MNRAS, 437, 1750
- McQuinn, M., Lidz, A., Zaldarriaga, M., et al. 2009, ApJ, 694, 842
- Michałowski, M. J., Dunlop, J. S., Cirasuolo, M., et al. 2012, A&A, 541, A85
- Michałowski, M. J., Hayward, C. C., Dunlop, J. S., et al. 2014, ArXiv e-prints 1405.2335
- Mitchell, P. D., Lacey, C. G., Baugh, C. M., & Cole, S. 2013, MNRAS, 435, 87
- Murali, C., Katz, N., Hernquist, L., Weinberg, D. H., & Davé, R. 2002, ApJ, 571, 1
- Nelson, D., Vogelsberger, M., Genel, S., et al. 2013, MNRAS, 429, 3353
- Ocvirk, P., Pichon, C., & Teyssier, R. 2008, MNRAS, 390, 1326
- Ono, Y., Ouchi, M., Shimasaku, K., et al. 2010, ApJ, 724, 1524
- Oppenheimer, B. D., Davé, R., Kereš, D., et al. 2010, MNRAS, 406, 2325
- Overzier, R., Lemson, G., Angulo, R. E., et al. 2013, MNRAS, 428, 778
- Papovich, C., Dickinson, M., & Ferguson, H. C. 2001, ApJ, 559, 620
- Pedrosa, S. E., Tissera, P. B., & De Rossi, M. E. 2014, A&A, 567, A47
- Price, S. H., Kriek, M., Brammer, G. B., et al. 2014, ApJ, 788, 86
- Quinn, P. J. 1984, ApJ, 279, 596
- Robertson, B., Li, Y., Cox, T. J., Hernquist, L., & Hopkins, P. F. 2007, ApJ, 667, 60
- Sales, L. V., Navarro, J. F., Theuns, T., et al. 2012, MNRAS, 423, 1544
- Scannapieco, C., Gadotti, D. A., Jonsson, P., & White, S. D. M. 2010, MNRAS, 407, L41
- Schaerer, D., & de Barros, S. 2009, A&A, 502, 423
- Schaye, J., Dalla Vecchia, C., Booth, C. M., et al. 2010, MNRAS, 402, 1536
- Shen, S., Madau, P., Guedes, J., et al. 2013, ApJ, 765, 89
- Sijacki, D., & Springel, V. 2006, MNRAS, 366, 397
- Sijacki, D., Springel, V., Di Matteo, T., & Hernquist, L. 2007, MNRAS, 380, 877
- Snyder, G. F., Cox, T. J., Hayward, C. C., Hernquist, L., & Jonsson, P. 2011, ApJ, 741, 77
- Snyder, G. F., Hayward, C. C., Sajina, A., et al. 2013, ApJ, 768, 168
- Snyder, G. F., Lotz, J., Moody, C., et al. 2014, ArXiv e-prints
- Somerville, R. S., Gilmore, R. C., Primack, J. R., & Domínguez, A. 2012, MNRAS, 423, 1992
- Springel, V. 2005, MNRAS, 364, 1105
- —. 2010, MNRAS, 401, 791
- Springel, V., Di Matteo, T., & Hernquist, L. 2005a, MNRAS, 361, 776
- Springel, V., & Hernquist, L. 2003a, MNRAS, 339, 289
- —. 2003b, MNRAS, 339, 312
- Springel, V., White, M., & Hernquist, L. 2001, ApJ, 549, 681
- Springel, V., White, S. D. M., Jenkins, A., et al. 2005b, Nature, 435, 629
- Stark, D. P., Schenker, M. A., Ellis, R., et al. 2013, ApJ, 763, 129
- Stinson, G. S., Brook, C., Macciò, A. V., et al. 2013, MNRAS, 428, 129
- Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253
- Teyssier, R., Pires, S., Prunet, S., et al. 2009, A&A, 497, 335
- Theuns, T., Leonard, A., Efstathiou, G., Pearce, F. R., & Thomas, P. A. 1998, MNRAS, 301, 478
- Thomas, D., Maraston, C., & Bender, R. 2003, MNRAS, 339, 897
- Tinsley, B. M. 1972, ApJ, 178, 319
- Torrey, P., Vogelsberger, M., Genel, S., et al. 2014, MNRAS, 438, 1985
- Torrey, P., Vogelsberger, M., Sijacki, D., Springel, V., & Hernquist, L. 2012, MNRAS, 427, 2224
- van de Voort, F., & Schaye, J. 2012, MNRAS, 423, 2991
- van de Voort, F., Schaye, J., Booth, C. M., Haas, M. R., & Dalla Vecchia, C. 2011, MNRAS, 414, 2458
- Vázquez, G. A., & Leitherer, C. 2005, ApJ, 621, 695
- Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, MNRAS, 436, 3031
- Vogelsberger, M., Sijacki, D., Kereš, D., Springel, V., & Hernquist, L. 2012, MNRAS, 425, 3024
- Vogelsberger, M., Genel, S., Springel, V., et al. 2014a, ArXiv e-prints 1405.2921
- —. 2014b, Nature, 509, 177
- Weinberg, D. H., Hernquist, L., & Katz, N. 1997, ApJ, 477, 8
- Wiersma, R. P. C., Schaye, J., & Smith, B. D. 2009a, MNRAS, 393, 99
- Wiersma, R. P. C., Schaye, J., Theuns, T., Dalla Vecchia, C., & Tornatore, L. 2009b, MNRAS, 399, 574
- Wild, V., Walcher, C. J., Johansson, P. H., et al. 2009, MNRAS, 395, 144
- Worthey, G. 1994, ApJS, 95, 107
- Wuyts, S., Cox, T. J., Hayward, C. C., et al. 2010, ApJ, 722, 1666
- Wuyts, S., Franx, M., Cox, T. J., et al. 2009a, ApJ, 696, 348
- Wuyts, S., Labbé, I., Franx, M., et al. 2007, ApJ, 655, 51
- Wuyts, S., Franx, M., Cox, T. J., et al. 2009b, ApJ, 700, 799
- Wuyts, S., Förster Schreiber, N. M., Genzel, R., et al. 2012, ApJ, 753, 114
- Zackrisson, E., Bergvall, N., & Leitet, E. 2008, ApJ, 676, L9
- Zhang, Y., Anninos, P., & Norman, M. L. 1995, ApJ, 453, L57