Modeling Type Ia supernova explosions

Despite their astrophysical significanceas a major contributor to cosmic nucleosynthesis and as distance indicators in observational cosmologyType Ia supernovae lack theoretical explanation. Not only is the explosion mechanism complex due to the interaction of (potentially turbulent) hydrodynamics and nuclear reactions, but even the initial conditions for the explosion are unknown. Various progenitor scenarios have been proposed. After summarizing some general aspects of Type Ia supernova modeling, recent simulations of our group are discussed. With a sequence of modeling starting (in some cases) from the progenitor evolution and following the explosion hydrodynamics and nucleosynthesis we connect to the formation of the observables through radiation transport in the ejecta cloud. This allows us to analyze several models and to compare their outcomes with observations. While pure deflagrations of Chandrasekhar-mass white dwarfs and violent mergers of two white dwarfs lead to peculiar events (that may, however, find their correspondence in the observed sample of SNe Ia), only delayed detonations in Chandrasekhar-mass white dwarfs or sub-Chandrasekhar-mass explosions remain promising candidates for explaining normal Type Ia supernovae. © 2011 Elsevier B.V. All rights reserved.

GIANT LOBES OF CENTAURUS A RADIO GALAXY OBSERVED WITH THE SUZAKU X-RAY SATELLITE

We report on Suzaku observations of selected regions within the southern giant lobe of the radio galaxy Centaurus A. In our analysis we focus on distinct X-ray features detected with the X-ray Imaging Spectrometer within the range 0.5-10 keV, some of which are likely associated with fine structure of the lobe revealed by recent high-quality radio intensity and polarization maps. With the available photon statistics, we find that the spectral properties of the detected X-ray features are equally consistent with thermal emission from hot gas with temperatures kT > 1 keV, or with a power-law radiation continuum characterized by photon indices Gamma similar to 2.0 +/- 0.5. However, the plasma parameters implied by these different models favor a synchrotron origin for the analyzed X-ray spots, indicating that a very efficient acceleration of electrons up to greater than or similar to 10 TeV energies is taking place within the giant structure of Centaurus A, albeit only in isolated and compact regions associated with extended and highly polarized radio filaments. We also present a detailed analysis of the diffuse X-ray emission filling the whole field of view of the instrument, resulting in a tentative detection of a soft excess component best fitted by a thermal model with a temperature of kT similar to 0.5 keV. The exact origin of the observed excess remains uncertain, although energetic considerations point to thermal gas filling the bulk of the volume of the lobe and mixed with the non-thermal plasma, rather than to the alternative scenario involving a condensation of the hot intergalactic medium around the edges of the expanding radio structure. If correct, this would be the first detection of the thermal content of the extended lobes of a radio galaxy in X-rays. The corresponding number density of the thermal gas in such a case is n(g) similar to 10(-4) cm(-3), while its pressure appears to be in almost exact equipartition with the volume-averaged non-thermal pressure provided by the radio-emitting electrons and the lobes' magnetic field. A prominent large-scale fluctuation of the Galactic foreground emission, resulting in excess foreground X-ray emission aligned with the lobe, cannot be ruled out. Although tentative, our findings potentially imply that the structure of the extended lobes in active galaxies is likely to be highly inhomogeneous and non-uniform, with magnetic reconnection and turbulent acceleration processes continuously converting magnetic energy to internal energy of the plasma particles, leading to possibly significant spatial and temporal variations in the plasma beta parameter around the volume-averaged equilibrium condition beta similar to 1.

Self-organized electromagnetic field structures in laser-produced counter-streaming plasmas

Self-organization(1,2) occurs in plasmas when energy progressively transfers from smaller to larger scales in an inverse cascade(3). Global structures that emerge from turbulent plasmas can be found in the laboratory(4) and in astrophysical settings; for example, the cosmic magnetic field(5,6,) collisionless shocks in supernova remnants(7) and the internal structures of newly formed stars known as Herbig-Haro objects(8). Here we show that large, stable electromagnetic field structures can also arise within counter-streaming supersonic plasmas in the laboratory. These surprising structures, formed by a yet unexplained mechanism, are predominantly oriented transverse to the primary flow direction, extend for much larger distances than the intrinsic plasma spatial scales and persist for much longer than the plasma kinetic timescales. Our results challenge existing models of counter-streaming plasmas and can be used to better understand large-scale and long-time plasma self-organization.

COMPUTATION OF SYNTHETIC SPECTRA FROM SIMULATIONS OF RELATIVISTIC SHOCKS

