Tidal coupling constants with MESA/GYRE

General discussion of all things GYRE-related (e.g., results, talks, ideas, tips)
Post Reply
alekseygenerozov
Posts: 2
Joined: Wed Aug 01, 2018 11:04 am

Tidal coupling constants with MESA/GYRE

Post by alekseygenerozov » Wed Aug 01, 2018 12:36 pm

Hello. I am interested in tidal interactions due to close encounters between stars in clusters.

As described in Press&Teukolsky 1977 (PT77) the energy deposited into a star following a close encounter can be encapsulated by a “tidal coupling constant” (really a function of the encounter distance and stellar model). I am currently working on tabulating tidal coupling constants for a broad range of stellar models (see the github repository here).

For each star the process is: (a) construct a MESA model (b) calculate the spectrum of normal modes with GYRE (c) calculate the overlap integrals for each mode. The overlap integral for a particular mode is

$Q= \int_0^1 r^2 \rho \ell r^{\ell-1} [ \xi^R + (\ell+1) \xi^H] dr$,
where rho is the stellar density, and xi^R (xi^H) are the radial (tangential) components of the mode (renormalized following PT77). Computing the tidal coupling constants involves a summation over different modes. So far I just have the quadrupole order (l=2) coupling constants. (Apologies I could not figure out how to render the above equation inline)

I am using MESA v. 9793, GYRE v. 5.0. I have attached sample in-lists for both.

Hopefully these tables will eventually be incorporated into the FEWBODY code!

I would appreciate any feedback on this project, particularly on the following two points.

(1) As I mention in the README here the overlap integrals of massive stars do not vary smoothly between consecutive, high order g-modes. I have tried playing with some of numerical parameters in the MESA/GYRE inlists (e.g. grid resolution), but this does not change the result. The experiments are described in the linked repo. Overall, the behavior of the overlap integrals for high order g-modes looks fish to me. But maybe it is physical? Do you have any insight on this?

(2) We would be interested in doing this for red giants also. I was wondering if any pit-falls are immediately apparent. (A collaborator expressed a concern that calculating the tidal coupling constants for red giants would be difficult to do robustly due to dense oscillation frequencies and small amplitudes).

Apologies if the above is somewhat open-ended...
Attachments
gyre_20Msun.in
(886 Bytes) Downloaded 10 times

fffeynman
Posts: 3
Joined: Wed Mar 18, 2015 2:12 pm
Location: Penn State University
Contact:

Re: Tidal coupling constants with MESA/GYRE

Post by fffeynman » Thu Aug 09, 2018 2:58 pm

Hi, Aleksey?

Calculating the overlap integral Q_nl directly from the displacement eigenfunctions will inevitably suffer from the numerical inaccuracies, especially for
high order g-modes. One way to circumvent the problem is to directly use the Eulerian perturbation of gravitation potential (see eq. 9 in Burkart et al 2012, MNRAS, 421, 983), this is is also one of the outputs that GYRE can offer, without adopting the Cowling approximation of course.
Hope it helps.

best regards,
Zhao

alekseygenerozov
Posts: 2
Joined: Wed Aug 01, 2018 11:04 am

Re: Tidal coupling constants with MESA/GYRE

Post by alekseygenerozov » Sun Aug 12, 2018 10:24 pm

This is very helpful. Thank you! It is definitely a more straightforward way to get the overlap integrals. Although, the results are similar to what I get using the displacement eigenfunction directly. The exception is high order p-modes in low mass stars.

l=2 , 20 Msun
comp_0.3_2.png
comp_0.3_2.png (22.15 KiB) Viewed 61 times
l=2, 0.3 Msun
The attachment comp_0.3_2.png is no longer available

Post Reply