Showing papers in "arXiv: Astrophysics in 1993"

Using FU Orionis Outbursts to Constrain Self-Regulated Protostellar Disk Models

Abstract: One dimensional, convective, vertical structure models and one dimensional, time dependent, radial diffusion models are combined to create a self-consistent picture in which FU~Orionis outbursts occur in young stellar objects (YSOs) as the result of a large scale, self-regulated, thermal ionization instability in the surrounding protostellar accretion disk. By fitting the results of time dependent disk models to observed time scales of FU~Orionis events, we estimate the magnitude of the effeciency of the effective viscous stress in the inner disk (r < 1 AU) to be, in accordance with the ad hoc ``alpha'' prescription, 10^{-4} where hydrogen is neutral and 10^{-3} where hydrogen is ionized. We hypothesize that all YSOs receive infall onto their outer disks which is steady (or slowly declining with time) and that FU~Orionis outbursts are self-regulated, disk outbursts which occur only in systems which transport matter inward at a rate sufficiently high to cause hydrogen to be ionized in the inner disk. We estimate a critical mass flux of 5x10^{-7} solar masses / year {\it independent of the magnitude of alpha} for systems with one solar mass, three solar radius central objects. Infall accretion rates in the range of (1-10)x10^{-6} solar masses per year produce observed FU~Orionis time scales consistent with estimates of spherical molecular cloud core collapse rates.

Element Diffusion in the Solar Interior

Abstract: We study the diffusion of helium and other heavy elements in the solar interior by solving exactly the set of flow equations developed by Burgers for a multi-component fluid, including the residual heat-flow terms. No approximation is made concerning the relative concentrations and no restriction is placed on the number of elements considered. We give improved diffusion velocities for hydrogen, helium, oxygen and iron, in the analytic form derived previously by Bahcall and Loeb. These expressions for the diffusion velocities are simple to program in stellar evolution codes and are expected to be accurate to $\sim 15\%$. Our complete treatment of element diffusion can be directly incorporated in a standard stellar evolution code by means of an exportable subroutine, but, for convenience, we also give simple analytical fits to our numerical results.

Biasing and Hierarchical Statistics in Large-scale Structure

Abstract: In the current paradigm there is a non-trivial bias expected in the process of galaxy formation. Thus, the observed statistical properties of the galaxy distribution do not necessarily extend to the underlying matter distribution. Gravitational evolution of initially Gaussian seed fluctuations predicts that the connected moments of the matter fluctuations exhibit a hierarchical structure, at least in the limit of small dispersion. This same hierarchical structure has been found in the galaxy distribution, but it is not clear to what extent it reflects properties of the matter distribution or properties of a galaxy formation bias. In this paper we consider the consequences of an arbitrary, effectively local biasing transformation of a hierarchical underlying matter distribution. We show that a general form of such a transformation preserves the hierarchical properties and the shape of the dispersion in the limit of small fluctuations, i.e. on large scales, although the values of the

Ultraviolet Radiation from Evolved Stellar Populations -- I. Models

