\section{Introduction}
\label{sec:fuv:intro}

Horizontal branch (HB) stars are core helium-burning stars that have evolved off the red giant branch (RGB). In globular clusters (GCs) they are so named because of their appearance on optical colour-magnitude diagrams, where they trace a horizontal path blueward of the RGB. These stars are thought to have similar core masses of $\sim$0.5$M_{\odot}$. Their location along the HB is therefore determined primarily by the mass and opacity of their thin stellar envelopes. To explain the observed HB stars in GCs, it is thought that substantial mass loss must occur during the star's RGB phase \citep[e.g.][]{Rood73}. The HB morphologies of different Galactic GCs are known to vary significantly. This suggests that the cluster properties have a significant effect on the evolution of stars onto the HB. The metallicity of a GC has long been proposed as the `first parameter' related to the morphology of a GC's HB, with metal poor clusters often having bluer HB stars than metal rich GCs \citep[e.g.][]{Sandage60,Dorman95}. However, the metallicity of a cluster alone is insufficient to explain the HB morphology in all Galactic GCs. Firstly, the clusters which host the bluest HBs actually have intermediate metallicities \citep[e.g.][]{OConnell99}. Also, some clusters are observed to have very different HB morphologies, but similar metallicities \citep[e.g. NGC~288 and NGC~362:][]{Bellazzini01}. This has led to the notorious search for a `second parameter' to describe HB morphology. Many different parameters have been proposed as a second parameter, all of which may play some role. These include the cluster's age \citep[e.g.][]{Demarque88,Lee94}, its helium abundance \citep[e.g.][]{Sweigart97,Lee05}, the stellar core rotation \citep[e.g.][]{Mengel76,Peterson85}, mass loss on the RGB \citep[e.g.][]{Rood73} and cluster core density \citep{Buonanno97}. Despite decades of research, the relative effects of these different parameters are yet to be fully understood \citep[for a discussion of proposed second parameters, see e.g. ][]{Fusi_Pecci97,Catelan09}.

