Bright Planetary Nebulae and their Progenitors in Galaxies Without Star Formation

Abstract: We present chemical abundances for planetary nebulae in M32, NGC 185, and NGC 205 based upon spectroscopy obtained at the Canada-France-Hawaii Telescope using the Multi-Object Spectrograph. From these and similar data compiled from the literature for other Local Group galaxies, we consider the origin and evolution of the stellar progenitors of bright planetary nebulae in galaxies where star formation ceased long ago. The ratio of neon to oxygen abundances in bright planetary nebulae is either identical to that measured in the interstellar medium of star-forming dwarf galaxies or at most changed by a few percent, indicating that neither abundance is significantly altered as a result of the evolution of their stellar progenitors. Several planetary nebulae appear to have dredged up oxygen, but these are the exception, not the rule. The progenitors of bright planetary nebulae typically enhance their original helium abundances by less than 50%. In contrast, nitrogen enhancements can reach factors of 100. However, nitrogen often shows little or no enhancement, suggesting that nitrogen enrichment is a random process. The helium, oxygen, and neon abundances argue that the typical bright planetary nebulae in all of the galaxies considered here are the progeny of stars with initial masses of approximately 1.5 Msun or less, based upon the nucleosynthesis predictions of current theoretical models. These models, however, are unable to explain the nitrogen enrichment or its scatter. Similar conclusions hold for the bright planetary nebulae in galaxies with ongoing star formation. Thus, though composition varies significantly, there is unity in the sense that the progenitors of typical bright planetary nebulae appear to have undergone similar physical processes. (Abridged)

\section{Introduction}

The planetary nebula luminosity function (PNLF) has proven remarkably successful as a distance indicator \citep[e.g.,][]{ciardulloetal2002}.  This success is based upon the constancy of the peak luminosity as a function of galaxy morphology, which is presumably a proxy for the mean age of the stellar progenitors of the brightest planetary nebulae \citep[e.g.,][]{marigoetal2004}.  The constancy of the peak luminosity suggests that it is either degenerate to a variety of stellar evolution parameters, or that some subset of planetary nebulae populations are common to the majority of galaxies.

Chemical abundances derived from spectroscopy of the brightest planetary nebulae in different galaxies provide a diagnostic tool to determine how similar these planetary nebulae are.  Such studies now exist for many galaxies within the Local Group.  Generally, the abundances in bright planetary nebulae trace those in the stellar populations of their host galaxy.  Host metallicities span a range of about an order of magnitude, as judged from oxygen abundances in planetary nebulae or from integrated light (see the references in \S \ref{other_galaxies} and Table \ref{table_pn_abun}).  Therefore, bright planetary nebulae do not arise from some stellar population that is common to all or most galaxies, at least as regards chemical composition.  These studies can be taken a step further, though, and used to investigate whether the stellar progenitors of these bright planetary nebulae underwent similar evolutionary processes.

It is well-known that the stellar progenitors of planetary nebulae modify their initial chemical composition through the course of their evolution.  Furthermore, the chemical abundances that are modified and the extent to which they are modified depend upon the star's initial mass and metallicity as well as the assumptions adopted concerning convective overshooting, at the lower boundaries of both the convective envelope and the pulse-driven convective zone during a thermal pulse \citep[e.g., ][]{forestinicharbonnel1997, marigo2001, herwig2005, karakaslattanzio2007}.  Typically, the matter that is returned to the interstellar medium via the nebular shell is enriched in helium and often in nitrogen, carbon, and s-process elements \citep[e.g., ][]{forestinicharbonnel1997, bussoetal2001, marigo2001}.  Recently, theoretical work has questioned whether oxygen is also produced \citep[e.g.,]{herwig2000} and some examples where this occurs have been found \citep[e.g., He 2-436 in the Sagittarius dwarf galaxy or the PN in Sex A;][]{pequignotetal2000,magrinietal2005}.

Here, we focus on the general problem of which elemental abundances are modified during the evolution of the progenitors of the brightest planetary nebulae.  The brightest planetary nebulae are of greatest interest because they are the ones that are being used to draw inferences about the chemical evolution of their hosts.  In particular, we argue that the great majority had progenitors that were not sufficiently massive to dredge up significant quantities of oxygen, which sets an upper limit to their initial mass, or that convective overshooting at the lower boundary of the pulse-driven convective zone is inefficient.  Primarily, we focus upon the planetary nebulae in galaxies where star formation ceased long ago, but then also consider the planetary nebulae in galaxies (or galactic components) with ongoing star formation.  In Section 2, we present new spectroscopic observations of planetary nebulae in M32, NGC 185, and NGC 205 as well as the reduction and analysis of these data.   In Section 3, we compile similar data for the planetary nebulae in other galaxies of the Local Group.  In Section 4, we determine what trends exist among elemental abundances and then proceed to interpret these in terms of the evolution of the progenitors of bright planetary nebulae in Section 5.  In Section 6, we present our conclusions.

\section{Observations and Analysis}

\subsection{Observations and Data Reductions}

The observations were obtained at the \facility{Canada-France-Hawaii Telescope} using the Multi-Object Spectrograph \citep[MOS;][]{lefevreetal1994} on UT 17-19 August 1998.  The MOS is a multi-object, imaging, grism spectrograph that uses focal plane masks constructed from previously acquired images.  The object slits were chosen to be 1\arcsec\ wide, but of varying lengths to accommodate the objects of interest.  In all cases, however, the slit lengths exceeded 12\arcsec.  We used a 600\,l\,mm$^{-1}$ grism that gave a dispersion of 105\,\AA\,mm$^{-1}$ and a central wavelength of 4950\,\AA.  The detector  was the STIS2 $2048 \times 2048$ CCD.  The pixel size was 21\,microns, yielding a spatial scale at the detector of 0\farcs 43\,pixel$^{-1}$.  The spectral coverage varied depending upon the object's position within the field of view, but usually covered the 3700-6700\,\AA\ interval.  Likewise, the dispersion was typically 2.1\,\AA\,pixel$^{-1}$, but varied depending upon the object's position within the field of view (from 1.9\,\AA\,pixel$^{-1}$ to 2.2\,\AA\,pixel$^{-1}$).  Images in the light of [\ion{O}{3}]$\lambda$5007 and the continuum at 5500\AA\ were acquired on UT 17 August 1998 and used to select objects for observation.  Spectra were obtained along columns as detailed in the observing log in Table \ref{table_observing_log}.  The photometric calibration was achieved using observations of BD+17$^\circ$4708, BD+25$^\circ$3941, BD+26$^\circ$2606, BD+33$^\circ$2642, and HD 19445 obtained over the course of the three nights (at least two standard stars were observed per night).  Spectra of the mercury, neon, and argon lamps were used to calibrate in wavelength.  Pixel-to-pixel variations were removed using spectra of the halogen lamp taken through the masks created for the standard stars and the three galaxies.

The spectra were reduced using the specred package of the Image Reduction and Analysis Facility\footnote{IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the Associated Universities for Research in Astronomy, Inc., under contract to the National Science Foundation.} (IRAF).  The procedure for data reduction followed approximately that outlined in \citet{masseyetal1992}.  Briefly, the mean of the overscan section was subtracted from all images.  A zero-correction image was constructed from overscan-subtracted bias images and was subtracted from all images.  The fit1d task was used to fit the columns of the overscan- and zero-correction-subtracted flat field images with a many-piece linear spline to remove the shape of the halogen lamp to form normalized flat field images.  Images of the standard stars and the galaxies were divided by the appropriate normalized flat field image.  The spectra of the planetary nebulae were then extracted and calibrated in wavelength using the spectra of the arc lamps.  Next, the spectra of the standard stars were used to calibrate in flux.  Finally, flux- and wavelength-calibrated spectra were co-added to produce the final spectrum for each object.

\subsection{Line Intensities}

For each planetary nebula, Tables \ref{table_n185_int}-\ref{table_m32b_int} present the raw line intensities normalized to H$\beta$ and their uncertainties ($1\,\sigma$) on a scale where the H$\beta$ line has an intensity of 100.  The line intensities were measured using the software described in \citet{mccalletal1985}.  This software simultaneously fits a sampled gaussian profile to the emission line(s) and a straight line to the continuum.  The quoted errors include the uncertainties from the fit to the line itself, the fit to the reference line (H$\beta$), and the noise in the continuum for both lines.  For those lines where no uncertainty is quoted, the intensity given is an upper limit ($2\,\sigma$).

Tables \ref{table_n185_int}-\ref{table_m32b_int} also include the reddening determined from $\mathrm H\alpha/\mathrm H\beta$ and $\mathrm H\gamma/\mathrm H\beta$, assuming intrinsic ratios appropriate for the temperature and density observed.  The temperature and density finally adopted for each object is given in Tables \ref{table_n185_abun}-\ref{table_m32_abun} (see \S 2.3 for details).  Normally, both estimates of the reddening agree within errors, implying that there are no severe errors in the flux calibration.  The \citet{fitzpatrick1999} reddening law was used, parametrized with a total-to-selective extinction of 3.041.  This parametrization delivers a true ratio of total-to-selective extinction of 3.07 when integrated over the spectrum of Vega \citep{mccall2004}, which is the average value for the diffuse component of the interstellar medium of the Milky Way \citep{mccallarmour2000}.  The uncertainties quoted for reddenings, temperatures, and densities are derived from the maximum and minimum line ratios allowed considering the uncertainties in the line intensities.

The line intensities were corrected for reddening according to

\begin{equation}
\log \frac{F(\lambda)}{F(\mathrm H\beta)} = \log \frac{I(\lambda)}{I(\mathrm H\beta)} -0.4 E(B-V)(A_1(\lambda)-A_1(\mathrm H\beta))
\end{equation}

\noindent where $F(\lambda)/F(\mathrm H\beta)$ and $I(\lambda)/I(\mathrm H\beta)$ are the observed and reddening-corrected line intensity ratios, respectively, $E(B-V)$ is the reddening determined from $F(\mathrm H\alpha)/F(\mathrm H\beta)$, when available, and $A_1(\lambda)$ is the extinction in magnitudes for $E(B-V)=1$\,mag, i.e., $A_1(\lambda) = A(\lambda)/E(B-V)$, $A(\lambda)$ being the reddening law \citep{fitzpatrick1999}.  The values of $A_1(\lambda)$ used for all lines are given in column 3 of Tables \ref{table_n185_int}-\ref{table_m32b_int}.  The optical depth at $1 \, \rm \mu m$, which is a much clearer descriptor of the amount of obscuration, can be computed by multiplying $E(B-V)$ by 1.054 \citep[see][]{mccall2004}.

