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

Robert Farmer rjfarmer at asu.edu
Wed Apr 20 13:01:53 EDT 2016


Hi
So i can't tell how well your data matches as that is missing from the the
email. The luminosity is expected (and did change for  MESA3, but that was
left out of the final paper)  to change as you change the network size as
larger nets allow you to complete more of the rp process and thus release
more energy. As flux is a distance dependent quantity we simply rescaled
each model to have the same peak flux, as the rise and decay time are the
more interesting quantities than the peak flux.

Rob

On Wed, Apr 20, 2016 at 4:02 AM, 박병찬 (물리 및 응용수학과) <allmest at unist.ac.kr>
wrote:

> 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.
>
>
> 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.
>
>
> 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
>
>
> ------------------------------------------------------------------------------
> Find and fix application performance issues faster with Applications
> Manager
> Applications Manager provides deep performance insights into multiple
> tiers of
> your business applications. It resolves application problems quickly and
> reduces your MTTR. Get your free trial!
> https://ad.doubleclick.net/ddm/clk/302982198;130105516;z
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/caddc409/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 46345 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/caddc409/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pastedImage.png
Type: image/png
Size: 20210 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160420/caddc409/attachment-0001.png>


More information about the Mesa-users mailing list