The bluest HB stars observed in GCs are the `blue tail' or `extreme-HB' (EHB) stars. These are analogues to the subdwarf~B (sdB) stars observed in the Galactic field and are considered here to to be those stars with T$_{eff}>$20,000K \citep{Brown01}. These stars have become increasingly important as the leading contenders in explaining the `ultraviolet excess' (UVX) observed in elliptical galaxies \citep[e.g.][]{Code69,Dorman95,OConnell99}. Because of their high effective temperatures, they are relatively faint at optical wavelengths but make an important, often dominating, contribution to the FUV emission from old stellar populations. EHB stars are observed in both metal rich and metal poor clusters \citep[e.g.][]{DCruz96,Rich97} and it is likely that the metallicity may have little direct effect on the formation of these stars \citep[e.g.][]{Heber86,Dorman93,OConnell99}. Due to their very low envelope masses \citep[$\la$0.02M$_{\odot}$][]{Heber86}, it is likely that their formation and colour is more strongly related to the efficiency of mass loss mechanisms than the metallicity. At the hottest and optically faintest end of the EHB population, some Galactic clusters are observed to host a population of `blue hook' stars \citep[BHk;][]{Brown01,Dieball09}. These stars have similar UV colours to EHB stars, but are slightly fainter, lying just below the zero age horizontal branch population in UV colour magnitude diagrams (CMDs). Although BHk stars appear to be relatively rare, if a cluster hosts a significant population of these stars, then they are likely to make a contribution to the FUV luminosity of the cluster.

In this study, we investigate the effect of a cluster's core density on its HB population. A relationship has previously been proposed between core density and the extent of the HB population in the Milky Way's GCs \citep{Fusi_Pecci93,Buonanno97}. There are several reasons to suspect that dense stellar environments will enhance the formation of EHB stars. Firstly, mass loss may be enhanced through tidal interactions in close encounters between stars. Close encounters may also increase mixing in RGB stars and lead to helium enrichment, which will make the HB stars bluer \citep{Suda07}. It has also been proposed that EHB stars may be formed in binary systems either via white dwarf (WD)-WD mergers or from mass loss through either Roche lobe overflow or common envelope mechanisms \citep{Han02,Han03,Han07}. This theory is supported by observations of Galactic sdB stars, a large fraction of which are found to be binary systems \citep[e.g.][]{Maxted01,Reed04}. Interestingly, current observations suggest that the binary fraction among EHB stars may be much lower in GCs than in the field \citep{Moni_Bidin09}. This may be due to the relative contributions of different formation mechanisms evolving with time, with WD-WD mergers becoming increasingly import for old dense stellar populations, such as GCs \citep{Moni_Bidin08,Han08}. If EHB stars do have a binary origin, then a relationship between the core density of a cluster and the population of EHB stars it hosts may be expected. This is because of dynamical formation of these binary systems in the cores of these clusters.

Only a few of M31's GCs have CMDs which allow direct observation of their HB stars, and these are generally incomplete \citep[e.g.][]{Rich05}. Even in the Milky Way's GCs, studying EHB/BHk stars based on optical observations is difficult due to the severe crowding of sources and because these stars are optically faint. They are better identified in FUV CMDs, although such observations of Galactic clusters are severely affected by extinction. This means that only a few Galactic GCs have been studied from UV CMDs \citep[see e.g.][]{Brown01,Knigge02,Dieball05,Dieball10}. In the absence of direct stellar counts for M31's clusters, we consider their integrated FUV and NUV colours. This photometry provides a probe into the hot stellar populations in these GCs. Main sequence and RGB stars are too cool to have significant emission at these wavelengths and the clusters' integrated FUV and NUV magnitudes are likely to be dominated by objects hotter than $\sim$10,000K. If a cluster hosts a population of blue-HB (BHB), EHB or BHk stars, then these will likely dominate the cluster luminosity at these wavelengths. In the Milky Way, the integrated FUV-V colour is found to correlate with the HB morphology of the cluster \citep[e.g. fig. 3 of ][]{Catelan09}.

Recent \textit{GALEX} FUV and NUV observations of M31 have been presented by \citet{Rey05,Rey07}. In their work, the expected relationship between the clusters metallicity and both its FUV-V and NUV-V colour were shown. However, there is significant scatter in the FUV relationship. This confirms the need for a second parameter to describe the FUV luminosity of clusters in M31. \citet{Rey07} also identify three metal rich clusters which have a FUV excess, compared with both stellar models and other metal rich clusters in the galaxy. To help explain their observations, they consider the effects of age and helium abundance on the FUV luminosities of M31's GCs. Helium abundance is also used by \citet{Sohn06} and \citet{Kaviraj07a} to explain the FUV observations of M87's GCs. In this study, we consider correlations between the luminosity and core density of a GC and its FUV properties.

\section{M31 globular cluster data}
\label{sec:fuv:data}

In the following analysis, we use data from the recent catalogue of the M31's GCs from \citet{Peacock10}. This catalogue provides updated classifications, \textit{ugriz} and K-band photometry and structural parameters for objects listed in the Revised Bologna Catalogue of GCs in M31 \citep[RBC:][]{Galleti04}. We use this catalogue to identify 416 confirmed old clusters in the galaxy, of which 213 have structural parameters. The clusters with estimated structural parameters are spatially limited to the disk region of M31 (minus the very central region) and cover the entire GC luminosity function \citep[although they are thought to be relatively unreliable for the faintest clusters, as discussed in][]{Peacock09}. As the primary aim of this work is to investigate the effect of stellar density on the UV properties of the clusters, we restrict our primary sample to those clusters with structural parameters.

We combine this catalogue with the spectroscopic metallicities collated by \citet{Fan08} from the studies of \citet{Huchra91,Barmby00,Perrett02}. Where these metallicities are available from more than one study, the values from the larger and more recent catalogue of \citet{Perrett02} are chosen over \citet{Barmby00}. We reject the cluster B235 because of a large discrepancy ($\Delta$[Fe/H]=0.6) between the metallicities available from \citet{Barmby00} and \citet{Perrett02}. We also remove those clusters with large errors ($\sigma_{\rm{[Fe/H]}}>$0.6) on their metallicities (this removes the clusters B214, B229, B352 and BA11 from our analysis). This gives spectroscopically estimated metallicities for 219/416 confirmed old clusters in our catalogue. \citet{Fan08} also present reddening estimates for these clusters based on this spectroscopy and on photometry from the RBC. We use these estimates of $E(B-V)$ to deredden our photometry using the extinction curves of \citet{Cardelli89}. The relative extinction in the different bands were evaluated from these curves by \citet{Schlegel98} (for the \textit{ugriz} filters) and by \citet{Rey07} (for the \textit{GALEX} FUV, $R_{\rm{FUV}}$=8.16 and NUV, $R_{\rm{NUV}}$=8.90). These values are relatively uncertain. This is because the reddening curve in the UV is known to vary along different lines of sight in the Milky Way. It is also not guaranteed that M31 has the same extinction curve as the Milky Way, although \citet{Barmby00} found no evidence of significant differences in the optical bands. To minimise the effects of this on our UV colours, we limit our analysis to the clusters with relatively low extinction \citep[$E(B-V)<$0.16, consistent with the limit adopted by][]{Rey07}.

\begin{figure}
 \centering
 \includegraphics[height=84mm,angle=270]{figures/gi_g.eps}
 \caption{Optical colour magnitude diagram for all M31 GCs with data from \citet{Peacock10} and E(B-V)$<$0.16 from \citet{Fan08}. Small grey points indicate all GCs, while open black circles indicate those clusters detected in the \textit{GALEX} FUV observations.}
 \label{fig:mag_colour}
\end{figure}

\begin{figure*}
 \includegraphics[width=100mm,angle=270]{figures/uv_gk_ZF08_rho0_inter.eps}
 \caption{The ultraviolet colours of M31's GCs as a function of both \textit{g}-K colour (\textit{top}) and metallicity (\textit{bottom}). Colours are dereddened using the reddening estimates of \citet{Fan08}. The grey points/arrows show those clusters detected in the NUV but not the FUV. Filled circles indicate clusters with intermediate metallicities (-1.8$<$[Fe/H]$<$-1.2).}
 \label{fig:UV_Z}
\end{figure*}

\begin{figure*}
 \includegraphics[width=100mm,angle=270]{figures/uv_lumin_rho0_inter.eps}
 \caption{The ultraviolet colours of M31's GCs as a function of core density (\textit{top}) and K-band luminosity (\textit{bottom}). The bottom figures include 15 additional clusters, not covered by the survey of \citet{Peacock09} but, whose K-band luminosities are available from 2MASS observations (squares). As in figure \ref{fig:UV_Z}, filled symbols indicate intermediate metallicity clusters. }
 \label{fig:UV_lumin_rho0}
\end{figure*}

The FUV and NUV magnitudes of these clusters are taken from \citet{Rey05,Rey07}. These data provide integrated FUV and NUV magnitudes for 104 and 210 of the 416 confirmed old clusters in \citet{Peacock10}. For full details of these UV observations, we refer the reader to \citet{Rey05,Rey07}. The \textit{GALEX} FUV (1344-1786$\rm{\AA}$) and NUV (1771-2831$\rm{\AA}$) filters have an effective wavelength of 1528$\rm{\AA}$ and 2271$\rm{\AA}$, respectively \citep{Morrissey05}. The photometry is on the standard \textit{GALEX} ABmag photometric system, where:
\begin{equation}
 \label{eq:fuv_flux}
 \rm{FUV_{AB}} = -2.5\times\rm{log_{10}(F_{FUV}/1.40\times10^{-15}}) + 18.82
\end{equation}
\begin{equation}
 \label{eq:nuv_flux}
 \rm{NUV_{AB}} = -2.5\times\rm{log_{10}(F_{NUV}/2.06\times10^{-16}}) + 20.08
\end{equation}
Here, F is the flux of the source (erg $\rm{sec^{-1}cm^{-2}\AA^{-1}}$). The sample of clusters is magnitude limited to $\sim$22.5 in both the FUV and NUV. As highlighted by figure \ref{fig:mag_colour}, this does not correspond to a homogeneous optical magnitude limit. It can be seen that the blue (metal poor) clusters are detected across the GC luminosity function. However, very few of the red (metal rich) clusters are detected in the FUV. This is discussed by \citet{Rey07} and highlights the effect of metallicity on the HB properties, and hence FUV brightness, of a cluster. This plot also identifies the three metal rich clusters that are detected due to proposed excess FUV emission, as discussed by \citet{Rey07}.

Our final dataset contains 51 confirmed old clusters which have $E(B-V)<$0.16, NUV and ugriz photometry, estimates of their core densities and spectroscopic metallicities. Of these clusters 29 have FUV photometry.

\section{UV properties of M31 clusters}

\subsection{Metallicity relationships}

Figure \ref{fig:UV_Z} shows the FUV-\textit{g} and NUV-\textit{g} colours of M31's clusters as a function of metallicity and \textit{g}-K colour. It can be seen that the NUV-\textit{g} colour of a cluster is strongly correlated with its metallicity. Although relatively weak, the FUV-\textit{g} colour is also correlated with the metallicity. Considering the clusters with only upper limits available for their FUV magnitude confirms this relationship, demonstrating a lack of FUV bright, metal rich, clusters. These relationships are to be expected and highlight the effect of metallicity on the HB population in a cluster. This result was already identified by \citet{Rey07}. It is included here for completeness, but for a more detailed discussion, we refer the reader to this work.

The bottom right panel of figure \ref{fig:UV_Z} shows the metallicity as a function of FUV-NUV colour. It can be seen that the relationship is much less clear in this colour. A Spearman rank test actually shows that \textit{metal rich} clusters appear to be slightly bluer (with number of clusters, N=29, Spearman rank correlation coefficient, $\rho_{s}$=-0.39 and P-value for non-correlation, P=0.03). There is no obvious reason for such a relationship and it should be noted that selection effects, which bias us against FUV faint, metal rich clusters, could result in such a relationship. Also, this relationship is dominated by the two metal rich clusters in our sample. The FUV luminosity of these metal rich clusters may follow a different relationship to metal poor and intermediate metallicity clusters. However, it can be seen from these data that, unlike the UV-optical colour, there is no evidence for the FUV-NUV colour increasing with the metallicity of a cluster. This suggests that the effect of metallicity is much weaker on the UV colour of these clusters.

The two metal rich clusters in our sample have very blue FUV-NUV colours. This could result from selection effects, as these bias us against all but the FUV brightest metal rich clusters. However, the blue colours of these clusters could indicate that the FUV output of metal rich clusters (which generally have redder horizontal branches and few BHB stars) may be dominated by a population of EHB stars. EHB stars will have similar (or brighter) FUV luminosities to BHB stars. However, they have significantly bluer FUV-NUV colours \citep[as shown by the UV CMD of e.g. NGC~2808 and M80:][]{Brown01,Dieball10}.

\subsection{Stellar density$/$ Luminosity relationships}

As discussed in section \ref{sec:fuv:intro}, one may expect that the density of a cluster could effect its FUV luminosity. Figure \ref{fig:UV_lumin_rho0} shows the FUV-\textit{g}, NUV-\textit{g} and FUV-NUV colours of M31's GCs as a function of K-band luminosity and core density.

It can be seen from the top right panel of figure \ref{fig:UV_lumin_rho0} that the cluster core density appears to correlate with the FUV-NUV colour of the clusters. A Spearman rank test demonstrates a significant relationship between these parameters, with higher density clusters having bluer FUV-NUV colours (with N=29, $\rho_{s}$=-0.65 and P=$1.4\times10^{-4}$). Since EHB stars are relatively bright at FUV wavelengths, bluer UV colours could be produced by a larger fraction of EHB stars with respect to BHB stars. The observed relationship could thus be suggestive of a population of EHB stars, the formation of which is related to the core density of the cluster. This relationship is shown in more detail in figure \ref{fig:UV_rho0_limits}. In this figure we have highlighted the densest 50$\%$ of clusters in our sample. It can be seen that the densest group of clusters do appear to be offset to bluer colours. A Kolmogorov-Smirnov test between the high and low density clusters indicates that there is a low probability of 0.8$\%$ that they are drawn from the same population. The lack of dense, FUV faint, clusters can not be explained by selection effects. This is because the denser clusters are also the brighter clusters, and hence generally detected. It can be seen that there is significant scatter in the observed relationship. This is to be expected for two reasons. Firstly, if there is a relationship between the extent of the HB in a cluster and its core density, this is very unlikely to be the only parameter involved. Secondly, the UV colour is an indirect measure of the distribution of HB stars in a cluster. Because of this, other hot populations may weaken genuine trends.

\begin{figure}
 \includegraphics[width=84mm,angle=0]{figures/FUV_NUV_rho0_wHist.eps}
 \caption{\textit{Bottom panel:} The dereddened FUV-NUV colour of M31's GCs as a function of core density. Solid circles show the densest 50$\%$ of clusters. \textit{Top panel:} compares the FUV-NUV colour of all clusters in our sample (open) with the densest 50$\%$ of clusters in our sample (solid). It can be seen that the denser clusters appear to have bluer UV colours.}
 \label{fig:UV_rho0_limits}
\end{figure}

The bottom panels of figure \ref{fig:UV_lumin_rho0} show the UV colours of the clusters as a function of their total luminosity. To increase the number of clusters studied, we also include in this figure 15 additional clusters, that were not covered by the survey of \citet{Peacock09}. The K-band luminosities of these clusters were taken from the 2MASS survey, as presented by \citet{Galleti04}. A relationship is observed between the total K-band luminosity of a cluster and its UV colour, with the brightest clusters having blue UV colours (N=44, $\rho_{s}$=0.39 and P=0.009). This trend could be due to the same processes as the core density relationship, because of the strong correlation between the density and luminosity of a cluster \citep[as shown in e.g.][]{Peacock09}. However, the relationship may also be suggestive of self enrichment. The greater potential of the most massive clusters may result in them retaining more of the helium they produce. Since higher helium fractions are known to produce bluer HB stars, this could potentially produce the observed trend. Because of the relatively close relationship between the luminosity and density of a cluster, and the weakness of the observed correlations, it is very difficult to distinguish between the relative effects of these parameters.

As discussed previously, the metallicity of a cluster is expected (and observed) to effect its UV luminosity. For this reason, we also consider a sub-sample of intermediate metallicity clusters (with -1.8$<$[Fe/H]$<$-1.2; highlighted as filled symbols in figures \ref{fig:UV_Z} and \ref{fig:UV_lumin_rho0}). For the Galactic GCs, this method is found to highlight the effects of other parameters on HB stars \citep[since it reduces the influence of metallicity; e.g.][]{Recio_Blanco06}. In this sample, we observe the same core density and luminosity correlations found previously. However, these have a lower significance. This is likely due to the smaller sample size combined with the fact that some of the brightest and densest clusters (in which such an effect should be most evident) lie outside this metallicity range.

No significant trends are observed between the clusters K-band luminosities, or their core densities, and their FUV-\textit{g} or NUV-\textit{g} colours. This suggests that the UV-optical colours of these clusters are dominated by other effects, such as the cluster metallicity. However, since the effect of metallicity, on the FUV luminosity of a cluster, is strongest between the metal rich and metal poor clusters (as demonstrated by the low number of metal rich clusters detected in our data, see figure \ref{fig:mag_colour}), it is interesting to investigate relationships between only the metal poor clusters in our sample. We do this by removing the two clusters with [Fe/H]$>$-0.8. For this sample of metal poor/ intermediate metallicity clusters, a significant relationship is now identified between the FUV-\textit{g} colour and the cluster core density (with N=27, $\rho_{s}$=-0.60 and P=0.001). In contrast, the relationship between FUV-\textit{g} and metallicity is relatively weak for this sample (with N=27, $\rho_{s}$=0.32 and P=0.10). Given this observed correlation between the cluster core density and FUV-\textit{g}, we investigate whether this effect can help to explain some of the scatter that is observed in the metallicity relationship. To do this, we assume a planar relationship among a cluster's FUV-\textit{g} colour, its density and its metallicity. These data are fit to the function:
\begin{equation}
 \label{eq:fuv_g_fit}
 (\rm{FUV}-\textit{g})_{o} = \alpha\rm{[Fe/H]} + \beta\rm{log(\rho_{c})} + c
\end{equation}
From this we obtain $\alpha$=0.349 and $\beta$=-0.213. Using this relationship, we detrend our data for the effects of cluster core density and investigate whether this produces a tighter correlation between the clusters FUV-\textit{g} colour and its metallicity. This is plotted in figure \ref{fig:fit_FUVg_Z_rho0}. It can be seen that, accounting for the derived core density effect, the correlation between the metallicity and FUV-\textit{g} colour is stronger. This is confirmed by a Spearman rank test, which shows that the correlation with metallicity is stronger having removed the proposed core density effect (with N=27, $\rho_{s}$=0.57 and P=0.002). This suggests that effects due to stellar density may help to explain some of the scatter observed in the relationship between metallicity and FUV luminosity of the metal poor clusters.

\begin{figure}
 \centering
 \includegraphics[height=66mm,angle=270]{figures/fit_FUVg_Z_rho0.eps}
 \caption{FUV-\textit{g} colour as a function of core density (\textit{top}) and metallicity (\textit{middle}) for the metal-poor clusters (with [Fe/H]$<$-0.8). The \textit{bottom} panel shows the FUV-\textit{g} - metallicity relationship, detrended for the derived core density relationship. }
 \label{fig:fit_FUVg_Z_rho0}
\end{figure}

These data suggest that the cluster core density does have a significant effect on the FUV emission from these clusters. We caution that, while the above relationships are formally significant, they are based on a relatively small sample of clusters. However, we know of no reasons why selection effects should bias our results in such a way as to produce the observed relationship.

\subsection{Milky Way globular cluster data}

\begin{figure}
 \includegraphics[width=86mm,angle=270]{figures/MWGC_UV_newdata.eps}
 \caption{The UV colour of 24 of the Milky Way's GCs. UV colours are taken from \textit{ANS}, \textit{IUE} and \textit{OAO} observations as presented in \citet{Dorman95} and \citet{Sohn06}. Colours are transformed into the \textit{GALEX} filters using equation \ref{eq:AB_ST_mag}. In the top left plot, grey points indicate all Galactic GCs listed in the Harris catalogue. Open circles show all clusters with FUV and NUV photometry, filled circles highlight the clusters with intermediate metallicities. }
 \label{fig:UV_MWGC}
\end{figure}

\begin{table*}
 \footnotesize
 \centering
 \begin{minipage}{150mm}
  \caption{FUV luminosity of selected Galactic LMXBs \label{tab:FUV_MW_LMXBs}}
\end{minipage}
 \begin{minipage}{140mm}
  \centering
  \begin{tabular}{@{}lrrrr@{}}
  \hline
  \hline
  Name               & d [kpc]    & A$_{\rm{FUV}}$ & F$_{1400\rm{\AA}}$ [erg/cm$^{2}$/s/$\rm{\AA}$]  & L$_{1400\rm{\AA}}$ [10$^{32}$erg/s/$\rm{\AA}$] \\
  \hline
  Sco~X-1$^{1}$                &  2.8  & 2.430 & 2$\times10^{-13}$    &   18.0 \\
  Cyg~X-2$^{2}$                &  8.0  & 3.645 & 5$\times10^{-15}$    &   11.0 \\
  Her~X-1$^{3}$                &  5.5  & 0.405 & 7$\times10^{-14}$    &    3.7 \\
  AC211~(M15)$^{4}$            & 10.3  & 0.810 & 1.0$\times10^{-14}$  &    2.4 \\
  4U~1820-30~(NGC~6624)$^{5}$    & 7.9   & -     & 1.8$\times10^{-14}$  &    1.3 \\
  $[$ M31 GC-B022              & 780   & 0.324 & 4.0$\times10^{-17}$  &   40 $]$ \\
  $[$ M31 GC-B338              & 780   & 1.134 & 1.1$\times10^{-15}$  & 2000 $]$ \\
  $[$ M31 GC-B193              & 780   & 0.891 & 7.8$\times10^{-17}$  &  110 $]$ \\
  $[$ M31 GC-B225              & 780   & 0.810 & 4.0$\times10^{-16}$  &  550 $]$ \\
  \hline
  \end{tabular}\\
\end{minipage}
\begin{minipage}{150mm}
\footnotesize
The estimated distance to, extinction (assuming $R_{\rm{FUV}}$=8.1), FUV flux and (extinction corrected) FUV luminosity of selected Galactic LMXBs. For comparison, the integrated FUV luminosity of the the faintest and brightest M31 GCs detected by \citet{Rey07} are listed. Also listed are the FUV bright metal rich clusters in M31, which are known to host LMXBs. The Galactic LMXB data are taken from: $^{1}$ \citet{Vrtilek91}; $^{2}$ \citet{Vrtilek90}; $^{3}$ \citet{Boroson00}; $^{4}$ \citet{Dieball07}; $^{5}$ \citet{King93}.
\end{minipage}
\end{table*}

Given the relatively low number of clusters in our sample, and the observed (and expected) subtlety of any relationships, it is also interesting to consider the UV properties of the Milky Way's clusters. The distances, metallicities and structural parameters of the Milky Way GCs were taken from the Harris catalogue \citep{Harris96}. Ultraviolet observations of the Milky Way's GCs are complicated due to their large angular size, contamination (from foreground stars and background active galactic nuclei) and the high (and often variable) extinction in the direction of some of these clusters. This means that data are currently only available for a small sample of clusters. Integrated FUV and NUV photometry of 24 clusters were presented by \citet{Dorman95} using a combination of observations from the \textit{Astronomical Netherlands Satellite (ANS)}, \textit{Orbiting Astronomical Observatory (OAO)} and \textit{International Ultraviolet Explorer (IUE)}. These data were taken from table~1 of \citet[][for the \textit{IUE} data]{Dorman95} and table~6 of \citet[][which makes use of updated extinction estimates for these data]{Sohn06}. We exclude clusters whose photometry is considered unreliable by \citet{Dorman95} and the cluster NGC~6441, whose colour may be contaminated by a forground star. These data are on the STmag system and can be transformed to the \textit{GALEX} ABmag system via \citep{Sohn06}:

\begin{equation}
 \label{eq:AB_ST_mag}
 \rm{(FUV-NUV)_{AB} = (FUV-NUV)_{STmag}+0.854}
\end{equation}

Figure \ref{fig:UV_MWGC} shows the UV colour of the Milky Way's GCs as a function metallicity, luminosity, and core density. It can be seen that the FUV-NUV colour does not get redder with increasing metallicity (with N=24, $\rho_{s}$=-0.14, P=0.51). Interestingly, there appears to be a lack of FUV bright metal poor clusters, with intermediate metallicity clusters having the bluest colours. The lack of a correlation between FUV-NUV and metallicity is in agreement with the observations of M31's clusters and adds support to the idea that metallicity has less of an influence on the UV colour than it does on the UV-optical colour.

We also find no significant correlations between the FUV-NUV colour of the clusters and either their luminosity or core density. Limiting the sample to intermediate metallicity clusters (-2.0$<$[Fe/H]$<$-1.2, solid points in figure \ref{fig:UV_MWGC}), does identify a weak correlation between the FUV-NUV colour and core density. Unlike our M31 sample, a few of the Milky Way's low density clusters are found to have very blue FUV-NUV colours. We note that the bluest of these is $\omega$~Cen. This is the most massive cluster in our sample, so its blue colours may be suggestive of a mass effect. However, this cluster is also likely to have had a more complex star formation history than a typical cluster. As such, its UV colour may be atypical.

The relative weakness of correlations observed in the Milky Way's clusters (compared with those observed in our M31 data) may also be due to the reliability of their UV colours. The data available for the Milky Way's clusters is less homogeneous than our \textit{GALEX} observations and their UV colours may have increased errors due to the difficulty of obtaining accurate UV photometry in the Galaxy. Indeed a previous study by \citet{Djorgovski92} used a similar \textit{ANS} and \textit{IUE} dataset to demonstrate that ``more luminous and more concentrated clusters have bluer FUV colours". It is unclear why we do not identify this relationship in figure \ref{fig:UV_MWGC}, since the data are from similar surveys. However, there are several differences between the datasets. Firstly, \citet{Djorgovski92} have a larger, but possibly more contaminated dataset, because they include eight clusters that are excluded by \citet[][and hence our study]{Dorman95} because their photometry is potentially unreliable. Slightly different UV colours are used by \citet{Djorgovski92} and they include more colours from \textit{IUE} data (which are obtained through a smaller aperture than the \textit{ANS} data, hence excluding the outer regions of some clusters but possibly suffering from less contamination from non-cluster sources). Finally, both datasets are relatively inhomogeneous, being taken from multiple studies using multiple telescopes. We believe this discrepancy highlights the difficulty of studying FUV correlations in the Milky Way's clusters (particularly considering the small sample currently available). We do not attempt to draw any new conclusions from these Milky Way data, other than to suggest that it does not disagree with our conclusions based on M31's clusters. The need for higher quality UV photometry of the Galactic GCs is clear, and may become available from analysis of future/archived \textit{SWIFT/GALEX} observations (although such work will still have to contend with contamination and extinction issues).

Data for a larger sample of the Milky Way's clusters are available from studying optical CMDs. Using HST observations of 54 clusters, \citet{Recio_Blanco06} have identified that the more massive globular clusters do tend to host bluer HB stars (consistent with our findings in M31). This relationship is particularly strong for intermediate metallicity clusters (where the influence of metallicity is relatively small; see their figure 4). A relationship has also been identified between the length of a cluster's `blue tail' and its core density \citep{Fusi_Pecci93,Buonanno97}. Both of these findings are in agreement with our conclusions based on integrated UV observations of M31's clusters. Interestingly, such a core density relationship was \textit{not} found in the study of \citet{Recio_Blanco06}. The reason for this discrepancy is not clear. However, it possibly highlights the difficulty in accurately studying the blue extremes of a clusters HB population using optical CMDs.

