Kinematics of local and high-z galaxies through 3D modeling of emission-line datacubes

The kinematics is a fundamental tool to infer the dynamical structure of galaxies and to understand their formation and evolution. Spectroscopic observations of gas emission lines are often used to derive rotation curves and velocity dispersions. It is however difficult to disentangle these two quantities in low spatial-resolution data because of beam smearing. In this thesis, we present 3D-Barolo, a new software to derive the gas kinematics of disk galaxies from emission-line data-cubes. The code builds tilted-ring models in the 3D observational space and compares them with the actual data-cubes. 3D-Barolo works with data at a wide range of spatial resolutions without being affected by instrumental biases. We use 3D-Barolo to derive rotation curves and velocity dispersions of several galaxies in both the local and the high-redshift Universe. We run our code on HI observations of nearby galaxies and we compare our results with 2D traditional approaches. We show that a 3D approach to the derivation of the gas kinematics has to be preferred to a 2D approach whenever a galaxy is resolved with less than about 20 elements across the disk. We moreover analyze a sample of galaxies at z~1, observed in the H-alpha line with the KMOS/VLT spectrograph. Our 3D modeling reveals that the kinematics of these high-z systems is comparable to that of local disk galaxies, with steeply-rising rotation curves followed by a flat part and H-alpha velocity dispersions of 15-40 km/s over the whole disks. This evidence suggests that disk galaxies were already fully settled about 7-8 billion years ago. In summary, 3D-Barolo is a powerful and robust tool to separate physical and instrumental effects and to derive a reliable kinematics. The analysis of large samples of galaxies at different redshifts with 3D-Barolo will provide new insights on how galaxies assemble and evolve throughout cosmic time.

X-ray emission from clusters and groups of galaxies

Recent major advances in x-ray imaging and spectroscopy of clusters have allowed the determination of their mass and mass profile out to ≈1/2 the virial radius. In rich clusters, most of the baryonic mass is in the gas phase, and the ratio of mass in gas/stars varies by a factor of 2–4. The baryonic fractions vary by a factor of ≈3 from cluster to cluster and almost always exceed 0.09 h50−[3/2] and thus are in fundamental conflict with the assumption of Ω = 1 and the results of big bang nucleosynthesis. The derived Fe abundances are 0.2–0.45 solar, and the abundances of O and Si for low redshift systems are 0.6–1.0 solar. This distribution is consistent with an origin in pure type II supernova. The amount of light and energy produced by these supernovae is very large, indicating their importance in influencing the formation of clusters and galaxies. The lack of evolution of Fe to a redshift of z ≈ 0.4 argues for very early enrichment of the cluster gas. Groups show a wide range of abundances, 0.1–0.5 solar. The results of an x-ray survey indicate that the contribution of groups to the mass density of the universe is likely to be larger than 0.1 h50−2. Many of the very poor groups have large x-ray halos and are filled with small galaxies whose velocity dispersion is a good match to the x-ray temperatures.

First science with SHARDS: emission line galaxies

SHARDS is an unbiased ultra-deep spectro-photometric survey with GTC@OSIRIS aimed at selecting and studying massive passively evolving galaxies at z=1.0-2.3 using a set of 24 medium-band filters (FWHM~17 nm) at 500-950 nm in GOODS-N. Our observing strategy is optimized to detect at z>1 the prominent Mg absorption feature at rest-frame ~ 280 nm, a distinctive, necessary, and sufficient feature of evolved stellar populations. Nonetheless, the data quality allow a plethora of studies on galaxy populations, including Emission Lines Galaxies (ELGs) about which we have started our first science verification project presented in this contribution.

Contribution from normal and starburst galaxies to the extragalactic gamma-ray background (EGRB)

The extragalactic diffuse emission at gamma-ray energies has interesting cosmological implications since these photons suffer little or no attenuation during their propagation from the site of origin. The emission could originate from either truly diffuse processes or from unresolved point sources such as AGNs, normal galaxies and starburst galaxies. Here, we examine the unresolved point source origin of the extragalactic gamma-ray background emission from normal galaxies and starburst galaxies. gamma-ray emission from normal galaxies is primarily coming from cosmic-ray interactions with interstellar matter and radiation (similar to 90%) along with a small contribution from discrete point sources (similar to 10%). Starburst galaxies are expected to have enhanced supernovae activity which leads to higher cosmic-ray densities, making starburst galaxies sufficiently luminous at gamma-ray energies to be detected by the current gamma-ray mission(Fermi Gamma-ray Space Telescope).