Abstract: This series of papers comprises a systematic exploration of the hypothesis that the far ultraviolet radiation from star clusters and elliptical galaxies originates from extremely hot horizontal-branch (HB) stars and their post-HB progeny. This first paper presents an extensive grid of calculations of stellar models from the Zero Age Horizontal Branch through to a point late in post-HB evolution or a point on the white dwarf cooling track. We use the term `Extreme Horizontal Branch' (EHB) to refer to HB sequences of constant mass that do not reach the thermally-pulsing stage on the AGB. These models evolve after core helium exhaustion

Dipole Anisotropy in the COBE DMR First-Year Sky Maps

Abstract: We present a determination of the cosmic microwave background dipole amplitude and direction from the COBE Differential Microwave Radiometers (DMR) first year of data. Data from the six DMR channels are consistent with a Doppler-shifted Planck function of dipole amplitude Delta T = 3.365 +/-0.027 mK toward direction (l,b) = (264.4 +/- 0.3 deg, 48.4 +/- 0.5 deg). The implied velocity of the Local Group with respect to the CMB rest frame is 627 +/- 22 km/s toward (l,b) = (276 +/- 3 deg, 30 +/- 3 deg). DMR has also mapped the dipole anisotropy resulting from the Earth's orbital motion about the Solar system barycenter, yielding a measurement of the monopole CMB temperature at 31.5, 53, and 90 GHz, to be 2.75 +/- 0.05 K.

A Family of Models for Spherical Stellar Systems

Abstract: We describe a one-parameter family of models of stable spherical stellar systems in which the phase-space distribution function depends only on energy. The models have similar density profiles in their outer parts ($\rho\propto r^{-4}$) and central power-law density cusps, $\rho\propto r^{3-\eta}$, $0<\eta\le 3$. The family contains the Jaffe (1983) and Hernquist (1990) models as special cases. We evaluate the surface brightness profile, the line-of-sight velocity dispersion profile, and the distribution function, and discuss analogs of King's core-fitting formula for determining mass-to-light ratio. We also generalize the models to a two-parameter family, in which the galaxy contains a central black hole; the second parameter is the mass of the black hole. Our models can be used to estimate the detectability of central black holes and the velocity-dispersion profiles of galaxies that contain central cusps, with or without a central black hole.

Nonlinear Instabilities in Shock-Bounded Slabs

Abstract: (substantial changes to section 3.2, otherwise minor) We present an analysis of the hydrodynamic stability of a cold slab bounded by two accretion shocks. Previous numerical work has shown that when the Mach number of the shock is large the slab is unstable. Here we show that to linear order both the bending and breathing modes of such a slab are stable. However, nonlinear effects will tend to soften the restoring forces for bending modes, and when the slab displacement is comparable to its thickness this gives rise to a nonlinear instability. The growth rate of the instability, above this threshold but for small bending angles, is $\sim c_sk (k\eta)^{1/2}$, where $\eta$ is the slab displacement. When the bending angle is large the slab will contain a local vorticity comparable to $c_s/L$, where $L$ is the slab thickness. We discuss the implications of this work for gravitational instabilities of slabs. Finally, we examine the cases of a decelerating slab bounded by a single shock and a stationary slab bounded on one side by thermal pressure. The latter case is stable, but appears to be a special case. The former case is subject to a nonlinear overstability driven by deceleration effects. We conclude that shock bounded slabs with a high density compression ratio generically produce substructure with a strong local shear, a bulk velocity dispersion like the sound speed in the cold layer and a characteristic scale comparable to the slab thickness.

Very High Energy Gamma-Rays from AGN: Cascading on the Cosmic Background Radiation Fields and the Formation of Pair Halos

Abstract: Recent high energy gamma-ray observations (E>100 MeV) of blazar AGN show emission spectra with no clear upper energy cutoff. AGN, considered to be possible sources for the highest energy cosmic rays, may have emission extending well into the VHE (very high energy, E>100 GeV) domain. Because VHE gamma-rays are absorbed by pair production on the intergalactic background radiation, much of this emission is not visible. Electromagnetic cascades initiated by absorbed VHE gamma-rays, however, may be observable. Since the velocities of e+e- pairs produced in a cascade are likely to be isotropized by an ambient random magnetic field, extended ``halos'' (R>1 Mpc) of pairs will be formed around AGN with VHE emission. Cascade radiation from these pair halos is emitted isotropically and is observable below a few TeV. The halo radiation can be distinguished by its characteristic variation in spectrum and intensity with angular distance from the central source. This variation depends weakly on the details of the central source model, e.g., the orientation and beaming/opening angle of an emitting jet. Limiting or determining the intensity of the pair halo can thus serve as a model-independent bound on or measure of the VHE power of AGN. Next-generation Cherenkov telescopes may be able to image a pair halo.

Radio Transients from Gamma-Ray Bursters

Abstract: The rapid time variability of gamma-ray bursts implies the sources are very compact, and the peak luminosities are so high that some matter must be ejected at ultra-relativistic speeds. The very large Lorentz factors of the bulk flow are also indicated by the very broad and hard spectra. It is natural to expect that when the relativistic ejecta interact with the interstellar matter a strong synchrotron radio emission is generated, as is the case with supernova remnants and radio galaxies. We estimate that the strongest gamma-ray bursts may be followed by radio transients with peak fluxes as high as 0.1 Jy. The time of peak radio emission depends on the distance scale; it is less than a minute if the bursts are in the galactic halo, and about a week if the bursts are at cosmological distances.

Standard and Non-Standard Plasma Neutrino Emission Revisited

Abstract: On the basis of Braaten and Segel's representation of the electromagnetic dispersion relations in a QED plasma we check the numerical accuracy of several published analytic approximations to the plasma neutrino emission rates. As we find none of them satisfactory we derive a new analytic approximation which is accurate to within 4\%\ where the plasma process dominates. The correct emission rates in the parameter regime relevant for the red giant branch in globular clusters are larger by about $10-20$\% than those of previous stellar evolution calculations. Therefore, the core mass of red giants at the He flash is larger by about $0.005\M_\odot$ or 1\% than previously thought. Our bounds on neutrino magnetic dipole moments remain virtually unchanged.

The Effect of Dissipation on the Shapes of Dark halos

Abstract: The dissipative infall of gas during the formation of a galaxy modifies the density profile and shape of the dark halo. Gas dissipates energy radiatively and sinks to the center of the dark halo forming the luminous part of a galaxy. The resulting central density enhancement can alter the halo's orbital distribution. We simulate dissipative infall inside of an initially triaxial N-body dark halo by slowly growing a potential in the center of the particle distribution. The dark halo transforms from a prolate-triaxial halo ($T \sim 0.8$) to an oblate-triaxial halo ($T \sim 0.5$) while approximately preserving the flattening ($c/a \sim 0.5$). The main implication is that dark halos are rounder and more oblate than previous predictions of purely collisionless simulations with the new constraint that $b/a \gapp 0.7$. At the same time, the distribution of intrinsic flattenings of dark halos ($\langle c/a \rangle = 0.5$, $\langle (c/a)^2 \rangle ^{1/2} = 0.15$) is preserved during the period of baryonic infall. The oval distortions of disk galaxies should therefore be slightly less than original predictions from collisionless dark halos. The predicted distribution of shapes for dark halos from cosmological N-body simulations is in better agreement with distribution of ellipticities of elliptical galaxies if we assume all halos and galaxies become oblate-triaxial in response to baryonic dissipation. The observed distribution of kinematic misalignment angles is also consistent with this shape distribution if

On the Relationship between BL-Lacertae Objects and FRI Radio Galaxies

Abstract: The idea that Fanaroff-Riley Class~I Radio Galaxies contain core jets with Lorentz factors of the order of a few and are the parent population for BL~Lac Objects is examined, with particular reference to the data on two FRI Radio Galaxies, NGC~315 and NGC~6251. Conservation laws for an entraining relativistic jet are derived and are used to determine a relationship between Mach number and velocity for initially relativistic jets. One immediate consequence of this relationship is that, if an initially relativistic jet is decelerated to a transonic Mach number between, say 0.5 and 1.5, then its velocity at this point is between 0.3 and 0.7 times the speed of light. Analysis of the milliarcsecond and arcsecond data on NGC~315 and NGC~6251 shows that the Lorentz factors of the core jets in these galaxies can be as high as 2-4 provided that the jet pressures are not much more than an order of magnitude higher than the minimum pressures. Such Lorentz factors are more plausible if the jets are inclined at $\sim 30$ degrees to the line of sight, although larger angles would also be possible. For core Lorentz factors in the range 2-4, the jets are mildly relativistic on the large scale and their velocities are consistent with the jet brightness asymmetries being caused by Doppler beaming. The Lorentz factors required for significant beaming of core jets may be reconciled with subluminal motions of knots in the core if these are reverse shocks advected by the jet. Moderate jet Lorentz factors are possible if the shocks are normal and higher Lorentz factors are possible if the shocks are oblique. This paper provides support for the idea that the transition from two-sided to one-sided jets across the FRI/II break is due to the transition from

Structure Formation with Cold + Hot Dark Matter

Abstract: We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an $\Omega=1$ cosmological model with a mixture of Cold plus Hot Dark Matter (C+HDM) having $\Omega_{\rm cold}=0.6$, $\Omega_ u=0.3$, and $\Omega_{\rm baryon}=0.1$. We present analytic fits to the C+HDM power spectra for both cold and hot ($ u$) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included six times as many neutrino particles as cold particles. Our simulation boxes were 14, 50, and 200~Mpc cubes (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$); we also did comparison simulations for Cold Dark Matter (CDM) in a 50~Mpc box. C+HDM with linear bias factor $b=1.5$ is consistent both with the COBE data and with the galaxy correlations we calculate. We find the number of halos as a function of mass and redshift in our simulations; our results for both CDM and C+HDM are well fit by a Press-Schechter model. The number density of galaxy-mass halos is smaller than for CDM, especially at redshift $z>2$, but the numbers of cluster-mass halos are comparable. We also find that on galaxy scales the neutrino velocities and flatter power spectrum in C+HDM result in galaxy pairwise velocities that are in good agreement with the data, and about 30\% smaller than in CDM with the same biasing factor. On scales of several tens of Mpc, the C+HDM streaming velocities are considerably larger than CDM. Thus C+HDM looks promising as a model of structure formation.

Oxygen Abundances in Low Surface Brightness Disk Galaxies

Abstract: The oxygen abundances in the \HII regions of a sample of low surface brightness (LSB) disk galaxies are presented. In general, LSB galaxies are found to be metal poor ($Z < {1 \over 3} Z\solar$). Indeed, some LSB galaxies rival the lowest abundance extragalactic objects known. These low metallicities indicate that LSB galaxies evolve slowly, forming relatively few stars over a Hubble time. The low metallicities of LSB galaxies occur even though many are comparable in size and mass to the prominent spirals which define the Hubble sequence. As well as being low in surface brightness, these galaxies tend to be isolated. This suggests that surface mass density and environment are more relevant to galaxy evolution than gross size. Despite the low surface brightness of the disks, massive ($M > 60 M\solar$) stars are inferred to be present and no abnormality of the IMF is indicated. Many low excitation \HII regions exist at low metallicity in LSB galaxies, and the ionization parameter is not tightly correlated with metallicity. However, there does seem to be a significant envelope of maximum ionization at a given metallicity.

The Environment of Lyman a Absorbers in the Sightline toward 3C273

Abstract: We present new ground-based data following up on the HST discovery of low-redshift Lya absorption in the sight-line to the quasar 3C273. Narrow-band filter observations show that there are no H II regions within a 12 kpc radius of the line-of-sight to the quasar, at the velocities of three of the absorbers. Broad-band imaging shows that there are no dwarf galaxies at Virgo distances with absolute magnitude above MB~-13.5 and within a radius of 40 kpc. We present fiber spectroscopy of galaxies within a radius of 1 deg, down to an apparent magnitude of B~19. We show that the absorbers are definitely not distributed at random with respect to the galaxies, but also that the absorber-galaxy correlation function is not as strong as the galaxy-galaxy correlation function on large scales. Our data are consistent with the hypothesis that all galaxies more luminous than 1/10 L* have effective cross-sections (for association with absorbers with Log(NH)>13.0), of between 0.5 and 1 Mpc. We also show a clear case of a Lya absorber which has no galaxy brighter than MB=-18 within a projected distance of 4.8 Mpc, and discuss the possibility that Lya absorbers are destroyed in a rich galaxy environment.

Dynamics with a non-standard inertia-acceleration relation: an alternative to dark matter

Abstract: We investigate particle laws of motion derived from nonstandard kinetic actions of a special form. We are guided by a phenomenological scheme--the modified dynamics (MOND)--that imputes the mass discrepancy observed in galactic systems to a departure from Newtonian dynamics below a certain acceleration scale a0. In the limit a0 goes to 0 the theory goes to Newtonian dynamics. In the opposite limit the action becomes proportional to 1/a0. Galilei-invariant such theories must be strongly non-local; this is a blessing, as such theories need not suffer from the illnesses that are endemic to higher-derivative theories. We discuss the possibility that such a modified law of motion is an effective theory resulting from the elimination of degrees of freedom pertaining to the universe at large (the near equality a0=cH0 being a trace of that connection). A virial relation for bounded trajectories is derived. Exact solutions are obtained for circular orbits, which pertain to rotation curves of disk galaxies. We also explore, in passing, theories that depart from the conventional Newtonian dynamics for very low frequencies. The present formulation contrasts, on important accounts, with earlier, modified-gravity formulations.

Weakly Non-Linear Gaussian Fluctuations and the Edgeworth Expansion

Abstract: We calculate the cosmological evolution of the 1-point probability distribution function (PDF), using an analytic approximation that combines gravitational perturbation theory with the Edgeworth expansion of the PDF. Our method applies directly to a smoothed mass density field or to the divergence of a smoothed peculiar velocity field, provided that rms fluctuations are small compared to unity on the smoothing scale, and that the primordial fluctuations that seed the growth of structure are Gaussian. We use this `Edgeworth approximation' to compute the evolution of $ $ and $ $; these measures are similar to the skewness and kurtosis of the density field, but they are less sensitive to tails of the probability distribution, so they may be more accurately estimated from surveys of limited volume. We compare our analytic calculations to cosmological N-body simulations in order to assess their range of validity. When $\sigma \ll 1$, the numerical simulations and perturbation theory agree precisely, demonstrating that the N-body method can yield accurate results in the regime of weakly non-linear clustering. We show analytically that `biased' galaxy formation preserves the relation $ \propto ^2$ predicted by second-order perturbation theory, provided that the galaxy density is a local function of the underlying mass density.

