Nonadiabatic charged spherical evolution in the postquasistatic approximation

We apply the postquasistatic approximation, an iterative method for the evolution of self-gravitating spheres of matter, to study the evolution of dissipative and electrically charged distributions in General Relativity. We evolve nonadiabatic distributions assuming an equation of state that accounts for the anisotropy induced by the electric charge. Dissipation is described by streaming out or diffusion approximations. We match the interior solution, in noncomoving coordinates, with the Vaidya-Reissner-Nordstr\"om exterior solution. Two models are considered: i) a Schwarzschild-like shell in the diffusion limit; ii) a Schwarzschild-like interior in the free streaming limit. These toy models tell us something about the nature of the dissipative and electrically charged collapse. Diffusion stabilizes the gravitational collapse producing a spherical shell whose contraction is halted in a short characteristic hydrodynamic time. The streaming out radiation provides a more efficient mechanism for emission of energy, redistributing the electric charge on the whole sphere, while the distribution collapses indefinitely with a longer hydrodynamic time scale.Comment: 11 pages, 16 Figures. Accepted for publication in Phys Rev

Disks Surviving the Radiation Pressure of Radio Pulsars

The radiation pressure of a radio pulsar does not necessarily disrupt a surrounding disk. The position of the inner radius of a thin disk around a neutron star can be estimated by comparing the electromagnetic energy density generated by the neutron star with the kinetic energy density of the disk. Inside the light cylinder, the near zone electromagnetic field is essentially the dipole magnetic field, and the inner radius is the conventional Alfven radius. Far outside the light cylinder, in the radiation zone, $E=B$ and the electromagnetic energy density is $/c \propto 1/r^2$ where $S$ is the Poynting vector. Shvartsman (1970) argued that a stable equilibrium can not be found in the radiative zone because the electromagnetic energy density dominates over the kinetic energy density, with the relative strength of the electromagnetic stresses increasing with radius. In order to check whether this is true also near the light cylinder, we employ global electromagnetic field solutions for rotating oblique magnetic dipoles (Deutsch 1955). Near the light cylinder the electromagnetic energy density increases steeply enough with decreasing $r$ to balance the kinetic energy density at a stable equilibrium. The transition from the near zone to the radiation zone is broad. The radiation pressure of the pulsar can not disrupt the disk for values of the inner radius up to about twice the light cylinder radius if the rotation axis and the magnetic axis are orthogonal. This allowed range beyond the light cylinder extends much further for small inclination angles. We discuss implications of this result for accretion driven millisecond pulsars and young neutron stars with fallback disks.Comment: Accepted by Astrophysical Journal, final version with a minor correctio

The Power of General Relativity

We study the cosmological and weak-field properties of theories of gravity derived by extending general relativity by means of a Lagrangian proportional to $R^{1+\delta}$. This scale-free extension reduces to general relativity when $\delta \to 0$. In order to constrain generalisations of general relativity of this power class we analyse the behaviour of the perfect-fluid Friedmann universes and isolate the physically relevant models of zero curvature. A stable matter-dominated period of evolution requires $\delta >0$ or $\delta <-1/4$. The stable attractors of the evolution are found. By considering the synthesis of light elements (helium-4, deuterium and lithium-7) we obtain the bound $-0.017<\delta <0.0012.$ We evaluate the effect on the power spectrum of clustering via the shift in the epoch of matter-radiation equality. The horizon size at matter--radiation equality will be shifted by $\sim 1$ for a value of $\delta \sim 0.0005.$ We study the stable extensions of the Schwarzschild solution in these theories and calculate the timelike and null geodesics. No significant bounds arise from null geodesic effects but the perihelion precession observations lead to the strong bound $\delta =2.7\pm 4.5\times 10^{-19}$ assuming that Mercury follows a timelike geodesic. The combination of these observational constraints leads to the overall bound $0\leq \delta <7.2\times 10^{-19}$ on theories of this type.Comment: 26 pages and 5 figures. Published versio

Stellar-Mass Black Holes in the Solar Neighborhood

We search for nearby, isolated, accreting, ``stellar-mass'' (3 to $100M_\odot$) black holes. Models suggest a synchrotron spectrum in visible wavelengths and some emission in X-ray wavelengths. Of 3.7 million objects in the Sloan Digital Sky Survey Early Data Release, about 150,000 objects have colors and properties consistent with such a spectrum, and 87 of these objects are X-ray sources from the ROSAT All Sky Survey. Thirty-two of these have been confirmed not to be black-holes using optical spectra. We give the positions and colors of these 55 black-hole candidates, and quantitatively rank them on their likelihood to be black holes. We discuss uncertainties the expected number of sources, and the contribution of blackholes to local dark matter.Comment: Replaced with version accepted by ApJ. 40 pages, 8 figure

A variational approach to the stochastic aspects of cellular signal transduction

