A Two-Dimensional MagnetoHydrodynamics Scheme for General Unstructured Grids

We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of div(\bB) by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant.Comment: 22 pages, including 11 figures; Accepted to the Astrophysical Journal. Higher resolution figures available at http://zenith.as.arizona.edu/~burrows/mhd-code

North-South Neutrino Heating Asymmetry in Strongly Magnetized and Rotating Stellar Cores

We perform a series of two-dimensional magnetohydrodynamic simulations of supernova cores. Since the distributions of the angular momentum and the magnetic fields of strongly magnetized stars are quite uncertain, we systematically change the combinations of the strength of the angular momentum, the rotations law, the degree of differential rotation, and the profiles of the magnetic fields to construct the initial conditions. By so doing, we estimate how the rotation-induced anisotropic neutrino heating are affected by the strong magnetic fields through parity-violating effects and first investigate how the north-south asymmetry of the neutrino heating in a strongly magnetized supernova core could be. As for the microphysics, we employ a realistic equation of state based on the relativistic mean field theory and take into account electron captures and the neutrino transport via the neutrino leakage scheme. With these computations, we find that the parity-violating corrections reduce $\lesssim 0.5$ of the neutrino heating rate than that without the magnetic fields in the vicinity of the north pole of a star, on the other hand, enhance about $\lesssim 0.5$ in the vicinity of the south pole. If the global asymmetry of the neutrino heating in the both of the poles develops in the later phases, the newly born neutron star might be kicked toward the north pole in the subsequent time.Comment: 25 pages, 6 figures, ApJ in press. A paper with higher-resolution figures available at http://www-utap.phys.s.u-tokyo.ac.jp/~kkotake/lonbun.htm

Dirac neutrino magnetic moment and the shock wave revival in a supernova explosion

The process of the two-step conversion of the neutrino helicity, $\nu_L \to \nu_R \to \nu_L$, is analysed in the supernova conditions, where the first stage is realized due to the interaction of the neutrino magnetic moment with the plasma electrons and protons in the supernova core. The second stage is caused by the neutrino resonant spin-flip in a magnetic field of the supernova envelope. Given the neutrino magnetic moment within the interval $10^{-13} \mu_{\rm B} < \mu_\nu < 10^{-12} \mu_{\rm B}$, and with the existence of the magnetic field at the scale $\sim 10^{13}$ G between the neutrinosphere and the shock-wave stagnation region, it is shown that an additional energy of the order of $10^{51}$ erg can be injected into this region during the typical time of the shock-wave stagnation. This energy could be sufficient for stumulation of the damped shock wave.Comment: 6 pages, LaTeX, 2 PS figures, based on the talk presented by N.V. Mikheev at the XV International Seminar Quarks'2008, Sergiev Posad, Moscow Region, May 23-29, 2008, to appear in the Proceeding

Observational Constraints on the Angular and Spectral Distributions of Photons in Gamma-Ray Burst Sources

The typical spectra of gamma-ray bursts (GRBs) are discussed in the context of the compactness problem for GRB sources and how it is resolved in the popular fireball model. In particular, observational (model-independent) constraints on the collimation of the gamma-rays and the dependence of the collimation angle on the photon energy are considered. The fact that the threshold for the creation of $e^{-}e^{+}$ pairs depends on the angle between the momenta of the annihilating photons in the GRB source provides an alternative solution to the compactness problem. A new approach to explaining GRBs, taking into account the angular dependence for pair creation, is proposed, and the main features of a scenario describing a GRB source with a total (photon) energy smaller or of the order of $10^{49}$ erg are laid out. Thus, we are dealing with an alternative to an ultra-relativistic fireball, if it turns out (as follows from observations) that all "long" GRBs are associated with normal (not peculiar) core-collapse supernovae. The effects of radiation pressure and the formation of jets as a consequence of even a small amount of anisotropy in the total radiation field in a (compact) GRB source are examined in this alternative model. Possible energy release mechanisms acting in regions smaller or of the order of $10^{8}$ cm in size (a compact model for a GRB) are discussed. New observational evidence for such compact energy release in the burst source is considered.Comment: 15 pages, no figures, no table

How Massive Single Stars End their Life