Skewness and Kurtosis in Large-Scale Cosmic Fields

Abstract: In this paper, I present the calculation of the third and fourth moments of both the distribution function of the large--scale density and the large--scale divergence of the velocity field, $\theta$. These calculations are made by the mean of perturbative calculations assuming Gaussian initial conditions and are expected to be valid in the linear or quasi linear regime. The moments are derived for a top--hat window function and for any cosmological parameters $\Omega$ and $\Lambda$. It turns out that the dependence with $\Lambda$ is always very weak whereas the moments of the distribution function of the divergence are strongly dependent on $\Omega$. A method to measure $\Omega$ using the skewness of this field has already been presented by Bernardeau et al. (1993). I show here that the simultaneous measurement of the skewness and the kurtosis allows to test the validity of the gravitational instability scenario hypothesis. Indeed there is a combination of the first three moments of $\theta$ that is almost independent of the cosmological parameters $\Omega$ and $\Lambda$, $${( -3 ^2) \over ^2}\approx 1.5,$$ (the value quoted is valid when the index of the power spectrum at the filtering scale is close to -1) so that any cosmic velocity field created by gravitational instabilities should verify such a property.

Evolution of One-Point Distributions from Gaussian Initial Fluctuations

Abstract: We study the quasilinear evolution of the one-point probability density functions (PDFs) of the smoothed density and velocity fields in a cosmological gravitating system beginning with Gaussian initial fluctuations. Our analytic results are based on the Zel'dovich approximation and laminar flow. A numerical analysis extends the results into the multistreaming regime using the smoothed fields of a CDM N-body simulation. We find that the PDF of velocity, both Lagrangian and Eulerian, remains Gaussian under the laminar Zel'dovich approximation, and it is almost indistinguishable from Gaussian in the simulations. The PDF of mass density deviates from a normal distribution early in the quasilinear regime and it develops a shape remarkably similar to a lognormal distribution with one parameter, the \rms density fluctuation $\sigma$. Applying these results to currently available data we find that the PDFs of the velocity and density fields, as recovered by the \pot\ procedure from observed velocities assuming $\Omega=1$, or as deduced from a redshift survey of \iras\ galaxies assuming that galaxies trace mass, are consistent with Gaussian initial fluctuations.

