<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_default"><div class="gmail_default" style="font-size:large">Dear Cactus users,  </div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">I am trying to simulate a precessing BBH system. When I analysed the mass of the black hole (qlm_mass), I found that the spin evolution (qlm_spin) is abnormal.</div><div class="gmail_default" style="font-size:large"><img src="cid:ii_kntpgi333" alt="qlm_spin.png" width="432" height="288"><br></div><div class="gmail_default" style="font-size:large">The coordinate-base spin vector gives a value closer to what I set.</div><div class="gmail_default" style="font-size:large"><img src="cid:ii_kntoza9a2" alt="qlm_coordspin_norm.png" width="432" height="288"><br></div><div class="gmail_default" style="font-size:large">This is the initial parameters.</div><div class="gmail_default" style="font-size:large"><pre style="box-sizing:unset;font-size:14px;border:none;margin-top:0px;margin-bottom:0px;padding:0px;overflow:auto;word-break:break-all;word-wrap:break-word;white-space:pre-wrap">ADMMass = 0.9898872694943947
Mass Ratio = 1.0
separation = 10.4
spin1 = [-0.125  0.     0.   ]
spin2 = [0.125 0.    0.   ]
chi1 = [-0.5  0.   0. ]
chi2 = [0.5 0.  0. ]</pre></div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">I tried different methods for finding the killing vector field (i.e. QuasiLocalMeasures::killing_vector_method), but it didn&#39;t work.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">If the small black hole does not spin (spin2 = [0, 0, 0]), things will get better.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">The following is the difference between my parfile and GW150914.rpar</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">$ diff GW150914.rpar BBH000.rpar</div><div class="gmail_default" style="font-size:large">32,33c32,33</div><div class="gmail_default" style="font-size:large">&lt; D = 10.0                    # Separation</div><div class="gmail_default" style="font-size:large">&lt; q = 36.0/29.0               # Mass ratio: q = mp/mm &gt;= 1</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; D = 10.4                    # Separation</div><div class="gmail_default" style="font-size:large">&gt; q = 1              # Mass ratio: q = mp/mm &gt;= 1</div><div class="gmail_default" style="font-size:large">35,38c35,38</div><div class="gmail_default" style="font-size:large">&lt; chip = [0, 0,  0.31]        # Dimensionsless spin of + BH (x0 &gt; 0, more massive)</div><div class="gmail_default" style="font-size:large">&lt; chim = [0, 0, -0.46]        # Dimensionsless spin of - BH (x0 &lt; 0, less massive)</div><div class="gmail_default" style="font-size:large">&lt; Pr = -0.00084541526517121;  # Radial linear momentum</div><div class="gmail_default" style="font-size:large">&lt; Pphi = 0.09530152296974252; # Azimuthal linear momentum</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; chip = [-0.5, 0, 0]        # Dimensionsless spin of + BH (x0 &gt; 0, more massive)</div><div class="gmail_default" style="font-size:large">&gt; chim = [0.5, 0, 0]        # Dimensionsless spin of - BH (x0 &lt; 0, less massive)</div><div class="gmail_default" style="font-size:large">&gt; Pr = -0.000878483532190008</div><div class="gmail_default" style="font-size:large">&gt; Pphi = 0.0927554222753844</div><div class="gmail_default" style="font-size:large">261,265c261,265</div><div class="gmail_default" style="font-size:large">&lt; CoordinatesSymmetry::reflection_z       = yes</div><div class="gmail_default" style="font-size:large">&lt; CoordinatesSymmetry::stagger            = no</div><div class="gmail_default" style="font-size:large">&lt; Coordinates::symmetry                   = &quot;+z bitant&quot;</div><div class="gmail_default" style="font-size:large">&lt; Coordinates::additional_symmetry_size   = 1</div><div class="gmail_default" style="font-size:large">&lt; Coordinates::verbose                    = no</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; ##CoordinatesSymmetry::reflection_z       = yes</div><div class="gmail_default" style="font-size:large">&gt; ##CoordinatesSymmetry::stagger            = no</div><div class="gmail_default" style="font-size:large">&gt; ##Coordinates::symmetry                   = &quot;+z bitant&quot;</div><div class="gmail_default" style="font-size:large">&gt; ##Coordinates::additional_symmetry_size   = 1</div><div class="gmail_default" style="font-size:large">&gt; ##Coordinates::verbose                    = no</div><div class="gmail_default" style="font-size:large">552,553c552,553</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::maxntheta             = 66</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::maxnphi               = 124</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::maxntheta             = 161</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::maxnphi               = 320</div><div class="gmail_default" style="font-size:large">559,560c559,560</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::ntheta            [2] = 41</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::nphi              [2] = 80</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::ntheta            [2] = 161</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::nphi              [2] = 320</div><div class="gmail_default" style="font-size:large">565,566c565,566</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::ntheta            [3] = 41</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::nphi              [3] = 80</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::ntheta            [3] = 161</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::nphi              [3] = 320</div><div class="gmail_default" style="font-size:large">571,572c571,572</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::ntheta            [4] = 41</div><div class="gmail_default" style="font-size:large">&lt; SphericalSurface::nphi              [4] = 80</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::ntheta            [4] = 161</div><div class="gmail_default" style="font-size:large">&gt; SphericalSurface::nphi              [4] = 320</div><div class="gmail_default" style="font-size:large">633,634c633,638</div><div class="gmail_default" style="font-size:large">&lt; IOScalar::outScalar_reductions          = &quot;minimum maximum average&quot;</div><div class="gmail_default" style="font-size:large">&lt; IOScalar::outScalar_vars                = &quot;SystemStatistics::process_memory_mb&quot;</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; IOScalar::outScalar_reductions          = &quot;minimum maximum average norm1 norm2&quot;</div><div class="gmail_default" style="font-size:large">&gt; IOScalar::outScalar_vars                = &quot;</div><div class="gmail_default" style="font-size:large">&gt;   SystemStatistics::process_memory_mb</div><div class="gmail_default" style="font-size:large">&gt;   ML_ADMConstraints::ML_Ham</div><div class="gmail_default" style="font-size:large">&gt;   ML_ADMConstraints::ML_mom</div><div class="gmail_default" style="font-size:large">&gt; &quot;</div><div class="gmail_default" style="font-size:large">670c674</div><div class="gmail_default" style="font-size:large">&lt; IOHDF5::out2D_every                     = $out2d_every</div><div class="gmail_default" style="font-size:large">---</div><div class="gmail_default" style="font-size:large">&gt; IOHDF5::out2D_every                     = 0</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Thank You</div><div class="gmail_default" style="font-size:large">Best Regards</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">---------------------------</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Yu Liu</div><div class="gmail_default" style="font-size:large">PhD Student</div><div class="gmail_default" style="font-size:large">Department of Astronomy</div><div class="gmail_default" style="font-size:large">Huazhong University of Science and Technology</div></div></div></div></div></div></div></div></div></div></div>