Primordial black holes, formed from rare peaks in the primordial fluctuations $\zeta$, are non-rotating at zeroth order in $\zeta_{\rm rms}$. We show that the spin also vanishes at first order in $\zeta_{\rm rms}$, suggesting the dimensionless spin parameter $a_{\rm rms} \sim \zeta_{\rm rms}^2$. We identify one quadratic contribution to the spin by calculating (and extrapolating to the formation time) the torque on a black hole due to ambient acoustic waves. For a reasonable density of primordial black holes this implies a percent level spin parameter.

Nonlinear radio waves modulate the plasma, scatter on the modulations, and develop an intermittent power spectrum -- perhaps. Rudiments of theory, numerical simulations, and qualitative modeling of nonlinear scattering are presented. Nonlinear radio waves in electron-positron plasma are shown to be unstable, the instability growth rate is calculated. The long-term outcome of the instability is studied numerically; an intermittent power spectrum is found. The frequency coherence scale of the intermittent power spectrum is interpreted as a measure of the nonlinear scattering screen.

An over-dense wake is created by a gravitating object moving through a gaseous medium, and this wake pulls back on the object and slows it down. This is conventional dynamical friction in a gaseous medium. We argue that if the object drives a sufficiently powerful outflow, the wake is destroyed and instead an extended under-dense region is created behind the object. In this case the overall gravitational force is applied in the direction of the object's motion, producing a negative dynamical friction (NDF). Black holes in dense gas drive powerful outflows and may experience the NDF, although extensive numerical work is probably needed to demonstrate or refute this conclusively. NDF may be important for stellar-mass black holes and neutron stars inside `common envelopes' in binary systems, for stellar mass black holes inside AGN discs, or for massive black holes growing through super-Eddington accretion in early Universe.

Faraday Rotation-Conversion is the simultaneous rotation of all three Stokes polarization parameters $Q$, $U$, $V$ as an electromagnetic wave propagates through a magnetized plasma. In this regime the Faraday plasma screen is characterized by more than just a Rotation Measure. We define the Conversion Measure that characterizes the wavelength-dependent conversion between the linear and circular polarization. In a cold plasma, the conversion occurs at the localized regions along the wave's path, where the large-scale magnetic field is perpendicular to the propagation direction. The number of these regions along the line of sight through the screen, and their individual contributions to the conversion measure, can be inferred from the polarization measurements. We argue that the simultaneous measurement of wavelength-dependent linear and circular polarization might give an important insight into the magnetic-field geometry of the Faraday screen in FRB121102 and other repeating Fast Radio Bursts.

A sufficient condition for maser instability in a weakly magnetized relativistic plasma with an isotropic particle distribution function is given. The maser growth rates and polarizations are computed starting from the exact dielectric permittivity tensor of a magnetized plasma. For very weak magnetic fields, our results confirm the approximate validity of the 'standard maser theory', which is based on the Einstein coefficients method, with one significant exception. For inclined propagation and realistic (small but finite) field, the growth rates of the two (nearly circular) polarizations differ significantly, while the standard theory predicts two (nearly circular) polarizations with similar growth rates. We show that this deviation is due to circularly polarized synchrotron emission, which is neglected in the standard theory. The maser is shown to grow slower than Langmuir waves. Nevertheless, significant generation of EM waves is seen in (highly simplified) direct numerical simulations. We study the nonlinear saturation of the maser instability and find that it offers a mechanism for the conversion of a significant fraction of the plasma energy into radio waves. We briefly discuss the conditions under which the maser instability may operate in astrophysical sources, and provide rough estimates that may be used as a guidance when studying particular astrophysical sources/phenomena.

Cosmic strings and primordial black holes (PBHs) commonly and naturally form in many scenarios describing the early universe. Here we show that if both cosmic strings and PBHs are present, their interaction leads to a range of interesting consequences. At the time of their formation, the PBHs get attached to the strings and influence their evolution, leading to the formation of black-hole-string networks and commonly to the suppression of loop production in a range of redshifts. Subsequently, reconnections within the network give rise to small nets made of several black holes and connecting strings. The number of black holes in the network as well as the stability of the nets depend on the topological properties of the strings. The nets oscillate and shrink exponentially due to the emission of gravitational waves. This leads to potentially observable string-driven mergers of PBHs. The strings can keep PBHs from galactic halos, making the current bounds on PBHs not generally applicable. Alternatively, heavy PBHs can drag low-tension strings into the centers of galaxies. The superconducting strings can appear as radio filaments pointing towards supermassive black holes.

Auger Collaboration has reported a large-scale anisotropy in the arrival directions of cosmic rays above 8 EeV. The dipole direction, at 125$^\circ$ from the Galactic center, is taken as an indication of an extragalactic origin of these cosmic rays. We show, both theoretically and by direct numerical simulations, that this is not necessarily true. Intermediate mass nuclei originating in the Galaxy and diffusing in the Galactic magnetic field can have a dipole anisotropy pointing away from the Galactic center. Our theory predicts a quadrupole anisotropy of the same order of magnitude as the dipole.

Thermalized fireballs should be created by cusp events on superconducting cosmic strings. This simple notion allows to reliably estimate particle emission from the cusps in a given background magnetic field. With plausible assumptions about intergalactic magnetic fields, the cusp events can produce observable fluxes of high-energy photons and neutrinos with unique signatures.

The saturation mechanism for the fastest-growing instability of massive scalar field in Kerr metric is identified, assuming saturation by cubic or quartic nonlinearities of the field potential. The resulting spindown rate of the black hole is calculated. The (rather involved) saturation scenario is confirmed by numerical simulations.

It is proposed that there exists a class of pulsars, called weak pulsars, for which the large-scale magnetosphere, and hence the gamma-ray emission, are independent of the detailed pattern of plasma production. The weak pulsar magnetosphere and its gamma-ray emission are uniquely determined by just three parameters: spin, dipole, and the spin-dipole angle. We calculate this supposedly unique pulsar magnetosphere in the axisymmetric case. The magnetosphere is found to be very close to (although interestingly not fully identical with) the magnetosphere we have previously calculated, explaining the phenomenological success of the old calculation. We offer only a highly tentative proof of this "Pulsar No-Hair Theorem". Our analytics, while convincing in its non-triviality, is incomplete, and counts only as a plausibility argument. Our numerics, while complete, is dubious. The plasma flow in the weak pulsar magnetosphere turns out to be even more intricate than what we have previously proposed: some particles, after being created near the star, move beyond the light cylinder and then return to the star.