Structural Characteristics and Stellar Composition of Low Surface Brightness Disk Galaxies

Abstract: We present UBVI surface photometry of a sample of low surface brightness (LSB) disk galaxies. LSB disk galaxies are fairly well described as exponential disks with no preferred value for either scale length, central surface brightness, or rotational velocity. Indeed, the distribution of scale lengths is indistinguishable from that of high surface brightness spirals, indicating that dynamically similar galaxies (e.g., those with comparable Rv^2) exist over a large range in surface density. These LSB galaxies are strikingly blue. The complete lack of correlation between central surface brightness and color rules out any fading scenario. Similarly, the oxygen abundances inferred from HII region spectra are uncorrelated with color so the low metallicities are not the primary cause of the blue colors. While these are difficult to interpret in the absence of significant star formation, the most plausible scenario is a stellar population with a young mean age stemming from late formation and subsequent slow evolution. These properties suggest that LSB disks formed from low initial overdensities with correspondingly late collapse times.

Moments of the Counts Distribution in the 1.2 Jy IRAS Redshift Survey

Abstract: We derive the volume-averaged 2, 3, 4, and 5-point correlation functions from the moments of the Count probability distribution function of a redshift survey of IRAS galaxies, and find them all to be reasonably well-described by power laws. Weak systematic effects with the sample size provide evidence for stronger clustering of galaxies of higher luminosity on small scales. Nevertheless, remarkably tight relationships hold between the correlation functions of different order. In particular, the ``normalized" skewness defined by the ratio $S_3\equiv \bar{\xi_3} / \bar{\xi_2}^2$ varies at most weakly with scale in the range $0.1 < \bar{\xi_2} < 10$. That is, $S_3$ is close to constant ($=1.5\pm 0.5$) from weakly to strongly non-linear scales. Furthermore, we find that the void probability function obeys a scaling relation with density to great precision, in accord with the scale-invariance hypothesis ($\bar{\xi_N}\propto\bar{\xi_2}^{N-1}$).

