  The Galactic globular cluster NGC 1851 has raised much interest since HST
photometry revealed that it hosts a double subgiant branch. Here we report on
our homogeneous study into the cyanogen (CN) bandstrengths in the RGB
population (17 stars) and AGB population (21 stars) using AAOmega/2dF spectra
with R $\sim 3000$. We discover that NGC 1851 hosts a quadrimodal distribution
of CN bandstrengths in its RGB and AGB populations. This result supports the
merger formation scenario proposed for this cluster, such that the CN
quadrimodality could be explained by the superposition of two `normal' bimodal
populations. A small sample overlap with an abundance catalogue allowed us to
tentatively explore the relationship between our CN populations and a range of
elemental abundances. We found a striking correlation between CN and [O/Na]. We
also found that the four CN peaks may be paired -- the two CN-weaker
populations being associated with low Ba and the two CN-stronger populations
with high Ba. If true then s-process abundances would be a good diagnostic for
disentangling the two original clusters in the merger scenario. More
observations are needed to confirm the quadrimodality, and also the
relationship between the subpopulations. We also report CN results for NGC 288
as a comparison. Our relatively large samples of AGB stars show that both
clusters have a bias towards CN-weak AGB populations.
