Environments of Acitive Galaxies

This dissertation presents studies on the environments of active galaxies. Paper I is a case study of a cluster of galaxies containing BL Lac object RGB 1745+398. We measured the velocity dispersion, mass, and richness of the cluster. This was one of the most thorough studies of the environments of a BL Lac object. Methods used in the paper could be used in the future for studying other clusters as well. In Paper II we studied the environments of nearby quasars in the Sloan Digital Sky Survey (SDSS). We found that quasars have less neighboring galaxies than luminous inactive galaxies. In the large-scale structure, quasars are usually located at the edges of superclusters or even in void regions. We concluded that these low-redshift quasars may have become active only recently because the galaxies in low-density environments evolve later to the phase where quasar activity can be triggered. In Paper III we extended the analysis of Paper II to other types of AGN besides quasars. We found that different types of AGN have different large-scale environments. Radio galaxies are more concentrated in superclusters, while quasars and Seyfert galaxies prefer low-density environments. Different environments indicate that AGN have different roles in galaxy evolution. Our results suggest that activity of galaxies may depend on their environment on the large scale. Our results in Paper III raised questions of the cause of the environment-dependency in the evolution of galaxies. Because high-density large-scale environments contain richer groups and clusters than the underdense environments, our results could reflect smaller-scale effects. In Paper IV we addressed this problem by studying the group and supercluster scale environments of galaxies together. We compared the galaxy populations in groups of different richnesses in different large-scale environments. We found that the large-scale environment affects the galaxies independently of the group richness. Galaxies in low-density environments on the large scale are more likely to be star-forming than those in superclusters even if they are in groups with the same richness. Based on these studies, the conclusion of this dissertation is that the large-scale environment affects the evolution of galaxies. This may be caused by different “speed” of galaxy evolution in low and high-density environments: galaxies in dense environments reach certain phases of evolution earlier than galaxies in underdense environments. As a result, the low-density regions at low redshifts are populated by galaxies in earlier phases of evolution than galaxies in high-density regions.

Can retired galaxies mimic active galaxies? Clues from the Sloan Digital Sky Survey