A few of the final spectra of the planetary nebulae in M32 and NGC 205 had anomalously low ratios of $\mathrm H\alpha/\mathrm H\beta$.  In all cases, these spectra were from slits close to the edges of the field of view (in the spectral direction).  The affected spectra are those of PN1 and PN3 in NGC 205 and PN7 and PN30 in M32.  These spectra may be identified in Tables \ref{table_n185_int}-\ref{table_m32b_int} since a measurement is given for H$\alpha$, but no reddening is computed from the $\mathrm H\alpha/\mathrm H\beta$ ratio.  In those cases, the line intensity for H$\alpha$ was scaled to the value expected given the reddening found from the $\mathrm H\gamma/\mathrm H\beta$ ratio.  Also, the line intensities for lines in the red, $\lambda \ge 5876$\AA, were calculated relative to H$\alpha$ and scaled by the same factor as H$\alpha$.  In principle, this effect should be due to some undocumented vignetting within the spectrograph.  The principal uncertainty resulting from our correction scheme will be regarding the abundance of He$^{\circ}$ for the affected objects.  N$^+$ abundances should not be affected due to the proximity of the [\ion{N}{2}] lines to H$\alpha$.

All of the planetary nebulae observed here were previously known from the list of \citet{ciardulloetal1989}, save the object \# 3273 from the list of \cite{merrettetal2006}.  The finding chart and coordinates for this object are given in Fig. \ref{fig_new_PN}.

\subsection{Physical Conditions and Chemical Abundances}\label{chemical_abundances}

Tables \ref{table_n185_abun}-\ref{table_m32_abun} present the electron temperatures as well as the ionic and elemental abundances derived for each object.  The atomic data employed for ions of N, O, and Ne are listed in Table \ref{table_atomic_data}.  For H$^{\circ}$ and He$^+$, the emissivities of \citet{storeyhummer1995} were used.  For He$^{\circ}$, we used an extended list of the emissivities from \citet{porteretal2005} that was kindly provided by R. Porter.

The ionic abundances were derived using the SNAP software package \citep{krawchuketal1997}.  The electron temperature adopted is always based upon the [\ion{O}{3}]$\lambda\lambda$4363/5007 ratio, the upper limit being used when [\ion{O}{3}]$\lambda$4363 was not detected.  An electron density of 2000\,cm$^{-3}$ was adopted for all objects, based upon the distribution of densities in Galactic planetary nebulae found by \citet{riesgolopez2006}.  When detected, the [\ion{S}{2}]$\lambda\lambda$6716,6731 lines are usually compatible with this density, within uncertainties.  If only an upper limit is available for the lines of a given ion, that limit is used.

Using a lower limit for the electron temperature will underestimate the abundances of oxygen, nitrogen, and neon ions (collisionally-excited lines), while it will slightly overestimate the abundances of helium ions (recombination lines).  If limits are available for both an abundance diagnostic and electron temperature, it is not strictly possible to know whether the result is a limit to the ionic abundance.  However, since a metal diagnostic varies exponentially with temperature, but linearly with abundance, it is likely that the result is still a lower limit to the ionic abundance.

The elemental abundances were calculated from the ionic abundances using the ionization correction factors (ICFs) proposed by \citet{kingsburghbarlow1994}.  In this scheme, the $\mathrm N/\mathrm O$ and $\mathrm{Ne}/\mathrm O$ abundance ratios that are tabulated in Tables \ref{table_n185_abun}-\ref{table_m32_abun} are simply the $\mathrm N^+/\mathrm O^+$ and $\mathrm{Ne}^{2+}/\mathrm O^{2+}$ ionic ratios, respectively.  We are interested primarily in the helium and oxygen abundances and in the $\mathrm N/\mathrm O$ and $\mathrm{Ne}/\mathrm O$ abundance ratios, so the ICFs matter most for oxygen.  The ICF for oxygen is a function of the He$^+$ and He$^{2+}$ abundances, and so of the intensities of \ion{He}{1}\,$\lambda$5876 and \ion{He}{2}\,$\lambda$4686.  We almost always measure \ion{He}{1}\,$\lambda$5876 (PN30 in M32 is the only exception), but often only have an upper limit to \ion{He}{2}\,$\lambda$4686.  When we only have upper limits to the helium lines, our ICF for oxygen will be slightly overestimated.  Nonetheless, since the ICFs are usually near unity (PN30 in M32 is the only exception), any error in the oxygen abundance should also be small.  Since the ICFs required to compute the nitrogen and neon abundances are linear functions of the oxygen abundance, any overestimate in the oxygen abundance translates into similar overestimates of the nitrogen and neon abundances, though the $\mathrm N/\mathrm O$ and $\mathrm{Ne}/\mathrm O$ ratios are unaffected.

In Tables \ref{table_n185_abun}-\ref{table_m32_abun}, the uncertainties quoted for all quantities account for the uncertainties in the line intensities involved in their derivation, including the uncertainties in the reddening.  However, the quoted uncertainties in the elemental abundances do not include the uncertainties in the ICFs.

In the calculation of the elemental abundances, we made an effort to adopt a scheme that treated all objects as homogeneously as possible.  Occasionally, this may have resulted in less than optimal elemental abundances for some objects.  However, it has the advantage of uniformity, i.e., any drawbacks are common to all objects.  For example, it is feasible to use singlet lines of \ion{He}{1} to compute $\mathrm{He}^+/\mathrm H$ in some objects, but not in all.  To avoid the risk of introducing spurious differences, we chose to use common lines for all objects.  The one exception to this rule is the $\mathrm O^+/\mathrm H$ ionic abundance.  Generally, the ionic abundance based upon [\ion{O}{2}]$\lambda$3727 was used in favor of that based upon [\ion{O}{2}]$\lambda\lambda$7319,7331, but the latter lines were used when [\ion{O}{2}]$\lambda$3727 was not available.

\subsection{Reliability of the Line Intensities and Elemental Abundances}

The uncertainties tabulated for our line intensity measurements are the formal uncertainties associated with the analysis of the data.  Given the difficulty of obtaining the data with 4-m-class telescopes, it is helpful to compare with independent measurements of the same objects to gauge the absolute uncertainties.  Most valuable to this cause are measurements by \citet{richermccall1995} and \citet{richeretal1999}.  In Fig. \ref{fig_m32_comp_int}, we compare our line intensities for planetary nebulae in M32 in common with \citet{richeretal1999}.  Generally, the agreement is reasonable.  Lines stronger than H$\beta$ from both studies agree to within $\pm 10$\%, while lines of order $0.1 I(\mathrm H\beta)$ have relative uncertainties of $\sim 100$\%.  A similar trend is found when comparing the current data for planetary nebulae in NGC 185 and NGC 205 with those presented in \citet{richermccall1995} though fewer data are involved and the measurements of \citet{richermccall1995} are considerably more uncertain than those of \citet{richeretal1999}.  \citet{richermccall2007a} discuss this issue further.

Oxygen abundances for some of the planetary nebulae in M32 were previously presented by \citet{richeretal1999}.  We compare the oxygen abundances found here with that study in Fig. \ref{fig_m32_comp_o}.  Considering that many of the previous abundances were lower limits due the availability of only upper limits for the electron temperatures (i.e., [\ion{O}{3}]$\lambda$4363 was often not previously detected), we find excellent agreement.  In the two cases where  there were previously measured temperatures, the oxygen abundance agrees in one case and is apparently discrepant in the other.  However, the uncertainties quoted here include the uncertainty in reddening whereas the uncertainties quoted by \citet{richeretal1999} do not.  Therefore, this discrepancy is smaller than it appears.  Where there were previously temperature limits and now measured temperatures, the current oxygen abundances are higher, as expected.  Finally, where there is currently still only a limit to the temperature, the oxygen abundance is higher, again as expected given our better sensitivity (and a lower temperature limit).

\section{Data for Planetary Nebulae in Other Galaxies}\label{other_galaxies}

In order to construct a more representative sample of bright planetary nebulae in galaxies without star formation, we compiled the data available for planetary nebulae in the bulge of M31, M32, NGC 147, Sagittarius, and Fornax from the literature \citep{walshetal1997, jacobyciardullo1999, richeretal1999, rothetal2004, kniazevetal2006, zijlstraaetal2006, goncalvesetal2006}.  For comparison with bright planetary nebulae in galaxies with star formation, we adopt the data compilation from \citet{richermccall2007} augmented with data for planetary nebulae in the disk of M31 from \citet{jacobyford1986} and \citet{jacobyciardullo1999} as well as the three planetary nebulae reported  here (PN4, PN17, and M3273 in the field of M32).  Table \ref{table_pn_abun} summarizes our abundances for the planetary nebulae in all of these galaxies.  The total sample consists of 84 and 82 bright planetary nebulae in galaxies with and without star formation, respectively.

The data compiled from the literature were analyzed in the same fashion as described in \S \ref{chemical_abundances} for our own data.  One difference, however, is the absence of upper limits to undetected lines of [\ion{O}{3}]$\lambda$4363 and \ion{He}{2}\,$\lambda$4686.  When a temperature or its limit could not be determined, no abundances were computed.  When \ion{He}{2}\,$\lambda$4686 was not detected, the ICF for oxygen is a lower limit (unity), leading to lower limits to the abundances of oxygen, nitrogen, and neon, though $\mathrm N/\mathrm O$ and $\mathrm{Ne}/\mathrm O$ are not affected.

