[Mesa-users] Custom Reaction Rates in MESA

[Brian Lutz]

Here are the rate files for the three pp1 reactions, and then the rate_list



#Random rates for 1h(p,e+nu)2h

15


0.01  2.11E-10

0.02  3.12E-10

0.03  1.11E-10

0.04  9.33E-10

0.05  8.32E-10

0.06  3.92E-10

0.07  2.76E-10

0.08  3.65E-10

0.09  1.87E-10

0.10  9.76E-9

0.11  8.61E-9

0.12  7.90E-8

0.13  4.89E-8

0.14  5.77E-7

0.15  6.43E-7



#Random rates for 2h(p,g)3he

15


0.01  8.71E-10

0.02  8.32E-10

0.03  8.21E-10

0.04  1.23E-10

0.05  1.62E-10

0.06  1.82E-10

0.07  7.06E-10

0.08  9.05E-10

0.09  6.17E-10

0.10  8.26E-9

0.11  1.41E-9

0.12  1.95E-8

0.13  8.84E-8

0.14  5.70E-7

0.15  3.49E-7



#Random rates for 3he(3he,2p)4he

15


0.01  2.00E-10

0.02  1.37E-10

0.03  3.87E-10

0.04  3.98E-10

0.05  5.67E-10

0.06  8.06E-10

0.07  9.00E-10

0.08  1.65E-10

0.09  1.91E-10

0.10  5.76E-9

0.11  9.43E-9

0.12  9.18E-8

0.13  7.90E-8

0.14  3.22E-7

0.15  1.99E-7




! this is an example of a rates list file for use with mesa/rates


! the mesa/data/net_data/rates directory has sample rate files


! pairs of rate name and rate file


r_h1_h1_wk_h2                   'toy_pp.txt'

r_h2_pg_he3                     'toy_dp.txt'

r_he3_he3_to_h1_h1_he4          'toy_2he3.txt'

On Fri, Oct 21, 2022 at 1:39 AM Rob Farmer <robert.j.farmer37 at gmail.com>
wrote:

> Hi
>
> What does your toy_pp.txt and your rate_list.txt file look like?
>
> Rob
>
> On Fri, 21 Oct 2022 at 08:44, Brian Lutz via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> Hello,
>>
>> I am working on creating a 1 solar mass model with custom reaction rates
>> (currently, the reaction rates are completely random) for the pp1 chain and
>> have been encountering a problem.
>> Inside my work directory, called Custom_PP_Rates, I have a directory
>> called customRates which includes the three rate tables (toy_pp.txt,
>> toy_dp.txt, toy_2he3.txt) as well as rate_list.txt.  I then call
>> rates_tables_dir = 'customRates' inside the &star_job section in my
>> inlist_project.  When I do ./mk and then ./rn I get the following error
>> message:
>>
>>
>> add_reaction_for_handle failed in reaclib_parse_handle r_h1_h1_wk_h2
>> 'toy_pp.txt'
>>
>>  invalid rate file r_h1_h1_wk_h2 'toy_pp.txt'
>>
>>  rates_init failed in init_raw_rates_records
>>
>>  read_raw_rates_records failed in net_setup_tables
>>
>>  set_net failed in net_tables
>>
>>  failed in set_net
>>
>>  star_create_pre_ms_model ierr          -1
>>
>>  do_load1_star ierr          -1
>>
>>  before_evolve_loop ierr          -1
>>
>> DATE: 2022-10-21
>>
>> TIME: 00:00:57
>>
>> The &star_job section of my inlist project is:
>>
>> &star_job
>>
>>   ! see star/defaults/star_job.defaults
>>
>>
>>   ! begin with a pre-main sequence model
>>
>>     create_pre_main_sequence_model = .true.
>>
>>
>>   ! save a model at the end of the run
>>
>>     save_model_when_terminate = .false.
>>
>>     save_model_filename = '1M.mod'
>>
>>
>>   ! nuclear reactions
>>
>>     rate_tables_dir = 'customRates'
>>
>>
>>
>>     change_initial_net = .true.
>>
>>     new_net_name = 'pp_and_cno_extras.net'
>>
>>
>>
>>   ! pause before terminate
>>
>>     pause_before_terminate = .true.
>>
>>
>>
>>   ! display on-screen plots
>>
>>     pgstar_flag = .true.
>>
>>
>>
>> / ! end of star_job namelist
>>
>> I would appreciate some support on how to properly add custom reaction
>> rates, and use them in my simulations.
>>
>> Cheers,
>>
>> Brian Lutz
>>
>> Physics Major
>>
>> South Dakota School of Mines & Technology
>>
>> Phone: 303-350-8767
>>
>> Email: Brian.Lutz at mines.sdsmt.edu
>> _______________________________________________
>> 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/20221021/254c83a0/attachment.htm>