[Mesa-users] the problem of making light curve

[游凱安]

Dear Jared,



I have tried make 13msun and 40msun model in make_pre_ccsn for three weeks,



but it still failed.



I have sent the problem of making 13Msun and 40Msun supernovae to you.



Can you tell me your current progress and status?



I think that the 40Msun massive star Finally becomes a black hole



And maybe the 13Msn becomes a white dwarf.



If  you need more information, please let me know.



Thanks for your patience and help.



~Kai-An


游凱安 <timyu930218 at gmail.com> 於 2021年6月6日 週日 下午1:18寫道：

> Dear Jared,
>
> Thanks for your helping.
>
> however, I have tried make 13msun and 40msun model in make_pre_ccsn for
> two weeks,
>
> but it still failed.
>
> I have sent the problem of making 13Msun and 40Msun supernovae to you.
>
> Can you tell me your situation or process?
>
> If  you need more information, please let me know.
>
> ~Kai-An
>
> Jared Goldberg <goldberg at physics.ucsb.edu> 於 2021年6月6日 週日 下午12:21寫道：
>
>> Hi Kai-an,
>>
>> It looks like my suspicion was correct, and your rn script for ccsn_IIp
>> is skipping steps and loading standard models. Add MESA_RUN_OPTIONAL=t to
>> your rn script or comment out the if statements and you should be good to
>> go.
>>
>> ~Jared
>>
>> On Sat, Jun 5, 2021 at 6:07 PM 游凱安 <timyu930218 at gmail.com> wrote:
>>
>>> Dear Jared,
>>>
>>> Below are inlist_mass_Z_wind_rotation and rn  in make_pre_ccsn,
>>>
>>> inlist_edep and rn in ccsn_llp.
>>>
>>> If you need more information, please let me know.
>>>
>>> ~Kai-An
>>>
>>>
>>> Jared Goldberg <goldberg at ucsb.edu> 於 2021年6月6日 週日 上午12:46寫道：
>>>
>>>> Hi Kai-an,
>>>>
>>>> For some reason, Stella is almost certainly reading the same input file
>>>> for all of these models. It’s particularly apparent because if there were
>>>> different Nickels it would necessarily affect the late time behavior.
>>>> However, we cant know exactly why without more information. Can you send
>>>> your inlists and your rn script?
>>>>
>>>> Without looking at those files, my guess is this:
>>>> Look at the difference between the old rn script and the new rn script.
>>>> The newest MESA revisions skip steps in the explosion and instead load
>>>> in a “standard” model. You might just be seeing the explosion of that
>>>> standard model.
>>>>
>>>> To fix this add
>>>> export MESA_RUN_OPTIONAL=False
>>>> To the rn script (or to your .bashrc file or both)
>>>>
>>>> Or Comment out the lines containing the if statements.
>>>>
>>>> And try the run again.
>>>>
>>>> One other thing I do with 15140 to make sure this doesn’t happen when I
>>>> copy a test_suite case is move the standard_*.mod files to a new directory,
>>>> just so MESA will crash if it tries to read one of them in.
>>>>
>>>> Also, this is something you may be able to debug on your own before
>>>> coming to the users list, if my guess is wrong. Compare the output at every
>>>> step of the simulation, see where the output starts to be the same, and
>>>> that will tell you which of the many steps things are failing.
>>>>
>>>> Hope this helps!
>>>>
>>>> Cheers,
>>>> ~Jared
>>>>
>>>>
>>>>
>>>> On Jun 5, 2021, at 2:22 AM, 游凱安 via Mesa-users <
>>>> mesa-users at lists.mesastar.org> wrote:
>>>>
>>>> 
>>>> Dear professors,
>>>>
>>>> I have issues with making light curve.
>>>>
>>>> I have several independent variable that control the development of
>>>> supernova,
>>>>
>>>> including *initial mass*
>>>> (13msun,15msun,20msun,25msun,30msun,35msun,40msun),
>>>>
>>>> *metallicity* (z=0.02,z=0.002,z=0.0002, the factor is new_Z in
>>>> inlist_mass_Z_wind_rotation of make_pre_ccsn),
>>>>
>>>> *explosion energy*(1.2foe(ergx1), 2.4foe(ergx2), 3.6foe(ergx3),
>>>> 4.8foe(ergx4), the factor is inject_until_reach_model_with_total_energy in
>>>> inlist_edep of ccsn_llp),
>>>>
>>>> *nickel mass*(the factor is x_ctrl(12) in inlist_shock_part5 of
>>>> ccsn_llp).
>>>>
>>>> as you can see, there are many models in it.
>>>>
>>>> BUT THE RESULTS ARE VERY STRANGE.
>>>>
>>>> below are the sample
>>>> <15:0.02:1:0.01.png>
>>>>
>>>> <15:0.002:1:0.01.png>
>>>>
>>>> <15:0.0002:1:0.01.png>
>>>>
>>>> <15:0.02:2:0.01.png>
>>>>
>>>> <15:0.002:2:0.01.png>
>>>>
>>>> <15:0.0002:2:0.01.png>
>>>>
>>>> <15:0.02:3:0.01.png>
>>>>
>>>> <15:0.0002:3:0.01.png>
>>>>
>>>> as you can see, *it is very strange that meltallicity and explosion
>>>> energy have no effect on the light curve*, but it make no sense.
>>>>
>>>> I was thinking that maybe I put the wrong model to plot or something,
>>>>
>>>> but every different masses model has such situation,
>>>>
>>>> the mesa version is r15140,
>>>>
>>>> *and I have used r-12778 before, I follow the steps mentioned and the
>>>> light curves made from it could change with the manipulation of nickel mass
>>>> and explosion energy*.(p.s. I have not tried metallicity before)
>>>>
>>>> how can I fix the problem?
>>>>
>>>> Below are my 15msun models(due to space storage, I only can transfer
>>>> one mass model to you)
>>>>
>>>>
>>>> https://drive.google.com/drive/folders/1247h2HCGoWgxvKY0EWTBHYW5NtziOrIn?usp=sharing
>>>>
>>>> If you need more information, please let me know.
>>>>
>>>> Thanks for your patience and help.
>>>>
>>>> ~Kai-An
>>>>
>>>> _______________________________________________
>>>> 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/20210611/6d45ca80/attachment.htm>