Detection of Signature Consistent with Cosmological Time Dilation in Gamma-Ray Bursts

Abstract: If gamma-ray bursters are at cosmological distances - as suggested by their isotropic distribution on the sky and by their number-intensity relation - then the burst profiles will be stretched in time, by an amount proportional to the redshift, 1 + $z$. We have tested data from the {\it Compton} Gamma Ray Observatory's Burst and Transient Source Experiment (BATSE) for such time dilation. Our measures of time scale are constructed to avoid selection effects arising from intensity differences by rescale all bursts to fiducial levels of peak intensity and noise bias. The three tests involved total count rate above background, wavelet decomposition, and alignment of the highest peaks. In all three tests, the dim bursts are stretched by a factor of about two relative to the bright ones, over seven octaves of time scale. We calibrated the measurements by dilating synthetic bursts that approximate the temporal characteristics of bright BATSE bursts. Results are consistent with bursts of BATSE's peak-flux completeness limit being at cosmological distances corresponding to $z \sim 1$, and thus with independent cosmological interpretations of the BATSE number-intensity relation.

Skewness Induced by Gravity

Abstract: Large-scale structures, observed today, are generally believed to have grown from random, small-amplitude inhomogeneities, present in the early Universe. We investigate how gravitational instability drives the distribution of these fluctuations away from the initial state, assumed to be Gaussian. Using second order perturbation theory, we calculate the skewness factor, $S_3 \equiv \langle \delta^3 \rangle \, /\, \langle \delta^2 \rangle^2$. Here the brackets, $\langle \ldots \rangle$, denote an ensemble average, and $\delta$ is the density contrast field, smoothed with a low pass spatial filter. We show that $S_3$ decreases with the slope of the fluctuation power spectrum; it depends only weakly on $\Omega$, the cosmological density parameter. We compare perturbative calculations with N-body experiments and find excellent agreement over a wide dynamic range. If galaxies trace the mass, measurements of $S_3$ can be used to distinguish models with Gaussian initial conditions from their non-Gaussian alternatives.