Generation of rotationally dominated galaxies by mergers of pressure-supported progenitors

Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure-supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depending on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on an orbit with high angular momentum. Systems with these hybrid properties have recently been observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with these observations and demonstrate that even mergers composed of non rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii.

An HI study of three long-tailed irregular galaxies in the cluster Abell 1367

We present the results on the distribution and kinematics of HI gas with higher sensitivity and in one case of higher spectral resolution as well than reported earlier, of three irregular galaxies CGCG 097073, 097079 and 097087 (UGC 06697) in the cluster Abell 1367. These galaxies are known to exhibit long (50 - 75 kpc) tails of radio continuum and optical emission lines (H alpha) pointing away from the cluster centre and arcs of starformation on the opposite sides of the tails, These features as well as the HI properties, with two of the galaxies (CGCG 097073 and 097079) exhibiting sharper gradients in HI intensity on the side of the tails, are consistent with the HI gas being affected by the ram pressure of the intracluster medium. However the HI emission in all the three galaxies extends to much smaller distances than the radio-continuum and H alpha tails, and are possibly still bound to the parent galaxies. Approximately 20 - 30 per cent of the HI mass is seen to accumulate on the downstream side due to the effects of ram pressure.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

Tidally compressed gas in centers of early-type and ultraluminous galaxies

In this paper we propose that the compressive tidal held in the centers of flat-core early-type galaxies and ultraluminous galaxies compresses molecular clouds producing dense gas observed in the centers of these galaxies. The effect of galactic tidal fields is usually considered disruptive in the literature. However, for some galaxies, the mass profile flattens toward the center and the resulting galactic tidal field is not disruptive, but instead it is compressive within the flat-core region. We have used the virial theorem to determine the minimum density of a molecular cloud to be stable and gravitationally bound within the tidally compressive region of a galaxy. We have applied the mechanism to determine the mean molecular cloud densities in the centers of a sample of flat-core, early-type galaxies and ultraluminous galaxies. For early-type galaxies with a core-type luminosity profile, the tidal held of the galaxy is compressive within half the core radius. We have calculated the mean gas densities for molecular gas in a sample of early-type galaxies which have already been detected in CO emission, and we obtain mean densities of [n] similar to 10(3)-10(6) cm(-3) within the central 100 pc radius. We also use our model to calculate the molecular cloud densities in the inner few hundred parsecs of a sample of ultraluminous galaxies. From the observed rotation curves of these galaxies we show that they have a compressive core within their nuclear region. Our model predicts minimum molecular gas densities in the range 10(2)-10(4) cm(-3) in the nuclear gas disks; the smaller values are applicable typically for galaxies with larger core radii. The resulting density values agree well with the observed range. Also, for large core radii, even fairly low-density gas (similar to 10(2) cm(-3)) can remain bound and stable close to the galactic center.

Search for sub-mm, mm and radio continuum emission from extremely red objects

We present the results of sub-mm, mm (850 mum, 450 mum and 1250 mum) and radio (1.4 and 4.8 GHz) continuum observations of a sample of 27 K-selected Extremely Red Objects, or EROs, (14 of which form a complete sample with K < 20 and I - K > 5) aimed at detecting dusty starbursts, deriving the fraction of UltraLuminous Infrared Galaxies (ULIGs) in ERO samples, and constraining their redshifts using the radio-FIR correlation. One ERO was tentatively detected at 1250 mum and two were detected at 1.4 GHz, one of which has a less secure identification as an ERO counterpart. Limits on their redshifts and their star forming properties are derived and discussed. We stacked the observations of the undetected objects at 850 mum, 1250 mum and 4.8 GHz in order to search for possible statistical emission from the ERO population as a whole, but no significant detections were derived either for the whole sample or as a function of the average NIR colours. These results strongly suggest that the dominant population of EROs with K < 20 is not comprised of ULIGs like HR 10, but is probably made of radio-quiet ellipticals and weaker starburst galaxies with L < 10(12) L . and SFR < few hundred M. yr(-1).

Luminosity profiles of advanced mergers of galaxies using 2MASS data