We restrict our samples of planetary nebulae to those that are intrinsically brightest, within 2\,mag of the brightest in each galaxy.  In principle, these should be the objects with the spectra of the highest quality within each galactic sample.  We note, however, that all of the planetary nebulae in NGC 147 have intrinsic luminosities considerably fainter than those included from the other galaxies \citep{corradietal2005}.  Selecting objects bright in [\ion{O}{3}]$\lambda$5007 affects the planetary nebula sample in two known ways (and perhaps in others yet unknown).  First, as Fig. \ref{fig_lo3} illustrates, at least for oxygen abundances below that of the interstellar medium (ISM) in the SMC, high [\ion{O}{3}]$\lambda$5007 luminosity favors the most oxygen-rich planetary nebulae in each galaxy \citep{dopitaetal1992, richermccall1995}, and so may preferentially select objects derived from recent star formation.  Second, high [\ion{O}{3}]$\lambda$5007 luminosity also favors planetary nebulae early in their evolution \citep[e.g., ][]{jacoby1989, stasinskaetal1998, schonberneretal2007}.

\section{Abundance trends}

Figure \ref{fig_ne_o} presents the correlation between neon and oxygen abundances for bright planetary nebulae in different galaxies.  To help improve the clarity of this and the following figures, error bars are only shown for the data presented here.  For the planetary nebulae in other galaxies, those in dwarf irregulars have somewhat lower uncertainties (especially in the case of Leo A and Sextans A and B), those in Fornax and Sagittarius have much smaller uncertainties, and those in M31 have slightly larger uncertainties.   The dashed line shows the Ne-O relation for bright planetary nebulae while the solid line is the relation between these two elements in the interstellar medium in emission-line galaxies \citep{izotovetal2006}.  A fit of the Ne-O relation for the bright planetary nebulae in galaxies without ongoing star formation (M31 bulge, M32, dwarf spheroidals) yields

\begin{equation}
12 + \log(\mathrm{Ne}/\mathrm H) = (1.016\pm 0.051) X + (-0.87\pm 0.43),
\end{equation}

\noindent where $X=12+\log(\mathrm O/\mathrm H)$.  Within uncertainties, this fit is identical to that found by \citet{richermccall2007} for the planetary nebulae in dwarf irregular galaxies.  Formally, the slope and intercept differ from the values found by \citet{izotovetal2006}, but the difference between these relations is minimal, never exceeding 0.05\,dex, over the abundance range in common, $7.5\,\mathrm{dex}<12+\log(\mathrm O/\mathrm H)<8.5\,\mathrm{dex}$.

Figure \ref{fig_no_o} plots the $\mathrm N/\mathrm O$ ratio as a function of oxygen abundance.  The solid line is the limiting value of $\mathrm N/\mathrm O$ found in star-forming galaxies \citep{izotovetal2006}.  Although the uncertainties are substantial, since $\mathrm N/\mathrm O$ is determined from the $\mathrm N^+/\mathrm O^+$ ratio, both of which are minority species, the scatter does appear to be real.  In particular, no anomalously low values are found, suggesting confidence in the $\mathrm N/\mathrm O$ ratios.

Figure \ref{fig_he_o} presents the helium abundance as a function of the oxygen abundance.  The solid line is a fit to these elemental abundances in emission-line galaxies from \citet{oliveetal1997}.  This relation is a reasonable approximation to a lower envelope for the abundances in planetary nebulae, demonstrating that many planetary nebulae do not modify their original helium content.  Note that this relation is only calibrated to $12+\log(\mathrm O/\mathrm H)\sim 8.2$\,dex.
Plotting helium abundances as a function of the $\mathrm N/\mathrm O$ ratio yields no new information.  In particular, no significant correlation is found, though this may be partly due to the large uncertainties in $\mathrm N/\mathrm O$.

In Fig. \ref{fig_nne_neo}, we attempt to discern cases where oxygen might have been produced by the stellar progenitors by plotting $\mathrm N/\mathrm{Ne}$ versus $\mathrm{Ne}/\mathrm O$.  Planetary nebulae whose progenitors produced oxygen should appear to the left in this diagram.  Their vertical position will depend upon whether nitrogen production also took place.  We normalize relative to neon since neon is not expected to be modified by the stellar progenitors of most planetary nebulae \citep{karakaslattanzio2003}.  Only PN6 and perhaps PN4 in NGC 205 appear to be possible examples of planetary nebulae whose progenitors produced oxygen.  From Fig. \ref{fig_ne_o}, it would appear that PN56 and FJCHP57 in M31's bulge and disk, respectively, might also have produced oxygen.  Unfortunately, neither nitrogen abundances nor uncertainties for the oxygen and neon abundances are available for these two objects.

\section{Nucleosynthesis in the Progenitors of Bright Planetary Nebulae}

\subsection{Oxygen and Neon}

On cosmic scales, the production of oxygen and neon is dominated by core collapse supernovae \citep{clayton2003, herwig2004}.  Nonetheless, models of the progenitors of planetary nebulae predict that such stars potentially produce both.  At very low metallicities, e.g., below one hundredth of the solar metallicity, the progenitors of planetary nebulae are capable of synthesizing the most abundant isotope of neon, $^{20}$Ne \citep{herwig2004, karakaslattanzio2007}.  At all metallicities, stars with masses in the range $\sim 2-4\,M_{\sun}$ produce $^{22}$Ne in quantities that can significantly change the observed total neon abundance in planetary nebulae, e.g., change the total abundance by at least $+0.1$\,dex.  Finally, depending upon the physics adopted to model semiconvection and convective overshoot in AGB stars, stars with masses of $\sim 2-3\,M_{\sun}$ and metallicities comparable to those in the Magellanic Clouds \emph{may} dredge up significant quantities of $^{16}$O \citep[our qualitative interpolation of][]{herwig2000,herwig2004}.  Note that these models achieve significant oxygen dredge-up because of the semi-convection attributed to the pulse-driven convective zone that arises with each helium shell flash and mixes out carbon and oxygen from the outer part of the helium-free core \citep{herwig2000}.  If this convective overshooting is efficient, oxygen is expected to be dredged up along with carbon during the third dredge-up; otherwise no oxygen is dredged up \citep{marigo2001, karakaslattanzio2007}.  Note that all of these models indicate that oxygen destruction should occur for the progenitors of planetary nebulae with the highest masses.  Below 2\,$M_{\sun}$, current models predict that changes in oxygen and neon abundances should be slight.

Nonetheless, as \citet{karakaslattanzio2003} concluded, the progenitors of planetary nebulae are unlikely to mimic the Ne-O relation found in star-forming galaxies.  Any oxygen production in low- or intermediate-mass stars would proceed through the same channel as in high-mass stars, the endpoint of helium burning, via the $^{12}\mathrm C(\alpha,\gamma)^{16}\mathrm O$ reaction.  In high-mass stars, neon is a product of carbon burning \citep{clayton2003}.  In low- or intermediate-mass stars, any $^{22}$Ne produced occurs during He-shell burning of CNO-processed material \citep[two alpha captures on $^{14}$N, e.g.,][]{herwig2004} while (at very low metallicities) any $^{20}$Ne arises primarily from $\alpha$ particle captures on $^{16}$O, $^{16}\mathrm O(\alpha,\gamma)^{20}\mathrm{Ne}$ \citep[e.g.,][]{herwig2004}.  Thus, while oxygen production will result in the dominant isotope in low- or intermediate-mass stars, most of the neon production would yield a minority isotope.  It would therefore be highly coincidental if oxygen and neon production in low- or intermediate-mass stars produced an oxygen-to-neon ratio identical to that resulting from the nucleosynthesis of high-mass stars.

Figures \ref{fig_ne_o} and \ref{fig_nne_neo} argue that oxygen and neon abundances are coupled in the majority of bright planetary nebulae.  Since bright planetary nebulae follow the same Ne-O relation as emission line galaxies, where the nucleosynthetic yields of type II supernovae fix the Ne-O relation, the progenitors of bright planetary nebulae either modify both abundances to the same extent or they do not modify either abundance.  As we have just argued, the available evidence indicates that the former is highly unlikely.  In accord with \citet{karakaslattanzio2003}, we conclude that the progenitors of bright planetary nebulae typically do not modify either their oxygen or neon abundances.

Although significant oxygen production is uncommon, both Figs. \ref{fig_ne_o} and \ref{fig_nne_neo} indicate that, occasionally, oxygen is produced during the evolution of the stellar progenitors in galaxies where star formation ceased long ago.  This conclusion was also reached by \citet{richermccall2007} for the progenitors of bright planetary nebulae in dwarf irregulars, and we extend it here with the example of FJCHP57 in the disk of M31.  In general, however, these objects are the exception, not the rule.  Note that \emph{all these objects are characterized by unusual Ne/O ratios}.

Based upon luminosity considerations, \citet{ciardulloetal2005} have argued that bright planetary nebulae should descend from stars of about $2\,M_{\sun}$.  (See \citet{richeretal1997} for a counter-example.)  If bright planetary nebulae do indeed descend from progenitors of about $2\,M_{\sun}$, many of the bright planetary nebulae in M32, NGC 185, NGC 205, Fornax, Sagittarius, and the bulge of M31 might be expected to be enriched in oxygen, since many of them have oxygen abundances intermediate between those of the SMC and LMC.  These objects, however, have neon and oxygen abundances in accord with the ratio observed in star-forming galaxies, a result that suggests that their oxygen content has not been modified.  Therefore, either these bright planetary nebulae descend from progenitors of less massive stars, and hence are less susceptible to oxygen enrichment, or convective overshooting is not as efficient as some models suggest.

\citet{pequignotetal2000} have argued in favor of oxygen production in He 2-436 in Sagittarius, with its progenitor having increased its oxygen content by $7\pm 4$\%.  Such a small increase would not be detectable here, in general, so we cannot rule out production at this level from our abundance measurements.  Interestingly, \citet{dudziaketal2000} suggest an initial mass of $\sim 1.2\,M_{\sun}$ for the progenitor of He 2-436.  If this is typical, oxygen production is expected to be slight, regardless of the assumptions concerning convective overshooting \citep{marigo2001, karakaslattanzio2007}.

While our individual oxygen and neon abundances are not sufficiently precise to detect small changes in oxygen abundances, the large sample at our disposal should allow us to detect small systematic differences with respect to the neon and oxygen abundances in emission-line galaxies.  The slope of our Ne-O relation differs slightly from that found by \citet{izotovetal2006} (Fig. \ref{fig_ne_o}).  Our slope is consistent with unity, as has been found for bright planetary nebulae in dwarf irregular galaxies \citep{richermccall2007}.  \citet{izotovetal2006} found a slope for the Ne-O relation in star-forming galaxies that is significantly greater than unity.  This difference is small, but formally significant, and may be interpreted in a variety of ways.

