abstract: Although Natural Language Processing (NLP) is at the core of many tools young people use in their everyday life, high school curricula (in Italy) do not include any computational linguistics education. This lack of exposure makes the use of such tools less responsible than it could be and makes choosing computational linguistics as a university degree unlikely. To raise awareness, curiosity, and longer-term interest in young people, we have developed an interactive workshop designed to illustrate the basic principles of NLP and computational linguistics to high school Italian students aged between 13 and 18 years. The workshop takes the form of a game in which participants play the role of machines needing to solve some of the most common problems a computer faces in understanding language: from voice recognition to Markov chains to syntactic parsing. Participants are guided through the workshop with the help of instructors, who present the activities and explain core concepts from computational linguistics. The workshop was presented at numerous outlets in Italy between 2019 and 2021, both face-to-face and online.
author: Ludovica Pannitto
University of Trento
`firstname.lastname@example.com`
Lucia Busso
Aston University
`firstname.lastname@example.com`
Claudia Roberta Combei
University of Bologna
`firstname.lastname@example.com   ` Lucio Messina
Independent Researcher
`firstname.lastname@example.com`
Alessio Miaschi
University of Pisa
`email@example.com   `
Gabriele Sarti
   University of Trieste
`   email@example.com`
Malvina Nissim
University of Groningen
`email@example.com`
bibliography: anthology.bib; custom.bib
title: Teaching NLP with Bracelets and Restaurant Menus:
  An Interactive Workshop for Italian Students

# Introduction

*Have you used Google this week?* This question would kick off the activity that we describe in this paper every time we delivered it. And a number of follow-up comments would generally appear. *What for? Translating, getting some help for homework, looking for info, writing collaboratively - and getting spelling correction!*

In our workshops, we talk to groups of teenagers – even if someone has not personally used any of those tools on a daily basis, it is utmost unlikely that they have never interacted with a vocal assistant, wondered how their email spam filter works, used text predictions, or spoken to a chat-bot. Also, applications that do not require a proactive role of the user are growing: most of us, for example, are subject to targeted advertising, profiled on the content we produce and share on social media.

Natural Language Processing (NLP) has grown at an incredibly fast pace, and it is at the core of many of the tools we use every day.[^1] At the same time, though, awareness of its underlying mechanisms and the scientific discussion that has led to such innovations, and even NLP's very existence as a scientific discipline is generally much less widespread and is basically unknown to the general public .

A concurrent cause to this lack of awareness resides in the fact that in more traditional high-school formal education systems, such as the Italian one, "young disciplines" such as Linguistics and Computer Science tend to be overlooked. Grammar, that in a high-school setting is the closest field to Linguistics, is rarely taught as a descriptive discipline; oftentimes, it is presented as a set of norms that one should follow in order to speak and write correctly in a given language. While this approach has its benefits, it is particularly misleading when it comes to what actual linguistic research is about. Similarly, Computer Science is often misread by the general public as an activity that deals with computers, while aspects concerning information technology and language processing are often neglected. This often leads to two important consequences. First, despite the overwhelming amount of NLP applications, students and citizens at large lack the basic notions that would allow them to fully understand technology and interact with it in a responsible and critical way. Second, high-school students might not be aware of Computational Linguistics as an option for their university degree. Oftentimes, students that enrol in Humanities degrees are mainly interested in literature and they only get acquainted with linguistics as discipline at university. At the same time, most Computer Science curricula in Italian higher education rarely focus on natural language-based applications. As a result, Computational Linguistics as such is practically never taught before graduate studies.

As members of the Italian Association for Computational Linguistics (AILC, www.ai-lc.it) we have long felt the need to bridge this knowledge gap, and made dissemination a core goal of the Association. As a step in this direction, we have developed a dissemination activity that covers the basic aspects of what it means to process and analyze language computationally. This is the first activity of its kind developed and promoted by AILC, and to the best of our knowledge, among the first in Italy at large.

This contribution describes the activity itself, the way it was implemented as a workshop for high school students in the context of several dissemination events, and how it can serve as a blueprint to develop similar activities for yet new languages.

# Genesis and Goals

We set to develop an activity whose main aim would be to provide a broad overview of language modeling, and, most importantly, to highlight the open challenges in language understanding and generation.

Without any ambition to present and explain the actual NLP techniques to students, we rather focused on showing how language, which is usually conceptualized by the layperson as a simple and monolithic object, is instead a complex stratification of interconnected layers that need to be disentangled in order to provide a suitable formalization.

In developing our activity, we took inspiration from the *word salad* Linguistic Puzzle, as published in :

> Charlie and Jane had been passing notes in class, when suddenly their teacher Mr. Johnson saw what was going on. He rushed to the back of the class, took the note Charlie had just passed Jane, and ripped it up, dropping the pieces on the floor. Jane noticed that he had managed to rip each word of the message onto a separate piece of paper. The pieces of paper were, in alphabetical order, as follows: *dog, in, is, my, school, the*. Most likely, what did Charlie's note originally say?

The problem includes a number of follow up questions that encourage the student to reflect upon the boundaries of sentence structure. In particular, we found that the *word salad* puzzle would give us the opportunity to introduce some of the core aspects of Computational Linguistics' research. Approaching the problem with no previous knowledge helps raising some crucial questions, such as: *what are the building blocks of our linguistic ability that allow us to perform such a task?*, *how much knowledge can we extract from text alone?*, *what does linguistic knowledge look like?*

Since the workshop here presented is the first activity of this kind in the Italian context, we took inspiration from games and problems such as those outlined in and used for the North American Computational Linguistics Olympiads, similar to the ones described in and . Particularly, we were inspired by the institution of (Computational) Linguistic Olympiads in making our workshop a problem-solving game with different activities, each related to a different aspect of computational language processing. Linguistic Olympiads are now an established annual event in many parts of the world since they first took place in Moscow in 1965. In these competitions students (generally of high-school age) are faced with linguistic problems of varying nature, that require participants to use problem-solving abilities to uncover underlying patterns or rules in the data. For an in-depth discussion of the history and diffusion of Linguistic Olympiads in the world, see and .

```latex
\begin{table*}[ht!]\center{
\begin{tabular}{lcp{8cm}}
\hline
\multicolumn{1}{c}{\textbf{Activity}}               &
\StopWatchStart &
\multicolumn{1}{c}{\textbf{Aim}} \\
\hline
1. Get to know a (computational) linguist  & 10'   & collaboratively build a definition for \textit{linguistics} as a study field  \\ \hline
2. Are computers able to \textit{hear}?    & 15'   & familiarize with the concept of \textit{simulation} of humans' speech perception abilities \\ \hline
3. Are computers able to \textit{read}?    & 30'   & introduce corpora as sources of linguistic knowledge and statistical patterns as structural aspects of language \\ \hline
4. Do computers \textit{know} grammar?     & 30'   & introduce human annotation and meta-linguistic knowledge as powerful research tools \\ \hline
5. Becoming a computational linguist       & 5'    & evaluate pros and cons of the two presented approaches, future directions and discuss about what's needed to become a computational linguist \\
\hline
\end{tabular}
}
\caption{Sections of the activity, with their planned duration and a broad aim for each of them.}
\label{tab:table-times}
\end{table*}
```

In the choice of algorithms to include in our dissemination activity, we decided to leave aside neural networks and instead focus on traditional statistical approaches, both for historical reasons and for the fact that these convey more clearly the distinction between different layers of linguistic information and their roles in language modeling. A brief discussion of most recent NLP technologies, including the application of neural networks, is included in the final part of the workshop (Sec.  ).

