format: text

text:
# The diffuse radio emission from the Coma cluster at 2.675 GHz and 4.85 GHz

## 1 Introduction

The diffuse, extended radio emission of the Coma cluster of galaxies (Abell cluster number 1656) was detected by Large et al. (large59 (1959)) and first investigated by Willson (willson70 (1970)). The prominent diffuse, extended radio source Coma C is the prototype of a cluster–wide radio halo. It was thoroughly investigated during the past years (e.g. Schlickeiser et al. schlickeiser87 (1987) (hereafter referred to as SST), Kim et al. kim90 (1990), Giovannini et al. gio93 (1993), Deiss et al. deiss97 (1997)). In the frequency range below 1.4 GHz there are many observations available that suggest a consistent power–law shape of the radio spectrum. SST present measurements of the halo source using the 100–m Effelsberg telescope111The 100–m telescope at Effelsberg is operated by the Max–Planck–Institut für Radioastronomie (MPIfR) on behalf of the Max–Planck–Gesellschaft (MPG). at 2.7 GHz. They found an integrated flux density significantly below that expected from the extrapolation from lower frequencies. The authors claim the existence of a spectral steepening of the radio spectrum of Coma C at frequencies above 1 GHz. Comparing predictions of the spectral shape of the three basic models for the radio halo formation (primary, secondary electron model, in–situ acceleration model) SST favour the latter for the Coma cluster. Observations at 1.4 GHz by Deiss et al. (deiss97 (1997)) with the Effelsberg telescope show the halo emission to be very extended. The subtraction of the contribution of point sources from the map results in a weak, diffuse radio component extended by more than 80$`\mathrm{}`$ in east–western direction, and by about 45$`\mathrm{}`$ north–south. The integrated flux density of this structure fits much better to the low–frequency part of the spectrum than to the 2.7 GHz point obtained by SST. Deiss et al. deemed the extremely sharp spectral break above 1.4 GHz as suggested by their own and the SST data to be unrealistic. The authors assume that either the 2.7 GHz map suffers from incorrect baselines, or that the mapped area was too small to cover the whole halo source. The shape of the radio spectrum provides constraints on the physical background of the halo formation. The work of Deiss et al. stimulated new sensitive observations of the diffuse emission in the Coma cluster at high frequencies ($`>`$1.4 GHz) in order to be able to derive constraints on the theory. The Effelsberg telescope is especially powerful to detect weak, extended emission.

The Coma cluster of galaxies is located at a redshift $`z=0.0232`$ (Struble & Rood struble91 (1991)), which implies a distance to the target of $`90h_{75}^1`$ Mpc and an angular scale of $`26h_{75}^1`$ kpc/arcmin (we adopt $`H_0=75`$ km/sec/Mpc throughout the paper).

## 2 Observations and data analysis

The observations reported here were carried out with the Effelsberg 100–m radio telescope between July 1998 and May 1999. We mapped the Coma cluster of galaxies at the frequencies of 2.675 GHz ($`\lambda 11.2`$cm) and 4.85 GHz ($`\lambda 6.2`$cm). Both receiver systems are installed in the secondary focus of the telescope. Their main properties can be found in Table 1.

We scanned the telescope in azimuth at both frequencies. The three Stokes parameters $`I`$, $`U`$ and $`Q`$ were recorded simultaneously. We obtained several coverages of the cluster area. At 2.675 GHz we mapped a large area, covering the central region and the extended source 1253+275. For the central region we added a couple of additional coverages. At 4.85 GHz we observed the center of the Coma cluster. The data reduction followed standard procedures for continuum observations with the 100–m telescope. For details of the observations consult Table 2. Table 3 lists the calibration sources together with the adopted parameters.

In Fig. 1 the radio map at 2.675 GHz showing the total intensity and the polarized emission is displayed. Fig. 2 zooms the center of the cluster. This map results from 33 coverages of the area, yielding an rms noise of 1.1 mJy/beam. The cluster center does not show any significant linear polarization at 2.675 GHz.

