Detailed models of observed interacting galaxies suffer from the extended parameter space. Here, we present results from our code MINGA which couples an evolutionary optimization strategy (a genetic algorithm) with a fast N-body method. MINGA allows for an automatic search of the optimal region(s) in parameter space within a few hours to a few days of CPU time on a modern PC by investigating of the order of 10^5 models. We demonstrate its applicability by modelling the HI intensity and velocity maps of the interacting system M51 and NGC 5195. We get a good fit for the HI intensity map and we can reproduce the counter-rotation feature of the HI arm. Our result corroborates the results of Salo & Laurikainen (2000) who favour a model with multiple passages through M51's disk.

The main problems in modeling interacting galaxies are the extended parameter space and the fairly high CPU costs of self-consistent N-body simulations. Therefore, traditional modeling techniques suffer from either extreme CPU demands or trapping in local optima (or both). A very promising alternative approach are evolutionary algorithms which mimic natural adaptation in order to optimize the numerical models. One main advantage is their very weak dependence on starting points which makes them much less prone to trapping in local optima. We present a Genetic Algorithm (GA) coupled with a fast (but not self-consistent) restricted N-body solver. This combination allows us to identify interesting regions of parameter space within only a few CPU hours on a standard PC or a few CPU minutes on a parallel computer. Especially, we demonstrate here the ability of GA-based fitting procedures to analyse observational data automatically, provided the data are sufficiently accurate.