Two-Dimensional Hydrodynamics of Pre-Core Collapse: Oxygen Shell Burning

By direct hydrodynamic simulation, using the Piecewise Parabolic Method (PPM) code PROMETHEUS, we study the properties of a convective oxygen burning shell in a SN 1987A progenitor star prior to collapse. The convection is too heterogeneous and dynamic to be well approximated by one-dimensional diffusion-like algorithms which have previously been used for this epoch. Qualitatively new phenomena are seen. The simulations are two-dimensional, with good resolution in radius and angle, and use a large (90-degree) slice centered at the equator. The microphysics and the initial model were carefully treated. Many of the qualitative features of previous multi-dimensional simulations of convection are seen, including large kinetic and acoustic energy fluxes, which are not accounted for by mixing length theory. Small but significant amounts of carbon-12 are mixed non-uniformly into the oxygen burning convection zone, resulting in hot spots of nuclear energy production which are more than an order of magnitude more energetic than the oxygen flame itself. Density perturbations (up to 8%) occur at the `edges' of the convective zone and are the result of gravity waves generated by interaction of penetrating flows into the stable region. Perturbations of temperature and electron fraction at the base of the convective zone are of sufficient magnitude to create angular inhomogeneities in explosive nucleosynthesis products, and need to be included in quantitative estimates of yields. Combined with the plume-like velocity structure arising from convection, the perturbations will contribute to the mixing of nickel-56 throughout supernovae envelopes. Runs of different resolution, and angular extent, were performed to test the robustness of theseComment: For mpeg movies of these simulations, see http://www.astrophysics.arizona.edu/movies.html Submitted to the Astrophysical Journa

Re-test of Rhinocyllus conicus host specificity, and the prediction of ecological risk in biological control

Biological control is proposed as an ecological strategy to manage the threat of invasive plants, especially in natural areas. To pursue this strategy, we need to know that the host specificity criteria used to evaluate ecological risk with deliberate introduction of an exotic insect for biocontrol are sufficient to predict potential impact on native species. Host specificity is defined by adult feeding and oviposition preferences and larval development. One way to evaluate the criteria is to re-examine case histories where ecological effects are recorded, such as that of Rhinocyllus conicus Frölich. This flower head weevil, released in North America in 1968 to control exotic thistles like Musk thistle (Carduus nutans L), is now reducing seed production by multiple native North American thistle species (Cirsium spp.), and local population density of Platte thistle (Cirsium canescens Nutt.). We hypothesized that host specificity of R. conicus has changed since pre-release testing, providing an explanation for the unexpected magnitude of the documented ecological effects. Instead, when we re-tested host specificity of weevils naturalized over 28 generations, we found that host specificity has not changed. Naturalized adults of R. conicus showed strong feeding and oviposition preference for Musk thistle over Platte thistle. In addition, larval development by these weevils was faster and more successful (to larger size) on Musk thistle than on Platte thistle. Thus, our results indicate that a change in host specificity cannot explain the unexpectedly large build-up of R. conicus and significant ecological effect on Platte thistle. We conclude that accurate prediction of the potential level of impact on native host plants in the field requires further ecological information in addition to host specificity

Evidence for a Mid-Atomic-Number Atmosphere in the Neutron Star 1E1207.4-5209

Recently Sanwal et al. (2002) reported the first clear detection of absorption features in an isolated neutron star, 1E1207.4-5209. Remarkably their spectral modeling demonstrates that the atmosphere cannot be Hydrogen. They speculated that the neutron star atmosphere is indicative of ionized Helium in an ultra-strong (~1.5x10^{14} G) magnetic field. We have applied our recently developed atomic model (Mori & Hailey 2002) for strongly-magnetized neutron star atmospheres to this problem. We find that this model, along with some simp le atomic physics arguments, severely constrains the possible composition of the atmosphere. In particular we find that the absorption features are naturally associated with He-like Oxygen or Neon in a magnetic field of ~10^{12} G, comparable to the magnetic field derived from the spin parameters of the neutron star. This interpretation is consistent with the relative line strengths and widths and is robust. Our model predicts possible substructure in the spectral features, which has now been reported by XMM-Newton (Mereghetti et al. 2002). However we show the Mereghetti et al. claim that the atmosphere is Iron or some comparable high-Z element at ~ 10^{12} G is easily ruled out by the Chandra and XMM-Newton data.Comment: 5 pages, AASTeX, Revised version. Accepted for publication in ApJ Letter

A Finite Difference Representation of Neutrino Radiation Hydrodynamics in Spherically Symmetric General Relativistic Space-Time

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.Comment: reformatted to 63 pages, 24 figures, to be published in ApJ

Collapsars - Gamma-Ray Bursts and Explosions in "Failed Supernovae"

Using a two-dimensional hydrodynamics code (PROMETHEUS), we study the continued evolution of rotating massive helium stars whose iron core collapse does not produce a successful outgoing shock, but instead forms a black hole. We study the formation of a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter, alpha ~ 0.001 and 0.1. For the standard 14 solar mass model the average accretion rate for 15 s is 0.07 solar masses per second and the total energy deposited along the rotational axes by neutrino annihilation is (1 - 14) x 10**51 erg, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of this energy in the polar regions results in strong relativistic outflow - jets beamed to about 1.5% of the sky. The jets remain highly focused, and are capable of penetrating the star in 5 - 10 s. After the jet breaks through the surface of the star, highly relativistic flow can commence. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that all GRBs longer than a few seconds will make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to make in this model.Comment: Latex, 66 pages including 27 figures (9 color), Submitted to The Astrophysical Journal, latex uses aaspp4.sty. Figures also available at http://www.ucolick.org/~andre

Modeling core collapse supernovae in 2 and 3 dimensions with spectral neutrino transport

The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion.Comment: 10 pages, 3 figure

Development of the SAMEX vector magnetograph at the Marshall Space Flight Center

A breadboard design to prove the operational feasibility of SAMEX Vector Magnetograph is being developed. Although the breadboard design will not include all of the elements of the original design concept, critical elements such as the large detector array and the high resolution polarimeter will be important parts of the breadboard design to study the data analysis and compression techniques that will be needed in a SAMEX instrument, to study the calibration techniques for systemmatic errors in the polarimeter, and to obtain high resolution vector magnetograms during the next solar maximum. Although the SAMEX polarimeter is not optimum for a ground-based patrol instrument, the design concept can be confirmed with ground-based measurements and direct comparisons with the existing vector magnetograph. The extension of the scientific objectives for this breadboard design is possible if a tunable filter can be acquired

Collisional Dark Matter and the Origin of Massive Black Holes

If the cosmological dark matter is primarily in the form of an elementary particle which has cross section and mass for self-interaction having a ratio similar to that of ordinary nuclear matter, then seed black holes (formed in stellar collapse) will grow in a Hubble time, due to accretion of the dark matter, to a mass range 10^6 - 10^9 solar masses. Furthermore, the dependence of the final black hole mass on the galaxy velocity dispersion will be approximately as observed and the growth rate will show a time dependence consistent with observations. Other astrophysical consequences of collisional dark matter and tests of the idea are noted.Comment: 7 pages, no figures, LaTeX2e, Accepted for publication in Phys. Rev. Lett. Changed conten

Fabrication and characterization of Ag- and Ga-doped mesoporous glass-coated scaffolds based on natural marine sponges with improved mechanical properties

Natural marine sponges were used as sacrificial template for the fabrication of bioactive glassbased scaffolds. After sintering at 1050 ºC, the resulting samples were additionally coated with a sol silicate solution containing biologically active ions (Ag and Ga), well-known for their antibacterial properties in comparison with standard scaffolds made by PU foam templates. The produced scaffolds were characterized by superior mechanical properties (maximum compressive strength of 4 MPa) and total porosity of ~80%. Direct cell culture tests performed on the uncoated and coated samples showed positive results in terms of adhesion, proliferation, and differentiation of MC3T3-E1 cells. Moreover, vascular endothelial growth factor (VEGF) secretion from cells in contact with scaffold dissolution products was measured after 7 and 10 days of incubation, showing promising angiogenic results for bone tissue engineering applications. The antibacterial potential of the produced samples was assessed by performing agar diffusion tests against both Gram-positive and Gram-negative bacteria.EU Horizon 2020 project COACH 64255