Gamma-Ray Bursts: Multiwaveband Spectral Predictions for Blast Wave Models

Abstract: In almost any scenario for 'cosmological' gamma-ray bursts (and in many models where they originate in our own Galaxy), the initial energy density is so large that the resulting relativistic plasma expands with $v\sim c$ producing a blast wave ahead of it and a reverse shock moving into the ejecta, as it ploughs into the external medium. We evaluate the radiation expected from these shocks,for both cosmological and galactic bursts, for various assumptions about the strength of the magnetic field and the particle acceleration mechanisms in the shocks. The spectra are evaluated over the whole range from the IR to $>$ GeV, and are compared with the variety of spectral behavior reported by BATSE, and with the X-ray and optical constraints. For bursts of duration $\simg 1\s$ acceptable $\gamma$-ray spectra and $L_x/L_\gamma$ ratios are readily obtained for 'cosmological' models. Blast waves in galactic models can produce bursts of similar gamma-ray fluence and duration, but they violate the X-ray paucity constraint, except for the shorter bursts ($\siml 1\s$). We discuss the prospects for using O/UV and X-ray observations to discriminate among alternative models.

Simulation of Thick Accretion Disks with Standing Shocks by Smoothed Particle Hydrodynamics

Abstract: We present results of numerical simulation of inviscid thick accretion disks and wind flows around black holes. We use Smoothed Particle Hydrodynamics (SPH) technique for this purpose. Formation of thick disks are found to be preceded by shock waves travelling away from the centrifugal barrier. For a large range of the parameter space, the travelling shock settles at a distance close to the location obtained by a one-and-a-half dimensional model of inviscid accretion disks. Occasionally, it is observed that accretion processes are aided by the formation of oblique shock waves, particularly in the initial transient phase. The post-shock region (where infall velocity suddenly becomes very small) resembles that of the usual model of thick accretion disk discussed in the literature, though they have considerable turbulence. The flow subsequently becomes supersonic before falling into the black hole. In a large number of cases which we simulate, we find the formation of strong winds which are hot and subsonic when originated from the disk surface very close to the black hole but become supersonic within a few tens of the Schwarzschild radius of the blackhole. In the case of accretion of high angular momentum flow, very little amount of matter is accreted directly onto the black hole. Most of the matter is, however, first squeezed to a small volume close to the black hole, and subsequently expands and is expelled as a strong wind. It is quite possible that this expulsion of matter and the formation of cosmic radio jets is aided by the shock heating in the inner parts of the accretion disks.