The map of the cluster at 4.85 GHz is shown in Fig. 3. One can clearly see the polarized emission of the central galaxies NGC 4869 ($`\alpha =12^\mathrm{h}56^\mathrm{m}56\stackrel{s}{.}5,\delta =28\mathrm{°}10^{}50^{\prime \prime }`$) and NGC 4874 ($`\alpha =12^\mathrm{h}57^\mathrm{m}10\stackrel{s}{.}7,\delta =28\mathrm{°}13^{}46^{\prime \prime }`$). Note that no extended emission above the 3$`\sigma `$ level is seen. Nevertheless the map reveals slightly extended emission at a level of 0.67 to 1.34 mJy/beam ($`\stackrel{}{=}1\mathrm{}2\sigma `$). In order to make it visible, the central part of the map is displayed in Fig. 3b down to zero flux. One can see extended structure around the two central galaxies.

## 3 Results

### 3.1 Morphology

#### The emission at 2.675 GHz

The 2.675 GHz map shows various point sources, emission by galaxies, extended diffuse emission in the center of the cluster, and the extended diffuse the source 1253+275 (see Sect. 5) southwest of the cluster center. There is no evidence for the bridge of low–brightness emission detected by Kim et al. (kim89 (1989)). Its intensity is obviously below the noise level of the 2.675 GHz map.

The central structure of the cluster emission is consistent with observations presented by Wielebinski (wieleb78 (1978)) and by SST. The map reveals a couple of point sources, mostly with low–level polarized emission. At the center one finds the radiation of the two galaxies NGC 4869 and NGC 4874 convolved to an extended source by the finite telescope beam. Diffuse emission surrounds the central maximum out to about $`20\mathrm{}`$. The highest sidelobes of the antenna pattern appear at about 19dB at a distance of about $`7\stackrel{}{.}5`$ from the center. Hence, this diffuse structure around the center is not created by sidelobes. Note that we could not apply the CLEAN algorithm (see e.g. Klein & Mack klein95 (1995)) because of the lack of a deep antenna pattern at this frequency due to the confusion by background sources. Comparing our map with that of the diffuse halo emission presented by Deiss et al. (deiss97 (1997)), as done in Figure 4, one clearly sees the smaller extension of the diffuse emission at 2.675 GHz. Our map is large enough to cover the whole halo source and to ensure emission–free area for baseline fitting. Our observations exhibit a large, frequency dependent extension of the halo source Coma C. Such a behaviour is to be expected in view of the spectral index map of the central region of the Coma cluster of galaxies presented by Giovannini et al. (gio93 (1993)). Those authors found an inner plateau of about $`8\mathrm{}`$ radius with a spectral index $`\alpha =0.8`$, whereas outside that central part $`\alpha `$ strongly steepens and reaches values larger than $`1.8`$. A comparison of both maps considering the different frequencies and beamsizes and the spectral index shows consistency of the intensity levels of the extended emission in the 2.675 GHz map with the 1.4 GHz map.

#### The emission at 4.85 GHz

Figure 3b shows the central part of the Coma cluster of galaxies at 4.85 GHz. Away from the point sources one finds only marginally extended low–brightness emission surrounding the two NGC–galaxies that are resolved at this frequency. In particular, no extended emission above the 3$`\sigma `$ level is detected.

### 3.2 Integrated diffuse radio flux

In order to estimate the integrated flux of the halo we first integrated the whole flux surrounding the center down to the zero level using polygons. In case of the 2.675 GHz map we included the bridges of emission towards the east and north of the center, as seen in Fig. 2, at the thinnest parts. We subsequently subtracted the fluxes of the point sources lying within the integration area. We calculated the point source flux densities on the basis of the spectral indices provided by Kim (kim94 (1994)). For those sources without any spectral index information we extrapolated the flux densities based on an assumed uniform spectral index of $`\alpha =0.8`$. All the extrapolated values are below 1$`\mathrm{}`$2 mJy, so that these faint sources sum up to only about 4% of the contaminating flux at 2.675 GHz, and about 5% at 4.85 GHz, respectively. The results are summarized in Table 4.

Most of the spectral indices in the Kim list are calculated based on low–frequency ($`<`$1.6 GHz) measurements. We expect that the spectral index increases with frequency ($`S\nu ^\alpha `$), so that our calculated point source fluxes will be overestimated. Such a behaviour is evident in case of the two central galaxies NGC 4869 and NGC 4874 (numbers 187 and 194 in the Kim list). The catalog contains 4.835 GHz fluxes for both sources given by observations. As shown in Table 5 their flux densities are much lower than those calculated with the derived spectral indices (resulting from best fits to the data). In assuming a constant spectral index one overestimates the point source flux and subsequently underestimates the flux density of the diffuse component.