The activity is targeted at students in their last year of middle school (13 years of age) or older. While we believe 13 is a good entry point, there isn't an actual upper-bound, since the activity can be enjoyed by people of any older age (though in practice participation was mainly offered to schools, with the oldest students being 18-19). We thought this would be the appropriate target audience of this workshop for two main reasons. On the one hand, we believe that coming to the activity with a richer metalinguistic background, typically acquired during the first Italian school cycle, would be beneficial for the attendees to better grasp the differences between the scientific approach to language and the more prescriptive approach they are exposed to in school. On the other hand we also conceived the activity as a way of helping students in their future study choices: we therefore included both students attending their last year of middle school and therefore about to choose a high school curriculum as well as high school students, the latter in order to provide them with more options for their university choices.

# The activity

We planned our dissemination activity for a 90 minutes time slot, divided into five parts, as detailed in Table .

## Get to know a (computational) linguist

The first 15 minutes of the workshop are organized both as an ice-breaker activity for the attendees, and as a brief introduction to linguistics and computational linguistics more specifically.

During the introduction we tried to demystify some common misconceptions about linguistics, (i.e., *a linguist knows many languages*, *linguists get sometimes confused with speech therapists*, *a linguist will correct my grammar*, etc.): we presented participants with a list of possible definitions, and they had to identify appropriate ones. We broadly defined linguistics as the study of language as a biological, psychological, and cultural object. Computational Linguistics was then introduced both as a commercial and engineering-oriented field, as well as a purely scientific research discipline.[^2]

We also briefly discussed *linguistic questions* with participants as an example of the kind of problems that a linguist tries to approach during their research activity. These included: "How many ways of pronouncing *n* do we have in Italian?", "Do numbers from one to ten exist in all languages?".

For the following parts of the activity, the introduction to each sub-part was dedicated to a reflection upon what it means for us humans to *hear*, *read* and *understand grammar*, and whether there is a difference when the same tasks are performed by computers.

## Are computers able to *hear*?

As vocal assistants such as Alexa, Siri, Google Home etc. have become increasingly popular, we chose them as tangible examples of NLP technologies to kick-off the games. The aim of this section of the workshop was to demonstrate two points:

-   computers do not necessarily solve linguistic tasks the way we solve them; they are therefore *simulating* our abilities without *replicating* them;

-   consequently, the concepts of *easy* and *difficult* tasks have to be carefully revised when applied to language models.

We introduced the McGurk effect , to clarify how hearing language is a complex task in itself, involving a broad set of aspects beyond simple sound perception, such as the visual system as well as the expectations regarding the upcoming input. Computers on the other hand *hear* on the basis of an audio signal that is processed (Figure ), at least in the most traditional and well-established architectures, without access to higher level linguistic knowledge or information from the communicative context.

We then briefly presented speech recognition as a direct optimization task (i.e., given an audio signal, find the word in a given database that maximises the probability of being associated to that signal) and introduced one of the major challenges that speech recognition models are still facing, namely the ability to adapt to different speech styles (i.e., speakers of different language varieties and dialects, non-native speakers, speakers with impairments etc.).

In order to further demonstrate this, we tested the attendees' ability to adapt their hearing skills to different speech varieties by making them hear conversations in various Italian regional accents.[^3] While we asked attendees to guess the name of the region of the speakers, the actual aim was to show how we easily adapt to understand speech, differently from speech recognition systems.

## Are computers able to *read*?

From this moment on, the attendees worked on written text. The activity described in this section was aimed at showing how salient aspects related to language structure can be derived from the statistical properties of language.

Following the "Word Salad" puzzle , the ability of *reading* was presented as follows: given a set of words, are we able to rearrange them in a plausible sentence-like order? We demonstrated how this is an easy task for human beings, when one deals with a language they are familiar with (Figure ), while, generally one may not be able to perform the same task in case of unknown languages (Figure ), where each possible ordering seems equally plausible.

We therefore gave the attendees a deck of unknown, mysterious tokens (left card in Figure ) and asked them to come up with the most plausible sentence that contained all of them. The cards represented either Italian or English words (participants were divided into two teams, each one dealing with a different *masked* language) which had however been transliterated into an unknown alphabet. While this was obviously an impossible task to solve, it gave us the opportunity to introduce a notion of probability in the linguistic realm. We cast it as the expectation that we humans bear on the appearance of a specific linguistic sequence, and the subsequent need for a source of linguistic knowledge to implement the same notion.

When asked to perform the same task on the words of Figure , participants produced sentences in a quite consistent order, and the most prototypical sentences (e.g., *il mio cane è nel giardino*, en. *my dog is in the garden*) were usually elicited before some less typical ones (e.g., *il cane nel giardino è mio*, en. *the dog in the garden is mine*), while agrammatical sentences (i.e., random permutations) were never produced.

We justified their responses by explaining that humans accumulate a great amount of linguistic knowledge throughout their lifetime that helps them refine these expectations, while machines are instead in principle unbiased towards having any specific preference. This observation allowed us to introduce participants to the notion of *corpus* as a large collection of linguistic data that mimics the amount of data we are exposed to as humans.

Each team was then provided with a corpus written in an unknown language (approximately 60 sentences hand-crafted by transliterating a portion of the "Snow White" tale into a mysterious alphabet Figure ). Concurrently, we introduced a simple algorithm to tell apart sentence-like orderings of the provided tokens from the random ones.

The algorithm, which we called *The bracelet method* (Figure ), is based on the Markov Chain procedure: in a scenario similar to that of lining pearls up to form a bracelet, participants could decompose the task of forming up a sentence into smaller tasks. To make the operation more concrete, we equipped each card with a button loop as shown in Figure : this way cards could be physically threaded together to form a sentence.

The first step consisted in choosing the first word, for the beginning of the sentence. Since participants were facing a language they didn't know, they were not aware of language structures nor of the meaning of the tokens. In such a situation, they could decide whether it was plausible to use a given word at the beginning of a sentence just by looking up in the corpus sentences that began the same way. If they found at least one sentence that began with the same word, it meant that that was a licit position and they could use it to start the sentence.

The activity continued as follows: sticking to the bracelet metaphor, participants needed to select and insert the following "pearl" into the thread: ideally the pearl should pair well with the previous one, as we might want to avoid colour mismatch (e.g., it is well known that blue and green do not go well together). The metaphor highlights therefore a core aspect that holds true for language as well, namely that we can condition our choice on a variable number of previous choices, and this will influence the complexity of the obtained pattern.

The activity resulted in a number of sentence-bracelets, as shown in Figure , that were then kept aside to be translated at the end of the workshop.

## Do computers *know* grammar?

While the previous game highlighted the importance of statistical information in NLP, in accordance with the overall aim of the activity, we also wanted to introduce some of the algorithms that are more deeply rooted in the linguistic tradition.

In order to do so, we introduced the notion of grammar as a descriptive abstraction over a set of examples. Out of the linguistic context, to exemplify this notion of grammar metaphorically, we presented participants with a set of possible restaurant menus (Figure ), and encouraged them to come up with the general rule that the restaurant owner must have had in mind when choosing those combinations. All menus were built as a traditional Italian full meal, composed by two main dishes and a dessert. We perpetuated the metaphor showing how, once a set of rules is defined, these can be used both to decide if a new menu is likely to be part of the same set (Figure ) and also to generate new meals (Figure ).

This metaphor, which was readily grasped by most participants, was useful to show how different components can be combined together in a meaningful way, as it happens in grammar. Before moving back to the corpus, we showed them what a formal grammar of the menus could look like.