The classification of galaxies as star forming or active is generally done in the ([O III]/H beta, [N II]/H alpha) plane. The Sloan Digital Sky Survey (SDSS) has revealed that, in this plane, the distribution of galaxies looks like the two wings of a seagull. Galaxies in the right wing are referred to as Seyfert/LINERs, leading to the idea that non-stellar activity in galaxies is a very common phenomenon. Here, we argue that a large fraction of the systems in the right wing could actually be galaxies which stopped forming stars. The ionization in these `retired` galaxies would be produced by hot post-asymptotic giant branch stars and white dwarfs. Our argumentation is based on a stellar population analysis of the galaxies via our STARLIGHT code and on photoionization models using the Lyman continuum radiation predicted for this population. The proportion of LINER galaxies that can be explained in such a way is, however, uncertain. We further show how observational selection effects account for the shape of the right wing. Our study suggests that nuclear activity may not be as common as thought. If retired galaxies do explain a large part of the seagull`s right wing, some of the work concerning nuclear activity in galaxies, as inferred from SDSS data, will have to be revised.

Discovery of "Warm Dust" galaxies in clusters at z ~ 0.3: evidence for stripping of cool dust in the dense environment?

Using far-infrared imaging from the "Herschel Lensing Survey," we derive dust properties of spectroscopically confirmed cluster member galaxies within two massive systems at z ~ 0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30 K). Several sub-luminous infrared galaxy (LIRG; L_IR < 10^11 L_☉) Bullet Cluster members are much warmer (T_dust > 37 K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (>3σ), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e., preferentially located in the cluster periphery, although the galaxy hosts tend toward lower stellar masses (M_* < 10^10 M_☉). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30%-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover an SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR luminosity or dust mass than the currently observed normal star-forming population.

Local luminous infrared galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs, L_IR=10^11-10 L_⨀) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a guaranteed time observation (GTO) Spitzer InfraRed Spectrograph (IRS) program aimed to obtain spectral mapping of a sample of 14 local d<76Mpc LIRGs. The data cubes map, at least, the central 20arcsec X 20arcsec to 30 arcsec X 30 arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38 μ m spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [Ne II] 12.81 μ m, [Ne III]15.56 μ m, [S III] 18.71 μ m, H_2 at 17 μ m), continuum, the 6.2 and 11.3 μ m polycyclic aromatic hydrocarbon (PAH) features, and the 9.7 μ m silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7 μ m silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of active galactic nuclei (AGN) indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6 9.7 μ m . We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

A sub-kiloparsec scale view of star formation in M31

Thesis (Ph.D.)--University of Washington, 2016-08

Nuclei of Nearby Active Galaxies: NIR and Sub-mm Views and Finalizing the LINC-NIRVANA Fringe and Flexure Tracking System: Flexure and Temperature Behavior

The work presented in my thesis addresses the two cornerstones of modern astronomy: Observation and Instrumentation. Part I deals with the observation of two nearby active galaxies, the Seyfert 2 galaxy NGC 1433 and the Seyfert 1 galaxy NGC 1566, both at a distance of $\sim10$ Mpc, which are part of the Nuclei of Galaxies (NUGA) sample. It is well established that every galaxy harbors a super massive black hole (SMBH) at its center. Furthermore, there seems to be a fundamental correlation between the stellar bulge and SMBH masses. Simulations show that massive feedback, e.g., powerful outflows, in Quasi Stellar Objects (QSOs) has an impact on the mutual growth of bulge and SMBH. Nearby galaxies follow this relation but accrete mass at much lower rates. This gives rise to the following questions: Which mechanisms allow feeding of nearby Active Galactic Nuclei (AGN)? Is this feeding triggered by events, e.g., star formation, nuclear spirals, outflows, on $\sim500$ pc scales around the AGN? Does feedback on these scales play a role in quenching the feeding process? Does it have an effect on the star formation close to the nucleus? To answer these questions I have carried out observations with the Spectrograph for INtegral Field Observation in the Near Infrared (SINFONI) at the Very Large Telescope (VLT) situated on Cerro Paranal in Chile. I have reduced and analyzed the recorded data, which contain spatial and spectral information in the H-band ($1.45 \mic-1.85 \mic$) and K-band ($1.95 \mic-2.45 \mic$) on the central $10\arcsec\times10\arcsec$ of the observed galaxies. Additionally, Atacama Large Millimeter/Sub-millimeter Array (ALMA) data at $350$ GHz ($\sim0.87$ mm) as well as optical high resolution Hubble Space Telescope (HST) images are used for the analysis. For NGC 1433 I deduce from comparison of the distributions of gas, dust, and intensity of highly ionized emission lines that the galaxy center lies $\sim70$ pc north-northwest of the prior estimate. A velocity gradient is observed at the new center, which I interpret as a bipolar outflow, a circum nuclear disk, or a combination of both. At least one dust and gas arm leads from a $r\sim200$ pc ring towards the nucleus and might feed the SMBH. Two bright warm H$_2$ gas spots are detected that indicate hidden star formation or a spiral arm-arm interaction. From the stellar velocity dispersion (SVD) I estimate a SMBH mass of $\sim1.74\times10^7$ \msol. For NGC 1566 I observe a nuclear gas disk of $\sim150$ pc in radius with a spiral structure. I estimate the total mass of this disk to be $\sim5.4\times10^7$ \msol. What mechanisms excite the gas in the disk is not clear. Neither can the existence of outflows be proven nor is star formation detected over the whole disk. On one side of the spiral structure I detect a star forming region with an estimated star formation rate of $\sim2.6\times10^{-3}$ \msol\ yr$^{-1}$. From broad Br$\gamma$ emission and SVD I estimate a mean SMBH mass of $\sim5.3\times10^6$ \msol\ with an Eddington ratio of $\sim2\times10^{-3}$. Part II deals with the final tests of the Fringe and Flexure Tracker (FFTS) for LBT INterferometric Camera and the NIR/Visible Adaptive iNterferometer for Astronomy (LINC-NIRVANA) at the Large Binocular Telescope (LBT) in Arizona, USA, which I conducted. The FFTS is the subsystem that combines the two separate beams of the LBT and enables near-infrared interferometry with a significantly large field of view. The FFTS has a cryogenic system and an ambient temperature system which are separated by the baffle system. I redesigned this baffle to guarantee the functionality of the system after the final tests in the Cologne cryostat. The redesign did not affect any scientific performance of LINC-NIRVANA. I show in the final cooldown tests that the baffle fulfills the temperature requirement and stays $<110$ K whereas the moving stages in the ambient system stay $>273$ K, which was not given for the old baffle design. Additionally, I test the tilting flexure of the whole FFTS and show that accurate positioning of the detector and the tracking during observation can be guaranteed.

Star formation history of Canis Major R1 I. Wide-Field X-ray study of the young stellar population

Aims. The CMa R1 star-forming region contains several compact clusters as well as many young early-B stars. It is associated with a well-known bright rimmed nebula, the nature of which is unclear (fossil HII region or supernova remnant). To help elucidate the nature of the nebula, our goal was to reconstruct the star-formation history of the CMa R1 region, including the previously unknown older, fainter low-mass stellar population, using X-rays. Methods. We analyzed images obtained with the ROSAT satellite, covering similar to 5 sq. deg. Complementary VRI photometry was performed with the Gemini South telescope. Colour-magnitude and colour-colour diagrams were used in conjunction with pre-main sequence evolutionary tracks to derive the masses and ages of the X-ray sources. Results. The ROSAT images show two distinct clusters. One is associated with the known optical clusters near Z CMa, to which similar to 40 members are added. The other, which we name the ""GU CMa"" cluster, is new, and contains similar to 60 members. The ROSAT sources are young stars with masses down to M(star) similar to 0.5 M(circle dot), and ages up to 10 Myr. The mass functions of the two clusters are similar, but the GU CMa cluster is older than the cluster around Z CMa by at least a few Myr. Also, the GU CMa cluster is away from any molecular cloud, implying that star formation must have ceased; on the contrary (as already known), star formation is very active in the Z CMa region.

Ultra-high-energy cosmic rays from centaurus A: jet interaction with gaseous shells

Ultra-high-energy cosmic rays (UHECRs), with energies above similar to 6 x 10(19) eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model, many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.

Substructure and Dynamics of the Fornax Cluster

We present the first dynamical analysis of a galaxy cluster to include a large fraction of dwarf galaxies. Our sample of 108 Fornax Cluster members measured with the UK Schmidt Telescope FLAIR-II spectrograph contains 55 dwarf galaxies (15.5 > b(j) > 18.0 or -16 > M-B > -13.5). H alpha emission shows that of the dwarfs are star forming, twice the fraction implied by morphological classifications. The total sample has a mean velocity of 1493 +/- 36 kms s(-1) and a velocity dispersion of 374 +/- 26 km s(-1). The dwarf galaxies form a distinct population: their velocity dispersion (429 +/- 41 km s(-1)) is larger than that of the giants () at the 98% confidence level. This suggests that the dwarf population is dominated by infalling objects whereas the giants are virialized. The Fornax system has two components, the main Fornax Cluster centered on NGC 1399 with cz = 1478 km s(-1) and sigma (cz) = 370 km s(-1) and a subcluster centered 3 degrees to the southwest including NGC 1316 with cz = 1583 km s(-1) and sigma (cz) = 377 km s(-1). This partition is preferred over a single cluster at the 99% confidence level. The subcluster, a site of intense star formation, is bound to Fornax and probably infalling toward the cluster core for the first time. We discuss the implications of this substructure for distance estimates of the Fornax Cluster. We determine the cluster mass profile using the method of Diaferio, which does not assume a virialized sample. The mass within a projected radius of 1.4 Mpc is (7 +/- 2) x 10(13) M-., and the mass-to-light ratio is 300 +/- 100 M-./L-.. The mass is consistent with values derived from the projected mass virial estimator and X-ray measurements at smaller radii.

The extraordinarily bright optical afterglow of GRB 991208 and its host galaxy

Observations of the extraordinarily bright optical afterglow (OA) of GRB 991208 started 2.1 d after the event. The flux decay constant of the OA in the R-band is -2.30 +/- 0.07 up to 5 d, which is very likely due to the jet effect, and after that it is followed by a much steeper decay with constant -3.2 +/- 0.2, the fastest one ever seen in a GRB OA. A negative detection in several all-sky films taken simultaneously to the event implies either a previous additional break prior to 2 d after the occurrence of the GRB (as expected from the jet effect). The existence of a second break might indicate a steepening in the electron spectrum or the superposition of two events. Once the afterglow emission vanished, contribution of a bright underlying SN is found, but the light curve is not sufficiently well sampled to rule out a dust echo explanation. Our determination of z = 0.706 indicates that GRB 991208 is at 3.7 Gpc, implying an isotropic energy release of 1.15 x 10E53 erg which may be relaxed by beaming by a factor > 100. Precise astrometry indicates that the GRB coincides within 0.2' with the host galaxy, thus given support to a massive star origin. The absolute magnitude is M_B = -18.2, well below the knee of the galaxy luminosity function and we derive a star-forming rate of 11.5 +/- 7.1 Mo/yr. The quasi-simultaneous broad-band photometric spectral energy distribution of the afterglow is determined 3.5 day after the burst (Dec 12.0) implying a cooling frequency below the optical band, i.e. supporting a jet model with p = -2.30 as the index of the power-law electron distribution.

NGC 2782: a merger remnant with young stars in its gaseous tidal tail

We have searched for young star-forming regions around the merger remnant NGC 2782. By using Galaxy Evolution Explorer far-ultraviolet and near-ultraviolet imaging and H i data we found seven ultraviolet sources, located at distances greater than 26 kpc from the centre of NGC 2782, and coinciding with its western H i tidal tail. These regions were resolved in several smaller systems when Gemini/Gemini multi-object spectrograph (GMOS) r-band images were used. We compared the observed colours to stellar population synthesis models and found that these objects have ages of similar to 1 to 11 Myr and masses ranging from 103.9 to 104.6 M circle dot. By using Gemini/GMOS spectroscopic data we confirm memberships and derive high metallicities for three of the young regions in the tail (12+log(O/H) = 8.74 +/- 0.20, 8.81 +/- 0.20 and 8.78 +/- 0.20). These metallicities are similar to the value presented by the nuclear region of NGC 2782 and also similar to the value presented for an object located close to the main body of NGC 2782. The high metallicities measured for the star-forming regions in the gaseous tidal tail of NGC 2782 could be explained if they were formed out of highly enriched gas which was once expelled from the centre of the merging galaxies when the system collided. An additional possibility is that the tail has been a nursery of a few generations of young stellar systems which ultimately polluted this medium with metals, further enriching the already pre-enriched gas ejected to the tail when the galaxies collided.

Millimeter characteristics of Galactic and extragalactic starbursts

The goal of the thesis was to compare Galactic and extragalactic star-forming environments, to understand whether we can derive scaling relations for the extragalactic star formation from our understanding of the Galactic star formation; and to analyze the effect of the angular resolution of the observations and the molecular tracer used in extragalactic studies. It is therefore necessary to perform observations at the same linear resolution in both cases. With this in mind we have started a program aimed at comparing properties of mini-starburst regions in our Galaxy and starbursts in nearby galaxies at similar linear resolutions.

Investigating evidence for different black hole accretion modes since redshift z ∼ 1

Chandra data in the COSMOS, AEGIS-XD and 4 Ms Chandra Deep Field South are combined with multiwavelength photometry available in those fields to determine the rest-frame U − V versus V − J colours of X-ray AGN hosts in the redshift intervals 0.1 < z < 0.6 (mean z¯=0.40) and 0.6 < z < 1.2 (mean z¯=0.85). This combination of colours provides an effective and least model-dependent means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasizes differences between AGN populations split by their U − V versus V − J colours. AGN in quiescent galaxies consist almost exclusively of bulges, while star-forming hosts are equally split between early- and late-type hosts. The position of AGN hosts on the U − V versusV − J diagram is then used to set limits on the accretion density of the Universe associated with evolved and star-forming systems independent of dust induced biases. It is found that most of the black hole growth at z≈ 0.40 and 0.85 is associated with star-forming hosts. Nevertheless, a non-negligible fraction of the X-ray luminosity density, about 15–20 per cent, at both z¯=0.40 and 0.85, is taking place in galaxies in the quiescent region of the U − V versus V − J diagram. For the low-redshift sub-sample, 0.1 < z < 0.6, we also find tentative evidence, significant at the 2σ level, that AGN split by their U − V and V − J colours have different Eddington ratio distributions. AGN in blue star-forming hosts dominate at relatively high Eddington ratios. In contrast, AGN in red quiescent hosts become increasingly important as a fraction of the total population towards low Eddington ratios. At higher redshift, z > 0.6, such differences are significant at the 2σ level only for sources with Eddington ratios ≳ 10^− 3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of supermassive black holes at the centres of galaxies. We compare these results with the predictions of theGALFORM semi-analytic model for the cosmological evolution of AGN and galaxies. This model postulates two black hole fuelling modes, the first is linked to star formation events and the second takes place in passive galaxies. GALFORM predicts that a substantial fraction of the black hole growth at z < 1 is associated with quiescent galaxies, in apparent conflict with the observations. Relaxing the strong assumption of the model that passive AGN hosts have zero star formation rate could bring those predictions in better agreement with the data.

[C II] and 12^CO(1-0) emission maps in HLSJ091828.6+514223: a strongly lensed interacting system AT z = 5.24

We present Submillimeter Array [C II] 158 μm and Karl G. Jansky Very Large Array 12^CO(1-0) line emission maps for the bright, lensed, submillimeter source at z = 5.2430 behind A 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple 12^CO rotational transitions, [C I], water, and [N II], providing some of the best constraints on the properties of the interstellar medium in a galaxy at z > 5. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8–1000 μm) = (1.6 ± 0.1) × 10^14 L_☉ μ^–1, where the total magnification μ_total = 8.9 ± 1.9, via a new lens model from the [C II] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components that correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source-plane separation <1 kpc. The reddest dominates the continuum map (demagnified L_FIR, component = (1.1 ± 0.2) × 10^13 L_☉) and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (ΔV ~ 500 km s^–1) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region and photodissociation region analysis suggests that it is lower density, cooler, and forming stars less vigorously than the other components. Strikingly, it has very strong [N II] emission, which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

Supernova 2013fc in a circumnuclear ring of a luminous infrared galaxy: The big brother of SN 1998S

We present photometric and spectroscopic observations of SN 2013fc, a bright type II supernova (SN) in a circumnuclear star-forming ring in the luminous infrared galaxy ESO 154-G010, observed as part of the Public ESO Spectroscopic Survey of Transient Objects. SN 2013fc is both photometrically and spectroscopically similar to the well-studied type IIn SN 1998S and to the bright type II-L SN 1979C. It exhibits an initial linear decline, followed by a short plateau phase and a tail phase with a decline too fast for ⁵⁶Co decay with full γ -ray trapping. Initially, the spectrum was blue and featureless. Later on, a strong broad (~8000 km s^-1) H α emission profile became prominent. We apply a STARLIGHT stellar population model fit to the SN location (observed when the SN had faded) to estimate a high extinction of A_V = 2.9 ± 0.2 mag and an age of 10⁺³ _-2 Myr for the underlying cluster.We compare the SN to SNe 1998S and 1979C and discuss its possible progenitor star considering the similarities to these events. With a peak brightness of B = -20.46 ± 0.21 mag, SN 2013fc is 0.9 mag brighter than SN 1998S and of comparable brightness to SN 1979C.We suggest that SN 2013fc was consistent with a massive red supergiant (RSG) progenitor. Recent mass loss probably due to a strong RSG wind created the circumstellar matter illuminated through its interaction with the SN ejecta. We also observe a near-infrared excess, possibly due to newly condensed dust.