\subsection{Low mass X-ray binaries and other FUV-bright sources}

It is generally assumed that, if a cluster hosts a significant population of BHB, EHB or BHk stars, then these will dominate the integrated FUV luminosity of the cluster. This is because most other stars in a cluster will be too cool to emit significantly in the FUV. However, there are other FUV bright sources in GCs. We consider the influence of these on the clusters integrated luminosity below.

Low mass X-ray binaries (LMXBs) are found in the cores of many GCs. These objects can have very high X-ray luminosities which can irradiate the systems accretion disks and$/$or donor stars and heat them to very high temperatures. The result of these large, hot, accretion disks is that LMXBs can be very bright in the FUV. Indeed, the LMXB 4U~1820-30 in the Galactic GC NGC~6624 dominates the total FUV luminosity of the cluster \citep{King93}. Even in the GC M15, which hosts a large population of FUV bright HB stars, the LMXB AC211 is the brightest FUV source in the cluster core \citep{Dieball07}. In table \ref{tab:FUV_MW_LMXBs} we list the FUV luminosity of these two Galactic GC LMXBs and three other Galactic LMXBs. The fluxes of these LMXBs are highly variable. This table therefore provides only an estimate of the highest FUV fluxes reached during these observations. It can be seen that the bright LMXBs Sco~X-1 and Cyg~X-2 are also very bright FUV sources. However, in the context of GCs, it should be noted that GCs are unlikely to host more than one or two LMXBs in outburst. This means that, for clusters with a significant BHB population, a single bright LMXB is unlikely to contribute more than $\sim10\%$ to the integrated FUV luminosity of the clusters. However, for metal rich GCs (which generally lack BHB stars), LMXBs could dominate the emission. This effect will be enhanced by that fact that metal rich clusters preferentially form LMXBs \citep[see e.g.][for the first clear evidence of this effect]{Kundu02}. It can be seen from table \ref{tab:FUV_MW_LMXBs} that a single LMXB is unlikely to explain the FUV bright metal rich clusters identified by \citet{Rey07}. However, it is interesting to note that two of the three metal rich clusters they identified do host LMXBs (B225 and B193). In particular, the cluster B225 hosts an LMXB with $L_{x}>10^{38}$ergs/s \citep{Trudolyubov04}. Such a source will be very bright in the FUV.

