Search

Bipolar magnetic regions (BMRs) are the cornerstone of solar variability. They are tracers of the large-scale magnetic processes that give rise to the solar cycle, shapers of the solar corona, building blocks of the large-scale solar magnetic field, and significant contributors to the free-energetic budget that gives rise to flares and coronal mass ejections. Surprisingly, no homogeneous catalog of BMRs exists today, in spite of the existence of systematic measurements of the magnetic field since the early 1970's. The purpose of this work is to address this deficiency by creating a homogenous catalog of BMRs from the 1970's until the present. For this purpose, in this paper we discuss the strengths and weaknesses of the automatic and manual detection of BMRs and how both methods can be combined to form the basis of our Bipolar Active Region Detection (BARD) code and its supporting human supervision module. At present, the BARD catalog contains more than 10,000 unique BMRs tracked and characterized during every day of their observation. Here we also discuss our future plans for the creation of an extended multi-scale magnetic catalog combining the SWAMIS and BARD catalogs.

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., Comment: 13 pages, 12 figures

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}.

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., Comment: 6 pages, 4 figures; supplemental information included (14 pages, 10 figures)

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., Comment: 12 pages, 9 figures

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.

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., Comment: Updated Author and Reviewer Information, 23 pages 17 figures

Three dimensional FLASH magneto-hydrodynamics(MHD) modeling is carried out to interpret the OMEGA laser experiments of strongly magnetized, highly collimated jets driven by a ring of 20 OMEGA beams. The predicted optical Thomson scattering spectra and proton images are in good agreement with a subset of the experimental data. Magnetic fields generated via the Biermann battery term are amplified at the boundary between the core and the surrounding of the jet. The simulation predicts multiple axially aligned magnetic flux ropes with alternating poloidal component. Future applications of the hollow ring configuration in laboratory astrophysics are discussed.

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., Comment: Main paper pp. 1-29; appendices pp. 30-84. 24 figures, 2 tables

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., Comment: Accepted for publication on Physics of Plasmas, 15 pages 12 figures

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 $\mu$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., Comment: Submitted for publication, 5 pages 4 figures with supplementary information (25 pages 9 figures)

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., Comment: 6 pages, 5 figures, 1 table; accepted by Physical Review Letters

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., Comment: Position Paper

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.

We study the magnetic properties of small-scale transients in coronal hole. We found all brightening events are associated with bipolar regions and caused by magnetic flux emergence followed by cancellation with the pre-existing and newly emerging magnetic flux. In the coronal hole, 19 of 22 events have a single stable polarity which does not change its position in time. In eleven cases this is the dominant polarity. The dominant flux of the coronal hole form the largest concentration of magnetic flux in terms of size while the opposite polarity is distributed in small concentrations. In the coronal hole the number of magnetic elements of the dominant polarity is four times higher than the non-dominant one. The supergranulation configuration appears to preserve its general shape during approximately nine hours of observations although the large concentrations in the network did evolve and were slightly displaced, and their strength either increased or decreased. The emission fluctuations seen in the X-ray bright points are associated with reoccurring magnetic cancellation in the footpoints. Unique observations of an X-ray jet reveal similar magnetic behaviour in the footpoints, i.e. cancellation of the opposite polarity magnetic flux. We found that the magnetic flux cancellation rate during the jet is much higher than in bright points. Not all magnetic cancellations result in an X-ray enhancement, suggesting that there is a threshold of the amount of magnetic flux involved in a cancellation above which brightening would occur at X-ray temperatures. Our study demonstrates that the magnetic flux in coronal holes is continuously recycled through magnetic reconnection which is responsible for the formation of numerous small-scale transient events. The open magnetic flux forming the coronal-hole phenomenon is largely involved in these transient features., Comment: 19 pages, 18 figures, A&A in press

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., Comment: 22 pages, 8 figures; ApJ

