[Mesa-users] the problem of making light curve

[游凱安]

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/20210606/14208fd4/attachment.htm>