We had previously annotated the corpus with syntactic and morpho-syntactic information (i.e., part of speech), as shown in Figures  and . We therefore asked participants to extract from the corpus a possible grammar, namely a set of attested rules and use them to generate a new sentence.

In order to write the grammar, participants were given several physical materials: felt strips reproducing the colors of the annotation, a deck of cards with numbers (identifying parts of speech) and a deck of "$=$" symbols (Figure ).

With a new deck of words (Figure , right panel), not all of which present in the corpus, participants had to generate a sentence using the previously composed rules.

## Becoming a computational linguist

At this point, participants had created a number of sentences by means of the two techniques described above. It is now time to discover that the mysterious languages they worked on were actually English and Italian. This was achieved in practice by superimposing a Plexiglas frame on the A3 corpus pages (Figure ): the true nature of the corpora was this way revealed as the participants could see the original texts (in Italian and English) and translate the sentences they had created previously.

The outcome of the activity stimulated discussion amongst the participants, highlighting pros and cons of each approach and how could they be integrated into real-life technologies. Our workshop ended with a brief description of more recent NLP technologies and their commercial applications, such as recommender systems, automatic translation, text completion, etc.

Since the target audience consisted mostly of middle- and high-school students, we offered an overview of what it takes to become a computational linguist and where one could study computational linguistics in Italy.

# The workshop in action

The activity – here outlined in its complete and original form – was presented in various outlets during the last year and a half.

It was initially designed for the 2019 edition of the "Bergamo Scienza" Science Festival[^4], where it was presented live to over 450 participants in the course of two weeks. A simplified "print-and-play" version of the workshop was also presented at the 2020 edition of the SISSA[^5] Student Day.

Due to the Covid-19 pandemic, all other presentations of the activity had to be moved online. Transposing the workshop crucially meant striving to maintaining the interactive nature of the activities without the possibility of meeting face to face. To do so, we integrated our original presentation on Google slides with the interactive presentation software Mentimeter[^6] - which we used for questions, polling and quizzes. The corpus and tokens were presented via a web interface created especially for this purpose[^7].

The interface presented the masked corpus, complete with POS tags and syntactic annotations. For the bracelet activity participants were automatically assigned a number of tokens which they could use to build a sentence by simply dragging and dropping them. (Figures  and  ).

This online version was crafted in the first place to be presented at "Festival della Scienza"[^8] (Figure ), a science festival primarily aimed at school students held each year in Genoa, where multiple 45' sessions of the workshop were run over the course of two days. The fourth activity (Section ) involved "bootstrapping" syntactic rules based only on our color-based annotation. To simplify online interaction, we only used the unmasked Italian version of the corpus and participants played collectively, helping each other to build sentences and grammatical rules: rules were collected through a Mentimeter poll, while a sentence was generated in a guided demonstration by the presenter of the workshop.

Overall, the workshop transposed really well online, and was extremely well-received in this version as well. The online modality also allowed us to present it to a more vast and varied audience than just students: a version for the general public was presented at *European Researchers' Night* (*Bright Night*[^9]) at the University of Pisa, thanks to a collaboration with the Computational Linguistics Laboratory[^10] and ILC-CNR[^11] in November 2020, a dedicated session was run for the High school ITS Buzzi[^12] (Prato, Italy) in December 2020, and during the second edition of the Science Web Festival[^13] in April 2021.

# Reusability: this activity as a blueprint

As mentioned in Section , our activity was inspired by a collection of English-language problems created for students participating in the Computational Linguistics Olympiad .

We adapted the original activity to the Italian language and context. While we tried to choose widely shared linguistic principles to communicate, the operation of adapting the game to a different language obviously requires language specific choices and details, which would have to be re-evaluated when porting the activity to yet new languages. Particularly, since the materials have been developed for Italian, it might be the case that transposition is not straightforward for some languages. However, we believe that the general structure of our workshop can serve as a blueprint for extension to new languages. For this reason all the relevant materials are made available in a dedicated repository. The repository includes both scripts to reproduce our activity as well as a general set of insights/recommendations regarding the structure and principles of the workshop.

All of the scripts necessary to produce the materials used in the game's workflow in a different language are made available in our open-access repository[^14]. To get the activity into production using the scripts provided, one only needs to create an annotated corpus in the target language.

Our specific choice of "Snow White" as a corpus is motivated by the fact that the story can be phrased in a fairly repetitive formulation, with two advantages. One is that enough bigrams are present that enable the generation of new sentences with the *bracelet method*. The other one is that the story contains recognizable characters (e.g., the dwarfs, the evil queen etc.), so that, when the unmasked text is revealed, the process results intuitively transparent. Such characteristics are desirable for the activity, and should be kept in mind when choosing a new text for a new language.

In addition to sharing scripts, we are sharing here the core structure and principles the workshop relies on, which can be reproduced when replicating the activity also for a different language.

#### Parts 1&2

At the beginning of the workshop (see Section ) we show the limits of current technologies, in particular in terms of adaptability. To this end, we exploited diatopic variations, such as Italian regional accents, since this is an aspect readily available to Italian students. In the context of other languages, the same concept could be however shown by different means, such as the influence of jargon or minority languages, or simply differences in accuracy depending on, age or other socio-demographic and socio-cultural variables.

#### Part 3

The next part of the activity (Section ) is the most easily reproducible in a different language as it only exploits statistical co-occurrences as a cue for linguistic structure. The only pre-requisite here is the availability of a sufficiently standardized tokenization for the language of interest. As described above, we masked the language through a transliteration system: this was achieved either substituting words with sequences of symbols, or with non-words. In the case of non-words, these were generated by manually defining a series of phonotactic constraints for Italian, which should be adapted to the target language.

The main message to be conveyed through this activity is that language is a complex system; we did this disentangling semantics from the purely symbolic tier, which is the one computers most commonly manipulate.

#### Part 4

The following part of the activity (Section ) focuses on the expressive power residing in the definition of auxiliary categories as descriptors of linguistic evidence. We achieved this by simplifying a constituent-based annotation for our corpus: having continuous constituents easily allowed for the physical implementation of re-writing rules (i.e., participants had some physical constituents and tokens that could be used to simulate the generation process). We built categories in order to extract from the text a simple *regular grammar*, and we were especially careful about the fact that both the Italian and the English corpus showed a similar structure and therefore complexity level. More specifically, we restricted to five types of higher-level categories that acted as phrases (i.e., sentence, noun phrase, verb phrase, prepositional phrase, subordinate clause): this is of course a huge simplification and, for the sake of the activity, we overlooked some relevant linguistic phenomena.

We reckon that this approach might not be portable to languages that exhibit a flexible word order, so alternative solutions should be sought.

#### Part 5

The last section of the activity was dedicated to a discussion on the presented methods for language generation. After that and depending on the audience, we presented some options to pursue studies in Computational Linguistics in Italy, which would of course have to be adapted to the target social context. Lay publications concerned with Computational Linguistics are also unfortunately not very common in Italy, therefore we took the opportunity to provide participants with some suggestions for further readings .

# Looking back and ahead

In the previous sections we described an interactive workshop designed to illustrate the basic principles of Computational Linguistics to young Italian students. It is the first activity of its kind to be designed by the Italian Association for Computational Linguistics and among the first dissemination activities in Italy for Computational Linguistics directed to young students.

The activity had the broad aim of increasing awareness towards applications based on language technology, and introduce students to the study of language as a scientific discipline.

