Mass estimation in the outer regions of galaxy clusters

We present a technique for estimating the mass in the outskirts of galaxy clusters where the usual assumption of dynamical equilibrium is not valid. The method assumes that clusters form through hierarchical clustering and requires only galaxy redshifts and positions on the sky. We apply the method to dissipationless cosmological N-body simulations where galaxies form and evolve according to semi-analytic modelling. The method recovers the actual cluster mass profile within a factor of two to several megaparsecs from the cluster centre. This error originates from projection effects, sparse sampling, and contamination by foreground and background galaxies. In the absence of velocity biases, this method can provide an estimate of the mass-to-light ratio on scales ~1-10 Mpc/h where this quantity is still poorly known.Comment: 14 pages, 7 figures, MN LaTeX style, MNRAS, in pres

No Pulsar Kicks from Deformed Neutrinospheres

In a supernova core, magnetic fields cause a directional variation of the neutrino refractive index so that resonant flavor oscillations would lead to a deformation of the "neutrinosphere" for, say, tau neutrinos. The associated anisotropic neutrino emission was proposed as a possible origin of the observed pulsar proper motions. We argue that this effect was vastly overestimated because the variation of the temperature over the deformed neutrinosphere is not an adequate measure for the anisotropy of neutrino emission. The neutrino flux is generated inside the neutron star core and is transported through the atmosphere at a constant luminosity, forcing the temperature gradient in the atmosphere to adjust to the inflow of energy from below. Therefore, no emission anisotropy is caused by a deformation of the neutrinosphere to lowest order. An estimate of the higher-order corrections must take into account the modified atmospheric temperature profile in response to the deformation of the neutrinosphere and the corresponding feedback on the core. We go through this exercise in the framework of a simplified model which can be solved analytically.Comment: Final version with minor corrections, to be published in PRD. Includes a "Note Added" in response to astro-ph/981114

Dynamics of Primordial Black Hole Formation

We present a numerical investigation of the gravitational collapse of horizon-size density fluctuations to primordial black holes (PBHs) during the radiation-dominated phase of the Early Universe. The collapse dynamics of three different families of initial perturbation shapes, imposed at the time of horizon crossing, is computed. The perturbation threshold for black hole formation, needed for estimations of the cosmological PBH mass function, is found to be $\delta_{\rm c} \approx 0.7$ rather than the generally employed $\delta_{\rm c} \approx 1/3$, if $\delta$ is defined as \Delta M/\mh, the relative excess mass within the initial horizon volume. In order to study the accretion onto the newly formed black holes, we use a numerical scheme that allows us to follow the evolution for long times after formation of the event horizon. In general, small black holes (compared to the horizon mass at the onset of the collapse) give rise to a fluid bounce that effectively shuts off accretion onto the black hole, while large ones do not. In both cases, the growth of the black hole mass owing to accretion is insignificant. Furthermore, the scaling of black hole mass with distance from the formation threshold, known to occur in near-critical gravitational collapse, is demonstrated to apply to primordial black hole formation.Comment: 10 pages, 8 figures, revtex style, submitted to PR

Three-dimensional simulations of type Ia supernovae

We present the results of three-dimensional hydrodynamical simulations of the subsonic thermonuclear burning phase in type Ia supernovae. The burning front model contains no adjustable parameters so that variations of the explosion outcome can be linked directly to changes in the initial conditions. In particular, we investigate the influence of the initial flame geometry on the explosion energy and find that it appears to be weaker than in 2D. Most importantly, our models predict global properties such as the produced nickel masses and ejecta velocities within their observed ranges without any fine tuning.Comment: 7 pages, 5 figures, accepted by A&

The Correlation Function of Clusters of Galaxies and the Amplitude of Mass Fluctuations in the Universe

