strange behavior of l=1 ridge for Procyon

[mjoyce]

Hi All / Rich,
I've been trying to produce synthetic frequency spectra of Procyon over a coarse parameter space (something I know is definitely possible) where I've used DSEP to generate the gonge files. Across all parameter combinations, I'm getting a weird ridge crossing for L = 1 in the Echelle diagrams. To first order, I'm trying to sort out whether this is a DSEP issue or some problem with my model.in controls for GYRE.

I've attached the model.in file and the gonge files. It seems the file uploader won't accept my GYRE output, so I'm pasting it in line below this message...sorry for the ugliness.

Does anyone have any idea what could be going wrong here? Tim White and I have been stuck on this for a while now...any insight is appreciated

-------------------------------------------------------------------------------------------------

1 2 3
M_star R_star freq_units
0.2906075100000000E+034 0.1447774470000000E+012UHZ
1 2 3 4 5 6 7
l n_pg n_p n_g Re(freq) Im(freq) E_norm
0 9 9 0 0.5123709074425948E+003 0.0000000000000000E+000 0.1638005975048221E-006
0 10 10 0 0.5677200563646284E+003 0.0000000000000000E+000 0.1101053363480010E-006
0 11 11 0 0.6238573133424256E+003 0.0000000000000000E+000 0.7569105626923405E-007
0 12 12 0 0.6805952602803055E+003 0.0000000000000000E+000 0.4918306377543802E-007
0 13 13 0 0.7368121555359837E+003 0.0000000000000000E+000 0.3025544999004832E-007
0 14 14 0 0.7922868713680920E+003 0.0000000000000000E+000 0.1869514928567667E-007
0 15 15 0 0.8472858054805251E+003 0.0000000000000000E+000 0.1210015353075965E-007
0 16 16 0 0.9019270379310273E+003 0.0000000000000000E+000 0.8289128904047418E-008
0 17 17 0 0.9568398844355487E+003 0.0000000000000000E+000 0.6078388747108554E-008
0 18 18 0 0.1012504023951596E+004 0.0000000000000000E+000 0.4552238957900554E-008
0 19 19 0 0.1068855305140205E+004 0.0000000000000000E+000 0.3240439435136320E-008
0 20 20 0 0.1125393822426251E+004 0.0000000000000000E+000 0.2173363583355066E-008
0 21 21 0 0.1181923318060494E+004 0.0000000000000000E+000 0.1367740158255394E-008
0 22 22 0 0.1238537434673724E+004 0.0000000000000000E+000 0.7969061339004060E-009
0 23 23 0 0.1294850985539960E+004 0.0000000000000000E+000 0.4260204251263952E-009
0 24 24 0 0.1351016630212948E+004 0.0000000000000000E+000 0.2041375091811279E-009
0 25 25 0 0.1407178625740893E+004 0.0000000000000000E+000 0.8197353915584875E-010
0 26 26 0 0.1463254779193073E+004 0.0000000000000000E+000 0.2348712997467410E-010
0 27 27 0 0.1518824788782096E+004 0.0000000000000000E+000 0.6099599714800318E-011
0 28 28 0 0.1522870204611406E+004 0.0000000000000000E+000 0.5947660361204929E-011
0 29 29 0 0.1575469603024135E+004 0.0000000000000000E+000 0.1378232164386426E-010
0 30 30 0 0.1631725549699907E+004 0.0000000000000000E+000 0.3619465516376979E-010
0 31 31 0 0.1688067739753007E+004 0.0000000000000000E+000 0.6593880543732452E-010
0 32 32 0 0.1744470111601784E+004 0.0000000000000000E+000 0.9715071653551950E-010
0 33 33 0 0.1801084250740879E+004 0.0000000000000000E+000 0.1256434944248514E-009
0 34 34 0 0.1857641363427997E+004 0.0000000000000000E+000 0.1475189150921327E-009
0 35 35 0 0.1914254141103571E+004 0.0000000000000000E+000 0.1614741316396465E-009
0 36 36 0 0.1970885013648154E+004 0.0000000000000000E+000 0.1683230447415073E-009

