Research Article  Open Access
Numerical Experiments for Nuclear Flashes toward Superbursts in an Accreting Neutron Star
Abstract
We show that the superburst would be originated from thermonuclear burning ignited by accumulated fuels in the deep layers compared to normal Xray bursts. Two cases are investigated for models related to superbursts by following thermal evolution of a realistic neutron star: helium flash and carbon flash accompanied with many normal bursts. For a helium flash, the burst shows the long duration when the accretion rate is low compared with the observation. The flash could become a superburst if the burning develops to the deflagration and/or detonation. For a carbon flash accompanied with many normal bursts, after successive 2786 normal bursts during 1.81 × 10^{9} s, the temperature reaches the deflagration temperature. This is due to the produced carbon which amount reaches to ≈0.1 in the mass fraction. The flash will develop to dynamical phenomena of the deflagration and/or detonation, which may lead to a superburst.
1. Introduction
Type I Xray bursts have been identified to the thermonuclear explosions on the surface region of accreting neutron (compact) stars. As a consequence, the phenomenon has been studied from both nuclear reactions and nuclear structure inside the compact stars. However, there still remain many uncertainties concerning the elementary processes associated with the bursts [1–3]. Superbursts have been detected from 13 Xray bursters by BeppoSAX and RXTE (see, e.g., Table 2 in [4]). In particular, 4U 1636536 exhibited four superbursts, where the shortest recurrence time is 1.5 years (http://www.astronomerstelegram.org/?read=2140) [5, 6]. Clearly, the light curve consists of a fast rise and slower power lawlike decay [2, 7]. The spectrum hardens during the rise phase to the maximum in luminosity, whereas it softens in the decay phase. This is also reflected in the spectral fits to the timeresolved preburst subtracted from Xray spectra. Each burst has energy of ergs and duration of a few hours. They are usually best described in terms of a blackbody model. The effective temperature increases and decreases during the rise and decay phase, respectively. These superbursts are 1000 times luminous and 1000 times long in the duration compared with the normal bursts though the spectral evolution is similar.
Even now, quantitative explanation and/or numerical simulation of superbursts using the stellar evolution code are limited. For example, Keek and Heger [8] do not selfconsistently produce the carbon from hydrogen/helium burning in their calculation but instead accrete the carbon directly onto the neutron star, bypassing the hydrogen/helium burning stages. The superbursts last too long and their energy release is too much to be explained in terms of unstable burning of hydrogen/helium so far considered [9]. Moreover, regular normal Xray bursts are observed before the occurrence of the superburst [10] that includes the precursor burst [11]. The long rise and decay times of superbursts are consistent with the model of unstable burning in the deep layer below the hydrogen/helium burning region. Therefore, it has been suggested that unstable burning of carbon is the origin of the superbursts [9, 12].
If the accreted material onto the neutron star is pure helium, carbon can be produced when helium is burned stably [9]. Unstable helium burning often involves alpha captures on carbon and carbon does not remain much. This would apply to the helium accretor 4U 182030 that shows long periods of high Xray intensity during which no burst occurs, which is consistent with a period of stable helium burning. Note that unstable carbon burning only reproduces the observed feature in superbursts when we take into account neutrino losses and significant heat flux transported from deeper into the accreting layer of the neutron star [9]. Cumming et al. [2] also show that the observed superburst energy is around ergs and more or less independent of ignition depth because of neutrino emission of the excess energy. While Cumming [13] expected recurrence times of the order of 12 years, Strohmayer and Brown [9] obtained a recurrence time of about 10 years.
If the accreted material onto the neutron star is a mixture of hydrogen and helium, either unstable or stable burning of hydrogen/helium can produce carbon. While the amount obtained by numerical calculations has been only limited [14, 15] after bursts, carbon is much more readily produced in stable burning [1]. In the observed sources of superbursts, normal XRBs have been observed with a mean rate of about 3 times per day during the period of the observation [16, 17]. This indicates that at least some amounts of the accreted material have been burned stably before a superburst. Furthermore, the superburst from 4U 125469 indicates that much of the accreted fuel burns stably [18]. Detection of superbursts at near Eddington accretion rate would reveal the relation between the recurrence time and remained nuclear fuels, which is studied by using the parameter [19].
Cumming and Bildsten [12] suggested that a small amount of carbon could be enough to trigger a thermonuclear runaway with energy comparable to the superburst if carbon resides in a bath of heavy elements. These heavy elements are the products of unstable burning through the rpprocess during the mixed hydrogen/helium burning of XRBs [14, 20]. In this case, the superburst recurrence time would depend on accretion rates, being in the order of a few decades, a year to a decade, or a week to a month, according to the accretion rate which is about 0.1, 0.3, or 1 times the Eddington accretion rate. On the other hand, it is proposed from the analysis of the photospheric radius expansion that an accretion rate may change during the burst by some factors [21]. Changes in accretion rates are also considered from recent study of the outburst of a transient Xray binary [22]. Therefore, it is reasonable to include accretion rate variations in numerical calculations because the observations indicate that the luminosity of accreting neutron stars is variable. It has been suggested that the high temperature reached during a superburst induced the photodisintegration [23]. As a consequence, they got energy generations comparable to those due to the carbon induced superburst.
Another scenario was proposed by Kuulkers et al. [17]. They suggested that hydrogen left after the burning in the hydrogen/helium layer is reignited by the electron capture that is followed by successive captures of neutrons by heavy nuclei occurred in the deeper layer. The recurrence time of superbursts is estimated to be less than one year [10, 17]. However, large amounts of hydrogen should be remained after bursts to explain the energy release in superbursts. Recent calculations have revealed that hydrogen is completely depleted after the hydrogen/helium burning [14, 15]. On the other hand, an exotic process of the diquark pair formation was proposed to explain superbursts [24].
On the other hand, significant progress has been done so far concerning the construction of model to study the phenomena of neutron star; two dimensional hydrodynamical model calculations of Xray bursts [25, 26], study of the propagation of deflagration wave of a rapidly rotating neutron star [27], magnetorotational study [28, 29], and rotating hotspot model examination of rotating neutron star [30] concerning burst ocillations. Unfortunately, it would be insufficient to study Xray bursts in details beyond spherically symmetric model. Since there is no model in the literature that selfconsistently calculates the production of carbon in hydrogen/helium and its subsequent ignition as superburst, we perform evolutionary calculations adopting spherically symmetric models of accreting neutron stars with important physical processes included [31].
In Section 2 physical inputs and our evolutionary code are explained. A helium flush model is presented in Section 3 related to superbursts. In Section 4, we present a model of carbon flash accompanied with many normal bursts of H/He combined burnings and show the possibility of a superburst. Concluding remarks are given in Section 5.
2. Evolution Code of a Neutron Star
The general relativistic evolutionary equations of spherical stars in hydrostatic equilibrium as formulated by Throne [32] are written as where The basic quantities are defined as follows: : rest mass density, : temperature, : abundance of the th particle, : proper mass inside the radius , : Schwarzschild time coordinate (proper time at a distant observer), : total mass inside the radius , : gravitational potential, : local luminosity, : pressure, : total nongravitational massenergy density in mass units, : specific entropy, : heating rate by nuclear burning, : cooling rate by escaping neutrinos, : nuclear reaction rate for the th particle, and (): the radiative (adiabatic) temperature gradient. In the accretion layer, the Eulerian coordinate (the mass fraction coordinate with changing mass ) is used, which is the most suitable method for computations of stellar structure when stellar mass varies [33]. In our calculation, the mixing length theory of convection with the mixing length equal to the pressure scale height has been used [34]. The radiative zero boundary conditions are imposed at the outer boundary. An outermost mesh point, which is regarded as the photosphere, is given at [35].
The above set of general relativistic equations for the evolution of spherically symmetric stars has been solved by using a Henyeytype numerical scheme of implicit method. We adopt the evolution code of a spherically symmetric neutron star [35, 36]. The star is divided into 266 meshes of the Lagrange masscoordinate. The gravitational mass and radius of the neutron star are initially set to be 1.3 and 8.1 km, respectively. Then we have log for the gravitational acceleration at the surface of the neutron star . The depth of the accretion layer is ~10 m and consists of 170 meshes. In the region of the combined hydrogen/helium burning, the interval between meshes is typically 56 cm. Physical inputs and an approximation network are almost the same as adopted by Fujimoto et al. [31] except for the changes of reaction rates [37] and alphanetwork used in Section 3. It should be noted that our approximation network has been constructed based on the onezone calculation with use of a large network for the bursts [38, 39]. Our network includes 16 nuclei: , , , , , , , , , , , , , , , and . This network can be used till the hydrogen decreases until around 10% in mass fractions [39].
Although input physics has been described in detail [35], we briefly write it again. This is because input physics is crucial to determine the neutron star properties and it is closely related to the study of elementary processes of high density matter. As for the equation of state of outer part of the neutron star , an ideal gas plus radiation is assumed with the electron degeneracy and the Coulomb liquid correction included [40]. For the inner part, the equation of state has been constructed by Richardson et al. [41] based on Canuto [42, 43]. Neutrino emissivities include bremsstrahlung of nucleons [44] and electronion [45], and electronpositron pair, photo, and plasmon processes [46]. Opacities for iron and neutronrich material include the radiative opacities for by Malone [47] and the thermal conductivities [48–50]. For the opacities of lighter elements the analytical approximations by Iben [51] are adopted to the radiative ones [52, 53], to the thermal conductivity for a nonrelativistic electron gas of [54], and to that for a relativistic one [55], respectively. Screening effects of nuclear reactions are taken from Dewitt et al. [56]. In the present paper, we include the crustal heating [57], where is the number of the reaction, is the mass accretion rate in units of , and is the number of the effective heat energy per nucleon in MeV for the th reaction which is tabulated in their paper.
We assume that material is accreted with the same entropy at the stellar surface, neglecting the surface effects caused by the accretion flow. Within the framework of the spherical symmetry, the kinetic energy of the falling material has little influence on the structure in the layer as deep as the burning shell, since the radial motion will be dissipated in the surface layers [35].
Initial models have been constructed through the continuous accretion ( constant) without nuclear burning until the steady state is achieved, where the nonhomologous part of the gravitational energy release vanishes [35].
3. Simple Features of Models Concerning Superbursts—Helium Flash Model
Superbursts observed in 4U 182030 are considered to be accreted by pure helium, while 4U 1636536, 4U 173544, and 4U 1254690 are accreted by hydrogen/helium [58]. The difference of the accretion matter may affect the mechanism of superbursts. Therefore, we first examine rather simplified models of superbursts, which are triggered by a helium flash. We note that deep helium ignition towards superbursts was discussed by Kuulkers et al. [59] on the Xray binary 4U 0614+091. For the carbon flash, we also investigate realistic situation accompanied with many normal bursts in Section 4.
The superburst of 4U 182030 was observed in the duration of 2.5 hr, where the burst energy is ergs [9], and the peak luminosity is [17]. To follow a pure helium flash, we use an alphanetwork consisted of alphanuclei up to , where the nuclear reaction rates are taken from Angulo et al. [37]. This network includes not only reactions but also reactions; the latter in competition with the former is assumed to follow reactions instantaneously. From Figure 1, we can see that although is consistent with the observation for a high accretion rate, , and intermediate rate, , the duration of the burst is too short to explain the burst energy. It should be noted that attains the Eddington limit as described in Hanawa and Fujimoto [36]. On the other hand, for a low accretion rate, , we have obtained the burst energy around ergs which lasts more than 3 hr (see Figure 1). It is noted that the ignition pressures in units of dyn for the above accretion rates are , , and , respectively. Since of the superburst in 4U 182030 is estimated to be several times [13], we recognize that simple models based on the single burst with an accretion rate assumed for the accreting neutron star are inconsistent with the observation (see Section 4). However, since the ultracompact source 4U 0614+91 likely accretes helium at with a superburst [59], we could carefully study the helium flash.
Figure 2 shows the temperature distribution for against the density of the initial state (dotted line) and that of the stage at the maximum nuclear luminosity (solid line), where denotes the nuclear luminosity. The dashed line indicates the ignition curve of the reaction for the helium mass fraction , with the nuclear energy generation rate of reaction and the radiative energy loss rate [14]. Although this criterion should be carefully reconsidered [60], we adopt the present one for simplicity. The dotted line is the deflagration temperature defined by equating the dynamical time scale and the nuclear heating time scale () with where is the pressure scale height, is the sonic velocity, and is the specific heat under the constant pressure. It is remarkable that the temperature in the layers of log exceeds log for : the flash may become the deflagration. It needs to perform a dynamical calculation to elucidate how the deflagration develops inside the accretion layers. Although we cannot represent the proper consistent with the observation of , the helium flash in low accretion rate could become a possible site of the superburst.
4. Model Accompanied with Normal Bursts of Combined Hydrogen and Helium Burnings
In the previous section, we have shown that simple helium burst models cannot explain the observed superbursts. New model is needed to produce enough carbon for a superburst to occur. Therefore, we present a sequence of calculations until the amount of carbon increases enough. We first adopt an accretion rate of . And to acquire more carbon, we only change it to , which save the computational time. To increase the temperature in the distribution, we change the accretion rate to the first one. Furthermore, we also raise the crustal heating by a factor of ten to save the computational time.
4.1. Gross Features toward the Carbon Flash
Let us make an initial model to simulate a superburst with the accretion rate of . This accretion rate is considered to be in the reasonable range of with for the observed sources [61]. Note that superbursts at nearEddington rate are suggested for GX 17+ 2 [19]. With the nuclear burning suppressed, we construct an initial temperature distribution by continuous accretion with this specified . Mass fractions in the accretion matter are assumed to be H (73.0), (25.0), (0.7), and (1.3). The bottom of the accretion layer is set to be pure , which is equivalent to in the present purpose. The steady state concerning the accretion rate is assumed; that is, ; this means that accreted matter increases the mass of the layer of pure ; that is, accreting matter becomes promptly and omits nuclear burning processes. At the end, we obtain the isothermal temperature distribution of log .
We repeatedly calculate normal bursts (combined hydrogen/helium burning) with use of the evolutionary code which includes the approximation reaction network. Though hydrogen consumption after a flash with use of this network might have been underestimated by ~10% for log [14], remaining hydrogen is less than ~1% in the bottom layer of normal bursts due to the convection. Therefore, in the deep region related to a superburst (log ), where rpprocess does not work anymore, our network can also be used except for detailed abundance distribution. Figure 3 shows the bursts from the beginning of accretion to s with the accretion rate . In Figure 4, we show the twelve bursts till s which should be compared with those of Woosley et al. [15] having the solar initial composition and . In view of the fact that accretion rates are different from each other, both cases produce the regular bursts. In our case, these regular bursts continue s and then the recurrence time becomes longer. Since there exists significant amount of produced materials inside the deep region, the heat transported to the inner part of the neutron star results in the lengthened recurrence intervals of bursts.
The totally calculated evolutionary time is s and the total number of normal bursts amounts to 2786. The time sequence of the bursts is illustrated in Figure 5. The time interval and the number of bursts are given in Table 1 for the individual period specified in terms of . Figures 6–8 show the light curve during the intervals (3–11), (11–15), and s. Critical bursts are marked by 1–7, whose profiles are discussed in Section 4.2. Important epochs leading to a superburst are specified as follows. The epoch is 930 s after the burst at s during “period 1,” and the epoch is 854 s after the burst at s during “period 2.” The epoch corresponds to the end of the last burst 7 and is just before the ignition which is 5500 s after the epoch .


(a)
(b)
Figure 9 shows temperature profiles against density and pressure, respectively, at the epochs . The ignition curve and deflagration line of the reaction are also shown in Figure 9. We note that in the context of detonations in superbursts, comparing to the dynamical time is discussed by Weinberg and Bildsten [62]. The two convexities near log are due to the effect of the unstable hydrogen/helium burning. The energy generation rate is shown in Figure 10. HotCNO cycle, (, ) (p, ) (, ) (p, ) (p, ) (p, ), has produced the energy of for log [38]. Around log , the energy generation rates of , , and before the ignition become small because there remains small amount of fuel as the result of both the rpprocess and steady burning of the reaction. It should be noted that the curve of for falls down steeply around log because of the significant decrease in abundances due to the sudden development of the convection as illustrated in Figure 14.
For “period 1,” the accretion rate is set to be from the beginning of the accretion to s. The number of bursts occurred in this interval is 1237. The temperature distribution inside the deep accreted layers (log and log ) remains isothermal of log (see the temperature distribution at s in Figure 9). We can recognize that unstable combined hydrogen/helium burning has been generated for log and log . Figure 11 shows the composition distribution at s (). The range of the pressure is equal to that in Figure 9. Around log , the rpprocess produces both and . Though there is no remained fuel of hydrogen [14], remained helium produces for log (see the solid line in Figure 11) due to the steady burning. The increase in is ascribed to the numerical diffusion of the initial distribution and convective mixing at the beginning of the accretion for log ; thus, we can consider that the mass fraction of in should be added to . For “period 1,” the carbon burning is stable (log , log and ) in the sense that the increase in does not reach the ignition curve. In the bottom of the accretion layer (log ), carbon sometimes burns intensively and the temperature of the layer has increased to log . The amount of carbon is too small () to trigger the nuclear flash. The nuclear energy generation for is supplied by helium burning for log , while it is supplied by reactions such as and for log (Figure 10). The crustal heating help to increase again for log .
After s, the accretion rate is artificially changed to . This results in the increase in abundance. We keep this accretion rate in the interval of s, and the 1416 bursts occur during this interval (period 2). The temperature distribution at s for this accretion rate is shown by the dotted line in Figure 9 (denoted by ). Figure 12 shows the composition distribution at this stage . For “period 2,” hydrogen is consumed completely in the region of log and the remained helium increases appreciably. The rpprocess first produces and Zn and afterwards also and . The appreciable decrease in corresponds to the dip of the composition distribution around log in Figure 12.
We have changed again the accretion rate to and keep it for s (period 3). The number of bursts is 133 in this interval. Though the temperature distribution has become lower compared with “period 1,” if we continue the calculations of normal bursts, the temperature distribution should be recovered to that in “period 1.” Therefore, we increase the crustal heating by a factor of 10 to save the computational time. Although this artificial change appears to be unreal, from the point of nuclear physics, physical process concerning the crustal heating has been rather uncertain; the rate is often multiplied some factors (see, e.g., [8, 63]). We note that the heating rates (ergs/g) change up to 20 times for the pressures shown in tables of Haensel and Zdunik [57], and we have implemented the crustal heating using the tables as the heating source.
Figure 13 shows the composition distribution at the stage (). As in the case of “period 1,” hydrogen is consumed in the bottom region of the hydrogen/helium burning and the remained helium forms C continuously below the region. After 5500 s () from the end () of the last burst (7), we obtain the carbon flash. Finally, the reaction is successfully ignited in log at the last epoch ( s). The composition distributes uniformly for log due to the convection as shown in Figure 14. We note that the flash does not yet develop to the peak of the total nuclear energy generation rate (Figure 19).
4.2. Characteristic Features of Key Bursts Associated with Many Normal Bursts
Let us discuss the profiles and/or sequences of the bursts for each period. In Figure 6(b), there is a significant interval at s. We call it the “long interval 1.” We pay attention to five bursts (1–5) before the “long interval 1” and the first burst (6) after the interval.
Five bursts 1–5 extracted from Figure 6 are shown in Figure 15. “Long interval 1” begins at s. The light curve of the burst 1 is shown in Figures 16 and 17, where we set the time to be zero when the burst 1 starts. The roman numerals specify the epochs: I: just before the convection, II: peak of the total nuclear energy generation rate , and III: 1/10 of the maximum luminosity . In Figure 16, the epoch “IV” for the burst 5 corresponds to 220 s after the onset of the burst.
The shapes of the light curves and the products for the bursts 1–4 and 6 are almost same. The burst 5 is the longest burst compared with other bursts of 1–4 (Figure 16) and 6. At “IV” of the burst 5, the hydrogen burns in log , and is produced where the convection occurs for log that leads to the exhaustion of hydrogen. The heat from this burning region of the convection is the reason for the long burst as seen in Figure 16, and the consumption of fuels leads to the long interval after 5 (Figures 6(b) and 15). This feature can be seen in other sequences of bursts. At the epoch “IV” of the burst 3 (2240 s after the onset of the burst), the convection also occurs between log where hydrogen is completely consumed. Therefore, as seen in Figure 15 the interval to the next burst 4 becomes long compared to other intervals except for the “long interval 1,” because the accumulated hydrogen has been consumed. The same situation also occurs for the rather long interval s in Figure 7.
The luminosity (Figure 17) in the burst of is a little small compared with the bursts 1–5 and 6 of . Though the products of , , , and are the same as those in the other bursts, the amount of heavy elements like and are produced more than those of the other bursts. The temperature distribution of the inner accretion layer (log ) is lower than that of the other bursts with (see Figure 9, where in the stage of for log is low due to the conduction). We can consider that due to the increased conduction the luminosity in the bursts decreases because the heat from the combined hydrogen/helium burning has flowed into the core. This phenomena can be understood from the idea of watershed introduced by Fujimoto et al. [35]. As illustrated in Figure 17, the last burst 7 before the carbon ignition is weaker in strength than that of the burst , because the hydrogen/helium burning is ignited at the deeper region (log ), which can be seen from the ash of the rpprocess in Figure 13, compared to the other bursts (log ) as seen in Figures 11 and 12.
The light curve and the total nuclear energy generation rate, , from the burst 7 to the onset of the carbon flash are shown in Figures 18 and 19. Just after the burst, , the convection occurs at the narrow region of log (Figure 13). At s in Figure 19, the convection begins to spread to lower pressure region, and it extends from log to 26 as shown in Figure 14. We note that at is not attained to the . However, since the temperature in log has increased beyond the deflagration temperature defined by as seen from Figure 9, we can infer that the flash should develop to a dynamical phenomena of deflagration. Finally, we stop the calculation of the flash due to the numerical difficulty in the assumption of the hydrostatical equilibrium; nevertheless, we insist that this flash could become a superburst.
5. Concluding Remarks
We have presented the two cases for the models of superbursts: single helium flash and carbon flash accompanied with many normal bursts. For the helium flash, the burst has a long duration time although the accretion rate is different from the observation. We suggest that the helium flash could originate the superburst if the burning develops to the deflagration and/or detonation. For a single carbon flash, the temperature does not reach the deflagration temperature due to the limit of heat conduction.
For the carbon flash accompanied with many normal bursts, we carried out the successive 2786 normal bursts up to the time of s using the observed accretion rates. We showed the profiles of the several normal bursts and the onset of the carbon flash after the normal bursts. Since there exists significant amount of produced elements in the deeper than the region where hydrogen/helium burning occurred, the heat transported to the core by them lengthens the recurrence intervals of bursts. We have also shown that the recurrence time becomes longer than other bursts when the little burning occurred for the bottom of the hydrogen/helium burning layer after the main burning had ceased. This leads to the consumption of nuclear fuels of H and He and increases the products of rpprocess. We have recognized that “long interval 1” is a remarkable case for the delay to the next burst due to the above reasons. Before the carbon flash, we find various profiles of the light curve and different intervals between bursts; Quantitative comparison with the observations should shed a new insight into the Xray burst research. Since normal bursts have been observed before the superbursts, our scenario is consistent with the observations. We conclude that a carbon flash accompanied with many normal bursts should trigger a superburst. It is interesting to apply the above idea to a superburst of a helium accretor like 4U182030.
For the helium flash shown in Section 3, accretion rate to produce a superburst is rather low compared with observations. If we repeat normal helium flashes with the observed as demonstrated in Section 4.2, a superburst could be triggered in a deep region of the accreting layers. Finally, we remark that to get a detailed history toward the superburst, it is needed to calculate many bursts using the code of the neutron star evolution with a large nuclear reaction network coupled.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 H. Schatz, L. Bildsten, A. Cumming, and M. Ouellette, “Nuclear physics in normal Xray bursts and superblasts,” Nuclear Physics A, vol. 718, pp. 247–254, 2003. View at: Publisher Site  Google Scholar
 A. Cumming, J. Macbeth, J. J. M. in't Zand, and D. Page, “Long type I xray bursts and neutron star interior physics,” The Astrophysical Journal, vol. 646, pp. 429–451, 2006. View at: Google Scholar
 A. Parikh, J. José, G. Sala, and C. Iliadis, “Nucleosynthesis in type I Xray bursts,” Progress in Particle and Nuclear Physics, vol. 69, no. 1, pp. 225–253, 2013. View at: Publisher Site  Google Scholar
 L. Keek, A. Heger, and J. J. M. in't Zand, “Superburst models for neutron stars with hydrogen and heliumrich atmospheres,” The Astrophysical Journal, vol. 752, no. 2, 150, 2012. View at: Publisher Site  Google Scholar
 R. Wijnands, “Recurrent very long type I Xray bursts in the lowmass Xray binary 4U 1636–53,” The Astrophysical Journal Letters, vol. 554, no. 1, article L59, 2001. View at: Publisher Site  Google Scholar
 T. E. Strohmayer and C. B. Markwardt, “Evidence for a millisecond pulsar in 4U 163653 during a superburst,” Astrophysical Journal Letters, vol. 577, no. 1, pp. 337–345, 2002. View at: Publisher Site  Google Scholar
 A. Cumming and J. Macbeth, “The thermal evolution following a superburst on an accreting neutron star,” Astrophysical Journal Letters, vol. 603, no. 1, pp. L37–L40, 2004. View at: Publisher Site  Google Scholar
 L. Keek and A. Heger, “Multizone models of superbursts from accreting neutron stars,” Astrophysical Journal, vol. 743, no. 2, article 189, 2011. View at: Publisher Site  Google Scholar
 T. E. Strohmayer and E. F. Brown, “A remarkable 3 hour thermonuclear burst from 4U 1820–30,” The Astrophysical Journal, vol. 566, pp. 1045–1059, 2002. View at: Google Scholar
 E. Kuulkers, “A superburst from GX 3+1,” Astronomy and Astrophysics, vol. 383, no. 1, pp. L5–L8, 2002. View at: Publisher Site  Google Scholar
 L. Keek, “Photospheric radius expansion in superburst precursors from neutron stars,” The Astrophysical Journal, vol. 756, article 130, 2012. View at: Publisher Site  Google Scholar
 A. Cumming and L. Bildsten, “Carbon flashes in the heavyelement ocean on accreting neutron stars,” Astrophysical Journal Letters, vol. 559, no. 2, pp. L127–L130, 2001. View at: Publisher Site  Google Scholar
 A. Cumming, “Models of type I XRAY bursts from 4U 182030,” The Astrophysical Journal, vol. 595, no. 2, pp. 1077–1085, 2003. View at: Publisher Site  Google Scholar
 O. Koike, M. Hashimoto, R. Kuromizu, and S. Fujimoto, “Final products of the rpprocess on accreting neutron stars,” The Astrophysical Journal, vol. 603, no. 1, pp. 242–251, 2004. View at: Publisher Site  Google Scholar
 S. E. Woosley, A. Heger, and A. Cumming, “Models for type I Xray bursts with improved nuclear physics,” The Astrophysical Journal Supplement Series, vol. 151, no. 1, p. 75, 2004. View at: Publisher Site  Google Scholar
 R. Cornelisse, E. Kuulkers, J. J. M. in't Zand, F. Verbunt, and J. Heise, “A fourhours long burst from Serpens X1,” Astronomy & Astrophysics, vol. 382, pp. 174–177, 2002. View at: Publisher Site  Google Scholar
 E. Kuulkers, J. J. M. In't Zand, and M. H. van Kerkwijk, “A halfaday long thermonuclear Xray burst from KS 1731260,” Astronomy and Astrophysics, vol. 382, no. 2, pp. 503–512, 2002. View at: Publisher Site  Google Scholar
 J. J. M. In't Zand, E. Kuulkers, F. Verbunt, J. Heise, and R. Cornelisse, “A superburst from 4U 125469,” Astronomy and Astrophysics, vol. 411, no. 3, pp. L487–L491, 2003. View at: Publisher Site  Google Scholar
 J. J. M. in't Zand, R. Cornelisse, and A. Cumming, “Superbursts at nearEddington mass accretion rates,” Astronomy & Astrophysics, vol. 426, no. 1, pp. 257–265, 2004. View at: Publisher Site  Google Scholar
 H. Schatz, L. Bildsten, A. Cumming, and M. Wiescher, “The rapid proton process ashes from stable nuclear burning on an accreting neutron star,” Astrophysical Journal Letters, vol. 524, no. 2, pp. 1014–1029, 1999. View at: Publisher Site  Google Scholar
 H. Worpel, D. K. Galloway, and D. J. Price, “Evidence for accretion rate change during type i Xray bursts,” Astrophysical Journal, vol. 772, no. 2, article 94, 2013. View at: Publisher Site  Google Scholar
 A. Bahramian, C. O. Heinke, and G. R. Sivakoff, “Discovery of the third transient Xray binary in the galactic globular cluster Terzan 5,” The Astrophysical Journal, vol. 780, no. 2, article 127, 2014. View at: Publisher Site  Google Scholar
 H. Schatz, L. Bildsten, and A. Cumming, “Photodisintegrationtriggered nuclear energy release in superbursts,” Astrophysical Journal Letters, vol. 583, no. 2, pp. L87–L90, 2003. View at: Publisher Site  Google Scholar
 M. Sinha, M. Dey, S. Ray, and J. Dey, “Super bursts and long bursts as surface phenomena of compact objects,” Monthly Notices of the Royal Astronomical Society, vol. 337, no. 4, pp. 1368–1372, 2002. View at: Publisher Site  Google Scholar
 C. M. Malone, A. Nonaka, A. S. Almgren, J. B. Bell, and M. Zingale, “Multidimensional modeling of type I Xray bursts. I. Twodimensional convection prior to the outburst of a pure ^{4}He accretor,” The Astrophysical Journal, vol. 728, no. 2, article 118, 2011. View at: Publisher Site  Google Scholar
 C. M. Malone, M. Zingale, A. Nonaka, A. S. Almgren, and J. B. Bell, “Multidimensional modeling of type I Xray bursts. II. Twodimensional convection in a mixed H/He accretor,” The Astrophysical Journal, vol. 788, article 115, 2014. View at: Google Scholar
 Y. Cavecchi, A. L. Watts, J. Braithwaite, and Y. Levin, “Flame propagation on the surfaces of rapidly rotating neutron stars during type I Xray bursts,” submitted, http://arxiv.org/abs/1212.2872. View at: Google Scholar
 R. V. E. Lovelace, A. K. Kulkarni, and M. M. Romanova, “Torsional magnetic oscillations in type I Xray bursts,” The Astrophysical Journal, vol. 656, no. 1, pp. 393–398, 2007. View at: Publisher Site  Google Scholar
 D. J. B. Payne and A. Melatos, “Magnetic burial and the harmonic content of millisecond oscillations in thermonuclear Xray bursts,” The Astrophysical Journal, vol. 652, no. 1, p. 597, 2014. View at: Publisher Site  Google Scholar
 R. Artigue, D. Barret, F. K. Lamb, K. H. Lo, and M. C. Miller, “Testing the rotating hotspot model using Xray burst oscillations from 4U 1636−536,” Monthly Notices of the Royal Astronomical Society, vol. 433, no. 1, pp. L64–L68, 2013. View at: Publisher Site  Google Scholar
 M. Fujimoto, T. Hanawa, I. Iben Jr., and M. B. Richardson, “Thermal evolution of accreting neutron stars. II—long Xray bursts as a probe into the interior,” The Astrophysical Journal, vol. 315, pp. 198–208, 1987. View at: Publisher Site  Google Scholar
 K. S. Thorne, “The relativistic equations of stellar structure and evolution,” Astrophysical Journal, vol. 212, pp. 825–831, 1977. View at: Publisher Site  Google Scholar
 D. Sugimoto, K. Nomoto, and Y. Eriguchi, “Stable numerical method in computation of stellar evolution,” Progress of Theoretical Physics Supplement, vol. 70, pp. 115–131, 1981. View at: Publisher Site  Google Scholar
 M. Y. Fujimoto and D. Sugimoto, “Helium shell flashes and evolution of accreting white dwarfs,” The Astrophysical Journal, vol. 257, pp. 291–302, 1982. View at: Google Scholar
 M. Fujimoto, T. Hanawa, I. Iben Jr., and M. B. Richardson, “Thermal evolution of accreting neutron stars,” The Astrophysical Journal, vol. 278, pp. 813–824, 1984. View at: Publisher Site  Google Scholar
 T. Hanawa and M. Y. Fujimoto, “Thermal response of neutron stars to shell flashes,” Publication of the Astronomical Society of Japan, vol. 36, pp. 199–214, 1984. View at: Google Scholar
 C. Angulo, M. Arnould, and M. Rayet, “A compilation of chargedparticle induced thermonuclear reaction rates,” Nuclear Physics A, vol. 656, no. 1, pp. 3–183, 1999. View at: Publisher Site  Google Scholar
 R. K. Wallace and S. E. Woosley, “Explosive hydrogen burning,” Astrophysical Journal, vol. 45, p. 389, 1981. View at: Publisher Site  Google Scholar
 T. Hanawa, D. Sugimoto, and M. A. Hashimoto, “Nucleosynthesis in explosive hydrogen burning and its implications in tenminute interval of Xray bursts,” Publications of the Astronomical Society of Japan, vol. 35, pp. 491–506, 1983. View at: Google Scholar
 W. L. Slattery, G. D. Doolen, and H. E. Dewitt, “Improved equation of state for the classical onecomponent plasma,” Physical Review A, vol. 21, no. 6, pp. 2087–2095, 1980. View at: Publisher Site  Google Scholar
 M. B. Richardson, H. M. Van Horn, K. F. Ratcliff, and R. C. Malone, “Neutron star evolutionary sequences,” The Astrophysical Journal, vol. 255, pp. 624–653, 1982. View at: Publisher Site  Google Scholar
 V. Canuto, “Equation of state at ultrahigh sensities,” Annual Review of Astronomy and Astrophysics, vol. 12, pp. 167–214, 1974. View at: Publisher Site  Google Scholar
 V. Canuto, “Equation of state at ultrahigh densities,” Annual Review of Astronomy and Astrophysics, vol. 13, pp. 335–380, 1975. View at: Publisher Site  Google Scholar
 B. L. Friman and O. V. Maxwell, “Neutrino emissivities of neutron stars,” The Astrophysical Journal, vol. 232, pp. 541–557, 1979. View at: Google Scholar
 G. G. Festa and M. A. Ruderman, “Neutrinopair bremsstrahlung from a degenerate electron gas,” Physical Review, vol. 180, no. 5, pp. 1227–1231, 1969. View at: Publisher Site  Google Scholar
 G. Beaudet, V. Petrosian, and E. E. Salpeter, “Energy losses due to neutrino processes,” The Astrophysical Journal, vol. 150, p. 979, 1967. View at: Publisher Site  Google Scholar
 R. C. Malone, Cooling of superfluid neutron stars [Ph.D thesis], Cornell University, Ithaca, NY, USA, 1974.
 G. Baym, C. Pethick, and D. Pines, “Electrical conductivity of neutron star matter,” Nature, vol. 224, pp. 674–675, 1969. View at: Publisher Site  Google Scholar
 E. Flowers and N. Itho, “Transport properties of dense matter,” The Astrophysical Journal, vol. 206, pp. 218–242, 1976. View at: Google Scholar
 N. Itoh, N. Matsumoto, M. Seki, and Y. Kohyama, “Neutrinopair bremsstrahlung in dense stars. II—crystalline lattice case,” The Astrophysical Journal, vol. 279, pp. 413–418, 1984. View at: Publisher Site  Google Scholar
 I. Iben Jr., “Thermal pulses; pcapture, alphacapture, sprocess nucleosynthesis; and convective mixing in a star of intermediate mass,” The Astrophysical Journal, vol. 196, part 1, pp. 525–547, 1975. View at: Publisher Site  Google Scholar
 A. N. Cox and J. N. Stewart, “Rosseland opacity tables for population I compositions,” The Astrophysical Journal Supplement, vol. 19, p. 243, 1970. View at: Publisher Site  Google Scholar
 A. N. Cox and J. N. Stewart, “Rosseland opacity tables for population II compositions,” The Astrophysical Journal, vol. 19, Supplement, p. 261, 1970. View at: Google Scholar
 W. B. Hubbard and M. Lampe, “Thermal conduction by electrons in stellar matter,” The Astrophysical Journal Supplement Series, vol. 18, p. 297, 1969. View at: Publisher Site  Google Scholar
 V. Canuto, “Electrical conductivity and conductive opacity of a relativistic electron gas,” The Astrophysical Journal, vol. 159, pp. 641–652, 1970. View at: Publisher Site  Google Scholar
 H. E. Dewitt, H. C. Gravoske, and M. S. Cooper, “Screening factors for nuclear reactions. I. General theory,” The Astrophysical Journal, vol. 181, pp. 439–456, 1973. View at: Google Scholar
 P. Haensel and J. L. Zdunik, “Nonequilibrium processes in the crust of an accreting neutron star,” Astronomy and Astrophysics, vol. 227, no. 2, pp. 431–436, 1990. View at: Google Scholar
 D. K. Galloway, M. P. Muno, J. M. Hartman, D. Psaltis, and D. Chakrabarty, “Thermonuclear (type i) Xray bursts observed by the rossi Xray timing explorer,” Astrophysical Journal, vol. 179, no. 2, pp. 360–422, 2008. View at: Publisher Site  Google Scholar
 E. Kuulkers, J. J. M. in't Zand, and J.L. Atteia, “What ignites on the neutron star of 4U 0614+091?” Astronomy & Astrophysics, vol. 514, article A65, 2010. View at: Publisher Site  Google Scholar
 L. Bildsten, Thermonuclear Burning on Rapidly Accreting Neutron Stars, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1998.
 E. P. J. van den Heuvel, in Proceedings of the 2nd BeppoSAX Meeting, E. P. J. van den Heuvel, J. J. M. In't Zand, and R. A. M. J. Wijers, Eds., 2003.
 N. N. Weinberg and L. Bildsten, “Carbon detonation and shocktriggered helium burning in neutron star superbursts,” Astrophysical Journal Letters, vol. 670, no. 2, pp. 1291–1300, 2007. View at: Publisher Site  Google Scholar
 Y. Matsuo, H. Tsujimoto, and T. Noda, “Effects of a new tripleα reaction on xray bursts of a heliumaccreting neutron star,” Progress of Theoretical Physics, vol. 126, article 1177, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Masaaki Hashimoto et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.