Properties of High-Redshift Lyman Alpha Clouds I. Statistical Analysis of the SSG Quasars

Abstract: Techniques for the statistical analysis of the \Lya\ forest in high redshift quasars are developed, and applied to the low resolution (25 \AA) spectra of 29 of the 33 quasars in the Schneider-Schmidt-Gunn (SSG) sample.We find that the mean absorption increases with $z$ approximately as a power law $(1+z)^{\gamma+1}$ with $\gamma = 2.46\pm 0.37$. The mean ratio of \Lya\ to Lyman $\beta$ absorption in the clouds is $0.476\pm 0.054$. We also detect, and obtain ratios, for Lyman $\beta$, $\gamma$, and possibly $\epsilon$. We are also able to quantify the fluctuations of the absorption around its mean, and find that these are comparable to, or perhaps slightly larger than, that expected from an uncorrelated distribution of clouds. The techniques in this paper, which include the use of bootstrap resampling of the quasar sample to obtain estimated errors and error covariances, and a mathematical treatment of absorption from a (possibly non-uniform) stochastic distribution of lines, should be applicable to future, more extensive, data sets.

The Evolution of the Globular Cluster System in a Triaxial Galaxy: Can a Galactic Nucleus Form by Globular Cluster Capture?

Abstract: Dynamical friction due to field stars and tidal disruption caused by a central nucleus are crucial in determining the evolution of the globular cluster system in an elliptical galaxy. In this paper I examine the possibility that some of galactic nuclei have been formed by frictionally decayied globular clusters moving in a triaxial potential. The initial rapid growth of the nucleus, due mainly to massive clusters on box orbits falling in a short time scale into the galactic centre, is found to be later slowed by tidal disruption induced by the nucleus itself on less massive clusters in the way described by Ostriker, Binney & Saha. The efficiency of dynamical friction is such to carry to the centre of the galaxy enough globular cluster mass available to form a compact nucleus, but the actual modes of its collisionless formation remains to be investigated. The mass of the resulting nucleus is determined by the mutual feedback of the two mentioned processes, together with the knowlegde of the initial spatial, velocity and mass distributions of the globular cluster family. Moreover, we discuss the possibility that the globular cluster fall to the galactic centre has been cause of primordial violent galactic activity. An application of the model to M31 is presented.