Particle-in-cell (PIC) simulations of relativistic shocks are in principle capable of predicting the spectra of photons that are radiated incoherently by the accelerated particles. The most direct method evaluates the spectrum using the fields given by the Lienard-Wiechart potentials. However, for relativistic particles this procedure is computationally expensive. Here we present an alternative method that uses the concept of the photon formation length. The algorithm is suitable for evaluating spectra both from particles moving in a specific realization of a turbulent electromagnetic field or from trajectories given as a finite, discrete time series by a PIC simulation. The main advantage of the method is that it identifies the intrinsic spectral features and filters out those that are artifacts of the limited time resolution and finite duration of input trajectories.

Stochastic particle acceleration in the lobes of giant radio galaxies

We investigate the acceleration of particles by Alfven waves via the second-order Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the accelerators of ultra-high-energy cosmic rays (UHECR). We focus on the nearby Fanaroff-Riley type I radio galaxy Centaurus A. This is motivated by the coincidence of its position with the arrival direction of several of the highest energy Auger events. The conditions necessary for consistency with the acceleration time-scales predicted by quasi-linear theory are reviewed. Test particle calculations are performed in fields which guarantee electric fields with no component parallel to the local magnetic field. The results of quasi-linear theory are, to an order of magnitude, found to be accurate at low turbulence levels for non-relativistic Alfven waves and at both low and high turbulence levels in the mildly relativistic case. We conclude that for pure stochastic acceleration via Alfven waves to be plausible as the generator of UHECR in Cen A, the baryon number density would need to be several orders of magnitude below currently held upper limits.

A current-driven instability in parallel, relativistic shocks

Recently, Bell ( 2004 Mon. Not. R. Astron. Soc. 353 550) has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments by Achterberg ( 1983 Astron. Astrophys. 119 274) and McKenzie and Volk ( 1982 Astron. Astrophys. 116 191), and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is also an important issue in the problem of the formation and structure of relativistic shock fronts, particularly in relation to models of gamma-ray bursts. We have therefore generalized the linear analysis to apply to this case, assuming a relativistic background plasma and a monoenergetic, unidirectional incoming proton beam. We find essentially the same non-resonant instability observed by Bell and show that also, under GRB conditions, it grows much faster than the resonant waves. We quantify the extent to which thermal effects in the background plasma limit the maximum growth rate.

STEADY-STATE SOLUTIONS IN NONLINEAR DIFFUSIVE SHOCK ACCELERATION

Stationary solutions to the equations of nonlinear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean free path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cutoff in energy space-particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for the substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven nonresonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off-a diagnostic that is in principle accessible to observation.

Dynamical symmetry breaking with optimal control: Reducing the number of pieces

We analyze the production of defects during the dynamical crossing of a mean-field phase transition with a real order parameter. When the parameter that brings the system across the critical point changes in time according to a power-law schedule, we recover the predictions dictated by the well-known Kibble-Zurek theory. For a fixed duration of the evolution, we show that the average number of defects can be drastically reduced for a very large but finite system, by optimizing the time dependence of the driving using optimal control techniques. Furthermore, the optimized protocol is robust against small fluctuations.

Multi-Dimensional Simulations of the Expanding Supernova Remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of 10 M-circle dot and an energy of 1.5 x 10(44) J. A termination shock in the progenitor's stellar wind at a distance of 0 ''.43-0 ''.51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0 ''.63 +/- 0 ''.01 and maximum particle number density of (7.11 +/- 1.78) x 10(7) m(-3) produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

MarcoPolo-R near earth asteroid sample return mission