The mechanism by which gamma-ray pulsars shine might be reproducible in a laboratory. This claim is supported by three observations: (i) properly focusing a few PW optical laser gives an electromagnetic field in the so-called Aristotelian regime, when a test electron is radiation-overdamped; (ii) the Goldreich-Julian number density of this electromagnetic field (the number density of elementary charges needed for a nearly full conversion of optical power into gamma-rays) is of order the electron number density in a solid; (iii) above about $50$PW, the external source of electrons is not needed -- charges will be created by a pair production avalanche.

Recent solution of a pulsar (although likely correct, as agreeing with phenomenology of weak pulsars without a single adjustable) is somewhat unsatisfactory on account of being heavily numerical. Here the main features of a pulsar -- (i) the existence of a force-free zone, (ii) bounded by a non-force-free radiation zone, (iii) where colliding Poynting fluxes annihilate into curvature radiation -- are reproduced as an exact solution of the same pulsar equations but in a much simpler geometry.

Assuming self-similarity of the first kind, we get three possible values p=1/2, 1, 3/2 for the exponent describing the density profile, rho r^{-p}, of a non-radiative (and hence quasi-spherical) accretion flow. The high and low p cases are known as Bondi and Convection-Dominated accretion flows. The 1-law flow we tentatively identify with the so-called Magnetically-Frustrated accretion flow. If our interpretation is correct, the accretion flow must be, roughly speaking, a collection of Prandtl's turbulent jets. The 1-law flow, being a first-kind self-similar solution, may actually occur in nature (in collisionless plasma).

Nuclei, if lifted into the radiation zone of a weak axisymmetric pulsar, are accelerated to quasi-monoenergetic spectrum, with more than half of the nuclear luminosity emitted nearly monoenergetically, at the energy E approximately = 4Ze(L/c)^(1/2), where L is the spin-down power. With a number of (not obviously wrong) assumptions, this turns Fermi (-like) weak pulsars into prominent accelerators at energies near the cosmic ray knee. The main uncertainty is purely theoretical: one needs to understand if nuclei are major positive charge carries in some weak pulsar magnetospheres.

