Duplicated Modes

Sometimes two oscillation modes with the same \(\numpg\) are found during a gyre calculation. This violates the expectation that \(\numpg\) be monotonic-increasing, and can happen for a few reasons.

Bad Stellar Model

If the input stellar doesn’t conserve mass properly, then one or more bogus (unphysical) modes can appear with the same \(\numpg\) as an existing mode, but a significantly different frequency. Such modes can also be spotted because their radial order is very different from the adjacent-in-frequency modes.

Bogus modes arise because the input stellar model doesn’t conserve mass. GYRE` assumes that the density \(\rho\) and interior mass \(M_{r}\) are related by equation 5. Given that there are many different ways to discretize this equation, there is a certain amount of numerical ‘slop’ that arises when going from discrete form assumed in the stellar evolution code that generated the model (e.g., MESA) through to the discrete form implicitly assumed by GYRE. This slop can sometimes produce bogus modes, and the fix is to re-create the model with a higher spatial resolution.

Another possibility arises if the input stellar model contains density discontinuities that aren’t properly marked using double points. Often, these discontinuities can be diagnosed by plotting the \(I_0\) or \(I_1\) first integrals (for radial and dipole modes, respectively) as a function of radius. These first integrals should vanish everywhere (as shown by Takata 2006a), but will typically show abrupt jumps to non-zero values at the location of unmarked discontinuities. The fix is to re-create the model with double points inserted as necessary; in the case of MESA models, this can be achieved using the add_double_points_to_pulse_data parameter of the &controls namelist group.

Non-adiabatic Effects

When non-adiabatic effects cause a mode to be mis-classified (as discussed in the Missing Modes section), often the incorrect \(\numpg\) value duplicates that of another mode. As before, a mis-classification must be fixed manually by determining which adiabatic mode the problematic non-adiabatic mode corresponds to.