If neither oxygen nor neon are affected by the evolution of the stellar progenitors and if neither oxygen nor neon production in supernovae changes with metallicity, a slope of unity is expected for the Ne-O relation.  \citet{izotovetal2006} attributed their larger slope to the depletion of oxygen onto dust grains in \ion{H}{2} regions at higher oxygen abundances.  However, a slope different from unity might also arise from changes in the reddening law associated with the dust in \ion{H}{2} regions or a dependence upon metallicity of the relative production of oxygen and neon in type II supernovae.  If the slope found by \citet{izotovetal2006} is due to either dust depletion or changes in the reddening law, the slightly shallower slopes found for bright planetary nebulae imply that there is no systematic oxygen production in their progenitors, even in small quantities.  If, instead, the slope found by \citet{izotovetal2006} reflects the true, changing, production of oxygen and neon by massive stars, the difference between emission-line galaxies and planetary nebulae points to a low level of self-enrichment in the progenitors of bright planetary nebulae proportional to the original oxygen abundance, of the order of a few percent.

\subsection{Helium}

Helium abundances span a significant range, even at a given oxygen abundance within a single galaxy (Fig. \ref{fig_he_o}).  However, the locus observed for \ion{H}{2} regions in star-forming dwarfs \citep{oliveetal1997} is a very reasonable fit to the lower boundary of the distribution.  Thus it is reasonable to presume that, at a given oxygen abundance, the minimum helium abundance was the initial helium abundance for the progenitors of these planetary nebulae and that the spread in helium abundances is the result of varying enrichments due to the evolution of the progenitors.

We define helium enrichment as the increase in helium abundance observed relative to the helium abundance expected at the same oxygen abundance for dwarf star-forming galaxies \citep{oliveetal1997}, i.e., $\mathrm{fractional\ enrichment} = ((\mathrm{He}/\mathrm H)_{PN,obs} -(\mathrm{He}/\mathrm H)_{ISM})/(\mathrm{He}/\mathrm H)_{ISM}$.  In Fig. \ref{fig_he_dist}, we plot histograms of the helium enrichment.  We separate bright planetary nebulae into three samples, those in dwarf irregulars, those in dwarf spheroidals and M32, and those in the bulge of M31.  The planetary nebulae in the first two samples have oxygen abundances that overlap the range studied by \citet{oliveetal1997} in the \ion{H}{2} regions in star-forming dwarfs.

Clearly, the planetary nebulae in all samples are enriched by less than 50\% on average, \emph{if} the relation between oxygen and helium abundances found by \citet{oliveetal1997} holds.  There is a tendency to observe larger helium enrichments in the planetary nebulae in galaxies where star formation ceased long ago, but a Kolmogorov-Smirnov test indicates that the difference between the samples in dwarf irregulars and in M32/dwarf spheroidals is not significant at even the 10\% level.  Considering theoretically-expected yields \citep[e.g., ][]{marigo2001}, the often low helium enrichments presumably imply that the progenitors of bright planetary nebulae are of rather low mass, probably below 1.5\,$M_{\sun}$, especially in those systems whose metallicities are below that of the LMC.  The bright planetary nebulae in the bulge of M31 appear to imply more enrichment at higher metallicities, but, at these metallicities, the \citet{oliveetal1997} relation must be extrapolated to compute the enrichment factor (see Fig. \ref{fig_he_o}), so this result is uncertain.

\subsection{Nitrogen}

The interpretation of nitrogen is complicated.  The history of star formation plays an important role, because nitrogen is produced by both type II supernovae and low- and intermediate-mass stars such as novae and AGB stars, as seen in studies of the ISM in star-forming galaxies \citep[e.g.,][]{vanzeeetal1998}.  Because of these two nitrogen sources, $\mathrm N/\mathrm O$ can vary with time in a complicated way in a galaxy's ISM.  For instance, for a constant rate of star formation, the $\mathrm N/\mathrm O$ ratio in a galaxy's ISM will initially be fixed at a low value as a result of type II supernovae.  As low- and intermediate-mass stars begin to contribute, the $\mathrm N/\mathrm O$ ratio will rise, eventually attaining some equilibrium value.  For a more complicated star formation rate, the $\mathrm N/\mathrm O$ ratio need not vary monotonically with time.  Furthermore, the $\mathrm N/\mathrm O$ ratio in the ISM is merely the initial value of this ratio in the stellar progenitor of a planetary nebula.  If this progenitor star synthesizes nitrogen, the resulting $\mathrm N/\mathrm O$ ratio will be even higher.  In principle, the $\mathrm N/\mathrm O$ ratio is expected to be a function of the mass of the stellar progenitor \citep[e.g.,][]{marigo2001}, so, for a constant star formation rate, one expects a trend of increasing $\mathrm N/\mathrm O$ with increasing oxygen abundance in each galaxy.

No such simple result is found.  In M32, planetary nebulae spanning a wide range of $\mathrm N/\mathrm O$ are found at all oxygen abundances.  In the bulge of M31, there are planetary nebulae with similar $\mathrm N/\mathrm O$ ratios whose oxygen abundances differ by an order of magnitude.  Clearly, then, the $\mathrm N/\mathrm O$ ratio is not a monotonic function of the progenitor mass, nor is it solely the result of the chemical evolution of the host galaxy.  The scatter of $\mathrm N/\mathrm O$ as a function of oxygen abundance is not entirely unexpected (see Fig. \ref{fig_no_o}), as it is seen in bright planetary nebulae in dwarf irregular galaxies \citep{richermccall2007}.  Likewise, planetary nebulae in the bulge of the Milky Way span a large range in $\mathrm N/\mathrm O$ \citep{cuisinieretal2000, escuderocosta2001, exteretal2004, gornyetal2004, escuderoetal2004}.

In galaxies without star formation, the minimum $\mathrm N/\mathrm O$ ratio is similar to the minimum value found in star-forming galaxies up to an oxygen abundance of $12+\log(\mathrm O/\mathrm H)\sim 8.5$\,dex (data are limited at higher $\mathrm O/\mathrm H$).  A simple interpretation of this result is that  the progenitors of their planetary nebulae may have had a similar initial $\mathrm N/\mathrm O$ ratio, i.e., the value in the ISM was dominated by the nitrogen production of type II supernovae.  If so, the factor by which the progenitor of any planetary nebulae enriched its nitrogen content may be found by comparing the observed $\mathrm N/\mathrm O$ ratio with the limiting value in emission line galaxies.  This is the maximum enrichment factor, since the progenitor's initial $\mathrm N/\mathrm O$ ratio is unknown and could have been higher.  Adopting this definition, the progenitors of the planetary nebulae in galaxies without star formation appear to have undergone nitrogen enrichment that varied by more than an order of magnitude.  Although  the amount of nitrogen enrichment appears to vary randomly, it could conceivably be related to some secondary characteristic of the progenitor stars, such as rotation or binarity.

Since oxygen consumption is not evident in either Figs. \ref{fig_ne_o} or \ref{fig_nne_neo}, nitrogen production presumably occurs exclusively at the expense of carbon for the bright planetary nebulae in galaxies without ongoing star formation.  If the largest nitrogen enrichments occur in progenitors with initial $\mathrm N/\mathrm O$ ratios similar to the limiting value in emission line galaxies, the enrichment may be up to a factor of 100.  If so, reprocessing of endogenous carbon within the progenitors would be necessary, e.g., via hot bottom burning or some similar process, since the carbon abundance in the interstellar medium is less than 10 times that of nitrogen \citep{garnett2002}.  Such reprocessing of carbon is not predicted in the mass range expected for the progenitors of bright planetary nebulae in systems where star formation stopped long ago.

For $12+\log(\mathrm O/\mathrm H)> 8.3$\,dex in Fig. \ref{fig_no_o}, the planetary nebulae from galaxies without star formation apparently dominate the lowest values of $\mathrm N/\mathrm O$.  If the effect is real, it may result from differences in the chemical evolution of their host galaxies rather than of differences in their stellar progenitors.  The higher $\mathrm N/\mathrm O$ ratio in planetary nebulae in galaxies with ongoing star formation does not necessarily indicate systematically greater enrichment in these planetary nebulae, since many of these planetary nebulae have $\mathrm N/\mathrm O$ ratios compatible with those in the interstellar medium of their host galaxies, i.e., compatible with no enrichment \citep{richermccall2007}.  Therefore, the stellar progenitors in star-forming galaxies need not have produced more nitrogen, but simply may have been born with a higher initial nitrogen content.

\subsection{More Generally}

A limitation of this entire discussion is our ignorance of the mass range of the progenitors of bright planetary nebulae.  In principle, models of nebular luminosities and the expected evolution of the nebular envelopes provide some restrictions, but these are so contradictory as to not be helpful \citep[e.g.,][]{richeretal1997, marigoetal2004, ciardulloetal2005, schonberneretal2007}.  Ideally, the chemical abundances observed in bright planetary nebulae would allow us to constrain the masses of their stellar progenitors.

In galaxies without star formation, planetary nebulae derived from the highest mass progenitors, $M_i > 3\,\mathrm M_{\sun}$, should be absent.  In this case, theory predicts little enrichment in nitrogen \citep[observations of red giants imply enrichment by a factor of up to $\sim 4$, e.g., ][]{grattonetal2000}, and up to a doubling of the helium content.  Large enhancements in the oxygen abundance (or none at all) could be common, depending upon the exact masses of the progenitor stars and what physics actually occurs \citep[e.g.,][]{marigo2001, karakaslattanzio2007}.  Likely, changes in neon abundances would be small.

Observations agree with the predictions for helium abundances, though the enhancement is usually less than 50\%\ (Fig. \ref{fig_he_o}).  Comparison of oxygen and neon abundances, however, imply that only rarely is oxygen synthesized in quantities exceeding a few percent of the initial abundance. Considering the helium, oxygen, and neon abundances, the bright planetary nebulae in galaxies without star formation are generally consistent with their progenitors having initial masses of about 1.5\,$M_{\sun}$ or less, given current theoretical models \citep[e.g.,][]{marigo2001,herwig2005,karakaslattanzio2007}.  The main argument against their progenitors having initial masses near 2\,$M_{\sun}$ is the rarity of planetary nebulae whose progenitors produced significant quantities of oxygen, unless convective overshoot is inefficient \citep{marigo2001, karakaslattanzio2007}.  Likewise, the rarity of unusual $\mathrm{Ne}/\mathrm O$ ratios implies masses below 2\,$M_{\sun}$.