How massive stars die -- what sort of explosion and remnant each produces -- depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are chiefly a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided single stars rotate rapidly enough at death, we speculate upon stellar populations that might produce gamma-ray bursts and jet-driven supernovae.Comment: 24 pages, 9 figues, submitted to Ap

A Dipole Vortex Model of Obscuring Tori in Active Galaxy Nuclei

The torus concept as an essential structural component of active galactic nuclei (AGN) is generally accepted. Here, the situation is discussed when the torus "twisting" by the radiation or wind transforms it into a dipole toroidal vortex which in turn can be a source of matter replenishing the accretion disk. Thus emerging instability which can be responsible for quasar radiation flares accompanied by matter outbursts is also discussed. The "Matreshka" scheme for an obscuring vortex torus structure capable of explaining the AGN variability and evolution is proposed. The model parameters estimated numerically for the luminosity close to the Eddington limit agree well with the observations.Comment: 17 pages, 11 figures, version of this paper is published in Astronomy Report

Physical Limits of Different Models of Cosmic Gamma-Ray Bursts

The present common view about GRB origin is related to cosmology, what is based on statistical analysis, and on measurements of the redshifts in the GRB optical afterglows of long GRB. No correlation is found between redshifts, GRB spectrum, and total GRB fluence. Comparison of KONUS and BATSE data about statistics and hard X-ray lines is done, and some differences are noted. Hard gamma-ray afterglows, prompt optical spectra, hard X-ray lines measurements could be important for farther insight into GRB origin. Possible possible connection of short GRB with soft gamma repeaters is discussed.Comment: 12 pages, 11 figures, Invited paper presented at the ``Frascati Workshop 2003'', Vulcano, Italy, 26-31 May, 200

The Proto-Magnetar Model for Gamma-Ray Bursts

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.Comment: 21 pages (plus 2 appendices), 21 figures, 1 table, now accepted to MNRA

A missense variant in CST3 exerts a recessive effect on susceptibility to age-related macular degeneration resembling its association with Alzheimer’s disease

Age-related macular degeneration (AMD) and Alzheimer’s disease (AD) are degenerative, multifactorial diseases involving age-related accumulation of extracellular deposits linked to dysregulation of protein homeostasis. Here, we strengthen the evidence that an nsSNP (p.Ala25Thr) in the cysteine proteinase inhibitor cystatin C gene CST3, previously confirmed by meta-analysis to be associated with AD, is associated with exudative AMD. To our knowledge, this is the first report highlighting a genetic variant that increases the risk of developing both AD and AMD. Furthermore, we demonstrate that the risk associated with the mutant allele follows a recessive model for both diseases. We perform an AMD-CST3 case–control study genotyping 350 exudative AMD Caucasian individuals. Bringing together our data with the previously reported AMD-CST3 association study, the evidence of a recessive effect on AMD risk is strengthened (OR = 1.89, P = 0.005). This effect closely resembles the AD-CST3 recessive effect (OR = 1.73, P = 0.005) previously established by meta-analysis. This resemblance is substantiated by the high correlation between CST3 genotype and effect size across the two diseases (R2 = 0.978). A recessive effect is in line with the known function of cystatin C, a potent enzyme inhibitor. Its potency means that, in heterozygous individuals, a single functional allele is sufficient to maintain its inhibitory function; only homozygous individuals will lack this form of proteolytic regulation. Our findings support the hypothesis that recessively acting variants account for some of the missing heritability of multifactorial diseases. Replacement therapy represents a translational opportunity for individuals homozygous for the mutant allele

Magnetars, Gamma-ray Bursts, and Very Close Binaries

We consider the possible existence of a common channel of evolution of binary systems, which results in a gamma-ray burst during the formation of a black hole or the birth of a magnetar during the formation of a neutron star. We assume that the rapid rotation of the core of a collapsing star can be explained by tidal synchronization in a very close binary. The calculated rate of formation of rapidly rotating neutron stars is qualitatively consistent with estimates of the formation rate of magnetars. However, our analysis of the binarity of newly-born compact objects with short rotational periods indicates that the fraction of binaries among them substantially exceeds the observational estimates. To bring this fraction into agreement with the statistics for magnetars, the additional velocity acquired by a magnetar during its formation must be primarily perpendicular to the orbital plane before the supernova explosion, and be large.Comment: 16 pages, 4 figure