Showing 1–50 of 50 results for author: Choi, E

The Origins of Gas Accreted by Supermassive Black Holes: the Importance of Recycled Gas

Authors: Ena Choi, Rachel S. Somerville, Jeremiah P. Ostriker, Michaela Hirschmann, Thorsten Naab

Abstract: We investigate the fueling mechanisms of supermassive black holes (SMBHs) by analyzing ten zoom-in cosmological simulations of massive galaxies, with stellar masses $10^{11-12} M_{\odot}$ and SMBH masses $10^{8.9-9.7}$ at $z=0$ and featuring various major and minor merger events. By tracing the gas history in these simulations, we categorize the gas accreted by the central SMBHs based on its origi… ▽ More We investigate the fueling mechanisms of supermassive black holes (SMBHs) by analyzing ten zoom-in cosmological simulations of massive galaxies, with stellar masses $10^{11-12} M_{\odot}$ and SMBH masses $10^{8.9-9.7}$ at $z=0$ and featuring various major and minor merger events. By tracing the gas history in these simulations, we categorize the gas accreted by the central SMBHs based on its origin. Gas that belonged to a different galaxy before accretion onto the BH is labeled as (1) ``external," while smoothly accreted cosmic gas is classified as (2) ``smooth." Gas produced within the primary halo through stellar evolution and subsequently accreted by the SMBH is classified as (3) ``recycled." Our analysis, which included stellar feedback, reveals that the primary fuel source for SMBHs is the recycled gas from dying stars. This recycled gas from stars in the inner region of the galaxy readily collapses toward the center, triggering starbursts, and simultaneously fueling the SMBH. Galaxy mergers also play a crucial role in fueling SMBHs in massive galaxies as SMBHs in massive halos tend to accrete a higher fraction of external gas from mergers compared to smoothly accreted gas. However, on average, it takes approximately 1.85 Gyr for external gas to enter the main galaxy and accrete onto the SMBH. Considering the presence of various other gas triggers for AGN activity alongside this time delay, the association between AGN and mergers may not always be obvious. △ Less

Submitted 19 February, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 14 pages, 7 figures, Accepted for publication in ApJ

On the Effects of Local Environment on Active Galactic Nucleus (AGN) in the Horizon Run 5 Simulation

Authors: Ankit Singh, Changbom Park, Ena Choi, Juhan Kim, Hyunsung Jun, Brad K. Gibson, Yonghwi Kim, Jaehyun Lee, Owain Snaith

Abstract: We use the Horizon Run 5 cosmological simulation to study the effect of galaxy intrinsic properties and the local environment on AGNs characterized by their threshold of the accretion rate. We select galaxies in the stellar mass range $10^{9.5} \le M^{}{*}/M^{}{\odot} \le 10^{10.5}$ in the snapshot at redshift $z$=0.625. Among various intrinsic properties, we find that the star formation rate of t… ▽ More We use the Horizon Run 5 cosmological simulation to study the effect of galaxy intrinsic properties and the local environment on AGNs characterized by their threshold of the accretion rate. We select galaxies in the stellar mass range $10^{9.5} \le M^{}{*}/M^{}{\odot} \le 10^{10.5}$ in the snapshot at redshift $z$=0.625. Among various intrinsic properties, we find that the star formation rate of the host galaxy is most correlated to the AGN activity. To quantify the environment, we use background galaxy number density (large-scale environment) and distance and morphological type of the nearest neighbors (small-scale environment), and study their relative effects on the AGN properties. We find that, compared to the background density, the nearest neighbor environment is the dominant quantity determining the bolometric luminosity, star formation rate, and kinematic properties of AGNs and better dictates the gas mass of the host galaxy. We show that the cold gas content in the host galaxies is crucial in triggering AGN activity. However, when the nearest neighbor environment effects start to act at the neighbor distance of less than about half the virial radius of the neighbor, the neighbor environmental effects are the most dominant factor for quasar activity. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 25 pages, 18 figures, 1 table, published in the the Astrophysical Journal (ApJ)

Emission-line properties of IllustrisTNG galaxies: from local diagnostic diagrams to high-redshift predictions for JWST

Authors: Michaela Hirschmann, Stephane Charlot, Anna Feltre, Emma Curtis-Lake, Rachel S. Somerville, Jacopo Chevallard, Ena Choi, Dylan Nelson, Christophe Morisset, Adele Plat, Alba Vidal-Garcia

Abstract: We compute synthetic, rest-frame optical and ultraviolet (UV) emission-line properties of galaxy populations at redshifts from z$\approx$0 to z=8 in a full cosmological framework. We achieve this by coupling, in post-processing, the cosmological IllustrisTNG simulations with new-generation nebular-emission models, accounting for line emission from young stars, post-asymptotic-giant-branch (PAGB) s… ▽ More We compute synthetic, rest-frame optical and ultraviolet (UV) emission-line properties of galaxy populations at redshifts from z$\approx$0 to z=8 in a full cosmological framework. We achieve this by coupling, in post-processing, the cosmological IllustrisTNG simulations with new-generation nebular-emission models, accounting for line emission from young stars, post-asymptotic-giant-branch (PAGB) stars, accreting black holes (BHs) and, for the first time, fast radiative shocks. The optical emission-line properties of simulated galaxies dominated by different ionizing sources are largely consistent with those expected from classical diagnostic diagrams and reflect the observed increase in [OIII]/H$β$ at fixed [NII]/H$α$ and the evolution of the H$α$, [OIII]$\lambda5007$ and [OII]$\lambda3727$ luminosity functions from z$\approx$0 to z$\sim$2. At higher redshift, we find that the emission-line galaxy population is dominated by star-forming and active galaxies, with negligible fractions of shock- and PAGB-dominated galaxies. We highlight 10 UV-diagnostic diagrams able to robustly identify the dominant ionizing sources in high-redshift galaxies. We also compute the evolution of several optical- and UV-line luminosity functions from z=4 to z=7, and the number of galaxies expected to be detectable per field of view in deep, medium-resolution spectroscopic observations with the NIRSpec instrument on board the James Webb Space Telescope. We find that 2-hour-long exposures are sufficient to achieve unbiased censuses of H$α$ and [OIII]$\lambda5007$ emitters, while at least 5 hours are required for H$β$, and even 10 hours will detect only progressively smaller fractions of [OII]$\lambda3727$, OIII]$\lambda1663$, CIII]$\lambda1908$, CIV$\lambda1550$, [NII]$\lambda6584$, SiIII]$\lambda1888$ and HeII$\lambda1640$ emitters, especially in the presence of dust. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 28 pages, 15 figures, submitted to MNRAS

Active galactic nucleus jet feedback in hydrostatic halos

Authors: Rainer Weinberger, Kung-Yi Su, Kristian Ehlert, Christoph Pfrommer, Lars Hernquist, Greg L. Bryan, Volker Springel, Yuan Li, Blakesley Burkhart, Ena Choi, Claude-André Faucher-Giguère

Abstract: Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct… ▽ More Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties. △ Less

Submitted 5 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

The Metallicity Distribution Function in Outer Halo Fields of Simulated Elliptical Galaxies Compared to Observations of NGC 5128

Authors: Ena Choi, Jeremiah P. Ostriker, Michaela Hirschmann, Rachel S. Somerville, Thorsten Naab