### 3.3 The spectrum of Coma C

All data available in the literature together with their references are summarized in Table 6. In contrast to earlier publications about the flux density of Coma C one will not find the work by Waldthausen (waldthausen1980 (1980)). His 5 GHz observation of the region only covers the central part of the cluster and was done to check the morphology of the central sources. In particular Waldthausen did not mention any 5 GHz flux value in his thesis. The often cited value ($`<`$52 mJy) is Waldthausen’s 2.7 GHz upper limit for the Coma C flux.

The resulting spectrum of the halo source Coma C is displayed in Fig. 5. The new measurements of this work are symbolized with filled squares. Despite the fact that our 2.675 GHz point is slightly higher than that of SST, the break in the spectrum, first mentioned by those authors, is visible. The first 4.85 GHz observation supports this spectral steepening. A more detailed discussion of the spectrum will follow in the next chapter.

## 4 Discussion

The goal of our new observations was the inspection of the spectrum of the halo source. Is the steepening of the spectrum real as claimed by SST? The plot (see Fig. 5) clearly shows a break towards high frequencies. Both, the SST observation and our measurements yield flux densities below those expected from a power–law extrapolation of the low–frequency range. In order to quantify this behaviour, we followed the same procedure as done by SST. We have fitted the different basic models for halo source formation to the available data points. A $`\chi ^2`$–test identifies the model that provides the best fit to the data, and which hence has to be favoured.

In our analysis we considered the primary electron models of Jaffe (jaffe77 (1977)) and Rephaeli (reph77 (1977), reph79 (1979)), the secondary electron model of Dennison (denn80 (1980)) and the in–situ acceleration model described by Jaffe (jaffe77 (1977)), Roland (rol81 (1981)) and SST. Table 7 gives a summary of the predicted spectra. The notation is identical to that used by SST.

We normalized the data at $`\nu _\mathrm{j}=0.43`$ GHz and determined the most probable values for the parameters using the $`\chi ^2`$ technique. Concerning the different measurements with heterogeneous errors we used a weight for the fit of $`(1+\mathrm{d}I/I)`$ instead of $`\mathrm{d}I^2`$ not to overestimate the errors. The last two rows of Table 7 display the obtained reduced minimized $`\chi ^2`$ (=$`\chi _{\mathrm{min}}^2`$ per degree of freedom) and the corresponding best values for the free parameters of the models.

Fig. 6 shows the comparison of the calculations. The models of Jaffe (jaffe77 (1977)) and Dennison (denn80 (1980)) yield the same predicted spectra and are displayed in the same plot. Our data set contains a couple of new measurements compared to the work of SST. Hence, we obtain different values for the free parameters, but the main conclusion is the same. The in–situ acceleration model gives the best fit to the available data and has to be favoured over the primary and secondary electron models. In particular, the shape of the spectrum reveals a steepening above frequencies of 1 GHz. This was claimed by SST and is proved by the present work.

#### Physical parameters

The derived shape of the spectrum from the in–situ fit and our map at 2.675 GHz can be used to calculate the magnetic field strength assuming equipartition between the energy density of the relativistic particles and the magnetic field. We model the source as an ellipsoid with axes of the lengths of about 770 kpc, 600 kpc and 600 kpc, yielding an emitting volume of about 1.5$`\times 10^8`$ kpc3. A filling factor of unity is assumed. As the spectral index we use $`\alpha =0.79`$, which describes the initial power–law in the in–situ fit. In order to avoid any underestimation of the energy of the relativistic particles we have to use a low–frequency flux density given by the fit (we take S$`_{20\mathrm{M}\mathrm{H}\mathrm{z}}`$=60 Jy). In order to estimate the total energy we integrate the energy spectrum of the particles from 300 MeV to infinity. This differs from the generally used minimum–energy method (integration over the frequency spectrum from 10 MHz to 10 GHz), but seems more appropriate (see e.g. Beck et al. beck96 (1996)). The calculation leads to an equipartition magnetic field of

$$B_{\mathrm{eq}}=0.57(1+k)^{1/(\alpha +3)}h_{75}^{1/(\alpha +3)}\mu \mathrm{G},$$

where k is the energy ratio of the positively and negatively charged particles and $`\alpha `$ the spectral index (here 0.79). The resulting total equipartition energy density is

