-
Insensitivity of a turbulent laser-plasma dynamo to initial conditions
Authors:
A. F. A. Bott,
L. Chen,
P. Tzeferacos,
C. A. J. Palmer,
A. R. Bell,
R. Bingham,
A. Birkel,
D. H. Froula,
J. Katz,
M. W. Kunz,
C. -K. Li,
H-S. Park,
R. Petrasso,
J. S. Ross,
B. Reville,
D. Ryu,
F. H. Séguin,
T. G. White,
A. A. Schekochihin,
D. Q. Lamb,
G. Gregori
Abstract:
It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetised laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one aris…
▽ More
It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetised laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one arising from two pre-magnetised plasma jets whose creation is identical save for the addition of a strong external magnetic field imposed by a pulsed magnetic field generator (`MIFEDS'). We investigate the differences between the two turbulent systems using a Thomson-scattering diagnostic, X-ray self-emission imaging and proton radiography. The Thomson-scattering spectra and X-ray images suggest that the presence of the external magnetic field has a limited effect on the plasma dynamics in the experiment. While the presence of the external magnetic field induces collimation of the flows in the colliding plasma jets and the initial strengths of the magnetic fields arising from the interaction between the colliding jets are significantly larger as a result of the external field, the energy and morphology of the stochastic magnetic fields post-amplification are indistinguishable. We conclude that, for turbulent laser-plasmas with super-critical magnetic Reynolds numbers, the dynamo-amplified magnetic fields are determined by the turbulent dynamics rather than the seed fields and modest changes in the initial flow dynamics of the plasma, a finding consistent with theoretical expectations and simulations of turbulent dynamos.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Strong suppression of heat conduction in a laboratory replica of galaxy-cluster turbulent plasmas
Authors:
J. Meinecke,
P. Tzeferacos,
J. S. Ross,
A. F. A. Bott,
S. Feister,
H. -S. Park,
A. R. Bell,
R. Blandford,
R. L. Berger,
R. Bingham,
A. Casner,
L. E. Chen,
J. Foster,
D. H. Froula,
C. Goyon,
D. Kalantar,
M. Koenig,
B. Lahmann,
C. -K. Li,
Y. Lu,
C. A. J. Palmer,
R. Petrasso,
H. Poole,
B. Remington,
B. Reville
, et al. (10 additional authors not shown)
Abstract:
Galaxy clusters are filled with hot, diffuse X-ray emitting plasma, with a stochastically tangled magnetic field whose energy is close to equipartition with the energy of the turbulent motions \cite{zweibel1997, Vacca}. In the cluster cores, the temperatures remain anomalously high compared to what might be expected considering that the radiative cooling time is short relative to the Hubble time \…
▽ More
Galaxy clusters are filled with hot, diffuse X-ray emitting plasma, with a stochastically tangled magnetic field whose energy is close to equipartition with the energy of the turbulent motions \cite{zweibel1997, Vacca}. In the cluster cores, the temperatures remain anomalously high compared to what might be expected considering that the radiative cooling time is short relative to the Hubble time \cite{cowie1977,fabian1994}. While feedback from the central active galactic nuclei (AGN) \cite{fabian2012,birzan2012,churazov2000} is believed to provide most of the heating, there has been a long debate as to whether conduction of heat from the bulk to the core can help the core to reach the observed temperatures \cite{narayan2001,ruszkowski2002,kunz2011}, given the presence of tangled magnetic fields. Interestingly, evidence of very sharp temperature gradients in structures like cold fronts implies a high degree of suppression of thermal conduction \cite{markevitch2007}. To address the problem of thermal conduction in a magnetized and turbulent plasma, we have created a replica of such a system in a laser laboratory experiment. Our data show a reduction of local heat transport by two orders of magnitude or more, leading to strong temperature variations on small spatial scales, as is seen in cluster plasmas \cite{markevitch2003}.
△ Less
Submitted 18 May, 2021;
originally announced May 2021.
-
Inefficient magnetic-field amplification in supersonic laser-plasma turbulence
Authors:
A. F. A. Bott,
L. Chen,
G. Boutoux,
T. Caillaud,
A. Duval,
M. Koenig,
B. Khiar,
I. Lantuéjoul,
L. Le-Deroff,
B. Reville,
R. Rosch,
D. Ryu,
C. Spindloe,
B. Vauzour,
B. Villette,
A. A. Schekochihin,
D. Q. Lamb,
P. Tzeferacos,
G. Gregori,
A. Casner
Abstract:
We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic plasma with a large magnetic Reynolds number ($\mathrm{Rm} \approx 45$) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller t…
▽ More
We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic plasma with a large magnetic Reynolds number ($\mathrm{Rm} \approx 45$) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller than the outer scale of the turbulent cascade was also observed. Our results support the notion that moderately supersonic, low-magnetic-Prandtl-number plasma turbulence is inefficient at amplifying magnetic fields.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
Time-resolved fast turbulent dynamo in a laser plasma
Authors:
A. F. A. Bott,
P. Tzeferacos,
L. Chen,
C. A. J. Palmer,
A. Rigby,
A. Bell,
R. Bingham,
A. Birkel,
C. Graziani,
D. H. Froula,
J. Katz,
M. Koenig,
M. W. Kunz,
C. K. Li,
J. Meinecke,
F. Miniati,
R. Petrasso,
H. -S. Park,
B. A. Remington,
B. Reville,
J. S. Ross,
D. Ryu,
D. Ryutov,
F. Séguin,
T. G. White
, et al. (3 additional authors not shown)
Abstract:
Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas ($\mathrm{Pm} < 1$). However, the…
▽ More
Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas ($\mathrm{Pm} < 1$). However, the same framework proposes that the fluctuation dynamo should operate differently when $\mathrm{Pm} \gtrsim 1$, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports a new experiment that creates a laboratory $\mathrm{Pm} \gtrsim 1$ plasma dynamo for the first time. We provide a time-resolved characterization of the plasma's evolution, measuring temperatures, densities, flow velocities and magnetic fields, which allows us to explore various stages of the fluctuation dynamo's operation. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude from its initial value and saturate dynamically. It is shown that the growth of these fields occurs exponentially at a rate that is much greater than the turnover rate of the driving-scale stochastic motions. Our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized MHD simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.
△ Less
Submitted 24 July, 2020;
originally announced July 2020.
-
Code-to-Code Comparison and Validation of the Radiation-Hydrodynamics Capabilities of the FLASH Code Using a Laboratory Astrophysical Jet
Authors:
Chris Orban,
Milad Fatenejad,
Don Q. Lamb
Abstract:
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available…
▽ More
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. (2008). That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. (2008) experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equation of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. (2008) experiment and discuss its relation to the Wan et al. (1997) experiment at the NOVA laser.
△ Less
Submitted 18 April, 2022; v1 submitted 17 June, 2020;
originally announced June 2020.
-
Retrieving fields from proton radiography without source profiles
Authors:
M. F. Kasim,
A. F. A. Bott,
P. Tzeferacos,
D. Q. Lamb,
G. Gregori,
S. M. Vinko
Abstract:
Proton radiography is a technique in high energy density science to diagnose magnetic and/or electric fields in a plasma by firing a proton beam and detecting its modulated intensity profile on a screen. Current approaches to retrieve the integrated field from the modulated intensity profile require the unmodulated beam intensity profile before the interaction, which is rarely available experiment…
▽ More
Proton radiography is a technique in high energy density science to diagnose magnetic and/or electric fields in a plasma by firing a proton beam and detecting its modulated intensity profile on a screen. Current approaches to retrieve the integrated field from the modulated intensity profile require the unmodulated beam intensity profile before the interaction, which is rarely available experimentally due to shot-to-shot variability. In this paper, we present a statistical method to retrieve the integrated field without needing to know the exact source profile. We apply our method to experimental data, showing the robustness of our approach. Our proposed technique allows not only for the retrieval of the path-integrated fields, but also of the statistical properties of the fields.
△ Less
Submitted 13 September, 2019; v1 submitted 30 May, 2019;
originally announced May 2019.
-
Transport of high-energy charged particles through spatially-intermittent turbulent magnetic fields
Authors:
L. E. Chen,
A. F. A. Bott,
P. Tzeferacos,
A. Rigby,
A. Bell,
R. Bingham,
C. Graziani,
J. Katz,
M. Koenig,
C. K. Li,
R. Petrasso,
H. -S. Park,
J. S. Ross,
D. Ryu,
T. G. White,
B. Reville,
J. Matthews,
J. Meinecke,
F. Miniati,
E. G. Zweibel,
S. Sarkar,
A. A. Schekochihin,
D. Q. Lamb,
D. H. Froula,
G. Gregori
Abstract:
Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results sho…
▽ More
Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.
△ Less
Submitted 30 March, 2020; v1 submitted 13 August, 2018;
originally announced August 2018.
-
Proton imaging of stochastic magnetic fields
Authors:
A. F. A. Bott,
C. Graziani,
P. Tzeferacos,
T. G. White,
D. Q. Lamb,
G. Gregori,
A. A. Schekochihin
Abstract:
Recent laser-plasma experiments report the existence of dynamically significant magnetic fields, whose statistical characterisation is essential for understanding the physical processes these experiments are attempting to investigate. In this paper, we show how a proton imaging diagnostic can be used to determine a range of relevant magnetic field statistics, including the magnetic-energy spectrum…
▽ More
Recent laser-plasma experiments report the existence of dynamically significant magnetic fields, whose statistical characterisation is essential for understanding the physical processes these experiments are attempting to investigate. In this paper, we show how a proton imaging diagnostic can be used to determine a range of relevant magnetic field statistics, including the magnetic-energy spectrum. To achieve this goal, we explore the properties of an analytic relation between a stochastic magnetic field and the image-flux distribution created upon imaging that field. We conclude that features of the beam's final image-flux distribution often display a universal character determined by a single, field-scale dependent parameter - the contrast parameter - which quantifies the relative size of the correlation length of the stochastic field, proton displacements due to magnetic deflections, and the image magnification. For stochastic magnetic fields, we establish the existence of four contrast regimes - linear, nonlinear injective, caustic and diffusive - under which proton-flux images relate to their parent fields in a qualitatively distinct manner. As a consequence, it is demonstrated that in the linear or nonlinear injective regimes, the path-integrated magnetic field experienced by the beam can be extracted uniquely, as can the magnetic-energy spectrum under a further statistical assumption of isotropy. This is no longer the case in the caustic or diffusive regimes. We also discuss complications to the contrast-regime characterisation arising for inhomogeneous, multi-scale stochastic fields, as well as limitations currently placed by experimental capabilities on extracting magnetic field statistics. The results presented in this paper provide a comprehensive description of proton images of stochastic magnetic fields, with applications for improved analysis of given proton-flux images.
△ Less
Submitted 5 August, 2017;
originally announced August 2017.
-
Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma
Authors:
P. Tzeferacos,
A. Rigby,
A. Bott,
A. R. Bell,
R. Bingham,
A. Casner,
F. Cattaneo,
E. M. Churazov,
J. Emig,
F. Fiuza,
C. B. Forest,
J. Foster,
C. Graziani,
J. Katz,
M. Koenig,
C. -K. Li,
J. Meinecke,
R. Petrasso,
H. -S. Park,
B. A. Remington,
J. S. Ross,
D. Ryu,
D. Ryutov,
T. G. White,
B. Reville
, et al. (5 additional authors not shown)
Abstract:
Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $μ$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are…
▽ More
Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $μ$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.
△ Less
Submitted 9 February, 2017;
originally announced February 2017.
-
Numerical modeling of laser-driven experiments aiming to demonstrate magnetic field amplification via turbulent dynamo
Authors:
P. Tzeferacos,
A. Rigby,
A. Bott,
A. R. Bell,
R. Bingham,
A. Casner,
F. Cattaneo,
E. M. Churazov,
J. Emig,
N. Flocke,
F. Fiuza,
C. B. Forest,
J. Foster,
C. Graziani,
J. Katz,
M. Koenig,
C. -K. Li,
J. Meinecke,
R. Petrasso,
H. -S. Park,
B. A. Remington,
J. S. Ross,
D. Ryu,
D. Ryutov,
K. Weide
, et al. (7 additional authors not shown)
Abstract:
The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and…
▽ More
The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.
△ Less
Submitted 9 February, 2017;
originally announced February 2017.
-
Inferring Morphology and Strength of Magnetic Fields From Proton Radiographs
Authors:
Carlo Graziani,
Petros Tzeferacos,
Donald Q. Lamb,
Chikang Li
Abstract:
Proton radiography is an important diagnostic method for laser plasma experiments, and is particularly important in the analysis of magnetized plasmas. The theory of radiographic image analysis has heretofore only permitted somewhat limited analysis of the radiographs of such plasmas. We furnish here a theory that remedies this deficiency. We show that to linear order in magnetic field gradients,…
▽ More
Proton radiography is an important diagnostic method for laser plasma experiments, and is particularly important in the analysis of magnetized plasmas. The theory of radiographic image analysis has heretofore only permitted somewhat limited analysis of the radiographs of such plasmas. We furnish here a theory that remedies this deficiency. We show that to linear order in magnetic field gradients, proton radiographs are projection images of the MHD current along the proton trajectories. We demonstrate that in the linear approximation, the full structure of the perpedicular magnetic field can be reconstructed by solving a steady-state inhomogeneous 2-dimensional diffusion equation sourced by the radiograph fluence contrast data. We explore limitations of the inversion method due to Poisson noise, to discretization errors, to radiograph edge effects, and to obstruction by laser target structures. We also provide a separate analysis that is well-suited to the inference of isotropic-homogeneous magnetic turbulence spectra. We discuss extension of these results to the nonlinear contrast regime.
△ Less
Submitted 28 March, 2016;
originally announced March 2016.
-
Creation of Magnetized Jet Using a Ring of Laser Beams
Authors:
Wen Fu,
Edison P. Liang,
Petros Tzeferacos,
Donald Q. Lamb
Abstract:
We propose a new way of generating magnetized supersonic jets using a ring laser to irradiate a flat surface target. Using 2D FLASH code simulations which include the Biermann Battery term, we demonstrate that strong toroidal fields can be generated and sustained downstream in the collimated jet outflow far from the target surface. The field strength can be controlled by varying the ring laser sep…
▽ More
We propose a new way of generating magnetized supersonic jets using a ring laser to irradiate a flat surface target. Using 2D FLASH code simulations which include the Biermann Battery term, we demonstrate that strong toroidal fields can be generated and sustained downstream in the collimated jet outflow far from the target surface. The field strength can be controlled by varying the ring laser separation, thereby providing a versatile laboratory platform for studying the effects of magnetic field in a variety of astrophysical settings.
△ Less
Submitted 10 December, 2015; v1 submitted 15 September, 2014;
originally announced September 2014.
-
The Biermann Catastrophe in Numerical MHD
Authors:
Carlo Graziani,
Petros Tzeferacos,
Dongwook Lee,
Donald Q. Lamb,
Klaus Weide,
Milad Fatenejad,
Joshua Miller
Abstract:
The Biermann Battery effect is frequently invoked in cosmic magnetogenesis and studied in High-Energy Density laboratory physics experiments. Generation of magnetic fields by the Biermann effect due to mis-aligned density and temperature gradients in smooth flow <i>behind</i> shocks is well known. We show that a Biermann-effect magnetic field is also generated <i>within</i> shocks. Direct implemen…
▽ More
The Biermann Battery effect is frequently invoked in cosmic magnetogenesis and studied in High-Energy Density laboratory physics experiments. Generation of magnetic fields by the Biermann effect due to mis-aligned density and temperature gradients in smooth flow <i>behind</i> shocks is well known. We show that a Biermann-effect magnetic field is also generated <i>within</i> shocks. Direct implementation of the Biermann effect in MHD codes does not capture this physical process, and worse, produces unphysical magnetic fields at shocks whose value does not converge with resolution. We show that this convergence breakdown is due to naive discretization, which fails to account for the fact that discretized irrotational vector fields have spurious solenoidal components that grow without bound near a discontinuity. We show that careful consideration of the kinetics of ion viscous shocks leads to a formulation of the Biermann effect that gives rise to a convergent algorithm. We note two novel physical effects: a <i>resistive magnetic precursor</i> in which Biermann-generated field in the shock "leaks" resistively upstream; and a <i>thermal magnetic precursor</i>, in which field is generated by the Biermann effect ahead of the shock front due to gradients created by the shock's electron thermal conduction precursor. Both effects appear to be potentially observable in experiments at laser facilities. We re-examine published studies of magnetogenesis in galaxy cluster formation, and conclude that the simulations in question had inadequate resolution to reliably estimate the field generation rate. Corrected estimates suggest primordial field values in the range $B\sim 10^{-22}$G --- $10^{-19}$G by $z=3$.
△ Less
Submitted 22 January, 2015; v1 submitted 18 August, 2014;
originally announced August 2014.
-
Experimental demonstration of an inertial collimation mechanism in nested outflows
Authors:
R. Yurchak,
A. Ravasio,
A. Pelka,
S. Pikuz Jr.,
E. Falize,
T. Vinci,
M. Koenig,
B. Loupias,
A. Benuzzi-Mounaix,
M. Fatenejad,
P. Tzeferacos,
D. Q. Lamb,
E. G. Blackman
Abstract:
Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-ene…
▽ More
Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-energy laser beam onto a solid target. A nested geometry was created by shaping the energy distribution at the focal spot with a dedicated phase plate. Optical and X-ray diagnostics were used to study the interacting flows. Experimental results and numerical hydrodynamic simulations indeed show the formation of strongly collimated jets. Our work experimentally confirms the "shock-focused inertial confinement" mechanism proposed in previous theoretical astrophysics investigations.
△ Less
Submitted 19 March, 2014;
originally announced March 2014.
-
Software Abstractions and Methodologies for HPC Simulation Codes on Future Architectures
Authors:
A. Dubey,
S. Brandt,
R. Brower,
M. Giles,
P. Hovland,
D. Q. Lamb,
F. Loffler,
B. Norris,
B. OShea,
C. Rebbi,
M. Snir,
R. Thakur
Abstract:
Large, complex, multi-scale, multi-physics simulation codes, running on high performance com-puting (HPC) platforms, have become essential to advancing science and engineering. These codes simulate multi-scale, multi-physics phenomena with unprecedented fidelity on petascale platforms, and are used by large communities. Continued ability of these codes to run on future platforms is as crucial to t…
▽ More
Large, complex, multi-scale, multi-physics simulation codes, running on high performance com-puting (HPC) platforms, have become essential to advancing science and engineering. These codes simulate multi-scale, multi-physics phenomena with unprecedented fidelity on petascale platforms, and are used by large communities. Continued ability of these codes to run on future platforms is as crucial to their communities as continued improvements in instruments and facilities are to experimental scientists. However, the ability of code developers to do these things faces a serious challenge with the paradigm shift underway in platform architecture. The complexity and uncertainty of the future platforms makes it essential to approach this challenge cooperatively as a community. We need to develop common abstractions, frameworks, programming models and software development methodologies that can be applied across a broad range of complex simulation codes, and common software infrastructure to support them. In this position paper we express and discuss our belief that such an infrastructure is critical to the deployment of existing and new large, multi-scale, multi-physics codes on future HPC platforms.
△ Less
Submitted 6 September, 2013;
originally announced September 2013.
-
Three-dimensional Simulations of Pure Deflagration Models for Thermonuclear Supernovae
Authors:
Min Long,
George C. Jordan IV,
Daniel R. van Rossum,
Benedikt Diemer,
Carlo Graziani,
Richard Kessler,
Bradley Meyer,
Paul Rich,
Don Q. Lamb
Abstract:
We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 d…
▽ More
We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy $E_{\mathrm{nuc}}$, have larger kinetic energies $E_{\mathrm{K}}$, and produce more $^{56}$Ni than those with many ignition points, and differ in the distribution of $^{56}$Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes $M_\mathrm{B}$ and light curves that rise and decline more quickly; their $M_\mathrm{B}$ and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.
△ Less
Submitted 13 May, 2014; v1 submitted 31 July, 2013;
originally announced July 2013.
-
Radiation Transport for Explosive Outflows: A Multigroup Hybrid Monte Carlo Method
Authors:
Ryan T. Wollaeger,
Daniel R. van Rossum,
Carlo Graziani,
Sean M. Couch,
George C. Jordan IV,
Donald Q. Lamb,
Gregory A. Moses
Abstract:
We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC…
▽ More
We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.
△ Less
Submitted 21 October, 2013; v1 submitted 24 June, 2013;
originally announced June 2013.
-
A Radiation-Hydrodynamics Code Comparison for Laser-Produced Plasmas: FLASH versus HYDRA and the Results of Validation Experiments
Authors:
Chris Orban,
Milad Fatenejad,
Sugreev Chawla,
Scott C. Wilks,
Donald Q. Lamb
Abstract:
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and expansion of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities recently added to the pu…
▽ More
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and expansion of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities recently added to the publicly available FLASH radiation-hydrodynamics code, we present detailed comparisons of FLASH results to new and previously published results from the HYDRA code used extensively at Lawrence Livermore National Laboratory. We focus on two very different problems of interest: (1) an Aluminum slab irradiated by 15.3 and 76.7 mJ of "pre-pulse" laser energy and (2) a mm-long triangular groove cut in an Aluminum target irradiated by a rectangular laser beam. Because this latter problem bears a resemblance to astrophysical jets, Grava et al., Phys. Rev. E, 78, (2008) performed this experiment and compared detailed x-ray interferometric measurements of electron number densities to HYDRA simulations. Thus, the former problem provides an opportunity for code-to-code comparison, while the latter provides an opportunity for both code-to-code comparison and validation. Despite radically different schemes for determining the computational mesh, and different equation of state and opacity models, the HYDRA and FLASH codes yield results that are in excellent agreement for both problems and with the experimental data for the latter. Having validated the FLASH code in this way, we use the code to further investigate the formation of the jet seen in the Grava et al. (2008) experiment and discuss its relation to the Wan et al. (1997) experiment at the NOVA laser.
△ Less
Submitted 6 June, 2013;
originally announced June 2013.
-
The Effect of the Pre-Detonation Stellar Internal Velocity Profile on the Nucleosynthetic Yields in Type Ia Supernova
Authors:
Yeunjin Kim,
G. C. Jordan IV,
Carlo Graziani,
B. S. Meyer,
D. Q. Lamb,
J. W. Truran
Abstract:
A common model of the explosion mechanism of Type Ia supernovae is based on a delayed detonation of a white dwarf. A variety of models differ primarily in the method by which the deflagration leads to a detonation. A common feature of the models, however, is that all of them involve the propagation of the detonation through a white dwarf that is either expanding or contracting, where the stellar i…
▽ More
A common model of the explosion mechanism of Type Ia supernovae is based on a delayed detonation of a white dwarf. A variety of models differ primarily in the method by which the deflagration leads to a detonation. A common feature of the models, however, is that all of them involve the propagation of the detonation through a white dwarf that is either expanding or contracting, where the stellar internal velocity profile depends on both time and space. In this work, we investigate the effects of the pre-detonation stellar internal velocity profile and the post-detonation velocity of expansion on the production of alpha-particle nuclei, including Ni56, which are the primary nuclei produced by the detonation wave. We perform one-dimensional hydrodynamic simulations of the explosion phase of the white dwarf for center and off-center detonations with five different stellar velocity profiles at the onset of the detonation. We observe two distinct post-detonation expansion phases: rarefaction and bulk expansion. Almost all the burning to Ni56 occurs only in the rarefaction phase, and its expansion time scale is influenced by pre-existing flow structure in the star, in particular by the pre-detonation stellar velocity profile. We find that the mass fractions of the alpha-particle nuclei, including Ni56, are tight functions of the empirical physical parameter rho_up/v_down, where rho_up is the mass density immediately upstream of the detonation wave front and v_down is the velocity of the flow immediately downstream of the detonation wave front. We also find that v_down depends on the pre-detonation flow velocity. We conclude that the properties of the pre-existing flow, in particular the internal stellar velocity profile, influence the final isotopic composition of burned matter produced by the detonation.
△ Less
Submitted 6 June, 2013;
originally announced June 2013.
-
Comparing the light curves of simulated Type Ia Supernovae with observations using data-driven models
Authors:
Benedikt Diemer,
Richard Kessler,
Carlo Graziani,
George C. Jordan IV,
Donald Q. Lamb,
Min Long,
Daniel R. van Rossum
Abstract:
We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a…
▽ More
We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a goodness-of-fit parameter. Each fit is performed multiple times with different choices of filter bands and epoch range in order to quantify the systematic uncertainty on the fitted parameters. We use a parametric population model for the distribution of observed SNIa parameters from large surveys, and extend it to represent red, dim, and bright outliers found in a low-redshift SNIa data set. We discuss the potential uncertainties of this population model and find it to be reliable given the current uncertainties on cosmological parameters. Using our population model, we assign each set of fitted parameters a likelihood of being observed in nature, and a figure of merit based on this likelihood. We define a second figure of merit based on the quality of the light curve fit, and combine the two measures into an overall figure of merit for each explosion model. We compute figures of merit for a variety of 1D, 2D and 3D explosion models and show that our evaluation method allows meaningful inferences across a wide range of light curve quality and fitted parameters.
△ Less
Submitted 2 August, 2013; v1 submitted 5 March, 2013;
originally announced March 2013.
-
Increase of the Density, Temperature and Velocity of Plasma Jets driven by a Ring of High Energy Laser Beams
Authors:
Wen Fu,
Edison P. Liang,
Milad Fatenejad,
Donald Q. Lamb,
Michael Grosskopf,
Hye-Sook Park,
Bruce Remington,
Anatoly Spitkovsky
Abstract:
Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platfor…
▽ More
Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic simulations, we demonstrate that the collimated jets can achieve much higher density, temperature and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, instead of focused onto a single spot. We carried out simulations with different ring radii and studied their effects on the jet properties. Implications for laboratory collisionless shock experiments are discussed.
△ Less
Submitted 15 April, 2013; v1 submitted 10 September, 2012;
originally announced September 2012.
-
The Detonation Mechanism of the Pulsationally-Assisted Gravitationally-Confined Detonation Model of Type Ia Supernovae
Authors:
G. C. Jordan IV,
C. Graziani,
R. T. Fisher,
D. M. Townsley,
C. Meakin,
K. Weide,
L. B. Reid,
J. Norris,
R. Hudson,
D. Q. Lamb
Abstract:
We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phas…
▽ More
We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final kinetic energy and nickel-56 yields conform better to observational values than is the case for the "classical" GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel), the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that are necessary for the Zel'dovich gradient mechanism to produce a detonation. We demonstrate feasibility of this scenario with three 3-dimensional (3D), full star simulations of this model using the FLASH code. We characterized the simulations by the energy released during the deflagration phase, which ranged from 38% to 78% of the white dwarf's binding energy. We show that the necessary conditions for detonation are achieved in all three of the models.
△ Less
Submitted 17 September, 2012; v1 submitted 17 February, 2012;
originally announced February 2012.
-
Spectral Lag Relations in GRB Pulses Detected with HETE-2
Authors:
M. Arimoto,
N. Kawai,
K. Asano,
K. Hurley,
M. Suzuki,
Y. E. Nakagawa,
T. Shimokawabe,
N. V. Pazmino,
R. Sato,
M. Matsuoka,
A. Yoshida,
T. Tamagawa,
Y. Shirasaki,
S. Sugita,
I. Takahashi,
J. -L. Atteia,
A. Pelangeon,
R. Vanderspek,
C. Graziani,
G. Prigozhin,
J. Villasenor,
J. G. Jernigan,
G. B. Crew,
T. Sakamoto,
G. R. Ricker
, et al. (14 additional authors not shown)
Abstract:
Using a pulse-fit method, we investigate the spectral lags between the traditional gamma-ray band (50-400 keV) and the X-ray band (6-25 keV) for 8 GRBs with known redshifts (GRB 010921, GRB 020124, GRB 020127, GRB 021211, GRB 030528, GRB 040924, GRB 041006, GRB 050408) detected with the WXM and FREGATE instruments aboard the HETE-2 satellite. We find several relations for the individual GRB puls…
▽ More
Using a pulse-fit method, we investigate the spectral lags between the traditional gamma-ray band (50-400 keV) and the X-ray band (6-25 keV) for 8 GRBs with known redshifts (GRB 010921, GRB 020124, GRB 020127, GRB 021211, GRB 030528, GRB 040924, GRB 041006, GRB 050408) detected with the WXM and FREGATE instruments aboard the HETE-2 satellite. We find several relations for the individual GRB pulses between the spectral lag and other observables, such as the luminosity, pulse duration, and peak energy (Epeak). The obtained results are consistent with those for BATSE, indicating that the BATSE correlations are still valid at lower energies (6-25 keV). Furthermore, we find that the photon energy dependence for the spectral lags can reconcile the simple curvature effect model. We discuss the implication of these results from various points of view.
△ Less
Submitted 26 February, 2010; v1 submitted 19 February, 2010;
originally announced February 2010.
-
Initiation of the detonation in the gravitationally confined detonation model of Type Ia supernovae
Authors:
Ivo R. Seitenzahl,
Casey A. Meakin,
Don Q. Lamb,
James W. Truran
Abstract:
We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have be…
▽ More
We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have been derived from 1D simulations that resolve the initiation of a detonation. By increasing the maximum resolution in a truncated cone that encompasses this region, beginning somewhat before initiation of the detonation occurs, we successfully simulate in situ the first gradient-initiated detonation in a whole-star simulation. The detonation emerges when a compression wave overruns a pocket of fuel situated in a Kelvin-Helmholtz cusp at the leading edge of the inwardly directed jet of burning carbon. The compression wave pre-conditions the temperature in the fuel in such a way that the Zel'dovich gradient mechanism can operate and a detonation ensues. We explore the dependence of the length scale of the temperature gradient on spatial resolution and discuss the implications for the robustness of this detonation mechanism. We find that the time and the location at which initiation of the detonation occurs varies with resolution. In particular, initiation of a detonation had not yet occurred in our highest resolution simulation by the time we ended the simulation because of the computational demand it required. We suggest that the turbulent shear layer surrounding the inwardly directed jet provides the most favorable physical conditions, and therefore the most likely location, for initiation of a detonation in the GCD model.
△ Less
Submitted 19 May, 2009;
originally announced May 2009.
-
Eulerian and Lagrangian Statistics from high resolution Numerical Simulations of weakly compressible turbulence
Authors:
R. Benzi,
L. Biferale,
R. Fisher,
D. Q. Lamb,
F. Toschi
Abstract:
We report a detailed study of Eulerian and Lagrangian statistics from high resolution Direct Numerical Simulations of isotropic weakly compressible turbulence. Reynolds number at the Taylor microscale is estimated to be around 600. Eulerian and Lagrangian statistics is evaluated over a huge data set, made by $1856^3$ spatial collocation points and by 16 million particles, followed for about one…
▽ More
We report a detailed study of Eulerian and Lagrangian statistics from high resolution Direct Numerical Simulations of isotropic weakly compressible turbulence. Reynolds number at the Taylor microscale is estimated to be around 600. Eulerian and Lagrangian statistics is evaluated over a huge data set, made by $1856^3$ spatial collocation points and by 16 million particles, followed for about one large-scale eddy turn over time. We present data for Eulerian and Lagrangian Structure functions up to ten order. We analyse the local scaling properties in the inertial range and in the viscous range. Eulerian results show a good superposition with previous data. Lagrangian statistics is different from existing experimental and numerical results, for moments of sixth order and higher. We interpret this in terms of a possible contamination from viscous scale affecting the estimate of the scaling properties in previous studies. We show that a simple bridge relation based on Multifractal theory is able to connect scaling properties of both Eulerian and Lagrangian observables, provided that the small differences between intermittency of transverse and longitudinal Eulerian structure functions are properly considered.
△ Less
Submitted 1 May, 2009;
originally announced May 2009.
-
Spontaneous Initiation of Detonations in White Dwarf Environments: Determination of Critical Sizes
Authors:
Ivo R. Seitenzahl,
Casey a. Meakin,
Dean M. Townsley,
Don Q. Lamb,
James W. Truran
Abstract:
Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations,…
▽ More
Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations, prompting the use of results of 1D simulations at higher resolution, such as the ones performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized non-uniformity of temperature -- a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in white dwarf matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia.
△ Less
Submitted 23 January, 2009;
originally announced January 2009.
-
Epeak estimator for Gamma-Ray Bursts Observed by the Swift Burst Alert Telescope
Authors:
T. Sakamoto,
G. Sato,
L. Barbier,
S. D. Barthelmy,
J. R. Cummings,
E. E. Fenimore,
N. Gehrels,
D. Hullinger,
H. A. Krimm,
D. Q. Lamb,
C. B. Markwardt,
D. M. Palmer,
A. M. Parsons,
M. Stamatikos,
J. Tueller,
T. N. Ukwatta
Abstract:
We report a correlation based on a spectral simulation study of the prompt emission spectra of gamma-ray bursts (GRBs) detected by the Swift Burst Alert Telescope (BAT). The correlation is between the Epeak energy, which is the peak energy in the νF_νspectrum, and the photon index (Γ) derived from a simple power-law model. The Epeak - Γrelation, assuming the typical smoothly broken power-law spe…
▽ More
We report a correlation based on a spectral simulation study of the prompt emission spectra of gamma-ray bursts (GRBs) detected by the Swift Burst Alert Telescope (BAT). The correlation is between the Epeak energy, which is the peak energy in the νF_νspectrum, and the photon index (Γ) derived from a simple power-law model. The Epeak - Γrelation, assuming the typical smoothly broken power-law spectrum of GRBs, is \log Epeak = 3.258 - 0.829Γ(1.3 < Γ< 2.3). We take into account not only a range of Epeak energies and fluences, but also distributions for both the low-energy photon index and the high-energy photon index in the smoothly broken power-law model. The distribution of burst durations in the BAT GRB sample is also included in the simulation. Our correlation is consistent with the index observed by BAT and Epeak measured by the BAT, and by other GRB instruments. Since about 85% of GRBs observed by the BAT are acceptably fit with the simple power-law model because of the relatively narrow energy range of the BAT, this relationship can be used to estimate Epeak when it is located within the BAT energy range.
△ Less
Submitted 20 November, 2008;
originally announced November 2008.
-
Intrinsic properties of a complete sample of HETE-2 gamma-ray bursts. A measure of the GRB rate in the Local Universe
Authors:
A. Pélangeon,
J. -L. Atteia,
Y. E. Nakagawa,
K. Hurley,
A. Yoshida,
R. Vanderspek,
M. Suzuki,
N. Kawai,
G. Pizzichini,
M. Boër,
J. Braga,
G. Crew,
T. Q. Donaghy,
J. P. Dezalay,
J. Doty,
E. E. Fenimore,
M. Galassi,
C. Graziani,
J. G. Jernigan,
D. Q. Lamb,
A. Levine,
J. Manchanda,
F. Martel,
M. Matsuoka,
J. -F. Olive
, et al. (10 additional authors not shown)
Abstract:
Aims: Taking advantage of the forthcoming Catalog of the HETE-2 mission, the aim of this paper is to evaluate the main properties of HETE-2 GRBs - the E_peak, the T_90 and the E_iso - in their source frames and to derive their unbiased distribution. Methods: We first construct a complete sample containing all the bursts localized by the WXM on-board HETE-2, which are selected with a uniform crit…
▽ More
Aims: Taking advantage of the forthcoming Catalog of the HETE-2 mission, the aim of this paper is to evaluate the main properties of HETE-2 GRBs - the E_peak, the T_90 and the E_iso - in their source frames and to derive their unbiased distribution. Methods: We first construct a complete sample containing all the bursts localized by the WXM on-board HETE-2, which are selected with a uniform criterion and whose observed parameters can be constrained. We then derive the intrinsic E_peak, T_90 and E_iso distributions using their redshift when it is available, or their pseudo-redshift otherwise. We finally compute the number of GRB (N_Vmax) within the visibility volume (V_max) of each GRB, in order to derive a weight for each detected burst accounting both for the detection significance and the star formation history of the universe. Results: The unbiased distributions obtained clearly show the predominence of X-ray flashes (XRFs) in the global GRB population. We also derive the rate of local GRBs: R0^H2 > 11 Gpc-3 yr-1, which is intermediate between the local rate obtained by considering only the high-luminosity bursts (~1 Gpc-3 yr-1) and that obtained by including the low-luminosity bursts (>200 Gpc-3 yr-1).
△ Less
Submitted 20 November, 2008;
originally announced November 2008.
-
Rings of Dark Matter in Collisions Between Clusters of Galaxies
Authors:
J. A. ZuHone,
D. Q. Lamb,
P. M. Ricker
Abstract:
Several lines of evidence suggest that the galaxy cluster Cl0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. Recent lensing observations suggest the presence of a ring-like dark matter structure, which has been interpreted as the result of such a collision. In this paper we present $N$-body simulations of cluster c…
▽ More
Several lines of evidence suggest that the galaxy cluster Cl0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. Recent lensing observations suggest the presence of a ring-like dark matter structure, which has been interpreted as the result of such a collision. In this paper we present $N$-body simulations of cluster collisions along the line of sight to investigate the detectability of such features. We use realistic dark matter density profiles as determined from cosmological simulations. Our simulations show a "shoulder" in the dark matter distribution after the collision, but no ring feature even when the initial particle velocity distribution is highly tangentially anisotropic ($σ_θ/σ_r >> 1$). Only when the initial particle velocity distribution is circular do our simulations show such a feature. Even modestly anisotropic velocity distributions are inconsistent with the halo velocity distributions seen in cosmological simulations, and would require highly fine-tuned initial conditions. Our investigation leaves us without an explanation for the dark matter ring-like feature in Cl 0024+17 suggested by lensing observations.
△ Less
Submitted 9 July, 2009; v1 submitted 18 September, 2008;
originally announced September 2008.
-
A Line-Of-Sight Galaxy Cluster Collision: Simulated X-Ray Observations
Authors:
J. A. Zuhone,
P. M. Ricker,
D. Q. Lamb,
H. -Y. Yang
Abstract:
Several lines of evidence have suggested that the the galaxy cluster Cl~0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. In this paper we present a high-resolution $N$-body/hydrodynamics simulation of such a collision. We have created mock X-ray observations of our simulated system using MARX, a program that simula…
▽ More
Several lines of evidence have suggested that the the galaxy cluster Cl~0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. In this paper we present a high-resolution $N$-body/hydrodynamics simulation of such a collision. We have created mock X-ray observations of our simulated system using MARX, a program that simulates the on-orbit performance of the Chandra X-ray Observatory. We analyze these simulated data to generate radial profiles of the surface brightness and temperature. At later times, $t = 2.0-3.0$ Gyr after the collision, the simulated surface brightness profiles are better fit by a superposition of two $β$-model profiles than a single profile, in agreement with the observations of Cl~0024+17. In general, due to projection effects, much of the post-collision density and temperature structure of the clusters is not seen in the observations. In particular, the observed temperatures from spectral fitting are much lower than the temperature of the hottest gas. We determine from our fitted profiles that if the system is modeled as a single cluster, the hydrostatic mass estimate is a factor $\sim$2-3 less than the actual mass, but if the system is modeled as two galaxy clusters in superposition, a hydrostatic mass estimation can be made which is accurate to within $\sim$10%. We examine some implications of these results for galaxy cluster X-ray surveys.
△ Less
Submitted 9 July, 2009; v1 submitted 6 August, 2008;
originally announced August 2008.
-
Multiple Component Analysis of Time Resolved Spectra of GRB041006: A Clue to the Nature of Underlying Soft Component of GRBs
Authors:
Yuji Shirasaki,
Atsumasa Yoshida,
Nobuyuki Kawai,
Toru Tamagawa,
Takanori Sakamoto,
Motoko Suzuki,
Yujin Nakagawam,
Akina Kobayashi,
Satoshi Sugita,
Ichiro Takahashi,
Makoto Arimoto,
Takashi Shimokawabe,
Nicolas Vasquez Pazmino,
Takuto Ishimura,
Rie Sato,
Masaru Matsuoka,
Edward E. Fenimore,
Mark Galassi,
Donald Q. Lamb,
Carlo Graziani,
Timothy Q. Donaghy,
Jean-Luc Atteia,
Alexandre Pelangeon,
Roland Vanderspek,
Geoffrey B. Crew
, et al. (8 additional authors not shown)
Abstract:
GRB 041006 was detected by HETE-2 at 12:18:08 UT on 06 October 2004. This GRB displays a soft X-ray emission, a precursor before the onset of the main event, and also a soft X-ray tail after the end of the main peak. The light curves in four different energy bands display different features; At higher energy bands several peaks are seen in the light curve, while at lower energy bands a single br…
▽ More
GRB 041006 was detected by HETE-2 at 12:18:08 UT on 06 October 2004. This GRB displays a soft X-ray emission, a precursor before the onset of the main event, and also a soft X-ray tail after the end of the main peak. The light curves in four different energy bands display different features; At higher energy bands several peaks are seen in the light curve, while at lower energy bands a single broader bump dominates. It is expected that these different features are the result of a mixture of several components each of which has different energetics and variability. To reveal the nature of each component, we analysed the time resolved spectra and they are successfully resolved into several components. We also found that these components can be classified into two distinct classes; One is a component which has an exponential decay of $E_{p}$ with a characteristic timescale shorter than $\sim$ 30 sec, and its spectrum is well represented by a broken power law function, which is frequently observed in many prompt GRB emissions, so it should have an internal-shock origin. Another is a component whose $E_{p}$ is almost unchanged with characteristic timescale longer than $\sim$ 60 sec, and shows a very soft emission and slower variability. The spectrum of the soft component is characterized by either a broken power law or a black body spectrum. This component might originate from a relatively wider and lower velocity jet or a photosphere of the fireball. By assuming that the soft component is a thermal emission, the radiation radius is initially $4.4 \times 10^{6}$ km, which is a typical radius of a blue supergiant, and its expansion velocity is $2.4 \times 10^{5}$ km/s in the source frame.
△ Less
Submitted 15 April, 2008; v1 submitted 12 February, 2008;
originally announced February 2008.
-
Global Properties of X-Ray Flashes and X-Ray-Rich Gamma-Ray Bursts Observed by Swift
Authors:
T. Sakamoto,
D. Hullinger,
G. Sato,
R. Yamazaki,
L. Barbier,
S. D. Barthelmy,
J. R. Cummings,
E. E. Fenimore,
N. Gehrels,
H. A. Krimm,
D. Q. Lamb,
C. B. Markwardt,
J. P. Osborne,
D. M. Palmer,
A. M. Parsons,
M. Stamatikos,
J. Tueller
Abstract:
We describe and discuss the spectral and temporal characteristics of the prompt emission and X-ray afterglow emission of X-ray flashes (XRFs) and X-ray-rich gamma-ray bursts (XRRs) detected and observed by Swift between December 2004 and September 2006. We compare these characteristics to a sample of conventional classical gamma-ray bursts (C-GRBs) observed during the same period. We confirm the…
▽ More
We describe and discuss the spectral and temporal characteristics of the prompt emission and X-ray afterglow emission of X-ray flashes (XRFs) and X-ray-rich gamma-ray bursts (XRRs) detected and observed by Swift between December 2004 and September 2006. We compare these characteristics to a sample of conventional classical gamma-ray bursts (C-GRBs) observed during the same period. We confirm the correlation between Epeak_obs and fluence noted by others and find further evidence that XRFs, XRRs and C-GRBs form a continuum. We also confirm that our known redshift sample is consistent with the correlation between the peak energy in the GRB rest frame (Epeak_src) and the isotropic radiated energy (Eiso), so called the Epeak_src-Eiso relation. The spectral properties of X-ray afterglows of XRFs and C-GRBs are similar, but the temporal properties of XRFs and C-GRBs are quite different. We found that the light curves of C-GRB afterglows show a break to steeper indices (shallow-to-steep break) at much earlier times than do XRF afterglows. Moreover, the overall luminosity of XRF X-ray afterglows is systematically smaller by a factor of two or more compared to that of C-GRBs. These distinct differences between the X-ray afterglows of XRFs and C-GRBs may be the key to understanding not only the mysterious shallow-to-steep break in X-ray afterglow light curves, but also the unique nature of XRFs.
△ Less
Submitted 28 January, 2008;
originally announced January 2008.
-
Flame Evolution During Type Ia Supernovae and the Deflagration Phase in the Gravitationally Confined Detonation Scenario
Authors:
D. M. Townsley,
A. C. Calder,
S. M. Asida,
I. R. Seitenzahl,
F. Peng,
N. Vladimirova,
D. Q. Lamb,
J. W. Truran
Abstract:
We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order t…
▽ More
We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.
△ Less
Submitted 7 June, 2007;
originally announced June 2007.
-
Three-Dimensional Simulations of the Deflagration Phase of the Gravitationally Confined Detonation Model of Type Ia Supernovae
Authors:
G C Jordan IV,
R T Fisher,
D M Townsley,
A C Calder,
C Graziani,
S Asida,
D Q Lamb,
J W Truran
Abstract:
We report the results of a series of three-dimensional (3-D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite breakout on the…
▽ More
We report the results of a series of three-dimensional (3-D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3-D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly-directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly-directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of $^{56}$Ni and small amounts of intermediate mass elements. These particular simulations are therefore relevant only to high luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.
△ Less
Submitted 16 March, 2008; v1 submitted 21 March, 2007;
originally announced March 2007.
-
HETE-2 Observations of the X-Ray Flash XRF 040916
Authors:
Makoto Arimoto,
Nobuyuki Kawai,
Atsumasa Yoshida,
Toru Tamagawa,
Yuji Shirasaki,
Motoko Suzuki,
Masaru Matsuoka,
Jun'ichi Kotoku,
Rie Sato,
Takashi Shimokawabe,
Nicolas Vasquez Pazmino,
Takuto Ishimura,
Yujin Nakagawa,
Nobuyuki Ishikawa,
Akina Kobayashi,
Satoshi Sugita,
Ichiro Takahashi,
Makoto Kuwahara,
Makoto Yamauchi,
Kunio Takagishi,
Isamu Hatsukade,
Jean-Luc Atteia,
Alexandre Pelangeon,
Roland Vanderspek,
Carlo Graziani
, et al. (21 additional authors not shown)
Abstract:
A long X-ray flash was detected and localized by the instruments aboard the High Energy Transient Explorer II (HETE-2) at 00:03:30 UT on 2004 September 16. The position was reported to the GRB Coordinates Network (GCN) approximately 2 hours after the burst. This burst consists of two peaks separated by 200 s, with durations of 110 s and 60 s. We have analyzed the energy spectra of the 1st and 2n…
▽ More
A long X-ray flash was detected and localized by the instruments aboard the High Energy Transient Explorer II (HETE-2) at 00:03:30 UT on 2004 September 16. The position was reported to the GRB Coordinates Network (GCN) approximately 2 hours after the burst. This burst consists of two peaks separated by 200 s, with durations of 110 s and 60 s. We have analyzed the energy spectra of the 1st and 2nd peaks observed with the Wide Field X-Ray Monitor (WXM) and the French Gamma Telescope (FREGATE). We discuss the origin of the 2nd peak in terms of flux variabilities and timescales. We find that it is most likely part of the prompt emission, and is explained by the long-acting engine model. This feature is similar to some bright X-ray flares detected in the early afterglow phase of bursts observed by the Swift satellite.
△ Less
Submitted 7 March, 2007; v1 submitted 2 March, 2007;
originally announced March 2007.
-
Discovery of a New X-Ray Burst/Millisecond Accreting Pulsar HETE J1900.1-2455
Authors:
M. Suzuki,
N. Kawai,
T. Tamagawa,
A. Yoshida,
Y. E. Nakagawa,
K. Tanaka,
Y. Shirasaki,
M. Matsuoka,
G. R. Ricker,
R. Vanderspek,
N. Butler,
D. Q. Lamb,
C. Graziani,
G. Pizzichini,
R. Sato,
M. Arimoto,
J. Kotoku,
M. Maetou,
M. Yamauchi
Abstract:
A class of low-mass X-ray binary sources are known to be both X-ray burst sources and millisecond pulsars at the same time. A new source of this class was discovered by High Energy Transient Explorer 2 (HETE-2) on 14 June 2005 as a source of type-I X-ray bursts, which was named HETE J1900.1-2455. Five X-ray bursts from HETE J1900.1-2455 were observed during the summer of 2005. The time resolved…
▽ More
A class of low-mass X-ray binary sources are known to be both X-ray burst sources and millisecond pulsars at the same time. A new source of this class was discovered by High Energy Transient Explorer 2 (HETE-2) on 14 June 2005 as a source of type-I X-ray bursts, which was named HETE J1900.1-2455. Five X-ray bursts from HETE J1900.1-2455 were observed during the summer of 2005. The time resolved spectral analysis of these bursts have revealed that their spectra are consistent with the blackbody radiation throughout the bursts. The bursts show the indication of radius expansion. The bolometric flux remains almost constant during the photospheric radius expansion while blackbody temperature dropped during the same period. Assuming that the flux reached to the Eddington limit on a standard 1.4 solar mass neutron star with a helium atmosphere, we estimate the distance to the source to be $\sim$ 4 kpc.
△ Less
Submitted 11 December, 2006;
originally announced December 2006.
-
Swift Discovery of Gamma-Ray Bursts without Jet Break Feature in their X-Ray Afterglows
Authors:
G. Sato,
R. Yamazaki,
K. Ioka,
T. Sakamoto,
T. Takahashi,
K. Nakazawa,
T. Nakamura,
K. Toma,
D. Hullinger,
M. Tashiro,
A. M. Parsons,
H. A. Krimm,
S. D. Barthelmy,
N. Gehrels,
D. N. Burrows,
P. T. O'Brien,
J. P. Osborne,
G. Chincarini,
D. Q. Lamb
Abstract:
We analyze Swift gamma-ray bursts (GRBs) and X-ray afterglows for three GRBs with spectroscopic redshift determinations -- GRB 050401, XRF 050416a, and GRB 050525a. We find that the relation between spectral peak energy and isotropic energy of prompt emissions (the Amati relation) is consistent with that for the bursts observed in pre-Swift era. However, we find that the X-ray afterglow lightcur…
▽ More
We analyze Swift gamma-ray bursts (GRBs) and X-ray afterglows for three GRBs with spectroscopic redshift determinations -- GRB 050401, XRF 050416a, and GRB 050525a. We find that the relation between spectral peak energy and isotropic energy of prompt emissions (the Amati relation) is consistent with that for the bursts observed in pre-Swift era. However, we find that the X-ray afterglow lightcurves, which extend up to 10-70 days, show no sign of the jet break that is expected in the standard framework of collimated outflows. We do so by showing that none of the X-ray afterglow lightcurves in our sample satisfies the relation between the spectral and temporal indices that is predicted for the phase after jet break. The jet break time can be predicted by inverting the tight empirical relation between the peak energy of the spectrum and the collimation-corrected energy of the prompt emission (the Ghirlanda relation). We find that there are no temporal breaks within the predicted time intervals in X-ray band. This requires either that the Ghirlanda relation has a larger scatter than previously thought, that the temporal break in X-rays is masked by some additional source of X-ray emission, or that it does not happen because of some unknown reason.
△ Less
Submitted 9 November, 2006; v1 submitted 4 November, 2006;
originally announced November 2006.
-
Capturing the Fire: Flame Energetics and Neutronizaton for Type Ia Supernova Simulations
Authors:
A. C. Calder,
D. M. Townsley,
I. R. Seitenzahl,
F. Peng,
O. E. B. Messer,
N. Vladimirova,
E. F. Brown,
J. W. Truran,
D. Q. Lamb
Abstract:
We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release…
▽ More
We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.
△ Less
Submitted 1 November, 2006;
originally announced November 2006.
-
Two Additions to the New Class of Low Accretion-Rate Magnetic Binaries
Authors:
Gary D. Schmidt,
Paula Szkody,
Arne Henden,
Scott F. Anderson,
Don Q. Lamb,
Bruce Margon,
Donald P. Schneider
Abstract:
Two new magnetic white dwarf accretion binaries with extremely low mass-transfer rates have been discovered in the course of the Sloan Digital Sky Survey. Measured magnetic fields are 42 MG and 57 MG, and one system orbits with a period of just 82 min. The new systems therefore significantly expand the range in properties exhibited by the small class. The measured accretion rates are very low, 0…
▽ More
Two new magnetic white dwarf accretion binaries with extremely low mass-transfer rates have been discovered in the course of the Sloan Digital Sky Survey. Measured magnetic fields are 42 MG and 57 MG, and one system orbits with a period of just 82 min. The new systems therefore significantly expand the range in properties exhibited by the small class. The measured accretion rates are very low, 0.6-5x10^{-13} Msun/yr, and multiple visits spanning more than a year confirm that this is not a short-lived characteristic. It is becoming increasingly clear that the low-mdot magnetic white dwarf binaries accrete by nearly complete magnetic capture of the stellar wind from the secondary star rather than by Roche lobe overflow. The accretion rates therefore provide some of the first realistic estimates of the total wind loss rates from M dwarfs. Although one or more of the eight systems known to date may be interrupted or possibly even extinct Polars, several lines of evidence suggest that most are pre-Polars whose evolution has not yet brought the secondaries into contact with their Roche surfaces. Considering the difficulties of identifying binaries over a wide range in field strength and accretion rate, it is quite possible that the space density of wind-accreting magnetic binaries exceeds that of the classical X-ray emitting, Roche-lobe overflow Polars.
△ Less
Submitted 27 October, 2006;
originally announced October 2006.
-
X-ray Flashes or soft Gamma-ray Bursts? The case of the likely distant XRF 040912
Authors:
G. Stratta,
S. Basa,
N. Butler,
J. L. Atteia,
B. Gendre,
A. Pelangeon,
F. Malacrino,
Y. Mellier,
D. A. Kann,
S. Klose,
A. Zeh,
N. Masetti,
E. Palazzi,
J. Gorosabel,
A. J. Castro-Tirado,
A. de Ugarte Postigo,
M. Jelinek,
J. Cepa,
H. Castaneda,
D. Martinez-Delgado,
M. Boer,
J. Braga,
G. Crew,
T. Q. Donaghy,
J. -P. Dezalay
, et al. (27 additional authors not shown)
Abstract:
In this work, we present a multi-wavelength study of XRF 040912, aimed at measuring its distance scale and the intrinsic burst properties. We performed a detailed spectral and temporal analysis of both the prompt and the afterglow emission and we estimated the distance scale of the likely host galaxy. We then used the currently available sample of XRFs with known distance to discuss the connecti…
▽ More
In this work, we present a multi-wavelength study of XRF 040912, aimed at measuring its distance scale and the intrinsic burst properties. We performed a detailed spectral and temporal analysis of both the prompt and the afterglow emission and we estimated the distance scale of the likely host galaxy. We then used the currently available sample of XRFs with known distance to discuss the connection between XRFs and classical Gamma-ray Bursts (GRBs). We found that the prompt emission properties unambiguously identify this burst as an XRF, with an observed peak energy of E_p=17+/-13 keV and a burst fluence ratio S(2-30keV)/S(30-400keV)>1. A non-fading optical source with R~24 mag and with an apparently extended morphology is spatially consistent with the X-ray afterglow, likely the host galaxy. XRF 040912 is a very dark burst since no afterglow optical counterpart is detected down to R>25 mag (3 sigma limiting magnitude) at 13.6 hours after the burst. The host galaxy spectrum detected from 3800A to 10000A, shows a single emission line at 9552A. The lack of any other strong emission lines blue-ward of the detected one and the absence of the Ly alpha cut-off down to 3800A are consistent with the hypothesis of the [OII] line at redshift z=1.563+/-0.001. The intrinsic spectral properties rank this XRF among the soft GRBs in the E_peak-E_iso diagram. Similar results were obtained for most XRFs at known redshift. Only XRF 060218 and XRF 020903 represent a good example of instrinsic XRF(i-XRF) and are possibly associated with a different progenitor population. This scenario may calls for a new definition of XRFs.
△ Less
Submitted 5 October, 2006;
originally announced October 2006.
-
HETE-2 Localizations and Observations of Four Short Gamma-Ray Bursts: GRBs 010326B, 040802, 051211 and 060121
Authors:
T. Q. Donaghy,
D. Q. Lamb,
T. Sakamoto,
J. P. Norris,
Y. Nakagawa,
J. Villasenor,
J. -L. Atteia,
R. Vanderspek,
C. Graziani,
N. Kawai,
G. R. Ricker,
G. B. Crew,
J. Doty,
G. Prigozhin,
J. G. Jernigan,
Y. Shirasaki,
M. Suzuki,
N. Butler,
K. Hurley,
T. Tamagawa,
A. Yoshida,
M. Matsuoka,
E. E. Fenimore,
M. Galassi,
M. Boer
, et al. (12 additional authors not shown)
Abstract:
Here we report the localizations and properties of four short-duration GRBs localized by the High Energy Transient Explorer 2 satellite (HETE-2): GRBs 010326B, 040802, 051211 and 060121, all of which were detected by the French Gamma Telescope (Fregate) and localized with the Wide-field X-ray Monitor (WXM) and/or Soft X-ray Camera (SXC) instruments. We discuss eight possible criteria for determi…
▽ More
Here we report the localizations and properties of four short-duration GRBs localized by the High Energy Transient Explorer 2 satellite (HETE-2): GRBs 010326B, 040802, 051211 and 060121, all of which were detected by the French Gamma Telescope (Fregate) and localized with the Wide-field X-ray Monitor (WXM) and/or Soft X-ray Camera (SXC) instruments. We discuss eight possible criteria for determining whether these GRBs are "short population bursts" (SPBs) or "long population bursts" (LPBs). These criteria are (1) duration, (2) pulse widths, (3) spectral hardness, (4) spectral lag, (5) energy Egamma radiated in gamma rays (or equivalently, the kinetic energy E_KE of the GRB jet), (6) existence of a long, soft bump following the burst, (7) location of the burst in the host galaxy, and (8) type of host galaxy. In particular, we have developed a likelihood method for determining the probability that a burst is an SPB or a LPB on the basis of its T90 duration alone. A striking feature of the resulting probability distribution is that the T90 duration at which a burst has an equal probability of being a SPB or a LPB is T90 = 5 s, not T90 = 2 s, as is often used. All four short-duration bursts discussed in detail in this paper have T90 durations in the Fregate 30-400 keV energy band of 1.90, 2.31, 4.25, and 1.97 sec, respectively, yielding probabilities P(S|T90) = 0.97, 0.91, 0.60, and 0.95 that these bursts are SPBs on the basis of their T90 durations alone. All four bursts also have spectral lags consistent with zero. These results provide strong evidence that all four GRBs are SPBs (abstract continues).
△ Less
Submitted 3 June, 2006; v1 submitted 22 May, 2006;
originally announced May 2006.
-
GRB 060121: Implications of a Short/Intermediate Duration Gamma-Ray Burst at High Redshift
Authors:
A. de Ugarte Postigo,
A. J. Castro-Tirado,
S. Guziy,
J. Gorosabel,
G. Johannesson,
M. A. Aloy,
S. McBreen,
D. Q. Lamb,
N. Benitez,
M. Jelinek,
S. B. Pandey,
D. Coe,
M. D. Perez-Ramirez,
F. J. Aceituno,
M. Alises,
J. A. Acosta-Pulido,
G. Gomez,
R. Lopez,
T. Q. Donaghy,
Y. E. Nakagawa,
T. Sakamoto,
G. R. Ricker,
F. R. Hearty,
M. Bayliss,
G. Gyuk
, et al. (1 additional authors not shown)
Abstract:
Since the discovery of the first short-hard gamma-ray burst afterglows in 2005, the handful of observed events have been found to be embedded in nearby (z < 1), bright underlying galaxies. We present multiwavelength observations of the short-duration burst GRB 060121, which is the first observed to clearly outshine its host galaxy (by a factor >10^2). A photometric redshift for this event places…
▽ More
Since the discovery of the first short-hard gamma-ray burst afterglows in 2005, the handful of observed events have been found to be embedded in nearby (z < 1), bright underlying galaxies. We present multiwavelength observations of the short-duration burst GRB 060121, which is the first observed to clearly outshine its host galaxy (by a factor >10^2). A photometric redshift for this event places the progenitor at a most probable redshift of z = 4.6, with a less probable scenario of z = 1.7. In either case, GRB 060121 could be the farthermost short-duration GRB detected to date and implies an isotropic-equivalent energy release in gamma-rays comparable to that seen in long-duration bursts. We discuss the implications of the released energy on the nature of the progenitor. These results suggest that GRB 060121 may belong to a family of energetic short-duration events, lying at z > 1 and whose optical afterglows would outshine their host galaxies, unlike the first short-duration GRBs observed in 2005. The possibility of GRB 060121 being an intermediate duration burst is also discussed.
△ Less
Submitted 19 May, 2006;
originally announced May 2006.
-
An Optically Dark GRB Observed by HETE-2: GRB 051022
Authors:
Yujin E. Nakagawa,
Atsumasa Yoshida,
Satoshi Sugita,
Kaoru Tanaka,
Nobuyuki Ishikawa,
Toru Tamagawa,
Motoko Suzuki,
Yuji Shirasaki,
Nobuyuki Kawai,
Masaru Matsuoka,
Jean-Luc Atteia,
Alexandre Pelangeon,
Roland Vanderspek,
Geoff B. Crew,
Joel S. Villasenor,
Nat Butler,
John Doty,
George R. Ricker,
Graziella Pizzichini,
Timothy Q. Donaghy,
Donald Q. Lamb,
Carlo Graziani,
Rie Sato,
Miki Maetou,
Makoto Arimoto
, et al. (11 additional authors not shown)
Abstract:
GRB 051022 was detected at 13:07:58 on 22 October 2005 by HETE-2. The location of GRB 051022 was determined immediately by the flight localization system. This burst contains multiple pulses and has a rather long duration of about 190 seconds. The detections of candidate X-ray and radio afterglows were reported, whereas no optical afterglow was found. The optical spectroscopic observations of th…
▽ More
GRB 051022 was detected at 13:07:58 on 22 October 2005 by HETE-2. The location of GRB 051022 was determined immediately by the flight localization system. This burst contains multiple pulses and has a rather long duration of about 190 seconds. The detections of candidate X-ray and radio afterglows were reported, whereas no optical afterglow was found. The optical spectroscopic observations of the host galaxy revealed the redshift z = 0.8. Using the data derived by HETE-2 observation of the prompt emission, we found the absorption N_H = 8.8 -2.9/+3.1 x 10^22 cm^-2 and the visual extinction A_V = 49 -16/+17 mag in the host galaxy. If this is the case, no detection of any optical transient would be quite reasonable. The absorption derived by the Swift XRT observations of the afterglow is fully consistent with those obtained from the early HETE-2 observation of the prompt emission. Our analysis implies an interpretation that the absorbing medium could be outside external shock at R ~ 10^16 cm, which may be a dusty molecular cloud.
△ Less
Submitted 6 July, 2006; v1 submitted 17 May, 2006;
originally announced May 2006.
-
A Study of Compact Object Mergers as Short Gamma-ray Burst Progenitors
Authors:
Krzysztof Belczynski,
Rosalba Perna,
Tomasz Bulik,
Vassiliki Kalogera,
Natalia Ivanova,
Donald Q. Lamb
Abstract:
We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for di…
▽ More
We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B)
△ Less
Submitted 26 May, 2006; v1 submitted 19 January, 2006;
originally announced January 2006.
-
An improved redshift indicator for Gamma-Ray Bursts, based on the prompt emission
Authors:
Alexandre Pélangeon,
Jean-Luc Atteia,
Don Q. Lamb,
George R. Ricker,
the the HETE Science Team Collaboration
Abstract:
We propose an improved version of the redshift indicator developed by Atteia (2003), which gets rid of the dependence on the burst duration and provides better estimates for high-redshift GRBs. We present the derivation and the definition of this redshift indicator, then its calibration with 17 GRBs with known redshifts detected by HETE-2 and 2 more detected by Konus-Wind. We also provide an est…
▽ More
We propose an improved version of the redshift indicator developed by Atteia (2003), which gets rid of the dependence on the burst duration and provides better estimates for high-redshift GRBs. We present the derivation and the definition of this redshift indicator, then its calibration with 17 GRBs with known redshifts detected by HETE-2 and 2 more detected by Konus-Wind. We also provide an estimation of the redshift for 59 bursts, and we finally discuss the redshift distribution of HETE-bursts and the possible other applications of this redshift indicator.
△ Less
Submitted 7 January, 2006;
originally announced January 2006.
-
Discovery of the short gamma-ray burst GRB 050709
Authors:
J. S. Villasenor,
D. Q. Lamb,
G. R. Ricker,
J. -L. Atteia,
N. Kawai,
N. Butler,
Y. Nakagawa,
J. G. Jernigan,
M. Boer,
G. B. Crew,
T. Q. Donaghy,
J. Doty,
E. E. Fenimore,
M. Galassi,
C. Graziani,
K. Hurley,
A. Levine,
F. Martel,
M. Matsuoka,
J. -F. Olive,
G. Prigozhin,
T. Sakamoto,
Y. Shirasaki,
M. Suzuki,
T. Tamagawa
, et al. (8 additional authors not shown)
Abstract:
Gamma-Ray Bursts (GRBs) fall into two classes: short-hard and long-soft bursts. The latter are now known to have X-ray and optical afterglows, to occur at cosmological distances in star-forming galaxies, and to be associated with the explosion of massive stars. In contrast, the distance scale, the energy scale, and the progenitors of short bursts have remained a mystery. Here we report the disco…
▽ More
Gamma-Ray Bursts (GRBs) fall into two classes: short-hard and long-soft bursts. The latter are now known to have X-ray and optical afterglows, to occur at cosmological distances in star-forming galaxies, and to be associated with the explosion of massive stars. In contrast, the distance scale, the energy scale, and the progenitors of short bursts have remained a mystery. Here we report the discovery of a short-hard burst whose accurate localization has led to follow-up observations that have identified the X-ray afterglow and (for the first time) the optical afterglow of a short-hard burst. These, in turn, have led to identification of the host galaxy of the burst as a late-type galaxy at z=0.16 showing that at least some short-hard bursts occur at cosmological distances in the outskirts of galaxies, and are likely to be due to the merging of compact binaries.
△ Less
Submitted 6 October, 2005;
originally announced October 2005.
-
A Gamma-Ray Burst Mission to Investigate the Properties of Dark Energy
Authors:
D. Q. Lamb,
G. R. Ricker,
D. Lazzati,
G. Ghirlanda,
G. Ghisellini,
C. Firmani,
L. Amati,
J. -L. Atteia,
V. Avila-Reese,
S. Burles,
N. Butler,
H. -W. Chen,
E. Costa,
J. Doty,
F. Frontera,
A. Fruchter,
P. Garnavich,
C. Graziani,
J. G. Jernigan,
N. Kawai,
P. Mazzali,
P. Meszaros,
L. Piro,
T. Sakamoto,
K. Stanek
, et al. (4 additional authors not shown)
Abstract:
[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae exce…
▽ More
[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae excellent standard candles for cosmology, so too, the relations between Epk, Eiso, and Egamma point toward a methodology for using GRBs as excellent standard candles for cosmology. In addition, GRBs occur over the broad redshift range from z=0.1 to at least z=4.5, and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great promise as cosmological ``yardsticks'' to measure the rate of expansion of the universe over time, and therefore the properties of dark energy (i.e., Omega_M, Omega_Lambda, w_0, and w_a).
We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the redshift range 0.1 >~ z <~ 10 during a 2-year mission. This burst sample would enable both Omega_M and w_0 to be determined to +/- 0.07 and +/- 0.06 (68% CL), respectively, and w_a to be significantly constrained.
△ Less
Submitted 14 July, 2005;
originally announced July 2005.
-
Gamma-Ray Burst Jet Profiles And Their Signatures
Authors:
C. Graziani,
D. Q. Lamb,
T. Q. Donaghy
Abstract:
HETE-II and BeppoSAX have produced a sample of GRBs and XRFs with known redshifts and $E_{pk}$. This sample provides four important empirical constraints on the nature of the source jets: Log $E_{iso}$ is approximately uniformly distributed over several orders of magnitude; the inferred prompt energy Log $E_γ$ is narrowly distributed; the Amati relation holds between $E_{iso}$ and $E_{pk}$; and…
▽ More
HETE-II and BeppoSAX have produced a sample of GRBs and XRFs with known redshifts and $E_{pk}$. This sample provides four important empirical constraints on the nature of the source jets: Log $E_{iso}$ is approximately uniformly distributed over several orders of magnitude; the inferred prompt energy Log $E_γ$ is narrowly distributed; the Amati relation holds between $E_{iso}$ and $E_{pk}$; and the Ghirlanda relation holds between $E_γ$ and $E_{pk}$.
We explore the implications of these constraints for GRB jet structure during the prompt emission phase. We infer the underlying angular profiles from the first two of the above constraints assuming all jets have the same profile and total energy, and show that such ``universal jet'' models cannot satisfy both constraints.
We introduce a general and efficient method for calculating relativistic emission distributions and $E_{pk}$ distributions from jets with arbitrary (smooth) angular jet profiles. We also exhibit explicit analytical formulas for emission from top-hat jets (which are not smooth). We use these methods to exhibit $E_{pk}$ and $E_{iso}$ as a function of viewing angle, for several interesting families of GRB jet profiles. We use the same methods to calculate expected frequency distributions of $E_{iso}$ and $E_γ$ for the same families of models.
We then proceed to explore the behavior of universal jet models under a range of profile shapes and parameters, to map the extent to which these models can conform to the above four empirical constraints.
△ Less
Submitted 31 May, 2005;
originally announced May 2005.
-
Early-Time Chromatic Variations in the Wind-Swept Medium of GRB 021211 and the Faintness of its Afterglow
Authors:
M. C. Nysewander,
D. E. Reichart,
H. -S. Park,
G. G. Williams,
K. Kinugasa,
D. Q. Lamb,
A. A. Henden,
S. Klose,
T. Kato,
A. Harper,
H. Yamaoka,
C. Laws,
K. Torii,
D. G. York,
J. C. Barentine,
J. Dembicky,
R. J. McMillan,
J. A. Moran,
D. H. Hartmann,
B. Ketzeback,
M. B. Bayliss,
J. W. Bartelme,
J. A. Crain,
A. C. Foster,
M. Schwartz
, et al. (6 additional authors not shown)
Abstract:
We present Follow-Up Network for Gamma-Ray Bursts (FUN GRB) Collaboration observations of the optical afterglow of GRB 021211 made between 143 seconds and 102 days after the burst. Our unique data set includes the earliest filtered detections and color information for an afterglow in the pre-Swift era. We find that the afterglow is best described by (1) propagation through a wind-swept medium, (…
▽ More
We present Follow-Up Network for Gamma-Ray Bursts (FUN GRB) Collaboration observations of the optical afterglow of GRB 021211 made between 143 seconds and 102 days after the burst. Our unique data set includes the earliest filtered detections and color information for an afterglow in the pre-Swift era. We find that the afterglow is best described by (1) propagation through a wind-swept medium, (2) a cooling break that is blueward of the observed optical frequencies, and (3) a hard electron energy distribution. However, superimposed on this "standard model" behavior we find significant chromatic variations within the first few hours after the burst. We consider possible reasons for these variations, including the possibility that they are due to a dust echo. Finally, we constrain physical parameters that describe the afterglow and surrounding medium for a variety of scenarios and find that GRB 021211's afterglow is faint for a combination of 3-4 reasons: (1) a low fraction of energy in relativistic electrons, (2) a low density for the wind-swept medium, implying either a low mass-loss rate and/or a high wind velocity for the progenitor, (3) a wide opening/viewing angle for the jet, and possibly (4) moderate source frame extinction. The jet appears to be significantly far from equipartition and magnetically dominated. More extreme versions of this might explain the darkness of many afterglows in the Swift era.
△ Less
Submitted 16 June, 2006; v1 submitted 23 May, 2005;
originally announced May 2005.
-
Jet Models of X-Ray Flashes
Authors:
D. Q. Lamb,
T. Q. Donaghy,
C. Graziani
Abstract:
One third of all HETE-2--localized bursts are X-Ray Flashes (XRFs), a class of events first identified by Heise in which the fluence in the 2-30 keV energy band exceeds that in the 30-400 keV energy band. We summarize recent HETE-2 and other results on the properties of XRFs. These results show that the properties of XRFs, X-ray-rich gamma-ray bursts (GRBs), and GRBs form a continuum, and thus p…
▽ More
One third of all HETE-2--localized bursts are X-Ray Flashes (XRFs), a class of events first identified by Heise in which the fluence in the 2-30 keV energy band exceeds that in the 30-400 keV energy band. We summarize recent HETE-2 and other results on the properties of XRFs. These results show that the properties of XRFs, X-ray-rich gamma-ray bursts (GRBs), and GRBs form a continuum, and thus provide evidence that all three kinds of bursts are closely related phenomena. As the most extreme burst population, XRFs provide severe constraints on burst models and unique insights into the structure of GRB jets, the GRB rate, and the nature of Type Ib/Ic supernovae. We briefly mention a number of the physical models that have been proposed to explain XRFs. We then consider two fundamentally different classes of phenomenological jet models: universal jet models, in which it is posited that all GRBs jets are identical and that differences in the observed properties of the bursts are due entirely to differences in the viewing angle; and variable-opening angle jet models, in which it is posited that GRB jets have a distribution of jet opening angles and that differences in the observed properties of the bursts are due to differences in the emissivity and spectra of jets having different opening angles. We consider three shapes for the emissivity as a function of the viewing angle theta_v from the axis of the jet: power-law, top hat (or uniform), and Gaussian (or Fisher). We then discuss the effect of relativistic beaming on each of these models. We show that observations can distinguish between these various models.
△ Less
Submitted 8 May, 2005;
originally announced May 2005.