Context. More than 60 planets have been discovered so far in systems that harbour two stars, some of which have binary semi-major axes as small as 20 au. It is well known that the formation of planets in such systems is strongly influenced by the stellar components, since the protoplanetary disc and the particles within are exposed to the gravitational influence of the binary. However, the question on how self-gravitating protoplanetary bodies affect the evolution of a radiative, circumprimary disc is still open. Aims. We present our 2D hydrodynamical GPU-CPU code and study the interaction of several thousands of self-gravitating particles with a viscous and radiative circumprimary disc within a binary star system. To our knowledge this program is the only one at the moment that is capable to handle this many particles and to calculate their influence on each other and on the disc. Methods. We performed hydrodynamical simulations of a circumstellar disc assuming the binary system to be coplanar. Our gridbased staggered mesh code relies on ideas from ZEUS-2D, where we implemented the FARGO algorithm and an additional energy equation for the radiative cooling according to opacity tables. To treat particle motion we used a parallelised version of the precise Bulirsch - Stoer algorithm. Four models in total where computed taking into account (i) only N-body interaction, (ii) N-body and disc interaction, (iii) the influence of computational parameters (especially smoothing) on N-body interaction, and (iv) the influence of a quiet low-eccentricity disc while running model (ii). The impact velocities where measured at two different time intervals and were compared. Results. We show that the combination of disc- and N-body self-gravity can have a significant influence on the orbit evolution of roughly Moon sized protoplanets.

We studied the impact of the revisited values for the LSR circular velocity of the Milky Way (Reid et al. 2004) on the formation of the Magellanic Stream. The LSR circular velocity was varied within its observational uncertainties as a free parameter of the interaction between the Large (LMC) and the Small (SMC) Magellanic Clouds and the Galaxy. We have shown that the large-scale morphology and kinematics of the Magellanic Stream may be reproduced as tidal features, assuming the recent values of the proper motions of the Magellanic Clouds (Kallivayalil et al. 2006). Automated exploration of the entire parameter space for the interaction was performed to identify all parameter combinations that allow for modeling the Magellanic Stream. Satisfactory models exist for the dynamical mass of the Milky Way within a wide range of 0.6*10^12Msun to 3.0*10^12Msun and over the entire 1-sigma errors of the proper motions of the Clouds. However, the successful models share a common interaction scenario. The Magellanic Clouds are satellites of the Milky Way, and in all cases two close LMC-SMC encounters occurred within the last 4Gyr at t<-2.5Gyr and t approx. -150Myr, triggering the formation of the Stream and of the Magellanic Bridge, respectively. The latter encounter is encoded in the observed proper motions and inevitable in any model of the interaction. We conclude that the tidal origin of the Magellanic Stream implies the previously introduced LMC/SMC orbital history, unless the parameters of the interaction are revised substantially.

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval of ~20 Myr after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence.

The dark matter content of galaxies is usually determined from galaxies in dynamical equilibrium, mainly from rotationally supported galactic components. Such determinations restrict measurements to special regions in galaxies, e.g. the galactic plane(s), whereas other regions are not probed at all. Interacting galaxies offer an alternative, because extended tidal tails often probe outer or off-plane regions of galaxies. However, these systems are neither in dynamical equilibrium nor simple, because they are composed of two or more galaxies, by this increasing the associated parameter space.We present our genetic algorithm based modeling tool which allows to investigate the extended parameter space of interacting galaxies. From these studies, we derive the dynamical history of (well observed) galaxies. Among other parameters we constrain the dark matter content of the involved galaxies. We demonstrate the applicability of this strategy with examples ranging from stellar streams around theMilkyWay to extended tidal tails, from proto-typical binary galaxies (like M51 or the Antennae system) to small group of galaxies.

