[Mesa-users] nucleosynthesis of PISN

[Rob Farmer]

>And it should be after the explosive burning, right?

What happens if you let the model run for longer? Do you get a different
nucleosynthetic signature?

Rob

On Fri, 14 Jul 2023, 10:29 辛文宇, <xinwenyu16 at mails.ucas.ac.cn> wrote:

> Hi, Rob,
>
>
> > Sounds like an interesting science experiment for you to do to see what
> physics you need to change in a MESA model
>
> > to reproduce the results of Takahashi et al 2018.
>
> Yes, I am trying to reproduce the results of Takahashi et al 2018.
>
> Here is the output in the terminal at the final.
>
> =====================================================================================================================================================
>
>
> 4071   7.825243  6.570E+04   9.038694  11.441619  1.020E+02  0.000E+00
>  0.000000   0.000000   1.000000  -6.903022   2574     18
> -7.5315E+00   7.825237   2.141591 -43.681218  11.062596 -99.000000
> 102.000000   0.000000   0.000007  -0.99E+02   0.568267      3
>  2.6365E+05   0.954152   8.508267  11.439898 -21.444051 -11.742023
> 95.570074   0.000002   0.000000   0.504892  0.132E+03    lgRho_cntr
>                        log_rel_run_E_err       -6.6458850423420524
>                                rel_E_err    5.5001182058632353D-14
>
>
>  k, q, energy_removed_layers before adjustment is        2574
>  7.8177274553515981E-008   4.0017822825775093E+052
>  k, q, energy_removed_layers after adjustment is        2574
>  7.8177274553515981E-008   4.0017822825775093E+052
>  Layers above q=   7.8177274553515981E-008 will be removed
>  checking for conditions inside q=  0.98999999999999999      of material
> that will remain
>  check time left   1010.0351562500000
>  max vel inside fraction of cutoff  -1.7976931348623158E+303
>  max c/cs inside fraction of cutoff  -1.7976931348623157E+308
>  average v/vesc outside cutoff   2.4437965679800371
>  Kinetic plus potential energy outside cutoff   4.0017822825775093E+052
>  mass inside cutoff        2574   7.9740820044586302E-006
>  relax counter           0
>  log max v [km/s]=   5.3905238713708199
>  save LOGS/pisn_prof.data for model 4071
>  Entire star is expanding above the escape velocity, PISN!
>
>
>                  runtime (minutes), retries, steps      326.92        18
>     4071
>
>
> =====================================================================================================================================================
>
>
> It tells me: Entire star is expanding above the escape velocity, PISN!
>
> The entire star is unbound means it has exploded as a PISN.
>
> And it should be after the explosive burning, right?
>
> But why the production of 56Ni and 56Fe are so different to Takahashi's.
>
>
> BTW, is there someone who has calculated the nucleosynthesis of PISN with
> MESA?
>
>
> Thanks,
>
> Wenyu
>
>
>
> -----原始邮件-----
> *发件人:*"Rob Farmer" <robert.j.farmer37 at gmail.com>
> *发送时间:*2023-07-14 16:04:39 (星期五)
> *收件人:* "辛文宇" <xinwenyu16 at mails.ucas.ac.cn>
> *抄送:* mesa-users at lists.mesastar.org
> *主题:* Re: [Mesa-users] nucleosynthesis of PISN
>
> Takahashi et al 2018 uses a different code to MESA, so it's not surprising
> there are differences.
>
> Sounds like an interesting science experiment for you to do to see what
> physics you need to change in a MESA model to reproduce the results of
> Takahashi et al 2018
>
> Look at the terminal output when the model ends. That will tell you when
> and why the model stopped.
>
> Rob
>
> On Thu, 13 Jul 2023 at 16:28, 辛文宇 <xinwenyu16 at mails.ucas.ac.cn> wrote:
>
>> Hi, everyone,
>>
>>
>> I am trying to calculate the nucleosynthesis yields of PISN with the
>> "ppisn" package of mesa-r23.05.1.
>>
>> the initial He star mass is 102 Msun, Z=0 and I use the network of
>> mesa_128.net.
>>
>> More detail please check the inlist0.zip.
>>
>> I also attached the T-rho diagram and the abundance pattern of this test.
>>
>>
>> The problem is:
>>
>>
>> Compared with the model in Takahashi 2018 ApJ.875:111,
>> We have a higher [Si, S, Ar, Ca/Fe], about 1 dex, but lower
>> [Cl, K, Sc/Fe] also about 1 dex. Specially, we have extremely
>> higher Co/Fe and Ni/Fe in this model.
>> I also note that the total mass of 56Fe and Ni56 are 7E-8 and 3.3 Msun
>> in this model.
>>
>> It seems that this abundance pattern is before explosion and not that
>>
>>
>> after PISN explosion. Right? If right, how can we get the nucleosynthesis
>> yields
>>
>> after explosion?
>>
>>
>> Thanks,
>>
>> Wenyu
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230714/04940471/attachment.htm>