$$u_{\mathrm{tot}}=2.6\times 10^{14}(1+k)^{2/(\alpha +3)}h_{75}^{2/(\alpha +3)}\mathrm{erg}/\mathrm{cm}^3.$$

Note that changing k from 1 to 100 changes $`B_{\mathrm{eq}}`$ from 0.68 to 1.9 $`\mu `$G. These values are in agreement with the equipartition estimate by Giovannini et al. (gio93 (1993)) as well as with the work of Kim et al. (kim90 (1990)), who used rotation measures of radio sources projected on the Coma cluster. Nevertheless one has to keep in mind that an estimate of numbers based on the equipartition assumption only provides order of magnitude values.

With the radio spectrum now well known, we can calculate whether there is a balance between the internal (relativistic gas) and external (thermal gas) pressure. The latter comes from X–ray observations. The pressure of the relativistic component is given by

$$P_{\mathrm{rel}}=(\gamma 1)\times P_{\mathrm{part}}+P_{\mathrm{mag}},$$

where $`P_{\mathrm{part}}`$ and $`P_{\mathrm{mag}}`$ are the pressures of particles and fields, respectively, and $`\gamma `$ the ratio of specific heats. For relativistic particles, $`\gamma =4/3`$, hence $`P_{\mathrm{rel}}=4/3\times P_{\mathrm{mag}}`$ in case of energy equipartition between particles and fields. Inserting numbers, we obtain

$$P_{\mathrm{rel}}=1.7\times 10^{14}(1+k)^{2/(\alpha +3)}h_{75}^{2/(\alpha +3)}\mathrm{dyn}/\mathrm{cm}^2.$$

The thermal gas pressure in clusters of galaxies is given by

$$P(r)=P_0\left(1+\frac{r^2}{r_\mathrm{c}^2}\right)^{\frac{3\beta }{2}}$$

where $`P_0=2n_0kT`$ is the central pressure, $`r_\mathrm{c}`$ is the core radius, $`n_0`$ is the central electron density, $`\beta `$ is the ratio of galaxy to gas temperature, and $`T`$ the gas temperature (e.g. Feretti et al. feretti92 (1992)). At the core radius we thus obtain a thermal gas pressure

$$P_{\mathrm{th}}(r=r_\mathrm{c})=4.0\times 10^{11}h_{75}^{1/2}\mathrm{dyn}/\mathrm{cm}^2$$

using the intracluster plasma properties in Coma as derived by Briel et al. (briel92 (1992)):

$`n_0`$ $`=`$ $`(3.54\pm 0.05)\times 10^3h_{75}^{1/2}\mathrm{cm}^3`$
$`T`$ $`=`$ $`7.8_{1.8}^{+3.8}\mathrm{keV}`$
$`r_\mathrm{c}`$ $`=`$ $`(10.5\pm 0.6)\mathrm{arcmin}`$
$`\beta `$ $`=`$ $`0.75\pm 0.03.`$
The above calculation shows that the thermal pressure exceeds the relativistic one by a factor of about 500, even if we allow for protons and other nuclei in the relativistic gas, (in which case k = 100). It is thus obvious that the pool of relativistic particles forming the radio halo is strongly confined by the ambient thermal gas. This confinement was already reported for tailed and low–luminosity radio sources by Morganti et al. (mor88 (1988)) and by Feretti et al. (feretti92 (1992)).

Recently, magnetic field strengths in galaxy clusters have been determined by Clarke et al. (cla01 (2001)) using rotation measures of background sources. An average intracluster magnetic field strength of $`B5(\mathrm{}/10\mathrm{kpc})^{1/2}h_{75}^{1/2}\mu \mathrm{G}`$ has been obtained by them, using a simple tangled–cell model for the magnetic field with a constant coherence length $`\mathrm{}=10h_{75}^1\mathrm{kpc}`$. Such a high value appears to be in conflict with what is obtained if the extreme–ultraviolett (EUV) radiation and high energy X–ray (HEX) excess seen in Coma are due to inverse Compton (IC) scattering of background photons (see e.g. Enßlin & Biermann ens98a (1998)). In the case of the Coma halo this would imply magnetic field strengths of a fraction of a $`\mu G`$, much lower than what is implied by the Faraday rotation experiments, and also much lower than the values we obtain from the radio synchrotron emission.