We present hydrodynamical models of the grand design spiral M51 (NGC 5194), and its interaction with its companion NGC 5195. Despite the simplicity of our models, our simulations capture the present day spiral structure of M51 remarkably well, and even reproduce details such as a kink along one spiral arm, and spiral arm bifurcations. We investigate the offset between the stellar and gaseous spiral arms, and find at most times (including the present day) there is no offset between the stars and gas to within our error bars. We also compare our simulations with recent observational analysis of M51. We compute the pattern speed versus radius, and like the observations, find no single global pattern speed. We also show that the spiral arms cannot be fitted well by logarithmic spirals. We interpret these findings as evidence that M51 does not exhibit a quasi-steady density wave, as would be predicted by density wave theory. The internal structure of M51 derives from the complicated and dynamical interaction with its companion, resulting in spiral arms showing considerable structure in the form of short-lived kinks and bifurcations. Rather than trying to model such galaxies in terms of global spiral modes with fixed pattern speeds, it is more realistic to start from a picture in which the spiral arms, while not being simple material arms, are the result of tidally induced kinematic density `waves' or density patterns, which wind up slowly over time.

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, depending on the shape of the potential, tides can also become fully compressive. When that is the case, they might trigger or strengthen the formation of galactic substructures (star clusters, tidal dwarf galaxies), instead of destroying them. We perform N-body simulations of interacting galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a galaxy merger, independently of the specific choice of parameterization. With a model tailored to the Antennae galaxies, we show that the distribution of compressive tides matches the locations and timescales of observed substructures. After extending our study to a broad range of parameters, we conclude that neither the importance of the compressive tides (~15% of the stellar mass) nor their duration (~ 10 Myr) are strongly affected by changes in the progenitors' configurations and orbits. Moreover, we show that individual clumps of matter can enter compressive regions several times in the course of a simulation. We speculate that this may spawn multiple star formation episodes in some star clusters, through e.g., enhanced gas retention.

We study the stability properties of isolated star forming dwarf galaxies which undergo dynamically driven starbursts induced by stellar feedback. Here we focus on the impact of the adopted ISM model, i.e. either a diffuse or a clumpy ISM. We apply a one-zone model extended for active dynamical evolution. We found two major types of repetitive star bursts: one set (type A) of quasi-periodic starbursts is related to the dynamical timescale of the galaxy. In that case, the star formation follows the variations of the gas density induced by decaying virial oscillations. The second set (type B) of starbursts is characterized by long quiescence periods given by the sum of the dynamical and the dissipative timescale: after a first burst, the inserted energy leads to a substantial expansion of the system, by this stopping any significant star formation activity. A next burst might occur, when the gas reaches high densities again, i.e. after the gas recollapsed and the energy injected by stellar feedback is dissipated. In case of a diffuse ISM model, type A bursts are the most common type due to the high efficiency of radiative cooling (no type B bursts are found). Bursts occur then mainly during an initial transitory phase. In case of a clumpy ISM model (i.e. dissipation by inelastic cloud-cloud collisions), the dissipative timescale is of the order of the dynamical time or longer. This allows for both, type A and type B bursts.

The M51 galaxy is one of best examples of grand design structure in galaxies. The most likely origin of the spiral arms in M51 is the ongoing interaction with NGC 5195. Here we report recent calculations which model the orbit of M51 and NGC 5195. These calculations, which for the first time focus on the gas dynamics, are able to reproduce even detailed features in M51, as well as demonstrate the nature of the spiral structure of M51.

Galaxy interactions are a common phenomenon in clusters of galaxies. Especially major mergers are of particular importance, because they can change the morphological type of galaxies. They have an impact on the mass function of galaxies and they trigger star formation - the main driver of the Galactic Matter Cycle. Therefore, we conducted a study of major mergers by means of a multi-phase ISM code. This code is based on a TREE-SPH-code combined with a sticky particle method allowing for star formation controlled by the properties of a multi-phase ISM. This is in contrast to the usually implemented Schmidt law depending mainly on the gas density. Previously, this code was used on isolated galaxies. Since our star formation recipe is not restricted to a special type of galaxy, it is interesting to apply it to interacting galaxies, too. Our study on major mergers includes a research of global properties of the interacting system, namely the star formation rate and the star formation efficiency, the evaporation and condensation rates, as well as the mass exchange of distinct components, namely stars, diffuse ISM, and clouds. Investigating these properties provides insight to interrelations between various physical processes. The results indicate that the star formation efficiency as well as the evaporation and condensation rates are influenced by the interaction.

Physics lectures always refer to the tides as a disruptive effect. However, tides can also be compressive. When the potential of two galaxies overlap, as happens during a merger, fully compressive tides can develop and have a strong impact on the dynamics of substructures such as star clusters or tidal dwarf galaxies. Using N-body simulations of a large set of mergers, we noticed the importance of these tidal modes at cluster scale. With a model of the Antennae galaxies, we conclude that the positions and timescales of these tidal modes match the actual distribution of young clusters. A detailed study of the statistics of the compressive tides shows a stunning correlation between this purely gravitational effect and the observed properties of the star clusters. In this contribution, we introduce the concept of compressive tide and show its relevance in the special case of the Antennae galaxies. We extend our conclusions to a broad range of parameters and discuss their implications on several critical points such as the infant mortality, multiple star formation epochs in clusters or the age distribution.

Starbursts and substantial variations in the star formation histories are a common phenomenon in galaxies. We study the stability properties of isolated star-forming dwarf galaxies with the aim of identifying starburst modes. The impact of the stellar birth function, the initial mass function (IMF), the stellar feedback and the interstellar medium (ISM) model are investigated. We apply a one-zone model for a star-gas system coupled by mass and energy transfer. Additionally, we extend the network for active dynamical evolution. This allows for a coupling between the dynamical state of the galaxy and its internal properties. While the influence of the dynamics on the total star formation rate is strong, the coupling of the internal properties (gas temperature) on the dynamics is rather limited, because radiative cooling keeps the gas temperature well below the virial temperature. Because of short cooling and feedback timescales, the star formation rate is close to the equilibrium star formation rates. Quasi-periodic starbursts occur, because star formation follows the variations in the gas density induced by decaying virial oscillations. This behaviour is quite insensitive to the nature and the details of the stellar birth description, viz. whether spontaneous or induced star formation is considered or the IMF is varied. A second type of burst is found as an instability operating when the cooling may drop at very low densities with increasing temperature. Bursts of star formation occur during transitory phases, when dynamical equilibrium is established. Then they are quasi-periodic on the dynamical timescale. Because of short heating and cooling timescales, the star formation rate follows the equilibrium star formation rate corresponding to the actual gas density.

The dwarf satellite galaxies in the Local Group are generally considered to be hosted in dark matter subhalos that survived the disruptive processes during infall onto their host halos. It has recently been argued that if the majority of satellites entered the Milky Way halo in a group rather than individually, this could explain the spatial and dynamical peculiarities of its satellite distribution. Such groups were identified as dwarf galaxy associations that are found in the nearby Universe. In this paper we address the question whether galaxies in such associations can be the progenitors of the Milky Way satellite galaxies. We find that the dwarf associations are much more extended than would be required to explain the disk-like distribution of the Milky Way and Andromeda satellite galaxies. We further identify a possible minor filamentary structure, perpendicular to the supergalactic plane, in which the dwarf associations are located, that might be related to the direction of infall of a progenitor galaxy of the Milky Way satellites, if they are of tidal origin.

We expand our pervious numerical study of the properties of the stellar velocity distribution within the disk of a two-armed spiral galaxy by considering spiral stellar density waves with different global Fourier amplitudes, C_2. We confirm our previous conclusion that the ratio \sigma_1:\sigma_2 of smallest versus largest principal axes of the stellar velocity ellipsoid becomes abnormally small near the outer edges of the stellar spiral arms. The extent to which the stellar velocity ellipsoid is elongated (as compared to the unperturbed value typical for the axisymmetric disk) increases with the strength of the spiral density wave. In particular, the C_2=0.06 spiral can decrease the unperturbed value of \sigma_1:\sigma_2 by 20%, while the C_2=0.13 spiral can decrease the unperturbed \sigma_1:\sigma_2 by a factor of 3. The abnormally small values of the \sigma_1:\sigma_2 ratio can potentially be used to track the position of {\it stellar} spiral density waves. The \sigma_{\phi\phi}:\sigma_{rr} ratio is characterized by a more complex behaviour and exhibits less definite minima near the outer edges of the spiral arms. We find that the epicycle approximation is violated near the spiral arms and cannot be used in spiral galaxies with C_2 >= 0.05-0.06 or in galaxies with the amplitude of the spiral stellar density wave (relative to the unperturbed background) of order 0.1 or greater.

Motivated by the closest major merger, the Antennae Galaxies (NGC4038/4039), we want to improve our genetic algorithm based modeling code Minga (Theis 1999). The aim is to reveal the major interaction and galaxy parameters, e.g. orbital information and halo properties of such an equal mass merger system. Together with the sophisticated search strategy of Minga, one needs fast and reliable models in order to investigate the high dimensional parameter space of this problem. Therefore we use a restricted N-body code which is based on the approach by Toomre & Toomre (1972), however with some refinements like consistent orbits of extended dark matter halos. Recently also dynamical friction was included to this code (Petsch 2007). While a good description for dynamical friction was found for mass ratios up to q = 1/3 (Petsch & Theis 2008), major merger systems were only imperfectly remodeled. Here we show recent improvements for a major merger system by including mass-loss and using NFW halos.

The tidal field of galaxies is known generally to be disruptive. However, in the case of galaxy mergers, a compressive mode of tidal wave may develop and last long enough to cocoon the formation of star clusters. Using an N-body simulation of the Antennae galaxies, we derive the positions of these compressive regions and the statistics of their duration. Excellent agreement between the spatial distribution of tides and observed young clusters is found, while the characteristic e-folding times of 10 to 30 Myrs derived for the tidal field compare well with cluster formation time-scales.

Recently, Kallivayalil et al. derived new values of the proper motion for the Large and Small Magellanic Clouds (LMC and SMC, respectively). The spatial velocities of both Clouds are unexpectedly higher than their previous values resulting from agreement between the available theoretical models of the Magellanic System and the observations of neutral hydrogen (HI) associated with the LMC and the SMC. Such proper motion estimates are likely to be at odds with the scenarios for creation of the large-scale structures in the Magellanic System suggested so far. We investigated this hypothesis for the pure tidal models, as they were the first ones devised to explain the evolution of the Magellanic System, and the tidal stripping is intrinsically involved in every model assuming the gravitational interaction. The parameter space for the Milky Way (MW)-LMC-SMC interaction was analyzed by a robust search algorithm (genetic algorithm) combined with a fast restricted N-body model of the interaction. Our method extended the known variety of evolutionary scenarios satisfying the observed kinematics and morphology of the Magellanic large-scale structures. Nevertheless, assuming the tidal interaction, no satisfactory reproduction of the HI data available for the Magellanic Clouds was achieved with the new proper motions. We conclude that for the proper motion data by Kallivayalil et al., within their 1-sigma errors, the dynamical evolution of the Magellanic System with the currently accepted total mass of the MW cannot be explained in the framework of pure tidal models. The optimal value for the western component of the LMC proper motion was found to be pm_w(LMC) > -1.3 mas/yr in case of tidal models. It corresponds to the reduction of the Kallivayalil et al. value for pm_w(LMC) by approx. 40% in its magnitude.

Modeling interacting galaxies to reproduce observed systems is still a challenge due to the extended parameter space (among other problems). Orbit and basic galaxy parameters can be tackled by fast simulation techniques like the restricted N-body method, applied in the fundamental work by Toomre & Toomre (1972). This approach allows today for the study of millions of models in a short time. One difficulty for the classical restricted N-body method is the missing orbital decay, not allowing for galaxy mergers. Here we present an extension of the restricted N-body method including dynamical friction. This treatment has been developed by a quantitative comparison with a set of self-consistent merger simulations. By varying the dynamical friction (formalism, strength and direction), we selected the best-fitting parameters for a set of more than 250000 simulations. We show that our treatment reliably reproduces the orbital decay and tidal features of merging disk galaxies for mass ratios up to q=1/3 between host and satellite. We implemented this technique into our genetic algorithm based modeling code MINGA and present first results.

In the framework of hierarchical structure formation ellipticals can form from merging of smaller disk galaxies. The nearby interacting 'Antennae' galaxy pair (NGC 4038/39) is one of the best-studied local systems of merging spirals, thus presenting us with an ideal laboratory for the study of galaxy evolution models. The Antennae are believed to be in a state prior to their final encounter with rapid subsequent merging, which puts them in the first position of the Toomre (1977) merger sequence. Here we present first numerical high-resolution, self-consistent, smoothed particle hydrodynamics (SPH) simulations of the Antennae system, including star formation and stellar feedback, and compare our results to VLA HI observations by Hibbard et al. (2001). We are able to obtain a close, but not yet perfect match to the observed morphology and kinematics of the system.

Gravitational tides are widely understood to strip and destroy galactic substructures. In the course of a galaxy merger, however, transient totally compressive tides may develop and prevent star forming regions from dissolving, after they condensed to form clusters of stars. We study the statistics of such compressive modes in an N-body model of the galaxy merger NGC 4038/39 (the Antennae) and show that ~15% of the disc material undergoes compressive tides at pericentre. The spatial distribution of observed young clusters in the overlap and nuclear regions of the Antennae matches surprisingly well the location of compressive tides obtained from simulation data. Furthermore, the statistics of time intervals spent by individual particles embedded in a compressive tide yields a log-normal distribution of characteristic time ~10 Myr, comparable to star cluster formation timescales. We argue that this generic process is operative in galaxy mergers at all redshifts and possibly enhances the formation of star clusters. We show with a model calculation that this process will prevent the dissolution of a star cluster during the formation phase, even for a star formation efficiency as low as ~10%. The transient nature of compressive tides implies that clusters may dissolve rapidly once the tidal field switches to the usual disruptive mode.

We present numerical studies of the properties of the stellar velocity distribution in galactic disks which have developed a saturated, two-armed spiral structure. In previous papers we used the Boltzmann moment equations (BME) up to second order for our studies of the velocity structure in self-gravitating stellar disks. A key assumption of our BME approach is the zero-heat flux approximation, i.e. the neglection of third order velocity terms. We tested this assumption by performing test particle simulations for stars in a disk galaxy subject to a rotating spiral perturbation. As a result we corroborated qualitatively the complex velocity structure found in the BME approach. It turned out that an equilibrium configuration in velocity space is only slowly established on a typical timescale of 5 Gyrs or more. Since many dynamical processes in galaxies (like the growth of spirals or bars)act on shorter timescales, pure equilibrium models might not be fully appropriate for a detailed comparison with observations like the local Galactic velocity distribution. Third order velocity moments were typically small and uncorrelated over almost all of the disk with the exception of the 4:1 resonance region (UHR). Near the UHR (normalized) fourth and fifth order velocity moments are still of the same order as the second and third order terms. Thus, at the UHR higher order terms are not negligible.

We present a numerical study of the properties of the stellar velocity distribution in stellar discs which have developed a saturated, two-armed spiral structure. We follow the growth of the spiral structure deeply into the non-linear regime by solving the Boltzmann moment equations up to second order. By adopting the thin-disc approximation, we restrict our study of the stellar velocity distribution to the plane of the stellar disc. We find that the outer (convex) edges of stellar spiral arms are characterized by peculiar properties of the stellar velocity ellipsoids, which make them distinct from most other galactic regions. In particular, the ratio \sigma_1:\sigma_2 of the smallest versus largest principal axes of the stellar velocity ellipsoid can become abnormally small (as compared to the rest of the disc) near the outer edges of spiral arms. Moreover, the epicycle approximation fails to reproduce the ratio of the tangential versus radial velocity dispersions in these regions. These peculiar properties of the stellar velocity distribution are caused by large-scale non-circular motions of stars, which in turn are triggered by the non-axisymmetric gravitational field of stellar spiral arms. The magnitude of the vertex deviation appears to correlate globally with the amplitude of the spiral stellar density perturbations. However, locally there is no simple correlation between the vertex deviation and the density perturbations. (Abstract abridged).

Context: Dwarf galaxies can arise from self-gravitating structures emerging from tidal tails. What fraction of the known dwarf galaxies in the Local Universe can have this origin is still a matter of debate. Aims: In our effort to understand the origin and evolution of tidal dwarf galaxies and their correspondence with local objects, the first step is to understand how these galaxies (which are supposed to have a limited amount of dark matter) react to the feedback of the ongoing star formation. Methods: We make use of 2-D chemodynamical calculations in order to study the early evolution of isolated, dark matter-free dwarf galaxies. We present models in which feedback parameters are varied. We also compare the results with dark matter-dominated dwarf galaxy models. Results: All the considered models show that the star formation proceeds for more than 300 Myr, therefore dwarf galaxies without large dark matter halos are not necessarily quickly destroyed. The chemical evolution of these objects is consistent with the main chemical properties of the dSphs of the Local Group. Models with large dark matter halos show results consistent with models free of dark matter, indicating that the distribution of gas is more important than the depth of the potential well in determining the global behaviour of dSph-sized dwarf galaxies.

Our aim is to study the evolution of tidal dwarf galaxies. The first step is to understand whether a model galaxy without Dark Matter can sustain the feedback of the ongoing star formation. We present tests of the evolution of models in which star formation efficiency, temperature threshold, initial distribution of gas and infall are varied. We conclude that it is feasible to keep a fraction of gas bound for several hundreds of Myr and that the development of galactic winds does not necessarily stop continuous star formation.

We present the Boltzmann moment equation approach for the dynamics of stars (BEADS-2D), which is a finite-difference Eulerian numerical code designed for the modelling of anisotropic and non-axisymmetric flat stellar disks. The BEADS-2D code solves the Boltzmann moment equations up to second order in the thin-disk approximation. This allows us to obtain the anisotropy of the velocity ellipsoid and the vertex deviation in the plane of the disk. We study the time-dependent evolution of exponential stellar disks in the linear regime and beyond. The disks are initially characterized by different values of the Toomre parameter Q_s and are embedded in a dark matter halo, yielding a rotation curve composed of a rigid central part and a flat outer region. Starting from a near equilibrium state, several unstable modes develop in the disk. In the early linear phase, the very centre and the large scales are characterized by growing one-armed and bisymmetric positive density perturbations, respectively. This is in agreement with expectations from the swing amplification mechanism of short-wavelength trailing disturbances, propagating through the disk centre. In the late linear phase, the overall appearance is dominated by a two-armed spiral structure localized within the outer Lindblad resonance (OLR). During the non-linear evolutionary phase, radial mass redistribution due to the gravitational torques of spiral arms produces an outflow of mass, which forms a ring at the OLR, and an inflow of mass, which forms a transient central bar. This process of mass redistribution is self-regulatory and it terminates when spiral arms diminish due to a shortage of matter. Finally, a compact central disk and a diffuse ring at the OLR are formed (the abstract is abridged).

We performed an extended analysis of the parameter space for the interaction of the Magellanic System with the Milky Way (MW). The varied parameters cover the phase space parameters, the masses, the structure, and the orientation of both Magellanic Clouds, as well as the flattening of the dark matter halo of the MW. The analysis was done by a specially adopted optimization code searching for the best match between numerical models and the detailed HI map of the Magellanic System by Bruens et al. (2005). The applied search algorithm is a genetic algorithm combined with a code based on the fast, but approximative restricted N-body method. By this, we were able to analyze more than 10^6 models, which makes this study one of the most extended ones for the Magellanic System. Here we focus on the flattening q of the axially symmetric MW dark matter halo potential, that is studied within the range 0.74<=q<=1.20. We show that creation of a trailing tail (Magellanic Stream) and a leading stream (Leading Arm) is quite a common feature of the Magellanic System-MW interaction, and such structures were modeled across the entire range of halo flattening values. However, important differences exist between the models, concerning density distribution and kinematics of HI, and also the dynamical evolution of the Magellanic System. Detailed analysis of the overall agreement between modeled and observed distribution of neutral hydrogen shows that the models assuming an oblate (q<1.0) dark matter halo of the Galaxy allow for better satisfaction of HI observations than models with other halo configurations.

We present a new particle code for modelling the evolution of galaxies. The code is based on a multi-phase description for the interstellar medium (ISM). We included star formation (SF), stellar feedback by massive stars and planetary nebulae, phase transitions and interactions between gas clouds and ambient diffuse gas, namely condensation, evaporation, drag and energy dissipation. The latter is realised by radiative cooling and inelastic cloud-cloud collisions. We present new schemes for SF and stellar feedback. They include a consistent calculation of the star formation efficiency (SFE) based on ISM properties as well as a detailed redistribution of the feedback energy into the different ISM phases. As a first test example we show a model of the evolution of a present day Milky-Way-type galaxy. Though the model exhibits a quasi-stationary behaviour in global properties like mass fractions or surface densities, the evolution of the ISM is locally strongly variable depending on the local SF and stellar feedback. We start only with two distinct phases, but a three-phase ISM is formed soon consisting of cold molecular clouds, a warm gas disk and a hot gaseous halo. Hot gas is also found in bubbles in the disk accompanied by type II supernovae explosions. The star formation rate (SFR) is ~1.6 M_sun/year on average decreasing slowly with time due to gas consumption. In order to maintain a constant SFR gas replenishment, e.g. by infall, of the order 1 M_sun/year is required. Our model is in fair agreement with Kennicutt's (1998) SF law including the cut-off at ~10 M_sun/pc^2. Models with a constant SFE, i.e. no feedback on the SF, fail to reproduce Kennicutt's law.

To investigate the stability properties of polar disks we performed two-dimensional hydrodynamical simulations for flat polytropic gaseous self-gravitating disks which were perturbed by a central S0-like component. Our disk was constructed to resemble that of the proto-typical galaxy NGC 4650A. This central perturbation induces initially a stationary two-armed tightly-wound leading spiral in the polar disk. For a hot disk (Toomre parameter Q>1.7), the structure does not change over the simulation time of 4.5 Gyr. In case of colder disks, the self-gravity of the spiral becomes dominant, it decouples from the central perturbation and grows, until reaching a saturation stage in which an open trailing spiral is formed, rather similar to that observed in NGC4650A. The timescale for developing non-linear structures is 1-2 Gyr; saturation is reached within 2-3 Gyr. The main parameter controlling the structure formation is the Toomre parameter. The results are surprisingly insensitive to the properties of the central component. If the polar disk is much less massive than that in NGC4650A, it forms a weaker tightly-wound spiral, similar to that seen in dust absorption in the dust disk of NGC2787. Our results are derived for a polytropic equation of state, but appear to be generic as the adiabatic exponent is varied between \gamma = 1 (isothermal) and \gamma = 2 (very stiff).

We present a parameter study of the possibility of tidally triggered disk instability. Using a restricted N-body model which allows for a survey of an extended parameter space, we show that a passing dwarf star with a mass between 0.1 and 1 M_sun can probably induce gravitational instabilities in the pre-planetary solar disk for prograde passages with minimum separations below 80-170 AU for isothermal or adiabatic disks. Inclined and retrograde encounters lead to similar results but require slightly closer passages. Such encounter distances are quite likely in young moderately massive star clusters (Scally & Clarke 2001; Bonnell et al. 2001). The induced gravitational instabilities may lead to enhanced planetesimal formation in the outer regions of the protoplanetary disk, and could therefore be relevant for the existence of Uranus and Neptune, whose formation timescale of about 100 Myr (Wuchterl, Guillot & Lissauer 2000) is inconsistent with the disk lifetimes of about a few Myr according to observational data by Haisch, Lada & Lada (2001). The relatively small gas/solid ratio in Uranus and Neptune can be matched if the perturbing fly-by occurred after early gas depletion of the solar system, i.e. when the solar system was older than about 5 Myr. We also confirm earlier results by Heller (1993) that the observed 7 degree tilt of the solar equatorial plane relative to the ecliptic plane could be the consequence of such a close encounter.

We show that the shape of the observed distribution of Milky Way (MW) satellites is inconsistent with being drawn from a cosmological sub-structure population with a confidence of 99.5 per cent. Most of the MW satellites therefore cannot be related to dark-matter dominated satellites.

Dynamically cold components are well known to destabilize hotter, even much more massive components. In this paper we studied the dynamical influence of a cold dust component on the gaseous phase in the central regions of galactic disks. We performed two-dimensional hydrodynamical simulations for flat multi-component disks embedded in a combined static stellar and dark matter potential. The pressure-free dust component is coupled to the gas by a drag force depending on their velocity difference. It turned out that the most unstable regions are those with either a low or near to minimum Toomre parameter or with rigid rotation, i.e. the central area. In that regions the dust-free disks become most unstable for high azimuthal modes (m~8), whereas in dusty disks all modes have a similar amplitude resulting in a patchy appearance. The structures in the dust have a larger contrast between arm and inter-arm regions than those of the gas. The dust peaks are frequently correlated with peaks of the gas distribution, but they do not necessarily coincide with them. Therefore, a large scatter in the dust-to-gas ratios is expected. The appearance of the dust is more cellular (i.e. sometimes connecting different spiral features), whereas the gas is organized in a multi-armed spiral structure. An admixture of 2% dust destabilizes gaseous disks substantially, whereas dust-to-gas ratios below 1% have no influence on the evolution of the gaseous disk. For a high dust-to-gas ratio of 10% the instabilities reach a saturation level already after 30 Myr.

We studied the dynamical influence of a dust component on the gaseous phase in central regions of galactic disks. Therefore, we performed two-dimensional hydrodynamical simulations for flat multi-component disks embedded in a stellar and dark matter potential. The pressure-free dust component is coupled to the gas by a drag force depending on their velocity difference. It turned out that the most unstable regions are those with either a low or near to minimum Toomre parameter or with rigid rotation, i.e. the central area. In that regions the dust-free disks become most unstable for a small range of high azimuthal modes (m \sim 8), whereas in dusty disks all modes have similar amplitudes resulting in a patchy appearance. The structures in the dust have a larger contrast between arm and inter-arm regions than those of the gas. The dust peaks are frequently correlated with peaks of the gas distribution, but they do not necessarily coincide with them. This leads to a large scatter in the dust-to-gas ratios. The appearance of the dust is more cellular (i.e. sometimes connecting different spiral features), whereas the gas is organized in a multi-armed spiral structure. We found that an admixture of 2% dust (relative to the mass of the gas) destabilizes gaseous disks substantially, whereas dust-to-gas ratios below 1% have no influence on the evolution of the gaseous disk. For a high dust-to-gas ratio of 10% the instabilities reach the saturation level already after 30 Myr.

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.

We present a modified TREESPH code to model galaxies in 3d. The model includes a multi-phase description of the interstellar medium which combines two numerical techniques. A diffuse warm/hot gas phase is modelled by SPH while a sticky particle scheme is used to represent a cloudy medium. Interaction processes, such as star formation and feedback, cooling and mixing by condensation and evaporation, are taken into account. Here we apply our model to the evolution of a Milky Way type galaxy. After an initial stage, a quasi-equilibrium state is reached. It is characterised by a star formation rate of ~1 M_sun/year. Condensation and evaporation rates are in balance at 0.1-1 M_sun/year.

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.

We present a 3d code for the dynamical evolution of a multi-phase interstellar medium (ISM) coupled to stars via star formation (SF) and feedback processes. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation processes. The cycle of matter is completed by SF and feedback by SNe and PNe. A SF scheme based on a variable SF efficiency as proposed by Elmegreen & Efremov (1997) is presented. For a Milky Way type galaxy we get a SF rate of ~1 M_sun/yr with an average SF efficiency of ~5%.

Chemo-dynamical models have been introduced in the late eighties and are a generally accepted tool for understanding galaxy evolution. They have been successfully applied to one-dimensional problems, e.g. the evolution of non-rotating galaxies, and two-dimensional problems, e.g. the evolution of disk galaxies. Recently, also three-dimensional chemo-dynamical models have become available. In these models the dynamics of different components, i.e. dark matter, stars and a multi-phase interstellar medium, are treated in a self-consistent way and several processes allow for an exchange of matter, energy and momentum between the components or different gas phases. Some results of chemo-dynamical models and their comparison with observations of chemical abundances or star formation histories will be reviewed.

In this paper we present some test results of our newly developed Multi-Phase Chemo-Dynamical Smoothed Particle Hydrodynamics (MP- CD-SPH) code for galaxy evolution. At first, we present a test of the ``pure'' hydro SPH part of the code. Then we describe and test the multi-phase description of the gaseous components of the interstellar matter. In this second part we also compare our condensation and evaporation description with the results of a previous 2d multi-phase hydrodynamic mesh code.

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless $N$-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example here we present the dynamical, chemical and photometric evolution of a star forming dwarf galaxy with a total baryonic mass of $2 \times 10^9 {\rm M}_\odot$.

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. This code follows the evolution of different galactic components like stars, dark matter and different components of the interstellar medium (ISM), i.e. a diffuse gaseous phase and the molecular clouds. Stars and dark matter are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse gas and a sticky particle scheme for the molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example we show the dynamical and chemical evolution of a star forming dwarf galaxy with a total baryonic mass of $2 \cdot 10^9 M_\odot$. After a moderate collapse phase the stars and the molecular clouds follow an exponential radial distribution, whereas the diffuse gas shows a central depression as a result of stellar feedback. The metallicities of the galactic components behave quite differently with respect to their temporal evolution as well as their radial distribution. Especially, the ISM is at no stage well mixed.

We present a new particle based code with a multi-phase description of the ISM implemented in order to follow the chemo-dynamical evolution of galaxies. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): Exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation. Based on time scales we show that in Milky-Way-like galaxies the drag force is for molecular clouds only important, if their relative velocities exceed 100 km/s. For the mass exchange we find that clouds evaporate only if the temperature of the ambient gas is higher than one million Kelvin. At lower temperatures condensation takes place at time scales of the order of 1-10 Gyr.

Most formation scenarios of globular clusters assume a molecular cloud as the progenitor of the stellar system. However, it is still unclear, how this cloud is transformed into a star cluster, i.e. how the destructive processes related to gas removal or low star formation efficiency can be avoided. Here a scheme of supernova (SN) induced cluster formation is studied. According to this scenario an expanding SN shell accumulates the mass of the cloud. This is accompanied by fragmentation resulting in star formation in the shell. Provided the stellar shell expands sufficiently slow, its self-gravity stops the expansion and the shell recollapses, by this forming a stellar system. I present N-body simulations of collapsing shells which move in a galactic potential on circular and elliptic orbits. It is shown that typical shells (10^5 Msun, 30 pc) evolve to twin clusters over a large range of galactocentric distances. Outside this range single stellar systems are formed, whereas at small galactocentric distances the shells are tidally disrupted. In that case many small fragments formed during the collapse survive as single bound entities. About 1/3 of the twin cluster systems formed on circular orbits merge within 400 Myr. On elliptic orbits the merger rate reduces to less than 4%. Thus, there could be a significant number of twin clusters even in our Galaxy, which, however, might be undetected as twins due to a large phase shift on their common orbit.

We examine the merger stage of IRAS 03158+4227, one of the most luminous ULIRGs from the IRAS 2 Jy sample. Deep optical images are presented along with high-resolution NIR images and optical low-resolution spectra. IRAS 03158+4227 is confirmed as a component of an equal-luminosity binary galaxy with a nuclei separation of 47 / (h/75) kpc. A long lopsided tail emanating from the companion, which harbours an active nucleus, is the most prominent peculiar feature of the binary. The results of numerical simulations admit the interpretation of this structure as a product of the tidal interaction between the two components. If the infrared-activity of IRAS 03158+4227 is also dynamically triggered by this process, this would implicate that the ULIRG phenomenon is not restricted to the final stage of a binary merger. Alternatively, the system may be a multiple merger where the partner(s) has/have already coalesced.

We performed two-dimensional non-linear hydrodynamical simulations of two-component gravitating disks aimed at studying stability properties of these systems. In agreement with previous analytical and numerical simulations, we find that the cold gas component strongly affects the growth rates of the unstable global spiral modes in the disk. Already a five percent admixture of cold gas increases approximately two-fold the growth rate of the most unstable global mode while a twenty percent gas admixture enhances the modal growth rate by a factor of four. The local stability properties of a two-component disk coupled by self-gravity are governed by a stability criterion similar to Toomre's Q-parameter derived for one-component systems. Using numerical simulations, we analyse the applicability of a two-component local stability criterion for the analysis of the stability properties of global modes. The comparison of non-linear simulations with the analytical stability criterion shows that the two-component disks can be globally unstable while being stable to the local perturbations. The minimum value of the local stability criterion provides, however, a rough estimate of the global stability properties of two-component systems. We find that two-component systems with a content of cold gas of ten percent or less are globally stable, if the minimum value of the stability parameter exceeds about 2.5.

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.

The formation of globular clusters is still an unsolved problem. Though most scenarios assume a massive molecular cloud as the progenitor, it is unclear, how the cloud is transformed into a star cluster. Here a scheme of supernova (SN) induced cluster formation is investigated. In this scenario the expanding SN shell accumulates the mass of the cloud. This is accompanied by fragmentation resulting in star formation in the shell. If this stellar shell expands sufficiently slow, its self-gravity leads to a recollapsing shell, by this forming one or several stellar clusters. In this paper N-body simulations of collapsing shells moving on circular orbits in a galactic potential are presented. It is shown that typical shells (10^5 Msun, 30 pc) evolve to twin clusters in the galactocentric distance range between 3 and 11 kpc. Their masses show a strong radial trend: on orbits inside 5 kpc both clusters have almost equal mass. Outside 5 kpc the more massive twin cluster contains about 55% of the shell's mass, whereas the mass of the smaller decreases linearily to 15% at 11 kpc. Outside 11 kpc the collapsing shells end up in a single cluster. Inside 3 kpc the shells are tidally disrupted and only fragments substantially less massive than the initial shell survive.

(abridged) We combined several N-body methods in order to investigate the interaction scenario between NGC 4449 and DDO 125, a close companion in projected space. In a first step fast restricted N-body models are used to confine a region in parameter space reproducing the main observational features. In a second step a genetic algorithm is applied for a uniqueness test of our preferred parameter set. We show that our genetic algorithm reliably recovers orbital parameters, provided that the data are sufficiently accurate, i.e. all the key features are included. In the third step the results of the restricted N-body models are compared with self-consistent N-body simulations. In the case of NGC 4449, the applicability of the simple restricted N-body calculations is demonstrated. Additionally, it is shown that the HI gas can be modeled here by a purely stellar dynamical approach. In a series of simulations, we demonstrate that the observed features of the extended HI disc can be explained by a gravitational interaction between NGC 4449 and DDO 125. According to these calculations the closest approach between both galaxies happened $\sim 4-6 \cdot 10^8$ yr ago at a minimum distance of $\sim 25$ kpc on a parabolic or slightly elliptic orbit. In the case of an encounter scenario, the dynamical mass of DDO 125 should not be smaller than 10% of NGC 4449's mass. Before the encounter, the observed HI gas was arranged in a disc with a radius of 35-40 kpc around the center of NGC 4449. It had the same orientation as the central ellipsoidal HI structure. The origin of this disc is still unclear, but it might have been caused by a previous interaction.

Though it is generally assumed that massive molecular clouds are the progenitors of globular clusters, their detailed formation mechanism is still unclear. Standard scenarios based on the collapse of a smooth matter distribution suffer from strong requirements with respect to cluster formation time scale, binding energy and star formation efficiency. An alternative model assuming cluster formation due to the recollapse of a supernova-induced, fragmented shell can relax these difficulties. In this paper the final collapse stages of the different scenarios are compared by N-body simulations for shells and spheres. It is shown that fragmentation is much more pronounced for shells. Taking a galactic tidal field into account shells preferably form twin (or multiple) systems, whereas spheres end up as single clusters. The twins are characterized by identical metallicities, and stellar mass functions; some of them show counter-rotating cores. Their orbital evolution can result in both, a final merger or well separated twins sharing a common galactic orbit.

We present a detailed study of the cluster pair SL 353 & SL 349. This candidate binary cluster is located at the northwestern rim of the LMC bar. Based on photometric data we find that both clusters are coeval with an age of 550 +/- 100 Myr. We use the Ca II triplet in the spectra of individual red giants to derive radial velocities. Both components of the binary cluster candidate show very similar mean velocities (approx. 274 +/- 10 km/sec for SL 349 and approx. 279 +/- 4 km/sec for SL 353) while the field stars show lower velocities (approx. 240 +/- 19 km/sec). These findings suggest a common origin of the two clusters from the same GMC. In this sense the cluster pair may constitute a true binary cluster. We furthermore investigate the stellar densities in and around the star clusters and compare them with isopleths created from artificial, interacting as well as non-interacting star clusters. Gravitational interaction leads to a distortion which can also be found in the observed pair.

Modeling of interacting galaxies suffers from an extended parameter space prohibiting traditional grid based search strategies. As an alternative approach a combination of a Genetic Algorithm (GA) with fast restricted N-body simulations can be applied. A typical fit takes about 3-6 CPU-hours on a PentiumII processor. Here we present a parallel implementation of our GA which reduces the CPU-requirement of a parameter determination to a few minutes on 100 nodes of a CRAY T3E.

(abridged) Several N-body methods were combined in order to develop a method for the determination of the parameters of interacting galaxies. This method has been applied to the HI distribution of NGC4449. In a first step the fast restricted N-body models were used to confine a region in parameter space reproducing the main observational features. In a second step a genetic algorithm is introduced which allows for an automated search in parameter space as well as for a uniqueness test of special parameter sets. Using artificial data I show that the genetic algorithm reliably recovers orbital parameters. In the third step the results of the restricted N-body models are compared with detailed self-consistent N-body simulations. In the case of NGC4449 the applicability of the simple restricted N-body calculations is demonstrated. NGC4449 is an active star-forming dwarf galaxy of Magellanic type. From radio observations van Woerden et al. (1975) found an extended HI-halo around NGC4449 which is at least a factor of 10 larger than the optical diameter. More recently, Bajaja et al. (1994) and Hunter et al. (1998) discerned details in the HI-halo: a disc-like feature around the center of NGC4449 and a lopsided arm structure. In a series of simulations it is demonstrated that the main features can be obtained by a gravitational interaction between NGC4449 and DDO125, another nearby dwarf galaxy. According to these calculations the closest approach between both galaxies happened 350 Myr ago at a minimum distance of 25 kpc on a nearly parabolic orbit. In case of the encounter scenario, the dynamical mass of DDO125 should not be smaller than 10% of NGC4449's mass. The origin of this disc is still unclear, but it might be caused by a previous interaction.

Using 2D nonlinear simulations, we study the generation and nonlinear evolution of spiral structure in a star-forming multi-component gravitating disk. We confirm in agreement with previous studies the destabilizing role of a cold gaseous component and extend this conclusion for multi-component star-forming disks exchanging mass and momentum between its components. We show that the spiral structure growing on a non-stationary multi-phase background reaches its saturation in a similar manner like the one-component disks. The spiral structure survives even if most of the gas is transformed into stellar remnants of larger velocity dispersion.

A database on shape evolution of direct N-body models formed out of cold, dissipationless collapse is generated using GRAPE and HARP special purpose computers. Such models are important to understand the formation of elliptical galaxies. Three dynamically distinct phases of shape evolution were found, first a fast dynamical collapse which gives rise to the radial orbit instability (ROI) and generates at its end the maximal triaxiality of the system. Subsequently, two phases of violent and two-body shape relaxation occur, which drive the system first towards axisymmetry, finally to spherical symmetry (the final state, however, is still much more concentrated than the initial model). In a sequence of models the influence of numerical and physical parameters, like particle number, softening, initial virial ratio, timestep choice, different N-body codes, are examined. We find that an improper combination of softening and particle number can produce erroneous results. Selected models were evolved on the secular timescale until they became spherically symmetric again. The secular shape relaxation time scale is shown to agree very well with the two-body relaxation time, if softening is properly taken into account for the latter. Finally, we argue, that the intermediate phase of violent shape relaxation after collapse is induced by strong core oscillations in the centre, which cause potential fluctuations, dampening out the triaxiality.

The network of interactions between hot gas, cool clouds, massive stars, and stellar remnants used in the chemodynamical modeling of the interstellar medium is investigated for the types of its solutions. In a physically consistent formulation for the energy transfer during condensation, oscillations due to a cyclic switching between condensation and evaporation never occur. A closed-box system evolves in a hierarchy of equilibria: thermal balance in the cloud gas, star-formation self-regulated due to heating of the clouds by massive stars, and the balance of condensation and evaporation. Except for an initial transitory phase, the evolution of the metallicity follows that of the Simple Model quite closely. For galaxies with a high initial density, or if the condensation rate is low, the metals produced by the stars may remain stored in the hot gas phase even in evolved systems with low gas fractions.

Conditions for the fragmentation of expanding shells due to gravitational instability are discussed. The self-similar analytical solution is compared with the results of 3-dimensional computer simulations for the expansion into homogeneous media. For realistic galactic disks we show that the formation of fragments is influenced by the amount of energy supply from the final number of young stars in an OB association, the value of the sound speed, the stratification and density of the ambient medium, the galactic differential rotation and the gravitational force perpendicular to the galactic plane. The typical size of gravitationally unstable shells is 1 kpc for an ambient gas density n=1 cm^-3. In thick disk galaxies the fragmentation occurs in nearly the whole shell while in thin disks it is restricted to the galactic equator. Unstable fragments may become molecular and trigger the formation of molecular clouds, and finally new star formation. We conclude that in dwarf galaxies the star formation may propagate in all directions turning the system into a star- burst. Contrary to that, the star formation in spiral galaxies propagates only in some directions in a thin strip near the symmetry plane, basically at the tips of the expanding shell.