strange behavior of l=1 ridge for Procyon

General discussion of all things GYRE-related (e.g., results, talks, ideas, tips)
Post Reply
mjoyce
Posts: 6
Joined: Tue Oct 30, 2018 6:42 pm

strange behavior of l=1 ridge for Procyon

Post by mjoyce » Mon Sep 09, 2019 5:47 am

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
Attachments
model.in
(2.04 KiB) Downloaded 16 times

warrick
Posts: 57
Joined: Wed Aug 28, 2013 2:47 am

Re: strange behavior of l=1 ridge for Procyon

Post by warrick » Mon Sep 30, 2019 4:59 am

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
Posts: 57
Joined: Wed Aug 28, 2013 2:47 am

Re: strange behavior of l=1 ridge for Procyon

Post by warrick » Mon Sep 30, 2019 6:13 am

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
Attachments
T_G1.png
T_G1.png (91.05 KiB) Viewed 523 times
x_G1.png
x_G1.png (94.46 KiB) Viewed 523 times

User avatar
rhtownsend
Site Admin
Posts: 315
Joined: Sun Mar 31, 2013 4:22 pm

Re: strange behavior of l=1 ridge for Procyon

Post by rhtownsend » Mon Sep 30, 2019 8:14 am

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

mjoyce
Posts: 6
Joined: Tue Oct 30, 2018 6:42 pm

Re: strange behavior of l=1 ridge for Procyon

Post by mjoyce » Thu Oct 03, 2019 10:46 pm

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
Posts: 6
Joined: Tue Oct 30, 2018 6:42 pm

Re: strange behavior of l=1 ridge for Procyon

Post by mjoyce » Thu Oct 17, 2019 10:00 pm

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

gamma_vs_x_DSEP_rescale5.png
gamma_vs_x_DSEP_rescale5.png (36.63 KiB) Viewed 4 times
The attachment gamma_vs_x_DSEP_rescale5.png is no longer available

mjoyce
Posts: 6
Joined: Tue Oct 30, 2018 6:42 pm

Re: strange behavior of l=1 ridge for Procyon

Post by mjoyce » Thu Oct 17, 2019 10:01 pm

DSEP_rescale5.png
DSEP_rescale5.png (67.58 KiB) Viewed 4 times

Post Reply