The discovery of the accelerated expansion of the universe using Type Ia supernovae (SNe Ia) has stimulated a tremendous amount of interest in the use of SNe Type Ia events as standard cosmological candles, and as a probe of the fundamental physics of dark energy. Recent observations of SNe Ia have indicated a significant population difference depending on the host galaxy. These observational findings are consistent with SNe Ia Ni-56 production in star-forming spiral galaxies some 0.1 solar masses higher - and therefore more luminous than in elliptical galaxies. We present recent full-star, 3D simulations of Type Ia supernovae which may help explain the nature of this systematic variation in SNe Ia luminosities, as well as the nature of the Ia explosion mechanism. These insights may in turn eventually shed light on the mystery of dark energy itself., Comment: 10 Pages, 3 Figures, Submitted to Proceedings of The Ninth Asia-Pacific International Conference on Gravitation and Astrophysics (ICGA9) held June 29 to July 3, 2009, at Huazhong University of Science & Technology in Wuhan, China.

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., Comment: 13 pages, 9 figures, accepted for the publication in PASJ (minor corrections)

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., Comment: 20 pages, 10 figures

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 \nu F_\nu spectrum, and the photon index (\Gamma) derived from a simple power-law model. The Epeak - \Gamma relation, assuming the typical smoothly broken power-law spectrum of GRBs, is \log Epeak = 3.258 - 0.829\Gamma (1.3 < \Gamma < 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., Comment: 27 pages, 31 figures, accepted for publication in ApJ

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)., Comment: 20 pages, 14 figures, 3 tables (editor version: 15 pages, 14 figures, tables 2 and 3 only available in electronic form at http://www.aanda.org)

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 ($\sigma_\theta/\sigma_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., Comment: 7 pages (emulateapj), 9 figures. Expanded figures and text to match accepted version

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 $\beta$-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., Comment: 14 pages (emulateapj), 13 figures. Updated text to match accepted version

We present a collection of eight data sets, from state-of-the-art experiments and numerical simulations on turbulent velocity statistics along particle trajectories obtained in different flows with Reynolds numbers in the range $R_\lambda \in [120:740]$. Lagrangian structure functions from all data sets are found to collapse onto each other on a wide range of time lags, revealing a universal statistics, and calling for a unified theoretical description. Parisi-Frisch Multifractal theory, suitable extended to the dissipative scales and to the Lagrangian domain, is found to capture intermittency of velocity statistics over the whole three decades of temporal scales here investigated., Comment: 5 pages, 1 figure; content changed, references updated

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., Comment: 50 pages, 22 figures, Accepted for publication in ApJ

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., Comment: 14 pages, 10 figures, accepted to the Astrophysical Journal

Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work., Comment: In press for Astrophys. J. 2007

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., Comment: Accepted by ApJ

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., Comment: 7 pages, 4 figures

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., Comment: 24 pages, 6 figures, accepted for publication in ApJ; minor revision

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., Comment: 21 pages, 24 figures, to appear in ApJ

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., Comment: 10 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial., Comment: Matches accepted PRD version. SDSS data, likelihood code, Markov chains and ppt figures available at http://space.mit.edu/home/tegmark/sdss.html 36 journal pages, 25 figs. CosmoMC plugin at http://cosmologist.info/cosmomc/

Thanks to the rapid position notice and response by HETE-2 and Swift, the X-ray afterglow emissions have been found for four recent short gamma-ray bursts (GRBs; GRB 050509b, GRB 050709, GRB 050724, and GRB 050813). The positions of three out of four short GRBs are coincident with galaxies with no current or recent star formation. This discovery tightens the case for a different origin for short and long GRBs. On the other hand, from the prompt emission point of view, a short GRB shows a harder spectrum comparing to that of the long duration GRBs according to the BATSE observations. We investigate the prompt emission properties of four short GRBs observed by Swift/BAT. We found that the hardness of all four BAT short GRBs is in between the BATSE range for short and long GRBs. We will discuss the spectral properties of short GRBs including the short GRB sample of Konus-Wind and HETE-2 to understand the hard nature of the BATSE short GRBs., Comment: 5 pages, 8 figures. To appear in the proceedings of "Gamma Ray Bursts in the Swift Era", eds. S. Holt, N. Gehrels & J. Nousek

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)., Comment: 60 pages, 19 figures, submitted to ApJ; added new references and corrected typos

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., Comment: 4 pages, 4 figures. Submitted to ApJL

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., Comment: 19 pages, 4 figures. Submitted to Nature

In 2000, Lamb and Reichart predicted that gamma-ray bursts (GRBs) and their afterglows occur in sufficient numbers and at sufficient brightnesses at very high redshifts (z > 5) to eventually replace quasars as the preferred probe of element formation and reionization in the early universe and to be used to characterize the star-formation history of the early universe, perhaps back to when the first stars formed. Here we report the discovery of the afterglow of GRB 050904 and the identification of GRB 050904 as the first very high redshift GRB. We measure its redshift to be 6.39(+0.11,-0.12), which is consistent with the reported spectroscopic redshift (6.29 +/- 0.01). Furthermore, just redward of Ly-alpha the flux is suppressed by a factor of three on the first night, but returns to expected levels by the fourth night. We propose that this is due to absorption by molecular hydrogen that was excited to rovibrational states by the GRB's prompt emission, but was then overtaken by the jet. Now that very high redshift GRBs have been shown to exist, and at least in this case the afterglow was very bright, observing programs that are designed to capitalize on this science will likely drive a new era of study of the early universe, using GRBs as probes., Comment: Submitted to Nature, 11 pages, 3 figures, LaTeX

[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., Comment: White paper submitted to the Dark Energy Task Force. 37 pages

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0

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_{\gamma}$ is narrowly distributed; the Amati relation holds between $E_{iso}$ and $E_{pk}$; and the Ghirlanda relation holds between $E_{\gamma}$ 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_{\gamma}$ 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., Comment: 71 page, 33 figures. Submitted to ApJ

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., Comment: Accepted to The Astrophysical Journal, 30 pages, 5 Figures, LaTeX

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., Comment: 8 pages, 7 figures, 1 table. Invited review talk at the 4th Workshop Gamma-Ray Bursts in the Afterglow Era, Rome,18-22 October 2004. Editors: L. Piro, L. Amati, S. Covino, and B. Gendre. Il Nuovo Cimento, in press