Numerical Simulations of Equatorially-Asymmetric Magnetized Supernovae: Formation of Magnetars and Their Kicks

A series of numerical simulations on magnetorotational core-collapse supernovae are carried out. Dipole-like configurations which are offset northward are assumed for the initially strong magnetic fields together with rapid differential rotations. Aims of our study are to investigate effects of the offset magnetic field on magnetar kicks and on supernova dynamics. Note that we study a regime where the proto-neutron star formed after collapse has a large magnetic field strength approaching that of a ``magnetar'', a highly magnetized slowly rotating neutron star. As a result, equatorially-asymmetric explosions occur with a formation of the bipolar jets. Resultant magnetar's kick velocities are $\sim 300-1000$ km s$^{-1}$. We find that the acceleration is mainly due to the magnetic pressure while the somewhat weaker magnetic tension works toward the opposite direction, which is due to stronger magnetic field in the northern hemisphere. Noted that observations of magnetar's proper motions are very scarce, our results supply a prediction for future observations. Namely, magnetars possibly have large kick velocities, several hundred km s$^{-1}$, as ordinary neutron stars do, and in an extreme case they could have those up to 1000 km s$^{-1}$.Comment: 36 pages, 9 figures, accepted by the Astrophysical Journa

Gravitational waves from supernova matter

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.Comment: 10 pages, 6 figures, accepted, to be published in a Classical and Quantum Gravity special issue for MICRA200

Gravitational Waves from Core Collapse Supernovae

We present the gravitational wave signatures for a suite of axisymmetric core collapse supernova models with progenitors masses between 12 and 25 solar masses. These models are distinguished by the fact they explode and contain essential physics (in particular, multi-frequency neutrino transport and general relativity) needed for a more realistic description. Thus, we are able to compute complete waveforms (i.e., through explosion) based on non-parameterized, first-principles models. This is essential if the waveform amplitudes and time scales are to be computed more precisely. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models.Comment: 10 pages, 5 figure

Equilibrium Configurations of Strongly Magnetized Neutron Stars with Realistic Equations of State

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). Employing the Tomimura-Eriguchi scheme to construct the equilibrium configurations. we study the basic physical properties of the sequences in the framework of Newton gravity. In addition we newly take into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g., structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyze this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than $10^{15}\rm{g} \rm{cm}^{-3}$. The maximum baryon mass of the magnetized stars with axis ratio $q\sim 0.7$ increases about up to twenty percents for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly-born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.Comment: submitted to MNRA

Light Curves and Event Rates of Axion Instability Supernovae

It was recently proposed that exotic particles can trigger a new stellar instability which is analogous to the e-e+ pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ~400 keV--2 MeV and coupling with photons of g_{ag}~10^{-5} GeV^{-1}. It is found that the 56Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10--20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7--2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function.Comment: 9 pages, 6 figures, 1 table, submitted to Ap

A local grid refinement technique based upon Richardson extrapolation

A grid-embedding technique for the solution of two-dimensional incompressible flows governed by the Navier-Stokes equations is presented. A single coarse grid covers the whole domain, and local grid refinement B carried out in the regions of high gradients without changing the basic grid structure. A finite volume method with collocated primitive variables is employed, ensuring conservation at the interfaces of embedded grids, as well as global conservation. The method is applied to the simulation of a turbulent flow past a backward facing step, the flow over a square obstacle, and the flow in a sudden pipe expansion, and the predictions are compared with data published in the literature. They show that neither the convergence rate nor the stability of the method are affected by the presence of embedded grids. The grid-embedding technique yields significant savings in computing time to achieve the same accuracy obtained wing conventional grids. (C) 1997 by Elsevier Science Inc

General relativistic simulations of pasive-magneto-rotational core collapse with microphysics

This paper presents results from axisymmetric simulations of magneto-rotational stellar core collapse to neutron stars in general relativity using the passive field approximation for the magnetic field. These simulations are performed using a new general relativistic numerical code specifically designed to study this astrophysical scenario. The code is based on the conformally-flat approximation of Einstein's field equations and conservative formulations of the magneto-hydrodynamics equations. The code has been recently upgraded to incorporate a tabulated, microphysical equation of state and an approximate deleptonization scheme. This allows us to perform the most realistic simulations of magneto-rotational core collapse to date, which are compared with simulations employing a simplified (hybrid) equation of state, widely used in the relativistic core collapse community. Furthermore, state-of-the-art (unmagnetized) initial models from stellar evolution are used. In general, stellar evolution models predict weak magnetic fields in the progenitors, which justifies our simplification of performing the computations under the approach that we call the passive field approximation for the magnetic field. Our results show that for the core collapse models with microphysics the saturation of the magnetic field cannot be reached within dynamical time scales by winding up the poloidal magnetic field into a toroidal one. We estimate the effect of other amplification mechanisms including the magneto-rotational instability (MRI) and several types of dynamos.Comment: 25 pages, 15 figures, accepted for publication in Astronomy & Astrophysics July 31, 2007. Added 1 figure and a new subsectio

CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics

We present a formulation for multigroup radiation hydrodynamics that is correct to order $O(v/c)$ using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts, one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.Comment: accepted by ApJS, 27 pages, 20 figures, high-resolution version available at https://ccse.lbl.gov/Publications/wqzhang/castro3.pd

The patterned structure of galactoglucomannan suggests it may bind to cellulose in seed mucilage

FAPESP - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DE SÃO PAULOThe interaction between mannan polysaccharides and cellulose microfibrils contributes to cell wall properties in some vascular plants, but the molecular arrangement of mannan in the cell wall and the nature of the molecular bonding between mannan and cell178310111026FAPESP - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DE SÃO PAULOFAPESP - FUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DE SÃO PAULO2013/08293-7The authors would like to acknowledge Prof. George Lomonossoff (John Innes Centre, UK), who developed the pEAQ-HyperTrans expression system used in this study. Plant Bioscience Limited supplied the pEAQ-HT vector that was used in this wor