Single `UV bright' stars, such as post-AGB stars, can also contribute a significant fraction of the integrated luminosity of a cluster. However, clusters are also unlikely to host many of these sources due to their short lifetimes. It is thought that the contribution from such stars, to the integrated luminosity of a cluster, is unlikely to be more than 15$\%$ \citep{Moehler01}.

In addition to these objects, there are several other hot populations in GCs. Cataclysmic variables (CVs; white dwarfs accreting from Roche lobe overflowing companions) are relatively common in GCs. Because of their decreased bolometric luminosities, CVs are generally fainter than LMXBs in the FUV. However, they are more numerous and can reach FUV luminosities similar to a typical BHB star. Single white dwarf stars also emit strongly in the FUV. The average luminosity of both CVs and white dwarfs are thought to be too low to make a significant contribution to the integrated FUV luminosity of a cluster. The other notable FUV sources in GCs are blue straggler stars. These are very bright, blue objects and can also reach FUV luminosities similar to BHB stars, although blue stragglers this bright are likely to be rare. These stars are also unlikely to make a significant contribution to the integrated FUV luminosity of a cluster with a large population of blue HB stars. However, they may become increasingly important for clusters with redder HB populations.

\section{Conclusions}

We have considered the UV properties of M31's GCs. The previously identified relationships between the metallicity of a cluster and both its FUV-\textit{g} and NUV-\textit{g} colours are confirmed. Considering M31's metal poor/ intermediate metallicity clusters, we identify a relationship between the clusters core densities and their FUV-\textit{g} colours. This suggests that stellar density may influence the HB stars in these clusters. We demonstrate that the effects of stellar density on the properties of HB stars in a cluster may help to explain some of the scatter in the observed metallicity relationships. We also identify a weak anti-correlation between the FUV-NUV colour of a cluster and its metallicity. While this relationship is not very significant, it is in the opposite direction to the general relationship between metallicity and UV-optical colours. This suggests that the metallicity of a cluster may have little effect on the blue extremes of the HB distribution.

These data show a significant relationship between the FUV-NUV colours of these clusters and their core densities. We interpret this result in the context of a population of EHB stars, the production of which is enhanced in dense stellar environments. This may be due to enhanced mass loss or helium enhancement, in close encounters, or dynamical formation of tight binary systems (in which mass loss can occur, forming the EHB stars). We also find that more massive clusters tend to have bluer UV colours. This may be related to the same mechanisms suggested for the core density relationship, or it may point to an effect due to self enrichment. These trends are not confirmed in a smaller sample of integrated UV colours of the Milky Way's GCs. However, we note that this data is less homogeneous (and may suffer from more contamination) than our M31 GC data and that other studies of the Milky Way's clusters have suggested such trends (based on both UV and optical observations). It should also be noted that, if such relationships are present, they are expected to be relatively weak.

The relationship between cluster density and the population of EHB stars is potentially important in understanding the formation of these stars and the distribution of HB stars in clusters. Also, if a cluster's density does have a significant influence on its FUV emission, then it is important to consider this when using FUV observations to estimate the ages of globular clusters. Currently, such work does not take density effects into account.

We caution that, while these data suggest interesting correlations, the relationships observed are for a relatively small sample of M31's clusters. Obtaining new measurements of the structural parameters, metallicities and reddenings for more of M31's clusters would allow us to extend this work to a larger sample of clusters. Obtaining deeper UV data would also enable us to investigate similar relationships in the metal rich clusters in M31. However, there are no known reasons why the data currently available should bias us towards the observed relationships.

The FUV emission from other sources in a GC are also considered. These sources are unlikely to dominate the integrated FUV emission of a cluster if it contains a significant population of BHB, EHB or BHk stars (such as the FUV bright clusters considered in this study). However, other stellar populations, in particular LMXBs and blue stragglers, may contribute a significant fraction to the integrated FUV luminosity in clusters with only a red HB population.

\section*{Acknowledgements}

We would like to thank Steve Zepf and Arunav Kundu for very helpful discussions related to this work. We would also like to thank the reviewer of this paper for their helpful comments.

\bibliographystyle{mn2e}


\section*{Introduction}

Let $(M, \pi)$ be a Poisson manifold. A \textbf{symplectic
realization} of $(M, \pi)$ is a symplectic manifold $(S, \omega)$
together with a Poisson submersion
\[ \mu: (S, \omega) \longrightarrow (M, \pi).\]
Although the existence of symplectic realizations if a fundamental
result in Poisson geometry, the known proofs are rather involved.
Originally, the local result was proven in \cite{Wein} and a gluing
argument was provided in \cite{DWC}; the same procedure appears in
\cite{Kar}. The path approach to symplectic groupoids \cite{CaFe,
CrFe2} gives a different proof. Here we present a direct, global,
finite dimensional proof, based on the philosophy of
\textbf{contravariant geometry}: in Poisson geometry, the relevant
tangent directions come from the cotangent bundle $T^*M$ via the
bundle map
\[ \pi^{\sharp}: T^*M \longrightarrow TM\]
which is just $\pi$ converted into a linear map ($\beta(\pi^{\sharp}(\alpha))= \pi(\alpha, \beta)$). More on contravariant geometry can be found in the next section.
We will use a contravariant version of the notion of spray. In the following definition, for $t> 0$,
we denote by $m_t: T^*M\longrightarrow T^*M$ the fiberwise multiplication by $t$.

\begin{definition} A \textbf{Poisson spray} on the Poisson manifold $(M,\pi)$ is a vector field $\mathcal{V}_{\pi}$ on $T^*M$ satisfying the following two properties:
\begin{enumerate}
\item[(1)] $(dp)_{\xi}(\mathcal{V}_{\pi, \xi})= \pi^{\sharp}(\xi)$ for all $\xi\in T^*M$.
\item[(2)] $m_{t}^{*}(\mathcal{V}_{\pi}) = t\mathcal{V}_{\pi}$ for all $t> 0$.
\end{enumerate}
We denote by $\varphi_t$  the flow of $\mathcal{V}_{\pi}$.
\end{definition}

A short discussion on Poisson sprays-completely analogous the classical sprays \cite{Lang}- is given in the next section.
Condition (1) means that the integral curves of $\mathcal{V}_{\pi}$ are cotangent curves (see the next section);
it also appears in \cite{Wein2} under the name ``second order differential equation''.
Our main result is the following.

\begin{theorem} \label{main-theorem} Given the Poisson manifold $(M, \pi)$ and a contravariant spray $\mathcal{V}_{\pi}$,
there exists an open neighborhood $\mathcal{U}\subset T^*M$ of the
zero-section so that
\[ \omega:= \int_{0}^{1} (\varphi_t)^*\omega_{\textrm{can}} dt\]
is a symplectic structure on $\mathcal{U}$ and the canonical projection
$p: (\mathcal{U}, \omega)\longrightarrow (M, \pi)$
is a symplectic realization.
\end{theorem}

\begin{example} \rm \ When $M= U\subset \mathbb{R}^n$ open, denoting by $\pi_{i, j}$ the components of $\pi$, the simplest contravariant spray is
$\mathcal{V}_{\pi}(x, y)= \sum_{p, q} \pi_{p, q}(x) y_p\frac{\partial}{\partial x_q}$, where $x$ are the coordinates of $U$ and $(x, y)$ the induced coordinates
on $T^*U$. It is not difficult to see that the resulting $\omega$ coincides with the one constructed by A.Weinstein \cite{Wein}.
\end{example}

One may expect that the proof is ``just a computations''. Although
that is true in principle, the computation is more subtle then one
may believe. In particular, we will make use of the principle of
``contravariant geometry'' which is intrinsic to Poisson geometry.
The fact that the proof cannot be so trivial and hides some
interesting geometry behind was already observed in the local case
by A. Weinstein in \cite{Wein}: the notion of contravariant spray,
its existence, the formula for $\omega$ (giving a symplectic form on
an small enough $\mathcal{U}$)- they all make sense for any bivector
$\pi$, Poisson or not. But the fact that the push-down of (the
inverse of) $\omega$ is $\pi$ can only hold for Poisson bivectors.
Nowadays, with all the insight into symplectic groupoids, we can say
that we have the full geometric understanding of this theorem; in
particular, it can be derived from the path-approach to symplectic
groupoids of \cite{CaFe} and the resulting construction of local
symplectic groupoids \cite{CrFe2}. However, it is clearly worth
giving a more direct, global argument.

Let us start already with the first steps of the proof. Let's first
look at $\omega$ on vectors tangent to $T^*M$ at zero's $0_x\in
T^{*}_{x}M$ ($x\in M$). At such points one has a canonical
isomorphism $T_{0_x}(T^*M)\cong T_xM\oplus T_{x}^{*}M$ denoted
$v\mapsto (\overline{v}, \theta_v)$, and the canonical symplectic
form is
\begin{equation}\label{omega-can-0}
\omega_{\textrm{can}, 0_x}(v, w)= \langle \theta_{w},
\overline{v}\rangle -\langle \theta_{v}, \overline{w}\rangle.
\end{equation}
From the properties of $\mathcal{V}_{\pi}$ it follows that $\varphi_t(0_x)= 0_x$ for all $t$ and all $x$, hence $\varphi_{t}$ is well defined on a neighborhood of the zero-section, for all $t\in [0, 1]$.
From the same properties it also follows that
\[ (d\varphi_t)_{0_x} : T_{0_x}(T^*M)\longrightarrow T_{0_x}(T^*M) \]
is, in components,
\[ (\overline{v}, \theta_v)\mapsto (\overline{v}+ t\pi^{\sharp}\theta_{v}, \theta_{v}) .\]
From the definition of $\omega$ and the previous formula for $\omega_{\textrm{can}}$ we deduce that
\begin{equation}\label{omega-0}
\omega_{0_x}(v, w)= \langle \theta_{w}, \overline{v}\rangle- \langle
\theta_{v}, \overline{w}\rangle+ \pi(\theta_{v}, \theta_{w})
\end{equation}
for all $v, w\in T_{0_x}(T^*M)$. This implies that $\omega$ is nondegenerate at all zero's $0_x\in T^*M$. Hence we can find a neighborhood
$\mathcal{U}$ of the zero-section in $T^*M$ such that $\varphi_t$ is defined on $\mathcal{U}$ for all $t\in [0, 1]$ and $\omega|_{\mathcal{U}}$ is nondegenerate
(hence symplectic). Fixing such an $\mathcal{U}$, we still have to show that the map
\[ (dp)_{\xi}: T_{\xi}(\mathcal{U})\longrightarrow T_{p(\xi)}M \]
sends the bivector associated to $\omega$ to $\pi$, for all $\xi\in
\mathcal{U}$. The fact that this holds at all $\xi= 0_x$ follows
immediately from the previous expression for $\omega_{0_x}$. Our job
is to show that it holds at all $\xi$'s. Although it will not be
used in this paper, it is worth mentioning here Libermann's result
on symplectically complete foliation, concerning the following
question: given a symplectic manifold $(\mathcal{U}, \omega)$ and a
submersion $p: \mathcal{U}\longrightarrow M$, when can one push-down the
bivector associated to $\omega$ to a bivector on $M$? Libermann's
theorem (see e.g. Theorem 1.9.7 in \cite{DZ}) gives us the following
characterization: considering the involutive distribution
$\mathcal{F}(p)\subset T\mathcal{U}$ tangent to the fibers of $p$,
its symplectic orthogonal with respect to $\omega$,
$\mathcal{F}(p)^{\perp}\subset T\mathcal{U}$, must be involutive.
What happens in our case is the following:
\begin{equation}\label{main-lemma}
\mathcal{F}(p)^{\perp}= \mathcal{F}(p_1),
\end{equation}
where $\mathcal{F}(p_1)$ is the (involutive!) distribution tangent to the fibers of $p_1:= p\circ \varphi_1: \mathcal{U}\longrightarrow M$.
This will be proven in the last section. However, it turns out that the ingredients needed to prove this equality can be used
to show directly that $p$ is a Poisson map, without having
to appeal to Libermann's result.

\section{Contravariant geometry}

As we have already mentioned, the basic idea of contravariant geometry in Poisson geometry is that of replacing the tangent bundle $TM$ by $T^*M$. The two are related
by the bundle map $\pi^{\sharp}$.  But the main structure that makes everything work is the presence of a Lie bracket $[-, -]_{\pi}$ on $\Gamma(T^*M)$, which is the contravariant analogue of the Lie bracket on vector fields
(the two brackets being related via $\pi^{\sharp}$). It is uniquely determined by the condition
\[ [df, dg]_{\pi}=d\{f, g\}\]
and the Leibniz identity
\[ [\alpha, f\beta]_{\pi}= f[\alpha, \beta]_{\pi}+ L_{\alpha}(f) \beta \]
for all $\alpha, \beta\in \Omega^1(M)$, where $L_{\alpha}:= L_{\pi^{\sharp}\alpha}$ is the Lie derivative along the ordinary vector field $\pi^{\sharp}(\alpha)$
associated to $\alpha$. In other words, contravariant geometry is the geometry associated to the Lie algebroid $(T^*M, \pi^{\sharp}, [-, -]_{\pi})$.
Here are some examples of notions that are contravariant to the usual ones (see e.g. \cite{CrFe2, Fer}).

\vspace*{.1in}
A \textbf{contravariant connection} on a vector bundle $E$ over $M$ is a a bilinear map
\[\nabla: \Gamma(T^*M)\times \Gamma(E), \ \ (\alpha, s)\mapsto \nabla_{\alpha}(s)\]
satisfying
\[ \nabla_{f\alpha}(s)= f\nabla_{\alpha}(s), \ \ \nabla_{\alpha}(fs)= f\nabla_{\alpha}(s)+ L_{\alpha}(f) s\]
for all $f\in C^{\infty}(M)$, $\alpha\in \Gamma(T^*M)$, $s\in \Gamma(E)$. The standard operations with connections (duals, tensor products, etc) have an obvious contravariant
version.

\vspace*{.1in}
A \textbf{cotangent path} (or contravariant path) is a path $a: [0, 1]\longrightarrow T^*M$ sitting above some path $\gamma: [0, 1]\longrightarrow M$, such that
\[ \pi^{\sharp}(a(t))= \frac{d\gamma}{dt}(t).\]
Intuitively, the cotangent path is the pair $(a, \gamma)$ where $\gamma$ is a standard path and the role of $a$ is to encode ``the contravariant derivative of $\gamma$''.
The previous equation says that the contravariant derivative is related to the classical one via $\pi^{\sharp}$.

\vspace*{.1in} Given a contravariant connection $\nabla$ on a vector
bundle $E$, one has a well-defined notion of \textbf{derivative of
sections along cotangent paths}: given a cotangent path $(a,
\gamma)$ and a path $u: [0, 1]\longrightarrow E$ sitting above $\gamma$,
$\nabla_{a}(u)$ is a new path in $E$ sitting above $\gamma$. Writing
$u(t)= s_t(\gamma(t))$ for some time dependent section of $E$,
\[\nabla_{a}(u)= \nabla_{a}(s_t)(x)+ \frac{d  s_t}{dt}(x), \ \ \ \textrm{at}\ x= \gamma(t).\]

\vspace*{.1in}
Given a contravariant connection $\nabla$ on $T^*M$, the \textbf{contravariant torsion} of $\nabla$ is the tensor $T_{\nabla}$ defined by
\[ T_{\nabla}(\alpha, \beta)= \nabla_{\alpha}(\beta)- \nabla_{\beta}(\alpha)- [\alpha, \beta]_{\pi}.\]

\vspace*{.1in} Given a metric $g$ on $T^*M$, one has an associated
\textbf{contravariant Levi-Civita connection}- the unique
contravariant metric connection $\nabla^{g}$ on $T^*M$ whose
contravariant torsion vanishes. The corresponding
\textbf{contravariant geodesics} are defined as usual, as the
(cotangent) curves $a$ satisfying $\nabla_{a}a= 0$. They are the
integral curves of a vector field $\mathcal{V}_{\pi}^{g}$ on $T^*M$,
called the \textbf{contravariant geodesic vector field}. In local
coordinates $(x, y)$ (where $x$ are the coordinates in $M$ and $y$
on the fiber),
\[ \mathcal{V}_{\pi}^g(x, y)= \sum_{p, q} \pi_{p, q}(x)y_p \frac{\partial}{\partial x_q}- \sum \Gamma_{p, q}^{r}(x) y_py_q \frac{\partial}{\partial y_r},\]
where $\Gamma_{p, q}^r$ are the coefficients in
$\nabla_{dx_p}(dx_q)=\sum \Gamma_{p,q}^{r} dx_r$. Geodesics and the
geodesic vector field are actually defined for any contravariant
connection $\nabla$ on $T^*M$, not necessarily of metric type. For
instance, any classical connection $\nabla$ on $T^*M$ induces a
contravariant connection with $\nabla_{\alpha}:=
\nabla_{\pi^{\sharp}\alpha}$ which, in general, is not of metric
type.

Back to our problem, the existence of contravariant sprays is now clear:

\begin{lemma}
Any contravariant geodesic vector field is a contravariant spray.
\end{lemma}

Recall also that (cf. e.g. \cite{CrFe2}) any classical connection $\nabla$ induces two contravariant connections, one on $TM$ and one on $T^*M$, both denoted by
$\overline{\nabla}$:
\[ \overline{\nabla}_{\alpha}(V)= \pi^{\sharp}\nabla_{V}(\alpha)+ [\pi^{\sharp}(\alpha), V], \ \overline{\nabla}_{\alpha}(\beta)= \nabla_{\pi^{\sharp}\beta}(\alpha)+ [\alpha, \beta]_{\pi}.\]
The two are related by the following formula, which follows
immediately from the fact that $\pi^{\sharp}$ is a Lie algebra map
from $(\Omega^1(M), [-, -]_{\pi})$ to the Lie algebra of vector
fields. Note also that this (and its consequences later on) is the
only place where we use that $\pi$ is Poisson.

\begin{lemma} For any classical connection $\nabla$,
\begin{equation}\label{EQ4}
\overline{\nabla}_{\alpha}(\pi^{\sharp}(\beta))=\pi^{\sharp}(\overline{\nabla}_{\alpha}(\beta)).
\end{equation}
\end{lemma}

In the next section we will be using a $\nabla$ which is torsion-free; this condition simplifies the computations
because of the following lemma.

\begin{lemma}\label{ConjConn2} If $\nabla$ is torsion-free then,
for any cotangent path $a$ with base path $\gamma$, for any smooth path $\theta$ in $T^*M$ above $\gamma$ and any smooth path $v$ in $TM$ above
$\gamma$,
\[\langle\overline{\nabla}_{a}(\theta),v\rangle+\langle \theta,\overline{\nabla}_{a}(v)\rangle=\frac{d}{dt}\langle \theta, v\rangle.\]
\end{lemma}

\begin{proof}
Choose a  time-dependent 1-form $A= A(t, x)$ such that $a(t)= A(t,
\gamma(t))$ and similarly a time-dependent 1-form $\Theta$
corresponding to $\theta$ and a time-dependent vector field
corresponding to $v$. Applying the definition of the derivatives
$\overline{\nabla}_{a}$ along cotangent paths and then the
definition of $\overline{\nabla}$ we find that the left hand side at
time $t$ coincides with the following expression on $(t, x)$
evaluated at $x= \gamma(t)$:
\begin{equation}\label{inter}
\langle \nabla_{\pi^{\sharp}\Theta}(A)+ [A, \Theta]_{\pi}+
\frac{d\Theta}{dt}, V\rangle+ \langle \Theta,
\pi^{\sharp}\nabla_{V}(A)+ [\pi^{\sharp}A, V]+ \frac{dV}{dt}\rangle
.
\end{equation}
For the two terms involving $\nabla$ we find
\begin{eqnarray*}
&& \langle \nabla_{\pi^{\sharp}\Theta}(A), V\rangle + \langle \Theta, \pi^{\sharp}\nabla_{V}(A)\rangle =\\
&&\phantom{12}= \langle \nabla_{\pi^{\sharp}\Theta}(A), V\rangle - \langle \nabla_{V}(A), \pi^{\sharp}\Theta \rangle =\\
&&\phantom{12}= L_{\pi^{\sharp}\Theta}\langle A, V \rangle- \langle A, \nabla_{\pi^{\sharp}\Theta}(V) \rangle- L_{V}\langle A, \pi^{\sharp}\Theta \rangle + \langle A, \nabla_{V}(\pi^{\sharp}\Theta)\rangle =\\
&&\phantom{12}= L_{\pi^{\sharp}\Theta}\langle A, V \rangle-
L_{V}\langle A, \pi^{\sharp}\Theta \rangle+ \langle A, [V,
\pi^{\sharp}\Theta]\rangle,
\end{eqnarray*}
where we have used the antisymmetry of $\pi$, then we passed from $\nabla$ on $T^*M$ to its dual on $TM$ and then
we used that $\nabla$ is torsion-free. For the term in (\ref{inter}) containing $[-, -]_{\pi}$, using the definition of this bracket we find
\begin{eqnarray*}
&&  \langle [A, \Theta]_{\pi}, V\rangle= \langle L_{\pi^{\sharp}A}(\Theta)-L_{\pi^{\sharp}\Theta}(A)- d\pi(A, \Theta), V\rangle = \\
&&\phantom{12}= L_{\pi^{\sharp}A}\langle \Theta, V\rangle- \langle
\Theta, [\pi^{\sharp}A, V]\rangle - L_{\pi^{\sharp}\Theta}\langle A,
V\rangle + \langle A, [\pi^{\sharp}\Theta, V]\rangle - L_V(\pi(A,
\Theta)).
\end{eqnarray*}
Plugging the last two expressions into (\ref{inter}) we find
\[ \langle \frac{d\Theta}{dt}, V\rangle + \langle \Theta, \frac{dV}{dt}\rangle+ L_{\pi^{\sharp}A}\langle \Theta, V\rangle .\]
As an expression on $(t, x)$, when evaluated at $x=\gamma(t)$, since
$\pi^{\sharp}A= \frac{d\gamma}{dt}$, we find precisely the right
hand side of the expression from the statement.
\end{proof}

\section{A different formula for $\omega$}

In this section we give another description of $\omega$. The
resulting formula is a generalization of the formula (\ref{omega-0})
from zero's $0_x$ to arbitrary $\xi$'s in $T^*M$. It will depend on
a connection $\nabla$ on $TM$ which is used in order to handle
tangent vectors to $T^*M$. Hence, from now on, we fix such a
connection which we assume to be torsion free. With respect to
$\nabla$, any tangent vector $v\in T_{\xi}(T^*M)$ is determined by
the tangent vector induced on $M$ and by its vertical component
\[ \overline{v}= (dp)_{\xi}(v)\in T_{p(\xi)}M , \ \ \theta_v= v- \textrm{hor}_{\xi}(\overline{v})\in T_{p(\xi)}^{*}M .\]
Of course, when $\xi= 0_x$, these coincide with the components mentioned in the introduction. The fact that
$\nabla$ is torsion-free ensures the following generalization of the formula (\ref{omega-0}) for $\omega_{\textrm{can}}$ at arbitrary $\xi$'s.

\begin{lemma}\label{torsion-free1} If $\nabla$ is torsion-free then, for any $v, w\in T_{\xi}(T^*M)$,
\begin{equation}\label{Omega_can}
\omega_{\mathrm{can}}(v, w)=\langle \theta_w, \overline{v}\rangle-\langle
\theta_v,\overline{w}\rangle,
\end{equation}
\end{lemma}

\begin{proof}  Since $\nabla$ is torsion free, it follows that the associated horizontal distribution $H\subset T(T^*M)$ is Lagrangian with respect to $\omega_{\textrm{can}}$ and then
the formula follows.
\end{proof}

To establish the generalization of (\ref{omega-0}) to arbitrary
$\xi$'s, we start with a tangent vector
\[ v_0\in T_{\xi}(T^*M),\]
with $\xi\in T^*M$ fixed with the property that $\varphi_t(\xi)$ is defined up to $t= 1$. Consider
\[ a: [0, 1]\longrightarrow T^*M, \ a(t)= \varphi_t(\xi) \]
which, from the properties of $\mathcal{V}_{\pi}$, is a cotangent
path; we denote by $\gamma= p\circ a$ its base path. Pushing $v$ by
$\varphi_t$ we obtain a path
\begin{equation}\label{v-t}
t\mapsto v_t:= (\varphi_t)_*(v_0)\in T_{a(t)}(T^*M).
\end{equation}
Taking the components with respect to $\nabla$, we obtain two paths above $\gamma$, one in $TM$ and one in $T^*M$
\[ t\mapsto \overline{v}_t\in T_{\gamma(t)}M, \  \ t\mapsto \theta_{v_t}\in T_{\gamma(t)}^{*}M .\]
We denote these paths by $\overline{v}$ and $\theta_v$. They are related in the following way:

\begin{lemma} $\overline{\nabla}_a \overline{v}=\pi^{\sharp}\theta_{v}$.
\end{lemma}

\begin{proof} We start with one remark on derivatives along vector fields. For any tangent vector $V$ to $\mathcal{U}$ along $a$, $t\mapsto V(t)\in T_{a(t)}\mathcal{U}$,
one has the Lie derivative of $V$ along $\mathcal{V}_{\pi}$, again a tangent vector along $a$, defined by
\[ L_{\mathcal{V}_{\pi}}(V)(t)= \frac{d}{ds}|_{s= 0}  (d\varphi_{-s})_{a(s+t)}(V(s+t)) \in T_{a(t)}\mathcal{U}.\]
We have the following two remarks:
\begin{enumerate}
\item For vertical $V$'s, i.e. coming from a 1-form $\theta$  on $M$ along $\gamma$
$dp(L_{\mathcal{V}_{\pi}}(V))= -\pi^{\sharp}(\theta)$. This follows immediately from the
first property of the spray (e.g. by a local computation).
\item For horizontal $V$'s, $(dp)(L_{\mathcal{V}_{\pi}}(V))= \overline{\nabla}_{a}(\overline{V})$, where $\overline{V}= (dp)(V)$ is a tangent vector to $M$ along $\gamma$. To check this, one may assume that
$V$ is a global vertical vector field on $M$ and one has to show that $(dp)_{\eta}(L_{\mathcal{V}_{\pi}}(V))= \overline{\nabla}_{\eta}(\overline{V})$ for all $\eta\in T^*\mathcal{U}$. Again, this follows immediately by a local computation.
\end{enumerate}
Hence, for an arbitrary $V$ (along $a$), using its components $(\overline{V}, \theta_V)$,
\[ dp(L_{\mathcal{V}_{\pi}}(V)) = -\pi^{\sharp}(\theta_V) + \overline{\nabla}_{a}(\overline{V}).\]
Finally, remark that for our $v$ from (\ref{v-t}), $L_{\mathcal{V}_{\pi}}(v)= 0$.
\end{proof}

We have the following version of (\ref{omega-0}) at arbitrary $\xi$'s in $\mathcal{U}$- a small enough neighborhood
of the zero-section in $T^*M$ on which $\omega$ is well-defined.

\begin{lemma}\label{frmla} Let $\xi\in \mathcal{U}$, $v_0, w_0 \in T_{\xi}(T^*M)$.
Let $v=v_t$ as before and let $\widetilde{\theta}_{v}$ be a path in $T^*M$, solution of the differential equation
\begin{equation}\label{diff-eq}
\nabla_{a}(\widetilde{\theta}_v)= \theta_v .
\end{equation}
Similarly, consider $w= w_t$ and $\widetilde{\theta}_{w}$ corresponding to $w_0$. Then
\[ \omega(v_0, w_0)= (\langle \widetilde{\theta}_w, \overline{v}\rangle - \langle \widetilde{\theta}_v, \overline{w}\rangle- \pi(\widetilde{\theta}_{v},\widetilde{\theta}_w) )|_{0}^{1} .\]
\end{lemma}

\begin{proof}
Since $\nabla$ is torsion-free, Lemma \ref{torsion-free1} implies that
\[ \omega(v_0, w_0)= \int_{0}^{1}(\langle \theta_w, \overline{v}\rangle - \langle \theta_v, \overline{w}\rangle ) dt.\]
Hence it suffices to show that
\[ \langle \theta_w, \overline{v}\rangle - \langle \theta_v, \overline{w}\rangle =
\frac{d}{dt}  (\langle \widetilde{\theta}_w, \overline{v}\rangle - \langle \widetilde{\theta}_v, \overline{w}\rangle- \pi(\widetilde{\theta}_{v},\widetilde{\theta}_w) ).\]
We start from the left hand side, in which we plug in $\theta_v= \nabla_{a}(\widetilde{\theta}_v)$ and the similar formula for $w$, followed by the use of Lemma \ref{ConjConn2}, then the previous lemma, then again the defining formula
for $\widetilde{\theta}_v$; we obtain
\begin{eqnarray*}
&& \frac{d}{dt}(\langle \widetilde{\theta}_w, \overline{v}\rangle - \langle \widetilde{\theta}_v, \overline{w}\rangle)- \langle \widetilde{\theta}_{w}, \overline{\nabla}_{a}(\overline{v})\rangle + \langle \widetilde{\theta}_{v}, \overline{\nabla}_{a}(\overline{w})\rangle = \\
&&\phantom{12}= \frac{d}{dt}(\langle \widetilde{\theta}_w, \overline{v}\rangle - \langle \widetilde{\theta}_v, \overline{w}\rangle)- \langle \widetilde{\theta}_{w}, \pi^{\sharp}(\theta_v)\rangle + \langle \widetilde{\theta}_{v}, \pi^{\sharp}(\theta_w)\rangle= \\
&&\phantom{12}= \frac{d}{dt}(\langle \widetilde{\theta}_w, \overline{v}\rangle - \langle \widetilde{\theta}_v, \overline{w}\rangle)- \langle \widetilde{\theta}_{w}, \pi^{\sharp}(\nabla_{a}(\widetilde{\theta}_v))\rangle + \langle \widetilde{\theta}_{v}, \pi^{\sharp}(\nabla_{a}(\widetilde{\theta}_w))\rangle .
\end{eqnarray*}
For the expression involving the last two term, using
$\pi^{\sharp}(\nabla_{a}(\widetilde{\theta}_v))=
\nabla_{a}(\pi^{\sharp}(\widetilde{\theta}_v))$ and the antisymmetry
of $\pi^{\sharp}$, we find
\[ - \langle \widetilde{\theta}_{w}, \nabla_{a}(\pi^{\sharp}(\widetilde{\theta}_v))\rangle - \langle \overline{\nabla}_a(\widetilde{\theta}_w), \pi^{\sharp}(\widetilde{\theta}_v) \rangle ,\]
which, by Lemma \ref{ConjConn2} again, equals to $-\frac{d}{dt}
\langle \widetilde{\theta}_w,
\pi^{\sharp}(\widetilde{\theta}_v)\rangle= -\frac{d}{dt}
\pi(\widetilde{\theta}_{v},\widetilde{\theta}_w)$. Plugging in the
previous formula, the desired equation follows.
\end{proof}

\section{The proof of the theorem}

We now return to the proof of Theorem \ref{main-theorem}. We start with the proof of the equality
(\ref{main-lemma}) from the introduction. By a dimension counting, it suffices to prove the reverse inclusion. Fix $\xi\in \mathcal{U}$.
We have to show that $\omega(v_0, w_0)= 0$ for all
\begin{equation}\label{cond}
v_0\in \mathcal{F}(p)_{\xi}, \ w_0\in \mathcal{F}(p_1)_{\xi}.
\end{equation}
These conditions are equivalent to $\overline{v}(0)= 0$, $\overline{w}(1)= 0$, where we use the notations from the previous section.
Remark that (\ref{diff-eq}), as an equation on  $\widetilde{\theta}_{v}$, is a linear ordinary differential equation; hence it has solutions
defined for all $t\in [0, 1]$, satisfying any given initial (or final) condition. Hence one may arrange that $\widetilde{\theta}_{v}(0)= 0$, $\widetilde{\theta}_{w}(1)=0$.
The formula from Lemma \ref{frmla} immediately implies that $\omega(v_0, w_0)= 0$.

\vspace*{.2in}

Finally, we show that $p$ is a Poisson map. We have to show that, for $\xi\in \mathcal{U}$ arbitrary,
$\theta\in T_{x}^{*}M$  ($x= p(\xi)$), the unique $v_0\in T_{\xi}\mathcal{U}$ satisfying
\begin{equation}\label{Poisson_map}
 p^{*}(\theta)_{\xi}= i_{v_0}(\omega)
\end{equation}
also satisfies $(dp)_{\xi}(v_0)= \pi^{\sharp}(\theta)$. From the
previous formula it immediately follows that $v_0$ is in
$\mathcal{F}(p)^{\perp}$, hence in $\mathcal{F}(p_1)$, hence
$\overline{v}(1)= 0$, where we start using the notations from Lemma
\ref{frmla}. Next, we evaluate (\ref{Poisson_map}) on an arbitrary
$w_0\in T_{\xi}\mathcal{U}$. We also use the formula for $\omega$
from Lemma \ref{frmla}, where $\widetilde{\theta}_v$ and
$\widetilde{\theta}_w$ are chosen so that $\widetilde{\theta}_v(1)=
0$ and $\widetilde{\theta}_w(0)= \eta\in T_{x}^{*}M$ is arbitrary.
We find:
\[ \theta(\overline{w}_0)= \langle \widetilde{\theta}_v(0), \overline{w}_0\rangle + \langle \eta, \pi^{\sharp}\widetilde{\theta}_v(0)- \overline{v}_0\rangle\]
Since this holds for all $w_0$ and all $\eta$, we deduce that
$\theta= \widetilde{\theta}_v(0)$,
$\pi^{\sharp}\widetilde{\theta}_v(0)= \overline{v}_0$. Hence
$\pi^{\sharp}(\theta)= \overline{v}_0= (dp)_{\xi}(v_0)$.

\section{Some remarks}

Here are some remarks on possible variations. First of all, regarding the notion of contravariant spray, the first condition means that, locally, $\mathcal{V}_{\pi}$ is
of type
\[ \mathcal{V}_{\pi}(x, y)= \sum_{p, q} \pi_{p, q}(x, y)y_p\frac{\partial}{\partial x_q}+ \sum_i \gamma^i(x, y) \frac{\partial}{\partial y_i}.\]
The second condition means that each $\gamma^i(x, y)$ is of type $\sum_{j, k} \gamma^{i}_{j, k}(x) y_j y_k$.
While the first condition has been heavily used in the paper, the second one was
only used to ensure that $\omega$ is well-defined and non-degenerate at elements $0_x\in T^{*}_{x}M$.

Another remark is that one can show that $\mathcal{U}$ can be made into a local symplectic groupoid, with source map $p$ and target map $p_1$;
see also \cite{Kar}.

Let us also point out where we used that $\pi$ is Poisson: it is
only for the compatibility relation (\ref{EQ4}) which, in turn, was
only used at the end of the proof of the Lemma \ref{frmla}. However,
it is easy to keep track of the extra-terms that show up for general
bivectors $\pi$: at the right hand side of (\ref{EQ4}) one has to
add the term $i_{\alpha\wedge \beta}(\chi_{\pi})$ where $\chi_{\pi}=
[\pi, \pi]$, while to the equation from Lemma \ref{frmla} the term
$\int_{0}^{1} \chi_{\pi}(a, \widetilde{\theta}_v,
\widetilde{\theta}_w) dt$. This is useful e.g. for handling various
twisted versions. E.g., for a  $\sigma$-twisted bivector $\pi$ on
$M$ in the sense of \cite{WeinSev} (i.e. satisfying $[\pi, \pi]=
\pi^{\sharp}(\sigma)$ where $\sigma$ is a given closed 3-form
$\sigma$ on $M$), the interesting (twisted symplectic) 2-form on
$\mathcal{U}$ is the previously defined $\omega$ to which we add the
new two-form $\omega_{\sigma}$ given by (compare with \cite{CaXu}):
\[ \omega_{\sigma} = \int_{0}^{1} \varphi_{t}^{*}(i_{\mathcal{V}_{\pi}}p^*(\sigma)) dt.\]

\bibliographystyle{amsplain}
\def{}


\section{Introduction}

The incidence of extrasolar planets is known to be a function of
both stellar mass and metallicity. If the present-day metallicity of
the star reflects that of its circumstellar material at early times,
then metal-rich discs would be expected to build planetary cores
more readily by grain coagulation. The stellar mass may enter if
more massive stars have higher-mass discs, as well as in factors
such as shorter dynamical times for grains to collide. In the
simplest approach, the mass of solids in a disc, i.e. the total mass
M multiplied by the solid fraction Z, may be the relevant threshold
for a successful outcome in forming planets or planetesimals
\citep{wyatt07,greaves07}. These two quantities M,Z can be readily
estimated using proxies. Millimetre-wavelength emission from dust is
rather optically thin \citep{andrews07a} and the derived disc masses
in present-day star formation regions are canonically multiplied by
a factor $\approx 100$ to include the gas content inherited from the
interstellar medium. These masses are proxies in the sense that
material already converted into larger bodies (from boulders up to
planets) has little emitting area, and so potentially only some of
the orbiting material is represented. A good proxy for the
refractory proportion of material in discs of non-solar metallicity
is the logarithmic [Fe/H] of the host star.

\citet{greaves07} showed that such a simple empirical model can
reproduce the properties of both planets and debris discs around
Sun-like stars. The presence of debris indicates collisions between
planetesimals, such as icy outer-system comets, treated in the model
as a less successful mode of planet formation ending only in small
bodies. Both the frequencies and metallicity-dependencies of these
two outcomes were well matched. In particular, forming planets
requires a high solid mass to build the core, and so high
metal-fractions are advantageous, explaining the strong metallicity
trend in the host stars \citep{wyatt07}. On the other hand, forming
comet belts requires only modest amounts of material, and many discs
meet this threshold regardless of metallicity, so it enters only
weakly \citep{greaves07}. More complex approaches such as the
simulations of \citet{johansen09} have confirmed that the outcomes
of planetesimal formation should depend strongly on metallicity.
\citet{raymond07} have used a model of disc-masses in solids to
investigate the formation of terrestrial planets, finding in
particular that habitable examples could be less likely for low-mass
stars if they had insubstantial discs. A recent summary of our
understanding of planetesimal growth is presented by
\citet{chiang10}

\begin{figure}
\label{fig1}
\includegraphics[height=85mm,width=70mm,angle=270]{KMestA.prn}
\includegraphics[height=85mm,width=70mm,angle=270]{KMestFGK.prn}
\includegraphics[height=85mm,width=70mm,angle=270]{KMestM.prn}
\caption{Distribution of disc masses (in gas plus dust) for A, FGK and
M stars. Disc detections and upper limits are shown by filled and
unfilled bars respectively; the red lines shows the fitted population
from survival analysis (see text). X-axis labels are for the
bin-centres. For reference, a solar mass is just over 1000 Jupiter
masses; the lowest bin comprises discs of $<1.3$~M$_{Jupiter}$ or
$<$4~M$_{Earth}$ in solids at solar metallicity.
}
\end{figure}

Here our empirical model is expanded to test the planet and debris
frequency expected for stars of higher and lower mass, namely A and M
dwarfs. Their planetary systems are now being discovered by radial
velocity surveys (of subgiant descendents in the case of A-stars), and
their debris discs seen as far-infrared excesses e.g. by the Spitzer
satellite \citep{su06,gautier07}. Here the M,Z distributions are
assessed, and while the [Fe/H] ranges among the stars are all similar,
the disc-mass distributions differ markedly, leading to strong changes
in incidence with stellar mass.

\section{Data}

Metallicity data were first compared for different stellar types.
Subgiants that are the descendents of A-stars of $\approx$~1.2-2.5
solar masses can be identified by absolute magnitude M$_V$ of 0.5-3.5
and B-V of 0.55-1.0 \citep{johnson06}, or similarly effective
temperature 5000-6300~K \citep{flower96}. The SPOCS study of nearby
stars \citep{valenti05} yields approximately 100 such subgiants with
[Fe/H] measurements. A (log-)normal fit to these values indicates mean
and standard deviation $\mu,\sigma$ of +0.02, 0.23. This is very
similar to values for nearby Sun-like dwarfs, with $\mu,\sigma$ of
-0.01,0.29 in the data of \citet{valenti05}. For M-stars, the
derivation of metallicities based on photometry has recently be
re-computed by \citet{johnson09}. Using data they present for M-stars
within 20 pc, and their relation for [Fe/H] derived from K-magnitude
height above the main-sequence, $\mu,\sigma$ are found to be
-0.06,0.24. These distributions\footnote{\citet{schlaufman10} argue
that the mean [Fe/H] for FGK and M stars are about 0.1 dex lower than
found here; such small shifts have negligible effect on the
M$_{solids}$ distribution.} are so similar that they can not
significantly affect planet-building potential around different classes
of star.

\begin{table*}
 \caption{Predicted and observed percentages of giant planets and
debris discs. The planet data are from the uniform-detectability
sample of \citet{johnson10} and debris frequencies are from
Spitzer 70~$\umu$m surveys \citep{gautier07,trilling08,beichman06,su06},
limited to data with signal-to-noise$> 3$. The debris frequency observed
for FGK stars is used to make the predictions for other spectral types.
The upper limit for M-stars is from $<$1/15 objects detected by Spitzer
at 70~$\umu$m in $\geq 3\sigma$ data \citep{gautier07}, while the bracketed
value is from 1/50 detections in the millimetre survey of \citet{lestrade09}.
Poisson errors are based on the number of planet and debris systems
detected, and the numbers of discs falling within the input mass-bounds
(see text). The parameters of Equation~1 are listed in the four columns
after the spectral types. The A-star mass marked with a $^*$-symbol is
actually the second most-massive value from \citet{natta04}, as
another disc with large distance and uncertain mass was removed.
The sample of \citet{alonso09} gives a comparable highest mass of
225 Jupiter masses among HAe stars within about 500~pc.
}
 \label{tab:list}
 \begin{tabular}{@{}ccccccccc}
  \hline
spectral type & max $M_{disc}$& min. [Fe/H]   & min. $M_{solids}$ & min. $M_{disc}$ & \multicolumn{2}{c}{f(planet) \%} & \multicolumn{2}{c}{f(debris) \%} \\
              &	(M$_{Jup}$)   &	(planet)      &	(M$_{\oplus}$, for planet) & (M$_{Jup}$, at [Fe/H]=0)   & predict& observe & predict & observe \\
  \hline
FGK	& 600 	& -0.7	& 400	& 120 	& 8 $\pm$ 3	& 8.5 $\pm$ 1.0 & --	& 	19 $\pm$ 3 \\
A	& 150$^*$& -0.36& 200 	& 60	& 20 $\pm$ 7	& 20 $\pm$ 3	& 55 $\pm$ 11	& 70 $\pm$ 11 \\
M	& 270	& -0.12	& 650 	& 200	& 1.9 $\pm$ 1.3	& 3.3 $\pm$ 1.5	& 14 $\pm$ 4	& $< 7$
($\sim 2$)\\
  \hline
 \end{tabular}
\end{table*}

The disc mass distributions were then constructed and found to be quite
different for A versus FGK and M stars (Figure~1). The masses for the
A-stars are from \citet{natta04}, in a millimetre-based compilation for
young Herbig Ae objects\footnote{The data presented by \citet{alonso09}
give a similar mass-distribution but with fewer HAe stars.}. Their plot
of M$_{disc}$/M$_{star}$ versus M$_{star}$ was used to derive
M$_{disc}$ for M$_{star}$ of 1.2-2.5 solar masses. Since 11/45 of the
disc-mass measurements were upper limits, a survival analysis was used
to find the underlying distribution, implemented with the {\sc asurv}
package \citep{lavalley92}. The resulting disc-mass distribution
appears rather log-normal in form, and is centred at a mean of 25
Jupiter masses of gas and dust. For comparison, this is just above the
Minimum Mass Solar Nebula of $\approx$20 Jupiter masses \citep{davis05}
needed to supply the cores of all the Sun's planets. For FGK and M
stars\footnote{Pre-main sequence M-stars have rather vertical Hayashi
tracks with modest change in effective temperature, so young and old
objects of M-type will have similar stellar masses. Other pre-main
sequence stars evolve towards hotter spectral types, but since most of
the observed T Tauri objects have K-classifications, they will still be
within roughly `Sun-like' classes at later times.}, the
millimetre-surveys of Tau/Aur and Oph by \citet{andrews05,andrews07b}
were used, comprising similar-aged stars to the HAe's of
\citet{natta04}. Since these authors used a conversion of flux to dust
mass yielding values a factor of three lower than in \citet{natta04},
the masses have been shifted upwards by one bin ($\delta$log-M = 0.5)
in Figure~1. The difference arises from using an opacity at 1~mm
wavelength of 0.03~cm$^2$/g rather than the 0.01~cm$^2$/g found in the
detailed grain model of \citet{natta04}, which also agrees with recent
literature values \citep{draine06}. Survival analysis was again used to
construct the base disc-mass populations for FGK and M stars, but more
of the points are now upper limits, as the discs are less massive than
for the HAe stars . There were 57 limits out of 111 measurements for
the FGK dwarfs, and 57/106 for the M-stars. The means are 4 and 3
Jupiter masses respectively. There is also a large tail of
insubstantial discs in both cases, which combined with the lower means,
suggests that these discs will have much less planet-forming capability
than those of A-stars. The low-mass tails include many Class~III T
Tauri stars, which may lose their discs early on \citep{luhman10}. All
this disc-mass distributions are very broad, whereas the metallicity
distributions would cover only about one bin on the scales of Figure~1.