A sample of 27 disturbed galaxies that show signs of interaction but have a single nucleus were selected from the Arp and the Arp-Madore catalogues. For these, the Ks band images from the Two Micron All Sky Survey (2MASS) are analysed to obtain their radial luminosity pro�les and other structural parameters. We �nd that in spite of their similar optical appearance, the sample galaxies vary in their dynamical properties, and fall into two distinct classes. The �rst class consists of galaxies which can be described by a single r1=4 law and the second class consists of galaxies that show an outer exponential disk. A few galaxies that have disturbed pro�les cannot be �t into either of the above classes. However, all the galaxies are similar in all other parameters such as the far-infrared colours, the molecular hydrogen content and the central velocity dispersion. Thus, the dynamical parameters of these sets seem to be determined by the ratio of the initial masses of the colliding galaxies. We propose that the galaxies in the �rst class result from a merger of spiral galaxies of equal masses whereas the second class of galaxies results from a merger of unequal mass galaxies. The few objects that do not fall into either category show a disturbed luminosity pro�le and a wandering centre, which is indicative of these being unrelaxed mergers. Of the 27 galaxies in our sample, 9 show elliptical-like pro�les and 13 show an outer exponential. Interestingly, Arp 224, the second oldest merger remnant of the Toomre sequence shows an exponential disk in the outer parts.

Origin of Disk Lopsidedness in Spiral Galaxies

In our work we have used the atomic hydrogen [HΙ] gas distribution in the HΙ 21-cm line emission to study the dark matter halo perturbations. For tHΙs analysis, the 2-D HΙ surface density and velocity maps (arcHΙval) of the galaxies in the Eridanus group (obtained using the GMRT) and in the Ursa Major group (obtained from WSRT) were used. In addition a few HΙckson Compact Groups of galaxies were also studied using the GMRT. The HΙ maps of these galaxies were Fourier analysed to estimate the asymmetry in the distribution and motion of gas. The average asymmetry parameter in the 1.5 to 2.5 K′-band scale lengths was found to be ~ 0.27 for the Eridanus group of galaxies wHΙle it was ~ 0.14 for the Ursa Major group of galaxies. The asymmetries in the distribution of HΙ as a function of Hubble type of galaxies were also studied and was found to be directly correlated with the compactness of the groups. In addition, the trend in the asymmetry as a function of the Hubble type of galaxies was opposite to that seen in the field galaxies, i.e., in the group galaxies, the early type galaxies showed more asymmetry than late type. These two aspects indicated that tidal interactions between the galaxies in a group environment to be the major cause of asymmetries. The observed asymmetry parameters were consistent with recent numerical simulations of asymmetries of gas disk caused by fly-by interactions. We have also estimated the perturbation of dark matter halo using the asymmetry parameter obtained from the Fourier series analysis of the surface density maps.

Long way to go: how outflows from large galaxies propagate through the hot halo gas