To illustrate the standard candle property of gamma-ray pulsars (and also to thereby confirm the recent first-principle calculation of pulsar gamma-ray emission), we "measure", via the lightcurve fitting, three distances and one moment of inertia of some weak pulsars. We are not sure what the three distances are good for, but the measurement of the moment of inertia must be of interest for nuclear physics. Although we must state that the quality of the numerical program which we use to calculate the lightcurves and efficiencies is inadequate (as are the author's qualifications as a numericist and data analyst), in good hands, and upon an easily doable extension to non-weak pulsars, the method's yield should be impressive.

Measurement of distances and moments of inertia of pulsars must be useful, for instance, for three-dimensional mapping of the dispersion and rotation measures, constraining the nuclear equation of state, etc. The distances and moments of inertia can be measured by fitting the gamma-ray lightcurves of pulsars, because the gamma-ray emission seems to be governed by easily calculable physics. The (first-principle) theoretical lightcurves have been computed only for weak pulsars (pair production near the light cylinder much smaller than Goldreich-Julian [1] per rotation), and at insufficient accuracy; but, since this computation has been done by a self-taught numericist, it must be possible to improve the theoretical accuracy for weak pulsars, and also to extend the computation to non-weak pulsars. To invite the computational effort of better-equipped researchers, we describe an (entirely obvious) procedure for measuring the distances and moments of inertia for weak pulsars.

Energy-resolved lightcurves are calculated for a weak pulsar.

Emission spectrum is calculated for a weak axisymmetric pulsar. Also calculated are the observed spectrum, efficiency, and the observed efficiency. The underlying flow of electrons and positrons turns out to be curiously intricate.

For a given solar magnetic field, the near-Sun (phase-space) density of cosmic ray electrons and positrons of energy above about 10GeV can be calculated from first principles, without any assumptions about the cosmic ray diffusion. This is because the sunlight Compton drag must be more important than diffusion. If the solar magnetic field has an appreciable dipole component, the electron/positron density should have a belt-like dent, perhaps extending to several solar radii. The belt structure appears because the quasi-bound orbits are depopulated by the sunlight Compton drag.

We define and calculate the "discharge mode" for a Schwarzschild black hole in massive electrodynamics. For small photon mass, the discharge mode describes the decay of the electric field of a charged star collapsing into a black hole. We argue that a similar "discharge of mass" occurs in massive gravity and leads to a strange process of black hole disappearance.

In Aristotelian Electrodynamics (AE), due to radiation over-damping, the velocity rather than the acceleration of a charge is determined by the local electromagnetic field. Treating electrons and positrons separately, AE seems to give a faithful description of the flow of charges in a pulsar magnetosphere. AE might allow to calculate the pulsar emission in full detail, at all energies except the radio. But AE is not a full theory -- the plasma production processes must be added. Here we show that the pulsar magnetosphere and the resulting emission are sensitive to the plasma production rate near the light cylinder. In strong pulsars (high production rate, due to photon-photon collisions), the number of field lines closing beyond the light cylinder decreases. This reduces both the pulsar efficiency (the ratio of the pulsed bolometric luminosity to the spin-down power) and the characteristic photon energy, in overall agreement with the Fermi data. For weak pulsars (low production rate), our numerical implementation of AE works stably and seems to be ready to calculate the lightcurves and spectra. But for strong pulsars, the code (included) needs improvement or better understanding.

Pulsar efficiency, defined as the ratio of the pulsed bolometric luminosity to the spin-down power, is calculated to be about 15% (averaged over the spin-dipole inclination angle, ranging between about 50% for the aligned and 10% for the orthogonal). We also estimate the characteristic photon energy and argue that our results agree with the Fermi pulsar catalog -- in a sense.

Electromagnetic field together with zero-mass charges moving in this field form a well-behaved semi-dissipative dynamical system -- Electrodynamics of Massless Charges (EMC). We give equations of EMC, argue that EMC is an adequate theory for calculating pulsar magnetospheres, give an illustrative numerical calculation (showing that bolometric luminosity of an aligned rotator is approximately equal to half the spin-down power). EMC looks like a portion of the full pulsar theory that will resolve the already calculated bolometric luminosity into light curves and spectra.

We have recently proposed that Force-Free Electrodynamics (FFE) does not apply to pulsars -- pulsars should be described by the high-conductivity limit of Strong-Field Electrodynamics (SFE), which predicts an order-unity damping of the Poynting flux, while FFE postulates zero damping. The strong damping result has not been accepted by several pulsar experts, who claim that FFE basically works and the Poynting flux damping can be arbitrarily small. Here we consider a thought experiment -- cylindrical periodic pulsar. We show that FFE is incapable of describing this object, while SFE predictions are physically plausible. The intrinsic breakdown of FFE should mean that the FFE description of the singular current layer (the only region of magnetosphere where FFE and the high-conductivity SFE differ) is incorrect. Then the high-conductivity SFE should be the right theory for real pulsars too, and the pure-FFE description of pulsars should be discarded.

Generically, massive gravity gives a non-unique gravitational field around a star. For a special family of massive gravity theories, we show that the stellar gravitational field is unique and observationally acceptable, that is close to Einsteinian. The black hole solutions in this family of theories are also studied and shown to be peculiar. Black holes have a near-horizon throat and the curvature diverging at the horizon. We show that there exists a sub-family of these massive gravity theories with non-singular at horizon black holes.

Ideal axisymmetric pulsar magnetosphere is calculated from the standard stationary force-free equation but with a new boundary condition at the equator. The new solution predicts Ohmic heating. About 50% of the Poynting power is dissipated in the equatorial current layer outside the light cylinder, with about 10% dissipated between 1 and 1.5 light cylinder radii. The Ohmic heat presumably goes into radiation, pair production, and acceleration of charges -- in an unknown proportion.

We argue that existing models of the ideal pulsar magnetosphere are incorrect because of the improper treatment of the singular current layer outside the light cylinder. We simulated the axisymmetric pulsar magnetosphere in the Force-Free limit of Strong-Field Electrodynamics. It turns out that even in the Force-Free limit, some field lines enter the singular current layer which lies beyond the light cylinder in the equatorial plane. As a result: (i) about 10% of the Poynting flux is dissipated between 1 and 1.5 light cylinder radii, (ii) there is no singular current layer within the light cylinder.

A formula is derived for the rate of thermal atmospheric escape, valid, and asymptotically exact, at low Knudsen number.

Stability and attractor property of free-floating axisymmetric magnetic bubbles in high-conductivity plasmas is (tentatively, numerically) demonstrated. The existence of compact non-singular axisymmetric magnetic equilibria is proved. Being attractors, the solitary magnetic bubbles should exist in nature.

We propose that generic magnetic equilibrium of an ideally conducting fluid contains a volume-filling set of singular current layers. Singular current layers should exist inside neutron stars. Residual dissipation in the singular current layers might be the main mechanism for the magnetic field decay. The slow decay of the field might be the clock responsible for triggering the magnetar flares.

We show that the first structures that form in the universe should spontaneously generate magnetic fields. No primordial seed field is required for this "first dynamo". Although the first dynamo starts with kinetic plasma instabilities, we argue that an adequate magnetohydrodynamic description might be possible via a simple trick. This should allow a numerical study of the effect of the first magnetic fields on the first baryonic objects.

Torsion oscillations of the neutron star crust are Landau damped by the Alfven continuum in the bulk. For strong magnetic fields (in magnetars), undamped Alfven eigenmodes appear.

Dissipative axisymmetric pulsar magnetosphere is calculated by a direct numerical simulation of the Strong-Field Electrodynamics equations. The magnetic separatrix disappears, it is replaced by a region of enhanced dissipation. With a better numerical scheme, one should be able to calculate the bolometric lightcurves for a given conductivity.

We predict a new class of very intense, short-duration AGN flares capable of accelerating the highest energy cosmic rays, resulting from the tidal disruption of a star or from a disk instability. The rate and power of these flares readily explains the observed flux and density statistics of UHECRs. The photon bursts produced by the predicted AGN flares are discussed; they may soon be detectable. Observations are shown to exclude that continuous jets of powerful Active Galactic Nuclei are the sole source of ultrahigh energy cosmic rays; the stringent requirements for Gamma Ray Bursts to be the source are delineated.

We describe a scenario for large-scale magnetic field generation and particle acceleration in a collisionless collision of cold plasma clouds. A first-principle (i.e. using particles) numerical simulation of this process might be possible. Our scenario is essentially 3D. We argue that {\it large-scale} magnetic fields are not generated in 2D, even in collisionless plasma. We calculate and numerically simulate magnetic field generation by relativistic collisionless Kelvin-Helmholtz instability in 2D. Collisionless tangential discontinuity might be more important than collisionless shock, because tangential discontinuity remains unstable even in the hydro limit, when the shock stabilizes.

Strong-Field Electrodynamics (SFE) is Maxwell theory with a certain Lorentz-covariant Ohm's law which uses only the electromagnetic degrees of freedom. We show that SFE is {\it semi-dissipative}: while the dissipation rate of the electromagnetic energy is non-negative, it can be exactly zero for non-trivial electromagnetic fields. It appears that SFE is well-defined for arbitrary electromagnetic fields. It should be possible to calculate the dissipative pulsar magnetosphere and resolve the magnetic separatrix using SFE. We show that SFE reduces to Force-Free Electrodynamics (FFE) in the large conductivity limit. In the regions where the ideal FFE 4-current is space-like, SFE predicts small dissipative corrections. In the regions where the ideal FFE 4-current is time-like, SFE predicts a zero correction. This indicates that bright pulsars radiate primarily from the magnetic separatrix.

We prove that static fluid stars can stably support magnetic fields (within the ideal MHD approximation).

A dissipative Lorentz-covariant Ohm's law which uses only the electromagnetic degrees of freedom is proposed. For large conductivity, Maxwell equations equipped with this Ohm's law reduce to the equations of Force-Free Electrodynamics (FFE) with small dissipative corrections, but only in the regions where the ideal FFE 4-current is space-like. This might indicate that the pulsar emission comes primarily from the magnetic separartrix.

This paper suffers from conceptual difficulties and unjustified approximations that render its conclusions dubious.

A spherically symmetrical ultra-relativistic blast wave is not an attractor of a generic asymmetric explosion. Spherical symmetry is reached only by the time the blast wave slows down to non-relativistic velocities, when the Sedov-Taylor-von Neumann attractor solution sets in. We show however, that a directed relativistic explosion, with the explosion momentum close to the explosion energy, produces a blast wave with a universal intermediate asymptotic -- a selfsimilar directed ultra-relativistic blast wave. This universality might be of interest for the astrophysics of gamma-ray burst afterglows.

Pulsar emission should primarily come from the magnetic separatrix... Combining theory and observations, we show that force-free electrodynamics (FFE) gives an accurate description of the large-scale electromagnetic field in the magnetospheres of Crab-like pulsars. A robust prediction of FFE is the existence and stability of a singular current layer on the magnetic separatrix. We argue that most of the observed pulsar emission comes from this singular current layer.

We show that causality constrains the sign of quartic Riemann corrections to the Einstein-Hilbert action. Our constraint constitutes a restriction on candidate theories of quantum gravity.

We argue that it is currently impossible to simulate X-ray clusters using correct equations, because even the MHD description is not applicable. But since fluid simulations actually reproduce observations quite well, one may try to improve the fluid codes by including molecular transport of heat and momentum. We calculate the effective molecular viscosity for the simplest model of magnetic field and obtain 1/5 of the Braginskii value, similar to 1/3 of Spitzer for the heat conduction. This is large enough to noticeably damp the X-ray cluster turbulence.

The non-linearity of Einstein's equations makes it possible for small-scale matter inhomogeneities to affect the Universe at cosmological distances. We study the size of such effects using a simple heuristic model that captures the most important backreaction effect due to nonrelativistc matter, as well as several exact solutions describing inhomogeneous and anisotropic expanding universes. We find that the effects are $O(H^2l^2/c^2)$ or smaller, where $H$ is the Hubble parameter and $l$ the typical size scale of inhomogeneities. For virialized structures this is of order $v^2/c^2$, where $v$ is the characteristic peculiar velocity.

Convection that develops behind the shock front during the first second of a core-collapse supernova explosion is believed to play a crucial role in the explosion mechanism. We demonstrate that the resulting turbulent density fluctuations may be directly observable in the neutrino signal starting at t>~ 3-4 s after the onset of the explosion. The effect comes from the modulation of the MSW flavor transformations by the turbulent density fluctuations. We derive a simple and general criterion for neutrino flavor depolarization in a Kolmogorov-type turbulence and apply it to the turbulence seen in modern numerical simulations. The turbulence casts a "shadow", by making other features, such as the shock front, unobservable in the density range covered by the turbulence.

We show that FFE (Force-Free Electrodynamics) describes pulsar magnetospheres -- even when radiation and particle wind power are comparable to the spin-down power. At the same time, we show that FFE is insufficient for calculating pulsar magnetospheres. This is because FFE admits a large family of stationary solutions with different currents flowing in the closed line region. To choose the actual solution one needs a model of current generation which goes beyond pure FFE. We calculate several FFE magnetospheres for the aligned rotator. The Poynting power of these solutions fills the range $c^{-3}\mu ^2\Omega^4~< ~L~ <~ 0.67 c^{-3}\mu ^2\Omega^4(c/\Omega R_s)^2$, for angular velocity $\Omega$ and magnetic field which is a pure dipole $\mu$ on the surface of the star of radius $R_s$. Anomalous braking indices of young pulsars might be explained by currents flowing in the closed-line region.

We formulate variation principle for force-free magnetosphere of an inclined pulsar: ${\cal E} +{\bf \Omega}\cdot {\bf M}$ ($\cal E$, ${\bf M}$ are electromagnetic energy and angular momentum, ${\bf \Omega}$ is the angular velocity of a star) is stationary under isotopological variations of magnetic field and arbitrary variations of electric field. The variation principle gives the reason for existence and proves local stability of current singular layers along magnetic separatrices. Magnetic field lines of inclined pulsar magnetosphere lie on magnetic surfaces, and do have magnetic separatrices. In the framework of the isotopological variation principle, inclined magnetospheres are expected to be simple deformations of the axisymmetric pulsar magnetosphere. A singular line should exist on the light cylinder, where inner separatrix terminates and outer separatrix emanates. The electromagnetic field should have an inverse square root singularity near the singular line inside the inner magnetic separatrix. Large distance asymptotic solution is calculated, and used to estimate the pulsar power, $L\approx c^{-3}\mu ^2\Omega^4$ for spin-dipole inclinations $\lesssim 30^{\circ}$ .

We speculate that stationary force-free pulsar magnetospheres are screw unstable, and the spin-down power of real pulsars is carried by turbulent Poynting jets. The turbulent jet entrains poloidal flux as it propagates away from the star. Due to poloidal flux entrainment, poloidal field inside the jet decreases at a slower rate than the dipole's $r^{-3}$, and the pulsar power increases accordingly: $L\sim (\mu ^2\Omega^4/c^3)(\Omega R/c)^{-\alpha}$, where $\mu$ is the magnetic dipole, $\Omega$ is the frequency, $R$ is the neutron star radius, and $\alpha \lesssim 1$ is an index which depends on spin-dipole angle. Our speculation is of interest because it seems to provide the only possible explanation (without fine-tuning) of observed pulsar braking indices.

Stationary force-free magnetosphere of an axisymmetric pulsar is shown to have a separatrix inclination angle of 77.3$^\circ$. The electromagnetic field has an $R^{-1/2}$ singularity inside the separatrix near the light cylinder. A numerical simulation of the magnetosphere which crudely reproduces these properties is presented. The numerical results are used to estimate the power of an axisymmetric pulsar: $L=(1\pm 0.1)\mu^2\Omega^4/c^3$. A need for a better numerical simulation is pointed out.

Single scalar field inflation with a generic, non-quadratic in derivatives, field Lagrangian is considered. It is shown that non-Gaussianity of curvature perturbations is characterized by two dimensionless amplitudes. One of these amplitudes can be expressed in terms of the usual inflationary observables -- the scalar power, the tensor power, and the tensor index. This consistency relation provides an observational test for the single scalar inflation.

Acceleration of the universe might be driven by a continuous elastic medium -- elastic dark energy (Bucher and Spergel 1999). Elastic dark energy can stably support equations of state with pressure to energy ratio w > -4/3. Stable expansion with w<-1 combined with an assumption of ``reheating of elastic energy'' leads to exotic possibilities such as Expanding Cyclic Universe -- an ever-expanding universe with periodically repeating inflationary epochs.

Inflation of a universe filled by an elastic continuous medium is considered. Elastic inflation, while capable of describing current observations, is qualitatively different from standard inflationary models. The scalar and tensor modes keep evolving after crossing the horizon. Due to this superhorizon evolution, the amplitude of the tensor mode (inflationary gravitational wave) is not simply proportional to the energy density of the universe at the time of horizon crossing. The spectral index of tensor modes can be positive.

It is noted that inflation driven by a three form with arbitrary self coupling satisfies a consistency relation $n_S-1=2n_T$, between spectral indices of scalar and tensor cosmological perturbations. The standard consistency relation for the ratio of the amplitudes of perturbations is only valid for quadratic self coupling of the three form.

Curvature perturbations generated after the end of inflation are estimated for a generic multi-field inflationary potential. The curvature perturbations are generated by conversion of inflationary isocurvature perturbations, similar to curvaton and inhomogeneous reheating scenarios. During inflation, the trajectory along which the fields are rolling down will fluctuate. When inflation ends, the fields enter the oscillatory phase along different directions in different parts of the universe. Different oscillatory trajectories can have different effective equations of state and reheating rates, so that different parts of the universe will evolve differently. This translates the isocurvature perturbations into the curvature perturbations. If all light fields make similar contributions to the potential energy, the postinflationary curvature perturbations are much smaller than the standard inflationary curvature perturbations. But if different light fields make very different contributions to the potential energy, the postinflationary curvature perturbations can exceed the inflationary perturbations. The characteristic feature of postinflationary perturbations is non-Gaussianity $\sim \pm 1$.

To describe a massive graviton in 4D Minkowski space-time one introduces a quadratic term in the Lagrangian. This term, however, can lead to a readjustment or instability of the background instead of describing a massive graviton on flat space. We show that for all local Lorentz-invariant mass terms Minkowski space is unstable. We start with the Pauli-Fierz (PF) term that is the only local mass term with no ghosts in the linearized approximation. We show that nonlinear completions of the PF Lagrangian give rise to instability of Minkowski space. We continue with the mass terms that are not of a PF type. Although these models are known to have ghosts in the linearized approximations, nonlinear interactions can lead to background change due to which the ghosts are eliminated. In the latter case, however, the graviton perturbations on the new background are not massive. We argue that a consistent theory of a massive graviton on flat space can be formulated in theories with extra dimensions. They require an infinite number of fields or non-local description from a 4D point of view.

Perturbations in a universe with matter and cosmological constant are computed in terms of the primordial scalar perturbation up to second order. The calculation is easy -- a generalization of the spherical dust solution. Our expression might allow accurate measurements of the primordial non-Gaussianity.

The Particle Data Group gives an upper bound on the photon mass $m < 2 \times 10^{-16}$eV from a laboratory experiment and lists, but does not adopt, an astronomical bound $m < 3 \times 10^{-27}$eV, both of which are based on the plausible assumption of large galactic vector-potential. We argue that the interpretations of these experiments should be changed, which alters significantly the bounds on $m$. If $m$ arises from a Higgs effect, both limits are invalid because the Proca vector-potential of the galactic magnetic field may be neutralized by vortices giving a large-scale magnetic field that is effectively Maxwellian. In this regime, experiments sensitive to the Proca potential do not yield a useful bound on $m$. As a by-product, the non-zero photon mass from Higgs effect predicts generation of a primordial magnetic field in the early universe. If, on the other hand, the galactic magnetic field is in the Proca regime, the very existence of the observed large-scale magnetic field gives $m^{-1}\gtrsim 1$kpc, or $m\lesssim 10^{-26}$eV.

We generalize a recently proposed mechanism for the origin of primordial metric perturbations in inflationary models. Quantum fluctuations of light scalar fields during inflation give rise to super-horizon fluctuations of masses and reaction rates of various particles. Reheating, freeze-out, and matter-domination processes become inhomogeneous and generate super-horizon metric perturbations. We also calculate the degree of non-Gaussianity $f_{nl}$ for this new model of cosmological perturbations. The precise value of $f_{nl}$ depends on the specific models, but $|f_{nl}|\sim$few is a natural lower bound for our mechanisms. This is much larger than the currently assumed theoretical value $f_{nl}\sim tilt \lesssim 0.05$, and is thought to be observable. In a particularly attractive model of inhomogeneous mass-domination, the non-Gaussianity of perturbations generated by our mechanism is simply $f_{nl}=5$, irrespective of the detailed structure of the underlying field theory.

We propose a new mechanism to generate density perturbations in inflationary models. Spatial fluctuations in the decay rate of the inflaton field to ordinary matter lead to fluctuations in the reheating temperature. We argue that in most realistic models of inflation the coupling of the inflaton to normal matter is determined by the vacuum expectation values of fields in the theory. If those fields are light during inflation (this is a generic situation in the minimal models of supersymmetric inflation) they will fluctuate leading to density perturbations through the proposed mechanism. We show that these fluctuations could easily dominate over the ones generated through the standard mechanism. The new scenario has several consequences for inflation model building and observations. The proposed mechanism allows to generate the observed level of density perturbations with a much lower scale of inflation and thus generically predicts a smaller level of gravitational waves. The relation between the slope of the spectrum of the produced density perturbations and the potential of the inflaton field is different from the standard relations obtained in the context of slow roll inflation. Because the field responsible for the fluctuations is not the inflaton, it can have significantly larger self couplings and thus density perturbations could be non-Gaussian. The non-Gaussianity can be large enough to be detectable by CMB and Large Scale Structure observations.

Propagation of a blast wave due to strong explosion in the center of a power-law-density ($\rho \propto r^{-\alpha}$) spherically symmetric atmosphere is studied. For adiabatic index of 5/3, the solution was known to be self-similar, (of type I) for $\alpha <3$, self-similar (of type II) for $\alpha >3.26$, and unknown in between. We find a self-similar solution for $3<\alpha <3.26$, and give a (tentative) numerical proof that this solution is indeed an asymptotic of the strong explosion. This self-similar solution is neither of type I (dimensional analysis does not work), nor of type II (the index of the solution is known without solving an eigenvalue problem).

Cosmologically motivated theories that explain small acceleration rate of the Universe via modification of gravity at very large, horizon or super-horizon distances, can be tested by precision gravitational measurements at much shorter scales, such as the Earth-Moon distance. Contrary to the naive expectation the predicted corrections to the Einsteinian metric near gravitating sources are so significant that fall within sensitivity of the proposed Lunar Ranging experiments. The key reason for such corrections is the van Dam-Veltman-Zakharov discontinuity present in linearized versions of all such theories, and its subsequent absence at the non-linear level ala Vainshtein.

We reexamine the transition magnetic moment solution to the solar neutrino problem. We argue that the absence of large time variations in the Super-Kamiokande rate provides strong evidence against spin-flavor flip in the solar convective zone. Spin-flavor flip could, however, occur in the primordial magnetic field in the radiative zone. We compute the longest lived toroidal mode for this field and show that spin-flavor flip in the radiative zone can account for all available solar data.

We discuss bounds on the strength of the magnetic fields that could be buried in the radiative zone of the Sun. The field profiles and decay times are computed for all axisymmetric toroidal Ohmic decay eigenmodes with lifetimes exceeding the age of the Sun. The measurements of the solar oblateness yield a bound <~ 7 MG on the strength of the field. A comparable bound is expected to come from the analysis of the splitting of the solar oscillation frequencies. The theoretical analysis of the double diffusive instability also yields a similar bound. The oblateness measurements at their present level of sensitivity are therefore not expected to measure a toroidal field contribution.

Effective heat conduction in a random variable magnetic field should be equal to one third of the Spitzer's value. Recent observations indicate that this heat conduction is sufficient to account for the bremsstrahlung in cooling X-ray clusters.

It has been suggested that the laws of thermodynamics are violated by what we have called a convection-dominated accretion flow (or a 1/2-law accretion flow) -- an accretion flow characterized by a constant outflow of energy. We show that both the 1/2-law flow and the Bondi flow (also known as ADAF, advection dominated accretion flow) are thermodynamically admissible.

It was suggested that observations of the solar system exclude massive gravity, in the sense that the graviton mass must be rigorously zero. This is because there is a discontinuity in the linearized gravity theory at graviton mass equal to zero. The linearized Schwarzschild metric is not recovered for infinitesimal graviton mass, contradicting observations on light deviation by the Sun and Mercury perihelion advance. It was then argued that non-perturbative effects make the massive gravity theory continuous in the graviton mass. Both the original suggestion and its refutation were based on a non-covariant and linearized action, and the physical interpretation of these results remained questionable. Here we use a covariant quasi-massive gravity theory that is known to be discontinuous in the graviton mass in the linear approximation. We show that non-perturbative effects do restore the continuity; the weak-field Schwarzschild solution is recovered in the limit of small graviton mass. We also show that weak-field Schwarzschild with matter is recovered for infinitesimal graviton mass. Thus: Observations of the solar system only put an upper limit on the graviton mass (in the context of the gravity theory that we use, inverse graviton mass, as measured at distances of order inverse graviton mass, is $\gtrsim 100$ Mpc). But graviton can be massive, with a cosmologically interesting mass.

A direct numerical simulation of a strong relativistic collisionless shock propagating into an unmagnetized medium has been performed in two spatial dimensions. It is found that: (i) collisionless shock exists, (ii) particle acceleration is insignificant, (iii) the shock does not generate ``quasi-static'' magnetic fields, contrary to recent claims by other authors. All our results agree with a simple plasma-physical model of collisionless shocks. On the other hand, astrophysical evidence indicates that shocks do accelerate charged particles, and GRB phenomenology indicates that shocks do generate strong magnetic fields. The conflict between the simplest plasma-physical model and astronomical observations of collisionless shocks is discussed.

A relativistic collisionless shock propagating into an unmagnetized medium leaves behind a strong large-scale magnetic field. This seems to follow from two assumptions: (i) GRB afterglows are explained by synchrotron emission of a relativistic shock, (ii) magnetic field can't exist on microscopic scales only, it would decay by phase space mixing. Assumption (i) is generally accepted because of an apparent success of the shock synchrotron phenomenological model of GRB afterglow. Assumption (ii) is confirmed in this work by a low-dimensional numerical simulation. One may hypothesize that relativistic shock velocities are not essential for the magnetic field generation, and that all collisionless shocks propagating into an unmagnetized medium generate strong large-scale magnetic fields. If this hypothesis is true, the first cosmical magnetic fields could have been generated in shocks of the first virialized objects.

At low densities, when radiation of an accreting gas is dynamically unimportant, quasi-spherical non-radiative accretion flow is formed around a compact object. We show that the density of the flow is n\propto r^{-1/2}. The gas compression is much weaker than the standard Bondi's, n\propto r^{-3/2}, because most of the gas does not reach the compact object. It can either circulate indefinitely long or it can be launched to infinity.

Linear eigenmodes of a spherically symmetrical ultra-relativistic blast wave (the Blandford-McKee, BMK, solution) are calculated. The BMK solution is shown to be stable and strongly non-universal. It is stable because all the eigenmodes decay. Non-universality of the BMK solution (BMK is not an intermediate asymptotic) is a consequence of causality. In terms of eigenmodes -- some eigenmodes decay too slowly. For each degree l, there exists an eigenmode which decays at the smallest rate. The amplitude of this eigenmode, defined as relative perturbation of energy behind the shock front, is constant at early time. Later, when the blast wave slows down to Gamma = l/10, an oscillatory decay commences, and the amplitude drops to less than 10% of the initial value at Gamma = l/100. The non-universality is surprisingly strong. Near the end of the ultra-relativistic stage, perturbations with l<200, that is more than 10,000 different harmonics, are suppressed by less than a factor of 10. Spherical symmetry is only reached by the time the blast wave slows down to non-relativistic velocities, when the Sedov/Taylor/von Neumann solution sets in.

We review five different derivations that demonstrate that the Salpeter formula for the plasma corrections to fusion rates is valid at the center of the sun with insignificant errors (~ percent). We point out errors in several recent papers that have obtained a variety of answers, some even with the wrong sign or the wrong functional dependence.

Two possible explanations for the low luminosity of the supermassive black hole at the center of our galaxy are (1) an accretion rate of order the canonical Bondi value (roughly 10^{-5} solar masses per year), but a very low radiative efficiency for the accreting gas or (2) an accretion rate much less than the Bondi rate. Both models can explain the broad-band spectrum of the Galactic Center. We show that they can be distinguished using the linear polarization of synchrotron radiation. Accretion at the Bondi rate predicts no linear polarization at any frequency due to Faraday depolarization. Low accretion rate models, on the other hand, have much lower gas densities and magnetic field strengths close to the black hole; polarization may therefore be observable at high frequencies. If confirmed, a recent detection of linear polarization from Sgr A$^*$ above 150 GHz argues for an accretion rate of order 10^{-8} solar masses per year, much less than the Bondi rate. This test can be applied to other low-luminosity galactic nuclei.

For Gamma-Ray Bursts occurring in dense star-forming regions, the X-ray afterglow behavior minutes to days after the trigger may be dominated by the small-angle scattering of the prompt X-ray emission off dust grains. We give a simple illustrative model for the X-ray light curves at different X-ray energies, and discuss possible implications. A bump followed by a steeper decay in soft X-rays is predicted for bursts which are heavily obscured in the optical.

We show how an idealized high-frequency VLBI can probe the metric of nearby supermassive black holes. No quantitative assumptions about plasma dynamics are used in the suggested diagnostic, which is based on strong-field lensing.

MeV seed photons produced in shocks in a variable ultra-relativistic outflow gain energy by the Fermi mechanism, because the photons Compton scatter off relativistically colliding shells. The Fermi-modified high-energy photon spectrum has a non-universal slope and a universal cutoff. A significant increase in the total radiative efficiency is possible. In some gamma ray bursts, most of the power might be emitted at the high-energy cutoff for this mechanism, which would be close to 100 MeV for outflows with a mean bulk Lorentz factor of 100.

Cold dark matter (CDM) models predict small-scale structure in excess of observations of the cores and abundance of dwarf galaxies. These problems might be solved, and the virtues of CDM models retained, even without postulating {\it ad hoc} dark matter particle or field interactions, if the dark matter is composed of ultra-light scalar particles (m ~ 10^{-22}eV), initially in a (cold) Bose-Einstein condensate, similar to axion dark matter models. The wave properties of the dark matter stabilize gravitational collapse providing halo cores and sharply suppressing small-scale linear power.

Two-temperature spherical accretion flows produce $\approx 100$ Mev gamma-rays from the decay of neutral pions created in proton-proton collisions close to the black hole; they also produce $\sim 10$ keV X-rays by bremsstrahlung emission at large radii. The gamma-ray to X-ray luminosity ratio is nearly independent of black hole mass and accretion rate. It does depend sensitively on the radial density profile of the accretion flow through the parameter a, where $n \propto r^{-a}$. For the canonical Bondi value of a = 3/2, the gamma-ray to X-ray luminosity ratio is $\approx 30$. We interpret a recent Chandra detection coincident with the massive black hole at the Galactic Center as being thermal bremsstrahlung emission from the accretion flow. With this normalization, the expected gamma-ray luminosity is $\approx 10^{35}$ ergs s$^{-1}$ if a = 3/2. This is nearly two orders of magnitude above the detection threshold of the GLAST telescope. For $a \approx 1/2$, however, (a value suggested by recent theoretical arguments), the expected gamma-ray luminosity is only $\approx 10^{29}$ ergs s$^{-1}$; GLAST should therefore provide an important probe of the true accretion rate and radial density profile of the accretion flow onto Sgr A*.

Non-radiating, advection-dominated, accretion flows are convectively unstable. We calculate the two-dimensional (r-theta) structure of such flows assuming that (1) convection transports angular momentum inwards, opposite to normal viscosity and (2) viscous transport by other mechanisms (e.g., magnetic fields) is weak (alpha << 1). Under such conditions convection dominates the dynamics of the accretion flow and leads to a steady state structure that is marginally stable to convection. We show that the marginally stable flow has a constant temperature and rotational velocity on spherical shells, a net flux of energy from small to large radii, zero net accretion rate, and a radial density profile proportional to r^{-1/2}, flatter than the r^{-3/2} profile characteristic of spherical accretion flows. This solution accurately describes the full two-dimensional structure of recent axisymmetric numerical simulations of advection-dominated accretion flows.

When a bright spot in a high-inclination disk orbits a black hole, its light is periodically lensed and Doppler shifted. If the spots occupy a narrow range of radii, a quasi-periodic oscillation (QPO) will result. I calculate the QPO frequency assuming that the spots appear near r_max - the radius of the maximal proper radiation flux from the disk. The calculated frequency depends on the black hole mass and spin. For the microquasar GRO J1655-40, the black hole mass is known and a 300 Hz QPO was observed. The inferred black hole spin is about 60% of the maximal. The orbit precession frequency at r_max is 20Hz, this QPO was also observed in GRO J1655-40. An approximate formula for r_max is given.

I have shown that the enhancement of fusion rates in stars due to electrostatic screening does not depend on the Gamow energy of the fusing nuclei (in the classical weak screening limit). Opher and Opher have recently published an ApJ paper claiming that my proof is wrong. I show here that their arguments are incorrect.

A Kerr black hole can flare quasi-periodically, if it is magnetically connected to the accretion disk. The flares extract the rotational energy of the hole. The flaring occurs because the hole twists the magnetic field lines beyond the threshold of the screw instability.

When a black hole accretes slowly, the radiative cooling of the infalling gas is weak and the accretion disk does not form. A hot collisionless quasi-spherical halo is formed instead. The properties of such halos are discussed. The rate of accretion, the radiative efficiency, and the temperature and density near the hole are evaluated.

The optical afterglows of the gamma ray bursts can be strongly polarized, in principle up to tens of percents, if: (i) the afterglow is synchrotron radiation from an ultra-relativistic blast, (ii) the blast is beamed during the afterglow phase, i.e. the shock propagates within a narrow jet, (iii) we observe at the right time from the right viewing angle, (iv) magnetic fields parallel and perpendicular to the jet have different proper strengths.

We show that flows of pure electromagnetic energy are subject to instabilities familiar from hydrodynamics and plasma physics---Kelvin-Helmholtz and screw instabilities. In the framework of force-free electrodynamics, it is found that electric-like tangential discontinuities are Kelvin-Helmholtz-unstable. All non-trivial cylindrically symmetrical magnetic configurations are screw-unstable. The Poynting jet of the Goldreich-Julian pulsar is screw-unstable if the current density in the polar cap region exceeds the charge density times the speed of light. This may be the process that determines the pulsar luminosity.

The rate of turbulent spherical accretion onto a compact object might be much smaller than the Bondi rate. It is suggested that the rate of accretion onto Sgr A-star is much smaller than the Bondi rate.

We calculate the proton distribution function due to heating by subsonic (Alfvenic) turbulence in a weakly magnetized collisionless plasma. The distribution function is nonthermal. For non-relativistic energies, it is an exponential of the magnitude of the proton velocity. For ultra-relativistic energies, it can be characterized as a power law with a momentum-dependent slope.

The observed frequency shifts of the high-degree solar fundamental mode are explained using a simple geometrical optics approximation. The predicted fractional frequency shift is proportional to the mean squared velocity of the convection flow at the surface of the Sun and inversely proportional to the group velocity of the mode squared. The data converge to the predicted behavior at very high degrees. The inferred root mean squared convection velocity is 1.3 km/s.

GRB afterglow polarization is discussed. We find an observable, up to 10%, polarization, if the magnetic field coherence length grows at about the speed of light after the field is generated at the shock front. Detection of a polarized afterglow would show that collisionless ultrarelativistic shocks can generate strong large scale magnetic fields and confirm the synchrotron afterglow model. Non-detection, at a 1% level, would imply that either the synchrotron emission model is incorrect, or that strong magnetic fields, after they are generated in the shock, somehow manage to stay un-dissipated at ``microscopic'', skin depth, scales. Analytic lightcurves of synchrotron emission from an ultrarelativistic self-similar blast wave are obtained for an arbitrary electron distribution function, taking into account the effects of synchrotron cooling. The peak synchrotron flux and the flux at frequencies much smaller than the peak frequency are insensitive to the details of the electron distribution function; hence their observational determination would provide strong constraints on blast wave parameters.

We review and analyze the available information for nuclear fusion cross sections that are most important for solar energy generation and solar neutrino production. We provide best values for the low-energy cross-section factors and, wherever possible, estimates of the uncertainties. We also describe the most important experiments and calculations that are required in order to improve our knowledge of solar fusion rates.

In a universe reionized in patches, the Doppler effect from Thomson scattering off free electrons generates secondary cosmic microwave background (CMB) anisotropies. For a simple model with small patches and late reionization, we analytically calculate the anisotropy power spectrum. Patchy reionization can, in principle, be the main source of anisotropies on arcminute scales. On larger angular scales, its contribution to the CMB power spectrum is a small fraction of the primary signal and is only barely detectable in the power spectrum with even an ideal, i.e. cosmic variance limited, experiment and an extreme model of reionization. Consequently patchy reionization is unlikely to affect cosmological parameter estimation from the acoustic peaks in the CMB. Its detection on small angles would help determine the ionization history of the universe, in particular the typical size of the ionized region and the duration of the reionization process.

We extend and reconcile recent work on turbulence and particle heating in advection-dominated accretion flows. For approximately equipartition magnetic fields, the turbulence primarily heats the electrons. For weaker magnetic fields, the protons are primarily heated. The division between electron and proton heating occurs between $\beta \sim 5$ and $\beta \sim 100$ (where $\beta$ is the ratio of the gas to the magnetic pressure), depending on unknown details of how Alfv\'en waves are converted into whistlers on scales of the proton Larmor radius. We also discuss the possibility that magnetic reconnection could be a significant source of electron heating.

We evaluate the effect of electrostatic screening by ions and electrons on low-Z thermonuclear reactions in the sun. We use a mean field formalism and calculate the electron density of the screening cloud using the appropriate density matrix equation of quantum statistical mechanics. Because of well understood physical effects that are included for the first time in our treatment, the calculated enhancement of reaction rates does not agree with the frequently used interpolation formulae. Our result does agree, within small uncertainties, with Salpeter's weak screening formula. If weak screening is used instead of the commonly employed screening prescription of Graboske et al., the predicted $^8$B neutrino flux is increased by 7% and the predicted chlorine rate is increased by 0.4 SNU.

Radiative efficiency of a slowly accreting black hole is estimated using a two-temperature model of accretion. The radiative efficiency depends on the magnetic field strength near the Schwarzschild radius. For weak magnetic fields (magnetic energy=equipartition/1000), the low efficiency 0.0001 assumed in some theoretical models might be achieved. For stronger fields, a significant fraction of viscous heat is dissipated by electrons and radiated away resulting in a larger efficiency. At equipartition magnetic fields, we estimate efficiency = of order 10%.

Perturbation theory is used to calculate the frequency shift of acoustic modes of an homogeneous turbulent fluid and the frequency shifts of solar modes due to turbulent convection. For sound waves in a random flow, the fractional frequency shift is +(11/30)M^2 in the long-wavelength limit, and the shift at short wavelengths is -(2/3)M^2, where M is the average Mach number of the flow. In our model of the solar convection, the low-degree f-mode shift is positive and unobservably small, whereas the fractional frequency shift at high degrees (l\gg 500) is -M^2, where M is the Mach number of convection near the solar surface.

For solar conditions, we numerically integrate the density matrix equation for a thermal electron in the field of a Be-7 ion and other plasma ions and smeared-out electrons. Our results are in agreement with previous calculations that are based on a different physical picture, a picture which postulates the existence of distinct continuum and bound state orbits for electrons. The density matrix calculation of the electron capture rate is independent of the nature of electron states in the solar plasma. To within a 1% accuracy, the effects of screening can be described at high temperatures by a Salpeter-like factor which can be derived from the density matrix equation. The total theoretical uncertainty in the electron capture rate is about 2%.

We show that there are no dynamic screening corrections to the Salpeter's enhancement factor in the weak-screening limit.

Generation of magnetic field energy, without mean field generation, is studied. Isotropic mirror-symmetric turbulence of a conducting fluid amplifies the energy of small-scale magnetic perturbations if the magnetic Reynolds number is high, and the dimensionality of space d satisfies 2.103 < d <8.765. The result does not depend on the model of turbulence, incompressibility and isotropy being the only requirements.