\section{Results and predictions}

\subsection{Planets}

Since the metallicity distributions are all narrow and similar, the
differences in M$_{disc}$ should be the dominant factor, and the
relevant thresholds for planet formation are now estimated.
\citet{greaves07} argued that if a minimum solid mass is needed to
grow a planetary core, then there will be observed examples of giant
planets arising from a range of discs with this same solid-mass, but
ranging between massive/low-metal and lower-mass/high-metal systems.
The lowest-metal system observed with a planet arises when the
solids criterion was met only because the disc was at the extreme
high-end of total mass. The solid mass can then be expressed using
the proxy data as
\begin{equation}
M^{min}_{solids} = M^{max}_{disc} 10^{\rm [Fe/H]_{min}}/100.
\end{equation}

The resulting parameters for all the spectral types are listed in
Table~1. The input data on the disc masses are from the references
in the previous section, and the metallicities of planet-hosts are
from http://exoplanet.eu/. The subsequent columns list the disc-mass
thresholds needed to form planets, and the predicted and observed
frequencies. For FGK stars, the predicted incidence of around 8~\%
is in excellent agreement with the 8.5~\% found in the systematic
planet sample of \citet{johnson10}. For A-stars the planet
agreement is also very good, with around 20~\% both predicted and
observed. For M-stars, the incidence is 2-3~\%, and the model
correctly predicts a low value because few of
the discs are massive (Figure~1).

