[mesa-users] Reproduce xray_burst results on MESA III paper.

Wed Apr 20 07:02:39 EDT 2016

Dear mesa-users

I'm trying to reproduce results of mesa for x-ray burst published on section 5 of 3rd MESA paper (MODULES FOR EXPERIMENTS IN STELLAR ASTROPHYSICS (MESA): BINARIES, PULSATIONS, AND EXPLOSIONS).

I used the inlist_xray of 'inlist_mesa3.zip' uploaded on 'Inlists of Published Results' in mesa community. The inlist says that it is for 'MESA r7484', but I used 'MESA r7499' because I couldn't find 'MESA r7484' on mesa sourceforge webpage.

[cid:cc67c97a-cbe4-4c1f-88eb-e9c67c75b156]

The inlist support to reproduce result for 4-type of network (see above table 5). I obtained logL vs star_age data from each simulation results and I draw "flux vs t`" data.

[cid:3f488b7f-8e40-4a56-a707-e291de323b0f]

The result of rp_53 is well matched up with corresponding result of 3rd MESA paper (see above figure 22), but result of rp_153 is obviously not matched up with the correseponding results of 3rd MESA paper. The peak value of rp_153 of 3rd MESA paper is about 28, but my result shows its peak value over 100!

I can't understand why I obtain different result. Here I attach text of my 'inlist_xray' for rp_153 simulation:

! Type 1 Xray bursts. MESA r7484
!
! How to run:
!
! Accretion rate:
! Set mass_change
!
!Accretion compostion:
! Set either accrete_given_mass_fractions = .true. and uncomment accretion_species_id section for the 2%N case
! or
! Set accrete_given_mass_fractions = .false. and uncomment accretion_h1 section for the 2% metals
!
! Choosing a net:
! RP nets:
!     change_initial_net=.true.
!     new_net_name = 'rp_53.net'
!
! Adaptive nets
!     enable_adaptive_network=.true.
!     max_Z_for_add = 52
!     max_N_for_add=82
!
!     min_x_for_keep = 1d-8
!     min_x_for_n = 1d-6
!     min_x_for_add = 1d-6
!

&star_job
      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = 'ns_relax.mod'

         change_initial_net = .true.
!         new_net_name = 'rp_53.net'
         new_net_name = 'rp_153.net'
!         new_net_name = 'rp_305.net'

!         enable_adaptive_network=.true.
!         max_Z_for_add = 52
!         max_N_for_add=82
!
!         min_x_for_keep = 1d-8
!         min_x_for_n = 1d-6
!         min_x_for_add = 1d-6

      set_rate_3a = 'FL87' ! Fushiki and Lamb, Apj, 317, 368-388, 1987
      set_rate_1212 = 'G05' ! Gasques, et al. Phys Review C, 72, 025806, 2005

      kappa_file_prefix = 'gs98'

! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
! Used when taking the NS mod file from ns_h test_suite to convert into ns_relax.mod.
! Dont try to change the net at the same time as using the relax_* options
!     replace_element=.true.
!     chem_name1='fe56'
!     chem_name2='ba138'
!     replace_element_nzlo=-1
!     replace_element_nzhi=-1
!
!     new_R_center=1.12d6
!     relax_initial_R_center = .true.
!
!     new_L_center=1.6d34
!     relax_initial_L_center = .true.
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !

/ ! end of star_job namelist

&controls

        fe_core_infall_limit = 1d99
        max_abar_for_burning = 137

        initial_mass = 1.4
        initial_z = 0.02d0

        use_Type2_opacities = .true.
        Zbase = 0.02d0

        photostep = 50
        profile_interval = 5
        history_interval = 1
        terminal_cnt = 10
        write_header_frequency = 10

        min_timestep_limit = 1d-9
        max_timestep_factor=2d0
        gamma_center_limit = 1d99

        varcontrol_target = 1d-3
        mesh_delta_coeff = 1.0

        mesh_min_dlnR = 1d-9
        merge_if_dlnR_too_small = .true.

        which_atm_option = 'grey_and_kap'
        Pextra_factor = 2

        RGB_wind_scheme = ''
        AGB_wind_scheme = ''
        RGB_to_AGB_wind_switch = 1d-4
        Reimers_wind_eta = 0.5d0
        Blocker_wind_eta = 0.1d0

        use_GR_factors = .true.

        accrete_same_as_surface = .false.
        accretion_dump_missing_metals_into_heaviest = .false.

        max_years_for_timestep = 1d-4
        set_min_D_mix = .true.
        min_D_mix = 1d6

        min_T_for_acceleration_limited_conv_velocity = 0.0

        okay_to_reduce_gradT_excess = .true.
        super_eddington_wind_eta = 1

        mixing_length_alpha = 1.5
        MLT_option = 'Henyey'

        use_Ledoux_criterion = .false.

        alpha_semiconvection = 0.1
        thermohaline_coeff = 2
        thermohaline_option = 'Kippenhahn'

        dX_nuc_drop_limit=7.5d-2
        dX_nuc_drop_limit_at_high_T=5d-3
        dX_nuc_drop_min_X_limit=7.5d-2
        dX_nuc_drop_max_A_limit=52

        newton_iterations_limit = 10

        mesh_dlog_pp_dlogP_extra = 1
        mesh_dlog_cno_dlogP_extra = 1

        mesh_dlog_3alf_dlogP_extra = 1
        mesh_dlog_burn_c_dlogP_extra = 1
        mesh_dlog_burn_n_dlogP_extra = 1
        mesh_dlog_burn_o_dlogP_extra = 1

        mesh_dlog_burn_ne_dlogP_extra = 1
        mesh_dlog_burn_na_dlogP_extra = 1
        mesh_dlog_burn_mg_dlogP_extra = 1

        mesh_dlog_cc_dlogP_extra = 1
        mesh_dlog_co_dlogP_extra = 1
        mesh_dlog_oo_dlogP_extra = 1

        mesh_dlog_burn_si_dlogP_extra = 1
        mesh_dlog_burn_s_dlogP_extra = 1
        mesh_dlog_burn_ar_dlogP_extra = 1
        mesh_dlog_burn_ca_dlogP_extra = 1
        mesh_dlog_burn_ti_dlogP_extra = 1
        mesh_dlog_burn_cr_dlogP_extra = 1
        mesh_dlog_burn_fe_dlogP_extra = 1

        mesh_dlog_pnhe4_dlogP_extra = 1
        mesh_dlog_other_dlogP_extra = 1
        mesh_dlog_photo_dlogP_extra = 1

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 2% N case
!      accrete_given_mass_fractions = .true.
!      num_accretion_species = 3
!      accretion_species_id(1) = 'h1'
!      accretion_species_xa(1) = 0.7048
!      accretion_species_id(2) = 'he4'
!      accretion_species_xa(2) = 0.2752
!      accretion_species_id(3) = 'n14'
!      accretion_species_xa(3) = 0.02
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 2% metals case
     accretion_h1 = 0.7048
     accretion_h2 = 0
     accretion_he3 = 0
     accretion_he4 = 0.2752
     accretion_zfracs = 3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Accretion rate
     mass_change = 3d-9
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/ ! end of controls namelist

