Effects of Rotation and Relativistic Charge Flow on Pulsar Magnetospheric Structure

We propose an analytical 3-D model of the open field-line region of a neutron star (NS) magnetosphere. We construct an explicit analytic solution for arbitrary obliquity (angle between the rotation and magnetic axes) incorporating the effects of magnetospheric rotation, relativistic flow of charges (e.g. primary electron beam) along the open field lines, and E X B drift of these charges. Our solution employs the space-charge-limited longitudinal current calculated in the electrodynamic model of Muslimov & Tsygan (1992) and is valid up to very high altitudes nearly approaching the light cylinder. We assume that in the innermost magnetosphere, the NS magnetic field can be well represented by a static magnetic dipole configuration. At high altitudes the open magnetic field lines significantly deviate from those of a static dipole and tend to focus into a cylindrical bundle, swept back in the direction opposite to the rotation, and also bent towards the rotational equator. We briefly discuss some implications of our study to spin-powered pulsars.Comment: 24 pages, 3 figures, accepted for publication in Ap

Plasma Modes Along the Open Field Lines of a Neutron Star

We consider electrostatic plasma modes along the open field lines of a rotating neutron star. Goldreich-Julian charge density in general relativity is analyzed for the neutron star with zero inclination. It is found that the charge density is maximum at the polar cap and it remains almost same in certain extended region of the pole. For a steady state Goldreich-Julian charge density we found the usual plasma oscillation along the field lines; plasma frequency resembles to the gravitational redshift close to the Schwarzschild radius. We study the nonlinear plasma mode along the field lines. From the system of equations under general relativity, a second order differential equation is derived. The equation contains a term which describes the growing plasma modes near Schwarzschild radius in a black hole environment. The term vanishes with the distance far away from the gravitating object. For initially zero potential and field on the surface of a neutron star, Goldreich-Julian charge density is found to create the plasma mode, which is enhanced and propagates almost without damping along the open field lines. We briefly outline our plan to extend the work for studying soliton propagation along the open field lines of strongly gravitating objects

Limits on the Halo White Dwarf Component of Baryonic Dark Matter from the {\em Hubble Deep Field}

The MACHO collaboration lensing event statistics suggest that a significant fraction of the dark galactic halo can be comprised of baryonic matter in the form of white dwarf stars with masses between 0.1 and 1.0 \Msun . Such a halo white dwarf population, in order to have escaped detection by those who observe the white dwarf luminosity function of the disk, must have formed from an old population. The observations indicate that the number of halo white dwarfs per cubic parsec per unit bolometric magnitude is less than $10^{-5}$ at $10^{-4.5}$\Lsun; the number must rise significantly at lower luminosities to provide the needed baryonic halo mass. Such white dwarfs may easily escape detection in most current and earlier surveys. Though it is limited in angular extent, the {\em Hubble Deep Field} (HDF) probes a sufficient volume of the galactic halo to provide interesting limits on the number of halo white dwarf stars, and on the fraction of the halo mass that they can make up. If the HDF field can be probed for stars down to $V=29.8$ then the MACHO result suggests that there could be up to 12 faint halo white dwarfs visible in the HDF. Finding (or not finding) these stars in turn places interesting constraints on star formation immediately following the formation of the galaxy.Comment: 10 pages, AASTEX, 1 table, no figures, accepted for publication in Ap.J. Letter

The Axisymmetric Pulsar Magnetosphere

We present, for the first time, the structure of the axisymmetric force-free magnetosphere of an aligned rotating magnetic dipole, in the case in which there exists a sufficiently large charge density (whose origin we do not question) to satisfy the ideal MHD condition, ${\bf E\cdot B}=0$, everywhere. The unique distribution of electric current along the open magnetic field lines which is required for the solution to be continuous and smooth is obtained numerically. With the geometry of the field lines thus determined we compute the dynamics of the associated MHD wind. The main result is that the relativistic outflow contained in the magnetosphere is not accelerated to the extremely relativistic energies required for the flow to generate gamma rays. We expect that our solution will be useful as the starting point for detailed studies of pulsar magnetospheres under more general conditions, namely when either the force-free and/or the ideal MHD condition ${\bf E\cdot B}=0$ are not valid in the entire magnetosphere. Based on our solution, we consider that the most likely positions of such an occurrence are the polar cap, the crossings of the zero space charge surface by open field lines, and the return current boundary, but not the light cylinder.Comment: 15 pages AAS Latex, 5 postscript figure

