Computational method for a fractional model of the helium burning network

Stellar cores may be considered as a nuclear reactor that play important role in injecting new synthesized elements in the interstellar medium. Helium burning is an important stage that contribute to the synthesis of key elements such as carbon, through the triple- α process, and oxygen. In the present paper, we introduce a computational method for the fractional model of the nuclear helium burning in stellar cores. The system of fractional differential equations is solved simultaneously using series expansion method. The calculations are performed in the sense of modified Riemann-Liouville fractional derivative. Analytic expressions are obtained for the abundance of each element as a function of time. Comparing the abundances calculated at the fractional parameter α=1, which represents the integer solution, with the numerical solution revealed a good agreement with maximum error ɛ=0.003. The product abundances are calculated at α=0.25,0.5,0.75 to declare the effects of changing the fractional parameters on the calculations.