.. _jdot_ml_check: ************* jdot_ml_check ************* .. tags:: binary, angular-momentum, orbital-angular-momentum, mass-loss, mass-transfer, roche-lobe-overflow, analytic-check, point-mass, low-mass This test case checks the evolution of one star plus a point mass, including mass transfer to the point mass. Using pre-specified efficiency options, verifies that the binary evolution, including mass transfer to the point mass, follows the analytical result of `Tauris & van den Heuvel (2006) `_. This test case has 1 part. Click to see a larger version of a plot. * Part 1 (``inlist_project``) loads a 1 Msun zero-age main sequence model and a 1.4 Msun point mass. The binary system is given an initial separation of 2.75 Rsun, and is evolved for 14 billion years. Along the way, the ``run_binary_extras.f90`` calculates the deviation from the analytical result of `Tauris & van den Heuvel (2006) `_: .. code-block:: console error in separation 7.0519985646485173D-05 Binary period and separation history: .. image:: ../../../binary/test_suite/jdot_ml_check/docs/grid_000891.svg :width: 100% Abundance and power profiles of star 1 at termination: .. image:: ../../../binary/test_suite/jdot_ml_check/docs/profile_s1_000891.svg :width: 100% pgstar commands used for the binary and star1 plots above: .. code-block:: console &pgstar file_white_on_black_flag = .true. ! white_on_black flags -- true means white foreground color on black background !file_device = 'png' ! png file_device = 'vcps' ! postscript pgstar_interval = 1 pgstar_left_yaxis_label_disp = 3.0 Text_Summary1_win_flag = .false. Text_Summary1_win_width = 10 Text_Summary1_win_aspect_ratio = 0.15 Text_Summary1_xleft = 0.01 Text_Summary1_xright = 0.99 Text_Summary1_ybot = 0.0 Text_Summary1_ytop = 1.0 Text_Summary1_txt_scale = 1.0 Text_Summary1_title = '' Text_Summary1_num_rows = 4 ! <= 20 Text_Summary1_num_cols = 8 ! <= 20 Text_Summary1_name(:,:) = '' Text_Summary1_name(1,1) = 'time_step' Text_Summary1_name(1,2) = 'donor_index' Text_Summary1_name(1,3) = 'binary_separation' Text_Summary1_name(1,4) = 'period_days' Text_Summary1_name(1,5) = 'eccentricity' Text_Summary1_name(1,6) = 'edot' Text_Summary1_name(1,7) = 'J_orb' Text_Summary1_name(1,8) = '' Text_Summary1_name(2,1) = 'sum_of_masses' Text_Summary1_name(2,2) = 'star_1_mass' Text_Summary1_name(2,3) = 'star_1_radius' Text_Summary1_name(2,4) = 'v_orb_1' Text_Summary1_name(2,5) = 'rl_1' Text_Summary1_name(2,6) = 'lg_mstar_dot_1' Text_Summary1_name(2,7) = 'J_spin_1' Text_Summary1_name(2,8) = '' Text_Summary1_name(3,1) = '' Text_Summary1_name(3,2) = 'star_2_mass' Text_Summary1_name(3,3) = 'star_2_radius' Text_Summary1_name(3,4) = 'v_orb_2' Text_Summary1_name(3,5) = 'rl_2' Text_Summary1_name(3,6) = 'lg_mstar_dot_2' Text_Summary1_name(3,7) = 'J_spin_2' Text_Summary1_name(3,8) = '' Text_Summary1_name(4,1) = 'Jdot' Text_Summary1_name(4,2) = 'jdot_gr' Text_Summary1_name(4,3) = 'jdot_mb' Text_Summary1_name(4,4) = 'jdot_ls' Text_Summary1_name(4,5) = 'jdot_ml' Text_Summary1_name(4,6) = '' Text_Summary1_name(4,7) = '' Text_Summary1_name(4,8) = '' History_Track1_win_flag = .false. History_Track1_win_width = 12 History_Track1_win_aspect_ratio = 0.75 History_Track1_txt_scale = 0.8 History_Track1_title = ' ' History_Track1_xname = 'model_number' History_Track1_yname = 'period_days' History_Track1_xaxis_label = 'Model Number' History_Track1_yaxis_label = 'Binary Period (day)' History_Track1_reverse_xaxis = .false. History_Track1_reverse_yaxis = .false. History_Track1_xmin = 0.0 History_Track1_xmax = 900.0 History_Track1_ymin = 0.3 History_Track1_ymax = 0.36 History_Track1_file_flag = .false. History_Track1_file_dir = 'pgstar_out' History_Track1_file_prefix = 'track1_s1_' History_Track1_file_interval = 10000 History_Track1_file_width = -1 History_Track1_file_aspect_ratio = -1 History_Track2_win_flag = .false. History_Track2_win_width = 12 History_Track2_win_aspect_ratio = 0.75 History_Track2_txt_scale = 0.8 History_Track2_title = '' History_Track2_xname = 'model_number' History_Track2_yname = 'binary_separation' History_Track2_xaxis_label = 'Model Number' History_Track2_yaxis_label = 'Binary Separation (R\d\(2281)\u)' History_Track2_reverse_xaxis = .false. History_Track2_reverse_yaxis = .false. History_Track2_xmin = 0.0 History_Track2_xmax = 900 History_Track2_ymin = 2.0 History_Track2_ymax = 3.0 History_Track2_file_flag = .false. History_Track2_file_dir = 'pgstar_out' History_Track2_file_prefix = 'track2_s1_' History_Track2_file_interval = 10000 History_Track2_file_width = -1 History_Track2_file_aspect_ratio = -1 Grid1_win_flag = .true. Grid1_win_width = 16 Grid1_win_aspect_ratio = 0.6 Grid1_title = 'jdot_ml_check' Grid1_txt_scale_factor(1) = 1.0 Grid1_num_cols = 2 ! divide plotting region into this many equal width cols Grid1_num_rows = 2 ! divide plotting region into this many equal height rows Grid1_num_plots = 3 ! <= 10 Grid1_plot_name(1) = 'Text_Summary1' Grid1_plot_row(1) = 1 ! number from 1 at top Grid1_plot_rowspan(1) = 1 ! plot spans this number of rows Grid1_plot_col(1) = 1 ! number from 1 at left Grid1_plot_colspan(1) = 2 ! plot spans this number of columns Grid1_plot_pad_left(1) = -0.10 ! fraction of full window width for padding on left Grid1_plot_pad_right(1) = -0.18 ! fraction of full window width for padding on right Grid1_plot_pad_top(1) = -0.00 ! fraction of full window height for padding at top Grid1_plot_pad_bot(1) = 0.32 ! fraction of full window height for padding at bottom Grid1_txt_scale_factor(1) = 0.70 ! multiply txt_scale for subplot by this Grid1_plot_name(2) = 'History_Track1' Grid1_plot_row(2) = 2 ! number from 1 at top Grid1_plot_rowspan(2) = 1 ! plot spans this number of rows Grid1_plot_col(2) = 1 ! number from 1 at left Grid1_plot_colspan(2) = 1 ! plot spans this number of columns Grid1_plot_pad_left(2) = -0.04 ! fraction of full window width for padding on left Grid1_plot_pad_right(2) = 0.04 ! fraction of full window width for padding on right Grid1_plot_pad_top(2) = -0.24 ! fraction of full window height for padding at top Grid1_plot_pad_bot(2) = 0.00 ! fraction of full window height for padding at bottom Grid1_txt_scale_factor(2) = 0.7 ! multiply txt_scale for subplot by this Grid1_plot_name(3) = 'History_Track2' Grid1_plot_row(3) = 2 ! number from 1 at top Grid1_plot_rowspan(3) = 1 ! plot spans this number of rows Grid1_plot_col(3) = 2 ! number from 1 at left Grid1_plot_colspan(3) = 1 ! plot spans this number of columns Grid1_plot_pad_left(3) = 0.04 ! fraction of full window width for padding on left Grid1_plot_pad_right(3) = -0.04 ! fraction of full window width for padding on right Grid1_plot_pad_top(3) = -0.24 ! fraction of full window height for padding at top Grid1_plot_pad_bot(3) = 0.00 ! fraction of full window height for padding at bottom Grid1_txt_scale_factor(3) = 0.7 ! multiply txt_scale for subplot by this Grid1_file_flag = .true. Grid1_file_dir = 'pgstar_out' Grid1_file_prefix = 'grid_' Grid1_file_interval = 10000 Grid1_file_width = -1 Grid1_file_aspect_ratio = -1 Profile_Panels2_win_flag = .true. Profile_Panels2_win_width = 10 Profile_Panels2_title = 'jdot_ml_check: star 1' Profile_Panels2_xaxis_name = 'mass' Profile_Panels2_xaxis_reversed = .false. Profile_Panels2_xmin = 0.0 Profile_Panels2_xmax = 1.0 Profile_Panels2_show_mix_regions_on_xaxis = .false. Profile_Panels2_xright = 0.92 Profile_Panels2_ytop = 0.92 num_abundance_line_labels = 5 Abundance_legend_max_cnt = 0 Profile_Panels2_yaxis_name(2) = 'Power' Profile_Panels2_ymin(2) = -10.0 Profile_Panels2_ymax(2) = 5.0 Profile_Panels2_file_flag = .true. Profile_Panels2_file_dir = 'pgstar_out' Profile_Panels2_file_prefix = 'profile_s1_' Profile_Panels2_file_interval = 10000 Profile_Panels2_file_width = -1 Profile_Panels2_file_aspect_ratio = -1 / ! end of pgstar namelist Last-Updated: 12Jul2021 (MESA 094ff71) by fxt.