We run the activity in both face-to-face and, due to COVID-19, online form: generally speaking, we received enthusiastic feedback both from younger participants and from the more general public. We adapted the activity to a variety of formats and time-slots, ranging from 30 to 90 minutes: the amount of time required to approach the game and get acquainted with the concepts is of course variable, and depends on the participants' background and on the level of engagement that is expected of them. Generally speaking 45 minutes are enough for a presentation including some interaction with the audience, especially in the online setting, but at least 90 minutes are needed for the participants to fully experiment with the hands-on activity.

We want to specifically stress how time is a central ingredient in the activity. While the game-related aspects remained engaging and fun even in the shortened, online versions, in order to fully grasp the mechanisms underlying the presented algorithms longer sessions would be needed. In fact, we often felt that more time would be beneficial for a deeper discussion about language as an object of study itself, and about language as data on which theories can be built. In particular, shorter time-slots or less guided activities enhance the risk of participants approaching the challenge as a puzzle that they can solve regardless of linguistic knowledge. This is because the approach to language we are presenting is entirely new to our audience, and not only to the younger students.

Although we did not have a formal system in place to collect systematic feedback, the overall response across venues and conditions has been extremely positive.[^15] Curiosity and engagement of participants remained high both onsite and online, and we received many questions on several aspects of Natural Language Processing and neural networks, as well as concerning its role in Artificial Intelligence at large.

The participants' enthusiastic questions gave us many ideas for future dissemination activities. In fact, the technological world is advancing fast and we firmly believe that it is necessary to spread more awareness on the inner workings of AI-based technologies, to develop a society-level conscience to approach them in a critical way.

The activity described in this paper was targeted at middle to high-school students as well as the general public. It would be interesting to engage with a younger audience as well, as communicating the study and (computational) modelling of language to them would raise awareness towards language studies as a scientific discipline.

For our activity, we took inspiration from one of the problems proposed in . Puzzles such as the ones presented at the Computational Linguistics Competitions are a great way to introduce both important challenges and the methodology to solve them, as they stimulate students to investigate linguistic aspects in a bottom-up fashion. Organizing the competition in Italy would represent a bigger project for our association, to be addressed in the coming years.

# Acknowledgements

The workshop has been developed with help and support from many parties. We would like to thank them all here. First and foremost, all the participants that enthusiastically played along and made the activity a success. Along with them, four tutors and science animators helped us greatly in putting our ideas into practice. We are also grateful to BergamoScienza, Festival della Scienza, Science Web Festival for hosting the activity during the festivals, to ILC-CNR "Antonio Zampolli" and ColingLab (University of Pisa) for hosting us during the European Researchers' Night, and to the Scuola Internazionale Superiore di Studi Avanzati (SISSA) and ITI Tullio Buzzi. We are also grateful to Dr. Mirko Lai, who has collaborated on the development of the web interface for the online versions of our activity. We would like to thank AILC (Italian Association of Computational Linguistics) for encouraging us to design the activity and supporting us throughout the process. Last but not least, we thank the true hero of our workshop: GingerCat[^16].

[^1]: In this discussion, and throughout the paper, we conflate the terms Natural Language Processing and Computational Linguistics and use them interchangeably.

[^2]: We relied on the definition reported in <https://www.thebritishacademy.ac.uk/blog/what-computational-linguistics/>.

[^3]: Conversations were extracted from corpus CLIPS .

[^4]: <https://festival.bergamoscienza.it/>

[^5]: <https://www.sissa.it/>

[^6]: <https://www.menti.com/>

[^7]: A demo of the interface can be found at <https://donde.altervista.org/>

[^8]: <http://www.festivalscienza.it>

[^9]: <https://www.bright-night.it/>

[^10]: <http://colinglab.humnet.unipi.it/>

[^11]: <http://www.ilc.cnr.it/>

[^12]: <https://www.tulliobuzzi.edu.it/>

[^13]: <https://www.sciencewebfestival.it/>

[^14]: <https://bitbucket.org/melfnt/malvisindi>

[^15]: Following one reviewer's suggestion, we have implemented such a system for our latest presentation at the Science Web Festival 2021 as a Google form that we circulated among participants at the end of the presentation.

[^16]: <https://icons8.com/illustrations/style--ginger-cat-1>


abstract: We present an analysis of compact star clusters in deep HST/WFPC2 images of NGC 1275. B and R band photometry of roughly 3000 clusters shows a bimodality in the B-R colors, suggesting that distinct old and young cluster populations are present. The small spread in the colors of the blue clusters is consistent with the hypothesis that they are a single age population, with an inferred age of 0.1 to 1 Gyr. The luminosity function shows increasing numbers of blue clusters to the limit of our photometry, which reaches several magnitudes past the turnover predicted if the cluster population were identical to current Galactic globulars seen at a younger age. The blue clusters have a spatial distribution which is more centrally peaked than that of the red clusters. The individual clusters are slightly resolved, with core radii $\, \buildrel < \over \sim \,0.75$ pc if they have modified Hubble profiles. We estimate the specific frequencies of the old and young populations and discuss the uncertainties in these estimates. We find that the specific frequency of the young population in NGC 1275 is currently larger than that of the old population and will remain so as the young population evolves, even if the majority of the low mass clusters are eventually destroyed. If the young population formed during a previous merger, this suggests that mergers can increase the specific frequency of globulars in a galaxy. However, the presently observed young population likely contains too few clusters to have a significant impact on the overall specific frequency as it will be observed in the future.
author: Matthew N. Carlson, Jon A. Holtzman, Alan M. Watson, Carl J. Grillmair, Jeremy R. Mould, Gilda E. Ballester, Christopher J. Burrows,
John T. Clarke, David Crisp, Robin W. Evans, John S. Gallagher III,
Richard E. Griffiths, J. Jeff Hester,
John G. Hoessel, Paul A. Scowen, Karl R. Stapelfeldt, John T. Trauger, and James A. Westphal
title: Deep HST Observations of Star Clusters in NGC 1275 [^1]

# Introduction

Early HST observations of NGC 1275, the central galaxy in the Perseus cluster, revealed a population of about 60 blue (V-R $\sim$ 0.3) star clusters surrounding the nucleus (Holtzman *et al.*  1992). The color of these clusters suggests an age of roughly 300 million years based on the models of Charlot and Bruzual (1991). Spectra of the brightest object (Zepf *et al.*  1995) also suggest an age of 0.1 to 0.9 Gyr, based on a comparison of the observed line widths to those predicted by the models of Bruzual and Charlot (1993). None of the clusters seem to have H$\alpha$ emission, with the exception of one object found by Shields & Filippenko (1990), which appears to be a much younger object. The blue clusters in the original WFPC1 images appear unresolved, suggesting sizes of less than 15 parsecs. The brightest object has V = 18.9, which corresponds to M$_V$ = -15.8 for $H_0$ = 75 km/s/Mpc and cz=5264 km/s (Strauss *et al.*  1992). The brightnesses of the blue clusters suggest masses of between 2 $\times$ 10$^4$ and 1 $\times$ 10$^8$ M$_{\odot}$, depending on the assumed age and distance, if one assumes that the objects are star clusters with a Salpeter IMF. The observed sizes, luminosities, and the estimated masses suggest that these objects may be young analogues of globular clusters.

