Second, the orbits used in Sect. 4 were all hyperbolic. In this section elliptical orbits will be considered as well. The elliptical orbits used will be fairly eccentric (e > 0.5) and will be chosen such that the discs are not severely damaged at pericentre passages.
Thus, the simulations to be described below constitute a greater challenge to the genetic algorithm. However, there is also an additional tool yet to be used, namely the radial velocity field. The use of the velocity field may be particularly useful for e.g. highly inclined systems for which the position data may yield ambiguous information.
The results from four runs with inclined discs are shown in Table 4. In these runs, both position and velocity data were used by the GA. As is evident from the table, the algorithm was able to find orbital parameters close to the correct ones in all cases. Note that the two first pairs shown in Table 4 correspond to elliptical orbits.
In order to test the importance of the velocity data, Run 7 was repeated using only position data, and the result from this run is shown in Table 4 as well. Clearly, the use of the velocity data improved the results, even though acceptable orbital parameters were found even in the case in which only position data were used.
For real systems, the inclinations and position angles are rarely known with high accuracy. In order to take this into account, the simulation program was modified to include the four angles and in the set of unknowns. Since the values of these angles can at least be estimated, they were allowed a range of variation of only degrees, in steps of 1 degree, thus increasing the size of the search space by a factor 214 = 194,481. Under these conditions, Run 7 was repeated once more and the results, including the inclinations and position angles found by the GA, are shown in the seventh row of Table 4. Due to the large increase in the size of the search space, the GA was allowed to run for 200 generations in this case, and, as can be seen from the table, acceptable orbital parameters were found in this case as well.
The observations corresponding to Runs 6 and 8 are shown in Fig. 4. Of particular interest is the fact that the algorithm is able to cope with the difficult case in which one of the discs is edge-on. Thus, the results reported in this section show that the GA method can successfully be applied to realistic interaction geometries.
Copyright The European Southern Observatory (ESO)