Using hydrodynamic simulations, we study the mass-loss due to supernova-driven outflows from Milky Way type disc galaxies, paying particular attention to the effect of the extended hot halo gas. We find that the total mass-loss at inner radii scales roughly linearly with total mass of stars formed, and that the mass loading factor at the virial radius can be several times its value at inner radii because of the swept up hot halo gas. The temperature distribution of the outflowing material in the inner region (similar to 10 kpc) is bimodal in nature, peaking at 10(5) K and 10(6.5) K, responsible for optical and X-ray emission, respectively. The contribution of cold/warm gas with temperature <= 10(5.5) K to the outflow rate within 10 kpc is approximate to 0.3-0.5. The warm mass loading factor, eta(3e5) (T <= 3 x 10(5) K) is related to the mass loading factor at the virial radius (eta(v)) as eta(v) approximate to 25 eta(3e5) (SFR/M-circle dot yr(-1))(-0.15) for a baryon fraction of 0.1 and a starburst period of 50 Myr. We also discuss the effect of multiple bursts that are separated by both short and long periods. The outflow speed at the virial radius is close to the sound speed in the hot halo, less than or similar to 200 km s(-1). We identify two `sequences' of outflowing cold gas at small scales: a fast (approximate to 500 km s(-1)) sequence, driven by the unshocked free-wind; and a slow sequence (approximate to +/- 100 km s(-1)) at the conical interface of the superwind and the hot halo.

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction (f(esc)) of ionizing photons from starburst galaxies. Using 2D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the centre of the disc along different angles through the superbubble and the gas disc. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of similar to 40 degrees, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scaleheights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed similar to1 - cos (1 rad)] = 0.5 from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time-and angle-averaged escape fraction on the mid-plane disc gas density (in the range n(0) = 0.15-50 cm(-3)) and the disc scaleheight (between z(0) = 10 and 600 pc). We find that the escape fraction is related to the disc parameters (the mid-plane disc density and scaleheight) roughly so that f(esc)(alpha)n(0)(2)z(0)(3) (with alpha approximate to 2.2) is a constant. For discs with a given warm neutral medium temperature, massive discs have lower escape fraction than low-mass galaxies. For Milky Way ISM parameters, we find f(esc) similar to 5 per cent, and it increases to approximate to 10 per cent for a galaxy 10 times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the Universe.

Detailed properties of high redshift galaxies

Galaxies evolve throughout the history of the universe from the first star-forming sources, through gas-rich asymmetric structures with rapid star formation rates, to the massive symmetrical stellar systems observed at the present day. Determining the physical processes which drive galaxy formation and evolution is one of the most important questions in observational astrophysics. This thesis presents four projects aimed at improving our understanding of galaxy evolution from detailed measurements of star forming galaxies at high redshift.

We use resolved spectroscopy of gravitationally lensed z ≃ 2 - 3 star forming galaxies to measure their kinematic and star formation properties. The combination of lensing with adaptive optics yields physical resolution of ≃ 100 pc, sufficient to resolve giant Hii regions. We find that ~ 70 % of galaxies in our sample display ordered rotation with high local velocity dispersion indicating turbulent thick disks. The rotating galaxies are gravitationally unstable and are expected to fragment into giant clumps. The size and dynamical mass of giant Hii regions are in agreement with predictions for such clumps indicating that gravitational instability drives the rapid star formation. The remainder of our sample is comprised of ongoing major mergers. Merging galaxies display similar star formation rate, morphology, and local velocity dispersion as isolated sources, but their velocity fields are more chaotic with no coherent rotation.

We measure resolved metallicity in four lensed galaxies at z = 2.0 − 2.4 from optical emission line diagnostics. Three rotating galaxies display radial gradients with higher metallicity at smaller radii, while the fourth is undergoing a merger and has an inverted gradient with lower metallicity at the center. Strong gradients in the rotating galaxies indicate that they are growing inside-out with star formation fueled by accretion of metal-poor gas at large radii. By comparing measured gradients with an appropriate comparison sample at z = 0, we demonstrate that metallicity gradients in isolated galaxies must flatten at later times. The amount of size growth inferred by the gradients is in rough agreement with direct measurements of massive galaxies. We develop a chemical evolution model to interpret these data and conclude that metallicity gradients are established by a gradient in the outflow mass loading factor, combined with radial inflow of metal-enriched gas.

We present the first rest-frame optical spectroscopic survey of a large sample of low-luminosity galaxies at high redshift (L < L*, 1.5 < z < 3.5). This population dominates the star formation density of the universe at high redshifts, yet such galaxies are normally too faint to be studied spectroscopically. We take advantage of strong gravitational lensing magnification to compile observations for a sample of 29 galaxies using modest integration times with the Keck and Palomar telescopes. Balmer emission lines confirm that the sample has a median SFR ∼ 10 M_sun yr^−1 and extends to lower SFR than has been probed by other surveys at similar redshift. We derive the metallicity, dust extinction, SFR, ionization parameter, and dynamical mass from the spectroscopic data, providing the first accurate characterization of the star-forming environment in low-luminosity galaxies at high redshift. For the first time, we directly test the proposal that the relation between galaxy stellar mass, star formation rate, and gas phase metallicity does not evolve. We find lower gas phase metallicity in the high redshift galaxies than in local sources with equivalent stellar mass and star formation rate, arguing against a time-invariant relation. While our result is preliminary and may be biased by measurement errors, this represents an important first measurement that will be further constrained by ongoing analysis of the full data set and by future observations.

We present a study of composite rest-frame ultraviolet spectra of Lyman break galaxies at z = 4 and discuss implications for the distribution of neutral outflowing gas in the circumgalactic medium. In general we find similar spectroscopic trends to those found at z = 3 by earlier surveys. In particular, absorption lines which trace neutral gas are weaker in less evolved galaxies with lower stellar masses, smaller radii, lower luminosity, less dust, and stronger Lyα emission. Typical galaxies are thus expected to have stronger Lyα emission and weaker low-ionization absorption at earlier times, and we indeed find somewhat weaker low-ionization absorption at higher redshifts. In conjunction with earlier results, we argue that the reduced low-ionization absorption is likely caused by lower covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. This result has important implications for the hypothesis that early galaxies were responsible for cosmic reionization. We additionally show that fine structure emission lines are sensitive to the spatial extent of neutral gas, and demonstrate that neutral gas is concentrated at smaller galactocentric radii in higher redshift galaxies.

The results of this thesis present a coherent picture of galaxy evolution at high redshifts 2 ≲ z ≲ 4. Roughly 1/3 of massive star forming galaxies at this period are undergoing major mergers, while the rest are growing inside-out with star formation occurring in gravitationally unstable thick disks. Star formation, stellar mass, and metallicity are limited by outflows which create a circumgalactic medium of metal-enriched material. We conclude by describing some remaining open questions and prospects for improving our understanding of galaxy evolution with future observations of gravitationally lensed galaxies.