In the past several years, massive young clusters have been observed in a variety of other galaxies. Lutz (1991) detected young globular cluster candidates in a ground-based study of the merger remnant NGC 3597, and these have been confirmed to be compact by HST observations (Holtzman *et al.*  1996). Candidate young globular clusters have also been found in other interacting systems, including NGC 7252 (Whitmore *et al.*  1993), NGC 4038/9 (Whitmore & Schweizer 1995), NGC 3921 (Schweizer *et al.*  1996), among others. A few massive clusters are seen in the starburst galaxy NGC 253 (Watson *et al.*  1996), while others have been detected in the ring galaxies NGC 1097 and NGC 6951 (Barth *et al.*  1995) and the dwarf galaxies NGC 1569, NGC 1705 (O'Connell *et al.*  1994), and He 2-10 (Conti & Vacca 1994). While no generally accepted picture has emerged as to what conditions lead to the formation of young clusters, galaxy interactions appear to be an important component.

It is particularly difficult to determine the precise mechanism responsible for the presence of these objects in NGC 1275 because of the myriad peculiarities of this galaxy, including the presence of a significant amount of dust, streamers of H$\alpha$ emission, an active nucleus, and the location of the galaxy at the center of a cooling flow in the Perseus cluster. Two hypotheses have been offered for the origin of these clusters, namely, that they were formed from the substantial mass deposition of the Perseus cluster cooling flow (200 M$_{\odot}$/yr) or that the cluster formation was triggered by a galaxy-galaxy interaction. Holtzman *et al.*  (1992) preferred the latter hypothesis based on the observed lack of spread in the WFPC1 colors, which implies a common age for the objects, and the appearance of a ripple in the galaxy light which suggests that a previous interaction may have occurred. Richer *et al.*  (1993) found larger color spreads in high resolution CFHT images and preferred the cooling flow hypothesis. However, Holtzman *et al.*  (1996) did not detect young clusters in three of a sample of four other cooling flow galaxies; the central cluster galaxy in Abell 1795 may have a few clusters, but it also has a peculiar morphology which suggests a previous interaction.

The merger hypothesis is particularly interesting in light of the theory that elliptical galaxies form through mergers and the observation that ellipticals have higher specific frequencies of globular clusters than spirals, where the specific frequency is a measure of the number of globular clusters per unit luminosity of the host galaxy. Ashman & Zepf (1992) investigate the proposition that clusters might form during mergers and suggest that, in such an event, cluster formation would occur in a brief burst, resulting in a set of newly formed clusters with a common age and color. This would lead to a cluster system with a bimodal distribution of cluster colors reflecting the difference in metallicity or age (or both) between the original and the newly formed clusters. They predict that the spatial distribution of the younger clusters would be more sharply peaked toward the center of the galaxy than that of the old cluster system, because the old globular cluster populations of the two progenitor galaxies would remain spatially extended and probably be dynamically heated during the merger, while the new clusters would be formed out of gas which becomes more centrally concentrated during a merger. However, as Harris *et al.*  (1995) note, an increase in specific frequency as a result of a merger requires not only that globular clusters form during such an event, but that they form *preferentially* over non-cluster stars as compared to the ratio of clusters to background stars in the progenitor galaxies.

Few estimates of the specific frequency of these young cluster systems exist. Watson *et al.*  (1996) suggest that the 4 young clusters in NGC 253 most likely have formed with a large specific frequency. In the merger remnant NGC 3921 Schweizer *et al.*  (1996) find that the blue clusters will increase the overall number of clusters enough that the galaxy will come to have the specific frequency of an elliptical within 7 Gyr. Miller *et al.*  (1997) find that the specific frequency of globular clusters in the merger remnant NGC 7252 will rise over the next 15 Gyr as the background population fades to resemble that of an elliptical. These results all suggest that interactions can increase the specific frequency in a galaxy.

Luminosity functions of old globular cluster systems have been well studied and used as a part of the extragalactic distance scale because of their uniformity from galaxy to galaxy; they are roughly Gaussian in shape with a peak near M$_V$ = -7.3 (Harris 1996). The luminosity functions of most of the recently discovered young cluster systems are poorly determined because of either small number statistics (few young clusters in the galaxy) or incompleteness. A notable exception is the young cluster system in NGC 4038/9, which shows no turnover to two magnitudes fainter than the turnover predicted for a typical old cluster system, even after allowing for the expected fading of the clusters based on stellar population models. Based on this, van den Bergh (1995) argues that this cluster system may be intrinsically different from old globular cluster systems. However, Mateo (1993) notes that, at least in the LMC, the observed increasing cluster luminosity function is well modelled by a combination of globular and open clusters. Also, several authors have recently suggested that substantial numbers of clusters may be destroyed over a Hubble time (Gnedin & Ostriker 1997; Elmegreen & Efremov 1997); mechanisms include evaporation and tidal disruption by galactic bulges and disks.

We have obtained WFPC2 images of NGC 1275 which go about 4 magnitudes deeper (in the red) than the WFPC1 observations of Holtzman et al. (1992). These observations provide more accurate colors than previous observations and probe the cluster population to approximately 2.5 magnitudes fainter than the turnover expected if this population is identical to the Galactic globular cluster system seen at an younger age.

Section 2 briefly summarizes the observations and the reductions. Section 3 discusses our analysis and potential sources of error. In section 4 we present our results including the photometry, luminosity function of the clusters, surface density distribution, an estimate of the specific frequency, and analysis of the sizes of the objects.

# Observations & Data Reduction

Observations of NGC 1275 were made with the WFPC2 on 1995 November 16 in the F450W and F702W filters. Exposure times in F450W were 200, 1000, and 3 $\times$ 1300 seconds. In F702W the exposures were 200, 900, 1000, and 2 $\times$ 1300 seconds. The F450W and F702W filters are roughly similar to broad band B and R filters.

The data were reduced using the procedures discussed by Holtzman *et al.*  (1995a). The individual exposures were all taken with a common pointing, and the frames in each color were combined. Cosmic rays were rejected in the averaging based on the expected variance from photon statistics and readout noise; an extra term was included in the expected variance which was proportional to the signal in the pixel to allow for small pointing differences between frames.

# Analysis

Combined color images of NGC 1275 in the PC and in the WFPC2 fields are shown in Figures 1a and 1b, respectively. The nucleus is in the center of the PC, and the other bright point source in the PC is a foreground F star (Hughes & Robson 1991). Almost all of the remaining point sources are candidate star clusters in NGC 1275.

To detect the compact clusters, we summed the two colors and subtracted a 5 $\times$ 5 boxcar-smoothed version of this image. We divided this by the square root of the smoothed image to provide an image whose units were proportional to the local noise (assuming the noise is dominated by the background, as is the case here). We then used IRAF's DAOFIND to detect objects on the resulting image, in order to provide a uniform detection threshold in signal-to-noise ratio across the image. The choice of detection threshold was determined by minimizing the number of spurious detections and maximizing the number of faint objects detected, judging from a visual inspection. Objects were filtered by shape based on roundness and sharpness criteria computed by DAOFIND; we required -1 $<$ *roundness $<$ 1 and 0.2 $<$ *sharpness $<$ 1. Several areas of the frame were masked out, including the galaxy nucleus, the bright saturated star in the PC, and a few bright stars and galaxies in the WFs.**

Aperture photometry was performed using a 1-pixel radius aperture in both the PC and the WFs. The small aperture was used to reduce background related errors, which are large in the central region of the galaxy. Sky values were determined by taking an estimate of the mode in a background annulus; in the PC, the annulus was 2 pixels wide with an inner radius of 12 pixels, while in each WF the corresponding annulus was 3 pixels wide with an inner radius of 6 pixels. We chose to use small sky annuli because the galaxy background is variable on small scales.

Corrections from 1-pixel aperture flux to flux within 0.5$^{\prime \prime}$ were determined from the nine brightest sources. The average aperture corrections were 1.03 and 1.15 magnitudes in the PC frame for F450W and F702W, respectively, and 0.63 and 0.71 magnitudes in the WFs. These are slightly larger than the corrections for bright point sources measured by Holtzman *et al.*  (1995a), but this is expected if these objects are slightly resolved. In the PC, measured aperture corrections varied on the order of 0.05 mag between different objects. Measured variations in the aperture corrections for the individual WFs were $\sim$ 0.02 mag, so an average aperture correction was assumed for all three. If there is an intrinsic spread in cluster sizes, using a single aperture correction for all objects would lead to systematic errors in the derived magnitudes. To first order, however, colors are unaffected by such a spread. The transformations from HST magnitudes (F450W, F702W) to Cousins magnitudes (B, R) were made using the transformations given by Holtzman *et al.*  (1995b).

Foreground Galactic extinction in the direction of NGC 1275 was taken to be $A_B$ = 0.70 (Burstein & Heiles 1984). We determined the extinction in the WFPC2 passbands by numerically integrating the product of an A star spectrum, the filter transmission, the system response, and an extinction curve from Cardelli *et al.*  (1989) derived assuming $R_V=3.1$. An A star was used because of its similarity to the observed spectrum of the bright clusters. This gave $A_{F450W} = 0.69$ and $A_{F702W} =
0.40$. These estimates do not include internal extinction, however, which may be important in NGC 1275 since patchy dust is evident; this is discussed further below.

Potential photometry errors include intrinsic Poisson error in the signal, error in the aperture correction, and error in the background level determination. Typical errors were estimated from the simulated clusters from our completeness tests and are shown in Figure 2. The dashed lines in color and magnitude indicate input magnitudes for the artificial clusters.

Completeness was estimated by generating simulated objects at each of 9 different input magnitudes (23.75 to 27.75 in B). The simulated objects were given the same noise characteristics as real objects. In the PC, 1000 objects were placed randomly in each of 3 separate annuli to better model the variation of completeness with distance from the center of the galaxy. Each set of artificial clusters was given a spatial distribution which mimicked the distribution of the brightest (B $<$ 25) blue and red clusters. In the WFs, the 1000 artificial objects were placed randomly in a uniform distribution. On each chip, a PSF of one of the isolated bright clusters was used to create the simulated objects (the use of other objects gave comparable results). Two sets of artificial objects were placed to estimate the difference in detection efficiency for objects of two different colors; the two sets of simulated clusters were given B-R=0.7 and B-R=1.6, comparable to the colors of the observed cluster populations (next section). The derived average completeness (weighted by the spatial distribution of the observed clusters) as a function of B magnitude is shown in Figure 3 for the whole image as well as for the section which fell in the PC.

We estimated the number of foreground stars based on the models of Bahcall and Soneira (1980) and expect $<$ 4 objects with B $<$ 27 in the PC. Since we expect such small numbers of foreground stars, we have neglected them.

# Results

## Photometry

We present the photometry of all detected objects in NGC 1275 by plotting color against apparent magnitude in Figure 4a. Extinction-corrected absolute magnitudes and colors are shown as well, using a distance of 70.2 Mpc determined by assuming pure Hubble flow with $H_0$ = 75 km/s/Mpc. We find 1181 blue objects (B-R $<$ 1.1) and 1855 red objects. Distinct red and blue populations are seen, although there is significant color scatter. Roughly 50% of the blue objects are seen in the PC, as compared to about 15% of the red objects.

These objects cannot be single stars because they are marginally resolved and because they are too bright. With the exception of the Shields-Filippenko object, the objects lack H$\alpha$ emission (Holtzman *et al.*  1992), which precludes their being emission line objects. Consequently, they are probably star clusters, as supported by spectroscopy of the brightest object (Zepf *et al.*  1995). They are brighter and more compact than open clusters in the Galaxy; the brightest object is 250 times brighter than the brightest Galactic open cluster. Their closest analog seems to be globular clusters, although the inferred ages are much smaller (see below).

To test whether the observed color scatter is likely to be related to internal extinction within NGC 1275 and photometric errors from errors in background subtraction, we show in Figure 4b only clusters in the southwestern half of the PC frame which have estimated photometric errors of less than 0.15 mag in both B and R. This region was chosen because it contains fewer obvious signatures of dust than the northeastern region. In fact, we find that the spread in color among the blue objects in Figure 4b is consistent with that associated with photometric errors (cf. Figure 2), supporting the conclusion that the observed color spread within the two populations is likely to be caused by the variation of internal extinction in the environment of the galaxy, rather than differences in age. We conclude that the blue clusters are a single color, single age population, with the spread in colors coming from Poisson errors and background subtraction errors. Figure 4b shows a distinct bimodality in the cluster population, consistent with the presence of two discrete populations.

The blue clusters have a typical dereddened (B-R)$_0 \sim$ 0.4. This is much bluer than typical old globular clusters, which have (B-R)$_0$ $\sim$ 1.2 (Harris 1996). The uniformity of the colors suggests that the clusters were formed in a relatively short duration single burst rather than an ongoing or continuous process. Figure 5 shows the magnitude and color evolution of a single burst population for a variety of metallicities, assuming a Salpeter initial mass function from 0.1 to 125 M$_\odot$ and a total mass of 1 M$_\odot$ (Bruzual & Charlot 1993). If we compare the measured colors of the blue clusters in the PC to the models, we infer ages of 10$^7$ to 10$^9$ years, allowing for small possible systematic errors in the models and the data as well as the unknown metallicity of the objects. We note that the clusters are far too blue to be explained by a metallicity effect alone. Ages younger than 10$^7$ years are ruled out by the absence of H$\alpha$ emission from the clusters. Spectra of the clusters (Zepf *et al.*  1995) argue for a several hundred Myr old population dominated in the blue by A type stars. It is also simply less probable that we that we should happen to observe these clusters at a younger age.

The LMC has a number of compact clusters with comparable colors (van den Bergh 1981). These objects are substantially fainter than the brightest objects in NGC 1275, implying that they are less massive. Still, observations of blue clusters in the LMC, where we can directly measure ages using the main sequence turnoff in color-magnitude diagrams, give us a less model dependent idea of the spread in colors that could be expected from a spread in ages. Figure 6 shows the distribution of B-V versus age for LMC clusters; ages are from Elson & Fall (1988) and Hodge (1981, 1983), and colors from van den Bergh (1981). Between 10$^8$ and 10$^9$ years, we see that the integrated colors of the clusters cover a range of 0.3 in B-V. This argues against a large spread in the ages for the NGC 1275 clusters, which show no such spread in colors. While our observations were in B and R, not B and V, the color spread in B-R for the same age range would be 0.8 magnitudes, assuming a V-R vs B-R relation from the population synthesis models.

Given ages, we can estimate masses for the clusters based on their luminosities. Age estimates of 10$^8$ and 10$^9$ years give mass estimates for the brightest object of 2 $\times$ 10$^7$ and 10$^8$ M$_\odot$, respectively. Blue clusters at our completeness limit of B $=$ 27 have masses of 2 $\times$ 10$^4$ to 10$^5$ M$_\odot$. These masses are very comparable with those of Milky Way globulars, which typically have masses between 10$^4$ and 3$\times$ 10$^6$ M$_\odot$, with a handful of clusters with M $<$ 10$^4$M$_{\odot}$ (Mandushev, Spassova, & Staneva 1991, Pryor and Meylan 1993). Our mass estimates depend on the assumption of a Salpeter IMF.

Two of the three brightest objects in the southwest half of the PC (Figure 4b) have much redder colors than the rest of the blue clusters (B-R of 0.95 and 1.23). While these objects are close to the nucleus, they are bright enough that errors in determining the galaxy background are not a large source of error. Possible explanations are reddening by dust, a younger age (slightly less than 10$^7$ years old) when the integrated light is dominated by red supergiants (Leitherer & Heckman , 1995), or the identification of these objects with late type foreground stars. For the redder object, the last interpretation is supported by our measurement of sizes ($\S 4.5$), since it appears to be unresolved while almost all other objects appear marginally resolved.

We suggest that the objects with B-R $\sim$ 1.6 ($\sim$ 1.3 dereddened) are the old population of globular clusters around NGC 1275. These are slightly redder than Galactic globular clusters which have an average (B-R)$_0$ $\sim$ 1.2. The brightest (M$_B$ = -10.7) has a magnitude similar to the brightest globulars seen in other central cluster galaxies.

## Luminosity Function

Figure 7a shows the B luminosity function for the blue clusters (B-R $<$ 1.1), and Figure 7b shows the B luminosity function for the red clusters (B-R $>$ 1.1). The dotted lines in each plot show the results after correction for incompleteness. The luminosity function of the blue clusters looks distinctly different from that expected for a typical old globular cluster population at a younger age. While the Galactic globular cluster luminosity function is Gaussian in shape with peak at M$_B$ $\sim$ -6.6 and $\sigma =$ 1.0 mag (Harris 1996), the luminosity function for the clusters in NGC 1275 is more closely fit by a exponential. If we were to observe the current Galactic globular cluster system at an age of 500 million years and at the distance of NGC 1275 (and the same Galactic extinction), we would observe a turnover in the luminosity function at B $\sim$ 24 ($H_0 =$ 75 km/s/Mpc). This prediction assumes 4 magnitudes (13 Gyr) of fading as suggested by the models of Bruzual & Charlot (1993). Our photometry is complete at the 50% level to 3 magnitudes fainter than this, yet there is no evidence of a turnover in the luminosity function of the clusters in NGC 1275.

The luminosity function of the blue clusters has the same trend as that of Galactic open clusters, namely, an increasing number of clusters at fainter magnitudes with no turnover (van den Bergh and LaFontaine 1984). However, the most luminous clusters in NGC 1275 are brighter than any observed Galactic open clusters, despite being older than the brightest of them. The compactness as inferred from the HST images is also different from that of open clusters ($\S 4.5$). Mateo (1993) finds that an increasing luminosity function for clusters in the LMC can be produced by a mixed population of true globulars and open clusters. In NGC 1275, a similar explanation would require a large number of open clusters which are far more massive than any Galactic open clusters. Consequently, as a system, these clusters are different from any cluster system seen in Local Group galaxies. As individual objects, however, the physical properties appear most similar to globular clusters.

In Figure 8, we compare the luminosity function of clusters in NGC 1275 with the luminosity functions in several of the other galaxies where reasonable numbers of blue clusters are found, as well as with the luminosity functions of Milky Way clusters. Dashed lines show rough completeness limits for the different observations, and dotted lines give the location of the expected turnover given the various age estimates of the clusters. There is no strong evidence in any of the young cluster systems for an intrinsic turnover in the luminosity function.

If these objects will evolve to look like the Galactic globular cluster system, some mechanism must preferentially destroy low luminosity, low mass clusters on timescales of 13 Gyr. Two possibilities are tidal disruption and evaporation (Gnedin & Ostriker 1997; Elmegreen & Efremov 1997). For tidal disruption to preferentially destroy low mass clusters, such clusters must have a lower mean density than more massive clusters. Evaporation through two-body relaxation can preferentially destroy lower mass clusters as their stars escape the cluster potential well and are captured by the gravitational potential well of the galaxy. Gnedin and Ostriker (1997) predict that from half to three quarters of the initial globular cluster population in a galaxy may be destroyed by a combination of these processes in a Hubble time.

## Surface Density Profile

Figure 9 shows the surface density distribution of red and blue clusters (open triangles and filled triangles, respectively). Corrections have been made for incomplete annuli at large radii, for masked areas, and for differing completeness levels with distance from the bright galaxy center. Error bars are simple Poisson errors in counts of clusters. The central area, within 18 pixels (0.8 arcsec) from the core of the galaxy, has been omitted from the surface density profile due to the bright nucleus. It is clear that the blue clusters are more centrally concentrated than the red ones.

We fit the surface density distribution of the blue and red cluster systems from 4$^{\prime \prime}$ to 120$^{\prime \prime}$ with power laws. For the blue clusters, we find that a power law slope of -1.3 provides a reasonable fit over the entire range. For the red clusters, we note a flattening in the inner portion of the distribution and fit the surface density with two power laws. Inside 0$^{\prime} \!$.5, we measure a slope of -0.6, while outside 0$^{\prime} \!$.5 we find a slope of -1.0. Flattening in the central part of the surface density distribution of a globular cluster system is also seen in other similar galaxies, such as M87 (McLaughlin 1995) and NGC 5128 (Harris *et al.*  1984). Our measurements for the distribution of old clusters are in good agreement with those of Kaisler *et al.*  (1996). However, they adopt a power law slope of -1.3, based on the observational relation between galaxy magnitude and power law slope of the surface density distribution of the globular cluster system (Harris *et al.*  1986, 1993). We note that our measured outer slope is sensitive to the choice of break point between the inner and outer regions; the data show a sign of steepening even more in the outermost parts. The spatial distribution of the clusters agrees with the galaxy merger model of Ashman and Zepf (1992), which predicts that the younger clusters will be peaked towards the center of the galaxy, while the older clusters will be more flattened in their distribution.

## Specific Frequency

The specific frequency of globular clusters in a galaxy is given by

$$S_N = N_T 10^{0.4(15+M_V)}$$

(Harris and van den Bergh 1981) where $N_T$ is the total number of clusters and $M_V$ is the absolute V magnitude of the galaxy. This is of interest in light of the hypothesis that elliptical galaxies form via mergers of spirals and the observation that ellipticals have a higher specific frequency of globular clusters than spirals. Consequently, the merger hypothesis requires that globular clusters form during mergers with a higher specific frequency than that of the progenitor population. Also, it requires that clusters form in numbers comparable to those of the old cluster population (Zepf & Ashman 1993). To address these questions, we calculate the specific frequency of the old population, the current specific frequency of the merger related population, and the predicted specific frequency of the merger related population in 13 Gyr.

To calculate specific frequencies, we need to estimate the total number of clusters associated with each population. These estimates come from integrating the surface density profiles of the red and blue clusters from the galaxy core to 100 kpc. This integration suggests that the total number of blue clusters is five times the number of clusters that we observe. For the red clusters, we correct this number for the clusters fainter than our completeness cutoff based on a Gaussian fit to the part of the magnitude distribution which we have observed. For the blue clusters, we make no correction for missing faint clusters and consider the resultant specific frequency to be a lower limit. To make a rough prediction of the number of clusters from the merger related population which survive to be 13 Gyr old, we assume that the cluster luminosity function will come to resemble the globular cluster luminosity function of the Galaxy after 4 magnitudes of fading and substantial cluster destruction. We assume that no clusters brighter than the predicted turnover of M$_B$ = -10.6 are destroyed, and that the remaining clusters fainter than the turnover are equal in number to those brighter than the turnover, forming a Gaussian luminosity function. We note that this would require the destruction of $\sim$ 90% of the young clusters in NGC 1275 seen within our field. If the young cluster population is older than 10$^8$ yrs, less cluster destruction would be required, and we would derive a higher specific frequency. Of course, if some fraction of the massive clusters are destroyed as well, we would derive a lower specific frequency. Our corrections yield estimates of 12,700 clusters for the old population, 5700 clusters for the young population, and 550 clusters remaining from the young cluster population after 13 Gyr.

We also need to estimate what fraction of the underlying stellar component is associated with each of the populations. NGC 1275 clearly has an excess blue population over what is normally observed in central cluster galaxies (Romanishin 1987). It is unknown whether this blue population is related to star formation in a recent merger or whether it is related to the presence of the cooling flow centered on NGC 1275. Assuming the entire excess blue light comes from a merger-related population, we derive a lower limit for the specific frequency of globular clusters in such a population. Romanishin (1987) estimates that 15% of the total B light comes from an excess blue population, and we adopt this contribution as an upper limit for the V band as well. With a total brightness of V=11.88 (deVaucouleurs *et al.*  1991), our adopted $A_V=0.53$, and $H_0=75$ km/s/Mpc, we derive $M_V^{tot}=-22.89$. Using a 15% contribution for the merger-related population gives $M_V^{young}=-20.83$ and $M_V^{old}=-22.71$. In order to estimate the specific frequency of the merger-related population *as it would appear in $\sim$ 13 Gyr*, we need to allow for a fading of the background population, and we adopt 4 magnitudes of fading based on the population synthesis models.

We combine the information on the total number of clusters with the total V magnitudes for the populations to calculate specific frequencies. These estimates yield a $S_N\sim 10$ for the old population, a $S_N\sim 27$ for the current young population (without corrections for background fading or for cluster destruction), and a $S_N\sim 102$ as an estimate of the specific frequency of the young population after 13 Gyr of cluster destruction and 4 magnitudes of background fading. Even without correction for clusters outside the WFPC2 field of view, we would predict 150 clusters surviving for 13 Gyr and a specific frequency of 28. All of these numbers, however, have significant uncertainties, as the corrections for completeness, cluster destruction, and fading are very substantial while at the same time being highly uncertain. From shallower ground-based observations over a larger field, Kaisler *et al.*  (1996) determine the S$_N$ for the old cluster population to be 4.3. The main sources of disagreement are that they have adopted a steeper surface density profile and a different total V magnitude.

The increase in the *overall specific frequency of the galaxy from this event, however, is small, as the number of these clusters predicted to survive 13 Gyr is small; only about 550 new clusters would be added to a population of 12,700. Whether the specific frequency increases significantly in a merger depends on the mass and gas content of the merging bodies. However, the numbers are intriguing because they suggest that, even allowing for substantial cluster destruction, the specific frequency of a merger-related population in 13 Gyr will be much higher than that typically observed in spiral galaxies. While this particular merger event will not have a large impact on the overall specific frequency in NGC 1275, it does show that a merger population can have a large specific frequency, which is the first requirement for increasing the overall specific frequency in a merger. Consequently, the observation that ellipticals have a higher specific frequency of globular clusters is not inconsistent with the hypothesis that these galaxies have formed via the merging of spiral systems.*

Our calculation indicates that an estimate of the future specific frequency is difficult to make, as several large corrections must be applied. The most uncertain of these is probably the estimate of the number of clusters which will be destroyed. Future work will need to carefully consider the probability of cluster destruction as a function of cluster mass and density if one wants to put serious constraints on the merger hypothesis for ellipticals from observations of their globular cluster systems.

## Sizes

Holtzman *et al.*  (1996) estimated the sizes of marginally resolved clusters in WFPC2 images using the difference between a 1 and a 2 pixel aperture magnitude. These measurements are compared with measurements of simulated objects made by convolving a model PSF with a modified Hubble profile. The comparison yields an estimate of the cluster core radius *if* the cluster shape is well modelled by a modified Hubble profile. Figure 10 shows measurements of $m_1-m_2$ in F702W for the NGC 1275 objects; note that values are near zero because $m_1-m_2$ is defined differentially with respect to a model point source placed at the same location on the frame as each cluster (see Holtzman *et al.*  1996 for details). Positive values of $m_1-m_2$ indicate that the objects are resolved, and core radii are estimated using Figure 4 of Holtzman *et al.*  (1996). For the brighter objects, we find that the objects are most likely resolved, with estimated core radii of $\sim 0.01-0.05$ pixels, corresponding to $\sim$ 0.15-0.75 pc using our adopted distance to NGC 1275. Despite the uncertainties in these measurements, it seems clear that these objects are compact, and much more comparable to globular clusters than to open clusters.

# Conclusions

From new deep observations of NGC 1275 with HST, we identify roughly 3000 objects which appear to be compact star clusters. The color distribution of the cluster system is bimodal, with a blue population which has (B-R)$_0$ $\sim$ 0.4 and a red population which has (B-R)$_0$ $\sim$ 1.3. We suggest that the red objects are members of the old globular cluster system and that the blue objects are members of a young globular cluster system. In an apparently dust-free region, the spread in the colors of the blue clusters is small enough that it is entirely attributable to scatter from errors in the photometry, mostly due to errors in sky subtraction. This suggests that the blue clusters are a single color, single age population. This argues against the hypothesis that clusters have been forming continuously from the cooling flow, and supports the hypothesis that their formation may have been triggered by a previous merger.

The luminosity function continues to rise to the limit of our observations and is inconsistent with a Gaussian globular cluster luminosity function peaking near M$_V$ = -7.3 after correction for evolutionary effects. A similar luminosity function observed in NGC 4038/9 (the Antennae) has been used by van den Bergh (1995) to argue that the young clusters observed there are not true globular clusters. However, the masses and sizes of the individual young clusters appear to be comparable to those of globulars. We suggest that either the luminosity function evolves, with fainter clusters being preferentially destroyed as time passes, or that the initial luminosity/mass function of the young cluster system is different from that of typical old globular cluster systems.

If cluster destruction is responsible for the difference in luminosity function, a destruction mechanism which preferentially destroys fainter, lower mass clusters must be invoked. In order to create a turnover in the luminosity function at the same mass as in old cluster systems, the destruction of over 90% of the clusters currently seen in NGC 1275 would be required, given our assumption of 4 magnitudes of fading.

Given these likely differences in the luminosity functions of different globular cluster systems, whether caused by evolutionary effects or by initial differences, it is clear that some caution must be exercised in using the cluster luminosity function as a distance indicator.

We have attempted to estimate the specific frequency of the young population seen in NGC 1275 to determine whether it is likely that the overall specific frequency as seen at a future time could be increased because of the proposed merger-related formation of globulars. We note that the overall specific frequency of NGC 1275 will only be slightly increased by the proposed merger event which formed the blue clusters, due to the small number of clusters expected to survive for a Hubble time. But we estimate a large specific frequency of the merger related population, which suggests that mergers are efficient in forming globular clusters. However, corrections for incompleteness at faint magnitudes, clusters residing outside the current field of view, uncertainties in the mass and luminosity of the stellar population formed during a merger event, and, especially, corrections for evolutionary effects in both the background and cluster populations are all significant effects for which we have only relatively crude estimates.

This work was supported in part by NASA under contract NAS7-918 to JPL and a grant to M. C. from the New Mexico Space Grant Consortium.

[^1]: Based on observations with the NASA/ESA *Hubble Space Telescope*, obtained at the Space Telescope Science Institute, operated by AURA Inc under contract to NASA