The solid-mass thresholds needed for forming giant planets are here
found to be 200, 400 and 650~M$_{\oplus}$ for A, FGK and M spectral
types respectively (Table~1). This trend qualitatively matches the
expectation that grains will coagulate more readily in discs around
more massive stars, as the dynamical (orbital) times are shorter.

Subsets of different types of planetary systems are not considered
in detail here. Rocky `super-Earth' planets form a small part of the
FGK-hosted population (e.g. $\la10$~M$_{\oplus}$ bodies contribute
around 0.5~\% incidence), and the requirements for building such
bodies are considered briefly below. Also, there are likely to be
systems with more gas giants orbiting further out -- for example,
\citet{gould10} estimate $\sim1/6$ for the frequency of two-giant
systems, from microlensing results -- but these are not included
here. The analysis is however self-consistent in using the
lowest-metal radial-velocity detection to predict the planet
frequency in this type of survey.

\subsection{Debris discs}

For the less-successful outcome of building planetesimals, the
threshold can not be straightforwardly derived. \citet{greaves07}
used the metallicity range of FGK debris-hosts to estimate the
mass-threshold for building comet belts, but there is little similar
data for A and M-stars. Hence the simplest procedure is to infer the
threshold from the data. In deep Spitzer surveys, approximately
19~\% of FGK stars have debris discs detected at 70~$\umu$m
\citep{greaves10}, associated with belts of colliding planetesimals.
Associating this outcome with discs of lower-mass than those that
build planets, this next 19~\% of the M$_{disc}$ distribution
corresponds to $>$50 Earth masses in solids or $>$15 Jupiter masses
in total at solar metallicity. This is eight times lower than the
value needed for building planets. Making the minimal assumption
that the efficiency of the planetesimal-building process is similar
in all systems, a similarly reduced value can be applied for the
other stellar masses. Discs of at least 25 and 80 Earth masses in
solids would then be needed for A and M stars to have debris,
respectively.

The model results (Table~1) show a predicted debris incidence
for A-stars that is slightly low, but only at the -1$\sigma$ bound
(error on observed$-$predicted of 15~\% in quadrature). For M-stars,
there is however an interesting discrepancy between the debris
prediction of around 14~\% and the observed incidence of as little
as 2~\% (Table~1). A few other M-host debris discs have been
discovered, but all around stars much younger than typical
main-sequence lifetimes \citep{liu04,forbrich08, plavchan09}. It is
possible that these M-discs evolve down to lower dustiness than seen
among A or FGK stars (Figure~2), but given that FGK and M discs are
similar at the T Tauri phase, an alternative explanation is
favoured. As M-dwarfs are of low luminosity, dust particles in
thermal equilibrium with the star will be cool, compared to a grain
on the same orbit around AFGK stars. In this case, the discs will
have faint emission shifted to long wavelengths, and this introduces
a survey bias \citep[e.g.]{matthews07}. There may be negligible
excess above the M-star photospheres at 70~$\umu$m, and exploratory
millimetre surveys may not yet have reached deep enough to detect
the discs.

There are some systems known with both gas giants and debris discs,
but the statistics suggest the two phenomena are not correlated
\citep{bryden09}; here in fact it is assumed that they are separate
outcomes. In the uncorrelated case for FGK stars, a proportion of
0.08 (planet-frequency) $\times$ 0.19 (debris frequency) or 1.5~\%
of stars would have planets and dust, so this should be a negligible
population. A small number of planets around A-hosts of debris discs
has also been found in direct imaging surveys. These cases are not
considered as a separate class here. The stability of debris discs
over several Gyr, in terms of size and mass and largest bodies
present, has been studied in detail by \citet{heng10}.

\begin{figure}
\label{fig2}
\includegraphics[width=88mm,angle=0]{AFGKMdustlum.eps}
\caption{Distributions of percentage of stars versus fractional dust
luminosities for A and FGK systems with debris detected at
24,70~$\umu$m \citep{su06,trilling08, beichman06}. The typical upper
limit for the M-stars observed by \citet{gautier07} is shown at the
top, and applies for dust at temperatures of $\approx$~40-100~K.
}
\end{figure}

\subsection{Planet-building efficiency}

The planetary systems of M-dwarfs are of particular interest because
these stars are so abundant in the Galaxy. The large solid-threshold
inferred for M-stars to form planets (Table~1) suggests the core
accretion process may be inefficient. A rough efficiency estimate
can be made, assuming the solid content in the giant planets in the
M-star sample of \citet{johnson10} is of order 20~M$_{Earth}$. For
example, these planets are of $\ga$~M$_{Saturn}$ and this is about
Saturn's solid content \citep{saumon04}. A lower-mass `super-Earth'
population is also being discovered of which the most massive
observed in transit has a 20~M$_{Earth}$ core plus $\approx$10-15~\%
atmospheric mass \citep{bean08}. If 20 Earth-masses is a rough
boundary above which planets got on to become gas giants, then we
could equate this outcome at $\sim3.3$~\% frequency with the same
fraction of top-end discs, which have $\ga 300$~M$_{\oplus}$ of
dust. In this case, the efficiency of incorporating disc solids (or
at least the millimetre-detected component) is $\sim 20/300$ or
7~\%.