Abstract: Stellar metallicity distribution functions (MDF) have been measured for resolved stellar populations in the outer halos of many galaxies in nearby groups. Among them, the MDF of NGC 5128, the central giant elliptical in the Centaurus group, provides essential constraints for theories of massive galaxy formation and hierarchical assembly. To investigate the formation and chemical evolution history… ▽ More Stellar metallicity distribution functions (MDF) have been measured for resolved stellar populations in the outer halos of many galaxies in nearby groups. Among them, the MDF of NGC 5128, the central giant elliptical in the Centaurus group, provides essential constraints for theories of massive galaxy formation and hierarchical assembly. To investigate the formation and chemical evolution history of the outer halo of giant elliptical galaxies, we examine the chemical properties of three zoom-in high resolution cosmological hydrodynamical simulations of an NGC 5128-like giant elliptical galaxy and compare their outer halo MDFs to the observed one of NGC 5128. Even though the simulated galaxies have different merging histories and age distributions, all predicted MDFs are in good qualitative agreement with the observed one. The median metallicity of the simulated galaxies is on average $\rm [M/H]=-0.41 \pm 0.06$ compared to the observed value of $\rm [M/H]=-0.38 \pm 0.02$ for NGC 5128, and the dispersion in metallicity is $\sim 0.77$ dex for both observed and simulated galaxies. We investigate the origin of the stars ending up in the outer halo field of simulated galaxies and show that most have an `accreted' origin, formed in other small galaxies and later accreted in mergers. Only $\sim 15$ percent of the stars are formed `in situ' within the main progenitor of galaxy and radially migrate outwards. We show that the contribution of metal-rich in situ stars is sub-dominant in the outer halos of our simulated galaxies, but can be prominent in the inner regions. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 16 pages, 9 figures, Accepted for publication in ApJ

IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust

Authors: ChangHoon Hahn, Tjitske K. Starkenburg, Daniel Anglés-Alcázar, Ena Choi, Romeel Davé, Claire Dickey, Kartheik G. Iyer, Ariyeh H. Maller, Rachel S. Somerville, Jeremy L. Tinker, L. Y. Aaron Yung

Abstract: We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitud… ▽ More We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 26 pages, 15 figures

The Connection between Mergers and AGN Activity in Simulated and Observed Massive Galaxies

Authors: Ray S. Sharma, Ena Choi, Rachel S. Somerville, Gregory F. Snyder, Dale D. Kocevski, Michaela Hirschmann, Benjamin P. Moster, Thorsten Naab, Desika Narayanan, Jeremiah P. Ostriker, David J. Rosario

Abstract: We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthe… ▽ More We analyze a suite of $30$ high resolution zoom-in cosmological hydrodynamic simulations of massive galaxies with stellar masses $M_{\ast} > 10^{10.9} M_\odot$, with the goal of better understanding merger activity in AGN, AGN activity in merging systems, SMBH growth during mergers, and the role of gas content. Using the radiative transfer code \textsc{Powderday}, we generate HST-WFC3 F160W synthetic observations of redshift $0.5 < z < 3$ central galaxies, add noise properties similar to the CANDELS survey, and measure morphological properties from the synthetic images using commonly adopted non-parametric statistics. We compare the distributions of morphological properties measured from the synthetic images with a sample of inactive galaxies and X-ray selected AGN hosts from CANDELS. We study the connection between mergers and AGN activity in the simulations, the synthetic images, and the observed CANDELS sample. We find that, in both the simulations and CANDELS, even the most luminous $(L_{\rm bol} > 10^{45}$ erg s$^{-1})$ AGN in our sample are no more likely than inactive galaxies $(L_{\rm bol} < 10^{43}$ erg s$^{-1})$ to be found in merging systems. We also find that AGN activity is not overall enhanced by mergers, nor enhanced at any specific time in the $1$ Gyr preceding and following a merger. Even gas rich major mergers (stellar mass ratio $>$1:4) do not necessarily enhance AGN activity or significantly grow the central SMBH. We conclude that in the simulated massive galaxies studied here, mergers are not the primary drivers of AGN. △ Less

Submitted 18 December, 2023; v1 submitted 5 January, 2021; originally announced January 2021.

IQ Collaboratory II: The Quiescent Fraction of Isolated, Low Mass Galaxies Across Simulations and Observations

Authors: Claire M Dickey, Tjitske K Starkenburg, Marla Geha, ChangHoon Hahn, Daniel Anglés-Alcázar, Ena Choi, Romeel Davé, Shy Genel, Kartheik G Iyer, Ariyeh H Maller, Nir Mandelker, Rachel S Somerville, L Y Aaron Yung

Abstract: We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each s… ▽ More We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range $\mathrm{M}_* = 10^{9-10} \ \mathrm{M}_\odot$, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below $\mathrm{M}_*=10^{9} \ \mathrm{M}_\odot$. We find that the low mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite apparent variation in the late time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path towards more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 19 pages, 8 figures. Figure 4 presents the main result. Code used in this work may be accessed at github.com/IQcollaboratory/orchard. Submitted to ApJ

The Impact of Outflows driven by Active Galactic Nuclei on Metals in and around Galaxies

Authors: Ena Choi, Ryan Brennan, Rachel S. Somerville, Jeremiah P. Ostriker, Michaela Hirschmann, Thorsten Naab