The in–situ fit provides both the spectral index of the initial power–law and the cutoff frequency. Using the latter ($`\nu _\mathrm{S}=0.44`$GHz) together with the determined magnetic field, we are able to determine the lifetime of the radiating electrons, $`\tau =5.1\times 10^7(1+k)^{3/(2(\alpha +3))}h_{75}^{3/(2(\alpha +3))}\mathrm{yr}`$, taking into account synchrotron and inverse Compton losses. Note that the magnetic field that is equivalent to the cosmic microwave background radiation ($`B_{\mathrm{CMB}}=4(1+z)^2\mu \mathrm{G}`$) is much larger than the equipartition magnetic field, hence the energy losses via comptonization dominate over the synchrotron losses.

## 5 The source 1253+275

### 5.1 Introduction

Our 2.675 GHz observation covers a large region of the Coma cluster including the source 1253+275. This structure was first mentioned by Jaffe & Rudnick (1979) and subsequently investigated by Andernach et al. (andernach84 (1984)) and Giovannini et al. (gio85 (1985)gio91 (1991)). 1253+275 is the prototype of a radio relic source. This kind of radio sources can be found at the periphery of galaxy clusters, and they have elongated shapes and often exhibit highly polarized emission. The relevant portion of the 2.675 GHz map that contains the source 1253+275 is enlarged in Fig. 7. The extended nature of the source is clearly visible. The map also reveals strong polarized emission within the structure. For the map parameters consult the first row of Table 2.

Since the source 1253+275 appears to lack any optical counterpart (Andernach et al. andernach84 (1984)), its nature is still a puzzle. Giovannini et al. (gio85 (1985)) initially deemed the cluster member IC3900 to be the source of the plasma, but later on (Giovannini et al. gio91 (1991)), led by the detection of a bridge of low brightness emission (Kim et al. kim89 (1989)), they interpreted 1253+275 as a part of the large–scale radio emission present in the cluster–wide region of the Coma cluster and Coma–A1367 supercluster. Enßlin et al. (ens98b (1998)) favoured NGC 4789 as the possible production site of the relativistic particles. Particles released from this galaxy should feel the large–scale matter flow, directed towards the Coma cluster center. At the contact surface with the dense intracluster medium a shock should be present, thus forcing the particles to emit detectable synchrotron radiation and forming the extended radio source.

### 5.2 Morphology

Our map of 1253+275 is in agreement with that of Andernach et al. (andernach84 (1984)), at the same frequency. The source is extended for $``$800 kpc along the major axis (position angle $`135\mathrm{°}`$), perpendicular to the direction towards the cluster center. The maximum extent along the minor axis is $``$450 kpc. Linear polarization is seen in the northern and southern part. Andernach et al. (andernach84 (1984)) obtained internal and foreground rotation measures of zero, so that rotation of the polarization vectors by 90$`\mathrm{°}`$ discloses the projected orientation and structure of the magnetic field. It can be seen in Fig. 7 that the magnetic field is aligned parallel to the major axis of 1253+275, which fits well into the picture of a shock in this region. Such a shock should be directed towards the cluster center. Upstream to the shock, the plasma is compressed and the magnetic field is aligned perpendicular to the shock direction. The sharp edge of the structure towards the south–west is an additional hint at a compression of the present plasma.

### 5.3 Integrated diffuse radio flux

The radio flux of the diffuse emission of source 1253+275 was obtained by integrating all emission within the 3$`\sigma `$ level. The contribution of point sources was calculated using the source list and the spectral index values provided by Giovannini et al. (gio91 (1991)). The resulting total flux of the extended structure at 2.675 GHz (point sources subtracted) is $`112\pm 10`$ mJy. Table 8 contains the relevant data, where the contributions have been calculated for the northern and southern part of the source (corresponding to regions of polarized emission at the 3$`\sigma `$ level). The polarization properties are also listed.

The flux density of 1253+275 obtained by us is lower than that of Andernach et al. (andernach84 (1984)) if one subtracts the same point source contribution from their number. Our new measurement fits well to the other data available in the literature. The spectrum of the source is plotted in Fig. 8, with the data listed in Table 9. The best fit to the data yields a spectral index of $`\alpha =1.18\pm 0.02`$. This is the same result as that of Giovannini et al. (gio91 (1991)), though with a smaller error.

