abstract: I report recent results on hidden charm spectroscopy from Belle. These include: observation of a near-threshold enhancement in the $\omega J/\psi$ invariant mass distribution for exclusive $B\rightarrow K \omega J/\psi$ decays; evidence for the decay $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$, where the $\pi^{+}\pi^{-}\pi^{0}$ invariant mass distribution has a strong peak between 750 MeV and the kinematic limit of 775 MeV, suggesting that the process is dominated by the sub-threshold decay $X\rightarrow\omega J/\psi$; and the observation of a peak near 3940 MeV in the $J/\psi$ recoil mass spectrum for the inclusive continuum process $e^{+} e^{-}\rightarrow J/\psi X$. The results are based on a study of a 287 fb$^{-1}$ sample $e^{+} e^{-}$ annihilation data collected at center-of-mass energies around the $\Upsilon(4S)$ in the Belle detector at the KEKB collider.
address: Department of Physics & Astronomy, University of Hawai'i at Manoa, Honolulu, HI 96822, USA
author: Stephen L. Olsen
title: New particles from Belle

.

# Introduction

The recent surge in activity in hadron spectroscopy and, I suppose, the main motivation for the formation of the Topical Group on Hadron Physics, is the result of renewed interest in a rather old question: *are there hadronic states with a more complex structure than the simple $q\bar{q}$ mesons and $qqq$ baryons of the original quark model*? This revival of interest has been driven by experimental reports of pentaquarks , the narrow $D_{sJ}$ states , and the $X(3872)$ .

In spite of considerable theoretical and experimental effort, the existence of non-$q\bar{q}$ mesons and/or non-$qqq$ baryons remains an open question. While the identification of a strangeness=+1 (or charm=-1) baryon would be definitive evidence for a non-$qqq$ baryon, the experimental situation regarding the existence of such states remains unsettled (and a major topic of discussion at this meeting ). On the other hand, while the $D_{sJ}$ and and $X(3872)$ are experimentally well established, the theoretical interpretation is not so clear. The $D_{sJ}$ states could be standard $P$-wave $c\bar{s}$ states and their narrowness is only surprising because the relativistic potential model calculations that predicted them to be heavier (and above $DK$ threshold) are wrong . Some theorists, including our opening speaker , remain hopeful that a $c\bar{c}$ charmonium assignment can be found for the $X(3872)$.

To sort this all out, I think that the so-called hidden charm mesons can and will play a decisive role for reasons that include:

-   the theory for these systems is well founded (and recently blessed by this year's Nobel Prize Committee) and has fewest ambiguities;

-   the experimental signatures tend to be clean;

-   $c\bar{c}$ meson states below open-charm threshold are narrow and do not overlap; and

-   lots of non-$c\bar{c}$-type mesons have been conjectured, including $D\bar{D^*}$ molecules  and $c\bar{c}$-$gluon$ hybrids .

Although the Belle detector  is specialized to studies of CP violation in $B$ meson decays, it has proven to be a useful device for studying particles containing $c\bar{c}$ pairs. Belle detects $c\bar{c}$ systems produced via weak decays of $b$ quarks —the $b\rightarrow c\bar{c}s$ process is a dominant $b$-quark decay mode— and the continuum production process $e^{+} e^{-}\rightarrow c\bar{c}c\bar{c}$, which has been found to be surprisingly large. The KEKB asymmetric energy $e^{+} e^{-}$ collider  operates at a center-of-mass (cms) energy corresponding to the $\Upsilon(4S)$ resonance and routinely delivers luminosities that are in excess of $10^{34}$cm$^{-2}$s$^{-1}$, thereby providing Belle with a huge data sample that contains about 300 million $B\bar{B}$ meson pair events and over one billion $e^{+} e^{-}\rightarrow q\bar{q}$ continuum annihilation events.

Belle results in the hidden charm meson sector include first observations of:

-   the $\eta_{c}^{\prime}$ via the sequence $B\rightarrow K\eta_{c}^{\prime}$, $\eta_{c}^{\prime}\rightarrow K_{S}K\pi$ ;

-   anomalously large cross sections for the exclusive process $e^{+} e^{-}\rightarrow J/\psi\eta_{c}$ and the inclusive process $e^{+} e^{-}\rightarrow J/\psi(c\bar{c})$ ;

-   the $X(3872)$ meson ;

-   a near-threshold $\omega J/\psi$ mass enhancement in exclusive $B\rightarrow K\omega J/\psi$ decays ; and

-   a peak at 3940 MeV in the $J/\psi$ recoil mass spectrum in the inclusive $e^{+} e^{-}\rightarrow J/\psi X$ process .

In this talk I will discuss the last two items as well as recent results on properties of the $X(3872)$. I will not have time to cover any of the many other Belle results on hadron spectroscopy, such as our many interesting results on charmed baryon spectroscopy , $D^{**}$  and $D_{sJ}$ mesons  and two-photon physics . In addition, I will not have time to report on Belle's lack of observation of pentaquarks  or the $D_{sJ}(2632)$ . All unpublished numbers reported here are preliminary.

# A near-threshold $\omega J/\psi$ mass enhancement in $B \rightarrow K \omega J/\psi$ decays

At the $\Upsilon(4S)$, $B\bar{B}$ meson pairs are produced with no accompanying particles. As a result, each $B$ meson has a total cms energy that is equal to $E_{\rm beam}$, the cms beam energy. We identify $B$ mesons using the beam-constrained $B$-meson mass $M_{\rm bc}=\sqrt{E_{\rm beam}^2 - p_B^2}$ and the energy difference $\Delta E= E_{\rm beam} - E_B$, where $p_B$ is the vector sum of the cms momenta of the $B$ meson decay products and $E_B$ is their cms energy sum. For the final states discussed here, the experimental resolutions for $M_{\rm bc}$ and $\Delta E$ are approximately 3 MeV and 13 MeV, respectively.

We select $B\rightarrow K\pi^{+}\pi^{-}\pi^0 J/\psi$ candidate events ($J/\psi\rightarrow\ell^{+}\ell^{-})$ track combinations with $M_{\rm bc}$ and $\Delta E$ values that are within $2.5\sigma$ of their nominal values. Figure  shows a scatterplot of $M(\pi^{+}\pi^{-}\pi^0J/\psi)$ (vertical) $versus$ $M(\pi^{+}\pi^{-}\pi^0)$ for selected events in the $\Delta E$-$M_{\rm bc}$ signal region. Here a distinct vertical band corresponding to $\omega\rightarrow\pi^{+}\pi^{-}\pi^0$ decays is evident near $M(\pi^{+}\pi^{-}\pi^0) = 0.782$ GeV.

We identify three-pion combinations with $M(\pi^{+}\pi^{-}\pi^0)$ within 25 MeV of $m_{\omega}$ as $\omega$ candidates and form the Dalitz plot of $M^2(\omega J/\psi)$ (vertical) $versus$ $M^2(\omega K)$ (horizontal) shown in Fig. . The clustering of events near the left side of the plot corresponds to $B\rightarrow K_XJ/\psi$; $K_X \rightarrow K\omega$ events, where $K_X$ denotes strange meson resonances such as $K_1(1270)$, $K_1(1400)$, and $K^*_2(1430)$ that are known to decay to $K\omega$. There is also a clustering of events with low $\omega J/\psi$ invariant masses near the bottom of the Dalitz plot. To study these, we suppress $K_X\rightarrow K\omega$ events by only looking at events in the region $M(K\omega)>1.6$ GeV, to the right of the dashed line in Fig. .

The $M_{\rm bc}$ and $\Delta E$ distributions of the selected events indicate that about half of the entries in the $M(\omega K)>1.6$ GeV Dalitz plot region are due to background. To perform a background subtraction and determine the level of $B\rightarrow K\omega J/\psi$ signal events, we separate the data into $40$ MeV-wide bins of $M(\omega J/\psi)$ and measure the $B$ meson signal levels in the $M_{\rm bc}$ and $\Delta E$ distributions. The histograms in Fig.  show the $M_{\rm bc}$ distributions for the twelve lowest $M(\omega J/\psi)$ mass bins, where strong peaks at $M_{\rm bc}=m_B$ are evident at low $\omega J/\psi$ masses, especially for the mass regions covered by Figs. (b) and (c). The corresponding $\Delta E$ distributions (not shown) show similar structure. We establish the $B\rightarrow K\omega J/\psi$ signal level for each $M(\omega J/\psi)$ mass bin by performing binned fits simultaneously to the $M_{\rm bc}$ and $\Delta E$ distributions with Gaussian functions for the signal and smooth background functions. The smooth curves in Fig.  indicate the fitted $M_{\rm bc}$ curves for each $\omega J/\psi$ mass bin.

The bin-by-bin signal yields are plotted $vs$ $M(\omega J/\psi)$ in Fig. . An enhancement is evident around $M(\omega J/\psi)=3940$ MeV. The curve in Fig.  is the result of a fit with a $S$-wave Breit Wigner function threshold function of the form $f(M) = A_0 q^*(M)$, where $q^*(M)$ is the momentum of the daughter particles in the $\omega J/\psi$ restframe. This functional form accurately reproduces the threshold behavior of Monte Carlo simulated $B\rightarrow K\omega J/\psi$ events that are generated uniformly distributed over phase-space.

The fit gives a Breit-Wigner signal yield of $58\pm 11$ events with a peak peak mass and total width of $$\begin{aligned}
M       &=&    3943 \pm 11 {\rm (stat)} \pm 13 {\rm (syst)}~{\rm MeV}\\
\Gamma  &=&     87  \pm 22 {\rm (stat)} \pm 26 {\rm (syst)}~{\rm MeV},
\end{aligned}$$ where the systematic errors are determined from variations in the values when different bin sizes, fitting shapes and selection criteria are used. The event yield translates into a product branching fraction (here we denote the enhancement as $Y(3940)$): $${\cal B}(B\rightarrow K Y(3940)){\cal B}(Y(3940)\rightarrow\omega J/\psi) =
(7.1 \pm 1.3 {\rm (stat)} \pm 3.1 {\rm (syst)})\times 10^{-5},$$ The statistical significance of the signal, determined from $\sqrt{-2\ln({\mathcal L}_0/{\mathcal L}_{\rm max})}$, where ${\mathcal L}_{\rm max}$ and ${\mathcal L}_0$ are the likelihood values for the best-fit and for zero-signal-yield, respectively, is $8.1\sigma$.

A $c\bar{c}$ charmonium meson a mass of 3943 MeV would dominantly decay to $D\bar{D}$ and/or $D\bar{D}^*$; hadronic charmonium transitions should have minuscule branching fractions. On the other hand, decays of $c\bar{c}$-$gluon$ hybrid charmonium to $D^{(*)}\bar{D^{(*)}}$ meson pairs are forbidden or suppressed, and the relevant "open charm" threshold is $m_D + m_{D^{**}} \simeq 4285$ MeV , where $D^{**}$ refers to the $J^P = (0,1,2)^+$ charmed mesons. Thus, a hybrid state with a mass equal to that of the peak we observe would have large branching fractions for decays to $J/\psi$ or $\psi '$ plus light hadrons . Moreover, lattice QCD calculations have indicated that partial widths for such decays can be comparable to the width that we measure . However, these calculations predict masses for these states that are between 4300 and 4500 MeV , substantially higher than our measured value.

# The $X(3872)$ with 253 fb$^{-1}$

The $X(3872)$ was discovered by Belle as a narrow $\pi^{+}\pi^{-}J/\psi$ mass peak in exclusive $B^-\rightarrow K^-\pi^{+}\pi^{-}J/\psi$ decays . Figure  shows the $X(3872)$ signal from a 253 fb$^{-1}$ data sample containing 275 million $B\bar{B}$ pairs. The observed mass and the narrow width are not compatible with expectations for any of the as-yet unobserved charmonium states . Moreover, the $\pi^{+}\pi^{-}$ invariant mass distribution, shown in Fig. , peaks near the upper kinematic limit of $M(\pi^{+}\pi^{-})=775$ MeV, and has a shape that is consistent with $\rho\rightarrow\pi^{+}\pi^{-}$ decays. Charmonium decays to $\rho J/\psi$ final states violate isospin and are expected to be suppressed. The $X(3872)$ and its above-listed properties were confirmed by the BaBar , CDF  and D0  experiments.

The $X(3872)$ mass ($3871.9\pm 0.5$ MeV ) is within errors of the $D^0\bar{D^{0*}}$ threshold ($3871.3\pm
1.0$ MeV ); the difference is $0.6\pm 1.1$ MeV. This has led to speculation that the $X$ might be a $D^0\bar{D^{0*}}$ bound state . According to ref. , the preferred quantum numbers for such a bound state would be either $J^{PC} = 0^{-+}$ or $1^{++}$. The decay of an $C=+1$ state to $\pi^{+}\pi^{-}J/\psi$ would proceed via an $I=1$ $\rho^0J/\psi$ intermediate state and produce the $\pi^{+}\pi^{-}$ mass spectrum like that we see. In this meson-meson bound state interpretation, the close proximity of the $X$ mass to $D^0\bar{D^{0*}}$ threshold compared to the $D^+ D^{-*}$-$D^0\bar{D^{0*}}$ mass splitting of 8.1 MeV produces a strong isospin violation.

Swanson made a dynamical model for the $X(3872)$ as a $D^0\bar{D^{0*}}$ hadronic resonance . In this model, $J^{PC}=1^{++}$ is strongly favored and the wave function has, in addition to $D^0\bar{D^{0*}}$, an appreciable admixture of $\omega J/\psi$ plus a small $\rho J/\psi$ component. The latter produces the $\pi^{+}\pi^{-}J/\psi$ decays that have been observed; the former gives rise to $\pi^{+}\pi^{-}\pi^{0}J/\psi$ decays via a virtual $\omega$ that are enhanced because of the large $\omega J/\psi$ component to the wavefunction. Swanson's model predicts that $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ decays should occur at about half the rate for $\pi^{+}\pi^{-}J/\psi$ and with a $\pi^{+}\pi^{-}\pi^{0}$ invariant mass spectrum that peaks near the upper kinematic boundary of 775 MeV (7.5 MeV below the $\omega$ peak).

$X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ decays would populate the horizontal band indicated by the horizontal lines in the scatterplot of Fig. . This corresponds to the $\pm 3\sigma$ band $|M(\pi^{+}\pi^{-}\pi^{0}J/\psi)-m_{X(3872)}|<16.5$ MeV.

Figure  shows the $M_{\rm bc}$ distributions for events that are in the $\Delta E$ and $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal regions for 25 MeV-wide $\pi^{+}\pi^{-}\pi^{0}$ invariant mass bins; Fig.  shows the corresponding $\Delta E$ distributions for events in the $M_{\rm bc}$ and $X$ signal regions. There are distinct $B$ meson signals in both the $M_{\rm bc}$ and $\Delta E$ distributions for the $M(\pi^{+}\pi^{-}\pi^{0})>750$ MeV bin and no evident signals for any of the other $3\pi$ mass bins. The curves in the figures are the results of binned likelihood fits that are applied simultaneously to the $M_{\rm bc}$ and $\Delta E$ distributions.

Figure  shows the fitted $B$-meson signal yields $vs$ $M(\pi^{+}\pi^{-}\pi^0)$. All of the fitted yields are consistent with zero except for the $M(\pi^{+}\pi^{-}\pi^{0})>750$ MeV bin, where the fit gives $12.4\pm 4.1$ events. The statistical significance of the signal in this one bin, determined from $\sqrt{-2\ln({\mathcal L}_0/{\mathcal L}_{\rm max})}$, where ${\mathcal L}_{\rm max}$ and ${\mathcal L}_0$ are the likelihood values for the best-fit and for zero-signal-yield, respectively, is $6.6\sigma$.

Figure  shows the $M(K\pi^{+}\pi^{-}\pi^{0})$ distribution for for the $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal events. The distribution is spread across the limited allowed kinematic region and there is no evident structure that might be producing the high mass peak in Fig.  by some sort of a kinematic reflection.

A possible background to the observed signal would be feed-across from the near-threshold $\omega J/\psi$ enhancement in $B\rightarrow K\omega J/\psi$ decays described above. Since the $\omega\rightarrow\pi^{+}\pi^{-}\pi^{0}$ resonance peak is at $m_{\omega}=782.5$ MeV, which is 7.5 MeV above the maximum possible $3\pi$ invariant mass value for $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ decays, there is no overlap between the centroids of the $\omega J/\psi$ and $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal bands in Fig. . However, there is some overlap in the tails of the kinematically allowed regions for the two processes that might result in some events from one signal feeding into the other.

We determine the level of signal cross-talk to be $0.75\pm0.14$ events from the integral of the fitted function in Fig.  over the region of overlap with the $X(3872)$ signal band. As an independent check, we refitted for the $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal yield with a tighter restriction on $M(\pi^{+}\pi^{-}\pi^{0}J/\psi)$, namely $m_X - 3\sigma < M(\pi^{+}\pi^{-}\pi^{0}J/\psi) < m_X + 1\sigma$, that has *no overlap* with the $\omega$ band. The $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal yield in the truncated region is $10.6\pm 3.6$ events. For a Gaussian signal distribution with no feed-across background, we expect the truncation of the signal region to reduce the signal by 2.1 events (16%); the observed reduction of 1.8 events is consistent with a feed-across level that is less than one event.

Another possible source of background to the $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal are non-resonant $B^-\rightarrow K^-\pi^{+}\pi^{-}\pi^{0}J/\psi$ decays. To determine the level of these, we looked for $B$-meson signals in the $M_{\rm bc}$-$\Delta E$ distributions for events in $M(\pi^{+}\pi^{-}\pi^{0}J/\psi)$ sidebands above and below the $X(3872)$ mass region. There is no evidence for significant signal yields in the $M_{\rm bc}$-$\Delta E$ distributions of either sideband. Fits gives an estimate of the non-resonant background in the $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal bin of $1.3\pm 1.0$ events.

To determine the branching fraction, we attribute all of the signal events with $M(\pi^{+}\pi^{-}\pi^{0})>750$ MeV to $X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ decay. We compute the ratio of $\pi^{+}\pi^{-}\pi^{0}J/\psi$ and $\pi^{+}\pi^{-}J/\psi$ branching fractions by comparing this to the number of $X\rightarrow\pi^{+}\pi^{-}J/\psi$ in the same data sample, corrected by MC-determined relative detection efficiencies. The ratio of branching fractions is $$\frac{{\cal B}(X\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi)}{{\cal B}(X\rightarrow\pi^{+}\pi^{-}J/\psi)}=
\frac{N_{ev}(\pi^{+}\pi^{-}\pi^{0}J/\psi)\varepsilon_{\pi^{+}\pi^{-}J/\psi}}
{N_{ev}(\pi^{+}\pi^{-}J/\psi)\varepsilon_{\pi^{+}\pi^{-}\pi^{0}J/\psi}}
= 1.1\pm0.4{\rm (stat)}\pm 0.3 {\rm (syst)},$$ where the systematic error reflects the uncertainty in the relative acceptance, the level of possible feed-across and nonresonant background, and possible event loss due to the $M(3\pi)>750$ MeV requirement, all added in quadrature. If we allow for cross-talk and non-resonant contributions at their maximum ($+1\sigma$) values, the statistical significance of the $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$ signal is reduced to $\simeq 4\sigma.$

# A new charmonium state in inclusive $e^{+} e^{-}\rightarrow J/\psi X$ annihilations.

Some of the biggest surprises from Belle have nothing to do with $B$-meson physics at all and have come, instead, from the inclusive $e^{+} e^{-}\rightarrow J/\psi X$ annihilation process. This is demonstrated in Fig. , which shows the distribution of masses for systems with more than two charged tracks that recoil against $J/\psi$ mesons produced in the $e^{+} e^{-}$ continuum at or near the $\Upsilon(4S)$ resonance. In this figure, which is based on a 280 fb$^{-1}$ data sample, the histogram indicates the background level derived from the $J/\psi\rightarrow\ell^{+}\ell^{-}$ mass sidebands.

The prominent peak at $M_{\rm recoil} \simeq 2.98$ GeV in Fig.  corresponds to the $\eta_{c}$. From the yield of events we determine a cross-section branching-fraction product  $$\sigma_{\rm Born}(e^{+} e^{-}\rightarrow J/\psi\eta_{c}){\cal{B}}(\eta_{c}\rightarrow>2tracks) = 25.6 \pm 4.4~{\rm fb.}$$ This is more than an order of magnitude higher than non-relativistic QCD (NRQCD) calculations of $\sim2$ fb$^{-1}$. There is no evident signal for any recoils with mass below $M_{\eta_{c}}$, which is also the $c\bar{c}$ mass threshold. Also contrary to NRQCD expectations, the four-charmed-quark process $e^{+} e^{-}\rightarrow c\bar{c}c\bar{c}$ dominates inclusive $J/\psi$ production. From the total number of charmonium states and charmed particles found in the recoil system, we determine the cross section ratio  $$\frac{\sigma(e^{+} e^{-}\rightarrow J/\psi(c\bar{c})}{\sigma(e^{+} e^{-}\rightarrow J/\psi X)} = 0.82\pm 0.21;$$ NRQCD predicts this ratio to be $\sim 0.1$ .

The second and third prominent peaks in Fig.  are at the masses of the $\chi_{c0}$ and $\eta_{c}^{\prime}$, respectively. The fourth peak is well fitted by a Gaussian function with a peak mass of $3940\pm 12$ MeV and a signal significance of $4.5\sigma$. The width of this state is consistent with experimental $M_{\rm recoil}$ resolution. Since this is rather poor, we can only derive an upper limit on the total width of $\Gamma < 96$ MeV (90% CL).

We investigated this peak further by studying events where a $D$ meson is identified in the $J/\psi$ recoil system, *e.g.* in events of the type $e^{+} e^{-}\rightarrow J/\psi D X$. Figure  shows the distribution of masses recoiling against the $J/\psi$ for $e^{+} e^{-}\rightarrow J/\psi D X$ events where $M_X = m_{D^*}$ (top) and $M_X = m_{D}$ (bottom). There is an evident $9.9\pm 3.3$ event signal for the $3940$ MeV state in the $D\bar{D^*}$ mass spectrum, with a statistical significance of $4.5\sigma$. The signal level is the $D\bar{D}$ mass spectrum is $4.1\pm 2.2$ events with a significance of only $2.1\sigma$.

This peak cannot be identified with any known charmonium state. An obvious guess is that it is either the $\chi_{c0}'$ or the $\eta_{c}''$. However, $\chi_{c0}'\rightarrow D\bar{D^*}$ is forbidden and, thus, ruled out. Likewise $\eta_{c}''\rightarrow D\bar{D}$ decays are also forbidden, but, since the $D\bar{D}$ "signal" is ambiguous, we can't use this to rule out this assignment. On the other hand, an $\eta_{c}''$ assignment to the observed peak would imply a $m_{\psi(3S)} - m_{\eta_{c}(3s)}$ mass splitting of $\sim100$ MeV, about twice as large as the measured splitting for the $2S$ states. This seems unlikely.

The mass of this fourth peak is very similar to that of the $\omega J/\psi$ peak seen in $B\rightarrow K\omega J/\psi$ and described above, and a search for it in the $\omega J/\psi$ decay channel is in progress. In addition, we are examining $B\rightarrow K D\bar{D^*}$ decays for a $D\bar{D^*}$ component of the $\omega J/\psi$ enhancement.

# Summary

We observe peaks near $3940$ MeV in the $\omega J/\psi$ mass distribution from $B\rightarrow K\omega J/\psi$ decays and in the recoil mass spectrum in the inclusive annihilation process $e^{+} e^{-}\rightarrow J/\psi X$. The latter peak is also seen in the exclusive process $e^{+} e^{-}\rightarrow J/\psi D\bar{D^*}$ and, thus, cannot be assigned to the $\chi_{c0}$ charmonium state. At this stage, we cannot tell whether or not the state seen in $B$ decays and the one seen in inclusive $J/\psi$ production are one and the same. Further investigation is in progress.

We observe a $\sim 4\sigma$ signal for $X(3872)\rightarrow\pi^{+}\pi^{-}\pi^{0}J/\psi$. This is the first measurement of an $X$ decay mode other than $\pi^{+}\pi^{-}J/\psi$. The $\pi^{+}\pi^{-}\pi^{0}$ invariant masses are strongly clustered above 750 MeV, near the upper kinematic boundary; this is suggestive of a sub-threshold decay via a virtual $\omega J/\psi$ intermediate state. Such a decay, at near the measured branching fraction, was predicted by Swanson based on a model where the $X(3872)$ is considered to be primarily a $D^0\bar{D^{0*}}$ hadronic resonance .

The presence of the $X(3872)\rightarrow\omega J/\psi$ decay process would establish the Charge-Conjugation quantum number of the $X(3872)$ as $C=+1$. This, in turn, would mean that the $\pi^{+}\pi^{-}$ system in $X\rightarrow\pi^{+}\pi^{-}J/\psi$ decay comes from the decay of a $\rho$ meson. The large isospin violation implied by the near equality of the $\rho J/\psi$ and $\omega J/\psi$ decay widths is difficult to accomodate in a $c\bar{c}$ charmonium interpretation of the $X$, but a natural consequence of the meson-meson bound state model point of view.

# Acknowledgements

I thank the organizers of the APS Topical Group on Hadron Physics for inviting me to give this talk. My Belle colleagues and I thank the KEKB group for the excellent operation of the accelerator, the KEK cryogenics group for the efficient operation of the solenoid, and the KEK computer group and the NII for valuable computing and Super-SINET network support. We acknowledge support from MEXT and JSPS (Japan); ARC and DEST (Australia); NSFC (contract No. 10175071, China); DST (India); the BK21 program of MOEHRD and the CHEP SRC program of KOSEF (Korea); KBN (contract No. 2P03B 01324, Poland); MIST (Russia); MESS (Slovenia); Swiss NSF; NSC and MOE (Taiwan); and DOE (USA).

# References


author: E. Puga; S. Hony; C. Neiner; A. Lenorzer; A.-M. Hubert; L.B.F.M. Waters; F. Cusano; V. Ripepi
bibliography: S284_Puga_bib.bib
date: Received ; accepted
subtitle: Searching for evidence of triggered star formation in an isolated region in the outer Milky Way
title: Spitzer/IRAC view of Sh 2-284 [^1]

Using Spitzer/IRAC observations of a region to be observed by the CoRoT satellite, we have unraveled a new complex star-forming region at low metallicity in the outer Galaxy. We perform a study of S284 in order to outline the chain of events in this star-forming region. We used four-band Spitzer/IRAC photometry as well as H$\alpha$ imaging obtained with INT/WFC. Combining these data with the optical photometry obtained in the frame of CoRoTs preparation and the 2MASS catalog we analysed the properties and distribution of young stellar objects (YSOs) associated with point-like sources. We also studied the SEDs of regions of extended emission, complementing our dataset with IRAS and MSX data. We find that S284 is unique in several ways: it is very isolated at the end of a spiral arm and both the diffuse dust and ionized emission are remarkably symmetric. We have partially resolved the central clusters of the three bubbles present in this region. Despite the different scales present in its multiple-bubble morphology, our study points to a very narrow spread of ages among the powering high-mass clusters. In contrast, the particular sawtooth structure of the extended emission at the rim of each ionized bubble harbours either small lower-mass clusters with a younger stellar population or individual young reddened protostars. In particular, triggered star formation is considered to be at work in these regions. \]

# Introduction

Several scenarios have been put forward to explain the role of expanding H  regions in the frame of triggered star formation . Observationally, it is difficult to find an isolated H  region in the Galactic plane where the impact of the expansion can be decoupled from severe crowding effects . However, triggered SF is commonly invoked when H  regions and bubbles are aligned and a multiple-bubble morphology is detected. Several examples of H  regions have been recently investigated a few degrees off the Galactic disk plane as an illustration of the specific hypothesis called: "collect and collapse mechanism" . The outer region of our Galaxy offers a benchmark to study this type of object . Away from the central spiral arms, the third quadrant is sparsely populated with the giant molecular clouds belonging to the Perseus and Norma-Cygnus arms. Moreover, lower metallicities are measured toward the outer Galaxy , guarantee of closer conditions to the pristine initial molecular clouds.

Sh 2-284 (hereafter S284) is a diffuse nebula observed for the Sharpless (1959) catalog. This region is located in the anti-centre of the Galactic disk (l=212, b=$-$`<!-- -->`{=html}1.3), and it is relatively isolated. It shows a very symmetric morphology in the red plate of the Palomar Observatory Sky Survey (POSS) with an approximated diameter of $80\arcmin$ in the sky. The kinematical study of the H$\alpha$ line in S284 revealed a two-component velocity profile . The radial velocities of these components are 43.0$\pm$`<!-- -->`{=html}0.5 and $-$`<!-- -->`{=html}20.1$\pm$`<!-- -->`{=html}0.5 km s$^{-1}$, while the linewidths are 48.6$\pm$`<!-- -->`{=html}0.5 and 26.4$\pm$`<!-- -->`{=html}0.8 km s$^{-1}$, respectively. Later studies of several recombination lines at 9 cm toward S284 encountered similar values for this first velocity component (V=45.2$\pm$`<!-- -->`{=html}1.5 km s$^{-1}$, $\Delta$V=28.2$\pm$`<!-- -->`{=html}3.6 km s$^{-1}$). The narrow linewidths of the ionized gas tracers toward S284 are compatible with a diffuse H  region origin. These estimates of the radial velocities agree with the values obtained from the observation of the CO emission . Radial velocities in the range V$_{\rm{LSR}}$ $\sim$`<!-- -->`{=html}42.0–46.8 km s$^{-1}$ have meen measured at various locations in the field, yielding heliocentric distances between 5.2 and 6.5 kpc.

S284 hosts a cluster of stars in its centre known as Cl Dolidze 25 ; aka C 0642+0.03 or OCL-537. The stellar population of this cluster was initially characterised through UBV–H$\alpha$ photometry by . They concluded that Dolidze 25 was a cluster of 10 O and B type stars and 1 F supergiant. However, in their study of a 23 $\times$ 16 field, the authors separated the cluster into two subclusters with different reddening. A northern group of foreground stars of E($B-V$)$\approx$`<!-- -->`{=html}15 at a distance of 0.75 kpc and a redder and more concentrated central cluster. have pinned down its reddening to E($B-V$)$\approx$`<!-- -->`{=html}0.83 and the distance to 5.0 kpc, allowing for a metallicity gradient in the Galactic rotation curve. These results cast doubt on the physical connection between the different clusters found in this field. Although much has been discussed about the actual members of the cluster , most authors agree on the distance to the central concentration of stars. recalculated this parameter using previous UBV photometry and spectroscopy yielding a distance of 6.03$\pm$`<!-- -->`{=html}1.16 kpc. However, the selected stars to derive this value do not belong to the central subcluster. A more recent spectroscopic study of the sources Dol 25–15, Dol 25–17, and Dol 25–22 within this central subcluster the authors determined: (i) earlier spectral types than prevously derived for these sources, (ii) the low metallicity of the cluster (Z$=$`<!-- -->`{=html}0.17 Z$_{\sun}$), (iii) a spectro-photometric distance of 5.5$\pm$`<!-- -->`{=html}0.5 kpc. We adopt a distance of 5 kpc throughout the present study. The age of the cluster (7 Myr) has been estimated through isochrone fitting of precision photometry in the 4 $\times$ 4 field around Dolidze 25 .

Water-maser emission has been detected toward a few locations around S284 , whereas no SiO, CS, or methanol maser emission (i.e. tracers of very dense ambient material) have been detected .

The objective of this paper is to determine the properties of the different regions identified toward this site and their physical interconnection. In particular, it focuses on scrutinizing the evidence for sequential star formation.

This article is structured as follows. In Sect.  we describe the observations of this star-forming region (SFR), obtained with Spitzer/IRAC (Sect. ). We also present ground-based H$\alpha$ observations obtained with INT/WFC (Sect. ), as well as additional optical data, which we cross-correlate with 2MASS and Spitzer/IRAC. Secondly, we focus on the stars and distribution of young stellar objects (Sect. ), followed by an analogous study of blobs of extended emission (Sect. ). In Sect.  we discuss the spectral types of the stellar content in S284 (Sect. ) and the relative ages and timescales we determine (Sect. ). In Sect.  we discuss the uniqueness of this nebula in terms of its particular environmental conditions. Finally, conclusions are drawn in Sect. .

# Observations and data reduction

## Spitzer/IRAC Imaging

We observed S284 with the IRAC (InfraRed Array Camera) onboard Spitzer during Cycle 4 (2005 March 28) under Programme ID 3340. The observations were conducted in the four IRAC filters at 3.6, 4.5, 5.8, and 8.0 $\mu$m . S284 spreads over an area of 1 $\times$ 1 with its centre at $\alpha_{2000}$$=$ 06$^h$ 45$^m$ 00,  $\delta_{2000}$$=$ $+$`<!-- -->`{=html}00 17 56 . We performed a raster map with IRAC mapping mode in each filter with a spacing between rows of $1/3^{\rm rd}$ of the camera field of view (5 $\times$ 5). The IRAC simultaneous imaging in two consecutive filters results in an offset for two of the mosaic fields of around 6 (in our case, in declination). A final area, covered in all filters, of approximately 0 $\times$ 1 was obtained. We obtained an average integration time of 80.4 seconds per sky pixel.

The Spitzer Science Center provided the Spitzer/IRAC basic calibrated data (BCD), using the pipeline version S11.4.0. Further treatment of the images with the contributed IRAC software was necessary to correct for the muxbleed and pulldown effects of the detector. Likewise, the difference in the background levels of overlapping frames present in the IRAC 8.0 $\mu$m channel was corrected with the post-BCD overlap correction software. Subsequently, we used the dedicated software package MOPEX (version 030106) to produce mosaics of the IRAC images and to perform PSF photometry .

MOPEX performs image interpolation and co-addition of the individual images onto a common grid. Moreover, it allows outlier rejection through four different algorithms. Considering the shallow coverage of the field, cosmic rays and other outliers were removed using the box outlier and dual outlier MOPEX modules. The mosaic final pixel scale corresponds to the original 1 /pix. An image of this field in three IRAC colours is shown in Fig. .

The source extraction is preceded by the construction of our own point response function (PRF i.e., an oversampled PSF). We selected 30 isolated bright stars by hand from the final mosaic for each filter. The postage stamp images of these sources were resampled to 1/4th. We normalized the PRF for a radius that encompasses the first Airy ring for each band, therefore an aperture correction needed to be applied for the extracted flux on each single channel. The presence of strong filamentary extended emission in the four bands required the use of a small-width median filter for a proper background subtraction (11for IRAC 3.6 $\mu$m and IRAC 4.5 $\mu$m and 6 for IRAC 5.8 $\mu$m and IRAC 8.0 $\mu$m). This aggressive filtering was accordingly applied to the channel PRF, producing an underestimate of the source flux. The correction factor was determined from the flux ratio between the unfiltered and filtered PRFs .

The source detection was performed in the single final coadded frame for each individual channel. To ensure a high reliability, we only selected sources with S/N $>$ 5. We finally obtained extraction lists of 124 657, 118 289, 19 445 and 15 856 sources in the IRAC 3.6, 4.5, 5.8, and 8.0 $\mu$m channels, respectively (hereafter referred to as IRAC 3.6, IRAC 4.5, IRAC 5.8, and IRAC 8.0). The final merged list of the four bands was obtained by matching the coordinates of the extracted sources in the four bands (using an optimally determined maximum distance of 1). To minimize the contamination by spurious detections in the source catalog we apply a filtering criterion based on simultaneous detection in multiple filters. We require detection of the object in either the two bluest (IRAC 3.6 and 4.5) or reddest bands (IRAC 5.8 and 8.0). This criterion accomodates faint blue sources and very red sources allowing them to survive into the source list. A total of 8 010 sources were extracted in the four filters. Since the IRAC/Spitzer observations were aimed at relatively faint sources, very bright objects appear saturated in our images. We determine the saturation limit of each band by comparison of the fluxes obtained with PRF fitting and aperture photometry (aperture radius of 3 pixels). The limit is set by eye at the beginning of the nonlinear regime between these two fluxes.

## Cross-correlation with 2MASS

In addition to the Spitzer/IRAC we used the 2MASS catalog. We cross-correlated the point sources found in the IRAC images with this catalog to obtain $J$-, $H$-, and $K_s$-band photometry. At least one of these colours was found in 2MASS for 23 734 point sources at a distance closer than 1.

## Cross-correlation with INT/WFC

Since the field containing S284 is a target field of the CoRoT satellite, it has been studied by CoRoT's exoplanet team led by M. Deleuil . As a preparation for the CoRoT mission, this team obtained 6187 images in CoRoT's fields with the WFC (Wide Field Camera) at the 2.5m INT (Isaac Newton telescope) in La Palma (Spain), in order to classify the stars and select the best candidates for planetary transits. They used the $B$ and $V$ Bessel filters and the *r'* and *i'* Sloan Gunn filters. See Deeg et al., in preparation for more details. The exoplanet team of CoRoT kindly put at our disposal the photometry they derived in these four filters from their observations for the field of S284. The mean seeing of the data for this field is about 1.5and the mean astrometric error is 0.06.

In addition, we have obtained specific observations of S284 in 2005 December 18 with the same instrumentation (WFC/INT) with a focus on the extended emission. We used the H$\alpha$ and ($\lambda$ 6568) Sloan Gunn *r'* filters to trace the recombination line and continuum, respectively. Two pointing positions were needed to cover the nearly 1 $\times$ 1 region. We obtained two exposures of 1200 s in H$\alpha$ and 300 s in *r'*. After bias, flatfield, and astrometry corrections, the 4 different CCDs were independently corrected for optical distortions with a $3^{\rm rd}$ order polynomial. The different exposures and pointings were combined with the IRAF routine MSCRED. Subsequently, the stellar continuum was subtracteed in order to obtain a line emission image. After exposure time synchronization, the *r'*-band mosaic was multiplied by the flux ratio of 23 stars in the field in both bands. Since we are interested in the extended emission, no absolute flux calibration was performed on the data. The H$\alpha$ (continuum subtracted) image is shown in the right panel of Fig. .

# Results

Large dust bubbles are ubiquitous in the Galactic plane toward the Galactic centre as shown by the GLIMPSE Survey . However, the study of these structures is being hampered by the difficulties imposed by crowding effects, the strong background toward the Galactic midplane, and the frequent asymmetries in the bubbles as a result of interactions with closeby regions.

The unprecedented view of S284 in the mid-IR provided by IRAC/Spitzer (see Fig.  *left panel*) has revealed in detail the presence of complex extended emission in this region, which was not previously detected in 2MASS observations. The emission is present in all four IRAC bands implying diffuse dust emission. This interpretation is corroborated by its detection (although at much poorer resolution) in IRAS maps, even at 100 $\mu$m. The mid-IR emission appears the strongest in the IRAC 8.0 band, since this band contains the PAH band emission at 7.7 $\mu$m and 8.6 $\mu$m. On the largest scale toward S284, the images reveal that the dust shells bound the bright H$\alpha$ emission quite well (see also Fig.  *right panel*), therefore tracing the photon-dominated regions (PDRs) and the dense swept-up shells. These two complementary maps allow us to identify three bubble-like structures shown in Fig. . In fact, the IR emission extends well beyond the 18-pc inner bubble radius of the largest bubble (hereafter B1). The dust emission almost fully surrounds the ionized gas although it appears less prominent to the southwest, likely due to a density gradient within the molecular cloud. This is also evident on the much larger scale traced by the Galactic dust extinction map (see Fig. ), where lower extinction values are detected in this direction. In the same direction, the H$\alpha$ map shows a faint extended emission component. This leak of ionizing photons suggests that the H  region may be density-bound to the southwest. The OB cluster Dolidze 25, as outlined by and , is located at the geometrical centre of B1. Interestingly, our H$\alpha$ image (Fig.  *right panel*) traces the bright rim of dense ionized gas, especially bright to the southeast of the ionized bubble.

The resolved filamentary structure of the dust allows differentiating among several other regions in the field of S284. Another two shells of dust enclose more compact H  regions within the field around S284. The bubble B2 is located to the northwest, sitting amidst the ionized emission of B1 and yet very close to its rim. This compact bubble (diameter$\sim$ 3.5 pc), seems associated with IRAS06422+0053 which is the brightest point-source at 100 $\mu$m. The low extinction observed toward this area of ionized gas, together with the higher brightness in H$\alpha$ indicate that it is probably located in the foreground of B1. We detect a stellar cluster in its centre labelled as Cl2 (see the left panel in Fig. ). Two of the sources of Cl2 belong to the trail of stars that has often been associated with Dolidze 25.

A third bubble of 4-pc diameter lies within the mid-IR emission and is located to the southeast of B1. The ionized gas is barely detected in the INT H$\alpha$ mosaic. This source coincides with IRAS06439-0000 and IRAS06439+0111. We report the detection of another stellar cluster (labelled as Cl3) at the centre of B3. Overall, these three bubbles show a remarkable symmetry despite the highly probable interaction among them.

Another prominent type of structure observed toward S284 are several cometary globules that are detected at the rim of bubble B1. These globules appear in emission in the IRAC images and as absorption silhouettes in the H$\alpha$ image. Moreover, they appear highlighted in the latter map by the presence of ionization rims at the head of the globules. All of them point toward the centre of B1 showing very diverse stellar content (see next section) and morphology with their location around the nebula. Globules to the northeast show wide tails and harbour resolved clusters of sources. They also seem to be located at various radial distances from Dolidze 25, possibly due to a projection effect (e.g. the region labelled as RN in Fig. ). The cometary globules or elephant trunks to the south and west possess longer and thinner tails and only contain a few unresolved sources at the head, in IRAC 3.6 and 4.5.

As measured by the shielded material behind the condensations, this distinctive morphology between cometary globules and elephant trunks supports our hypothesis of a density gradient in the ambient material. However, observations of the molecular gas are needed to confirm this statement.

Bubble B3 shows only one cometary globule at the eastern rim of the dust bubble. A resolved source appears in the head of this cometary structure, which points to cluster Cl3. Some of these globules and trunks are associated with isolated IRAS sources, indicating an enshrouded nature. Discussion on their stellar content follows in the next section.

## Stellar content: point-like sources

Our IRAC/Spitzer data has unveiled the presence of a rich stellar population in S284 in the form of various peripheral clusters. It also allows us to put constraints on the debated extent of the cluster members of Dolidze 25 through the link between cluster and larger dust structures. Another important point is the inferred presence of high-mass stars in clusters Cl2 and Cl3 through the detection of associated ionized regions, which requires the existence of at least one star of spectral type earlier than a B3V .

Our mid-IR IRAC data allows the classification of point-like sources based on the differentiation between sources that show an IR excess and those that are photosphere-dominated. Furthermore, it is also used as a young stellar object subclassification between Class I and Class II objects by comparison with models of proto-stellar envelopes and circumstellar disks . The effective temperature and luminosities of young stars can be used to infer their ages and thus the timescale for the different evolutionary stages of the circumstellar material. Although the literature comprises information on the spectral types of a small sub-sample of bright stars , the reported spectral types vary significantly; for instance, the spectral type clasiffication of the brightest source in Cl2 –an A0 star according to – is not sufficient to produce the observed H  region. As an alternative, a crude estimate of the spectral type can be derived from the stellar brightness in the IRAC 4.5.

### IRAC/Spitzer colour-colour diagram

Our IRAC observations open a new window on the stellar content in S284. This improvement does not come so much from peering through the extinction at longer wavelengths, but from the fact that the sensitivities achieved in IRAC 3.6 and IRAC 4.5 are higher than existing near-IR surveys. Our mid-IR colour-colour stellar classification is restricted by the limiting magnitudes in IRAC 5.8 and IRAC 8.0 (14 and 13.5 mag, respectively), since simultaneous detection in the four bands is needed for this classification method. In addition, to increase the reliability of our classification, we only consider those objects with photometric error smaller than 0.175 mag in each of the four bands.

The diagram using a total of 4046 sources is presented in Fig. . Bona fide candidate Class I objects are selected from the region marked in the diagram with the restriction that \[3.6\]$-$\[4.5\]$>$`<!-- -->`{=html}0.4 mag. The bulk of main sequence and giant stars is centred at a \[3.6\]$-$\[4.5\]$=$ $+$`<!-- -->`{=html}0.1 mag, thus on average slightly bluer than expected. As a consequence, main sequence and giant stars contaminate the area of the diagram where Class II objects lie according to the models . For this reason, we have designated as Class II sources those whose \[3.6\]$-$\[4.5\] colour is above the visible gap in this direction, namely \[3.6\]$-$\[4.5\]$>$ 0.3 mag. A particular group of sources lies in a region of the diagram where \[5.8\]$-$\[8.0\]$>$ 1.4 mag and \[3.6\]$-$\[4.5\]$<$ 0 mag. These points correspond to enhancements in the filamentary extended structure that MOPEX detected as point-like sources whose colours correspond to reflected emission, i.e., reflection nebulae.

The full lists of Class I and Class II candidates are available electronically in Tables 3 and 4.

### 2MASS diagrams

The complementary information provided by the 2MASS catalog is especially useful for the brightest stars that appear saturated in our IRAC images. This is the case for bright central sources located in the heads of a few cometary globules. It is important, however, to bear in mind the slightly lower spatial resolution of the 2MASS survey, especially toward the centre of the more dense clusters. The near-IR colour-magnitude and colour-colour diagrams of the sources within S284 are displayed in Fig. .

For the colour-magnitude diagram (CMD), we consider sources with relative photometric error below 20% in the three 2MASS bands. The main sequence (MS) at a distance of 5 kpc with A$_{V}=0$ and A$_{V}=5$ mag are displayed for reference. Note that the estimated average extinction toward Dolidze 25 corresponds to A$_{V}$ $=$ 1.98 mag . The M$_\mathrm{V}$ magnitudes of luminosity class V are obtained from . For the colours of the spectral types earlier than A0V we use the values from and for later spectral types we obtained the colours from . The near-IR extinction law adopted is from . We select the sources redder than the gap noticeable at $J$ $-$ $K_s$ $\sim$ $-$`<!-- -->`{=html}0.3 mag in our CMD and brighter than an accordingly reddened star B8V (points represented in red in Fig. ). These objects are either intermediate- and high-mass stars that are strongly reddened or more evolved stars that have reached the giant phase. These presumably reddened stars are plotted in the colour-colour diagram (CCD) shown in Fig.  (*right panel*), where many of them lie in the red giant locus.

In our CCD we consider that stars posses a near-IR excess when they appear redward of the region limited by the reddening vectors projected from the intrinsic colours of main sequence stars and giants and above the T Tauri locus . These objects (marked with triangles in the diagram) are considered to be YSO candidates . The same triangle-like symbols are used to highlight those objects back on the CMD.

### The distribution of YSOs

The spatial distribution of YSO candidates in S284 is shown in Fig. . The symbol sizes for the Class II candidates are scaled with the IRAC 4.5 brightness as a proxy for the luminosity of the YSO . This choice is motivated by the fact that the PRF is undersampled in the IRAC 3.6 band, while oversampled in the other bands; IRAC 4.5 thus offers a more reliable measure of the flux. Moreover, the extended emission appears least prominent in the 4.5 $\mu$m map. Analogously, the symbol sizes for Class I objects are proportional to the IRAC 5.8 brightness. Although IRAC 8.0 seems better suited to YSOs with a strong excess, the 7.7 $\mu$m PAH emission band could contribute significantly to the source flux estimate .

Most Class I and Class II candidates arise in clusters around the shell of the B1 dust bubble or at the locations of the powering clusters Dolidze 25, Cl2, and Cl3. Other YSOs around B1 seem to be more scattered, but for some cases their low brigthness make them dubious. Several bright Class I objects are located at the head of cometary globules or at elephant trunks around B1. In particular, the longest trunk of the region harbours several Class I candidates at various positions along the trunk (RW in Fig. ). Those candidates at the rim of the B1 bubble show a colour \[3.6\]$-$\[4.5\]$>$ 1.1 mag indicative of either strong reddening or high stellar temperatures . We also detect two clusters of Class I and Class II objects (RE in Fig.  and another concentration 15$\arcmin$ north of RE) a few parsec distant from the ionization front associated with B1. To the north of the IRAC field, we find one isolated Class I candidate at the centre of a region of strong radio continuum emission.

Surprisingly, Dolidze 25 comprises several Class I objects located in a 4$\arcmin$ projected diameter region. This result will be discussed further in Sect. . Cl2 and Cl3 seem to have a small population of YSOs dominated by Class II objects, whereas a Class I candidate is detected at the rim of each bubble. In conclusion, the location of the YSO candidates follows closely the morphology of the nebulae. There is a higher concentration of young objects along the IR bright rim of B1. However, a significant number of objects are detected as far as few tens of parsecs away from the rim of B1.

### Comments on the stellar population of individual regions

#### Dolidze 25

Most of the objects that have been traditionally considered as members of Dolidze 25 (e.g. \#17, \#15, \#16, \#35, \#22, \#18, \#19) do not show any circumstellar peculiarity at near- and mid-IR wavelengths (see Fig. ). A concentration of Class I candidates exists toward Dolidze 25, while low-luminosity Class II sources appear more scattered around the cluster. Most of the Class I candidates are also detected in the $J$ band ruling out the possibility of being reddened Class II objects.

#### Cl2

The area around Cl2, also including B2, shows a high concentration of Class II objects. The source located at the centre of the cluster (source \#9 in Table ) has been spectroscopically classified in the optical as an O9V star . The cluster core appears elongated at mid-IR wavelengths and starts to be resolved redward of 5.8 $\mu$m. The counterpart to O9V star exhibits the mid-IR colours of a Class II object, although this could be due to the contamination by the nearby object. Another source observed in the direction of the H  region (#8) appears extended in our IRAC 3.6 image and is barely resolved into two point-like sources in the redder IRAC maps. A bright Class I object is detected at the rim of the dust bubble, where the extended emission is the strongest. Another nearby object shows a strong near-IR excess (NYC1).

#### Cl3

The two brightest objects within the central cluster (Cl3\*1 and Cl3\*2) are detected in all four IRAC bands, but their SEDs only show a photospheric component. Their location in the near-IR CMD suggests that these are O-type stars. In particular, Cl3\*1 is closely surrounded by sources that show a mid-IR excess. NYC-2 is located at the head of the cometary globule east of bubble B3 and appears saturated in IRAC 3.6, IRAC 4.5, and IRAC 8.0, so it is not included in our analysis of Sect. . However, the overall shape of its SED corresponds to a Class I candidate.

#### RN

The field labelled as RN includes several cometary globules to the north of B1. Two of the brightest sources at the vortices of cometary globules appear saturated in one of our IRAC maps, and yet, bright Class I objects are unveiled within these clusters, pointing to the presence of various YSOs embedded in these regions. The detection of two Class II candidates (RN-13 and RN-3) is reinforced by the near-IR excess shown in the 2MASS photometry (sources NYC3 and NYC5). One of them coincides with the head of a non-saturated globule that hosts only few point-like objects.

#### RE

The field RE has two blobs of strong emission separated by a projected distance of $\sim$`<!-- -->`{=html}0.4 pc. The eastern blob is resolved into extended emission, while the western blob has several resolved point-like sources, the central one saturated. The 2MASS photometry toward this source reveals a near-IR excess at a slightly worse resolution than our IRAC images. A resolved cluster is detected to the south of the blobs. Seven bright Class I candidates are aligned in a trail, while several Class II sources are scattered around. This cluster is partially detected in the near-IR and optical, indicating that it is not deeply embedded within the dust; hence, we can rule out that the Class I candidates are reddened Class II objects. Moreover, these Class I candidates are the only objects not located at the rim of B1 that show a colour \[3.6\]$-$\[4.5\]$>$ 1.1 mag, pointing to relatively high luminosities , but not hot enough to generate an H  region.

#### RW and RS

The region RW encompasses the longest trunk at the rim of B1. Several Class I objects are located along the pillar. One of them is also detected by an associated near-IR excess.

Analogously, the RS region contains another two trunks that harbour Class I sources at the head of their cometary globules.

## Stellar content: embedded objects and blobs

The previous sections have focused on the stellar classification of point-like sources in S284. However, the filamentary high-density structures detected in our IRAC maps support the possible presence of embedded clusters in regions where the emission appears resolved.

We selected 11 blobs or regions of small extent (a sub-scale of bubbles B2 and B3) and high infrared brightness in order to search for young embedded clusters. The radius of these blobs was defined as 1.25 $\times$ HWHM in the IRAC 8.0 band and kept constant to calculate the integrated flux in other bands. The same radius was used to search for near-infrared counterparts in the 2MASS point-source database. When multiple sources were present, we summed their contribution to estimate the total near-infrared flux. For the objects that saturated an IRAC-band image, we considered the IRAC fluxes as a lower limit and MSX fluxes when available. We also searched for far-infrared counterparts in the IRAS point sources catalog and selected these that are compatible with the blobs as their main contibutors in the IRAS beam on the basis of HIRES images. As expected, the more extended blobs are more likely to have an IRAS counterpart.

All these data points were used to build the spectral energy distributions of the blobs presented in Fig. . We fitted the 2MASS, IRAC, MSX and IRAS flux distribution with a combination of blackbodies with temperature ranging from 30 to 1000 K to reproduce the shape of the SEDs and estimate their corresponding luminosities. Figure  displays the spectral energy distribution as well as the results of the fit. The blob characteristics are described in Table , while their location and radius are represented in Fig. .

We identify four different types of SEDs in column (6) of Table : (I) ClassI objects unresolved in IRAC 8.0 and with integrated near-infrared luminosities $<$ 2.0; (II) ClassII objects, also all unresolved in IRAC 8.0; (R) objects having reflection colours; and (Cl) objects falling in the Class I region in the IRAC colour-colour diagram with higher near-infrared luminosities. This last category also tends to have a steeper 3.6 to 8 $\mu$m slope than unresolved ClassI objects, nearly as steep as reflection objects for some of them. Still, their mid-infrared spectrum cannot be fitted with a single temperature blackbody, but on the contrary require an extended range of temperatures. Their luminosity, extent, and spectral distribution suggest that these are young clusters, either constituting an earlier cluster phase compared to bubbles B2 and B3 or lower-mass clusters.

# Discussion

## Environmental conditions

Before we consider the stellar content and the ages of the various distinct regions, it is instructive to debate the global morphology and what can be concluded from it. The shell thickness of dynamically blown bubbles is a function of the central star, the age of the bubble, and the density of the ambient ISM . The theory predicts that the ratio of shell thickness to radius increases with age and ambient interstellar density. Therefore, for relatively symmetric bubbles with the powering star(s) located at the geometrical centre (as is the case of S284), a smaller ratio is indicative of a lower ambient density. We observe that the ratio of shell thickness to radius is not constant around B1, decreasing along the southwest direction. A lower ambient density in this direction would explain this effect and support our interpretation of the lower dust extinction values inferred from the IRAS maps as caused by a local ambient density gradient in this direction.

The presence of elephant trunks suggests that the dynamics of the H  region within B1 are dominated by the thermal pressure of the gas at 10$^4$ K and not by shocked gas, either produced by stellar winds or supernovae. These sawtooth shapes have been successfully modelled as the result of dynamic instabilities in the ionization front of an expanding H  region, regardless of pre-existing density inhomogeneities, by . These models account for different cooling and background conditions and outline the effect of low-metallicities (final shell temperatures are higher) into thicker shells and wider and sparser trunks.

Some characteristics of this SFR may be determined by the combined effect of the lower than solar metallicity and a gas pressure that is lower than typically found in the inner MilkyWay. High-mass stars exhibit weaker winds at lower metallicity due to the less efficient momentum transfer of the stellar radiation on the outer envelope. In a study of young SFRs in the SMC, find that the morphology of the region is not shaped by wind interaction but more by the expansion of the ionized bubble, much like what is found for S284. A second effect concerns the dust shielding. Low-metallicity environments like the LMC and SMC show a deficit of the smallest dust grains , which may cause the UV radiation of the newly formed stars to permeate the surrounding molecular cloud and modify the conditions over great distances (e.g. dust composition, charge state of the smallest grains, photodissociation of molecular gas, Hony et al. in preparation). It is interesting to compare this with the results of . These authors study the star formation in a sample of late-type galaxies using IR and UV and find that the star formation in the outskirts of these galaxies is typified by less attenuation due to the lower metallicity and that these regions are underluminous in H$\alpha$. One possible explanation for the lack of prominent H$\alpha$ emission is that the lower ambient pressure makes the average star-forming cluster less massive, and as a result less very massive stars are formed. S284 might be an example of such a less massive star-formation cluster, with an O7V the most massive star detected.

## Spectral types and masses

The lack of detailed information on the spectral types of the complete stellar content precludes a thorough age analysis of S284. The use of spectral types has proven useful for determining age gradients in massive clusters like the Eagle nebula and probing the disk evolution toward several Galactic SF regions . However, we qualitatively discuss this issue with the information presently at hand.

Table  compiles the spectral type and luminosity class classification and membership for S284 from the literature. The various maps tracing the dust and ionized gas complement this existing information on S284 distinguishing among the different physical structures (clusters, H  regions and bubbles).

## Relative ages and timescales

```latex
\begin{table*}[!t]\caption{Characteristics of the selected embedded regions and
      comparison with Class I and II unresolved objects.}
    %
    \small
    \begin{tabular}{lccccccc}
      \hline\hline
      Blob ID & RA & DEC & Radius&Log(L/Lo)&Class&Region&IRAS assoc.\\
      &    &     &  ('') &         &     &      &           \\
      \hline
      18..............    & 06  45  40.8 & +00  22  27.1 &  5.73 &2.5 &Cl  &  &06430+0025\\
      22..............    & 06  46  29.6 & -00  03  29.2 &  6.06 &2.6 &Cl  &Cl3&06439-0000\\
      21..............    & 06  46  15.8 & +00  06  26.2 & 15.92 &3.0 &Cl  &RE&06437+0009\\
      12..............    & 06  45  16.0 & +00  22  26.9 &  8.24 &3.9 &Cl  &RN&06426+0025\\
      16..............    & 06  45  31.8 & +00  06  26.1 &  5.31 &2.0 &Cl  &RS&06430+0009\\
      11..............    & 06  45  14.1 & +00  25  18.6 & 10.86 &2.6 &Cl  &RN&06426+0028\\
      25..............    & 06  46  23.3 & +00  05  48.6 & 22.39 &2.9 &R   &RE&\\
      28..............    & 06  44  49.1 & +00  20  15.3 & 20.14 &2.7 &R   &Cl2&06422+0023\\
      26..............    & 06  44  31.8 & +00  23  23.6 & 13.22 &2.3 &R   &  &06419+0026\\
      27..............    & 06  45  18.9 & +00  03  49.2 & 12.80 &2.3 &R   &RS&\\
      14..............    & 06  45  17.5 & +00  03  30.3 &  5.04 &1.8 &R   &RS&\\
      15..............    & 06  45  22.0 & +00  04  49.4 &  U    &1.7 &I   &RS&\\
      19..............    & 06  45  58.8 & +00  13  52.9 &  U    &0.9 &I   &  &\\
      5................    & 06  44  36.5 & +00  13  51.2 &  U    &1.5 &I   &RW&\\
      13..............    & 06  45  17.3 & +00  23  45.6 &  U    &1.8 &II  &RN&\\
      17..............    & 06  45  37.4 & +00  19  49.8 &  U    &1.3 &II  &  &06431+0023\\
      4................    & 06  44  32.9 & +00  10  52.4 &  U    &1.2 &II  &RW&\\
      \hline
      \label{blob_table}
    \end{tabular}
\end{table*}
```

The different spatial scales of H  regions and their distribution around S284 is suggestive of different epochs of star formation. Nevertheless, a sequence of events and causality between them must be demonstrated to confirm a triggered star-formation scenario.

Low-mass Class I objects have an estimated age of a few 10$^5$ yr, while Class II objects range between 10$^6$ and 10$^7$ yr . We can therefore use these two classes of protostars present in the S284 field as a relative age scale to study the interaction of young clusters ($<$ 5 Myr) with their environment.

The three dust bubbles that bound the H  regions detected on our H$\alpha$ map (B1, B2, and B3) share a common morphology of cometary globules or elephant trunks. In the case of B1 and B3, they are clearly protruding into the ionized gas. The brighter young YSOs (Class I) and proto-clusters appear located at the vortices of these structures, pointing to a later formation than that of the central clusters in the H  regions. That no compact radio continuum emission, maser emission or strong CS have been detected at these locations supports the idea that these young protostars and proto-clusters are probably not massive. The morphology of this younger generation of protostars and proto-clusters is compatible with a process of dynamical instabilities of the ionization front . Because the condensations do not seem to share the expansion of the ionized gas, another plausible explanation would be the radiation driven implosion of pre-existing molecular condensation. The morphology of the region at the rim of B2 where several YSO candidates are found resembles the cases where the collect and collapse process has be reported to be at work . The velocity distributions of the ionized and molecular gas have been proven valuable when studying the association of condensations and YSOs with H  regions .

B1 shows a thick dust shell with a combination of Class I and older Class II candidates, both in clusters and in isolation, farther away from the ionization front. report a similar result in their study of the RCW 120 region, explaining it in terms of a leak of radiation into the PDR.

A question arises after discussion of the impact of the various clusters on their immediate surroundings: what is the interplay between the high-mass stellar populations in this region? are Cl3, and especially Cl2, produced by the expansion of the H  region powered by Dolidze 25?

The powering clusters Cl2 and Cl3 are dominated by Class II candidates. It is therefore puzzling to detect younger Class I candidates in the direction of the central cluster Dolidze 25, the presumed oldest cluster in this region. In the case that these YSO candidates were associated to the cluster, these results would challenge the triggered SF scenario initially suggested by the multiple bubble morphology of S284, in particular the close location of Cl2 to the B1 shell. The UV field produced by the several O-type stars of Dolidze 25 must contribute to the photo-evaporation of circumstellar material in nearby stars. However, YSO candidates have also been detected close to the young massive cluster NGC 6611 in the Eagle nebula . The authors conclude that circumstellar disks must be fairly robust to the cluster environment around NGC 6611. Another possible explanation would be that these Class I objects are located at the border of the bubble, either in the foreground or in the background. Further support for this hypothesis could come from the detection of several Class II candidates in Fig. . These Class II sources would trace the stellar population associated with a large emission filament that crosses B1 in the northeast direction in our IRAC images and appears in absorption in our H$\alpha$ image. The proximity *in projection* of this filament suggests that this or other structures in the line of sight could overlap with the cluster Dolidze 25. If the Class I objects detected near Dolidze 25 were not members, the cluster would be dominated by a stellar population of naked photospheres. A future radial velocity study of the stars around Dolidze 25 would aid in answering this question.

### Dynamical ages of the H  regions

|  |  |  |  |  |  |  |  |
|:---|:--:|:--:|:--:|:--:|:--:|:--:|:--:|
| ID | Classification |  |  |  |  |  |  |
| \(1\) | \(2\) | \(3\) | \(4\) | \(5\) | \(6\) | \(7\) | \(8\) |
| 1 | \(m\) | B0IV | B5(m?) |  | B5 |  | \(n\) |
| 2 | \(n\) |  | A5(–) |  | A5 |  | \(n\) |
| 3 | \(n\) |  | A0(–) |  | A0 |  | \(n\) |
| 4 | \(n\) |  | B2(–) |  | F5$^{\dag}$ |  | \(n\) |
| 5 | \(n\) |  | B5(–) |  | A0$^{\dag}$ |  | \(n\) |
| 6 | \(n\) |  | A5(–) |  | A5 |  | \(n\) |
| 7 | \(n\) |  | A0(–) |  | A0$^{\dag}$ |  | \(n\) |
| 8 | \(n\) |  | A5(–) |  | A2$^{\dag}$ |  | (Cl 2 m) |
| 9 | \(m\) | O9 | A0(m) |  | A0 |  | (Cl 2 m) |
| 10 | \(m\) | F2Iab | F0(m) |  | F5$^{\dag}$ |  | \(n\) |
| 11 | \(n\) |  | K2 |  | K2$^{\dag}$ |  | \(n\) |
| 12 | \(m\) | B0III | B2(m) |  | B2 |  | (Dol 25 m?) |
| 13 | \(m\) |  | B8(m) |  | B8 |  | (Dol 25 m?) |
| 14 | \(n\) |  | B8(–) |  | B8 |  | (Dol 25 m) |
| 15 | \(m\) | B0III | B8(m) | \(m\) | B8 | O6$^{\ddag}$ | (Dol 25 m) |
| 16a | \(x\) |  |  |  |  |  | (Dol 25 m) |
| 16b | \(m\) |  |  | \(m\) | B8 |  | (Dol 25 m) |
| 17 | \(m\) | B1V | B8(m) | \(m\) | B8 | O5.5–O6.5$^{\ddag}$ | (Dol 25 m) |
| 18 | \(m\) |  | F5(m) | \(m\) | F5 |  | (Dol 25 m) |
| 19 | \(m\) |  | A5(m) |  | A5 |  | (Dol 25 m) |
| 20 | \(n\) |  | F5(–) |  | F5 |  | (Dol 25 m) |
| 21 | \(n\) |  | F5(–) |  | F5 |  | (Dol 25 m) |
| 22 | \(m\) |  | A5(m) | \(m\) | A5 | O8–O9$^{\ddag}$ | (Dol 25 m) |
| 23 | \(n\) |  |  |  |  |  | \(n\) |
| 24 |  |  | B8(m?) |  |  |  | \(n\) |
| 25 |  |  |  |  |  |  | \(n\) |
| 26 |  |  | K0(–) |  |  |  | \(n\) |
| 27 |  |  | F5(m) |  |  |  | \(n\) |
| 28 |  |  | F5(m) |  |  |  | \(n\) |
| 31 |  |  | F0(m) |  |  |  | \(n\) |
| 32 |  |  | F0(m) |  |  |  | \(n\) |
| 33 |  |  | B8(m) |  |  |  | \(n\) |
| 34 |  |  | F5(m) |  |  |  | \(n\) |

Membership and Spectral Types of Dolidze 25. Column (1) indicates the star number following the notation of . (m) members, (n) non-members.

(2) , (3) , (4) , (5) , (6) WEBDA Open Cluster Database http://obswww.unige.ch/webda, (7) (8) This work, $^{\dag}$ Spectral type classification from , $^{\ddag}$ Spectral type determination following

The exceptional symmetry of S284 allows a rough calculation of the dynamical age of the H  region using a zero-order approximation (i.e. a simple pressure driven expansion), using the strong shock approximation and assuming a constant-density medium. We convert the spectral types of the most certain members from the literature (see Table ) into the corresponding Lyman continuum photon rate. Dolidze 25 and its associated H  region are considered separately from Cl2 and its compact ionized region. Adopting an isothermal sound speed $\sim$`<!-- -->`{=html}0.4 km s$^{-1}$, a recombination coefficient $\sim$`<!-- -->`{=html}3 $\times$`<!-- -->`{=html}10$^{-13}$ cm$^3$ s$^{-1}$ and the estimated electronic density n$_e$ $\sim$ 100 cm$^{-3}$ , we derive a dynamical age $<$ 3.8 Myr for B1. In the case of the compact H  region bound by B2, we assume n$_e$ $\sim$ 2$\times$`<!-- -->`{=html}10$^3$ cm$^{-3}$ and a sound speed $\sim$`<!-- -->`{=html}0.4 km s$^{-1}$, rendering a dynamical age $<$ 1.8 Myr.

### Comparison to isochrones and evolutionary tracks

We attempted a quantitative age analysis of the clusters Dolidze 25 and Cl2. For this purpose, we converted the photometry of the bright stellar population obtained with INT/WFC to the standard Johnson-Cousins photometric system (Deeg et al., private communication). We combined these data with the VIMOS photometry in the $R$ and $I$ bands of the faint population . The CMDs of both regions are presented in Fig. . The population around Dolidze 25 exhibits a main sequence, and even a hint of lower-mass pre-main sequence in the CMD. This is not the case for Cl2, where no distinctive features are detected in the CMD probably due to the poor statistics. The field around Dolidze 25 is not fully covered to the north of the cluster by the VIMOS observations. This could explain the void of sources detected at $I\sim$`<!-- -->`{=html}16 mag in the left panel of Fig. . The fit of a log(age\[yr\])$=$ 6.75 isochrone to the blue part of the main sequence renders an E($R-I$)$=$ 0.5 mag, considering a distance of 5 kpc.

The atmospheric parameter determination for source \#22 secures an age upper limit for the cluster of 6.2$\pm$`<!-- -->`{=html}0.6 Myr (assuming a coeval stellar population). This estimate was obtained via comparison to a 15 M$_{\sun}$-evolutionary track with Z=0.004 from , although source \#22 is located at the outskirts of Dolidze 25, while sources \#15 and \#17 belong to the high-stellar-density core of the cluster. Source \#15 is spectrally classified as a middle O-type dwarf star, while star \#17 appears as a more evolved hot star of spectral type III or Ib . Their positions in the CMD reveal the presence of differential extinction toward the centre of Dolidze 25, detected also in the $V$ vs. ($B-V$) diagram. Since the masses derived from atmospheric models are particularly uncertain for stars hotter than 35000 K, no comparison to an evolutionary track is attempted. Nevertheless, the presence of an O6V star (#15) confirms that Dolidze 25 is younger than 5 Myr.

The optimal isochrone fit to the stellar population of Cl2 in Fig.  indicates that it is also approximately 5 kpc away. We use the less certain spectral type classification of source \#9 as O9V and the comparison to a 15 M$_{\sun}$-evolutionary track to establish that the cluster must be younger than 1.0$\pm$`<!-- -->`{=html}1.0 Myr.

These results for Dolidze 25 and Cl2 are in good agreement with the dynamical ages estimated in the previous section from the correspondent H  regions. Moreover, in the case of Cl2, the detection of several Class II objects within the cluster is compatible with the estimated age.

The age estimates for the clusters Dolidze 25 and Cl2 are consistent with the inferred dynamical ages of their respective H  regions. This translates into an age spread of approximately 3 Myr between the two clusters. Considering the $\sim$`<!-- -->`{=html}20 pc separation between them, S284 lies on the age spread–linear size relationship described by for the LMC clusters. This correlation is explained as the result of the rapid formation of self gravitating sub-clumps of a larger cloud complex. In this context, these results hint that the formation of the massive clusters in S284 is not necessarily triggered.

# Conclusions

In this paper we have presented newly obtained IRAC/Spitzer and H$\alpha$/INT observations of S284. These observations show –at a very detailed level– a relatively undisturbed isolated site of high-mass star formation. We presented an overview of the young stellar content of this region and sketched a scenario for the evolution of this region, considering its morphology.

Our IRAC maps show the newly discovered dust diffuse component of the nebula, which allows the investigation of a multiple-bubble morphology. Since this sort of structure has recently been explained in terms of triggered star formation, we analyse the population and morphologycal features as proxies for their relative ages.

The dominant photospheric component traced by the filters IRAC 3.6 and IRAC 4.5 allows the detection of a cluster at the centre of the three H  regions shown by our H$\alpha$ map. These clusters must therefore harbour at least one high-mass star to blow such detectable ionized regions. Considering the remarkable symmetry of the bubbles and cluster centre positions, the morphologigal variations within the field are discussed in terms of the environmental conditions. Density variations along the field may play a major role in explaining the observed features.

Based on an IRAC colour-colour criterium on the photometry of point-like objects toward S284, we found a total of 155 Class I and 183 Class II candidates. Most of these young objects concentrate at the shells of the bubbles. The extrapolation of this analysis to the SEDs of integrated regions of extended emission yields the conclusion that six more embedded clusters are also located in the dust shell.

We have detected relative younger Class I candidates close to Dolidze 25 than the dominant Class II populations toward clusters Cl2 and Cl3. We considered a possible alignment of foreground or background, young YSOs in the line of sight of Dolidze 25, as well as a robust nature for low-mass disks, previously invoked for the case of NGC 6611 . We studied the dynamical ages inferred from the H  regions. These are in good agreement with the ages inferred from the older stellar population of Dolidze 25 and Cl2 through comparison of optical photometry with isochrones and evolution tracks. Based on the age-spread and linear-separation argument, we discard the expansion of the B1 bubble as triggering the formation of Cl2 and most probably of Cl3.

Reddened younger luminous objects (Class I) are only detected at the vortices of cometary globules or elephant trunks, while a mixed population of Class I and Class II objects is present both in a clustered and a more scattered mode farther away from the ionization front inside B1. This suggested age gradient is even more pronounced considering the dominant Class II populations within Cl2 and Cl3 and the presence of Class I candidates at their respective shells. As a result, we see indications of triggered star formation between the powering clusters of the H  regions and the physical structures at their respective shells.

Processes such as radiation-driven implosion of pre-existing dense clumps, dynamical instabilities of the ionization front, and the collapse of the collected layer may explain the variety of morphologies observed in this region. We conclude that the S284 fulfils the necessary conditions for an in-depth study of triggered star-formation by expansion of its individual H  regions.

[^1]: Based on data obtained with IRAC onboard Spitzer (Program ID 3340) and the WFC at Isaac Newton Telescope at La Palma.