We show that if a sample of galaxy clusters is complete above some mass threshold, then hierarchical clustering theories for structure formation predict its autocorrelation function to be determined purely by the cluster abundance and by the spectrum of linear density fluctuations. Thus if the shape of the initial fluctuation spectrum is known, its amplitude $\sigma_8$ can be estimated directly from the correlation length of a cluster sample in a way which is independent of the value of $\Omega_0$. If the cluster mass corresponding to the sample threshold is also known, it provides an independent estimate of the quantity $\sigma_8\Omega_0^{0.6}$. Thus cluster data should allow both $\sigma_8$ and $\Omega_0$ to be determined observationally. We explore these questions using N-body simulations together with a simple but accurate analytical model based on extensions of Press-Schechter theory. Applying our results to currently available data we find that if the linear fluctuation spectrum has a shape similar to that suggested by the APM galaxy survey, then a correlation length $r_0$ in excess of 20\mpch for Abell clusters would require $\sigma_8>1$, while r_0<15\mpch would require $\sigma_8<0.5$. With conventional estimates of the relevant mass threshold these imply \Omega_0\la 0.3 and \Omega_0\ga 1 respectively.Comment: Latex, 25 pages (including 8 PS figures). The PS-file of the paper is also available via anonymous ftp at: ftp://ibm-3.mpa-garching.mpg.de/pub/jing/xicc.ps . Submitted to MNRAS. In the replaced version, a typo in Eq.(1a) is fixe

On the Microlensing Optical Depth of the Galactic Bar

The microlensing probability (optical depth $\tau$) toward the Galactic center carries information about the mass distribution of the Galactic bulge/bar, so can be used to constrain the very uncertain shape parameters of the bar. We find $tau$ depends on the bar mass, radial profile, angle, axis scale lengths and boxyness by a few simple analytical formulae, which shows: (1) $\tau$ is proportional to the mass of the bar, $M$. (2) $\tau$ falls along the minor axis with a strong gradient. (3) An oblate bulge can have more optical depth than a triaxial bar if the bar angle $\alpha>45$ degress. (4) $\tau$ is the largest if the angle $\alpha$ and the axis ratio $y_0/x_0$ conspires so that $y_0/x_0=\tan \alpha$. (5) At a fixed field on the minor axis but away from the center, boxy bars with a flat density profile tend to give a larger optical depth than ellipsoidal bars with a steep profile. (6) Main sequence sources should have a significantly lower (20-50\% lower) optical depth than red clump giants if main sequence stars are not observed as deep as the bright clump giants. An application to four COBE-constrained models (Dwek et al. 1994) shows most models produce optical depth $2\sigma$ lower than MACHO and OGLE observed values even with both a massive bar $2.8\times 10^{10} M_\odot$ and a full disk. The high $\tau$ argues for a massive ($> 2\times 10^{10}M_\odot$) boxy bar with $y_0/x0\approx \tan\alpha$ and $\alpha<20$ deg and with a flat radial profile up to corotation.Comment: 28 pages including 6 postscript figures in uuencoded compressed tar file. Submitted to MNRA

Microwave polarization in the direction of galaxy clusters induced by the CMB quadrupole anisotropy

Electron scattering induces a polarization in the cosmic microwave background (CMB) signal measured in the direction of a galaxy cluster due to the presence of a quadrupole component in the CMB temperature distribution. Measuring the polarization towards distant clusters provides the unique opportunity to observe the evolution of the CMB quadrupole at moderate redshifts, z~0.5-3. We demonstrate that for the local cluster population the polarization degree will depend on the cluster celestial position. There are two extended regions in the sky, which are opposite to each other, where the polarization is maximal, 0.1(tau/0.02) microK in the Rayleigh-Jeans part of the CMB spectrum (tau being the Thomson optical depth across the cluster) exceeding the contribution from the cluster transverse peculiar motion if v_t<1300 km/s. One can hope to detect this small signal by measuring a large number of clusters, thereby effectively removing the systematic contribution from other polarization components produced in clusters. These polarization effects, which are of the order of (v_t/c)^2 tau, (v_t/c) tau^2 and (kT_e/m_ec^2) tau^2, as well as the polarization due to the CMB quadrupole, were previously calculated by Sunyaev and Zel'dovich for the Rayleigh-Jeans region. We fully confirm their earlier results and present exact frequency dependencies for all these effects. The polarization is considerably higher in the Wien region of the CMB spectrum.Comment: 8 pages, 5 figures, submitted to MNRA

Explosive Nucleosynthesis in GRB Jets Accompanied by Hypernovae