Abstract: Metals in the hot gaseous halos of galaxies encode the history of star formation as well as the feedback processes that eject metals from the galaxies. X-ray observations suggest that massive galaxies have extended distributions of metals in their gas halos. We present predictions for the metal properties of massive galaxies and their gaseous halos from recent high resolution zoom-in simulations t… ▽ More Metals in the hot gaseous halos of galaxies encode the history of star formation as well as the feedback processes that eject metals from the galaxies. X-ray observations suggest that massive galaxies have extended distributions of metals in their gas halos. We present predictions for the metal properties of massive galaxies and their gaseous halos from recent high resolution zoom-in simulations that include mechanical and radiation driven feedback from Active Galactic Nuclei (AGN). In these simulations, AGN launch high-velocity outflows, mimicking observed broad absorption line winds. By comparing two sets of simulations with and without AGN feedback, we show that our prescription for AGN feedback is capable of driving winds and enriching halo gas `inside-out' by spreading centrally enriched metals to the outskirts of galaxies, into the halo and beyond. The metal (iron) profiles of halos simulated with AGN feedback have a flatter slope than those without AGN feedback, consistent with recent X-ray observations. The predicted gas iron abundance of group scale galaxies simulated with AGN feedback is $Z_{\rm Fe} = 0.23$ $Z_{\rm Fe,\odot}$ at $0.5 r_{500}$, which is 2.5 times higher than that in simulations without AGN feedback. In these simulations, AGN winds are also important for the metal enrichment of the intergalactic medium, as the vast majority of metals ejected from the galaxy by AGN-driven winds end up beyond the halo virial radius. △ Less

Submitted 21 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ. 16 pages, 9 figures. Key figures are 8 & 9

First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations

Authors: Drummond B. Fielding, Stephanie Tonnesen, Daniel DeFelippis, Miao Li, Kung-Yi Su, Greg L. Bryan, Chang-Goo Kim, John C. Forbes, Rachel S. Somerville, Nicholas Battaglia, Evan E. Schneider, Yuan Li, Ena Choi, Christopher C. Hayward, Lars Hernquist

Abstract: We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolu… ▽ More We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

Powderday: Dust Radiative Transfer for Galaxy Simulations

Authors: Desika Narayanan, Matthew J. Turk, Thomas Robitaille, Ashley J. Kelly, B. Connor McClellan, Ray S. Sharma, Prerak Garg, Matthew Abruzzo, Ena Choi, Charlie Conroy, Benjamin D. Johnson, Benjamin Kimock, Qi Li, Christopher C. Lovell, Sidney Lower, George C. Privon, Jonathan Roberts, Snigdaa Sethuram, Gregory F. Snyder, Robert Thompson, John H. Wise

Abstract: We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time fle… ▽ More We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ; comments welcome

Ionized gas kinematics of massive elliptical galaxies in CALIFA and in cosmological zoom-in simulations

Authors: Jan Florian, Bodo Ziegler, Michaela Hirschmann, Polychronis Papaderos, Ena Choi, Matteo Frigo, Jean-Michel Gomes, Rachel S. Somerville

Abstract: (Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmologica… ▽ More (Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall. △ Less

Submitted 13 December, 2019; originally announced December 2019.

Comments: Figure resolution reduced due to size limitations. Submitted to A&A 02.08.2019, referee report received 31.08.2019, resubmitted 19.11.2019, accepted 13.12.2019

Journal ref: A&A 635, A41 (2020)

AGN-driven quenching of satellite galaxies

Authors: Gohar Dashyan, Ena Choi, Rachel S. Somerville, Thorsten Naab, Amanda C. N. Quirk, Michaela Hirschmann, Jeremiah P. Ostriker

Abstract: We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stell… ▽ More We explore the effect of active galactic nucleus (AGN) feedback from central galaxies on their satellites by comparing two sets of cosmological zoom-in runs of 27 halos with masses ranging from $10^{12}$ to $10^{13.4}$ solar masses at z=0, with (wAGN) and without (noAGN) AGN feedback. Both simulations include stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, AGB stars, radiative heating within Strömgren spheres and photoelectric heating. Our wAGN model is identical to the noAGN model except that it also includes a model for black hole seeding and accretion, as well as AGN feedback via high-velocity broad absorption line winds and Compton/photoionization heating. We show that the inclusion of AGN feedback from the central galaxy significantly affects the star formation history and the gas content of the satellite galaxies. AGN feedback starts to affect the gas content and the star formation of the satellites as early as z=2. The mean gas rich fraction of satellites at z=0 decreases from 15% in the noAGN simulation to 5% in the wAGN simulation. The difference between the two sets extends as far out as five times the virial radius of the central galaxy at z=1. We investigate the quenching mechanism by studying the physical conditions in the surroundings of pairs of satellites matched across the wAGN and noAGN simulations and find an increase in the temperature and relative velocity of the intergalactic gas. △ Less

Submitted 18 June, 2019; originally announced June 2019.

Comments: 15 pages, 12 figures; accepted for publication in MNRAS

MIND THE GAP: Is The Too Big To Fail Problem Resolved?

Authors: Jeremiah P. Ostriker, Ena Choi, Anthony Chow, Kundan Guha

Abstract: The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function… ▽ More The faintness of satellite systems in galaxy groups has contributed to the widely discussed "missing satellite" and "too big to fail" issues. Using techniques based on Tremaine & Richstone (1977), we show that there is no problem with the luminosity function computed from modern codes per se, but that the gap between first and second brightest systems is too big {\it given} the luminosity function, that the same large gap is found in modern, large scale baryonic $Λ$CDM simulations such as EAGLE and IllustrisTNG, is even greater in dark matter only simulations, and finally, that this is most likely due to gravitationally induced merging caused by classical dynamical friction. Quantitatively the gap is larger in the computed simulations than in the randomized ones by $1.79 \pm 1.04$, $1.51 \pm 0.93$, $3.43 \pm 1.44$ and $3.33 \pm 1.35$ magnitudes in the EAGLE, IllustrisTNG, and dark matter only simulations of EAGLE and IllustrisTNG respectively. Furthermore the anomalous gaps in the simulated systems are even larger than in the real data by over half a magnitude and are still larger in the dark matter only simulations. Briefly stated, $Λ$CDM does not have a problem with an absence of "too big to fail" galaxies. Statistically significant large gaps between first and second brightest galaxies are to be expected. △ Less

Submitted 16 July, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 9 pages, 8 figures, Accepted for publication in the Astrophysical Journal

RAMSES II - RAMan Search for Extragalactic Symbiotic Stars. Project concept, commissioning, and early results from the science verification phase

Authors: R. Angeloni, D. R. Gonçalves, S. Akras, G. Gimeno, R. Diaz, J. Scharwächter, N. E. Nuñez, G. J. M. Luna, H. W. Lee, J. E. Heo, A. B. Lucy, M. Jaque Arancibia, C. Moreno, E. Chirre, S. J. Goodsell, P. Soto King, J. L. Sokoloski, B. E. Choi, M. Dias Ribeiro

Abstract: Symbiotic stars (SySts) are long-period interacting binaries composed of a hot compact star, an evolved giant star, and a tangled network of gas and dust nebulae. They represent unique laboratories for studying a variety of important astrophysical problems, and have also been proposed as possible progenitors of SNIa. Presently, we know 257 SySts in the Milky Way and 69 in external galaxies. Howeve… ▽ More Symbiotic stars (SySts) are long-period interacting binaries composed of a hot compact star, an evolved giant star, and a tangled network of gas and dust nebulae. They represent unique laboratories for studying a variety of important astrophysical problems, and have also been proposed as possible progenitors of SNIa. Presently, we know 257 SySts in the Milky Way and 69 in external galaxies. However, these numbers are still in striking contrast with the predicted population of SySts in our Galaxy. Because of other astrophysical sources that mimic SySt colors, no photometric diagnostic tool has so far demonstrated the power to unambiguously identify a SySt, thus making the recourse to costly spectroscopic follow-up still inescapable. In this paper we present the concept, commissioning, and science verification phases, as well as the first scientific results, of RAMSES II - a Gemini Observatory Instrument Upgrade Project that has provided each GMOS instrument at both Gemini telescopes with a set of narrow-band filters centered on the Raman OVI 6830 A band. Continuum-subtracted images using these new filters clearly revealed known SySts with a range of Raman OVI line strengths, even in crowded fields. RAMSES II observations also produced the first detection of Raman OVI emission from the SySt LMC 1 and confirmed Hen 3-1768 as a new SySt - the first photometric confirmation of a SySt. Via Raman OVI narrow-band imaging, RAMSES II provides the astronomical community with the first purely photometric tool for hunting SySts in the local Universe. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: 23 pages, 16 figures, 5 tables, accepted for publication in AJ

Adding a Suite of Chemical Abundances to the MACER Code for the Evolution of Massive Elliptical Galaxies

Authors: Zhaoming Gan, Ena Choi, Jeremiah P. Ostriker, Luca Ciotti, Silvia Pellegrini

Abstract: We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size… ▽ More We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years. △ Less

Submitted 8 April, 2019; v1 submitted 25 December, 2018; originally announced December 2018.

Comments: 10 pages, 10 figures; accepted by ApJ

The impact of AGN on stellar kinematics and orbits in simulated massive galaxies

Authors: Matteo Frigo, Thorsten Naab, Michaela Hirschmann, Ena Choi, Rachel Somerville, Davor Krajnovic, Romeel Davé, Michele Cappellari

Abstract: We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulatio… ▽ More We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $ξ_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/σ$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $ξ_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Synthetic nebular emission from massive galaxies II: ultraviolet-line diagnostics of dominant ionizing sources

Authors: Michaela Hirschmann, Stéphane Charlot, Anna Feltre, Thorsten Naab, Rachel S. Somerville, Ena Choi

Abstract: We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which o… ▽ More We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$λ$5007/H$β$ versus [NII]$λ$6584/H$α$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$λ$1663) versus OIII]$λ$1663/HeII$λ$1640, CIII]$λ$1908/HeII$λ$1640 versus OIII]$λ$1663/HeII$λ$1640, and CIV$λ$1550/CIII]$λ$1908 versus CIII]$λ$1908/CII$λ$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$λ$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 23 pages, 10 figures, submitted to MNRAS

Optical observations of NEA 3200 Phaethon (1983 TB) during the 2017 apparition

Authors: Myung-Jin Kim, Hee-Jae Lee, Sang-Min Lee, Dong-Heun Kim, Fumi Yoshida, Przemyslaw Bartczak, Grzegorz Dudzinski, Jintae Park, Young-Jun Choi, Hong-Kyu Moon, Hong-Suh Yim, Jin Choi, Eun-Jung Choi, Joh-Na Yoon, Alexander Serebryanskiy, Maxim Krugov, Inna Reva, Kamoliddin E. Ergashev, Otabek Burkhonov, Shuhrat A. Ehgamberdiev, Yunus Turayev, Zhong-Yi Lin, Tomoko Arai, Katsuhito Ohtsuka, Takashi Ito , et al. (2 additional authors not shown)

Abstract: The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets o… ▽ More The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented. △ Less

Submitted 16 September, 2018; originally announced September 2018.

Comments: 14 pages, 4 figures, 1 figure in Appendix A. Accepted for publication in Astronomy & Astrophysics (A&A)

Journal ref: A&A 619, A123 (2018)

Linking galaxy structural properties and star formation activity to black hole activity with IllustrisTNG

Authors: Melanie Habouzit, Shy Genel, Rachel S. Somerville, Dale Kocevski, Michaela Hirschmann, Avishai Dekel, Ena Choi, Dylan Nelson, Annalisa Pillepich, Paul Torrey, Lars Hernquist, Mark Vogelsberger, Rainer Weinberger, Volker Springel

Abstract: We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good… ▽ More We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L_x ~ 10^{43-44} erg/s), and a lower number of bright AGN (L_x>10^{44} erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10^{9} Msun BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyze how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z=0 M_{star}>10^{11} Msun quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z~1.5-2), and the lowest among compact quiescent galaxies (6-10% at z~1.5-2). △ Less

Submitted 8 January, 2019; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: 35 pages, 22 figures, accepted for publication in MNRAS

The Role of Black Hole Feedback on Size and Structural Evolution in Massive Galaxies

Authors: Ena Choi, Rachel S. Somerville, Jeremiah P. Ostriker, Thorsten Naab, Michaela Hirschmann

Abstract: We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nuc… ▽ More We use cosmological hydrodynamical simulations to investigate the role of feedback from accreting black holes on the evolution of sizes, compactness, stellar core density and specific star-formation of massive galaxies with stellar masses of $ M_{star} > 10^{10.9} M_{\odot}$. We perform two sets of cosmological zoom-in simulations of 30 halos to z=0: (1) without black holes and Active Galactic Nucleus (AGN) feedback and (2) with AGN feedback arising from winds and X-ray radiation. We find that AGN feedback can alter the stellar density distribution, reduce the core density within the central 1 kpc by 0.3 dex from z=1, and enhance the size growth of massive galaxies. We also find that galaxies simulated with AGN feedback evolve along similar tracks to those characterized by observations in specific star formation versus compactness. We confirm that AGN feedback plays an important role in transforming galaxies from blue compact galaxies into red extended galaxies in two ways: (1) it effectively quenches the star formation, transforming blue compact galaxies into compact quiescent galaxies and (2) it also removes and prevents new accretion of cold gas, shutting down in-situ star formation and causing subsequent mergers to be gas-poor or mixed. Gas poor minor mergers then build up an extended stellar envelope. AGN feedback also puffs up the central region through the fast AGN driven winds as well as the slow expulsion of gas while the black hole is quiescent. Without AGN feedback, large amounts of gas accumulate in the central region, triggering star formation and leading to overly massive blue galaxies with dense stellar cores. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: 13 pages, 7 figures, Accepted for publication in ApJ

IQ-Collaboratory 1.1: the Star-Forming Sequence of Simulated Central Galaxies

Authors: ChangHoon Hahn, Tjitske K. Starkenburg, Ena Choi, Romeel Davé, Claire M. Dickey, Marla C. Geha, Shy Genel, Christopher C. Hayward, Ariyeh H. Maller, Nityasri Mandyam, Viraj Pandya, Gergö Popping, Mika Rafieferantsoa, Rachel S. Somerville, Jeremy L. Tinker

Abstract: A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this metho… ▽ More A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Comments: 28 pages, 15 figures

Momentum-driven Winds from Radiatively Efficient Black Hole Accretion and Their Impact on Galaxies

Authors: Ryan Brennan, Ena Choi, Rachel S. Somerville, Michaela Hirschmann, Thorsten Naab, Jeremiah P. Ostriker

Abstract: We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating,… ▽ More We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN' model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN' model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective' feedback - the outflows themselves which drive gas out of galaxies - and `preventative' feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: 27 pages, 20 figures, Accepted for publication in ApJ

