Normalized to: Barai, P.
[1]
oai:arXiv.org:1807.04768 [pdf] - 1901733
Intermediate-Mass Black Hole Growth and Feedback in Dwarf Galaxies at
High Redshifts
Submitted: 2018-07-12, last modified: 2019-06-10
Intermediate-mass black holes (IMBHs: masses between $100 - 10^{6}
M_{\odot}$) historically comprise of an elusive population compared to
stellar-mass and supermassive BHs. Recently IMBHs have started to be observed
at the centers of low-mass galaxies. We perform cosmological hydrodynamical
simulations of $(2 h^{-1} ~ {\rm Mpc})^3$ comoving boxes and investigate the
growth and feedback of central IMBHs in dwarf galaxies (DGs). The earliest BHs
appear at $z \sim 18 - 25$, and grow thereafter by accreting gas and by merger
with other BHs. We find that, starting from $10^{2} M_{\odot}$, it is possible
to build up IMBHs of a few$\times 10^{5} - 10^{6} M_{\odot}$ by $z = 5$, when
the BHs are seeded in halos less massive than $4 \times 10^{7} M_{\odot}$. The
BH accretion rates increase with time, and reaches $\dot{M}_{\rm BH} = (0.2 -
0.8) \dot{M}_{\rm Edd}$ for the massive IMBHs by $z = 4$. The star formation
rate density (SFRD) evolution of the DGs (stellar mass $10^{5} - 10^{8}
M_{\odot}$) has a peak plateau between $z = 4 - 6$. Star formation is quenched
between $z = 9 - 4$. The SFRD is reduced by factors up to $3$, when the BHs
have grown to a few times $10^5 M_{\odot}$. Even in the presence of stronger
SN-driven mass ejection, the BHs continue to grow up to $z \sim 6$, sustained
by gas inflows driven by galaxy mergers and interactions in a cosmological
environment. Our conclusions, based on numerical simulation results, support
the scenario that early feedback from IMBHs in gas-rich DGs at $z = 5 - 8$ can
potentially solve several anomalies in the DG mass range within the concordance
$\Lambda$CDM cosmological scenario (Silk 2017). Our results suggest that IMBHs
at DG centers grow faster than their host galaxies in the early Universe, and
the resulting BH feedback turns the DGs and the BHs dormant.
[3]
oai:arXiv.org:1903.04632 [pdf] - 1849249
Supermassive and Intermediate-Mass Black Hole Growth at Galaxy Centers
and Resulting Feedback using Cosmological Simulations
Submitted: 2019-03-11
Accretion of matter onto central Black Holes (BHs) in galaxies liberates
enormous amounts of feedback energy, which influence the formation and
evolution of structures, affecting the environment from pc to Mpc scales. These
BHs are usually Supermassive BHs (SMBHs: mass $\geq 10^6 M_{\odot}$) existing
at the centers of active galactic nuclei (AGN), which are widely observed
through their multi-wavelength emission at all cosmic epochs. The SMBH energy
output is often observed as powerful AGN outflows in a wide variety of forms.
Relatively recently, Intermediate-Mass BHs (IMBHs: mass = $100 - 10^6
M_{\odot}$) have started to be observed hosted in Dwarf Galaxy (DG) centers.
Some of the central IMBHs in DGs show signatures of activity in the form of
low-luminosity AGN. We have performed Cosmological Hydrodynamical Simulations
to probe SMBHs in high-z quasars (Barai et al. 2018), and IMBHs in DGs (Barai &
de Gouveia Dal Pino 2019). Our simulations employ the 3D TreePM SPH code
GADGET-3, and include metal cooling, star formation, chemical enrichment,
stellar evolution, supernova feedback, AGN accretion and feedback. Analyzing
the simulation output in post-processing, we investigate the growth of the
first IMBHs, and the growth of the first SMBHs, their impact on star-formation,
as well as their co-evolution with the respective host galaxies. We quantify
the impact of SMBHs and IMBHs on their host galaxies, especially the effects on
quenching star-formation. We also study the corresponding BH outflow
properties.
[4]
oai:arXiv.org:1811.06576 [pdf] - 2119801
Intermediate-Mass Black Hole Feedback in Dwarf Galaxies: a View from
Cosmological Simulations
Submitted: 2018-11-15
Black holes are usually observed to be of stellar-mass or supermassive. By
natural extension, there should be a population of Intermediate-Mass Black
Holes (IMBHs: with mass between $100$ to $10^6 M_{\odot}$) in the Universe;
which has started to been observed. An exciting claim has been made recently by
Silk (2017): that early feedback by IMBHs in gas-rich dwarf galaxies at
$z=5-8$, can potentially solve multiple dwarf galaxy problems within the
$\Lambda$-cold-dark-matter cosmology. We are performing Cosmological
Hydrodynamical Simulations of $(2 Mpc)^3$ volumes, starting from $z=100$, to
test the case for IMBHs in Dwarf Galaxies. Black holes of mass $1000 M_{\odot}$
are seeded inside halos when they reach a mass of $10^7 M_{\odot}$. The black
holes grow by accretion of gas from their surroundings and by merger with other
black holes, and consequently eject feedback energy. We analyze the simulation
output in post-processing to study the growth of the first IMBHs, and their
impact on star-formation. Our conclusions, based on numerical simulation
results, support the phenomenological ideas made by Silk (2017). IMBHs at the
centers of dwarf galaxies can be a strong source of feedback to quench
star-formation and generate outflows. At the same time, these IMBHs form the
missing link between stellar-mass and supermassive BHs.
[5]
oai:arXiv.org:1811.06025 [pdf] - 1791385
Large-Scale Diffuse Intergalactic Magnetic Fields Constraints with the
Cherenkov Telescope Array
Submitted: 2018-11-14
Magnetic fields of the order of $\mu$-Gauss are observationally detected in
galaxies and galaxy clusters, which can be (at least) in part originated by the
amplification of much weaker primordial seed fields. These fields should be
carried out by strong galactic outflows, magnetically enriching the
InterGalactic Medium (IGM). However direct observation of magnetic fields in
the IGM is scarce. This talk will give a review of how Intergalactic Magnetic
Field (IGMF) can be constrained using gamma-ray observations. High-energy TeV
photons emitted by distant blazars can interact with the cosmic extragalactic
optical/infrared/microwave background light, producing electron-positron pairs,
and initiating electromagnetic cascades in the IGM. The charged component of
these cascades is deflected by IGMFs, thereby reducing the observed point-like
TeV flux, and creating an extended image in the GeV energy range, which can
potentially be detected with $\gamma$-ray telescopes (Fermi-LAT, HESS, CTA).
Studies (e.g., Neronov & Vovk 2010, Dolag et al. 2011) have put lower limits on
the IGMF strength of the order of $10^{-16} - 10^{-15} G$, and filling factors
of $60\%$. This talk will describe the constraints which the Cherenkov
Telescope Array sensitivity is expected to give (CTA Consortium 2018).
[6]
oai:arXiv.org:1804.05566 [pdf] - 1779514
The dense molecular gas in the $\rm z\sim6$ QSO SDSS J231038.88+185519.7
resolved by ALMA
Feruglio, C.;
Fiore, F.;
Carniani, S.;
Maiolino, R.;
D'Odorico, V.;
Luminari, A.;
Barai, P.;
Bischetti, M.;
Bongiorno, A.;
Cristiani, S.;
Ferrara, A.;
Gallerani, S.;
Marconi, A.;
Pallottini, A.;
Piconcelli, E.;
Zappacosta, L.
Submitted: 2018-04-16, last modified: 2018-08-20
We present ALMA observations of the CO(6-5) and [CII] emission lines and the
sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS
J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic
beam that is ten times smaller in angular size, we have achieved ten times
better sensitivity in the CO(6-5) line, and two and half times better
sensitivity in the [CII] line, enabling us to resolve the molecular gas
emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and
of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is
thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO}
= 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of
$\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the
observed CO velocity gradient is due to an inclined rotating disk, we derive a
dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~
M_{\odot}$, which is a factor of approximately two smaller than the previously
reported estimate based on [CII]. Regarding the central black hole, we provide
a new estimate of the black hole mass based on the C~IV emission line detected
in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~
M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$,
where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio
$v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows
and/or complex kinematics due to a merger event. We estimate a global Toomre
parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare,
at the same angular resolution, the CO(6-5) and [CII] distributions, finding
that dense molecular gas is more centrally concentrated with respect to [CII].
We find that the current BH growth rate is similar to that of its host galaxy.
[7]
oai:arXiv.org:1706.06123 [pdf] - 1708260
Absorption systems at $z\sim 2$ as a probe of the circum galactic
medium: a probabilistic approach
Submitted: 2017-06-19, last modified: 2018-07-02
We characterize the properties of the intergalactic medium (IGM) around a
sample of galaxies extracted from state-of-the-art hydrodynamical simulations
of structure formation in a cosmological volume of 25 Mpc comoving at $z\sim
2$. The simulations are based on two different subresolution schemes for star
formation and supernova feedback: the MUlti-Phase Particle Integrator (MUPPI)
scheme and the Effective Model. We develop a quantitative and probabilistic
analysis based on the apparent optical depth method of the properties of the
absorbers as a function of impact parameter from their nearby galaxies: in such
a way we probe different environments from circumgalactic medium to low-density
filaments. Absorbers' properties are then compared with a spectroscopic
observational data set obtained from high-resolution quasar spectra. Our main
focus is on the N$_{\rm{ CIV}}$-N$_{\rm { HI}}$ relation around simulated
galaxies: the results obtained with MUPPI and the Effective model are
remarkably similar, with small differences only confined to regions at impact
parameters $b = [1-3] \times r_{\rm {vir}}$. Using $\mbox{C IV}$ as a tracer of
the metallicity, we obtain evidence that the observed metal absorption systems
have the highest probability to be confined in a region of 150-400 kpc around
galaxies. Near-filament environments have instead metallicities too low to be
probed by present-day telescopes, but could be probed by future spectroscopical
studies. Finally we compute $\mbox{C IV}$ covering fractions which are in
agreement with observational data.
[9]
oai:arXiv.org:1707.03014 [pdf] - 1598088
Quasar outflows at $z \geq 6$: the impact on the host galaxies
Submitted: 2017-07-10
We investigate quasar outflows at $z \geq 6$ by performing zoom-in
cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we
zoom in the $2 R_{200}$ region around a $2 \times 10^{12} M_{\odot}$ halo at $z
= 6$, inside a $(500 ~ {\rm Mpc})^3$ comoving volume. We compare the results of
our AGN runs with a control simulation in which only stellar/SN feedback is
considered. Seeding $10^5 M_{\odot}$ BHs at the centers of $10^{9} M_{\odot}$
halos, we find the following results. BHs accrete gas at the Eddington rate
over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{\rm BH} = 4
\times 10^9 M_{\odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($\dot{M}_{\rm
out} \sim 2000 M_{\odot}$/yr) outflows of shock-heated low-density gas form at
$z \sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over
$z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is
reduced by a factor of $5$ ($1000$). We analyse the relative contribution of
multiple physical process: (i) disrupting cosmic filamentary cold gas inflows,
(ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and
find that AGN feedback has the following effects at $z = 6$. The inflowing gas
mass fraction is reduced by $\sim 12 \%$, the high-density gas fraction is
lowered by $\sim 13 \%$, and $\sim 20 \%$ of the gas outflows at a speed larger
than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at
$z \geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN
feedback quenches their star formation dominantly by powerful outflows ejecting
gas out of the host galaxy halo.
[10]
oai:arXiv.org:1611.03805 [pdf] - 1532574
Exploring the Thermal State of the Low-Density Intergalactic Medium at
z=3 with an Ultra-High Signal-to-Noise QSO Spectrum
Rorai, A.;
Becker, G. D.;
Haehnelt, M. G.;
Carswell, R. F.;
Bolton, J. S.;
Cristiani, S.;
D'Odorico, V.;
Cupani, G.;
Barai, P.;
Calura, F.;
Kim, T. -S.;
Pomante, E.;
Tescari, E.;
Viel, M.
Submitted: 2016-11-11
At low densities the standard ionisation history of the intergalactic medium
(IGM) predicts a decreasing temperature of the IGM with decreasing density once
hydrogen (and helium) reionisation is complete. Heating the high-redshift,
low-density IGM above the temperature expected from photo-heating is difficult,
and previous claims of high/rising temperatures in low density regions of the
Universe based on the probability density function (PDF) of the opacity in
Lyman-$\alpha$ forest data at $2<z<4$ have been met with considerable
scepticism, particularly since they appear to be in tension with other
constraints on the temperature-density relation (TDR). We utilize here an
ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel
technique to study the low opacity part of the PDF. We show that there is
indeed evidence (at 90% confidence level) that a significant volume fraction of
the under-dense regions at $z \sim 3$ has temperatures as high or higher than
those at densities comparable to the mean and above. We further demonstrate
that this conclusion is nevertheless consistent with measurements of a slope of
the TDR in over-dense regions that imply a decreasing temperature with
decreasing density, as expected if photo-heating of ionised hydrogen is the
dominant heating process. We briefly discuss implications of our findings for
the need to invoke either spatial temperature fluctuations, as expected during
helium reionization, or additional processes that heat a significant volume
fraction of the low-density IGM.
[11]
oai:arXiv.org:1608.06116 [pdf] - 1525209
Metals in the z~3 intergalactic medium: results from an ultra-high
signal-to-noise ratio UVES quasar spectrum
D'Odorico, V.;
Cristiani, S.;
Pomante, E.;
Carswell, R. F.;
Viel, M.;
Barai, P.;
Becker, G. D.;
Calura, F.;
Cupani, G.;
Fontanot, F.;
Haehnelt, M. G.;
Kim, T-S.;
Miralda-Escude, J.;
Rorai, A.;
Tescari, E.;
Vanzella, E.
Submitted: 2016-08-22, last modified: 2016-09-13
In this work, we investigate the abundance and distribution of metals in the
intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the
analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar
HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$
to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60
per cent of the considered redshift range down to $\simeq11.1$. In our sample,
all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption.
In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an
associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates
drop to $<10$ per cent, possibly due to our sensitivity limits and not to an
actual variation of the gas abundance properties. In the range $\log N_{\rm HI}
\ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2
larger than the fraction of HI lines lying in the circum-galactic medium (CGM)
of relatively bright Lyman-break galaxies hosted by dark matter haloes with
$\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with
the output of a grid of photoionization models and of two cosmological
simulations implies that the volume filling factor of the IGM gas enriched to a
metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$
percent. In conclusion, our results favour a scenario in which metals are found
also outside the CGM of bright star-forming galaxies, possibly due to pollution
by lower mass objects and/or to an early enrichment by the first sources.
[12]
oai:arXiv.org:1605.08051 [pdf] - 1444016
Kinetic AGN Feedback Effects on Cluster Cool Cores Simulated using SPH
Submitted: 2016-05-25, last modified: 2016-06-09
We implement novel numerical models of AGN feedback in the SPH code GADGET-3,
where the energy from a supermassive black hole (BH) is coupled to the
surrounding gas in the kinetic form. Gas particles lying inside a bi-conical
volume around the BH are imparted a one-time velocity (10,000 km/s) increment.
We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h
M_sun), which is initially evolved to form a dense cool core, having central
T<10^6 K. A BH resides at the cluster center, and ejects energy. The
feedback-driven fast wind undergoes shock with the slower-moving gas, which
causes the imparted kinetic energy to be thermalized. Bipolar bubble-like
outflows form propagating radially outward to a distance of a few 100 kpc. The
radial profiles of median gas properties are influenced by BH feedback in the
inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the
feedback efficiency, depletes the inner cool gas and reduces the hot gas
content, such that the initial cool core of the cluster is heated up within a
time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and
the central entropy flattens. Our implementation of BH thermal feedback (using
the same efficiency as kinetic), within the star-formation model, cannot do
this heating, where the cool core remains. The inclusion of cold gas accretion
in the simulations produces naturally a duty cycle of the AGN with a
periodicity of 100 Myr.
[13]
oai:arXiv.org:1601.04660 [pdf] - 1359231
Triple-ionised carbon associated with the low-density neutral hydrogen
gas at 1.7 < z < 3.3: the integrated N(HI)-N(CIV) relation
Submitted: 2016-01-18
From the Voigt profile fitting analysis of 183 intervening CIV systems at 1.7
< z < 3.3 in 23 high-quality UVES/VLT and HIRES/Keck QSO spectra, we find that
a majority of CIV systems (~75%) display a well-characterised scaling relation
between integrated column densities of HI and CIV with a negligible redshift
evolution, when column densities of all the HI and CIV components are
integrated within a given (-150, +150) km/sec range centred at the CIV flux
minimum. The integrated CIV column density N(CIV, sys) increases with N(HI,
sys) at log N(HI, sys) = 14.0--15.5 and log N(CIV, sys) = 11.8--14.0, then
becomes almost independent of N(HI, sys) at log N(HI, sys) > 16, with a large
scatter: at log N(HI, sys) = 14--22, log N(CIV, sys) = C1 / (log(NHI, sys) +
C2) + C3, with C1 = -1.90+0.55, C2 = -14.11+0.19 and C3 = 14.76+0.17,
respectively. The steep (flat) part is dominated by SiIV-free (SiIV-enriched)
CIV systems. Extrapolating the N(HI, sys)-N(CIV, sys) relation implies that
most absorbers with log N(HI) < 14 are virtually CIV-free. The N(HI,
sys)-N(CIV, sys) relation does not hold for individual components, clumps or
the integration velocity range less than +-100 km/sec. It is expected if the
line-of-sight extent of CIV is smaller than HI and N(CIV, sys) decreases more
rapidly than N(HI, sys) at the larger impact parameter, regardless of the
location of the HI+CIV gas in the IGM filaments or in the intervening galactic
halos.
[14]
oai:arXiv.org:1211.6952 [pdf] - 922872
Black Hole spin dependence of general relativistic multi-transonic
accretion close to the horizon
Das, Tapas K.;
Nag, Sankhasubhra;
Hegde, Swathi;
Bhattacharya, Sourav;
Maity, Ishita;
Czerny, Bozena;
Barai, Paramita;
Wiita, Paul J.;
Karas, Vladimir;
Naskar, Tapan
Submitted: 2012-11-29, last modified: 2014-12-31
We introduce a novel formalism to investigate the role of the spin angular
momentum of astrophysical black holes in influencing the behaviour of low
angular momentum general relativistic accretion. We propose a metric
independent analysis of axisymmetric general relativistic flow, and
consequently formulate the space and time dependent equations describing the
general relativistic hydrodynamic accretion flow in the Kerr metric. The
associated stationary critical solutions for such flow equations are provided
and the stability of the stationary transonic configuration is examined using
an elegant linear perturbation technique. We examine the properties of
infalling material for both prograde and retrograde accretion as a function of
the Kerr parameter at extremely close proximity to the event horizon. Our
formalism can be used to identify a new spectral signature of black hole spin,
and has the potential of performing the black hole shadow imaging corresponding
to the low angular momentum accretion flow.
[15]
oai:arXiv.org:1411.1409 [pdf] - 1223025
Galactic outflow and diffuse gas properties at z>=1 using different
baryonic feedback models
Submitted: 2014-11-05
We measure and quantify properties of galactic outflows and diffuse gas at $z
\geq 1$ in cosmological hydrodynamical simulations. Our novel sub-resolution
model, MUPPI, implements supernova feedback using fully local gas properties,
where the wind velocity and mass loading are not given as input. We find the
following trends at $z = 2$ by analysing central galaxies having a stellar mass
higher than $10^{9} M_{\odot}$. The outflow velocity and mass outflow rate
($\dot{M}_{\rm out}$) exhibit positive correlations with galaxy mass and with
the star formation rate (SFR). However, most of the relations present a large
scatter. The outflow mass loading factor ($\eta$) is between $0.2 - 10$. The
comparison Effective model generates a constant outflow velocity, and a
negative correlation of $\eta$ with halo mass. The number fraction of galaxies
where outflow is detected decreases at lower redshifts, but remains more than
$80 \%$ over $z = 1 - 5$. High SF activity at $z \sim 2 - 4$ drives strong
outflows, causing the positive and steep correlations of velocity and
$\dot{M}_{\rm out}$ with SFR. The outflow velocity correlation with SFR becomes
flatter at $z = 1$, and $\eta$ displays a negative correlation with halo mass
in massive galaxies. Our study demonstrates that both the MUPPI and Effective
models produce significant outflows at $\sim 1 / 10$ of the virial radius; at
the same time shows that the properties of outflows generated can be different
from the input speed and mass loading in the Effective model. Our MUPPI model,
using local properties of gas in the sub-resolution recipe, is able to develop
galactic outflows whose properties correlate with global galaxy properties, and
consistent with observations.
[16]
oai:arXiv.org:1403.3063 [pdf] - 1208345
Major Cluster Mergers and the Location of the Brightest Cluster Galaxy
Submitted: 2014-03-12
Using a large N-body cosmological simulation combined with a subgrid
treatment of galaxy formation, we study the formation and evolution of the
galaxy and cluster population in a comoving volume (100 Mpc)^3 in a LCDM
universe. At z = 0, our computational volume contains 1788 clusters with mass
M_cl > 1.1x10^12 Msun, including 18 massive clusters with M_cl > 10^14 Msun. It
also contains 1 088 797 galaxies with mass M_gal > 2x10^9 Msun and luminosity L
> 9.5x10^5 Lsun. For each cluster, we identified the brightest cluster galaxy
(BCG). We then computed the fraction f_BNC of clusters in which the BCG is not
the closest galaxy to the center of the cluster in projection, and the ratio
Dv/s, where Dv is the difference in radial velocity between the BCG and the
whole cluster, and s is the radial velocity dispersion of the cluster. f_BNC
increases from 0.05 for low-mass clusters (M_cl ~ 10^12 Msun) to 0.5 for
high-mass ones (M_cl > 10^14 Msun), with no dependence on cluster redshift. The
values of Dv/s vary from 0 to 1.8. These results are consistent with previous
observational studies, and indicate that the central galaxy paradigm, which
states that the BCG should be at rest at the center of the cluster, is usually
valid, but exceptions are too common to be ignored. Analysis of the merger
trees for the 18 most massive clusters in the simulation reveals that 16 of
these clusters have experienced major mergers in the past. These mergers leave
each cluster in a non-equilibrium state, but eventually the cluster settles
into an equilibrium configuration, unless it is disturbed by another major
merger. We found evidence that these mergers are responsible for the off-center
positions and peculiar velocities of some BCGs. Our results thus support the
merging-group scenario, in which some clusters form by the merger of smaller
groups in which the galaxies have already formed.
[17]
oai:arXiv.org:1312.5310 [pdf] - 1202354
Simulated star formation rate functions at z~4-7, and the role of
feedback in high-z galaxies
Submitted: 2013-12-17, last modified: 2014-01-06
We study the role of feedback from supernovae and black holes in the
evolution of the star formation rate function (SFRF) of z~4-7 galaxies. We use
a new set of cosmological hydrodynamic simulations, ANGUS (AustraliaN GADGET-3
early Universe Simulations), run with a modified and improved version of the
parallel TreePM-smoothed particle hydrodynamics code GADGET-3 called
P-GADGET3(XXL), that includes a self-consistent implementation of stellar
evolution and metal enrichment. In our simulations both Supernova (SN) driven
galactic winds and Active Galactic Nuclei (AGN) act simultaneously in a complex
interplay. The SFRF is insensitive to feedback prescription at z>5, meaning
that it cannot be used to discriminate between feedback models during
reionisation. However, the SFRF is sensitive to the details of feedback
prescription at lower redshift. By exploring different SN driven wind
velocities and regimes for the AGN feedback, we find that the key factor for
reproducing the observed SFRFs is a combination of "strong" SN winds and early
AGN feedback in low mass galaxies. Conversely, we show that the choice of
initial mass function and inclusion of metal cooling have less impact on the
evolution of the SFRF. When variable winds are considered, we find that a
non-aggressive wind scaling is needed to reproduce the SFRFs at z>4. Otherwise,
the amount of objects with low SFRs is greatly suppressed and at the same time
winds are not effective enough in the most massive systems.
[18]
oai:arXiv.org:1307.5326 [pdf] - 1172861
Kinetic or thermal AGN feedback in simulations of isolated and merging
disk galaxies calibrated by the M-sigma relation
Submitted: 2013-07-19, last modified: 2013-10-14
(Abridged) We investigate two modes of coupling the feedback energy from a
central AGN to the neighboring gas in galaxy simulations: kinetic - velocity
boost, and thermal - heating. We formulate kinetic feedback models for
energy-driven wind (EDW) and momentum-driven wind (MDW), using two free
parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel
numerical algorithm is implemented in the SPH code GADGET-3, to prevent the
expansion of a hole in the gas distribution around the BH. We perform
simulations of isolated evolution and merger of disk galaxies, of Milky-Way
mass as well as lower and higher masses. We find that in the isolated galaxy BH
kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer
that current prescriptions for BH subgrid physics in galaxy simulations can
grow the BH to observed values even in an isolated disk galaxy. The BH growth
is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation
obtained in our simulations with observational data, we conclude that it is
possible to find parameter sets for a fit in all the models, except for the
case with MDW feedback in a galaxy merger, in which the BH is always too
massive. The BH thermal feedback implementation of Springel, Di Matteo &
Hernquist (2005) within the multiphase star-formation model is found to have
negligible impact on gas properties; and the effect claimed in all previous
studies is attributed to gas depletion around the BH by the creation of an
artificial hole. The BH mass accretion rate in our simulations exhibit heavy
fluctuations. The star formation rate is quenched with feedback by removal of
gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give
the best metallicity observational diagnostic to distinguish between BH models.
[19]
oai:arXiv.org:1210.3582 [pdf] - 1152085
Galactic Winds in Cosmological Simulations of the Circumgalactic Medium
Barai, Paramita;
Viel, Matteo;
Borgani, Stefano;
Tescari, Edoardo;
Tornatore, Luca;
Dolag, Klaus;
Killedar, Madhura;
Monaco, Pierluigi;
D'Odorico, Valentina;
Cristiani, Stefano
Submitted: 2012-10-12, last modified: 2013-01-21
(Abridged) We explore new observationally-constrained sub-resolution models
of galactic outflows and investigate their impact on the circumgalactic medium
(CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic
simulations, including star formation, chemical enrichment, and four cases of
SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind
(CW), a radially varying wind (RVWa) where the outflow velocity has a positive
correlation with galactocentric distance (r), and a RVW with additional
dependence on halo mass (RVWb). Overall, we find that the outflows expel
metal-enriched gas away from galaxies, significantly quench star formation, and
enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy
centers are most sensitive to the choice of the wind model for halo masses
(10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case,
except that it substantially enriches the CGM: the carbon metallicity (Z_C) is
10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 -
10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find
that the impact of models CW and RVWa are similar, with the following
differences. RVWa causes a greater suppression of star formation rate at z < 5,
and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K)
gas than in CW. Outflows in CW produce a higher and earlier enrichment of some
IGM phases than RVWa. By visual inspection, we note that the RVWa model shows
more pronounced bipolar outflows and galactic disks. We present fitting
formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a
function of r. We predict observational diagnostics to distinguish between
different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h
comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.
[20]
oai:arXiv.org:1207.6673 [pdf] - 1125157
The Fate of Dwarf Galaxies in Clusters and the Origin of Intracluster
Stars. II. Cosmological Simulations
Submitted: 2012-07-27
We combine a N-body simulation algorithm with a subgrid treatment of galaxy
formation, mergers, and tidal destruction, and an observed conditional
luminosity function Phi(L|M), to study the origin and evolution of galactic and
extragalactic light inside a cosmological volume of size (100 Mpc)^3, in a
concordance LCDM model. This algorithm simulates the growth of large-scale
structures and the formation of clusters, the evolution of the galaxy
population in clusters, the destruction of galaxies by mergers and tides, and
the evolution of the intracluster light. We find that destruction of galaxies
by mergers dominates over destruction by tides by about an order of magnitude
at all redshifts. However, tidal destruction is sufficient to produce
intracluster light fractions f_ICL that are sufficiently high to match
observations. The bulk of the intracluster light (60%) is provided by
intermediate galaxies of total masses 10^11 Msun-10^12 Msun and stellar masses
6x10^8 Msun-3x10^10 Msun that were tidally destroyed by even more massive
galaxies. The contribution of low-mass galaxies to the intracluster light is
small and the contribution of dwarf galaxies is negligible, even though, by
numbers, most galaxies that are tidally destroyed are dwarfs. Tracking clusters
back in time, we find that their values of f_ICL tend to increase over time,
but can experience sudden changes that are sometimes non-monotonic. These
changes occur during major mergers involving clusters of comparable masses but
very different intracluster luminosities. Most of the tidal destruction events
take place in the central regions of clusters. As a result, the intracluster
light is more centrally concentrated than the galactic light. Our results
support tidal destruction of intermediate-mass galaxies as a plausible scenario
for the origin of the intracluster light.
[21]
oai:arXiv.org:1112.5483 [pdf] - 562518
Multiphase, non-spherical gas accretion onto a black hole
Submitted: 2011-12-22, last modified: 2012-05-16
(Abridged) We investigate non-spherical behavior of gas accreting onto a
central supermassive black hole performing simulations using the SPH code
GADGET-3 including radiative cooling and heating by the central X-ray source.
As found in earlier 1D studies, our 3D simulations show that the accretion mode
depends on the X-ray luminosity (L_X) for a fixed density at infinity and
accretion efficiency. In the low L_X limit, gas accretes in a stable,
spherically symmetric fashion. In the high L_X limit, the inner gas is
significantly heated up and expands, reducing the central mass inflow rate. The
expanding gas can turn into a strong enough outflow capable of expelling most
of the gas at larger radii. For some intermediate L_X, the accretion flow
becomes unstable developing prominent non-spherical features, the key reason
for which is thermal instability (TI) as shown by our analyses. Small
perturbations of the initially spherically symmetric accretion flow that is
heated by the intermediate L_X quickly grow to form cold and dense clumps
surrounded by overheated low density regions. The cold clumps continue their
inward motion forming filamentary structures; while the hot infalling gas slows
down because of buoyancy and can even start outflowing through the channels in
between the filaments. We found that the ratio between the mass inflow rates of
the cold and hot gas is a dynamical quantity depending on several factors:
time, spatial location, and L_X; and ranges between 0 and 4. We briefly discuss
astrophysical implications of such TI-driven fragmentation of accreting gas on
the formation of clouds in narrow and broad line regions of AGN, the formation
of stars, and the observed variability of the AGN luminiosity.
[22]
oai:arXiv.org:1112.3525 [pdf] - 462749
Gas Accretion onto a Supermassive Black Hole: a step to model AGN
feedback
Submitted: 2011-12-15, last modified: 2012-01-11
We study the gas accretion onto a supermassive black hole (SMBH) using the 3D
SPH code GADGET-3 on scales of 0.1-200 pc. First we test our code with
spherically symmetric, adiabatic Bondi accretion problem. We find that our
simulation can reproduce the expected Bondi accretion flow very well for a
limited amount of time until the effect of outer boundary starts to be visible.
We also find artificial heating of gas near the inner accretion boundary due to
the artificial viscosity of SPH. Second, we implement radiative cooling and
heating due to X-rays, and examine the impact of thermal feedback by the
central X-ray source. The accretion flow roughly follows the Bondi solution for
low central X-ray luminosities, however, the flow starts to exhibit
non-spherical fragmentation due to thermal instability for a certain range of
central L_X, and a strong overall outflow develops for greater L_X. The cold
gas develops filamentary structures that fall into the central SMBH, whereas
the hot gas tries to escape through the channels in-between the cold filaments.
Such fragmentation of accreting gas can assist in the formation of clouds
around AGN, induce star-formation, and contribute to the observed variability
of narrow-line regions.
[23]
oai:arXiv.org:1102.3925 [pdf] - 463142
SPH Simulations of Black Hole Accretion: A Step to Model Black Hole
Feedback in Galaxies
Submitted: 2011-02-18, last modified: 2011-07-27
(Abridged) We test how accurately the smoothed particle hydrodynamics (SPH)
numerical technique can follow spherically-symmetric Bondi accretion. Using the
3D SPH code GADGET-3, we perform simulations of gas accretion onto a central
supermassive black hole (SMBH) of mass 10^8 M_sun within the radial range of
0.1 - 200 pc. We carry out simulations without and with radiative heating by a
central X-ray corona and radiative cooling. For an adiabatic case, the radial
profiles of hydrodynamical properties match the Bondi solution, except near the
inner and outer radius of the computational domain. We find that adiabatic
Bondi accretion can be reproduced for durations of a few dynamical times at the
Bondi radius, and for longer times if the outer radius is increased. With
radiative heating and cooling included, the spherically accreting gas takes a
longer time to reach a steady-state than the adiabatic Bondi accretion runs,
and in some cases does not reach a steady-state even within several hundred
dynamical times. We find that artificial viscosity in the GADGET code causes
excessive heating near the inner radius, making the thermal properties of the
gas inconsistent with a physical solution. This overheating occurs typically
only in the supersonic part of the flow, so that it does not affect the mass
accretion rate. We see that increasing the X-ray luminosity produces a lower
central mass inflow rate, implying that feedback due to radiative heating is
operational in our simulations. With a sufficiently high X-ray luminosity, the
inflowing gas is radiatively heated up, and an outflow develops. We conclude
that the SPH simulations can capture the gas dynamics needed to study radiative
feedback provided artificial viscosity alters only highly supersonic part of
the inflow.
[24]
oai:arXiv.org:1011.2509 [pdf] - 1041887
Anisotropic AGN Outflows and Enrichment of the Intergalactic Medium. II.
Metallicity
Submitted: 2010-11-10
We investigate the large-scale influence of outflows from AGNs in enriching
the IGM with metals in a cosmological context. We combine cosmological
simulations of large scale structure formation with a detailed model of metal
enrichment, in which outflows expand anisotropically along the direction of
least resistance, distributing metals into the IGM. The metals carried by the
outflows are generated by two separate stellar populations: stars located near
the central AGN, and stars located in the greater galaxy. Using this algorithm,
we performed a series of 5 simulations of the propagation of AGN-driven
outflows in a cosmological volume of size (128/h Mpc)^3 in a Lambda-CDM
universe, and analyze the resulting metal enrichment of the IGM. We found that
the metallicity induced in the IGM is greatly dominated by AGNs having
bolometric luminosity L > 10^9 L_sun, sources with 10^8 < L / L_sun < 10^9
having a negligible contribution. Our simulations produced an average IGM
metallicity of [O/H] = -5 at z = 5.5, which then rises gradually, and remains
relatively flat at a value [O/H] = -2.8 between z = 2 and z = 0. The ejection
of metals from AGN host galaxies by AGN-driven outflows is found to enrich the
IGM to > 10 - 20% of the observed values, the number dependent on redshift. The
enriched IGM volume fractions are small at z > 3, then rise rapidly to the
following values at z = 0: 6 - 10% of the volume enriched to [O/H] > -2.5, 14 -
24% volume to [O/H] > -3, and 34 - 45% volume to [O/H] > -4. At z > 2, there is
a gradient of the induced enrichment, the metallicity decreasing with
increasing IGM density, enriching the underdense IGM to higher metallicities, a
trend more prominent with increasing anisotropy of the outflows. This can
explain observations of metal-enriched low-density IGM at z = 3 - 4.
[25]
oai:arXiv.org:0905.3560 [pdf] - 24520
Anisotropic AGN Outflows and Enrichment of the Intergalactic Medium
Submitted: 2009-05-21
We investigate the cosmological-scale influence of outflows driven by AGNs on
metal enrichment of the intergalactic medium. AGNs are located in dense
cosmological structures which tend to be anisotropic. We designed a
semi-analytical model for anisotropic AGN outflows which expand away along the
direction of least resistance. This model was implemented into a cosmological
numerical simulation algorithm for simulating the growth of large-scale
structure in the universe. Using this modified algorithm, we perform a series
of 9 simulations inside cosmological volumes of size $(128 h^{-1}{\rm Mpc})^3$,
in a concordance $\Lambda$CDM universe, varying the opening angle of the
outflows, the lifetimes of the AGNs, their kinetic fractions, and their level
of clustering. For each simulation, we compute the volume fraction of the IGM
enriched in metals by the outflows. The resulting enriched volume fractions are
relatively small at $z \gtrsim 2.5$, and then grow rapidly afterward up to $z =
0$. We find that AGN outflows enrich from 65% to 100% of the entire universe at
the present epoch, for different values of the model parameters. The enriched
volume fraction depends weakly on the opening angle of the outflows. However,
increasingly anisotropic outflows preferentially enrich underdense regions, a
trend found more prominent at higher redshifts and decreasing at lower
redshifts. The enriched volume fraction increases with increasing kinetic
fraction and decreasing AGN lifetime and level of clustering.
[26]
oai:arXiv.org:0902.4093 [pdf] - 21709
Behaviour of relativistic black hole accretion sufficiently close to the
horizon
Submitted: 2009-02-24
This work introduces a novel formalism to investigate the role of the spin of
astrophysical black holes in determining the behaviour of matter falling onto
such accretors. Equations describing the general relativistic hydrodynamic
accretion flow in the Kerr metric are formulated, and stationary solutions for
such flow equations are provided. The accreting matter may become
multi-transonic, allowing a stationary shock to form for certain initial
boundary conditions. Such a shock determines the disc geometry and can drive
strong outflows. The properties of matter extremely close to the event horizon
are studied as a function of the Kerr parameter, leading to the possibility of
detecting a new spectral signature of black hole spin.
[27]
oai:arXiv.org:0809.0515 [pdf] - 15918
Anisotropic AGN Outflows Filling The Cosmological Volume
Submitted: 2008-09-02
We simulate anisotropic outflows of AGN, and investigate the large-scale
impact of the cosmological population of AGN outflows over the Hubble time by
performing N-body LambdaCDM simulations. Using the observed quasar luminosity
function to get the redshift and luminosity distribution, and analytical models
for the outflow expansion, AGNs are allowed to evolve in a cosmological volume.
By the present epoch, 13 - 25% of the total volume is found to be pervaded by
AGN outflows, with 10^{-9} G magnetic field.
[28]
oai:arXiv.org:0707.1533 [pdf] - 2968
The Fate of Dwarf Galaxies in Clusters and the Origin of Intracluster
Stars. I. Isolated Clusters
Submitted: 2007-07-10, last modified: 2008-07-02
The main goal of this paper is to compare the relative importance of
destruction by tides, vs. destruction by mergers, in order to assess if tidal
destruction of dwarf galaxies in clusters is a viable scenario for explaining
the origin of intracluster stars. We have designed a simple algorithm for
simulating the evolution of isolated clusters. The distribution of galaxies in
the cluster is evolved using a direct gravitational N-body algorithm combined
with a subgrid treatment of physical processes such as mergers, tidal
disruption, and galaxy harassment. Using this algorithm, we have performed a
total of 227 simulations. Our main results are (1) destruction of dwarf
galaxies by mergers dominates over destruction by tides, and (2) the
destruction of dwarf galaxies by tides is sufficient to explain the observed
intracluster light in clusters.
[29]
oai:arXiv.org:0806.1956 [pdf] - 13474
Large Scale Impact of the Cosmological Population of Expanding Radio
Galaxies
Submitted: 2008-06-11
We seek to compute the fraction of the volume of the Universe filled by
expanding cocoons of the cosmological population of radio galaxies over the
Hubble time as well as the magnetic field infused by them, in order to assess
their importance in the cosmic evolution of the Universe. Using N-body
$\Lambda$CDM simulations, radio galaxies distributed according to the observed
radio luminosity function are allowed to evolve in a cosmological volume as
using well defined prescriptions for their expansion. We find that the radio
galaxies permeate $10 - 30%$ of the total volume with $\sim 10^{-8}$ G magnetic
field by the present epoch.
[30]
oai:arXiv.org:0706.1221 [pdf] - 2028
Radial Velocities of Six OB Stars
Boyajian, T. S.;
Gies, D. R.;
Baines, E. K.;
Barai, P.;
Grundstrom, E. D.;
McSwain, M. V.;
Parks, J. R.;
Riddle, R. L.;
Ryle, W. T.;
Wingert, D. W.
Submitted: 2007-06-08
We present new results from a radial velocity study of six bright OB stars
with little or no prior measurements. One of these, HD 45314, may be a
long-period binary, but the velocity variations of this Be star may be related
to changes in its circumstellar disk. Significant velocity variations were also
found for HD 60848 (possibly related to nonradial pulsations) and HD 61827
(related to wind variations). The other three targets, HD 46150, HD 54879, and
HD 206183, are constant velocity objects, but we note that HD 54879 has
H$\alpha$ emission that may originate from a binary companion. We illustrate
the average red spectrum of each target.
[31]
oai:arXiv.org:astro-ph/0611689 [pdf] - 87075
Testing Models of Radio Galaxy Evolution and the Cosmological Impact of
FR II Radio Galaxies
Submitted: 2006-11-21
We investigate aspects of the cosmological evolution of FR II radio galaxies,
focusing first on the abilities of models to match data for linear-sizes, radio
powers, redshifts and spectral indices. Here we consider modifications to the
theoretical models we had treated earlier, primarily by accounting for the
growth of the radius of hotspots with source size. Better fits to the
distributions of most of the data in three low-frequency surveys can be found
with sensible choices of model parameters but no model yet considered gives a
good match to all of the survey data simultaneously, nor does any do a good job
of producing the spectral index distributions. The observational datasets are
too small to completely discriminate among the models. We calculate the volume
fraction of the ``relevant universe'' cumulatively occupied by the expanding
radio galaxy lobes over the quasar era, when these powerful radio galaxies were
much more common, and when they have been argued to play an important role in
triggering galaxy formation and spreading magnetic fields and metals. We found
the cumulative relevant volume filling factor of radio galaxies to be $\sim 5
%$, so we conclude that these impacts are smaller than previously estimated but
that they are still significant.
[32]
oai:arXiv.org:astro-ph/0510724 [pdf] - 77220
Testing Models of the Individual and Cosmological Evolutions of Powerful
Radio Galaxies
Submitted: 2005-10-25, last modified: 2006-07-19
We seek to develop an essentially analytical model for the evolution of
Fanaroff-Riley Class II radio galaxies as they age individually and as their
numbers vary with cosmological epoch. Such modeling is required in order to
probe in more detail the impact of radio galaxies on the growth of structures
in the universe, which appears likely to have been quite significant at z > 1.
In this first paper of a series we compare three rather sophisticated
analytical models for the evolution of linear size and lobe power of FR II
radio galaxies, those of Kaiser, Dennett-Thorpe & Alexander (1997), Blundell,
Rawlings & Willott (1999), and Manolakou & Kirk (2002). We perform
multi-dimensional Monte Carlo simulations in order to compare the predictions
of each model for radio powers, sizes, redshifts and spectral indices with
data. The observational samples used here are the low frequency radio surveys,
3CRR, 6CE and 7CRS, which are flux-limited and complete. We search for and
describe the best parameters for each model, after doing statistical tests on
them. We find that no existing model can give acceptable fits to all the
properties of the surveys considered, although the Kaiser, Dennett-Thorpe &
Alexander (1997) model gives better overall results than do the Manolakou &
Kirk (2002) or Blundell, Rawlings & Willott (1999) models for most of the tests
we performed. We suggest ways in which these models may be improved.
[33]
oai:arXiv.org:astro-ph/0411096 [pdf] - 68695
Expansion of radio galaxies in a cosmologically evolving medium:
Possible implications for the cosmic star-formation history
Submitted: 2004-11-03
We compare earlier estimates of the volumes filled by lobes of radio galaxies
during the quasar era based upon non-evolving ambient media with new ones
assuming a strong cosmological evolution of the ambient medium. If the sources
remain active for over 10^8 years the volumes filled by them are found to be
comparable for the two scenarios. This strengthens our earlier inference that
much of the cosmic web of gaseous filaments, the site of galaxy formation, was
probably permeated by radio lobes during the quasar era and this could have
triggered extensive star formation and made large contributions to the spread
of magnetic fields and metals through the universe by z ~ 2.
[34]
oai:arXiv.org:astro-ph/0408170 [pdf] - 66661
Dependence of General Relativistic Accretion on Black Hole Spin
Submitted: 2004-08-10
We propose a novel approach to the investigation of how black hole spin
influences the behavior of astrophysical accretion. We examine the terminal
behavior of general relativistic, multi-transonic, hydrodynamic advective
accretion in the Kerr metric, by analyzing various dynamical as well as
thermodynamic properties of accreting matter extremely close to the event
horizon, as a function of Kerr parameter $a$. The examination of the transonic
properties of such flows can be useful in self-consistently discriminating
between prograde and retrograde relativistic accretion, and in studying the
spectral signatures of black hole rotation.
[35]
oai:arXiv.org:astro-ph/0309734 [pdf] - 59561
Mass and Angular Momentum Transfer in the Massive Algol Binary RY Persei
Barai, P.;
Gies, D. R.;
Choi, E.;
Das, V.;
Deo, R.;
Huang, W.;
Marshall, K.;
McSwain, M. V.;
Ogden, C.;
Osterman, A.;
Riddle, R. L.;
Seymour, J. L.;
Wingert, D. W.;
Kaye, A. B.;
Peters, G. J.
Submitted: 2003-09-26, last modified: 2004-04-02
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.