The nitrogen abundances in bright planetary nebulae do not agree with the theoretical predictions.  Planetary nebulae in galaxies without star formation span a wide range of nitrogen enrichment, from none at all to enrichment exceeding that in red giants by an order of magnitude.   The largest nitrogen abundances agree with theoretical expectations only if the progenitors of these planetary nebulae had masses exceeding 3\,$M_{\sun}$ \citep{marigo2001,herwig2005}, a possibility excluded by the history of star formation in these galaxies.  Even if such masses were allowed, there would remain the problem of very large and very small enrichment at a given progenitor mass (oxygen abundance), e.g., M32 in Fig. \ref{fig_no_o}, since such randomness is not predicted at any mass.  Also, if such high masses were invoked to explain the nitrogen abundances, much higher helium abundances should be observed, and they are not.

Fig \ref{fig_nne_neo} further emphasizes the difficulty with nitrogen production.  All of the objects where oxygen has apparently been dredged up are highly enriched in nitrogen.  No models predict significant enrichment in both oxygen and nitrogen simultaneously.

If low mass stars really produce as much nitrogen as these observations suggest, these stars may be a much more important source of this element than has been previously considered.  Nitrogen abundances hinge upon the $\mathrm N^+/\mathrm O^+$ ionic ratio.  Both ions are minority species in bright planetary nebulae, so the ICF needed to derive the nitrogen abundance is large (see Tables \ref{table_n185_abun}-\ref{table_m32_abun}) and the resulting abundances uncertain.  If the nebulae are matter-bounded, the problem would be even more severe.  However, the uncertainties for our $\mathrm N/\mathrm O$ ratios are driven by uncertainties in line intensities.  Typically, the [\ion{O}{2}] and [\ion{N}{2}] lines are weak and so uncertain, as is also the case for [\ion{O}{3}]$\lambda$4363, which is needed to compute the temperatures.  Remedying this situation would require deeper spectroscopy.  The effort is worthwhile, however, since the low mass progenitors of planetary nebulae represent a substantial fraction of all stellar mass, e.g., the $1-1.5\, M_{\odot}$ mass range represents about 16\%\ of the total mass above 1\,$M_{\odot}$ for a \citet{salpeter1955} initial mass function and an upper mass limit of 100\,$M_{\odot}$.

Although the foregoing discussion focussed primarily upon the bright planetary nebulae in galaxies where star formation ceased long ago, the bright planetary nebulae in galaxies with ongoing star formation have very similar properties.  They too have modest enrichment in helium, negligible enrichment in oxygen, and a wide range of nitrogen enrichments.  Therefore, the progenitors of bright planetary nebulae in galaxies with ongoing star formation must also have low masses if current theory is to be trusted.

Thus, chemical abundances argue that the masses of the stellar progenitors of the brightest planetary nebulae are relatively low in all galaxies.  This is not necessarily incompatible with the high [\ion{O}{3}]$\lambda$5007 luminosities that they achieve, as argued by \citet{richeretal1997}.  In their (somewhat extreme) example of a $0.55 M_{\sun}$ central star, 13\%\ of the ionizing radiation need be converted to [\ion{O}{3}]$\lambda$5007 photons, which is similar to the maximum efficiency of 10\%\ expected theoretically \citep{dopitaetal1992,marigoetal2004}.  The unknown, as \citet{marigoetal2004} make clear, is the lack of quantitative knowledge of the evolution, particularly the time scales, of the planetary nebula system.  So, while we may not yet understand why the brightest planetary nebulae descend from these particular progenitor stars, the similarity of these stars from galaxy to galaxy is likely one of the underpinnings of the success of the planetary nebula luminosity function as a distance indicator.

The progenitors of bright extragalactic planetary nebulae span a wide range of chemical compositions, reflecting those of the progenitor stellar populations in their host galaxies.  That they enrich the resulting planetary nebulae similarly implies that they underwent similar physical processes during their evolution.  Presumably, then, the mass range of these progenitor stars is slightly different from galaxy to galaxy, reflecting the different metallicities of their stellar populations.  If so, the epoch of star formation that dominates the production of bright planetary nebulae should vary somewhat from galaxy to galaxy.

Based upon the foregoing, however, it would be a mistake to argue that the masses (or the mass range) of the progenitors of bright planetary nebulae in all galaxies are exactly the same.  As shown by \citet{stasinskaetal1998}, the distribution of \ion{He}{2}$\lambda$4686/H$\beta$ intensity ratios is dramatically different for bright planetary nebulae in the LMC and the bulge of M31.  This conclusion appears to hold generally between environments with and without star formation \citep{richer2006}, implying that the distribution of central star temperatures differs between the two environments.  Normally, this is interpreted in terms of the central star mass, but might also result from differences in envelope masses \citep{stanghellinirenzini2000}, or both.  As argued by both \citet{stasinskaetal1998} and \citet{richer2006}, it is likely that the ionization structure of the nebulae is also different in environments with and without star formation.  Differences in central star masses or in envelope evolution could easily arise from a difference in the mass range of the progenitor stars in different environments.

It would not be unusual if there were surprises lurking regarding the nucleosynthesis in the progenitors of fainter planetary nebulae, especially in galaxies with ongoing star formation.  Theory clearly predicts that the planetary nebulae descending from the progenitors of the highest masses should be depleted in oxygen because of hot bottom burning at the base of their convective envelopes \citep{marigo2001}.  Some theoretical work predicts that oxygen production should be common for planetary nebula progenitors with masses of order 2-3\,$M_{\odot}$ at metallicities comparable to or lower than the SMC's \citep{marigo2001,herwig2004}.  Dwarf irregular galaxies should, therefore, be the typical galactic hosts of planetary nebulae with unusual oxygen abundances.   Observationally, however, individual examples are not that common among the brightest planetary nebulae \citep{richermccall2007}, so presumably they are to be found among those of lower [\ion{O}{3}]$\lambda$5007 luminosity.  The exact luminosities and abundances where oxygen enrichment is common have yet to be identified.

\section{Conclusions}

We have obtained spectroscopy of a sample of bright planetary nebulae in M32, NGC 185, and NGC 205.  From these and similar data for planetary nebulae in other Local Group galaxies, we calculate their chemical composition.  These planetary nebulae are generally within 2\,mag of the peak of the luminosity function.

We find that helium abundances in bright planetary nebulae are typically changed by less than 50\% compared to the initial abundances in their progenitors.  Oxygen and neon abundances in bright planetary nebulae are not significantly modified, since they follow the relation between neon and oxygen abundances found in the interstellar medium in star-forming galaxies.  There is the possibility of a slight systematic increase in the oxygen abundance by a few percent, but this conclusion depends upon how one interprets the slight variation of the ratio of neon to oxygen as a function of oxygen abundance in star-forming galaxies.  In contrast, we often find very strong nitrogen enhancements in bright planetary nebulae, though very small enhancements are also found.  These nitrogen enhancements are not obviously related to the progenitor mass or to the chemical evolution of the host galaxy.  At present, they appear to be random.  In reality, they are likely to be linked to some secondary characteristic of the stellar progenitors, such as rotation or binarity.  These conclusions apply independently of whether we consider planetary nebulae in galaxies with or without ongoing star formation.

Considering the helium, oxygen, and neon abundances, bright extragalactic planetary nebulae in all galaxies appear to descend from relatively low mass progenitors, of approximately 1.5\,$M_{\sun}$ or less, based upon the calibrations of current theoretical models.  If their progenitors were significantly more massive, of order 2\,$M_{\sun}$ as has been argued on the basis of luminosity considerations, significant enrichment in oxygen could be relatively common.  A few examples of planetary nebulae with oxygen enrichment are found, but are rare.

The large nitrogen abundances often found suggest more massive stellar progenitors, {\emph i.e.}, exceeding 3\,$M_{\sun}$ (so as to allow the production of nitrogen via hot bottom burning).  Such masses are clearly excluded by the history of star formation in galaxies where star formation stopped long ago, and in any case don't explain the observation of both large and small nitrogen enhancements at a given oxygen abundance.  Also, we find no example of a bright extragalactic planetary nebula with very large abundances of both helium and nitrogen, which should be a clear signature of a high mass stellar progenitor.  This is the case even for galaxies where star formation is ongoing.  Thus, low- and intermediate-mass stars are a more important source of nitrogen than has been hitherto considered.  Despite their difficulty, deeper observations would be very worthwhile to derive more secure abundances.

Overall, the progenitors of all bright extragalactic planetary nebulae appear to have evolved similarly.  They were of sufficiently low mass that helium enrichment was modest and oxygen enrichment was minimal, if indeed it even generally takes place.  Nevertheless, despite this apparently low mass, these stars produced much greater quantities of nitrogen than are predicted by current models.  Presumably, then, the mass range of these stellar progenitors is similar in all galaxies, but displaced somewhat due to metallicity.  Consequently, there may be more unity to the properties of bright extragalactic planetary nebulae than has been apparent previously.  If this is the case, bright planetary nebulae also presumably probe a particular epoch during the evolution of their host galaxies.

\acknowledgments

We thank the time allocation committee of the Canada-France-Hawaii Telescope for granting us the opportunity to observe.  We thank R. Porter for providing us with a more extensive list of He$^{\circ}$ emissivities.  MGR acknowledges financial support from CONACyT through grant 43121 and from UNAM-DGAPA via grants IN112103, IN108406-2, and IN108506-2.  MLM thanks the Natural Sciences and Engineering Research Council of Canada for its continuing support.


Meson exchange models for Meson production