Synthetic nebular emission from massive galaxies I: origin of the cosmic evolution of optical emission-line ratios

Authors: Michaela Hirschmann, Stephane Charlot, Anna Feltre, Thorsten Naab, Ena Choi, Jeremiah P. Ostriker, Rachel S. Somerville

Abstract: Galaxies occupy different regions of the [OIII]$\lambda5007$/H$β$-versus-[NII]$\lambda6584$/H$α$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in… ▽ More Galaxies occupy different regions of the [OIII]$\lambda5007$/H$β$-versus-[NII]$\lambda6584$/H$α$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$β$, [NII]$\lambda6584$/H$α$, [SII]$λ\lambda6717,6731$/H$α$ and [OI]$\lambda6300$/H$α$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$β$ increases and [NII]/H$α$, [SII]/H$α$ and [OI]/H$α$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$β$ and [NII]/H$α$. For [SII]/H$α$ and [OI]/H$α$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios. △ Less

Submitted 31 August, 2017; v1 submitted 31 May, 2017; originally announced June 2017.

Comments: 31 pages, 14 figures, accepted for publication in MNRAS

Relativistic hydrodynamic jets in the intracluster medium

Authors: Eunwoo Choi

Abstract: We have performed the first three-dimensional relativistic hydrodynamic simulations of extragalactic jets of pure leptonic and baryonic plasma compositions propagating into a hydrostatic intracluster medium environment. The numerical simulations use a general equation of state for a multi-component relativistic gas which closely reproduces the Synge equation of state for a relativistic perfect gas… ▽ More We have performed the first three-dimensional relativistic hydrodynamic simulations of extragalactic jets of pure leptonic and baryonic plasma compositions propagating into a hydrostatic intracluster medium environment. The numerical simulations use a general equation of state for a multi-component relativistic gas which closely reproduces the Synge equation of state for a relativistic perfect gas. We find that morphological and dynamical differences between leptonic and baryonic jets are much less evident than those between hot and cold jets. In all these models, the jets first propagate with essentially constant velocities within the core radius of the intracluster medium and then accelerate progressively so as to increase the jet advance velocity by a factor of between 1.2 and 1.6 at the end of simulations, depending upon the models. The temporal evolution of the average cavity pressure is not consistent with that expected by the extended theoretical model even if the average cavity pressure decreases as a function of time with a power-law. Our simulations produce synthetic radio images that are dominated by bright hot spots and appear similar to observations of the extended radio galaxies with collimated radio jets. These bright radio lobes would be visible as dark regions in X-ray images and are morphologically similar to observed X-ray cavities in the intracluster medium. This supports the expectation that the bow shock surrounding the head of the jet is important mechanism for producing X-ray cavities in the intracluster medium. Although there are quantitative differences among the models the total radio and X-ray intensity curves show qualitatively similar trends in all of them. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS

