[Mesa-users] Nucleosynthesis using MESA

Mon Apr 19 04:15:54 UTC 2021

hi sudarshan,

there may not be anything wrong; you are using a different different network.
take a look at whole page http://cococubed.asu.edu/code_pages/burn_helium.shtml 
for your verification runs,

fxt

> On Apr 18, 2021, at 5:51 PM, Sudarshan Neopane <sneopane at umassd.edu> wrote:
> 
> 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
> 
> 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