Abstract: The production of mesons in nucleon-nucleon collisions is reviewed from the viewpoint of the meson-exchange picture. In the first part various possible production mechanisms and their relative importance are discussed. In addition, general features of meson production are described.In the second part special emphasis is put on pion production. Implications of chiral perturbation theory are discussed. Results based on a specific meson-exchange model for pion production in nucleon-nucleon collisions are presented. The model is utilized to calculate several spin-dependent observables of the reactions $pp \to pp\pi^0$ and $pp \to pn\pi^+$ such as spin-correlation parameters. This allows us to study the role of the Delta resonance at close to threshold energies.The talk closes with a brief discussion of the production of $\phi$ mesons as an example for the production of heavier mesons. A strategy is outlined that allows to extract information on the structure of the nucleon from the reaction $NN \to NN \phi$.

\section*{Introduction}

With the advent of new accelerator technology measurements with
very high accuracy became possible for meson production close to
 threshold \cite{DATA}. This improvement opens the possibility for the
investigation of a lot of interesting physics phenomena. It is the
goal of this talk to highlight a few of those.
This will we done using the meson exchange picture.

The largest part of the presentation will concentrate on pion production.
As this is the first and most important inelasticity of the
nucleon--nucleon system it shall deal as a guideline for models
of the production of heavier mesons. In addition, for this channel
there are not only analyzing powers but also double polarization
data available. A study of this system will therefore allow us
to get a feeling about what can be learned from meson production.
In particular it will be demonstrated that using polarization observables
will allow to study resonances that couple too weakly to
significantly influence the unpolarized observables -- in this example
it will be the $\Delta$ (1235) far below the resonance position.
However, the same kind of sensitivity
should be found for resonances that
have only a minor effect on the
total cross section even at their resonance position.
 In addition: in the unpolarised observables
almost all the effects of the Delta--isobar are driven by the
$^1D_2 \to ^5S_2$ $NN \to N\Delta$ transition leading to p--wave
$\pi$--production\footnote{This transition is prohibited for
the $pp \to pp\pi^0$ reaction and this is why the Delta
is less important for the neutral pion production (see below).}. Polarization observables
allow to study less dominant transition amplitudes selectively and
thus allow direct insight into the NN interaction.

There is one special feature to the production of pions: since they are
pseudo Goldstone bosons of the spontaneous broken chiral symmetry,
there is the possibility of studying pion production in
nucleon--nucleon collisions with chiral perturbation theory.
I will briefly compare the features of the meson exchange
approach to those of chiral perturbation theory.

As an example for the production of heavier mesons
I will continue with a short discussion of how one can use
data on the production of vector mesons close to threshold to
deduce information on the structure of the nucleon.

The presentation will close with a brief outlook.

\section*{General features of meson production in nucleon--nucleon
collisions}

To produce a meson in nucleon nucleon collisions the initial
kinetic energy needs to be large enough to put the outgoing meson on its mass
shell. The initial center of mass momentum required to produce
a meson of mass $m$ turns out to be at least
$
p = \sqrt{Mm+\frac{m^2}{4}} \ ,
$
where $M$ denotes the nucleon mass.
If the kinematics is chosen close to the threshold the outgoing
particles are (almost) at rest. Thus, large
characteristic for the meson production close to threshold.
This leads to a big momentum mismatch between the final
wave-function and the initial one leading to a suppression
of the direct
 production mechanisms (c.f. Fig. \ref{beitraege}a and d for pion production).
Therefore short range processes  dominate the production (c.f. Fig. \ref{beitraege}b
\footnote{For large space like momentum transfers, the $\pi N$ T--matrix enters
far off--shell and thus also the pion exchange contribution should
be considered as a short range process.} and c).
Note, however, once one moves away from threshold
the range of allowed final momenta increases. Thus moderate momentum transfers
become more and more likely and higher partial waves come into play.
Thus, one should expect very
different production mechanisms to be dominant as a function
of the excess energy.
The range of energies just described can be identified with
the range of $0 \leq \eta \leq 1$, where $\eta$ denotes the maximum meson
momentum in units of the meson mass. Therefore, in this presentation we
will mainly concentrate on this regime.

\begin{figure}[h]
\vspace{4.2cm}
\special{psfile=beitraege2.eps
hoffset=30 voffset=0 hscale=60 vscale=60}
\caption{Possible pion production mechanisms taken into account in our model:
(a) direct production; (b) pion rescattering; (c) contributions from pair
diagrams; (d) and (e) production involving the excitation of
the $\Delta (1232)$ resonance. Note that diagrams where the
$\Delta$ is excited after pion emission are also included.}
\label{beitraege}
\end{figure}

The goal of the meson exchange picture is to find a description
of meson production that is both transparent enough to be able
to identify the relevant physics of the process and still contains
all the relevant dynamics.
In this context
it is important to understand
the special role that the nucleon--nucleon system plays
in the initial as well as the final state.

As early as 1952 Watson pointed out the significance that the final state
nucleon--nucleon interaction should play in meson production processes
close to threshold \cite{watson}. To see this let us consider  meson
absorption on a two nucleon system, which is related to  meson production
through detailed balance. The cross section is a measure of the volume
in which all three particles must be simultaneously in order to make
the absorption process possible. However, since the elastic scattering
cross section for the two nucleons is very much larger than the production
cross section and the NN interaction close to threshold is strongly
attractive, a two step process is favored in which first the two
nucleons approach at close range -- driven by the NN interaction --
before the meson comes into play.
Therefore the energy dependence
of the absorption (production) matrix element should be given by the energy
dependence of the initial (final) nucleon--nucleon interaction.
If we neglect the effect of the meson--nucleon interaction
and use the effective range expansion for the nucleon--nucleon interaction
the energy dependence of the total production cross--section close to threshold
should be given by
\begin{equation}
\sigma_{NN \to NNx}(\eta ) \propto \int_0^{m_x \eta} dq' q'{}^2
\frac{p'}{1+a(a+r_o)p'{}^2} \ ,
\label{endep}
\end{equation}
where $a$($r_o$) denotes the nucleon--nucleon scattering length
(effective range) in the relevant channel and $p'$ and $q'$
denote the relative momenta of the two nucleons and of the meson
with respect to the two nucleons respectively.
It was shown in ref. \cite{MuS} that the inclusion of the
Coulomb effect is important, if
for a two proton final state.
The inclusion of the
Coulomb effect in eq. \ref{endep}
is straightforward \cite{meyer2}. It turns
out that the energy dependence of the total cross section
that follows from eq. \ref{endep}
is in nice agreement with almost all experimental evidence on
meson production near
threshold \cite{DATA}. The only exception is the production of $\eta$ mesons, since
the $\eta N$ interaction is too strong to be neglected (see
also ref. \cite{ulfneu}).

The above argument shows that it is important to include
the final state interactions of the produced particles
 in order to describe the
energy dependence of the total cross sections correctly. Reversing this
statement indicates, that it is possible to relate the energy dependence
of the near threshold production cross sections to low energy parameters
of the dominant interactions. Examples for this are given in ref.
\cite{janpkl} for the $\Lambda N$ interaction and in refs.
 \cite{Green,Grishina} for the $\eta N$ interaction.

Care should be taken, however, when trying to deduce
a general formula like eq. \ref{endep} not only to
 describe the energy dependence
but also to set the normalization of the effect of the FSI
independent of the production operator.
This statement is illustrated by different groups using
different prescriptions for the normalization
of the FSI factor that differ by more than an order of
magnitude, although both use the meson exchange picture to calculate
the production operator (c.f. e.g. refs. \cite{ulfneu} and \cite{sibir2}).
 Based on quite
general arguments one can show, that there is no model independent way
to fix the normalization introduced by the FSI factor \cite{fsicrit}
\footnote{The method applied in refs. \cite{ulfneu,sibir2} was also critizised
in ref. \cite{jounifsi} from a different point of view, namely, that the
inclusion of the FSI through using a plane wave that is modified in strength
according to the on--shell scattering data
misses the short range correlations introduced by the FSI and thus
leads to wrong relative weights of the individual contributions.}.

A similar reasoning that lead to the successful description of the
energy dependence of the total cross section leads to an expression
that describes the effect of the initial state interaction \cite{fsicrit}:
In order to allow the production to proceed, the two incoming nucleons have to come very close
to each other. There are two mechanisms that potentially prevent this,
namely elastic and inelastic scattering.
If we focus on the meson exchange picture, it can be argued that the
principal value integral that introduced the scheme dependence for the FSI
is actually small for high energies\footnote{Remember: there is no such thing as a
model independent splitting of the NN interaction and the production operator.}.
Based on this observation one can derive a rough estimate for
a suppression factor induced by the ISI
(for a discussion of the range of applicability see \cite{fsicrit})
\begin{equation}
\lambda = \frac{1}{4}(1+\eta_L)^2-\eta_L \sin ^2(\delta_L) \ ,
\end{equation}
where $\eta_L$ and $\delta_L$ denote the inelasticity and the phaseshifts
of the partial wave in the initial state relevant for the meson production
under consideration. Using the numbers given by the SAID database \cite{said}
we find for the $\eta$ and $\eta '$ production values of $\lambda = \frac{1}{5}$
and $\lambda = \frac{1}{3}$ respectively. Note, that the former number covers
the bulk effect of the initial state interaction calculated within a meson
exchange model \cite{lee}. Thus, the ISI leads to a non negligible suppression
of the total cross section and should be taken into account.

\section*{Pion production and the role of the $\Delta$ close to threshold}

As this is meant to be an overview presentation, this section starts
with a brief history of the field.

Essentially all recent theoretical investigations on pion production
near threshold are based on the model proposed by
Koltun and Reitan in 1966 \cite{KuR}. In this model two production
mechanisms are considered:  direct production
(Fig.~\ref{beitraege}a), and pion rescattering (Fig.~\ref{beitraege}b) --
the latter, however, in on--shell approximation, where the $\pi$N--T--matrix
is replaced by the scattering length. It is important to note that --
as a consequence of chiral symmetry -- the isoscalar scattering
length\footnote{The $\pi$N--T--matrix can be written as a linear combination
of an isoscalar and an isovector part.} almost vanishes.
This model was utilized by Miller and Sauer in 1991 \cite{MuS}
to analyze the first set of high precision data of the
reaction $pp \rightarrow pp\pi^0$ near threshold that became available
from IUCF \cite{IU1}.
It turned out that such a model grossly underestimates
the empirical cross section \cite{MuS}. Note, that isovector rescattering
is prohibited in this reaction channel.