AGN feedback, quiescence and CGM metal enrichment in early-type galaxies

Authors: Maximilian Eisenreich, Thorsten Naab, Ena Choi, Jeremiah P. Ostriker, Eric Emsellem

Abstract: We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the… ▽ More We present three-dimensional hydrodynamical simulations showing the effect of kinetic and radiative AGN feedback on a model galaxy representing a massive quiescent low-redshift early-type galaxy of $M_* = 8.41\times 10^{10} M_\odot$, harbouring a $M_\mathrm{BH} = 4\times 10^8 M_\odot $ black hole surrounded by a cooling gaseous halo. We show that, for a total baryon fraction of $\sim 20\%$ of the cosmological value, feedback from the AGN can keep the galaxy quiescent for about 4.35 Gyr and with properties consistent with black hole mass and X-ray luminosity scaling relations. However, this can only be achieved if the AGN feedback model includes both kinetic and radiative feedback modes. The simulation with only kinetic feedback fails to keep the model galaxy fully quiescent, while one with only radiative feedback leads to excessive black-hole growth. For higher baryon fractions (e.g. 50\% of the cosmological value), the X-ray luminosities exceed observed values by at least one order of magnitude, and rapid cooling results in a star-forming galaxy. The AGN plays a major role in keeping the circumgalactic gas at observed metallicities of $Z/Z_\odot \gtrsim 0.3 $ within the central $\sim 30$ kpc by venting nuclear gas enriched with metals from residual star formation activity. As indicated by previous cosmological simulations, our results are consistent with a model for which the black hole mass and the total baryon fraction are set at higher redshifts $z > 1$ and the AGN alone can keep the model galaxy on observed scaling relations. Models without AGN feedback violate both the quiescence criterion as well as CGM metallicity constraints. △ Less

Submitted 7 April, 2017; v1 submitted 22 February, 2017; originally announced February 2017.

Comments: 19 pages, 15 figures. Accepted for publication in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 468 (1): 751-768

Mitigating radiation damage of single photon detectors for space applications

Authors: Elena Anisimova, Brendon L. Higgins, Jean-Philippe Bourgoin, Miles Cranmer, Eric Choi, Danya Hudson, Louis P. Piche, Alan Scott, Vadim Makarov, Thomas Jennewein

Abstract: Single-photon detectors in space must retain useful performance characteristics despite being bombarded with sub-atomic particles. Mitigating the effects of this space radiation is vital to enabling new space applications which require high-fidelity single-photon detection. To this end, we conducted proton radiation tests of various models of avalanche photodiodes (APDs) and one model of photomult… ▽ More Single-photon detectors in space must retain useful performance characteristics despite being bombarded with sub-atomic particles. Mitigating the effects of this space radiation is vital to enabling new space applications which require high-fidelity single-photon detection. To this end, we conducted proton radiation tests of various models of avalanche photodiodes (APDs) and one model of photomultiplier tube potentially suitable for satellite-based quantum communications. The samples were irradiated with 106 MeV protons at doses approximately equivalent to lifetimes of 0.6 , 6, 12 and 24 months in a low-Earth polar orbit. Although most detection properties were preserved, including efficiency, timing jitter and afterpulsing probability, all APD samples demonstrated significant increases in dark count rate (DCR) due to radiation-induced damage, many orders of magnitude higher than the 200 counts per second (cps) required for ground-to-satellite quantum communications. We then successfully demonstrated the mitigation of this DCR degradation through the use of deep cooling, to as low as -86 degrees C. This achieved DCR below the required 200 cps over the 24 months orbit duration. DCR was further reduced by thermal annealing at temperatures of +50 to +100 degrees C. △ Less

Submitted 12 October, 2017; v1 submitted 3 February, 2017; originally announced February 2017.

Comments: The license has been corrected. Note that the license of v2 was incorrect and not valid. No other changes since v2

MSC Class: 81-05

Journal ref: EPJ Quantum Technol. 4, 10 (2017)

Modeling for Stellar Feedback in Galaxy Formation Simulations

Authors: Alejandro Núñez, Jeremiah P. Ostriker, Thorsten Naab, Ludwig Oser, Chia-Yu Hu, Ena Choi

Abstract: Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, e… ▽ More Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Strömgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses. △ Less

Submitted 4 January, 2017; originally announced January 2017.

Comments: Accepted for publication in the ApJ

The Nature of Massive Transition Galaxies in CANDELS, GAMA, and Cosmological Simulations

Authors: Viraj Pandya, Ryan Brennan, Rachel S. Somerville, Ena Choi, Guillermo Barro, Stijn Wuyts, Edward N. Taylor, Peter Behroozi, Allison Kirkpatrick, Sandra M. Faber, Joel Primack, David C. Koo, Daniel H. McIntosh, Dale Kocevski, Eric F. Bell, Avishai Dekel, Jerome J. Fang, Henry C. Ferguson, Norman Grogin, Anton M. Koekemoer, Yu Lu, Kameswara Mantha, Bahram Mobasher, Jeffrey Newman, Camilla Pacifici , et al. (3 additional authors not shown)

Abstract: We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom… ▽ More We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate Sérsic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$. △ Less

Submitted 14 July, 2017; v1 submitted 11 November, 2016; originally announced November 2016.

Comments: Re-submitted to MNRAS after referee revisions, the main change is that the paper has been significantly streamlined and is shorter (25 pages and 8 figures, without appendices)

Physics of Galactic Metals: Evolutionary Effects due to Production, Distribution, Feedback & Interaction with Black Holes

Authors: Ena Choi, Jeremiah P. Ostriker, Thorsten Naab, Rachel S. Somerville, Michaela Hirschmann, Alejandro Núñez, Chia-Yu Hu, Ludwig Oser

Abstract: We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, tu… ▽ More We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies. △ Less

Submitted 14 June, 2017; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: 16 pages, 11 figures, accepted for publication in ApJ

Radiation from Particles Accelerated in Relativistic Jet Shocks and Shear-flows

Authors: K. -I. Nishikawa, P. Hardee, I. Dutan, B. Zhang, A. Meli, E. J. Choi, K. Min, J. Niemiec, Y. Mizuno, M. Medvedev, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, D. Hartmann

Abstract: We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and… ▽ More We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants. △ Less

Submitted 22 December, 2014; originally announced December 2014.

Comments: 8 pages, 7 figures, 2014 Fermi Symposium proceedings - eConf C14102.1

A Study of the Early-stage Evolution of Relativistic Electron-Ion Shock using 3D PIC Simulations

Authors: E. J. Choi, K. Min, K. -I. Nishikawa, C. R. Choi

Abstract: We report the results of a 3D particle-in-cell (PIC) simulation carried out to study the early-stage evolution of the shock formed when an unmagnetized relativistic jet interacts with an ambient electron-ion plasma. Full-shock structures associated with the interaction are observed in the ambient frame. When open boundaries are employed in the direction of the jet; the forward shock is seen as a h… ▽ More We report the results of a 3D particle-in-cell (PIC) simulation carried out to study the early-stage evolution of the shock formed when an unmagnetized relativistic jet interacts with an ambient electron-ion plasma. Full-shock structures associated with the interaction are observed in the ambient frame. When open boundaries are employed in the direction of the jet; the forward shock is seen as a hybrid structure consisting of an electrostatic shock combined with a double layer, while the reverse shock is seen as a double layer. The ambient ions show two distinct features across the forward shock: a population penetrating into the shocked region from the precursor region and an accelerated population escaping from the shocked region into the precursor region. This behavior is a signature of a combination of an electrostatic shock and a double layer. Jet electrons are seen to be electrostatically trapped between the forward and reverse shock structures showing a ring-like distribution in a phase-space plot, while ambient electrons are thermalized and become essentially isotropic in the shocked region. The magnetic energy density grows to a few percent of the jet kinetic energy density at both the forward and the reverse shock transition layers in a rather short time scale. We see little disturbance of the jet ions over this time scale. △ Less