Cellular signaling networks have evolved to cope with intrinsic fluctuations, coming from the small numbers of constituents, and the environmental noise. Stochastic chemical kinetics equations govern the way biochemical networks process noisy signals. The essential difficulty associated with the master equation approach to solving the stochastic chemical kinetics problem is the enormous number of ordinary differential equations involved. In this work, we show how to achieve tremendous reduction in the dimensionality of specific reaction cascade dynamics by solving variationally an equivalent quantum field theoretic formulation of stochastic chemical kinetics. The present formulation avoids cumbersome commutator computations in the derivation of evolution equations, making more transparent the physical significance of the variational method. We propose novel time-dependent basis functions which work well over a wide range of rate parameters. We apply the new basis functions to describe stochastic signaling in several enzymatic cascades and compare the results so obtained with those from alternative solution techniques. The variational ansatz gives probability distributions that agree well with the exact ones, even when fluctuations are large and discreteness and nonlinearity are important. A numerical implementation of our technique is many orders of magnitude more efficient computationally compared with the traditional Monte Carlo simulation algorithms or the Langevin simulations.Comment: 15 pages, 11 figure

Three-dimensional MHD Simulations of Radiatively Inefficient Accretion Flows

We present three-dimensional MHD simulations of rotating radiatively inefficient accretion flows onto black holes. In the simulations, we continuously inject magnetized matter into the computational domain near the outer boundary, and we run the calculations long enough for the resulting accretion flow to reach a quasi-steady state. We have studied two limiting cases for the geometry of the injected magnetic field: pure toroidal field and pure poloidal field. In the case of toroidal field injection, the accreting matter forms a nearly axisymmetric, geometrically-thick, turbulent accretion disk. The disk resembles in many respects the convection-dominated accretion flows found in previous numerical and analytical investigations of viscous hydrodynamic flows. Models with poloidal field injection evolve through two distinct phases. In an initial transient phase, the flow forms a relatively flattened, quasi-Keplerian disk with a hot corona and a bipolar outflow. However, when the flow later achieves steady state, it changes in character completely. The magnetized accreting gas becomes two-phase, with most of the volume being dominated by a strong dipolar magnetic field from which a thermal low-density wind flows out. Accretion occurs mainly via narrow slowly-rotating radial streams which `diffuse' through the magnetic field with the help of magnetic reconnection events.Comment: 35 pages including 3 built-in plots and 14 separate jpg-plots; version accepted by Ap

Relic neutrino decoupling including flavour oscillations

In the early universe, neutrinos are slightly coupled when electron-positron pairs annihilate transferring their entropy to photons. This process originates non-thermal distortions on the neutrino spectra which depend on neutrino flavour, larger for nu_e than for nu_mu or nu_tau. We study the effect of three-neutrino flavour oscillations on the process of neutrino decoupling by solving the momentum-dependent kinetic equations for the neutrino spectra. We find that oscillations do not essentially modify the total change in the neutrino energy density, giving N_eff=3.046 in terms of the effective number of neutrinos, while the small effect over the production of primordial 4He is increased by O(20%), up to 2.1 x 10^{-4}. These results are stable within the presently favoured region of neutrino mixing parameters.Comment: 18 pages, 2 figure

On the structure of the burst and afterglow of Gamma-Ray Bursts I: the radial approximation

We have proposed three paradigms for the theoretical interpretation of gamma-ray bursts (GRBs). (1) The relative space-time transformation (RSTT) paradigm emphasizes how the knowledge of the entire world-line of the source from the moment of gravitational collapse is a necessary condition to interpret GRB data. (2) The interpretation of the burst structure (IBS) paradigm differentiates in all GRBs between an injector phase and a beam-target phase. (3) The GRB-supernova time sequence (GSTS) paradigm introduces the concept of induced supernova explosion in the supernovae-GRB association. These three paradigms are illustrated using our theory based on the vacuum polarization process occurring around an electromagnetic black hole (EMBH theory) and using GRB 991216 as a prototype. We illustrate the five fundamental eras of the EMBH theory: the self acceleration of the $e^+e^-$ pair-electromagnetic plasma (PEM pulse), its interaction with the baryonic remnant of the progenitor star (PEMB pulse). We then study the approach of the PEMB pulse to transparency, the emission of the proper GRB (P-GRB) and its relation to the ``short GRBs''. Finally the three different regimes of the afterglow are described within the fully radiative and radial approximations. The best fit of the theory leads to an unequivocal identification of the ``long GRBs'' as extended emission occurring at the afterglow peak (E-APE). The relative intensities, the time separation and the hardness ratio of the P-GRB and the E-APE are used as distinctive observational test of the EMBH theory and the excellent agreement between our theoretical predictions and the observations are documented. The afterglow power-law indexes in the EMBH theory are compared and contrasted with the ones in the literature, and no beaming process is found for GRB 991216.Comment: 96 pages, 40 figures, to appear on Int. Journ. Mod. Phys.

Effect of Neutrino Heating on Primordial Nucleosynthesis

We have modified the standard code for primordial nucleosynthesis to include the effect of the slight heating of neutrinos by $e^\pm$ annihilations. There is a small, systematic change in the $^4$He yield, $\Delta Y \simeq +1.5\times 10^{-4}$, which is insensitive to the value of the baryon-to-photon ratio $\eta$ for 10^{-10}\la \eta \la 10^{-9}. We also find that the baryon-to-photon ratio decreases by about 0.5\% less than the canonical factor of 4/11 because some of the entropy in $e^\pm$ pairs is transferred to neutrinos. These results are in accord with recent analytical estimates.Comment: 14 pages/4 Figs (upon request