MarcoPolo-R is a sample return mission to a primitive Near-Earth Asteroid (NEA) proposed in collaboration with NASA. It will rendezvous with a primitive NEA, scientifically characterize it at multiple scales,and return a unique sample to Earth unaltered by the atmospheric entry process or terrestrial weathering. MarcoPolo-R will return bulk samples (up to 2 kg) from an organic-rich binary asteroid to Earth for laboratory analyses, allowing us to: explore the origin of planetary materials and initial stages of habitable planet formation; identify and characterize the organics and volatiles in a primitive asteroid; understand the unique geomorphology, dynamics and evolution of a binaryNEA. This project is based on the previous Marco Polo mission study,which was selected for the Assessment Phase of the first round of Cosmic Vision. Its scientific rationale was highly ranked by ESA committees andit was not selected only because the estimated cost was higher than theallotted amount for an M class mission. The cost of Marco Polo-R will be reduced to within the ESA medium mission budget by collaboration withAPL (John Hopkins University) and JPL in the NASA program for coordination with ESA's Cosmic Vision Call. The baseline target is a binary asteroid (175706) 1996 FG3, which offers a very efficient operational and technical mission profile. A binary target also providesenhanced science return. The choice of this target will allow newinvestigations to be performed more easily than at a single object, andalso enables investigations of the fascinating geology and geophysics ofasteroids that are impossible at a single object. Several launch windows have been identified in the time-span 2020-2024. A number of otherpossible primitive single targets of high scientific interest have beenidentified covering a wide range of possible launch dates. The baselinemission scenario of Marco Polo-R to 1996 FG3 is as follows: a singleprimary spacecraft provided by ESA, carrying the Earth Re-entry Capsule, sample acquisition and transfer system provided by NASA, will be launched by a Soyuz-Fregat rocket from Kourou into GTO and using two space segment stages. Two similar missions with two launch windows, in 2021 and 2022 and for both sample return in 2029 (with mission durationof 7 and 8 years), have been defined. Earlier or later launches, in 2020 or 2024, also offer good opportunities. All manoeuvres are carried out by a chemical propulsion system. MarcoPolo-R takes advantage of three industrial studies completed as part of the previous Marco Polo mission (see ESA/SRE (2009)3, Marco Polo Yellow Book) and of the expertise of the consortium led by Dr. A.F. Cheng (PI of the NASA NEAR Shoemaker mission) of the JHU-APL, including JPL, NASA ARC, NASA LaRC, and MIT.

Sputtering of oxygen ice by low energy ions

Naturally occurring ices lie on both interstellar dust grains and on celestial objects, such as those in the outer Solar system. These ices are continuously subjected to irradiation by ions from the solar wind and/or cosmic rays, which modify their surfaces. As a result, new molecular species may form which can be sputtered off into space or planetary atmospheres. We determined the experimental values of sputtering yields for irradiation of oxygen ice at 10 K by singly (He+, C+, N+, O+ and Ar+) and doubly (C2 +, N2 + and O2 +) charged ions with 4 keV kinetic energy. In these laboratory experiments, oxygen ice was deposited and irradiated by ions in an ultra high vacuum chamber at low temperature to simulate the environment of space. The number of molecules removed by sputtering was observed by measurement of the ice thickness using laser interferometry. Preliminary mass spectra were taken of sputtered species and of molecules formed in the ice by temperature programmed desorption (TPD). We find that the experimental sputtering yields increase approximately linearly with the projectile ion mass (or momentum squared) for all ions studied. No difference was found between the sputtering yields for singly and doubly charged ions of the same atom within the experimental uncertainty, as expected for a process dominated by momentum transfer. The experimental sputter yields are in good agreement with values calculated using a theoretical model except in the case of oxygen ions. Preliminary studies have shown molecular oxygen as the dominant species sputtered and TPD measurements indicate ozone formation.

Turbulent amplification of magnetic fields in laboratory laser-produced shock waves

X-ray and radio observations of the supernova remnant Cassiopeia A reveal the presence of magnetic fields about 100 times stronger than those in the surrounding interstellar medium. Field coincident with the outer shock probably arises through a nonlinear feedback process involving cosmic rays. The origin of the large magnetic field in the interior of the remnant is less clear but it is presumably stretched and amplified by turbulent motions. Turbulence may be generated by hydrodynamic instability at the contact discontinuity between the supernova ejecta and the circumstellar gas9. However, optical observations of Cassiopeia A indicate that the ejecta are interacting with a highly inhomogeneous, dense circumstellar cloud bank formed before the supernova explosion. Here we investigate the possibility that turbulent amplification is induced when the outer shock overtakes dense clumps in the ambient medium. We report laboratory experiments that indicate the magnetic field is amplified when the shock interacts with a plastic grid. We show that our experimental results can explain the observed synchrotron emission in the interior of the remnant. The experiment also provides a laboratory example of magnetic field amplification by turbulence in plasmas, a physical process thought to occur in many astrophysical phenomena.

Astrochemistry