1 2 3
M_star R_star freq_units
0.2906075100000000E+034 0.1447774470000000E+012UHZ
1 2 3 4 5 6 7
l n_pg n_p n_g Re(freq) Im(freq) E_norm
1 7 7 1 0.5026400153479791E+003 0.0000000000000000E+000 0.1852436537228644E-006
1 8 8 1 0.5529483206869776E+003 0.0000000000000000E+000 0.1357384506439386E-006
1 9 9 1 0.5906931464531120E+003 0.0000000000000000E+000 0.1425588938256885E-006
1 10 10 1 0.6239918843560007E+003 0.0000000000000000E+000 0.8489784761127841E-007
1 11 10 0 0.6739020021825360E+003 0.0000000000000000E+000 0.5363505523976685E-007
1 12 11 0 0.7274941516597828E+003 0.0000000000000000E+000 0.3359688965484655E-007
1 13 12 0 0.7813023855155803E+003 0.0000000000000000E+000 0.2095720774360474E-007
1 14 13 0 0.8348586107055071E+003 0.0000000000000000E+000 0.1359775093498561E-007
1 15 14 0 0.8882469158857090E+003 0.0000000000000000E+000 0.9248776490998965E-008
1 16 15 0 0.9417631700825664E+003 0.0000000000000000E+000 0.6742831674889043E-008
1 17 16 0 0.9959392901983128E+003 0.0000000000000000E+000 0.5090482744066688E-008
1 18 17 0 0.1050984599843924E+004 0.0000000000000000E+000 0.3723430672937868E-008
1 19 18 0 0.1106278068168813E+004 0.0000000000000000E+000 0.2568111881097769E-008
1 20 19 0 0.1161706966610933E+004 0.0000000000000000E+000 0.1664050888716038E-008
1 21 20 0 0.1217235235505930E+004 0.0000000000000000E+000 0.1010780200732689E-008
1 22 21 0 0.1272627575566509E+004 0.0000000000000000E+000 0.5632598752754486E-009
1 23 22 0 0.1327917804670736E+004 0.0000000000000000E+000 0.2864678119013709E-009
1 24 23 0 0.1383091994159162E+004 0.0000000000000000E+000 0.1272927704307665E-009
1 25 24 0 0.1438404257964360E+004 0.0000000000000000E+000 0.4375656326611308E-010
1 26 25 0 0.1493603109246069E+004 0.0000000000000000E+000 0.1020568860167396E-010
1 27 26 0 0.1522319743910312E+004 0.0000000000000000E+000 0.5960943656598610E-011
1 28 27 0 0.1549003298826309E+004 0.0000000000000000E+000 0.7803621396186994E-011
1 29 28 0 0.1604440802556494E+004 0.0000000000000000E+000 0.2432423665950886E-010
1 30 29 0 0.1659996138035360E+004 0.0000000000000000E+000 0.5123339858591108E-010
1 31 30 0 0.1715754614063403E+004 0.0000000000000000E+000 0.8223168277529501E-010
1 32 31 0 0.1771605734010793E+004 0.0000000000000000E+000 0.1127807138473721E-009
1 33 32 0 0.1827613269590294E+004 0.0000000000000000E+000 0.1379645916781473E-009
1 34 33 0 0.1883670731662373E+004 0.0000000000000000E+000 0.1563002586072816E-009
1 35 34 0 0.1939696112410281E+004 0.0000000000000000E+000 0.1667136245998483E-009
1 36 35 0 0.1995911772017383E+004 0.0000000000000000E+000 0.1694556333791317E-009