Using the equipartition assumption we derive a magnetic field strength of $`B_{\mathrm{eq}}=0.56(1+k)^{1/(\alpha +3)}h_{75}^{1/(\alpha +3)}\mu \mathrm{G}`$ in 1253+275, which implies a total energy density of $`u_{\mathrm{tot}}=2.5\times 10^{14}(1+k)^{2/(\alpha +3)}h_{75}^{2/(\alpha +3)}\mathrm{erg}/\mathrm{cm}^3`$. Here we use a filling factor of 1, the spectral index of $`\alpha `$=1.18 and the flux density of our measurement at 2.675 GHz. The emitting volume of about 1.4$`\times 10^8`$ kpc3 results from treating the source as a flat cylinder of height 650′′ (280 kpc, orientated towards the cluster center) and diameter 1850′′ (800 kpc, parallel to the major axis). Our numbers differ only slightly from those used by Andernach (andernach84 (1984)) and Giovannini et al. (gio91 (1991)).

## 6 Summary

We presented new single–dish radio observations of the Coma cluster of galaxies using the Effelsberg 100–m telescope. Our investigations were focused on the diffuse extended radio sources in the cluster. The radio halo source Coma C was observed at two frequencies (2.675 and 4.85 GHz). After subtracting the contribution of point sources we obtained new determinations of the flux densities of the diffuse structures. These verify the strong steepening of the radio continuum spectrum of Coma C at frequencies above 1 GHz, as first noticed by Schlickeiser et al. (schlickeiser87 (1987)). In comparing the three basic models for radio halo formation (primary electron model, secondary electron model and in–situ acceleration model), similar to the work done in Schlickeiser et al., we find that the latter most probably reflects the process at work in the Coma cluster. Fits of the models to the available data using a $`\chi ^2`$ test yield the best result for the in–situ acceleration model. This describes the spectral steepening also quantitatively.

Our 2.675 GHz map is extended enough to cover the entire diffuse emission. Our map is in good agreement with that published by Schlickeiser et al., thus settling a long–standing uncertainty (see the discussion of Deiss et al. deiss97 (1997)). Our 4.85 GHz measurement corroborates the spectral steepening towards high frequencies. Part of the previous confusion results from a misquotation in a couple of papers concerning the flux density of Coma C at 5 GHz. The often cited Ph.D. thesis of Waldthausen (waldthausen1980 (1980)) does not provide any 5 GHz flux density of Coma C. The value of less than 52 mJy, used in the citations as the total flux density at 5 GHz, is Waldthausen’s 2.7 GHz upper limit to the flux density of the extended component. Our 4.85 GHz observation provides the first firm measurement of the flux density of Coma C at this frequency.

Based on the equipartition assumption we calculate a field strength $`B_{\mathrm{eq}}=0.57(1+k)^{1/(\alpha +3)}h_{75}^{1/(\alpha +3)}\mu \mathrm{G},`$ and a total energy density of $`u_{\mathrm{tot}}=2.6\times 10^{14}(1+k)^{2/(\alpha +3)}h_{75}^{2/(\alpha +3)}\mathrm{erg}/\mathrm{cm}^3`$ for in Coma C.

Our 2.675 GHz map includes the region in which Kim et al. (kim89 (1989)) detected a bridge of low–brightness emission. Our map does not show any significant feature to be identified with the bridge. Its brightness is probably too low at 2.675 GHz to be detected.

We derive a new 2.675 GHz flux density for the source 1253+275, an extended diffuse radio relic located at the south–western periphery of the Coma Cluster. Using this new measurement, we calculate a spectral index of $`\alpha `$=1.18$`\pm 0.02`$, a value that is more accurate than previous determinations. For 1253+275 we estimate an equipartition magnetic field strength $`B_{\mathrm{eq}}=0.56(1+k)^{1/(\alpha +3)}h_{75}^{1/(\alpha +3)}\mu \mathrm{G}`$ and a total energy density of $`u_{\mathrm{tot}}=2.5\times 10^{14}(1+k)^{2/(\alpha +3)}h_{75}^{2/(\alpha +3)}\mathrm{erg}/\mathrm{cm}^3`$.

###### Acknowledgements.

MTH thanks Dr. Ch. Nieten for his assistance in programming the fit procedure. Dr. R. Beck is acknowledged for critical reading of the manuscript.