<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div><br></div><div>Thanks for the help Barry, Ian and Erik.</div><div><br data-mce-bogus="1"></div><div>Cheers,</div><div>Francisco</div><div><br></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>From: </b>"Barry Wardell" &lt;barry.wardell@gmail.com&gt;<br><b>To: </b>"Ian Hinder" &lt;ian.hinder@aei.mpg.de&gt;<br><b>Cc: </b>"Einstein Toolkit Users" &lt;users@einsteintoolkit.org&gt;<br><b>Sent: </b>Wednesday, April 12, 2017 3:58:07 PM<br><b>Subject: </b>Re: [Users] units<br></div><br><div data-marker="__QUOTED_TEXT__"><div dir="ltr"><span style="font-size:12.800000190734863px">Hi Francisco,</span><br><div>To add a little more detail to what Ian and Erik already said, you might find it helpful to look at the Waveforms SimulationTools tutorial [<a href="http://simulationtools.org/examples/Waveforms.cdf" target="_blank">http://simulationtools.org/examples/Waveforms.cdf</a>] from&nbsp;the BBH gallery example page. This gives a short explanation of the units in that simulation along with demonstrations of how to interpret the data from the simulation in terms of those units.</div><br><div>Regards,</div><div>Barry</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 12 April 2017 at 21:21, Ian Hinder <span dir="ltr">&lt;<a href="mailto:ian.hinder@aei.mpg.de" target="_blank">ian.hinder@aei.mpg.de</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><br><div><span class=""><div>On 12 Apr 2017, at 22:12, Erik Schnetter &lt;<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>&gt; wrote:</div><br class="m_537190974548160429Apple-interchange-newline"><blockquote><div dir="ltr">Francisco<br><div>Without hydrodynamics, nothing fundamental sets the mass scale. If you assume that c = G = 1, then there is a rescaling freedom for M, the black hole mass. Apart from this, it is up to the designer of the parameter file what they call M. I think that the ADM mass is a good choice, but others might use a different convention.</div></div></blockquote><br></span><div>Hi Erik,</div><br><div>If you want the units of the simulation to be the ADM mass, then you need a mechanism to iterate the bare mass parameters of the punctures to make this the case.&nbsp; We don't have that in TwoPunctures; we iterate until we get the puncture masses equal to the targets.&nbsp; I think it is useful to compare two BBH systems where the BHs have the same total mass. This would be natural also if you were comparing with PN.&nbsp; I'm not sure it is so useful to compare systems with the same ADM mass; the ADM mass will consist of the BH masses, plus the energy content in the junk radiation, and the orbital energy.&nbsp; The junk radiation isn't very interesting, which is why I am rarely interested in M_ADM.</div><span class=""><br><blockquote><div dir="ltr"><div>For example, in &lt;<a href="http://einsteintoolkit.org/gallery/bbh/index.html" target="_blank">http://einsteintoolkit.org/gallery/bbh/index.html</a>&gt; I assume that M_ADM = 1, but I don't know whether this is taken from the initial conditions or whether this is measured after the initial junk radiation has left.</div></div></blockquote><br></span><div>No, M_ADM is not 1 for that parameter file.&nbsp; It's</div><br><div><div>initial-ADM-energy &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0.9899366929086094169</div><br><div>In the parameter file, we have</div></div><div><pre style="word-wrap:break-word;white-space:pre-wrap">q = 36.0/29.0               # Mass ratio: q = mp/mm &gt;= 1
M = 1.0                     # Total mass</pre></div><div><pre style="word-wrap:break-word;white-space:pre-wrap">mp = M * q/(1+q)           # Heavier, larger BH, AH1, SS 0
mm = M * 1/(1+q)           # Lighter, smaller BH, AH2, SS 1</pre></div><div><pre style="word-wrap:break-word;white-space:pre-wrap">TwoPunctures::target_M_plus             = $mp
TwoPunctures::target_M_minus            = $mm
</pre><br></div><div>so the sum of the puncture masses is 1 in the units of the simulation, so the units of the simulation are the sum of the puncture masses.</div></div><span class=""><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>--&nbsp;</div><div>Ian Hinder</div><div><a href="http://members.aei.mpg.de/ianhin" target="_blank">http://members.aei.mpg.de/ianhin</a></div></div></div></div></div>
</div>
<br></span></div><br>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>Users mailing list<br>Users@einsteintoolkit.org<br>http://lists.einsteintoolkit.org/mailman/listinfo/users<br></div></div></body></html>