Two-dimensional hydrodynamic simulations are performed to investigate explosive nucleosynthesis in a collapsar using the model of MacFadyen and Woosley (1999). It is shown that 56Ni is not produced in the jet of the collapsar sufficiently to explain the observed amount of a hypernova when the duration of the explosion is \sim 10 sec, which is considered to be the typical timescale of explosion in the collapsar model. Even though a considerable amount of 56Ni is synthesized if all explosion energy is deposited initially, the opening angles of the jets become too wide to realize highly relativistic outflows and gamma-ray bursts in such a case. From these results, it is concluded that the origin of 56Ni in hypernovae associated with GRBs is not the explosive nucleosynthesis in the jet. We consider that the idea that the origin is the explosive nucleosynthesis in the accretion disk is more promising. We also show that the explosion becomes bi-polar naturally due to the effect of the deformed progenitor. This fact suggests that the 56Ni synthesized in the accretion disk and conveyed as outflows are blown along to the rotation axis, which will explain the line features of SN 1998bw and double peaked line features of SN 2003jd. Some fraction of the gamma-ray lines from 56Ni decays in the jet will appear without losing their energies because the jet becomes optically thin before a considerable amount of 56Ni decays as long as the jet is a relativistic flow. We show that abundance of nuclei whose mass number \sim 40 in the ejecta depends sensitively on the energy deposition rate. So it may be determined by observations of chemical composition in metal poor stars which model is the proper one as a model of a gamma-ray burst accompanied by a hypernova.Comment: 29 pages with 16 figures. ApJ, accepte

The distribution of supermassive black holes in the nuclei of nearby galaxies

The growth of supermassive black holes by merging and accretion in hierarchical models of galaxy formation is studied by means of Monte Carlo simulations. A tight linear relation between masses of black holes and masses of bulges arises if if the mass accreted by supermassive black holes scales linearly with the mass forming stars and if the redshift evolution of mass accretion tracks closely that of star formation. Differences in redshift evolution between black hole accretion and star formation introduce considerable scatter in this relation. A non-linear relation between black hole accretion and star formation results in a non-linear relation between masses of remnant black holes and masses of bulges. The relation of black hole mass to bulge luminosity obseved in nearby galaxies and its scatter are reproduced reasonably well by models in which black hole accretion and star formation are linearly related but do not track each other in redshift. This suggests that a common mechanism determines the efficiency for black hole accretion and the efficiency for star formation, especially for bright bulges.Comment: 6 pages, 3 figures, submitted to MNRA

Diversity of chemistry and excitation conditions in the high-mass star forming complex W33

The object W33 is a giant molecular cloud that contains star forming regions at various evolutionary stages from quiescent clumps to developed H II regions. Since its star forming regions are located at the same distance and the primary material of the birth clouds is probably similar, we conducted a comparative chemical study to trace the chemical footprint of the different phases of evolution. We observed six clumps in W33 with the Atacama Pathfinder Experiment (APEX) telescope at 280 GHz and the Submillimeter Array (SMA) at 230 GHz. We detected 27 transitions of 10 different molecules in the APEX data and 52 transitions of 16 different molecules in the SMA data. The chemistry on scales larger than $\sim$0.2 pc, which are traced by the APEX data, becomes more complex and diverse the more evolved the star forming region is. On smaller scales traced by the SMA data, the chemical complexity and diversity increase up to the hot core stage. In the H II region phase, the SMA spectra resemble the spectra of the protostellar phase. Either these more complex molecules are destroyed or their emission is not compact enough to be detected with the SMA. Synthetic spectra modelling of the H$_{2}$CO transitions, as detected with the APEX telescope, shows that both a warm and a cold component are needed to obtain a good fit to the emission for all sources except for W33 Main1. The temperatures and column densities of the two components increase during the evolution of the star forming regions. The integrated intensity ratios N$_{2}$H$^{+}$(3$-$2)/CS(6$-$5) and N$_{2}$H$^{+}$(3$-$2)/H$_{2}$CO(4$_{2,2}$$-$3$_{2,1}$) show clear trends as a function of evolutionary stage, luminosity, luminosity-to-mass ratio, and H$_{2}$ peak column density of the clumps and might be usable as chemical clocks.Comment: 66 pages, 28 figures, 8 tables, accepted for publication at A&