[Mesa-users] the problem of making light curve

[Jared Goldberg]

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/20210605/74ac365b/attachment.htm>