In 1992 Niskanen extended this model by including the $\Delta$ (1235)
isobar (cf. Fig.~\ref{beitraege}d and e) \cite{Nis1}. Furthermore, he allowed for
an energy dependence in the s-wave rescattering term, however,
still keeping the on--shell approximation \cite{Nis2}. These
improvements roughly doubled the predicted $pp\rightarrow pp\pi^0$
cross section, but Niskanen's results still underestimate the IUCF data
by a factor of 3.6.

Another new production mechanism was introduced by Lee and Riska
in 1993 \cite{LuR}. These authors considered effects from meson-exchange
currents due to the exchange of heavy mesons
that excite a nucleon--antinucleon pair, as shown in Fig.~\ref{beitraege}c.
It was found that the resulting contributions (in particular the
one of the $\omega$ meson) enhance the pion production cross section
by a factor of 3-5 \cite{LuR,HMG} and thus eliminate most of the
under prediction found in earlier investigations.

However, in 1995 Hern\'andez and Oset presented an alternative
explanation for the missing strength in the $\pi^0$ production close
to threshold \cite{HO}. These authors took into account the
off-shell properties of the $\pi N$ amplitude in the
evaluation of the rescattering diagram.
It turned out, that this enhances the rescattering contribution sufficiently
to describe the empirical $pp \to pp\pi^0$ cross--section without
the inclusion of additional short ranged diagrams.
In the same year our group performed a calculation using a microscopic
model for the $\pi$N T--matrix \cite{Han1}. Although the effect
of the rescattering was not as large as reported in \cite{HO} its
enhancement due to the off--shell part of the T--matrix was quite significant.

\begin{figure}[t]
\vspace{9cm}
\special{psfile=csresults.eps
hoffset=40 voffset=-10 hscale=60 vscale=60}
\caption{\it{Total cross sections for the different
pion production channels.
The solid line represents the result
of the  full model of ref. \protect\cite{Han98a}, whereas for the dashed line the contributions
involving the $\Delta$ where omitted.
For a compilation of the experimental data see ref.
\protect\cite{DATA}.
}}
\label{csresults}
\end{figure}

Since the pion is the Goldstone Boson of the chiral symmetry one might
hope that the framework of chiral perturbation theory allows more insight
into this fundamental reaction. In the literature there is a series
of papers available, that use tree level chiral perturbation theory to calculate
the cross sections close to threshold for both neutral pion production
\cite{CP1,CP2,CP2b,CP3} and charged pion production \cite{HHH,CP4}.
In addition there are two calculations available that calculate the
$\pi^o$ production to one loop order \cite{moalem,fred} -- these, however,
use heavy baryon chiral perturbation theory which is inappropriate in
the present kinematics \cite{ulfneu}.
The advantage of chiral perturbation theory compared to the meson exchange
picture is, that one has an organizing principle at hand. Although the
expansion parameter is quite large in case of the $NN \to NN\pi$ reactions,
namely $Q=\sqrt{m/M} \simeq \frac{1}{3}$ \cite{CP2},
it is very illuminating to study the
lowest order contributions and to compare those to the diagrams depicted
in Fig. \ref{beitraege}.

Based on the counting rules employed in ref. \cite{CP2} the lowest order
contributions to pion production close to threshold (O($Q$)) are for the neutral
pion production the direct production from the nucleon as well as
from the delta (Fig. \ref{beitraege}a and d), whereas in case of the
charged pion production at the same order there appears a rescattering
diagram (Fig. \ref{beitraege}b) involving the isovector $\pi N$ interaction
in addition \cite{HHH,CP4}.
As it was argued above, the large momentum transfer characteristic of
meson production close to threshold leads to a suppression of the
direct production mechanisms in favor of rescattering and short
range processes. Therefore, one expects the production of charged pions to
be under better control than the one of neutral pions -- that is
indeed what we experience: the lowest order calculation for the $\pi^+$ production
deviates from the data by a factor of 2 \cite{HHH,CP4} whereas the one for $\pi^0$ production
deviates by more than an order of magnitude (the latter occurs because the
individual contributions are small and there is a cancelation between
the $\Delta$ contribution (Fig. \ref{beitraege}d) and the nucleonic one (Fig. \ref{beitraege}a)
 \cite{CP2}).

\begin{figure}[t]
\vspace{10cm}
\special{psfile=pp2pppi0_pol.ps
hoffset=50 voffset=-10 hscale=60 vscale=60}
\caption{\it{Spin correlation parameters for the reaction
$pp \to pp\pi^0$, where $A_\Sigma = A_{xx}+A_{yy}$ and
$A_\Delta = A_{xx}-A_{yy}$.
Curves are the same as in Fig. \protect\ref{csresults}.
The experimental data are taken from
Refs.~\protect\cite{Mey98,Mey99}.
}}
\label{Corr1}
\end{figure}

At the next order (O($Q^3$)) in both channels loops enter\footnote{This
is true only when using the counting scheme of ref. \cite{CP2}.} as well
as contact interactions (the effective field theory
analog of the short range contributions of the meson exchange picture) and
 rescattering diagrams. The strength of the contact interactions
is not constrained by chiral symmetry. Therefore the authors of refs.
\cite{CP2,CP2b,CP4} tried to estimate their strength by means of resonance
saturation: the contact interactions were identified with short ranged
diagrams (Fig. \ref{beitraege}c) that were evaluated using parameters
from realistic $NN$ potentials. However, it is not clear a priori, that
this procedure makes sense for the contact interactions.
The parameters relevant for the rescattering can be related to
$\pi N$ scattering \cite{bkm}.

Taking all this together it should not come as a surprise that the s--wave pion
production is not well under control theoretically. As it was argued
before the typical momentum transfers decrease as we move away from the
production threshold and it is the large momentum transfer that causes all
the trouble.
A more appropriate approach to pion production close to threshold thus seems
to be to first understand the higher partial waves and then to move down
in energy. This procedure should allow one to study the onset of new
physics in a more controlled way what
highlights the importance of looking at higher partial
waves. It should be noted that
one expects
 the chiral expansion to work better, when the pions are produced in  higher
partial waves\cite{biraprivate}. Therefore in what follows we will concentrate on polarization
observables for here the higher partial waves show up most clearly.

To be more concrete in what follows we will compare the results of a particular
meson exchange model \cite{Han98a} to the data.
 In the literature there is only one additional
model that considers all the different pion production channels as well as
higher partial waves, namely the one by the Osaka group.
This model was described in a different presentation in some detail  \cite{Tamura},
and thus we will not discuss it here.