The Nonlinear Evolution of Rare Events

Abstract: In this paper I consider the nonlinear evolution of a rare density fluctuation in a random density field with Gaussian fluctuations, and I rigorously show that it follows the spherical collapse dynamics applied to its mean initial profile. This result is valid for any cosmological model and is independent of the shape of the power spectrum. In the early stages of the dynamics the density contrast of the fluctuation is seen to follow with a good accuracy the form $$\delta=(1-\delta_L/1.5)^{-1.5}-1,$$ where $\delta_L$ is the linearly extrapolated overdensity. I then investigate the validity domain of the rare event approximation in terms of the parameter $ u=\delta_L/\sigma$ giving the initial overdensity scaled by the rms fluctuation at the same mass scale, and find that it depends critically on the shape of the power spectrum. When the power law index $n$ is lower than $-1$ the departure from the spherical collapse is expected to be small, at least in the early stages of the dynamics, and even for moderate values of $ u$ ($\vert u\vert\ge 2$). However, for $n\ge -1$ and whatever the value of $ u$, the dynamics seems to be dominated by the small--scale fluctuations and the subsequent evolution of the peak may not be necessarily correlated to its initial overdensity. I discuss the implications of these results for the nonlinear dynamics and the formation of astrophysical objects.

A Multiwavelength Study of Star Formation in the L1495E Cloud in Taurus

Abstract: We have carried out a deep (t=30000s) x-ray search of the eastern portion of the L1495 cloud centered on the well known weak line T Tauri star (WTTS) V410 Tau using the ROSAT PSPC. This deep exposure enabled a search for candidate pre-main sequence (PMS) objects in this cloud to a limit 20 times more sensitive than that typical of the fields examined with the Einstein searches. Despite assertions that the PMS population in Taurus-Auriga is nearly completely known, this x-ray survey revealed 8 new PMS objects in a region 50' in diameter, as compared to a previously known stellar population of 12 objects, including deeply embedded IRAS sources. Spectroscopic and photometric observations enable us to place these objects in the HR Diagram. The newly discovered objects are predominantly stars of spectral type M0 and later, and a large fraction (6/8) appear to be surrounded by circumstellar accretion disks as judged by their infrared excess and H alpha emission. We combined the data for these x-ray discovered objects with extant and new data for the previously identified PMS stars in this region to examine the history of star formation and the frequency distribution of stellar masses in this cloud. If the ``post ROSAT'' population is either complete or representative, we conclude (1) that star formation in L1495 East took place 1 x 10^6 yrs ago and that the spread in ages is small; (2) the frequency distribution of masses,

Shock Excited Maser Emission from the Supernova Remnant W28

Abstract: We present VLA observations in the hydroxyl line OH(1720 MHz) towards the supernova remnant W28. These high resolution images show a more extensive distribution of 1720 MHz emission than had previously been suspected from earlier single-dish observations. A total of 26 distinct maser spots were detected which lie on the edge of the radio continuum emission from W28 and are unresolved with our 12 arcsec resolution. Significant detection of linear and circular polarization in the maser spots is also made. The location of the maser spots, together with the physical conditions in the SNR and the adjacent molecular cloud, argue for a shock excitation mechanism. The possibilities of using the OH(1720 MHz) as a powerful diagnostic of the interaction of supernova remnants with molecular clouds is discussed.