Submitted 16 July, 2014; originally announced July 2014.

Comments: 11 pages, 5 figures, in press, in Physics of Plasmas

Journal ref: Phys. Plasmas 21, 072905 (2014)

The impact of mechanical AGN feedback on the formation of massive early-type galaxies

Authors: Ena Choi, Jeremiah P. Ostriker, Thorsten Naab, Ludwig Oser, Benjamin P. Moster

Abstract: We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run t… ▽ More We employ cosmological hydrodynamical simulations to investigate the effects of AGN feedback on the formation of massive galaxies with present-day stellar masses of $M_{stel} = 8.8 \times 10^{10} - 6.0 \times 10^{11} M_{sun}$. Using smoothed particle hydrodynamics simulations with a pressure-entropy formulation that allows an improved treatment of contact discontinuities and fluid mixing, we run three sets of simulations of 20 halos with different AGN feedback models: (1) no feedback, (2) thermal feedback, and (3) mechanical and radiation feedback. We assume that seed black holes are present at early cosmic epochs at the centre of emerging dark matter halos and trace their mass growth via gas accretion and mergers with other black holes. Both feedback models successfully recover the observed M_BH - sigma relation and black hole-to-stellar mass ratio for simulated central early-type galaxies. The baryonic conversion efficiencies are reduced by a factor of two compared to models without any AGN feedback at all halo masses. However, massive galaxies simulated with thermal AGN feedback show a factor of ~10-100 higher X-ray luminosities than observed. The mechanical/radiation feedback model reproduces the observed correlation between X-ray luminosities and velocity dispersion, e.g. for galaxies with sigma = 200 km/s, the X-ray luminosity is reduced from $10^{42}$ erg/s to $10^{40}$ erg/s. It also efficiently suppresses late time star formation, reducing the specific star formation rate from $10^{-10.5}$ $yr^{-1}$ to $10^{-14}$ $yr^{-1}$ on average and resulting in quiescent galaxies since z=2, whereas the thermal feedback model shows higher late time in-situ star formation rates than observed. △ Less

Submitted 13 March, 2015; v1 submitted 5 March, 2014; originally announced March 2014.

Comments: 13 pages, 11 figures, accepted for the publication in MNRAS

Consequences of Mechanical and Radiative Feedback from Black Holes in Disc Galaxy Mergers

Authors: Ena Choi, Thorsten Naab, Jeremiah P. Ostriker, Peter H. Johansson, Benjamin P. Moster

