Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALSimulation 6.2.0.1-b246709
•All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALSimIMRPhenomXHM_qnm.c
Go to the documentation of this file.
1
2/*
3 * Copyright (C) 2019 Marta Colleoni, Cecilio García Quirós
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 2 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with with program; see the file COPYING. If not, write to the
17 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
18 * MA 02110-1301 USA
19 */
20
21 /*
22 This file contains the functions for the fits of the ringdown and damping frequencies for each mode
23 and the fits for the real and imaginary part of the mixing coefficients (only for the 32 mode).
24 The fits build on data taken from https://pages.jh.edu/~eberti2/ringdown/, see also https://arxiv.org/abs/gr-qc/0512160 and https://arxiv.org/abs/1408.1860. Explicit expressions are given in the supplementary material, in IMRPhenomXHM_PhaseFits.nb, as well as in the dedicated folders QNMs and MixingCoefficients, see https://dcc.ligo.org/LIGO-P2000011 and https://arxiv.org/abs/2001.10914
25 */
26
27#include "LALSimIMRPhenomXHM_qnm.h"
28
29
30double evaluate_QNMfit_fring21(double finalDimlessSpin){
31
32 double return_val;
33
34 if (fabs(finalDimlessSpin) > 1.0) {
35 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fring21 function: |finalDimlessSpin| > 1.0 not supported");
36 }
37
38 double x2= finalDimlessSpin*finalDimlessSpin;
39 double x3= x2*finalDimlessSpin;
40 double x4= x2*x2;
41 double x5= x3*x2;
42
43 return_val = (0.059471695665734674 - 0.07585416297991414*finalDimlessSpin + 0.021967909664591865*x2 - 0.0018964744613388146*x3 + 0.001164879406179587*x4 - 0.0003387374454044957*x5)/(1 - 1.4437415542456158*finalDimlessSpin + 0.49246920313191234*x2);
44 return return_val;
45}
46
47double evaluate_QNMfit_fdamp21(double finalDimlessSpin){
48
49 double return_val;
50
51 if (fabs(finalDimlessSpin) > 1.0) {
52 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fdamp21 function: |finalDimlessSpin| > 1.0 not supported");
53 }
54
55 double x2= finalDimlessSpin*finalDimlessSpin;
56 double x3= x2*finalDimlessSpin;
57 double x4= x2*x2;
58 double x5= x3*x2;
59
60 return_val = (2.0696914454467294 - 3.1358071947583093*finalDimlessSpin + 0.14456081596393977*x2 + 1.2194717985037946*x3 - 0.2947372598589144*x4 + 0.002943057145913646*x5)/(146.1779212636481 - 219.81790388304876*finalDimlessSpin + 17.7141194900164*x2 + 75.90115083917898*x3 - 18.975287709794745*x4);
61 return return_val;
62}
63
64double evaluate_QNMfit_fring33(double finalDimlessSpin){
65
66 double return_val;
67
68 if (fabs(finalDimlessSpin) > 1.0) {
69 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fring33 function: |finalDimlessSpin| > 1.0 not supported");
70 }
71
72 double x2= finalDimlessSpin*finalDimlessSpin;
73 double x3= x2*finalDimlessSpin;
74 double x4= x2*x2;
75 double x5= x3*x2;
76 double x6= x3*x3;
77
78 return_val = (0.09540436245212061 - 0.22799517865876945*finalDimlessSpin + 0.13402916709362475*x2 + 0.03343753057911253*x3 - 0.030848060170259615*x4 - 0.006756504382964637*x5 + 0.0027301732074159835*x6)/(1 - 2.7265947806178334*finalDimlessSpin + 2.144070539525238*x2 - 0.4706873667569393*x4 + 0.05321818246993958*x6);
79 return return_val;
80}
81
82double evaluate_QNMfit_fdamp33(double finalDimlessSpin){
83
84 double return_val;
85
86 if (fabs(finalDimlessSpin) > 1.0) {
87 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fdamp33 function: |finalDimlessSpin| > 1.0 not supported");
88 }
89
90 double x2= finalDimlessSpin*finalDimlessSpin;
91 double x3= x2*finalDimlessSpin;
92 double x4= x2*x2;
93 double x5= x3*x2;
94
95 return_val = (0.014754148319335946 - 0.03124423610028678*finalDimlessSpin + 0.017192623913708124*x2 + 0.001034954865629645*x3 - 0.0015925124814622795*x4 - 0.0001414350555699256*x5)/(1 - 2.0963684630756894*finalDimlessSpin + 1.196809702382645*x2 - 0.09874113387889819*x4);
96 return return_val;
97}
98
99double evaluate_QNMfit_fring32(double finalDimlessSpin){
100
101 double return_val;
102
103 if (fabs(finalDimlessSpin) > 1.0) {
104 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fring32 function: |finalDimlessSpin| > 1.0 not supported");
105 }
106
107 double x2= finalDimlessSpin*finalDimlessSpin;
108 double x3= x2*finalDimlessSpin;
109 double x4= x2*x2;
110 double x5= x3*x2;
111 double x6= x3*x3;
112
113 return_val = (0.09540436245212061 - 0.13628306966373951*finalDimlessSpin + 0.030099881830507727*x2 - 0.000673589757007597*x3 + 0.0118277880067919*x4 + 0.0020533816327907334*x5 - 0.0015206141948469621*x6)/(1 - 1.6531854335715193*finalDimlessSpin + 0.5634705514193629*x2 + 0.12256204148002939*x4 - 0.027297817699401976*x6);
114 return return_val;
115}
116
117double evaluate_QNMfit_fdamp32(double finalDimlessSpin){
118
119 double return_val;
120
121 if (fabs(finalDimlessSpin) > 1.0) {
122 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fdamp32 function: |finalDimlessSpin| > 1.0 not supported");
123 }
124
125 double x2= finalDimlessSpin*finalDimlessSpin;
126 double x3= x2*finalDimlessSpin;
127 double x4= x2*x2;
128
129 return_val = (0.014754148319335946 - 0.03445752346074498*finalDimlessSpin + 0.02168855041940869*x2 + 0.0014945908223317514*x3 - 0.0034761714223258693*x4)/(1 - 2.320722660848874*finalDimlessSpin + 1.5096146036915865*x2 - 0.18791187563554512*x4);
130 return return_val;
131}
132
133double evaluate_QNMfit_fring44(double finalDimlessSpin){
134
135 double return_val;
136
137 if (fabs(finalDimlessSpin) > 1.0) {
138 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fring44 function: |finalDimlessSpin| > 1.0 not supported");
139 }
140
141 double x2= finalDimlessSpin*finalDimlessSpin;
142 double x3= x2*finalDimlessSpin;
143 double x4= x2*x2;
144 double x5= x3*x2;
145 double x6= x3*x3;
146
147 return_val = (0.1287821193485683 - 0.21224284094693793*finalDimlessSpin + 0.0710926778043916*x2 + 0.015487322972031054*x3 - 0.002795401084713644*x4 + 0.000045483523029172406*x5 + 0.00034775290179000503*x6)/(1 - 1.9931645124693607*finalDimlessSpin + 1.0593147376898773*x2 - 0.06378640753152783*x4);
148 return return_val;
149}
150
151double evaluate_QNMfit_fdamp44(double finalDimlessSpin){
152
153 double return_val;
154
155 if (fabs(finalDimlessSpin) > 1.0) {
156 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_fdamp44 function: |finalDimlessSpin| > 1.0 not supported");
157 }
158
159 double x2= finalDimlessSpin*finalDimlessSpin;
160 double x3= x2*finalDimlessSpin;
161 double x4= x2*x2;
162 double x5= x3*x2;
163 double x6= x3*x3;
164
165 return_val = (0.014986847152355699 - 0.01722587715950451*finalDimlessSpin - 0.0016734788189065538*x2 + 0.0002837322846047305*x3 + 0.002510528746148588*x4 + 0.00031983835498725354*x5 + 0.000812185411753066*x6)/(1 - 1.1350205970682399*finalDimlessSpin - 0.0500827971270845*x2 + 0.13983808071522857*x4 + 0.051876225199833995*x6);
166 return return_val;
167}
168
169
170double evaluate_QNMfit_re_l2m2lp2(double finalDimlessSpin){
171
172 double return_val;
173
174 if (fabs(finalDimlessSpin) > 1.0) {
175 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_re_l2m2lp2 function: |finalDimlessSpin| > 1.0 not supported");
176 }
177
178 double x2= finalDimlessSpin*finalDimlessSpin;
179 double x3= x2*finalDimlessSpin;
180 double x4= x2*x2;
181 double x5= x3*x2;
182 double x6= x3*x3;
183
184 return_val = (1 - 2.2956993576253635*finalDimlessSpin + 1.461988775298876*x2 + 0.0043296365593147035*x3 - 0.1695667458204109*x4 - 0.0006267849034466508*x5)/(1 - 2.2956977727459043*finalDimlessSpin + 1.4646339137818438*x2 - 0.16843226886562457*x4 - 0.00007150540890128118*x6);
185 return return_val;
186}
187
188double evaluate_QNMfit_im_l2m2lp2(double finalDimlessSpin){
189
190 double return_val;
191
192 if (fabs(finalDimlessSpin) > 1.0) {
193 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_im_l2m2lp2 function: |finalDimlessSpin| > 1.0 not supported");
194 }
195
196 double x2= finalDimlessSpin*finalDimlessSpin;
197 double x3= x2*finalDimlessSpin;
198 double x4= x2*x2;
199 double x5= x3*x2;
200 double x6= x3*x3;
201
202 return_val = (finalDimlessSpin*(0.3826673013161342 - 0.47531267226013896*finalDimlessSpin - 0.05898102880105067*x2 + 0.0724525431346487*x3 + 0.054714637311702986*x4 + 0.024544862718252784*x5))/(-38.70835035062785 + 69.82140084545878*finalDimlessSpin - 27.99036444363243*x2 - 4.152310472191899*x4 + 1.*x6);
203 return return_val;
204}
205
206double evaluate_QNMfit_re_l3m2lp2(double finalDimlessSpin){
207
208 double return_val;
209
210 if (fabs(finalDimlessSpin) > 1.0) {
211 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_re_l3m2lp2 function: |finalDimlessSpin| > 1.0 not supported");
212 }
213
214 double x2= finalDimlessSpin*finalDimlessSpin;
215 double x3= x2*finalDimlessSpin;
216 double x4= x2*x2;
217 double x5= x3*x2;
218
219 return_val = (finalDimlessSpin*(0.47513455283841244 - 0.9016636384605536*finalDimlessSpin + 0.3844811236426182*x2 + 0.0855565148647794*x3 - 0.03620067426672167*x4 - 0.006557249133752502*x5))/(-6.76894063440646 + 15.170831931186493*finalDimlessSpin - 9.406169787571082*x2 + 1.*x4);
220 return return_val;
221}
222
223double evaluate_QNMfit_im_l3m2lp2(double finalDimlessSpin){
224
225 double return_val;
226
227 if (fabs(finalDimlessSpin) > 1.0) {
228 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_im_l3m2lp2 function: |finalDimlessSpin| > 1.0 not supported");
229 }
230
231 double x2= finalDimlessSpin*finalDimlessSpin;
232 double x3= x2*finalDimlessSpin;
233 double x4= x2*x2;
234 double x5= x3*x2;
235 double x6= x3*x3;
236
237 return_val = (finalDimlessSpin*(-2.8704762147145533 + 4.436434016918535*finalDimlessSpin - 1.0115343326360486*x2 - 0.08965314412106505*x3 - 0.4236810894599512*x4 - 0.041787576033810676*x5))/(-171.80908957903395 + 272.362882450877*finalDimlessSpin - 76.68544453077854*x2 - 25.14197656531123*x4 + 1.*x6);
238 return return_val;
239}
240
241double evaluate_QNMfit_re_l2m2lp3(double finalDimlessSpin){
242
243 double return_val;
244
245 if (fabs(finalDimlessSpin) > 1.0) {
246 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_re_l2m2lp3 function: |finalDimlessSpin| > 1.0 not supported");
247 }
248
249 double x2= finalDimlessSpin*finalDimlessSpin;
250 double x3= x2*finalDimlessSpin;
251 double x4= x2*x2;
252 double x5= x3*x2;
253 double x6= x3*x3;
254
255 return_val = (finalDimlessSpin*(18.522563276099167 - 37.978140351289014*finalDimlessSpin + 19.030390708998894*x2 + 3.0355668591803386*x3 - 2.210028290847915*x4 - 0.37117112862247975*x5))/(164.52480238697507 - 377.9093045285145*finalDimlessSpin + 243.3353695550844*x2 - 30.79738566181734*x4 + 1.*x6);
256 return return_val;
257}
258
259double evaluate_QNMfit_im_l2m2lp3(double finalDimlessSpin){
260
261 double return_val;
262
263 if (fabs(finalDimlessSpin) > 1.0) {
264 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_im_l2m2lp3 function: |finalDimlessSpin| > 1.0 not supported");
265 }
266
267 double x2= finalDimlessSpin*finalDimlessSpin;
268 double x3= x2*finalDimlessSpin;
269 double x4= x2*x2;
270 double x5= x3*x2;
271 double x6= x3*x3;
272
273 return_val = (finalDimlessSpin*(-49.7688437256778 + 120.43773704442333*finalDimlessSpin - 82.95323455645332*x2 + 1.721453011852496*x3 + 11.540237244397877*x4 - 0.9819458637589314*x5))/(2858.5790831181725 - 6305.619505422591*finalDimlessSpin + 3825.6742092829054*x2 - 377.7822297815406*x4 + 1.*x6);
274 return return_val;
275}
276
277double evaluate_QNMfit_re_l3m2lp3(double finalDimlessSpin){
278
279 double return_val;
280
281 if (fabs(finalDimlessSpin) > 1.0) {
282 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_re_l3m2lp3 function: |finalDimlessSpin| > 1.0 not supported");
283 }
284
285 double x2= finalDimlessSpin*finalDimlessSpin;
286 double x3= x2*finalDimlessSpin;
287 double x4= x2*x2;
288 double x5= x3*x2;
289 double x6= x3*x3;
290
291 return_val = (1 - 2.107852425643677*finalDimlessSpin + 1.1906393634562715*x2 + 0.02244848864087732*x3 - 0.09593447799423722*x4 - 0.0021343381708933025*x5 - 0.005319515989331159*x6)/(1 - 2.1078515887706324*finalDimlessSpin + 1.2043484690080966*x2 - 0.08910191596778137*x4 - 0.005471749827809503*x6);
292 return return_val;
293}
294
295double evaluate_QNMfit_im_l3m2lp3(double finalDimlessSpin){
296
297 double return_val;
298
299 if (fabs(finalDimlessSpin) > 1.0) {
300 XLAL_ERROR(XLAL_EDOM, "PhenomXHM evaluate_QNMfit_im_l3m2lp3 function: |finalDimlessSpin| > 1.0 not supported");
301 }
302
303 double x2= finalDimlessSpin*finalDimlessSpin;
304 double x3= x2*finalDimlessSpin;
305 double x4= x2*x2;
306 double x5= x3*x2;
307 double x6= x3*x3;
308
309 return_val = (finalDimlessSpin*(12.45701482868677 - 29.398484595717147*finalDimlessSpin + 18.26221675782779*x2 + 1.9308599142669403*x3 - 3.159763242921214*x4 - 0.0910871567367674*x5))/(345.52914639836257 - 815.4349339779621*finalDimlessSpin + 538.3888932415709*x2 - 69.3840921447381*x4 + 1.*x6);
310 return return_val;
311}
evaluate_QNMfit_fdamp21
double evaluate_QNMfit_fdamp21(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:47
evaluate_QNMfit_re_l2m2lp3
double evaluate_QNMfit_re_l2m2lp3(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:241
evaluate_QNMfit_re_l2m2lp2
double evaluate_QNMfit_re_l2m2lp2(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:170
evaluate_QNMfit_fdamp44
double evaluate_QNMfit_fdamp44(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:151
evaluate_QNMfit_fring44
double evaluate_QNMfit_fring44(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:133
evaluate_QNMfit_fring33
double evaluate_QNMfit_fring33(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:64
evaluate_QNMfit_im_l2m2lp3
double evaluate_QNMfit_im_l2m2lp3(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:259
evaluate_QNMfit_fdamp32
double evaluate_QNMfit_fdamp32(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:117
evaluate_QNMfit_im_l3m2lp3
double evaluate_QNMfit_im_l3m2lp3(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:295
evaluate_QNMfit_im_l2m2lp2
double evaluate_QNMfit_im_l2m2lp2(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:188
evaluate_QNMfit_fring21
double evaluate_QNMfit_fring21(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:30
evaluate_QNMfit_re_l3m2lp3
double evaluate_QNMfit_re_l3m2lp3(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:277
evaluate_QNMfit_fdamp33
double evaluate_QNMfit_fdamp33(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:82
evaluate_QNMfit_fring32
double evaluate_QNMfit_fring32(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:99
evaluate_QNMfit_im_l3m2lp2
double evaluate_QNMfit_im_l3m2lp2(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:223
evaluate_QNMfit_re_l3m2lp2
double evaluate_QNMfit_re_l3m2lp2(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:206
LALSimIMRPhenomXHM_qnm.h
XLAL_ERROR
#define XLAL_ERROR(...)
XLAL_EDOM
XLAL_EDOM