&pgstar

! Useful plots to show during the run
        History_Track1_win_flag = .true.
        History_Track1_win_width = 12
        History_Track1_win_aspect_ratio = 0.75 ! aspect_ratio = height/width

        ! set default
        History_Track1_xname = 'star_age_hr'
        History_Track1_yname = 'log_L'
        History_Track1_xaxis_label = 'Age [hr]'
        History_Track1_yaxis_label = 'log L/L\d\(2281)'

        ! file output
        History_Track1_file_flag = .false.
        History_Track1_file_dir = 'png'
        History_Track1_file_prefix = 'track1_'
        History_Track1_file_cnt = 10 ! output when mod(model_number,History_Track1_file_cnt)==0
        History_Track1_file_width = -1 ! (inches) negative means use same value as for window
        History_Track1_file_aspect_ratio = -1 ! negative means use same value as for window

        !# TRho Profile window
        ! current model in T-Rho plane

        TRho_Profile_win_flag = .true.

        TRho_Profile_win_width = 12
        TRho_Profile_win_aspect_ratio = 0.75 ! aspect_ratio = height/width

        TRho_switch_to_Column_Depth = .true.
        ! if true, replace logRho for xaxis by log column depth (g/cm^2)

        show_TRho_Profile_legend = .true.
        TRho_Profile_legend_coord = 0.05
        TRho_Profile_legend_fjust = 0.0
        TRho_Profile_legend_disp1 = -2.0
        TRho_Profile_legend_del_disp = -1.3

        show_TRho_Profile_text_info = .false.
        TRho_Profile_text_info_xfac = 0.77 ! controls x location
        TRho_Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text
        TRho_Profile_text_info_yfac = 0.6 ! controls y location of 1st line
        TRho_Profile_text_info_dyfac = -0.04 ! controls line spacing

        show_TRho_Profile_mass_locs = .false.
        show_TRho_accretion_mesh_borders = .true.
        ! If true, plots the borders of the mesh and newly accreted material in the appropriate xcoord
        show_TRho_Profile_kap_regions = .false.
        show_TRho_Profile_eos_regions = .true.
        show_TRho_Profile_degeneracy_line = .true.
        show_TRho_Profile_Pgas_Prad_line = .true.
        show_TRho_Profile_burn_lines = .true.
        show_TRho_Profile_burn_labels = .true.
        show_TRho_Profile_logQ_limit = .true.

        show_TRho_Profile_annotation1 = .false.
        show_TRho_Profile_annotation2 = .false.
        show_TRho_Profile_annotation3 = .false.

        TRho_Profile_fname = '' ! file name for extra TRho profile data

        ! axis limits
        TRho_Profile_xmin = -5.0
        TRho_Profile_xmax = 10.2
        TRho_Profile_ymin = 2.6
        TRho_Profile_ymax = 10.2

        ! file output
        TRho_Profile_file_flag = .false.
        TRho_Profile_file_dir = 'png'
        TRho_Profile_file_prefix = 'trho_profile_'
        TRho_Profile_file_cnt = 5 ! output when mod(model_number,TRho_Profile_file_cnt)==0
        TRho_Profile_file_width = -1 ! (inches) negative means use same value as for window
        TRho_Profile_file_aspect_ratio = -1 ! negative means use same value as for window

/ ! end of pgstar namelist

I just changed the '!' sign from original inlist_xray.

Thanks

Park
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/aa758767/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 20210 bytes
Desc: pastedImage.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/aa758767/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 46345 bytes
Desc: pastedImage.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/aa758767/attachment-0001.png>