[Mesa-users] Nucleosynthesis using MESA

Mon Apr 19 00:51:20 UTC 2021

Dear Dr. Timmes,

I ran the constant density and temperature trajectory until 1e12 s to compare the results with the plot that you shared.
The yields that I am getting are orders of magnitude different than the yields presented in the plot, for example, ni56 in the plot is of the order 10^-(1), however, I am getting ni56 of order 10^(-27).
Any suggestions on what I am doing wrong and the changes that I should make would be greatly appreciated.
I have copied my inlist and the final abundance file below.

Sincerely,
Sudarshan Neopane

My inlist_one_zone_burn file is as follows:
-------------------------------------------------------------------------------
&one_zone

net_name = 'approx21.net'

data_filename = 'const_temp_dens_one_zone_burn.data'
data_heading_line = 'temp=3e9 rho=1e9'

burn_temp = 3e9
burn_rho = 1e9
burn_tend = 1e12 ! seconds

show_net_reactions_info = .true.

min_for_show_peak_abundances = 1d-14 ! output peak value if it is > this

save_final_abundances = .true.
final_abundances_filename = 'const_temp_dens_final.data'

num_isos_for_Xinit = 2
names_of_isos_for_Xinit(1) = 'c12'
values_for_Xinit(1) = 5.0d-01
names_of_isos_for_Xinit(2) = 'o16'
values_for_Xinit(2) = 5.0d-01

num_names_of_isos_to_show = 15
names_of_isos_to_show(1) = 'he4'
names_of_isos_to_show(2) = 'c12'
names_of_isos_to_show(3) = 'o16'
names_of_isos_to_show(4) = 'ne20'
names_of_isos_to_show(5) = 'mg24'
names_of_isos_to_show(6) = 'si28'
names_of_isos_to_show(7) = 's32'
names_of_isos_to_show(8) = 'ar36'
names_of_isos_to_show(9) = 'ca40'
names_of_isos_to_show(10) = 'ti44'
names_of_isos_to_show(11) = 'cr48'
names_of_isos_to_show(12) = 'fe52'
names_of_isos_to_show(13) = 'ni56'
names_of_isos_to_show(14) = 'fe54'
names_of_isos_to_show(15) = 'cr56'

burn_rtol = 1d-5
burn_atol = 1d-6
which_solver = 'ros2_solver'

decsol_switch = 50
! if current number of species <= switch,
! then use small_mtx_decsol,
! else use large_mtx_decsol.
small_mtx_decsol = 'lapack'
large_mtx_decsol = 'klu'

screening_mode = 2 ! extended_screening
theta_e_for_graboske_et_al = 1

which_rates_choice = 2
! rates_NACRE_if_available = 1
! rates_JR05_if_available = 2 << for jina reaclib

cache_suffix = '' ! use this only for reactions from rate_tables
/
---------------------------------------------------------------------------------

Also, the top four abundances from my run are:

fe56           1    9.9999051742102363D-01
fe54           2    9.4826266268704317D-06
neut           3    1.8172964571281977D-13
c12           4    1.4548950960774552D-14

-------------------------------------------------------------------------------

________________________________
From: Mesa-users <mesa-users-bounces at lists.mesastar.org> on behalf of Sudarshan Neopane via Mesa-users <mesa-users at lists.mesastar.org>
Sent: Monday, April 19, 2021 4:48 AM
To: Francis Timmes <fxt44 at mac.com>
Cc: mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] Nucleosynthesis using MESA

[EXTERNAL SENDER]
Dear Dr. Timmes,

Thanks for the prompt response and suggestion.
I will run the constant density and temperature trajectory for verification.

Sincerely,
Sudarshan Neopane

________________________________
From: Francis Timmes <fxt44 at mac.com>
Sent: Monday, April 19, 2021 4:16 AM
To: Sudarshan Neopane <sneopane at umassd.edu>
Cc: Frank Timmes <fxt44 at mac.com>; mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] Nucleosynthesis using MESA

[EXTERNAL SENDER]

hi sudarshan,

i suggest trying something simple for a verification run,
such as a trajectory that has a constant temperature (say 3e9 k),
constant density (say 1e9 g/cc), and starts from 50-50 carbon and oxygen.
compare the results to https://nam12.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcococubed.asu.edu%2Fimages%2Fnet_aprox13%2Faprox13_hydrostatic.svg&data=04%7C01%7Csneopane%40umassd.edu%7Cbbb10836befb4f8cc46a08d902b9c33c%7C328d6c0d0f2f4b7693109762ba1c3e2d%7C0%7C0%7C637543819132699970%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=PH4q1q3NHPa%2F%2FYt%2FhgKJ8VnVmQYLLp1B8NAicJ3prY4%3D&reserved=0<https://nam12.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcococubed.asu.edu%2Fimages%2Fnet_aprox13%2Faprox13_hydrostatic.svg&data=04%7C01%7Csneopane%40umassd.edu%7C4c6cbb276ffa47fc6a2708d902be4161%7C328d6c0d0f2f4b7693109762ba1c3e2d%7C0%7C1%7C637543838435409574%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=jLXniq4JnHbV8P2MFf7vQW4I%2BohnxS4yXlNAeNeAphg%3D&reserved=0>

fxt

> On Apr 18, 2021, at 3:19 PM, Sudarshan Neopane via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> Dear Mesa Users,
>
> I have been trying to run particle trajectories in MESA for the nucleosynthesis of type Ia supernovae events. To run the particle trajectories, I followed the tutorial on running trajectories from 2012 MESA summer school. I could get the trajectory running, but wanted to know if I am running everything properly or not.
>
> So, does anyone have a sample trajectory and inlist file, along with the nucleosynthetic result of the trajectory that they can share? The sample trajectory doesn't need to be for a supernova event.
>
> Sincerely,
> Sudarshan Neopane
> University of Massachusetts Dartmouth
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://nam12.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=04%7C01%7Csneopane%40umassd.edu%7Cbbb10836befb4f8cc46a08d902b9c33c%7C328d6c0d0f2f4b7693109762ba1c3e2d%7C0%7C0%7C637543819132699970%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=fz%2BO%2FAoyi2i7fzczQPq%2BsGJjep5EIHESFNU4qITN%2F4w%3D&reserved=0<https://nam12.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=04%7C01%7Csneopane%40umassd.edu%7C4c6cbb276ffa47fc6a2708d902be4161%7C328d6c0d0f2f4b7693109762ba1c3e2d%7C0%7C1%7C637543838435419567%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=3q%2BsDoCVAgB1v544hhNgIM3H3%2BwHeqA3ynZUAqamexI%3D&reserved=0>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210419/844448b2/attachment.htm>