Proceedings of the Workshop
"The Magellanic Clouds and Other Dwarf Galaxies"
of the Bonn/Bochum-Graduiertenkolleg

The Outer Halo of NGC 4449

Christian Theis1 and Sven Kohle2

1Institut für Theoretische Physik und Astrophysik, Universität Kiel,
Olshausenstr. 40, 24098 Kiel, Germany
2Radioastronomisches Institut der Universität Bonn,
Auf dem Hügel 71, 53121 Bonn, Germany

Received 20th March 1998
Abstract. NGC 4449 is an active star-forming dwarf galaxy of Magellanic type. From radio observations van Woerden et al. (1975) found an extended H I-Halo around NGC 4449 which is at least a factor of 10 larger than the optical diameter D25≅5.6 kpc. More recently, Bajaja et al. (1994) and Hunter et al. (1998) discerned details in the H I-halo: a disk-like feature around the center of NGC 4449 and a lopsided arm structure. In a series of simulations we demonstrated that the main features can be obtained by a gravitational interaction between NGC 4449 and another dwarf galaxy, DDO 125. According to these calculations the closest approach between both galaxies happened ∼3·108 yr ago at a minimum distance of 30-35 kpc on a nearly parabolic orbit. In that case, the dynamical mass of DDO 125 should not be smaller than 10% of NGC 4449. Before the encounter the observed H I was organized in a disk with a radius of 25-30 kpc around the center of NGC 4449 which had an inclination angle of 55° to 70°. The origin of this disk is still not clear, but it might have been a previous additional interaction.

1. Introduction

Comparing the optical images of NGC 4449 and the Large Magellanic Cloud, there are many similarities like the size or the assumed mass, though NGC 4449 has a more active star formation. Therefore, NGC 4449 has widely been used as a reference galaxy for Magellanic type irregulars. On the other hand, radio observations performed by van Woerden et al. (1975) exhibited a large H I-structure exceeding the Holmberg diameter of 11 kpc (we will assume a distance of 3.9 Mpc in this paper). Recent observations by Bajaja et al. (1994) with the Effelsberg telescope and, more detailed results by Hunter et al. (1998) with the VLA revealed a complex structure of several components (Fig. 1): An elongated ellipse of H I gas centered on NGC 4449 has a total mass of 1.1·109 Msun. This gas shows a rigid rotation inside 11 kpc reaching a level of 97 km s-1 in the deprojected speed. Outside 11 kpc the rotation curve is flat. Bajaja et al. (1994) derived a dynamical mass of 7.0·1010 Msun inside 32 kpc, and a H I gas mass fraction of 3.3%. Furthermore, the H I gas inside the optical part of NGC 4449 is counterrotating with respect to the outer regions.

South of the galactic center a streamlike structure of 25 kpc emanates which abruptly splits into two parts: A small spur (5 kpc) which points towards DDO 125, another close irregular galaxy, and a long extended stream. The latter first points 24 kpc to the north, then 49 kpc north-east and, finally, 27 kpc to the south-east, by this covering 180° around NGC 4449. Remarkable are the long straight sections and the abrupt changes of direction. The total mass of the streamers measured by the VLA observations is about 1·109 Msun. Hunter et al. (1998) pointed out that by comparison with Effelsberg data there is an additional diffuse H I mass not detected in the VLA observations which amounts to two-thirds of the H I in the extended structures.

In order to explain the outer distribution of H I there are two classes of explanations: The first one assumes that the H I around NGC 4449 is still in an ongoing galaxy formation stage and the detected H I might fall into the galaxy. However, it seems to be unlikely that such large-scale, lopsided and regular structure can be maintained over such a long period, as e.g. an orbital period of more than 1 Gyr for the outer streamer. Already, the high internal velocity dispersion of 10 km s-1 would destroy the streamers on a timescale of 2·108 yr (Hunter et al. 1998). The alternative model is an interaction between NGC 4449 and DDO 125 explaining the streamer as a tidal feature. However, it is not obvious how the abrupt angles can be produced or why there is no bridge between both galaxies. Additionally, the optical images of both galaxies show no sign of any interaction like tidal tails.

In this paper we investigate the latter scenario addressing the question whether the observed structures can be formed by tidal interaction at all and what are the constraints on the mass distribution or orbits for both galaxies.

[Click here to see Fig. 1!]

2. Numerical models

2.1. Methods

Since the gravitational force in an encounter of two galaxies is assumed to be the dominant force, we apply an N-body approach for the dynamics of the H I. Though the effects of neglecting gasdynamics are not clear, such an ansatz has been used successfully for e.g. M81 and NGC 3077 (Thomasson & Donner 1993).