1 2 3
M_star R_star freq_units
0.2906075100000000E+034 0.1447774470000000E+012UHZ
1 2 3 4 5 6 7
l n_pg n_p n_g Re(freq) Im(freq) E_norm
2 6 8 2 0.5357655492224557E+003 0.0000000000000000E+000 0.1375156168255531E-006
2 8 9 1 0.5718746672482990E+003 0.0000000000000000E+000 0.5295122028334099E-001
2 8 9 1 0.5909005931069748E+003 0.0000000000000000E+000 0.9418968606271006E-007
2 9 10 1 0.6469723890096532E+003 0.0000000000000000E+000 0.6383466316837104E-007
2 10 11 1 0.7033109079592767E+003 0.0000000000000000E+000 0.4026205787384602E-007
2 11 12 1 0.7590492573089913E+003 0.0000000000000000E+000 0.2466893332334293E-007
2 12 13 1 0.8140986373690509E+003 0.0000000000000000E+000 0.1551419035505534E-007
2 12 13 1 0.8569385416116281E+003 0.0000000000000000E+000 0.2060837419331491E-002
2 14 14 0 0.8688620048831914E+003 0.0000000000000000E+000 0.1030356231787981E-007
2 15 15 0 0.9234198162668071E+003 0.0000000000000000E+000 0.7221902241261261E-008
2 16 16 0 0.9785053114553768E+003 0.0000000000000000E+000 0.5374172796094853E-008
2 17 17 0 0.1034494608364638E+004 0.0000000000000000E+000 0.3970968631750999E-008
2 18 18 0 0.1090855397415137E+004 0.0000000000000000E+000 0.2761362359403025E-008
2 19 19 0 0.1147394560733466E+004 0.0000000000000000E+000 0.1805107599059839E-008
2 20 20 0 0.1203923889163793E+004 0.0000000000000000E+000 0.1107772466181238E-008
2 21 21 0 0.1260400169834894E+004 0.0000000000000000E+000 0.6235546576555902E-009
2 22 22 0 0.1316645312258335E+004 0.0000000000000000E+000 0.3214020810083703E-009
2 23 23 0 0.1372703273514133E+004 0.0000000000000000E+000 0.1459408483246729E-009
2 24 24 0 0.1428862289658499E+004 0.0000000000000000E+000 0.5261871208260846E-010
2 25 25 0 0.1484843499456604E+004 0.0000000000000000E+000 0.1288214324972144E-010
2 26 26 0 0.1522264123537352E+004 0.0000000000000000E+000 0.5962408953708333E-011
2 27 27 0 0.1541001333170352E+004 0.0000000000000000E+000 0.6712186041639706E-011
2 28 28 0 0.1597046307385264E+004 0.0000000000000000E+000 0.2092123600375919E-010
2 29 29 0 0.1653226472644892E+004 0.0000000000000000E+000 0.4651947693665262E-010
2 30 30 0 0.1709567300506241E+004 0.0000000000000000E+000 0.7700429241409206E-010
2 31 31 0 0.1765932595136182E+004 0.0000000000000000E+000 0.1074755237684011E-009
2 32 32 0 0.1822458423401692E+004 0.0000000000000000E+000 0.1331799412483206E-009
2 33 33 0 0.1878963364477321E+004 0.0000000000000000E+000 0.1521705223310330E-009
2 34 34 0 0.1935418544839363E+004 0.0000000000000000E+000 0.1632095449410302E-009
2 35 35 0 0.1992013205872309E+004 0.0000000000000000E+000 0.1668031738983728E-009

1 2 3
M_star R_star freq_units
0.2906075100000000E+034 0.1447774470000000E+012UHZ
1 2 3 4 5 6 7
l n_pg n_p n_g Re(freq) Im(freq) E_norm
3 4 7 3 0.5002748506923064E+003 0.0000000000000000E+000 0.1794106858808200E-006
3 5 7 2 0.5327402236844745E+003 0.0000000000000000E+000 0.4315993172211288E+003
3 6 8 2 0.5552820080403242E+003 0.0000000000000000E+000 0.1188058683776491E-006
3 7 9 2 0.6110887898287980E+003 0.0000000000000000E+000 0.8175790441175361E-007
3 8 10 2 0.6676694971526381E+003 0.0000000000000000E+000 0.5409108231124809E-007
3 9 10 1 0.6891226030632329E+003 0.0000000000000000E+000 0.1683359396630890E+001
3 10 11 1 0.7239883241987988E+003 0.0000000000000000E+000 0.3347495654733549E-007
3 11 12 1 0.7796254607957256E+003 0.0000000000000000E+000 0.2055355291806002E-007
3 12 13 1 0.8346719040480923E+003 0.0000000000000000E+000 0.1315331542699409E-007
3 13 14 1 0.8894055446139015E+003 0.0000000000000000E+000 0.8858331823131198E-008
3 14 15 1 0.9442214434983752E+003 0.0000000000000000E+000 0.6412480285597880E-008
3 15 15 0 0.9671268296674061E+003 0.0000000000000000E+000 0.5174714972232489E-002
3 16 16 0 0.9997070897602911E+003 0.0000000000000000E+000 0.4802904194477925E-008
3 17 17 0 0.1056046316626060E+004 0.0000000000000000E+000 0.3467468124279261E-008
3 18 18 0 0.1112586324602794E+004 0.0000000000000000E+000 0.2357097298503845E-008
3 19 19 0 0.1169193218925069E+004 0.0000000000000000E+000 0.1502563585237908E-008
3 20 20 0 0.1225875876658159E+004 0.0000000000000000E+000 0.8930245231941537E-009
3 21 21 0 0.1282305147141525E+004 0.0000000000000000E+000 0.4861794876544827E-009
3 22 22 0 0.1338591355561462E+004 0.0000000000000000E+000 0.2394625123134152E-009
3 23 23 0 0.1394769382391677E+004 0.0000000000000000E+000 0.1011451690342211E-009
3 24 24 0 0.1450959873830528E+004 0.0000000000000000E+000 0.3171994633908132E-010
3 25 25 0 0.1506995602281297E+004 0.0000000000000000E+000 0.7286173705468221E-011
3 26 26 0 0.1522421027115416E+004 0.0000000000000000E+000 0.5958313384906718E-011
3 27 27 0 0.1563244314449718E+004 0.0000000000000000E+000 0.1045669092912183E-010
3 28 28 0 0.1619531174003682E+004 0.0000000000000000E+000 0.3010296028336205E-010
3 29 29 0 0.1675856089340603E+004 0.0000000000000000E+000 0.5834584845253979E-010
3 30 30 0 0.1732297302747652E+004 0.0000000000000000E+000 0.8911718527095148E-010
3 31 31 0 0.1788882937604047E+004 0.0000000000000000E+000 0.1183055154113171E-009
3 32 32 0 0.1845456982620287E+004 0.0000000000000000E+000 0.1412506625060274E-009
3 33 33 0 0.1902094533302940E+004 0.0000000000000000E+000 0.1568389790386518E-009
3 34 34 0 0.1958670461639939E+004 0.0000000000000000E+000 0.1651596381174889E-009

