Unequal-mass Binaries

To construct a irrotational binary configuration with unequal-mass stars, you can use the ``setup2iq'' subroutine, which is called by setting the three letter code in sph.init to ``2iq''.

The particle positions and other important parameters are read in from files called ``pos.sph'' and ``pos2.sph'', which you will need to create. We have supplied the code to make executables called ``b2pos'' and ``b2pos2'', which will write these files for you. Based on responses to two questions, it can convert binary output files, in the form ``out012.sph'' (or any other 3-digit number) or ``restart.sph'' into the ASCII pos.sph and pos2.sph files. For more on file I/O, please see Sec. 5.

Note that the parameter ``N'' used to initialize the run is essentially ignored, since the number of particles is set by the configuration read in from the files.

Additionally, the run will need you to create input files called ``axes1.sph'' and ``axes2.sph''. These files, in ASCII format, each consist of three numbers representing the amount by which the spherical model is deformed in each of the three principal axes (x: the center of mass separation vector, y: the orbital velocity, z: the angular momentum axis) for the first and second star. Typically, all three values should be $\sim 1$, and more often than not, $a_1>a_2>a_3$. The axis ratios should differ between the two stars, especially when the mass ratio diverges from unity. These files will be created for you by the testinput executable discussed above.