For our models we use two different N-body schemes: a self-consistent N-body integration using the GRAPE3Af special purpose computer (Sugimoto et al. 1990) and a restricted N-body simulation (Pfleiderer & Siedentopf 1961; Toomre & Toomre 1972). The basic idea of the latter method is to reduce the N-body problem into N 1-body problems by assuming that the gravitational potential is formed by a simple relation. Here we assume that the gravitational forces on each particle are given by the superposition of the forces exerted by two point-like objects (the galaxies) moving on Keplerian orbits. Hence one needs only (test) particles at the interesting halo locations instead of a large number of particles which have to create a stable and self-consistent galactic configuration. Additionally, the computation of the forces is much faster than for a consistent model. In total, the speedup of the numerical calculation is a factor of 1000 or more. On the other hand, the simplified treatment of the galactic potential prohibits any effect of the tidal interaction on the orbits of both galaxies, e.g. no transfer of orbital angular momentum into galactic spin is possible (i.e. no merging). In order to overcome this problem, we used the restricted N-body method only for scanning the parameter range and checked representative models by detailed self-consistent simulations. However, no significant differences have been found in the parameter range of interest.

2.2. Initial conditions

The initial conditions for the simulations characterize both, the galactic orbits and the distribution of the H I gas. Since the difference in heliocentric velocity between NGC 4449 and DDO 125 is only 10 km s-1, we are looking almost face-on the orbital plane (except for the unlikely case that both galaxies are now in closest approach). Hence, the actual distance between both galaxies is 42 kpc. The mass of NGC 4449 is fixed to the observed value within 32 kpc, i.e. MNGC 4449 = 7·1010 Msun. Thus, the unknown orbital parameters are the eccentricity e of the orbit, the mass ratio q ≡ MDDO 125 / MNGC 4449 and the minimum distance dmin at closest approach.

The original configuration of the H I-distribution is much poorly constrained. Assuming that the central ellipsoidal H I feature is the tidally unaffected remnant of an originally rotationally supported H I disk around NGC 4449 and presuming a prograde encounter, we can determine the rotational sense of the disk (NE-part is pointing towards us) and roughly its orientation, i.e. position angle. However, the inclination i is less clear. Moreover, the initial size Rmax of the disk or the regularity of the H I distribution is completely unknown.

3. Results

In a series of simulations we varied the orbital parameter (minimum distance dmin, the eccentricity ε, the mass ratio q) as well as the disk parameter (size Rmax of the disk, its position angle α and an inclination i). Qualitatively, the main features found by Hunter et al. (1998) can be reproduced by a nearly parabolic orbit, a mass ratio q = 0.15 and a minimum distance dmin = 30 kpc (upper left panel of Fig. 2). The parameters for the disk orientation are chosen in agreement with Hunter et al.'s suggestion of i = 60°, α = 230°, whereas the disk radius was set to Rmax = 27 kpc ≤ dmin. Though the numerical model is not a fit of the data, the characteristic sizes and locations of the streamers as well as the unaffected disk can be found. The material seen between the streamer and the disk which is not found in Hunter's data probably corresponds to the extended H I-mass seen in the single dish observations by van Woerden et al. (1975).

In order to check the significance of the model parameters the remaining three plots in Fig. 2 show the influence of changing the inclination (i = 80°), the size of the H I-disk (Rmax = 20 kpc) or the mass ratio (q = 0.05). We find that the relative size of the streamers or their orientation with respect to DDO 125 strongly depend on these parameters.

[Click here to see Fig. 2!]

4. Discussion

The numerical models are basically able to reproduce the morphology of the streamers, at least if one identifies the outer edges of the particle distribution with the location of the streamers. A more detailed investigation of the surface densities in the N-body simulations exhibits clearly the large stream pointing to the North-West, and also the stream emanating from the centre, however, the stream starting at the tip of the central stream pointing to the North is less clear, if existing at all. Another feature seen in the numerical models is a strong mass concentration at the southwestern edge of the disk which seems not to be in the observational data. Finally, the smaller mass concentration detected in the South-West cannot be explained by the interaction at all. On the other hand, the sensitivity of the models to variations of the initial parameters, especially to the variation of the observed disk orientation makes it unlikely that the successfully reproduced structures are formed by chance and no interaction takes place at all. Moreover, we think that the initial H I-distribution is more complex and the model of a completely rotationally supported disk-like H I distribution which is very smoothly distributed might be too simple.

Another important constraint is the total mass of DDO 125. Tully et al. (1978) derived from single dish observations a total mass of 5·108 Msun (corrected for a distance of 3.9 Mpc). Such a low mass would rule out the interaction scenario, because it can not produce sufficient tidal features. However, Ebneter et al. (1987) reported VLA-observations which show a linear increase of the rotation curve out to the largest radii where H I has been detected, i.e. out to twice the optical radius. This already allows for a lower mass limit of DDO 125 which is a factor 8 larger than the value of Tully et al., i.e. 5.7% of NGC 4449. Additionally, no flattening of the rotation curve nor any decline has been observed. Thus, the dynamical mass of DDO 125 is probably at least 10% or more of the mass of NGC 4449 which is in agreement with the mass ratio derived from the numerical model.

Acknowledgments. The simulations were partly performed with the GRAPE3af special purpose computer in Kiel (DFG Sp345/5).

References


Links (back/forward) to:

First version: 12thAugust,1998
Last update: 14thNovember,1998

Jochen M. Braun   &   Tom Richtler
 (E-Mail: jbraun|richtler@astro.uni-bonn.de)