Abstract: We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN fe… ▽ More We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $ε_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - σ$) as well as X-ray luminosity ($L_X - σ$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim σ\lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model. △ Less

Submitted 10 March, 2014; v1 submitted 16 August, 2013; originally announced August 2013.

Comments: 16 pages, 10 figures, resubmitted to MNRAS

Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability

Authors: K. -I. Nishikawa, P. Hardee, B. Zhang, I. Dutan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, D. H. Hartmann

Abstract: We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath config… ▽ More We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($γ_{\rm j} = 1.5$) is larger than for a relativistic jet case ($γ_{\rm j} = 15$). △ Less

Submitted 10 July, 2013; originally announced July 2013.

Comments: 6 pages, 6 figures, presented at Dynamical processes in space plasmas II, Isradinamic 2012, in press, ANGEO. arXiv admin note: text overlap with arXiv:1303.2569

Journal ref: Ann.Geophys. 31 (2013) 1535

Radiation from accelerated particles in relativistic jets with shocks, shear-flow, and reconnection

Authors: K. -I. Nishikawa, P. Hardee, B. Zhang, I. Dutan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, D. H. Hartmann

Abstract: We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find… ▽ More We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants. △ Less

Submitted 11 March, 2013; originally announced March 2013.

Comments: 6 pages, 4 figures, 2012 Fermi Symposium proceedings - eConf C121028

Radiative and Momentum Based Mechanical AGN Feedback in a 3-Dimensional Galaxy Evolution Code

Authors: Ena Choi, Jeremiah P. Ostriker, Thorsten Naab, Peter H. Johansson

Abstract: We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and… ▽ More We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem. △ Less

Submitted 3 July, 2012; v1 submitted 9 May, 2012; originally announced May 2012.

Comments: 15 pages, 9 figures, ApJ in press

The Effect of Anisotropic Conduction on the Thermal Instability in the Interstellar Medium

Authors: Ena Choi, James M. Stone

Abstract: Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate… ▽ More Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM. △ Less

Submitted 20 December, 2011; originally announced December 2011.

Comments: 10 pages, 12 figures, Accepted for the publication in ApJ

Simulation of Relativistic Jets and Associated Self-consistent Radiation

Authors: K. -I. Nishikawa, E. J. Choi, K. Min, P. Hardee, Y. Mizuno, B. Zhang, J. Niemiec, M. Medvedev, A. Nordlund, J. Fredriksen, M. Pohl, H. Sol, D. H. Hartmann, G. J. Fishman

Abstract: lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially dev… ▽ More lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. △ Less

Submitted 15 November, 2011; originally announced November 2011.

Comments: 4 pages, 2 figures, 2011 Fermi Symposium proceedings - eConf C110509

A Multidimensional Relativistic Hydrodynamics Code with a General Equation of State

Authors: Eunwoo Choi, Paul J. Wiita

Abstract: The ideal gas equation of state with a constant adiabatic index, although commonly used in relativistic hydrodynamics, is a poor approximation for most relativistic astrophysical flows. Here we propose a new general equation of state for a multi-component relativistic gas which is consistent with the Synge equation of state for a relativistic perfect gas and is suitable for numerical (special) rel… ▽ More The ideal gas equation of state with a constant adiabatic index, although commonly used in relativistic hydrodynamics, is a poor approximation for most relativistic astrophysical flows. Here we propose a new general equation of state for a multi-component relativistic gas which is consistent with the Synge equation of state for a relativistic perfect gas and is suitable for numerical (special) relativistic hydrodynamics. We also present a multidimensional relativistic hydrodynamics code incorporating the proposed general equation of state, based on the HLL scheme, which does not make use of a full characteristic decomposition of the relativistic hydrodynamic equations. The accuracy and robustness of this code is demonstrated in multidimensional calculations through several highly relativistic test problems taking into account nonvanishing tangential velocities. Results from three-dimensional simulations of relativistic jets show that the morphology and dynamics of the relativistic jets are significantly influenced by the different equation of state and by different compositions of relativistic perfect gases. Our new numerical code, combined with our proposed equation of state is very efficient and robust, and unlike previous codes, it gives very accurate results for thermodynamic variables in relativistic astrophysical flows. △ Less

Submitted 20 September, 2010; originally announced September 2010.

Comments: 32 pages, 9 figures, accepted by ApJS

Momentum Driving: which physical processes dominate AGN feedback?

Authors: Jeremiah P. Ostriker, Ena Choi, Luca Ciotti, Gregory S. Novak, Daniel Proga

Abstract: The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $η=dot{M_outflow}/dot{M_accretion}= (2 ε_w c^2)/v_w^2$ where (… ▽ More The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $η=dot{M_outflow}/dot{M_accretion}= (2 ε_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $η$, and many prior treatments have essentially assumed that $η=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $ε_w \lesssim 10^{-4}$. △ Less

Submitted 13 July, 2010; v1 submitted 16 April, 2010; originally announced April 2010.

Comments: 10 pages, 8 figures, resubmitted to ApJ, added references

Tracing the Filamentary Structure of the Galaxy Distribution at z~0.8

Authors: Ena Choi, Nicholas A. Bond, Michael A. Strauss, Alison L. Coil, Marc Davis, Christopher N. A. Willmer

Abstract: We study filamentary structure in the galaxy distribution at z ~ 0.8 using data from the Deep Extragalactic Evolutionary Probe 2 (DEEP2) Redshift Survey and its evolution to z ~ 0.1 using data from the Sloan Digital Sky Survey (SDSS). We trace individual filaments for both surveys using the Smoothed Hessian Major Axis Filament Finder, an algorithm which employs the Hessian matrix of the galaxy de… ▽ More We study filamentary structure in the galaxy distribution at z ~ 0.8 using data from the Deep Extragalactic Evolutionary Probe 2 (DEEP2) Redshift Survey and its evolution to z ~ 0.1 using data from the Sloan Digital Sky Survey (SDSS). We trace individual filaments for both surveys using the Smoothed Hessian Major Axis Filament Finder, an algorithm which employs the Hessian matrix of the galaxy density field to trace the filamentary structures in the distribution of galaxies. We extract 33 subsamples from the SDSS data with a geometry similar to that of DEEP2. We find that the filament length distribution has not significantly changed since z ~ 0.8, as predicted in a previous study using a $\Lamda$CDM cosmological N-body simulation. However, the filament width distribution, which is sensitive to the non-linear growth of structure, broadens and shifts to smaller widths for smoothing length scales of 5-10 Mpc/h from z ~ 0.8 to z ~ 0.1, in accord with N-body simulations. △ Less

Submitted 16 March, 2010; originally announced March 2010.

Comments: 10 pages, 8 figures, accepted for the publication in MNRAS

An Explicit Scheme for Incorporating Ambipolar Diffusion in a Magnetohydrodynamics Code

Authors: Eunwoo Choi, Jongsoo Kim, Paul J. Wiita

Abstract: We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms are included by explicit finite difference operators in a fully unsplit way, maintaining second order accuracy. The divergence-free condition of magnetic fiel… ▽ More We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms are included by explicit finite difference operators in a fully unsplit way, maintaining second order accuracy. The divergence-free condition of magnetic fields is exactly ensured at all times by a flux-interpolated constrained transport scheme. The super time stepping method is used to accelerate the timestep in high resolution calculations and/or in strong ambipolar diffusion. We perform two test problems, the steady-state oblique C-type shocks and the decay of Alfvén waves, confirming the accuracy and robustness of our numerical approach. Results from the simulations of the compressible MHD turbulence with ambipolar diffusion show the flexibility of our method as well as its ability to follow complex MHD flows in the presence of ambipolar diffusion. These simulations show that the dissipation rate of MHD turbulence is strongly affected by the strength of ambipolar diffusion. △ Less

Submitted 19 December, 2008; originally announced December 2008.

Comments: 25 pages, 5 figures, ApJS accepted

Testing the AGB Scenario as the Origin of the Extreme-Helium Population in omega Centauri

Authors: Ena Choi, Sukyoung K. Yi

Abstract: The most massive Galactic globular cluster, Omega centauri, appears to have multiple populations. Its bluest main sequence and extended horizontal branch stars are suggested to have the common origin, that is, an extremely high helium abundance of Y ~ 0.4. The high helium abundance is most often attributed to asymptotic giant branch (AGB) stars. In this study we test the AGB hypothesis. We simul… ▽ More The most massive Galactic globular cluster, Omega centauri, appears to have multiple populations. Its bluest main sequence and extended horizontal branch stars are suggested to have the common origin, that is, an extremely high helium abundance of Y ~ 0.4. The high helium abundance is most often attributed to asymptotic giant branch (AGB) stars. In this study we test the AGB hypothesis. We simulate the maximum-AGB models where the impact of AGB stars is maximised by assuming that supernova explosions do not affect the chemical evolution of the proto cloud. We compare the enrichment history of helium, metals, carbon and nitrogen to the observed values. Even under the most generous condition, the maximum-AGB models fail to reproduce the large values of helium Y ~ 0.4 and helium enrichment parameter Delta Y / Delta Z ~ 70 which were deduced from the colour-magnitude diagram fits. They also fail to reproduce the C and N contents of the blue population spectroscopically determined. We conclude that the AGB scenario with the canonical stellar evolution theory cannot explain the observational constraints and that the self chemical enrichment does not provide a viable solution. Alternative processes are desperately called for. △ Less

Submitted 9 April, 2008; originally announced April 2008.

Comments: 6 pages, 5 figures, Accepted by MNRAS

First stars as a possible origin for the helium-rich population in omega Centauri

Authors: Ena Choi, Sukyoung K. Yi

Abstract: The most massive Galactic globular cluster omega Cen appears to have two, or perhaps more, distinct main sequences. Its bluest main sequence is at the centre of debate because it has been suggested to have an extremely high helium abundance of Y ~ 0.4. The same helium abundance is claimed to explain the presence of extreme horizontal branch stars of omega Cen as well. This demands a relative hel… ▽ More The most massive Galactic globular cluster omega Cen appears to have two, or perhaps more, distinct main sequences. Its bluest main sequence is at the centre of debate because it has been suggested to have an extremely high helium abundance of Y ~ 0.4. The same helium abundance is claimed to explain the presence of extreme horizontal branch stars of omega Cen as well. This demands a relative helium to metal enrichment of deltaY/deltaZ ~ 70; that is, more than one order of magnitude larger than the generally accepted value. Candidate solutions, namely, AGB stars, massive stars, and supernovae, have been suggested; but in this study, we show that none of them is a viable channel, in terms of reproducing the high value of deltaY/deltaZ for the constrained age difference between the red and blue populations. Essentially no populations with an ordinary initial mass function, including those candidates, can produce such a high deltaY/deltaZ because they all produce metals as well as helium. As an alternative, we investigate the possibility of the stochastic ``first star'' contamination to the gas from which the younger generation of omega Cen formed. This requires the assumption that Population III star formation episode overlaps with that of Population II. While the required condition appears extreme, very massive objects in the first star generation provide a solution that is at least as plausible as any other suggestions made before. △ Less

Submitted 6 November, 2006; originally announced November 2006.

Comments: 6 pages, 5 figures, accepted for publication of MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.375:L1-L5,2007

Hydrodynamic Interactions of Relativistic Extragalactic Jets with Dense Clouds

Authors: Eunwoo Choi, Paul J. Wiita, Dongsu Ryu

Abstract: We have studied three-dimensional hydrodynamic interactions of relativistic extragalactic jets with two-phase ambient media. These jets propagate through a denser homogeneous gas and then impact clouds with densities 100 to 1000 times higher than the initial beam density. The deflection angle of the jet is influenced more by the density contrast of the cloud than by the beam Mach number of the j… ▽ More We have studied three-dimensional hydrodynamic interactions of relativistic extragalactic jets with two-phase ambient media. These jets propagate through a denser homogeneous gas and then impact clouds with densities 100 to 1000 times higher than the initial beam density. The deflection angle of the jet is influenced more by the density contrast of the cloud than by the beam Mach number of the jet. A relativistic jet with low relativistic beam Mach number can eventually be slightly bent after it crosses the dense cloud; however, we have not seen permanently bent structures in the interaction of a high relativistic beam Mach number jet with a cloud. The relativistic jet impacts on dense clouds do not necessarily destroy the clouds completely, and much of the cloud body can survive as a coherent blob. This enhancement of cloud durability is partly due to the geometric influence of the off-axis collisions we consider and also arises from the lower rate of cloud fragmentation through the Kelvin-Helmholtz instability for relativistic jets. To compare our simulations with observed extragalactic radio jets, we have computed the approximate surface distributions of synchrotron emission at different viewing angles. These surface intensity maps show relativistic jets interacting with clouds can produce synchrotron emission knots similar to structures observed in many VLBI-scale radio sources. We find that the synchrotron emission increases steeply at the moment of impact and the emission peaks right before the jet passes through the cloud. △ Less

Submitted 16 October, 2006; originally announced October 2006.

Comments: 32 pages, 7 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.655:769-780,2007

Equation of State in Numerical Relativistic Hydrodynamics

Authors: Dongsu Ryu, Indranil Chattopadhyay, Eunwoo Choi

Abstract: Relativistic temperature of gas raises the issue of the equation of state (EoS) in relativistic hydrodynamics. We study the EoS for numerical relativistic hydrodynamics, and propose a new EoS that is simple and yet approximates very closely the EoS of the single-component perfect gas in relativistic regime. We also discuss the calculation of primitive variables from conservative ones for the EoS… ▽ More Relativistic temperature of gas raises the issue of the equation of state (EoS) in relativistic hydrodynamics. We study the EoS for numerical relativistic hydrodynamics, and propose a new EoS that is simple and yet approximates very closely the EoS of the single-component perfect gas in relativistic regime. We also discuss the calculation of primitive variables from conservative ones for the EoS's considered in the paper, and present the eigenstructure of relativistic hydrodynamics for a general EoS, in a way that they can be used to build numerical codes. Tests with a code based on the Total Variation Diminishing (TVD) scheme are presented to highlight the differences induced by different EoS's. △ Less

Submitted 22 May, 2006; originally announced May 2006.

Comments: To appear in the ApJS September 2006, v166n1 issue. Pdf with full resolution figures can be downloaded from http://canopus.cnu.ac.kr/ryu/ryuetal.pdf

Numerical Relativistic Hydrodynamics Based on the Total Variation Diminishing Scheme

Authors: Eunwoo Choi, Dongsu Ryu

Abstract: This paper describes a multidimensional hydrodynamic code which can be used for studies of relativistic astrophysical flows. The code solves the special relativistic hydrodynamic equations as a hyperbolic system of conservation laws based on the total variation diminishing (TVD) scheme. It uses a new set of conserved quantities and employs an analytic formula for transformation from the conserve… ▽ More This paper describes a multidimensional hydrodynamic code which can be used for studies of relativistic astrophysical flows. The code solves the special relativistic hydrodynamic equations as a hyperbolic system of conservation laws based on the total variation diminishing (TVD) scheme. It uses a new set of conserved quantities and employs an analytic formula for transformation from the conserved quantities in the reference frame to the physical quantities in the local rest frame. Several standard tests, including relativistic shock tubes, a relativistic wall shock, and a relativistic blast wave, are presented to demonstrate that the code captures discontinuities correctly and sharply in ultrarelativistic regimes. The robustness and flexibility of the code are demonstrated through test simulations of the relativistic Hawley-Zabusky shock and a relativistic extragalactic jet. △ Less

Submitted 5 July, 2005; originally announced July 2005.

Comments: To appear in New Astronomy. Pdf with full resolution figures can be downloaded from http://canopus.cnu.ac.kr/ryu/choiryu.pdf

Journal ref: New Astron.11:116-129,2005

Mass and Angular Momentum Transfer in the Massive Algol Binary RY Persei

Authors: P. Barai, D. R. Gies, E. Choi, V. Das, R. Deo, W. Huang, K. Marshall, M. V. McSwain, C. Ogden, A. Osterman, R. L. Riddle, J. L. Seymour Jr, D. W. Wingert, A. B. Kaye, G. J. Peters

Abstract: We present an investigation of H-alpha emission line variations observed in the massive Algol binary, RY Per. We give new radial velocity data for the secondary based upon our optical spectra and for the primary based upon high dispersion UV spectra. We present revised orbital elements and an estimate of the primary's projected rotational velocity (which indicates that the primary is rotating 7… ▽ More We present an investigation of H-alpha emission line variations observed in the massive Algol binary, RY Per. We give new radial velocity data for the secondary based upon our optical spectra and for the primary based upon high dispersion UV spectra. We present revised orbital elements and an estimate of the primary's projected rotational velocity (which indicates that the primary is rotating 7 times faster than synchronous). We use a Doppler tomography algorithm to reconstruct the individual primary and secondary spectra in the region of H-alpha, and we subtract the latter from each of our observations to obtain profiles of the primary and its disk alone. Our H-alpha observations of RY Per show that the mass gaining primary is surrounded by a persistent but time variable accretion disk. The profile that is observed outside-of-eclipse has weak, double-peaked emission flanking a deep central absorption, and we find that these properties can be reproduced by a disk model that includes the absorption of photospheric light by the band of the disk seen in projection against the face of the star. We developed a new method to reconstruct the disk surface density distribution from the ensemble of H-alpha profiles observed around the orbit, and this method accounts for the effects of disk occultation by the stellar components, the obscuration of the primary by the disk, and flux contributions from optically thick disk elements. The resulting surface density distribution is elongated along the axis joining the stars, in the same way as seen in hydrodynamical simulations of gas flows that strike the mass gainer near trailing edge of the star. This type of gas stream configuration is optimal for the transfer of angular momentum, and we show that rapid rotation is found in other Algols that have passed through a similar stage. △ Less

Submitted 2 April, 2004; v1 submitted 26 September, 2003; originally announced September 2003.

Comments: 39 pages, 12 figures, ApJ in press, 2004 June 20 issue

Journal ref: Astrophys.J.608:989-1000,2004

Baryonic fraction in the cold plus hot dark matter universe

Authors: Eunwoo Choi, Dongsu Ryu

Abstract: We report a study to constrain the fraction of baryonic matter in the cold plus hot dark matter (CHDM) universe by numerical simulations which include the hydrodynamics of baryonic matter as well as the particle dynamics of dark matter. Spatially flat, COBE-normalized CHDM models with the fraction of hot component $Ω_h\leq0.2$ are considered. We show that the models with $h/n/Ω_h=0.5/0.9/0.1$ an… ▽ More We report a study to constrain the fraction of baryonic matter in the cold plus hot dark matter (CHDM) universe by numerical simulations which include the hydrodynamics of baryonic matter as well as the particle dynamics of dark matter. Spatially flat, COBE-normalized CHDM models with the fraction of hot component $Ω_h\leq0.2$ are considered. We show that the models with $h/n/Ω_h=0.5/0.9/0.1$ and $0.5/0.9/0.2$ give a linear power spectrum which agrees well with observations. Here, $h$ is the Hubble constant in unit of $100~km/s/Mpc$ and $n$ is the spectral index of the initial power spectrum. Then, for the models with $h/n/Ω_h=0.5/0.9/0.2$ and baryonic fraction $Ω_b=0.05$ and 0.1 we calculate the properties of X-ray clusters, such as luminosity function, temperature distribution function, luminosity-temperature relation, histogram of gas to total mass ratio, and change of average temperature with redshift $z$. Comparison with the observed data of X-ray clusters indicates that the model with $Ω_b=0.05$ is preferred. The COBE-normalized CHDM model with $Ω_b>0.1$ may be ruled out by the present work, since it produces too many X-ray bright clusters. △ Less

Submitted 8 October, 1997; originally announced October 1997.

Comments: To appear in MN, using mn.sty and epsf.sty, 10 pages including ps figures, additional 2 gif figures