[warrick]

Hi Meridith,

Could you send the stellar model file? When oscillations look weird, I find the equilibrium model is more often at fault than GYRE (or any other oscillation code).

Also, what version of GYRE are you using? The attached input file has the old `op_type` control, so I guess you're using something from before 5. You could try upgrading to 5.x, though I have a feeling that won't help.

Cheers,
Warrick

[warrick]

Hi again,

I pestered Tim for a stellar model and he sent me a different model that shows a similar problem. I've been poking around and there are three places where the dimensionless Ledoux discriminant (A=N²r/g) is fishy. One is where the atmosphere is joined onto the model but I tried truncating the atmosphere to no effect, so I doubt that's to blame. There's also the convective core boundary, which is typically messy, since I presume this convective core has receded as the model evolved.

Finally, there's a third region where there are two points, seemingly in the middle of nowhere, around x=r/R_*=0.696, where the Ledoux discriminant jumps up for two points. Looking at the components of A, I see there's a sudden drop in the first adiabatic index. This might not be causing any problems in the oscillations but seems fairly nasty. Is this maybe where two EoS tables are linked? It happens to coincide with the point where the temperature is exactly 10^6 K. See if your problematic models show the same thing.

I was able to compute mode frequencies with GYRE 5.2 by the way. If you're getting something about "Non-monotonic data", check that the model is written in the correct order. i.e. from surface to centre. I think GYRE presumes this is the case.

Cheers,
Warrick

[rhtownsend]

That jump in A is definitely fishy -- I suggest the first place to start is to look into that.

[mjoyce]

Thanks both for the investigative efforts (and for bothering Tim W). It may simply be the case that DSEP is not equipped to operate in this parameter space, or that I've linked some external data tables incorrectly. Really appreciate the feedback about where the structure looks odd...!

[mjoyce]

I think we're doing a bit better...?

gamma_vs_x_DSEP_rescale5.png (36.63 KiB) Viewed 5734 times

The attachment gamma_vs_x_DSEP_rescale5.png is no longer available

[mjoyce]

DSEP_rescale5.png (67.58 KiB) Viewed 5734 times

[warrick]

Definitely! The Gamma_1 profile doesn't have that nasty glitch any more and the echelle diagram looks OK.

[rhtownsend]

Much better! How did you fix the discontinuity?

[mjoyce]

After much emailing, I learned that some extra work is required to force DSEP to push the envelope--interior matching point outwards. I don't have an automated solution yet, but a control file that forces several additional rescaling runs during the prems evolution seems to solve the problem about 50% of the time. More discoveries to come, I hope...