Magnetar Spindown, Hyper-Energetic Supernovae, and Gamma Ray Bursts

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength ($\sim10^{15}$ G) large scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch.Because the strong magnetic field forces the wind to co-rotate with the protoneutron star,this outflow can significantly effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber & Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown timescale is short with respect to the time required for the supernova shockwave to traverse the stellar progenitor, the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino-magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods ($P$) of $\sim1$ millisecond, that of order $10^{51}-10^{52}$ erg of rotational energy can be extracted in $\sim10$ seconds. If magnetars are born slowly rotating ($P\gtrsim10$ ms) they can spin down to periods of $\sim1$ second on the Kelvin-Helmholtz timescale.Comment: 16 pages, 5 figures, emulateap

Post-T Tauri stars: a false problem

We consider the problem of the apparent lack of old T Tauri stars in low-mass star forming regions in the framework of the standard model of low-mass star formation. We argue that the similarity between molecular cloud lifetime and ambipolar diffusion timescale implies that star formation does not take place instantaneously, nor at a constant rate. We conclude that the probability of finding a large population of old stars in a star forming region is intrinsically very small and that the post-T Tauri problem is by and large not existent.Comment: 6 pages (LaTeX), no Figures to be published in The Astrophysical Journal Letter

Nonlinear Criterion for the Stability of Molecular Clouds

Dynamically significant magnetic fields are routinely observed in molecular clouds, with mass-to-flux ratio lambda = (2 pi sqrt{G}) (Sigma/B) ~ 1 (here Sigma is the total column density and B is the field strength). It is widely believed that ``subcritical'' clouds with lambda < 1 cannot collapse, based on virial arguments by Mestel and Spitzer and a linear stability analysis by Nakano and Nakamura. Here we confirm, using high resolution numerical models that begin with a strongly supersonic velocity dispersion, that this criterion is a fully nonlinear stability condition. All the high-resolution models with lambda <= 0.95 form ``Spitzer sheets'' but collapse no further. All models with lambda >= 1.02 collapse to the maximum numerically resolvable density. We also investigate other factors determining the collapse time for supercritical models. We show that there is a strong stochastic element in the collapse time: models that differ only in details of their initial conditions can have collapse times that vary by as much as a factor of 3. The collapse time cannot be determined from just the velocity dispersion; it depends also on its distribution. Finally, we discuss the astrophysical implications of our results.Comment: 11 pages, 5 figure

Fragmentation Instability of Molecular Clouds: Numerical Simulations

We simulate fragmentation and gravitational collapse of cold, magnetized molecular clouds. We explore the nonlinear development of an instability mediated by ambipolar diffusion, in which the collapse rate is intermediate to fast gravitational collapse and slow quasistatic collapse. Initially uniform stable clouds fragment into elongated clumps with masses largely determined by the cloud temperature, but substantially larger than the thermal Jeans mass. The clumps are asymmetric, with significant rotation and vorticity, and lose magnetic flux as they collapse. The clump shapes, intermediate collapse rates, and infall profiles may help explain observations not easily fit by contemporary slow or rapid collapse models.Comment: 25pp, 20 small eps figures, in press ApJ, April 1, 200

Charged Condensate and Helium Dwarf Stars

White dwarf stars composed of carbon, oxygen or heavier elements are expected to crystallize as they cool down below certain temperatures. Yet, simple arguments suggest that the helium white dwarf cores may not solidify, mostly because of zero-point oscillations of the helium ions that would dissolve the crystalline structure. We argue that the interior of the helium dwarfs may instead form a macroscopic quantum state in which the charged helium-4 nuclei are in a Bose-Einstein condensate, while the relativistic electrons form a neutralizing degenerate Fermi liquid. We discuss the electric charge screening, and the spectrum of this substance, showing that the bosonic long-wavelength fluctuations exhibit a mass gap. Hence, there is a suppression at low temperatures of the boson contribution to the specific heat -- the latter being dominated by the specific heat of the electrons near the Fermi surface. This state of matter may have observational signatures.Comment: 10 pages; v2: to appear in JCAP, brief comments and section titles added, typos correcte

Optimal battery charge/discharge strategies for prosumers and suppliers

We discuss the application of classical variational methods to optimal charging/discharging strategies for a prosumer or storage supplier, where the price of electrical power is known in advance. We outline how a classical calculus of variations approach can be applied to two related problems: (i) how can a prosumer minimise the cost of charging/discharging a battery, when the price of electrical power is known throughout the charging/discharging period? and (ii) how can an electricity supplier incentivise desired prosumer/storage supplier behaviour by adjusting the price