In the J\"ulich model \cite{Han98a} all standard pion-production mechanisms
(direct production (Fig.~\ref{beitraege}a),
pion rescattering (Fig.~\ref{beitraege}b),
contributions from pair diagrams (Fig.~\ref{beitraege}c)
are considered -- the latter, however, are treated
as a parametrisation of all the missing short range mechanisms.
In addition, production mechanisms involving the
excitation of the $\Delta (1232)$ resonance
(cf. Fig.~\ref{beitraege}d,e) are taken into account explicitly.
All $NN$ partial waves up to orbital angular momenta $L_{NN} = 2$, and
all states with relative orbital angular momentum $l \leq 2$ between
the $NN$ system and the pion are considered in the final state.
Furthermore all $\pi N$ partial waves up to orbital angular momenta
$L_{\pi N} = 1$ are included in calculating the rescattering diagrams
in Fig.~\ref{beitraege}b,e. Thus, our model includes not only
s-wave pion rescattering but also contributions from
non resonant p-wave rescattering.
For details about the model we refer the reader to ref. \cite{Han98a}.

\begin{figure}[t]
\vspace{10cm}
\special{psfile=pp2pnpipl_pol.ps
hoffset=50 voffset=-10 hscale=60 vscale=60}
\caption{\it{Spin correlation parameters for the reaction
$pp \to pn\pi^+$. Same description of the curves as in
Fig.~\ref{csresults}. The experimental data are taken from
Ref.~\protect\cite{Sah99}.
}}
\label{Corr2}
\end{figure}

Results of this model for total cross sections for the reactions channels
$pp \rightarrow pp\pi^0$, $pp \rightarrow pn\pi^+$,
$pn \rightarrow pp\pi^-$, and $pp \rightarrow d\pi^+$ are displayed in
Fig. \ref{csresults}.
Results for
 analyzing powers
 were
presented in Refs.~\cite{Han98a,Han98b}. It was found that the model
yields a very good overall description of the data from the
threshold up to the $\Delta$ resonance region.
In fact, a nice quantitative agreement with basically all experimental
information (then available) was observed over a wide energy range.
Thus, this model is very well suited as a starting point for a detailed
analysis of the forthcoming spin-dependent observables of the reaction
$NN \rightarrow NN\pi$.

Looking at the results for the total cross sections one observes
that the role of the $\Delta$ is very different in the different
production channels: it dominates the $\pi^+$ production at
values of $\eta > 1$ but it appears much less relevant anywhere else.
As we will see in the following, this is no longer true for
polarization observables. Here the
$\Delta$ turns out to be very important in all
production channels in the whole energy range
if one wants to describe the spin observables as well.
This clearly demonstrates the power of polarization
observables to investigate resonances even if their effect
is too weak to show up in unpolarized observables!

Predictions for the spin correlation coefficient combinations
$A_\Sigma = A_{xx}+A_{yy}$, $A_\Delta = A_{xx}-A_{yy}$ and $A_{zz}$
are shown in Figs.~\ref{Corr1} (for $pp\rightarrow pp\pi^0$)
and \ref{Corr2} (for $pp\rightarrow pn\pi^+$). The polar integrals of
these observables are displayed in the left panels as a function of
$\eta$.
The other two panels contain the results for the different observables
as a function of the pion angle (middle
panel) and of the angle between the nucleons (right panel).
Data for these observables are currently being analyzed \cite{jan}.

Again the solid line shows the results for the full model, whereas the dashed line
are the results when all contributions involving the $\Delta$ are omitted.
As can be seen clearly from the figure, the inclusion of the isobar is
essential for all the polarization observables displayed even at energies
as low as $\eta = 0.4$!

Our model is not able to describe $A_{xx}-A_{yy}$ in the reaction $pp \to pp
\pi^0$ (c.f. Fig. \ref{Corr1}). The numerator of this observable is sensitive
to $Pp$ partial waves only \cite{Meyba} (Here capital letters denote the NN relative
 momentum whereas the small letters denote
the pion angular momentum with respect to the NN system).
Does figure \ref{Corr1} therefore tell
us that all but the $Pp$ piece is correctly reproduced by our model\footnote{This
was given as a possible explanation of this deviation in ref. \cite{Han98a}.}?
The answer is no. There is a possible other explanation, namely that it is not
the numerator of  $A_{xx}-A_{yy}$ but the denominator that causes the deviation.
Note, that our model underestimates the total cross section for $\pi^0$ production
by a factor of 2 at energies of $\eta = 1$ (c.f. Fig. \ref{csresults}).

As it was pointed out in ref. \cite{Mey99}, it is possible
to relate the integrated double polarization observables
to the spin cross sections (the cross sections one gets
for a given initial spin state). At moderate energies
these should be dominated by the lowest partial waves.
 A possible way to get better
insights could therefore be to look at the spin cross sections directly. However,
in their determination the total cross section enters besides the
double polarization observables.
The former, however, is badly known at energies around $\eta = 1$.
Better data is needed to
accurately constrain the spin cross sections in the energy region of interest.
Fortunately those data will soon be available \cite{Jozef}.

\section*{Vector meson production}

In this section we will briefly discuss, how the investigation of vector meson
production in nucleon--nucleon collisions can improve our knowledge on
the structure of the nucleon. In particular we will outline a method
on how to deduce information on the $NN\phi$ coupling constant from
data on differential cross sections recently measured by the DISTO collaboration \cite{Disto}.

This information can be regarded as complementary to
recent experiments on $\bar pp$ annihilation at rest, where
the unexpectedly large cross section
ratios $\sigma_{\bar pp \rightarrow \phi X} /
\sigma_{\bar pp \rightarrow \omega X}$
(cf. Ref.\cite{Ell95} for a compilation of data)
were interpreted by some
authors as a clear signal for an intrinsic
$\bar ss$ component in the nucleon\cite{Ell95}. However, in
an alternative approach based on two-step processes these data
were explained without introducing any strangeness
in the nucleon and any explicit violation of the OZI rule\cite{Loche}.

In this context $\phi$ production in nucleon-nucleon collisions is
of specific interest. Here one does not expect any significant
contributions from competing OZI-allowed two-step mechanisms.
Therefore cross section ratios
$\sigma_{pp\rightarrow pp \phi} / \sigma_{pp\rightarrow pp \omega}$
should provide a clear indication for a possible OZI violation and
the amount of hidden strangeness in the nucleon -- note, that
data recently presented by the DISTO collaboration indicate that this
ratio is about 8 times larger than the OZI estimate\cite{Disto}.
Based on a model calculation in this section it will be investigated, if the
observed enhancement over the OZI estimate in the cross section implies
a $g_{NN\phi}$ that is likewise enhanced and therefore
at variance with the OZI rule. It should be emphasized that the investigation
is by far not as ambitious as the study of pion production, for the
individual ingredients are much less known. The major goal of
this study is to extract constraints on $g_{NN\phi}$ from the $\omega$ to $\phi$ ratio.

As before
we also  describe the $pp \rightarrow pp\phi$ reaction within a
relativistic meson-exchange model, where
the transition amplitude is calculated in DWBA
in order to take the $NN$ final state interaction into account. (See
Ref.\cite{Nak1} for the details of the formalism.)
For the $NN$ interaction we employ the model Bonn B\cite{Mach89}.
The  effect of the ISI is accounted for via an appropriate adjustment of the
(phenomenological) form factors at the hadronic vertices.

In a previous study of the reaction $pp\rightarrow pp\omega$\cite{Nak1}
it was found that the dominant production mechanisms are the nucleonic and
$\omega\rho\pi$ mesonic currents, as depicted in Fig.~\ref{Fig1}.
Note, that in contrast to the pion production here we take the
Born diagram for the $\pi N \to \omega N$ T--matrix multiplied with a
phenomenological form--factor and not a microscopic T--matrix.
However, this is not a serious draw back for the current investigation,
since this formfactor is anyhow adjusted to the data.
The reason why it is still possible to extract important information
from the vector meson production is,
that the angular distribution of the produced vector
meson provides a unique and clear signature of the magnitude of these
currents, thus allowing one to disentangle these two reaction mechanisms.
Therefore, it is possible
to fix uniquely the magnitudes of the nucleonic and the meson-exchange
current by analyzing the angular
distribution of the $\phi$ meson measured by the DISTO
collaboration\cite{Disto}.
Furthermore, since the $NN\phi$ coupling constant enters only in the
nucleonic current it is possible to extract its value from such an
analysis. It is determined by the requirement of getting the
proper contribution of the nucleonic current needed to reproduce the
angular distributions.

\begin{figure}[t]
\vspace{4.0cm}
\special{psfile=Vprod.eps
hoffset=80 voffset=15 hscale=65 vscale=65}
\caption{$\phi$ and $\omega$-meson production currents included in
the present study: (a) nucleonic current (the diagrams with
the production after the meson exchange are also included),
 (b) meson exchange current.
$v = \omega, \phi$ and $M = \pi, \eta, \rho, \omega, \sigma, a_o$.}
\label{Fig1}
\end{figure}

The parameters of the model (coupling constants, cutoff masses of the
vertex form factors)
are mostly taken over from the employed $NN$ model. The $\phi\rho\pi$
coupling constant is obtained from the measured decay width of
$\phi \rightarrow \rho + \pi$. However, besides the $g_{NN\phi}$ that we
want to
extract from the analysis, there are still some more free parameters:
The cutoff mass of the $\phi\rho\pi$ vertex form factor
of the meson-exchange current, and the form factor and tensor- to
vector coupling constant ratio $\kappa_\phi \equiv f_{NN\phi}/g_{NN\phi}$
of the nucleonic current. It is possible to fix most of them by
performing a combined analysis of the $\omega$ and $\phi$ data
as well as other sources
 (cf. Ref.\cite{Nak2} for details).
The remaining unknown is the tensor coupling. Here
we assume that $\kappa_\phi =
\kappa_\omega$, as also suggested by SU(3) symmetry,
and $-0.5 \leq \kappa_\omega \leq 0.5$.

After the above considerations, we are now prepared to apply the model
to the reactions $pp\rightarrow pp\omega$ and $pp\rightarrow pp\phi$.
The angular distribution for $\phi$-meson production measured at
$T_{lab} = 2.85$ GeV is shown in Fig.~\ref{Fig2}.
We observe that the angular distribution is fairly flat. Recalling the
results we obtained for $\omega$ production\cite{Nak1} this
tells us that $\phi$-meson production should be almost entirely due to the
$\phi\rho\pi$ meson-exchange current. Only a very small contribution
of the nucleonic current is required if the angular
distribution drops at forward and backward angles, as indicated by
the data\footnote{Note, that since there are identical particles in
the initial state the angular distribution has to be symmetric around
90 degrees.}.

\begin{figure}[t]
\vspace{7.3cm}
\special{psfile=ang_dist_new.eps
hoffset=70 voffset=-15 hscale=50 vscale=50}
\caption{Angular distribution for the reaction $pp\rightarrow pp\phi$
at an incident energy of $T_{lab}=2.85$ GeV. The dashed-dotted curve corresponds
to the mesonic current contribution, the dashed curve to the nucleonic
current contribution. The solid curve is the total contribution. The experimental
data are from Ref.~\protect\cite{Disto}.}
\label{Fig2}
\end{figure}

Details about our strategy
for fixing the various parameters and information about the
data used in the analysis can be found
in Ref.\cite{Nak2}.
We get  a set of values
which range from
\begin{equation}
-1.4 \ \leq  \ g_{NN\phi}^{model} \ \leq \ -0.163 \ .
\end{equation}
Note,
that the spread is not an estimate of the theoretical uncertainty of the model
but reflects the remaining uncertainty within the model -- this will be
reduced sufficiently as soon as more experimental information is available.
Nevertheless, it is encouraging to see that the extracted values all lie within
fairly narrow bounds. This clearly indicates to us that the dependence on
the model parameters is not very strong, and that the magnitude of
$g_{NN\phi}$ is primarily determined by the experimental information used.

The values of $g_{NN\phi}$ obtained may be compared with those
resulting from SU(3) flavor symmetry considerations and imposition of the
OZI rule,
$$
g_{NN\phi}^{OZI} = - 3 g_{NN\rho} \sin(\alpha_v) \cong  -(0.60 \pm 0.15) \ ,
$$
where the factor $\sin(\alpha_v)$ is due to the deviation from the ideal
$\omega - \phi$ mixing. The numerical value is obtained using the values
of $g_{NN\rho} = 2.63 - 3.36$\cite{rhocoup} and $\alpha_v \cong 3.8^o$.
Comparing this value with the ones extracted from our model analysis, we
conclude that the preliminary data presently available can be described with
using a $NN\phi$ coupling constant that is compatible with the OZI rule.
This clearly indicates that a dynamical model is needed for drawing
any conclusion about the validity of the OZI rule.

\section*{Summary and Outlook}
The main points of the presentation can be summarized as
\begin{itemize}
\item it is important to take into account the nucleon--nucleon interaction
in the final as well as in the initial state
to get a quantitative understanding of meson production close to threshold,
\item the interesting energy regime for meson production is $0 \leq \eta \leq 1$, for
this regime includes both large and moderate momentum transfers,
\item polarization observables allow to detect resonances, even if their effect
is too small to show up in the total cross section. Note, that this statement
is true only if there is a spin dependent coupling of the resonance to the produced
meson.
\end{itemize}
It should be made clear that theory lacks far behind experiment. However,
the large amount of experimental information that is going to be available within
the next years will help us to get a better understanding of meson production
reactions and will therefore guide us to better models and insights into
the phenomenology of the nucleon--nucleon interaction.

\section*{Acknowledgments}

The author is very grateful to J. Haidenbauer for comments on the manuscript
and to G. A. Miller and U. van Kolck for useful discussions.
He also thanks K. Nakayama, J. W. Durso, J. Haidenbauer,
 J. Speth and O. Krehl for the collaboration that lead to the results presented here.