[Mesa-users] Evaluating rates been truncated when the lgT is above 10.000

[Bill Paxton]

> On Dec 26, 2020, at 5:44 PM, Francis Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> hi long,
> 
> published fits to tabulated thermonuclear reaction data are usually done
> temperatures less than 1e10k. evaluating those fits beyond this upper bound
> can produce silly numbers, hence the clipping and warning. as material at
> these temperatures is generally in nuclear statistical equilibrium, 
> such clipping is unlikely to significantly impact the evolution.
> 
> fxt

hi long,

Frank has it right of course.  One implementation detail to consider.  For such high T’s, you should be doing “op_split_burn” in mesa/star so that it can deal with the integration challenge.   For the test_suite cases that do core collapse, you’ll find things like this in the inlists:

      op_split_burn = .true.
      op_split_burn_min_T = 1d9
      burn_steps_limit = 150
      burn_steps_hard_limit = 250
      op_split_burn_eps = 1d-5
      op_split_burn_odescal = 1d-5

cheers,
bill



> 
> 
> 
>> On Dec 24, 2020, at 3:14 PM, 姜龙 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>> 
>> Dear all,
>> 
>>        When running to FeCCSN of Helium star, I encounter a list of warnings like:
>> 
>>        "evaluating rates with lgT=10.042050 which is above lgT=10.0000, rates have been truncated"
>> 
>>        Question: Is the truncation affect the results of evolution?  If Yes, how to avoid it or anything to do?    
>> 
>>        Best wishes.
>> 
>>        Long
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>> 
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>