One consequence is that to build an Earth-mass planet,
15~M$_{\oplus}$ of dust in the disc would be needed. This occurs for
around 1/3 of M-stars (Figure~1), so they are not all expected to
form Earth-analogues. Given that about 3/4 of main-sequence stars
are of M-type, at least a quarter of nearby stars could host an
Earth-mass planet under this simple scaling approach. However, the
disc material spreads far from the star \citep{andrews07b}, so these
planets may only be able to form in large orbits, particularly as
the Hill radius for accretion increases with radial distance.
Microlensing is particularly well suited to discovery of M-dwarf
systems hosting such low-mass planets in orbits out to around 10~AU
\citep{gould10}.

It is beyond the scope of this study to assess how many planets
might migrate inwards, to the the habitable zone around 0.1~AU where
temperatures are suitable for liquid water on the surface.
\citet{raymond07} have simulated the formation of terrestrial
planets, and find that habitable planets of $\ga0.3$~M$_{\oplus}$
(the boundary for plate tectonics and atmospheric retention) are
much less likely for M-dwarfs than G/K stars. Figure~1 also suggests
that discs substantial enough to form an Earth-mass planet are much
more common around A stars, although these discs may also form
giants, with implications for the supply of material and stability
for terrestrial bodies \citep[e.g.]{fogg09}.

\section{Conclusions}

The disc-mass distributions for young A and M stars have been
constructed and compared to those for Sun-like stars. Inferring the
mass-thresholds for forming extrasolar planets and debris discs, the
higher incidence of both phenomena among A-stars or their descendents
is well reproduced, as is the low incidence of giant planets around
M-stars, which rarely have substantial discs. However, only about 2~\%
of nearby M-stars have debris discovered so far, and several times
more are expected. There may be a population of cool debris discs
around red dwarfs largely missed in surveys so far. The new SCUBA-2
camera \citep{holland06} observing at 450 and 850~$\umu$m will be
ideally matched to cold dust, as is the highly sensitive Herschel
satellite \citep{pilbratt10}, now observing out to 500~$\umu$m. A
further prediction of the empirical calculations is that about a third
of M-stars have the capability to form a 1 Earth-mass planet --
however, such planets are likely to form far from the star, as the disc
material is very spread-out.

\section*{Acknowledgments}

I thank Sean Raymond for very useful discussions, an anonymous
referee for helpful comments that improved the clarity of the paper,
and STFC for a fellowship in support of this work.