[Mesa-users] Pitfalls of eosPTEH?

[Bill Paxton]

> On Dec 28, 2018, at 12:48 PM, Zhang, Zhaoxi M. (Michael) <mz at astro.caltech.edu> wrote:
> 
> Hi Bill,
> 
> I am reusing the old models.  They still work in version 10204 by setting use_eosPTEH_with_eosDT = .false. and use_PTEH_for_high_Z = .false. , but since version 10382, use_eosPTEH_with_eosDT is gone as an option and the old models no longer work.  Do you know how to restore the pre-10382 behavior?

         !### use_eosPTEH_for_low_density
         !### use_eosPTEH_for_high_Z
         !### Z_for_all_PTEH
         !### Z_for_any_PTEH
         !### logRho_min_for_all_OPAL
         !### logRho_min_for_any_OPAL
         
         ! overwrite the control of same name in the `Eos_General_Info` structure
         
         ! if use_eosPTEH_for_low_density is true, 
         ! for logT < 3.7, use PTEH for logRho < about 7.5
         ! for higher logT, use PTEH for logRho < logRho_min_for_any_OPAL
         ! and do blend of OPAL and PTEH for logRho < logRho_min_for_any_OPAL.
         
         ! if use_eosPTEH_for_high_Z is true.
         ! use PTEH in place of HELM for Z above OPAL/SCVH range
         ! Z >= Z_for_all_PTEH, then use pure PTEH with no blend
         ! Z < Z_for_any_PTEH, then do not use PTEH
         ! blend with PTEH for intermediate Z.
         ! if use_eosDT2, blend with DT2.
         ! else blend with OPAL/SCVH data.
         
      use_eosPTEH_for_low_density = .true.
      use_eosPTEH_for_high_Z = .true.
      Z_for_all_PTEH = 0.040d0
      Z_for_any_PTEH = 0.039d0
      logRho_min_for_all_OPAL = -9.1d0
      logRho_min_for_any_OPAL = -9.2d0 ! -9.6d0 

> 
> I can also try to generate new models, but the white dwarf models in star_data/white_dwarf_models have been deleted in the newest version, so I don't have an easy starting point.  Do you plan to regenerate them?

no.


> 
> 
> Thanks for your help,
> 
> Michael
> From: Bill Paxton <paxton at kitp.ucsb.edu <mailto:paxton at kitp.ucsb.edu>>
> Sent: Friday, December 28, 2018 7:35:49 AM
> To: Zhang, Zhaoxi M. (Michael)
> Cc: mesa-developers at lists.mesastar.org <mailto:mesa-developers at lists.mesastar.org>; mesa-users at lists.mesastar.org <mailto:mesa-users at lists.mesastar.org>
> Subject: Re: [Mesa-users] Pitfalls of eosPTEH?
>  
> Are you reusing old models or building new ones?  
> 
> For the cases you mention, the differences between the older eos and PTEH could be large, so it might not be possible to load a model created with the one and expect it to work with the other.
> 
> -b
> 
> 
> 
> 
>> On Dec 28, 2018, at 1:16 AM, Zhang, Zhaoxi M. (Michael) <mz at astro.caltech.edu <mailto:mz at astro.caltech.edu>> wrote:
>> 
>> Hello MESA users & developers,
>> 
>> I'm trying to reproduce the simulations I did for my paper (https://arxiv.org/abs/1812.08793 <https://arxiv.org/abs/1812.08793>) with the latest MESA code, but they fail immediately with a lot of retries and backups.  I managed to determine that version number 10,204--specifically, the introduction of eosPTEH--is what broke my simulations.  Turning off PTEH by setting use_eosPTEH_with_eosDT = .false. and use_PTEH_for_high_Z = .false. fixes the problem.
>> 
>> Is there a reason why PTEH isn't good to use in my scenario, namely iron enriched C/O white dwarfs with hot inflated envelopes?  Reading the latest release notes, it seems like low density and high metallicity environments are precisely what PTEH is meant for, so I have no idea why MESA can't even find a first model.  
>> 
>> Thanks for any help,
>> 
>> 
>> Michael
>> _______________________________________________
>> mesa-users at lists.mesastar.org <mailto:mesa-users at lists.mesastar.org>
>> https://lists.mesastar.org/mailman/listinfo/mesa-users <https://lists.mesastar.org/mailman/listinfo/mesa-users>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181228/2a80be87/attachment.html>