In the last 40 years a wide range of molecules, including neutrals, cations and anions, containing up to 13 atoms—in addition to detections of C₆₀ and C₇₀ — have been found in the harsh environment of the interstellar medium. The exquisite sensitivity and very high spectral and, more recently, spatial resolution, of modern telescopes has enabled the physics of star formation to be probed through rotational line emission. In this article, I review the basic properties of interstellar clouds and the processes that initiate the chemistry and generate chemical complexity, particularly in regions of star and planet formation. Our understanding of astrochemistry has evolved over the years. Before 1990, the general consensus was that molecules were formed in binary, gas-phase, or volume, reactions, most importantly ion-neutral reactions despite the very low ionization in clouds. Since then, observations have indicated unambiguously that there is also a contribution from surface processes, particularly on the icy mantles that form around refractory grain cores in cold, dense gas. The balance between these two processes depends on particular physical conditions and can vary during the life cycle of a particular volume of interstellar cloud.The complex chemistry that occurs in space is driven mostly through interaction of thegas with cosmic ray protons, a source of ionization that enables a rich ion-neutral chemistry. In addition, I show that the interaction between the gas and the dust in cold, dense regionsalso leads to additional chemical complexity through reactions that take place in ices at onlya few tens of degrees above absolute zero. Although densities are low compared to those in terrestrial environments, the extremely long life times of interstellar clouds and their enormous sizes, enable complex molecules to be synthesised and detected. I show that in some instances, particularly in reactions involving deuterium, the rotational populations of reactants, together with spin-selection rules, can determine the detailed abundances. Although the review is mainly focused on regions associated with star formation, I also consider chemistry in other interesting astronomical regions — in the early Universe and in the envelopes formed by mass loss during the final stages of stellar evolution.

Systematic Uncertainties Associated with the Cosmological Analysis of the First PAN-STARRS1 Type Ia Supernova Sample

We probe the systematic uncertainties from the 113 Type Ia supernovae (SN Ia) in the Pan-STARRS1 (PS1) sample along with 197 SN Ia from a combination of low-redshift surveys. The companion paper by Rest et al. describes the photometric measurements and cosmological inferences from the PS1 sample. The largest systematic uncertainty stems from the photometric calibration of the PS1 and low-z samples. We increase the sample of observed Calspec standards from 7 to 10 used to define the PS1 calibration system. The PS1 and SDSS-II calibration systems are compared and discrepancies up to ∼0.02 mag are recovered. We find uncertainties in the proper way to treat intrinsic colors and reddening produce differences in the recovered value of w up to 3%. We estimate masses of host galaxies of PS1 supernovae and detect an insignificant difference in distance residuals of the full sample of 0.037 ± 0.031 mag for host galaxies with high and low masses. Assuming flatness and including systematic uncertainties in our analysis of only SNe measurements, we find w = -1.120^+0.360_-0.206(Stat)^+0.269_-0.291(Sys). With additional constraints from Baryon acoustic oscillation, cosmic microwave background (CMB) (Planck) and H₀ measurements, we find w = -1.166^+0.072_-0.069 and Ω_m = 0.280^+0.013_-0.012 (statistical and systematic errors added in quadrature). The significance of the inconsistency with w = -1 depends on whether we use Planck or Wilkinson Microwave Anisotropy Probe measurements of the CMB: w_{BAO+H0+SN+WMAP} = -1.124^+0.083_-0.065.

Cosmological Constraints from Measurements of Type Ia Supernovae Discovered During the First 1.5 Yr of the Pan-STARRS1 Survey

We present griz_P1 light curves of 146 spectroscopically confirmed Type Ia supernovae (SNe Ia; 0.03 < z < 0.65) discovered during the first 1.5 yr of the Pan-STARRS1 Medium Deep Survey. The Pan-STARRS1 natural photometric system is determined by a combination of on-site measurements of the instrument response function and observations of spectrophotometric standard stars. We find that the systematic uncertainties in the photometric system are currently 1.2% without accounting for the uncertainty in the Hubble Space Telescope Calspec definition of the AB system. A Hubble diagram is constructed with a subset of 113 out of 146 SNe Ia that pass our light curve quality cuts. The cosmological fit to 310 SNe Ia (113 PS1 SNe Ia + 222 light curves from 197 low-z SNe Ia), using only supernovae (SNe) and assuming a constant dark energy equation of state and flatness, yields w = -1.120^+0.360_-0.206(Stat)^+0.269_0.291(Sys). When combined with BAO+CMB(Planck)+H₀, the analysis yields Ω_M = 0.280^+0.013_0.012 and w = -1.166^+0.072_-0.069 including all identified systematics. The value of w is inconsistent with the cosmological constant value of -1 at the 2.3σ level. Tension endures after removing either the baryon acoustic oscillation (BAO) or the H₀ constraint, though it is strongest when including the H₀ constraint. If we include WMAP9 cosmic microwave background (CMB) constraints instead of those from Planck, we find w = -1.124^+0.083_-0.065, which diminishes the discord to <2σ. We cannot conclude whether the tension with flat ΛCDM is a feature of dark energy, new physics, or a combination of chance and systematic errors. The full Pan-STARRS1 SN sample with ∼three times as many SNe should provide more conclusive results.