AnHST/ACS investigation of the spatial and chemical structure and sub-structure of NGC 891, a Milky Way analogue

We studied a sample of 58 edge-on spiral galaxies at redshifts $z \sim 1$ selected in the Hubble Ultra Deep Field. For all galaxies we analyzed the 2D brightness distributions in the $V_{606}$ and $i_{775}$ filters and measured the radial ($h_r$) and vertical ($h_z$) exponential scales of the brightness distribution. We obtained evidence that the relative thickness of the disks of distant galaxies, i.e., the ratio of the vertical scale height and radial scale length ($h_z/h_r$), on average, exceeds the relative thickness of the disks of nearby spiral galaxies. The vertical scale height $h_z$ of the stellar disks of galaxies shows no big changes at $z \leq 1$. The possibility of the evolution of the radial scale length $h_r$ for the brightness distribution with redshift is discussed.

The galaxy within which we reside, the Milky Way, offers perhaps the highest fidelity training ground for models of galaxy formation, providing insights into its history of formation and evolution on a star-by-star basis. Such insights are only truly useful for constraining disc galaxy formation models if we properly understand the wider context of the Milky Way among the plethora of galaxies in the Universe. This thesis aims to make progress toward answering the question of whether the Milky Way is a typical disc galaxy. Through the effort to answer this question, this thesis presents new measurements of the structure of our Galaxy and new insights into aspects of its history of assembly and evolution which have a strong influence on its α-element abundances. Providing a baseline constraint on future models for the formation of the Galaxy, I present a dissection of the Milky Way disc spatial structure as a function of [Fe/H], [α/Fe] and stellar ages (based on the surface abundances of Carbon and Nitrogen), as measured by the APOGEE survey. I measure the disc density profile, fitting the scale heights and lengths of mono-age, mono-[Fe/H] sub-populations in the high and low [α/Fe] disc. The fitted disc vertical scale height distribution is smooth when weighted by surface-mass density, suggesting that the high and low [α/Fe] populations are not vertically distinct at the solar radius, as would be expected if they were interchangeable with the geometric thin and thick disc components. I find that the surface density profile of low [α/Fe] mono-age, mono-[Fe/H] populations is best fit by a broken exponential, such that their density increases with R to a peak radius, and declines thereafter, whereas the high [α/Fe] populations are better described by a single, declining exponential within the range of Galactocentric radii observed by APOGEE. The trends of the density profile parameters as a function of age and metallicity provide insights into the structural evolution of the Galaxy, and provide strong constraints on future models of its formation and evolution. In particular, a main finding of this study is that the high and low [α/Fe] discs have a relatively distinct radial structure. To understand the origin of the structures found in the above study, and to generate novel and predictive models for the formation of the Galaxy, I perform an analysis of element abundances in Milky Way like galaxy discs in the EAGLE simulation. I concentrate on the abundance of α-elements in these galaxies, with a view to understanding the origin of the bimodality in [α/Fe] at fixed [Fe/H] which is apparent in both the Milky Way and in EAGLE Galactic analogues. I show that EAGLE reproduces broad expectations for the production of α-elements from simple chemical evolution models, namely that stars with enhanced [α/Fe] are born from gas in rapidly star forming environments at high density. These environments are conducive to forming α-enhanced stars because their gas consumption timescale is considerably shorter than the characteristic delay timescale for Type Ia supernovae. I further show that such environments are only achieved in EAGLE haloes when the rate of dark matter accretion is faster at early times than the majority of galaxies with similar stellar mass at z=0. The necessity for this atypical and rapid early assembly means that galaxies hosting discs with [α/Fe] bimodality are extremely rare, forming in ≈ 6% of galaxies at the Milky Way stellar mass range. Following the prediction of EAGLE that galaxies which host [α/Fe] bimodality like that of the Milky Way should have atypical histories of assembly, I then perform a study of Milky Way halo stars in common between APOGEE and Gaia DR2, with a view to placing constraints on the accretion history of the Galaxy. I present a detailed characterisation of the kinematics and abundances of the recently discovered Gaia-Enceladus association, which is proposed to be the debris of a singular and massive accretion event. By comparing the kinematics and abundances of Gaia-Enceladus to the debris of satellites accreted onto Milky Way analogues in the EAGLE simulations, I make a quantitative prediction of the stellar mass of the Milky Way satellite to be 10^8.5 < M* < 10^9 Msun and predict its earliest possible time of accretion to be z≈1.5. I also show that such mergers are uncommon in the simulations, in agreement with the prediction that the Milky Way assembly history should be atypical. The above results outline a new view of the Galaxy in which it is potentially not a good example of a `typical' disc galaxy, playing host to structural components which are linked to its chemistry in a complex manner, which may indicate that its history of assembly is not in common with other galaxies at the same stellar mass. In finding that the Galaxy is atypical in this way, this thesis has uncovered new aspects of the evolutionary history of the Milky Way, which pave the way for future work towards the goal of fully reconstructing the history of our Galaxy and using that understanding to formulate robust and general models for the formation of disc galaxies.

Detailed surface photometry of the edge-on spiral galaxy NGC 4565 was carried out by a computerized digital method. The luminosity distribution was analyzed in the standard manner and then decomposed into sub-components with the help of a spheroidal-shell model and an exponential-disk model. The results of the present analysis in the standard manner are in good agreement with those of previous investigations. The main results of the present model analysis are : (1) The spheroidal component of NGC 4565 consists of two spheroids, i.e., the central bulge and the extended halo. The indices are ;γ(bulge)~7.7 and ; γ(halo)=2.46 respectively, if the volume emissivity is expressed by the power law σκr-(γ+1) in an axisymmetric galaxy. The axial ratios are ~0.8 and 0.46 for the bulge and the halo respectively. (2) When the volume emissivity of the disk component is expressed by σ exp(–αr—β|z|), the radial scale length α−1 and the vertical scale height β−1 are 4.9 kpc and 0.05 kpc (distance 8.1 Mpc). (3) The face-on profile of NGC 4565, which was hypothetically generated by model calculation, is dominated by the disk component except for r≲1.′5. The B(0)c value in this model was found to be ~21.5 mag arcsec-2.

(abridged) We present an analysis of B-R and R-K_s color maps for 47 late-type, edge-on, unwarped, bulgeless disk galaxies spanning a wide range of mass. The color maps show that the thin disks of these galaxies are embedded within a low surface brightness red envelope that is substantially thicker than the thin disk (a/b~4:1 vs a/b>8:1), extends to at least 5 vertical disk scale heights above the galaxy midplane, and has a radial scale length that appears to be uncorrelated with that of the embedded thin disk. The color of the red envelope is similar from galaxy to galaxy and is consistent with a relatively old (>6Gyr) stellar population that is not particularly metal-poor. The color difference between the thin disk and the envelope varies systematically with rotation speed, indicating a younger thin disk relative to the red envelope in lower mass galaxies. The red stellar envelopes are similar to the MW thick disk, having common surface brightnesses, spatial distributions, mean ages, and metallicities. The ubiquity of the red stellar envelopes implies that the formation of the thick disk is a nearly universal feature of disk formation and need not be associated with bulge formation. Furthermore, our data suggest that the thick disk forms early, even in cases where the majority of star formation was recent. Finally, we find that our data and the observed properties of the MW thick disk argue in favor of a merger origin for the thick disk population. If so, then the age of the thick disk marks the end of the epoch of major merging, and the age difference between the thin and thick disks can become a strong constraint on cosmological constants and models of galaxy and/or structure formation.

We present a catalog of true edge-on disk galaxies automatically selected from the Seventh Data Release (DR7) of the Sloan Digital Sky Survey. A visual inspection of the $g$, $r$ and $i$ images of about 15000 galaxies allowed us to split the initial sample of edge-on galaxy candidates into 4768 (31.8% of the initial sample) genuine edge-on galaxies, 8350 (55.7%) non-edge-ons, and 1865 (12.5%) edge-on galaxies not suitable for simple automatic analysis because these objects show signs of interaction, warps, or nearby bright stars project on it. We added more candidate galaxies from RFGC, EFIGI, RC3, and Galaxy Zoo catalogs found in the SDSS footprints. Our final sample consists of 5747 genuine edge-on galaxies. We estimate the structural parameters of the stellar disks (the stellar disk thickness, radial scale length, and central surface brightness) in the galaxies by analyzing photometric profiles in each of the g, r, and i images. We also perform simplified 3-D modeling of the light distribution in the stellar disks of edge-on galaxies from our sample. Our large sample is intended to be used for studying scaling relations in the stellar disks and bulges and for estimating parameters of the thick disks in different types of galaxies via the image stacking. In this paper we present the sample selection procedure and general description of the sample.

\n To analyze the vertical structure of edge-on galaxies, we have used images\nof a large uniform sample of flat galaxies that have been taken during the\n2MASS all-sky survey. The photometric parameters, such as the radial scale\nlength, the vertical scale height, and the deprojected central surface\nbrightness of galactic disks have been obtained. We find a strong\ncorrelation between the central surface brightness and the ratio of the\nvertical scale height to the vertical scale length: the thinner the galaxy,\nthe lower the central surface brightness of its disk. The vertical scale \nheight does not increase systematically with the distance from the galaxy\ncenter in the frames of this sample.\n\n \n

There are stars in the halo of the Galaxy whose orbital properties are distinctly different from the majority of the halo population. One suggestion to account for these kinematically peculiar stars is that they have been gravitationally subsumed by the Milky Way through the breakup of nearby dwarf galaxies. One test of this hypothesis lies in determining the chemical composition of the deviant halo stars relative to that of normal (native) halo field stars. The possibly accreted stars are predicted to have different chemical enrichments due to a slower rate of star formation in nearby low-mass spheroidal or irregular galaxies. We have obtained echelle spectra of metal-poor halo dwarfs with unusual orbital properties: from the outer halo, from the halo, and on retrograde orbits. Our spectra are at high spectral resolution (35,000–48,000) and high signal-to-noise ratios (median value 140), primarily from the Keck I 10 m telescope with HIRES (53 stars), but also from the echelle spectrometer of the KPNO Mayall 4 m Telescope (three stars). The spectra cover a large spectral range: ~4000–6800 A. The strengths of approximately 9000 lines have been measured to determine the stellar parameters (Teff, log g, [Fe/H], and ξ) and the abundances of 10 elements, including Na, α-fusion products (Mg, Si, Ca, Ti), iron peak elements (Cr, Fe, Ni), and neutron capture elements (Y, Ba). A model atmosphere was computed for each star, which was used to determine the composition under the assumption of LTE. The comparison of the abundances with the kinematic properties revealed no special trends, except that the ratio of the mean of the α-fusion elements to Fe is weakly correlated with stellar apogalactica in the sense that the outermost stars have lower ratios of [α/Fe] than those with smaller apogalactica. The hallmark of an accreted halo star is presumed to be a deficiency (compared with normal stars) of the α-elements (e.g., O, Mg, Si, Ca, Ti) with respect to Fe, a consequence of sporadic bursts of star formation within the dwarf galaxies. However, we have found no new stars with such subsolar ratios of [α/Fe]. For BD +80°245 we derive [α/Fe] = -0.22 in agreement with earlier results. For all four of our α-element ratios, [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe], we find high values at low metallicities that decrease to near-solar values at high metallicities. The slope of the mean [α/Fe] with [Fe/H] is approximately -0.15. At a given metallicity, the outer halo stars have lower [α/Fe] than the inner halo stars. There is a spread in [α/Fe] beyond the individual measurement errors at both high and low [Fe/H]. It appears that the ratios of [Na/Fe] and [Mg/Fe] increase together. The ratio of [Ni/Fe] is persistently mildly deficient relative to solar. The ratio of [Ba/Fe] increases with [Fe/H]. The kinematically peculiar stars in our sample must have had their origins in localized star-forming regions far removed from the Galactic center. They do not carry the chemical signature of an accreted population. It seems clear that the chemical enrichment in our sample comes primarily from Type II supernovae with only a small component from Type Ia supernovae. Our stars would have formed at most ~1 Gyr after the beginning of a star formation episode in the remote halo. The general conclusion extracted from these data is that the formation of the nascent Milky Way was not dominated by the late accretion of dwarf galaxies like the ones that currently orbit the Galaxy. However, the assimilation of fragments early in the evolution of the Galaxy is a natural by-product of hierarchical models of structure formation and can explain many properties of the halo population.

The Milky Way’s stellar halo preserves a fossil record of smaller dwarf galaxies that merged with the Milky Way throughout its formation history. Currently, though, we lack reliable ways to identify which halo stars originated in which dwarf galaxies or even which stars were definitively accreted. Selecting stars with specific chemical signatures may provide a way forward. We investigate this theoretically and observationally for stars with r-process nucleosynthesis signatures. Theoretically, we combine high-resolution cosmological simulations with an empirically-motivated treatment of r-process enhancement. We find that around half of highly r-process-enhanced metal-poor halo stars may have originated in early ultra-faint dwarf galaxies that merged into the Milky Way during its formation. Observationally, we use Gaia DR2 to compare the kinematics of highly r-process-enhanced halo stars with those of normal halo stars. R-process-enhanced stars have higher galactocentric velocities than normal halo stars, suggesting an accretion origin. If r-process-enhanced stars largely originated in accreted ultra-faint dwarf galaxies, halo stars we observe today could play a key role in understanding the smallest building blocks of the Milky Way via this novel approach of chemical tagging

By means of N-body/SPH simulations we investigate the morphological properties of thick disks formed through minor mergers. We show that the vertical surface density profile of the post-merger thick disk follows a sech function and has an excess in the regions far from the disk mid-plane (z>2kpc). This stellar excess also follows a sech function with a larger scale height than the main thick disk component, and it is usually dominated by stars from the primary galaxy. Stars in the excess have a rotational velocity lower than that of stars in the thick disk, and they may thus be confused with stars in the inner galactic halo. The thick disk scale height increases with radius and the rate of its increase is smaller for more gas rich primary galaxies. On the contrary, the scale height of the stellar excess is independent of both radius and gas fraction. We also find that the post-merger thick disk has a radial scale length which is 10-50% larger than that of the thin disk. Two consecutive mergers have basically the same effect on heating the stellar disk as a single merger of the same total mass. To investigate how thick disks produced through secular processes may differ from those produced by minor mergers, we also simulated gravitationally unstable gas-rich disks. These disks do not produce either a stellar excess or a ratio of thick to thin disk scale lengths greater than one. Our simulation results are consistent with observations of the ratio of thick to thin disk scale lengths of the Milky Way and nearby galaxies, and with the Toomre diagram of the Milky Way. We conclude that minor mergers are a viable mechanism for the creation of galactic thick disks and investigating stars at several kpc above the mid-plane of the Milky Way and other galaxies may provide a quantitative method for studying the minor merger history of galaxies.

Aim.The vertical halo scale height is a crucial parameter to understand the transport of cosmic-ray electrons (CRE) and their energy loss mechanisms in spiral galaxies. Until now, the radio scale height could only be determined for a few edge-on galaxies because of missing sensitivity at high resolution.Methods.We developed a sophisticated method for the scale height determination of edge-on galaxies. With this we determined the scale heights and radial scale lengths for a sample of 13 galaxies from the CHANG-ES radio continuum survey in two frequency bands.Results.The sample average values for the radio scale heights of the halo are 1.1 ± 0.3 kpc inC-band and 1.4 ± 0.7 kpc inL-band. From the frequency dependence analysis of the halo scale heights we found that the wind velocities (estimated using the adiabatic loss time) are above the escape velocity. We found that the halo scale heights increase linearly with the radio diameters. In order to exclude the diameter dependence, we defined a normalized scale heighth˜which is quite similar for all sample galaxies at both frequency bands and does not depend on the star formation rate or the magnetic field strength. However,h˜shows a tight anticorrelation with the mass surface density.Conclusions.The sample galaxies with smaller scale lengths are more spherical in the radio emission, while those with larger scale lengths are flatter. The radio scale height depends mainly on the radio diameter of the galaxy. The sample galaxies are consistent with an escape-dominated radio halo with convective cosmic ray propagation, indicating that galactic winds are a widespread phenomenon in spiral galaxies. While a higher star formation rate or star formation surface density does not lead to a higher wind velocity, we found for the first time observational evidence of a gravitational deceleration of CRE outflow, e.g. a lowering of the wind velocity from the galactic disk.

It remains a mystery when our Milky Way first formed a stellar disk component that survived and maintained its disk structure from subsequent galaxy mergers. We present a study of the age-dependent structure and star formation rate of the Milky Way’s disk using high-α stars with substantial orbital angular momentum that have precise age determinations. Our results show that the radial scale length is nearly independent of age, whereas the vertical scale height experienced dramatic evolution. A disk-like geometry presents even for populations older than 13 Gyr, with the scale height-to-length ratio dropping below 0.5 for populations younger than 12.5 Gyr. We dub the oldest population that has maintained a disk geometry—apparently formed over 13 Gyr ago—PanGu. With an estimated present-day stellar mass of approximately 2 × 109 M⊙, PanGu is presumed to be a major stellar component of our Galaxy in the earliest epoch. The total present-day stellar mass of the whole high-α disk is 2 × 1010 M⊙, which was mostly formed during a distinct star formation rate peak of 11 M⊙ per year around 11 Gyr ago. A comparison with Milky Way analogues in the TNG50 simulations implies that our Galaxy has experienced an exceptionally quiescent dynamical history, even before the Gaia–Enceladus merger.

Dwarf Galaxies and Galaxy Halos

In order to test a possible evolutionary scenario of high-z compact quiescent galaxies (cQGs) that they can survive as local compact cores embedded in local massive galaxies with different morphology classes, we explore the star formation histories of local compact cores according to their spectral analysis. We build a sample of 182 massive galaxies with compact cores (M*, core > 1010.6 M⊙) at 0.02 ≤ z ≤ 0.06 from SDSS DR7 spectroscopic catalogue. The STARLIGHT package is used to analyze the median stacked spectra and derive the stellar ages and metallicities. Our main results show that local compact cores have the average age of about 12.1 ± 0.6 Gyr, indicating their early formation at z > 3, which is consistent with the formation redshifts of cQGs at 1 < z < 3. Together with previous studies, our result that local compact cores have similar formation redshifts as those of high-z cQGs, supports that local massive galaxies with compact cores are possible descendants of cQGs. Morphological study of local galaxies with compact cores suggests that there would be multiple possible evolutionary paths for high-z cQGs: most of them (> 80%) will evolve into local massive early-type galaxies according to dry minor merger, while some of them (∼ 15%) will build substantial stellar/gas discs according to the late-time gas accretion and sustaining star formation, and finally grow up into spiral galaxies.

The determination of the LSR is still a matter of debate. The classical value of the tangential peculiar motion of the Sun with respect to the LSR was challenged in recent years, claiming a significantly larger value. We show that the RAdial Velocity Experiment (RAVE) sample of dwarf stars is an excellent data set to derive tighter boundary conditions to chemodynamical evolution models of the extended solar neighbourhood. We present an improved Jeans analysis, which allows a better interpretation of the measured kinematics of stellar populations in the Milky Way disc. We propose an improved version of the Str\"omberg relation with the radial scalelengths as the only unknown. Binning RAVE stars in metallicity reveals a bigger asymmetric drift (corresponding to a smaller radial scalelength) for more metal-rich populations. With the standard assumption of velocity-dispersion independent radial scalelengths in each metallicity bin, we redetermine the LSR. The new Str\"omberg equation yields a joint LSR value of V_\sun=3.06 \pm 0.68 km/s, which is even smaller than the classical value based on Hipparcos data. The corresponding radial scalelength increases from 1.6 kpc for the metal-rich bin to 2.9 kpc for the metal-poor bin, with a trend of an even larger scalelength for young metal-poor stars. When adopting the recent Sch\"onrich value of V_\sun=12.24 km/s for the LSR, the new Str\"omberg equation yields much larger individual radial scalelengths of the RAVE subpopulations, which seem unphysical in part. The new Str\"omberg equation allows a cleaner interpretation of the kinematic data of disc stars in terms of radial scalelengths. Lifting the LSR value by a few km/s compared to the classical value results in strongly increased radial scalelengths with a trend of smaller values for larger velocity dispersions.

Context. The observation of total and linearly polarized synchrotron radiation of spiral galaxies in the radio continuum reveals the distribution and structure of their magnetic fields. By observing these, information about the proposed dynamo processes that preserve the large-scale magnetic fields in spiral galaxies can be gained. Additionally, by analyzing the synchrotron intensity, the transport processes of cosmic rays into the halo of edge-on spiral galaxies can be investigated. Aims. We analyze the magnetic field geometry and the transport processes of the cosmic rays of the edge-on spiral starburst galaxy NGC 4666 from CHANG-ES radio data in two frequencies; 6 GHz (C-band) and 1.5 GHz (L-band). Supplementary X-ray data are used to investigate the hot gas in NGC 4666. Methods. We determine the radio scale heights of total power emission at both frequencies for this galaxy. We show the magnetic field orientations derived from the polarization data. Using rotation measure (RM) synthesis we further study the behavior of the RM values along the disk in C-band to investigate the large-scale magnetic-field pattern. We use the revised equipartition formula to calculate a map of the magnetic field strength. Furthermore, we model the processes of cosmic-ray transport into the halo with the 1D SPINNAKER model. Results. The extended radio halo of NGC 4666 is box-shaped and is probably produced by the previously observed supernova-driven superwind. This is supported by our finding of an advective cosmic-ray transport such as that expected for a galactic wind. The scale-height analysis revealed an asymmetric halo above and below the disk as well as between the two sides of the major axis. A central point source as well as a bubble structure is seen in the radio data for the first time. Our X-ray data show a box-shaped hot halo around NGC 4666 and furthermore confirm the AGN nature of the central source. NGC 4666 has a large-scale X-shaped magnetic field in the halo, as has been observed in other edge-on galaxies. The analysis furthermore revealed that the disk of NGC 4666 shows hints of field reversals along its radius, which is the first detection of this phenomenon in an external galaxy.