{'fraction_non_alphanumeric': 0.03650611928987213, 'fraction_numerical': 0.050018543302642826, 'mean_word_length': 4.45837000528076, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A comprehensive study of young B stars in NGC 2264: I. Space photometry and asteroseismology

Space photometric time series of the most massive members of the young open cluster NGC 2264 allow us to study their different sources of variability down to the millimagnitude level and permits a search for Slowly Pulsating B (SPB) type pulsation among objects that are only a few million years old. Our goal is to conduct a homogeneous study of young B type stars in the cluster NGC 2264 using photometric time series from space in combination with high-resolution spectroscopy and spectropolarimetry obtained from the ground. The latter will be presented in a separate follow-up article. We performed frequency analyses for eleven B stars in the field of the young cluster NGC 2264 using photometric time series from the MOST, CoRoT and Spitzer space telescopes and the routines Period04 and SigSpec. We employ the MESA stellar evolution code in combination with the oscillation code GYRE to identify the pulsation modes for two SPB stars which exhibit short period spacing series. From our analysis we identify four objects that show SPB pulsations, five stars that show rotational modulation of their light curves caused by spots, one star that is identified to be a binary, and one object in the field of the cluster that is found to be a non-member Be star. In two SPB stars we detect a number of regularly spaced pulsation modes that are compatible with being members of a g mode period series. Despite NGC 2264's young age, our analysis illustrates that its B type members have already arrived on the zero-age main sequence (ZAMS). Our asteroseismic analysis yields masses between 4 and 6 Msun and ages between 1 and 6 million years, which agree well to the overall cluster age.


Introduction
Stellar rotation, pulsation, binarity and magnetic fields play important roles in massive stars influencing their structure, evolution and lifetimes across all evolutionary phases, from star formation until the latest stages (e.g., Langer 2012; Maeder & Meynet 2012;Aerts et al. 2010). Our theoretical concepts for the evolutionary properties of O and B stars have strong implications on the predictions for supernovae progenitors, galaxies and the chemical enrichment of the whole universe. Stellar evolution codes of single and binary stars include parametrized descriptions of phenomena crucial for massive stars, such as internal differential rotation, transport of angular momentum, convective core overshooting, chemical composition, and magnetic fields. However, our theoretical concepts mostly lack essential calibration by high-precision observations. First successful seismic studies revealed, for example, the internal rotation profile   or put tight constraints on core overshooting and diffusive mixing of a massive pulsating star based on observations conducted with the Kepler space telescope (Borucki et al. 2010), hence illustrating the importance and power of testing our theoretical predictions with precise observational data.
Slowly pulsating B (SPB) stars (Waelkens 1991) have spectral types between B3 and B9 corresponding to ∼ 11 000 to 22 000 K in effective temperature (T eff ) and masses ranging from 2.5 to 8 M (see e.g., Aerts et al. 2010, Chapter 2). They show non-radial, multi-periodic, high-order gravity (g) mode oscillations with periods between 0.5 to 3 days which are driven by the κ-mechanism connected to the iron-group elements opacity bump (Dziembowski et al. 1993). The SPB stars' g-modes of same degree and consecutive radial order n are expected to be nearly equally spaced in period. The deviations from this equal spacing carry information about the physical processes occurring near the core (Miglio et al. 2008).
This theoretical prediction was first investigated with observational data by Pápics et al. (2015), who detected a long series of g-modes with nearly equidistant spacings in period space for an SPB star which revealed clear signatures of chemical mixing and rotation. Moreover, asteroseismic modeling allowed ages for several SPB stars to be deduced from their central hydrogen mass fractions, X c (e.g., Moravveji et al. 2015Moravveji et al. , 2016, always under the assumption that the stars are in their main sequence evolutionary phases, i.e., burning hydrogen in their cores. Gruber et al. (2012) suggested for the first time the presence of SPB stars among pre-main sequence (pre-MS) members of the young open cluster NGC 2244 based on an asteroseismic analysis. Stars with masses higher than ∼6M do not have an optically observable pre-MS phase as the birthline intersects the ZAMS (Palla & Stahler 1990, 1993. B type stars with masses lower than ∼6M have a relatively short optically visible pre-MS evolutionary phase lasting at most a few million years. Generally, pre-MS stars mainly gain their energy from gravitational contraction until their cores become hot enough for hydrogen burning to ignite. As pre-MS SPB stars are already quite close to their arrival on the ZAMS, they are in a transition phase between the end of gravitational contraction as main energy source and the onset of hydrogen core burning. During this phase, stars undergo significant structural changes that shall be probed through asteroseismology. Only with the first observations of potential pre-MS SPB candidates from space using the MOST satellite (Walker et al. 2003), the search for and the investigation of this type of young pulsating objects could commence. For a full interpretation of the stars' photometric variability, for any asteroseismic interpretation and for tests of stellar evolution models, an accurate characterization of their fundamental parameters such as effective temperature, surface gravity, chemical abundances and magnetic fields is crucial.
The present study is the first of two articles prepared simultaneously that use the time series photometry for eleven B stars in the field of NGC 2264 in combination with high-resolution spectroscopy and spectropolarimetry for a comprehensive study of the young B stars in NGC 2264. In this paper (Paper I hereafter) we discuss the space photometry obtained with the MOST space telescope (Walker et al. 2003(Walker et al. ) in 2006(Walker et al. and 2011, with the CoRoT satellite (COnvection, ROtation and Planetary Transits; Baglin & Fridlund 2006) in 2008 and with the NASA Spitzer (Werner et al. 2004) satellite obtained during the CSI2264 (Coordinated Synoptic Investigation of NGC 2264) project (Cody et al. 2013) simultaneously to the 2011/12 MOST (and CoRoT) observations (Sect. 3). We characterize the different types of variability (Sect. 3), conduct an asteroseismic investigation for two selected young SPB stars showing indications for period spacings, and discuss the stars' evolutionary stage (Sect. 4).
The subsequent paper (Przybilla et al. 2017, in preparation; Paper II hereafter) will focus on the analysis of the highresolution spectroscopy and spectropolarimetry obtained with the Mc Donald Observatory 2.7 m telescope and the Tull spectrograph and with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) and describe the stars' fundamental parameters, atmospheric chemical abundances, and magnetic properties.

B stars in the field of NGC 2264
NGC 2264 (α 2000 = 06 h 41 m , δ 2000 = +09 • 35') is a very young (age ∼3 -8 Myr; e.g., Sung et al. 2004;Sagar et al. 1986) open cluster that is rich in young stellar objects and has been studied frequently in the past using various instruments in different wavelength ranges from space and from the ground. It has a diameter of ∼39 arcmin and is located in the Monoceros OB1 association about 30pc above the Galactic plane. The cluster distance is reported to be 759 ± 83 pc by Sung et al. (1997) corresponding to a distance modulus of 9.40 ± 0.25 mag. The proper motion values of NGC 2264 are −2.70 ± 0.25 mas/yr in right ascension and −3.50 ± 0.26 mas/yr in declination (Kharchenko 2001). The cluster reddening is rather low with E(B−V) = 0.071±0.033 mag and the differential reddening is negligible (Sung et al. 1997).
NGC 2264 is located in a region accessible to the space telescopes MOST, CoRoT, Spitzer and Chandra and was the target for the Coordinated Synoptic Investigation of NGC 2264 (CSI2264) project (Cody et al. 2013). The CSI2264 project lasted from the beginning of December 2011 until early January 2012 and involved simultaneous observations of NGC 2264 for more than one month using all four above mentioned space telescopes together with a coordinated ground-based observing campaign for complementary data. One of the main scientific topics addressed with these combined observations from space and from the ground was the discovery and analysis of pulsationally variable stars in the pre-and early main sequence stages. Zwintz et al. (2009) presented the first investigation of 68 stars in the field of NGC 2264 observed by the MOST satellite in 2006. They identified four pre-main sequence (pre-MS) δ Scuti type stars, one γ Doradus candidate, ten SPB candidates, and six stars showing rotational modulation caused by spots on their surfaces. A previous investigation of the two early B-type stars HD 47887 and HD 47777 in NGC 2264 revealed that both stars possess large-scale magnetic fields with average longitudinal field strengths of about 400 G Fossati et al. (2014). For completeness of the sample of B stars in NGC 2264 observed from space, we also include these two objects in the present study. In the following investigation we will refer to the results published in that paper and put both stars into context to the nine other stars.
The stars HD 47469, HD 48012, HD 261810, HD 261878, HD 47777, HD 47887, HD 261903, HD 261938, NGC 2264 137, HD 261054 and HD 47961 (numbered from 1 to 11 accordingly; see also Table 1) lie in the magnitude range between 7.1 and 9.9 mag in the V band, and -given their corresponding smallest B − V values -they belong to the hottest stars in NGC 2264. Note that the spectral type of "K1?" given by the SIMBAD data base for HD 261903 was derived from multi-color photometry and is likely a mis-identification: As HD 261903 is located in a relatively obscured region in NGC 2264 (which can also be identified in Figure 2), its colors would indeed point to a cool object, but its spectrum clearly indicates a spectral type B. The top panel of Fig. 1 shows the colour-magnitude diagram (CMD) of NGC 2264 in Johnson V and B photometry down to V = 16 mag, where the eleven B type stars are marked with black symbols. A zoom into the hottest part of the CMD (Fig. 1 bottom panel) makes it easier to identify the stars that are subject to this investigation.
All stars, except for NGC 2264 137, have been reported to be candidates for SPB type pulsations showing frequencies between 0.05d −1 and 6.5d −1 by Zwintz et al. (2009), but the authors did not describe the frequency content in detail. NGC 2264 137 (GSC 00750-01804) did not show clear variability in the MOST 2006 observations, hence was listed as a constant object (Zwintz et al. 2009). The star HD 48012 is listed in the CHAN-DRA Variable Guide Star Catalog as a pulsating variable star and UV source (Nichols et al. 2010). The stars HD 261810, HD 261903, HD 261938 and HD 47961 are reported to be X-ray sources from ROSAT measurements (Flaccomio et al. 2000) and HD 261810 and NGC 2264 137 are included in the X-ray observations by Dahm et al. (2007). The stars HD 47777 and HD 47887 have been reported to show variability caused by rotational modulation of a spotted surface and host large scale magnetic fields (Fossati et al. 2014).
For HD 261810, HD 261878, HD 261903, HD 261938, NGC 2264 137 and HD 47961 membership probabilities (see Table  1) ranging between 25% (for HD 261903) and 97% (for HD 261938) were taken from Vasilevskis et al. (1965) and Flaccomio et al. (2000). All stars were also included in the CSI2264 project where for HD 47469 and HD 261054 no membership information is given and for all other stars a very likely to possible cluster membership is reported (see Table 1). From radial velocity measurements, Fossati et al. (2014) concluded that HD 47777 and HD 47887 are cluster members.
Four additional criteria can be used to assess the cluster membership of stars: (i) their position on the sky with respect to the center of the cluster, (ii) their proper motion values in comparison to the cluster proper motion, (iii) the location of the stars in the cluster colour-magnitude diagram (CMD), and (iv) the stars' radial velocities. The location of the eleven stars on the sky is displayed in Fig. 2, which shows a one square degree field centered at the coordinates RA 2000 = 06:39:43.5 and DE 2000 = 09:40:50.3. It can clearly be seen that the two stars labeled #1 (HD 47469) and #10 (HD 261054) are located further away from the cluster center compared to the other nine objects. However, HD 47469's CMD position matches quite well those of the other stars with high membership probability (e.g., NGC 2264 137 with 96%) in its vicinity. Additionally, its proper motion agrees to the cluster values within the given errors (Table 1). Hence, we find more arguments to support HD 47469's membership to NGC 2264 than against it. On the other hand, HD 261054, is positioned with an offset in (B − V) in the cluster CMD ( Fig.  1) and has a proper motion that is significantly outside the errors of the cluster values. We therefore conclude, that HD 261054 is more likely to be a background field object. The radial velocities of all eleven stars have been determined homogeneously from our spectroscopic data and will be discussed in detail in Paper II (Przybilla et al. 2017 Hence, we conclude that in our sample ten stars are likely to be cluster members and only HD 261054's membership to NGC 2264 is unlikely.

Photometric analysis
The photometric time series presented in this study have been obtained using the three space telescopes MOST, CoRoT and Spitzer. Table 2 summarizes the observational material.

MOST observations and data reduction
The Canadian micro-satellite MOST carries a 15-cm Rumak Maksutov telescope which feeds a CCD photometer through a single custom broadband filter in the 3500-7500Å wavelength range. Three types of photometric data can be supplied simultaneously for multiple targets in a given field: Fabry Imaging, Direct Imaging and Guide Star Photometry data. Although the Canadian Space Agency ceased funding for MOST in 2015, the space telescope is still operational and can be used for dedicated Article number, page 3 of 19 A&A proofs: manuscript no. NGC2246Bstars_phot Table 1. Properties of the eleven B stars in the field of NGC 2264 and comparison to the cluster position and proper motions (Kharchenko 2001 Vasilevskis et al. (1965); reference # 2: Flaccomio et al. (2000) and membership as provided from the CSI2264 project through the NASA / IPAC Infrared Science Archive a (CSI2264), where 0 -no membership information, 1 -very likely cluster member, 2 -possible cluster member. Cluster proper motion values are taken from Kharchenko (2001). * : The spectral type of K1? for HD 261903 is likely a misidentification based on multi-color photometry of a star located in an obscured region of the cluster.  (Zwintz et al. 2009). The second MOST observing run on NGC 2264 lasted from December 5, 2011, to January 14, 2012, and was conducted as part of the CSI 2264 project.
Due to the magnitude range (7 < V < 12 mag) and the large number of targets, in both runs NGC 2264 was observed in the open field of the MOST Science CCD in Guide Star Photometry Mode. Because not all of the targets could be reached using a single pointing of the satellite, two fields of observations were chosen and observed in alternating halves of each 101-min orbit. Notes. Column headers are a running number as in Table 2 (#), identifier (Name), MOST guide star numbers and fields in 2006 (MOST 2006/12 (MOST 2011, availability of additional CoRoT or Spitzer data (other), and identified type of variability (var): slowly pulsating B type (SPB), rotational modulation (rot.mod), Be star (Be) and spectroscopic binary (SB).
Using this setting, MOST time series photometry was obtained for a total of 68 stars in the region of NGC 2264 in the 2006 run (Zwintz et al. 2009) and 67 in the 2011/12 run.
Data reduction of the MOST Guide Star photometry was conducted using the method developed by Hareter et al. (2008) and is described in detail in Zwintz et al. (2009). The corresponding MOST 2006 light curves of both fields -field A and field B -used for the analysis each have a time base of ∼22.7 d. The on-board exposure time was 1.5 s, 16 consecutive images were "stacked" on board resulting in an integration time of 24 s (see also Zwintz et al. 2009). The total length of the NGC 2264 data set from the 2011/12 run is 38.77 d for field 1 and 39.98 d for field 2. In 2011/12 on-board exposures were 3.01 s long (to satisfy the cadence of guide star ACS operations), and the integration time was 51.17 s as 17 consecutive images were "stacked" on board. An overview of the properties of the MOST data is given in Table 3. The corresponding spectral windows (top panels: 2006 data, bottom panels: 2011/12 data) are given in Fig.  3. Notes. Data set, information on the observations (obs. info), number of data points used for the analysis (N), time base (t base ), Rayleigh frequency resolution (1/T ), and Nyquist frequency (f Nyquist ).

CoRoT observations and data reduction
On December 27, 2006, the European satellite CoRoT was launched from Baikonur, Kazakhstan, into a circular polar orbit of about 896 km altitude. It carried a 27-cm telescope and observed stars inside two cones with a radius of 10 • centered around RA = 06:50 and RA = 18:50. The telescope's field-ofview is almost circular with a diameter of 3.8 • . The filter bandwidth ranges from 370 to 1000 nm. The CoRoT space telescope originally had two CCDs devoted to asteroseismology for stars with 5.7 < m V < 9.5 mag and two CCDs dedicated to the search for exoplanets where ∼6000 stars in the magnitude range from 10 to 16 mag in R per CCD were monitored simultaneously. Since a hardware failure in November 2012, CoRoT is not operational any more. The young cluster NGC 2264 was targeted by CoRoT for the first time during the Short Run SRa01 in 2008 for 23.4 days consecutively and for the second time during the Short Run SRa05 in 2011/12 as part of the previously mentioned CSI 2264 multisatellite campaign for a total time base of ∼39 days. In both observing runs, the complete cluster was placed in one exofield CCD and data were taken for all stars in the accessible magnitude range.
Therefore, the brightest stars in NGC 2264 (i.e., with V < 10 mag) have not been targets for the CoRoT exofield observations. But for selected stars data could be obtained using the method of the so-called CoRoT imagettes (for a detailed description see, e.g., Zwintz et al. 2011). NGC 2264 137 is the only B-type star in the field of the cluster that has been included as an imagette in the CoRoT SRa01 observations of NGC 2264. Due to the special observations of a star that is nominally too bright and saturates on the CCD, the light curve of NGC 2264 137 is only ∼ 19 days long (see Table 3).

Spitzer observations
During the CSI2264 campaign, the Spitzer satellite observed the young cluster NGC 2264 for roughly 30 days in the mid-infrared from the beginning of December 2011 to beginning of January 2012. Two 5.2' x 5.2' regions near the center of the cluster were observed twice each for about one day with Spitzer's high precision staring mode (< 1%) with a cadence of 0.1 or 1 minute (Cody et al. 2013). Data for the remaining targets were taken in mapping mode (1 -3% precision) with a cadence of 100 minutes (Cody et al. 2013).
HD 261903 is the only B type star that was observed in one of the two regions near the cluster center in Spitzer staring mode using the IRAC2 band at 4.5µm. The total time base of the IRAC2 light curve is 6.01 days where observations were taken in four parts of one day length each with two gaps of one day length each in between the observations. Hence, the resolution of the Spitzer light curve is sufficiently high to study variability on time scales of a few hours.
The IDL package Cluster Grinder (Gutermuth et al. 2009) was used to generate the light curves from the basic calibrated data (BCD) images released by the Spitzer Science Center (SSC) pipeline (version S19.1). Cluster Grinder provides mosaics for each Astronomical Observation Request (AOR), point source locations and photometric measurements from the mosaics. From both the BCDs and the mosaics, aperture photometry was computed using an aperture radius of 2.4" (2 pixels), and a sky annulus from 2.4" to 8.4" (2-6 pixels). A detailed description of the Spitzer data reduction can be found in Cody et al. (2014).

Frequency analysis
For the frequency analyses, we used the software package Pe-riod04 (Lenz & Breger 2005) that combines Fourier and leastsquares algorithms. Frequencies were then prewhitened and considered to be significant if their amplitudes exceeded four times the local noise level in the amplitude spectrum (Breger et al. 1993;Kuschnig et al. 1997).
We verified the analysis using the SigSpec software (Reegen 2007). SigSpec computes significance levels for amplitude spectra of time series with arbitrary time sampling. The probability density function of a given amplitude level is solved analytically and the solution includes dependences on the frequency and phase of the signal.

HD 47469
HD 47469 was observed by the MOST space telescope in 2006 and 2011/12 (Table 2) and found to be an SPB type star (Zwintz et al. 2009). Using both data sets, we could identify 12 pulsation frequencies (see Fig. 4 and Table 4). Ten of those were found in both data sets and two additional frequencies could be identified due to the higher frequency resolution only in the 2011/12 data set. We searched for linear combinations and identified F6 at 2.764d −1 to be three times F3 at 0.922d −1 . Apart from this, all other frequencies are independent from each other.
Frequency F11 at 1.226d −1 decreased its amplitude from 2006 to 2011/12 significantly by about 0.38 mmag. The amplitude of frequency F4 also decreased between the two observing runs with MOST, but this difference is likely a result of the different frequency resolution: in 2006 the two close frequencies F8 at 2.873 d −1 and F10 2.892 d −1 could not be resolved, hence the corresponding frequency of 2.892d −1 (F10 in 2006) includes the power of the two peaks that were identified separately in the 2011/12 data set. The amplitudes of F1 at 1.621d −1 , F5 at 1.377d −1 , and F6 at 2.764d −1 increased significantly between the 2006 and the 2011/12 observations. Amplitude variability has been reported previously, e.g., for δ Scuti type stars (e.g., Bowman & Kurtz 2014). With the present data sets we can only report about the observed change in amplitude for F1, F5 and F6, but cannot conduct a detailed study of their temporal behavior due to the relatively short MOST data sets available to us.  Longer photometric time series are needed for a more thorough interpretation of HD 47469's amplitude variability.

HD 48012
MOST took observations of HD 48012 only in 2006 (Table 2). Figure 5 shows the original amplitude spectrum of HD 48012 from the 2006 observations with the six identified pulsation frequencies (black solid lines) which are indicative for SPB type pulsation (see Table 5). It can clearly be seen that additional frequencies in the range between ∼2.5 and 3.0d −1 are buried in the noise.

HD261810
HD 261810 was observed by the MOST space telescope in 2006 and in 2011/12 (Table 2) and identified as SPB pulsator (Zwintz et al. 2009). Our frequency analysis revealed 13 intrinsic frequencies that can be divided in two groups: the lowest frequency range -hereafter named Group 1 -lower than 1d −1 , and the higher frequency range -hereafter named Group 2 -from about 1.8 to 5 d −1 . Figure 6 illustrates the original amplitude spectra for the 2006 data (pointing downwards) and the 2011/12 data (pointing upwards). The identified intrinsic frequencies are marked with black solid lines.  Notes. Pulsation frequencies, periods, amplitudes, signal-to-noise values and SigSpec significances identified from the MOST data sets from 2006 (amp 1 , S/N 1 ,sig 1 ) and 2011/12 (amp 2 , S/N 2 , sig 2 ) as well as the the corresponding amplitude differences (amp 2 -amp 1 ). The last-digit errors of the frequencies and periods are given in parentheses; frequency errors were computed according to Kallinger et al. (2008) and the errors in periods were derived using standard error propagation. Possible linear combinations are given in the last column (lin. combi). Notes. Pulsation frequencies, periods, amplitudes, signal-to-noise values and SigSpec significances identified from the MOST data set from 2006 (amp 1 , S/N 1 ,sig 1 ) The last-digit errors of the frequencies and periods are given in parentheses; frequency errors were computed according to Kallinger et al. (2008) and the errors in periods were derived using standard error propagation. Frequencies in Group 1 can all be explained as linear combinations of higher frequencies (see Table 6). We clearly identify frequencies in Group 2 to be caused by SPB type pulsations.

HD 261878
Observations for HD 261878 were taken in 2006 and in 2011/12 by the MOST satellite ( Table 2). The frequency analysis of the 2006 data set revealed no formally significant frequency while the 2011/12 data set shows six frequencies that can be attributed to SPB type pulsation. The reason for this difference between the two data sets is the rather high noise level in the 2006 data (see Fig. 7). The average noise level of the original 2006 data set was calculated to be 137 ppm in the range from 0 to 100d −1 while for the 2011/12 data set, the same average noise lies at 39 ppm. A comparison of the amplitude spectra for both MOST observing runs shows that the frequencies we identified from the 2011/12 data sets are buried in the noise of the 2006 data set.

HD 47777
The star HD 47777 was discussed in detail by Fossati et al. (2014). Here we summarize HD 47777's properties for completeness.
The MOST data for HD 47777 obtained in 2006 and 2011/12 revealed a single intrinsic frequency at 0.3787d −1 corresponding Article number, page 7 of 19 A&A proofs: manuscript no. NGC2246Bstars_phot Notes. Same format as in Table 4. Notes. Pulsation frequencies, periods, amplitudes, signal-to-noise values and SigSpec significances identified from the MOST data set from 2011 (amp 2 , S/N 2 ,sig 2 ) The last-digit errors of the frequencies and periods are given in parentheses; frequency errors were computed according to Kallinger et al. (2008) and the errors in periods were derived using standard error propagation.
to a rotation period of 2.641 days. The non-LTE analysis of highresolution spectroscopic data for HD 47777 yielded an effective temperature, T eff , of 22000±1000 K, a surface gravity, log g, of 4.2±0.1 and a projected rotational velocity, υ sin i, of 60 km s −1 (Fossati et al. 2014). Using spectropolarimetric data obtained with the ESPaDOnS spectropolarimeter of the Canada-France-Hawaii Telescope (CFHT) the authors reported the detection of a magnetic field with a strength of the longitudinal component of 469±87 G (Fossati et al. 2014).

HD 47887
The MOST space telescope observed HD 47887 only in 2006. A detailed discussion on HD 47887 is given in Fossati et al. (2014), hence, here we only give a brief summary of the star's properties. HD 47887 shows a single intrinsic frequency at 0.513d −1 corresponding to a period of 1.949 days. From highresolution spectroscopy, Fossati et al. (2014) found T eff to be 24000±1000 K, log g as 4.1±0.1 and υ sin i as 45 km s −1 . Spectropolarimetric measurements taken with the ESPaDOnS spectropolarimeter of the Canada-France-Hawaii Telescope (CFHT) revealed the presence of a magnetic field with a longitudinal component of 373±48 G (Fossati et al. 2014).

HD 261903
For HD 261903 we were able to investigate two MOST data sets and an infrared light curve obtained by the Spitzer space tele-scope (Table 2). Figure 8 shows the respective amplitude spectra from the MOST 2006 (upper panel pointing downwards) and 2011/12 (upper panel pointing upwards) data sets where the two frequencies F1 at 3.216 d −1 and F2 at 3.275 d −1 (see Table 8) are marked as solid lines. Figure 8 shows the amplitude spectrum of the IRAC2 light curve (bottom panel). It can clearly be seen that the frequency resolution is not sufficient to resolve both frequencies found in the MOST data. With only two very close frequencies in a rapidly rotating star (see Paper II), we interpret the variability to be caused by differential rotation of a spotted star.

HD 261938
The analysis of the two MOST data sets for HD 261938 (see Table 2) revealed two significant frequencies: F1 at 0.71845d −1 , and F2 at 6.50342d −1 (see Fig. 9). F1 corresponds to a period, P1 of about 1.392d −1 , which could be caused by pulsations, binarity or rotational modulation of a spotted surface. F2 on the other hand is much too high to be explained as, e.g., SPB type pulsation. But a closer inspection revealed that F2 is nine times F1 within the given errors. With the present observational material we cannot decide whether F2 is indeed the ninth harmonic of F1 and the "missing" harmonics are buried in the noise or whether F2 has a different origin. Hence, we cannot draw a final conclusion of HD 261938's variability. We can only speculate that we might either see rotational modulation of a spotted surface or that we might observe a binary system with tidally induced pulsations.

NGC 2264 137
NGC 2264 137 was observed by MOST in 2006 and 2011/12; additionally it was also included as imagette in the CoRoT Short Run SRa01 in 2008.
For NGC 2264 137 we identified five formally significant frequencies from the CoRoT data set where four can be explained as multiples of the frequency with the highest amplitude (F1) at 0.471d −1 as illustrated in the bottom panel of Fig. 10. F1 also appears in the MOST data set obtained in 2011/12 (top panel in Fig. 10). The MOST 2006 data set unfortunately has a significantly higher noise level which makes it impossible to identify a frequency with an amplitude at the millimagnitude level.  Notes. Intrinsic frequencies, amplitudes and signal-to-noise values identified from the MOST data sets from 2006 (amp 1 , ,phase 1 S/N 1 ) and 2011/12 (amp 2 , phase 2 , S/N 2 ) and the Spitzer IRAC2 data sets from 2011/12 (amp 3 , phase 3 , S/N 3 ). Frequency errors were computed according to Kallinger et al. (2008). Notes. Pulsation frequencies, periods, amplitudes, signal-to-noise values and SigSpec significances identified from the MOST data set from 2006 (amp 1 , S/N 1 ,sig 1 ) and 2011/12 (amp 2 , S/N 2 ,sig 2 ). The last-digit errors of the frequencies and periods are given in parentheses; frequency errors were computed according to Kallinger et al. (2008) and the errors in periods were derived using standard error propagation.
face. The rotational period would then be the inverse of F1, i.e., 2.123 days.

HD 261054
MOST time series photometry was taken for HD 261054 only in 2006 (Table 2). The frequency analysis revealed only three significant intrinsic frequencies (see Table 11). Although F3 has a value of nearly exactly 2d −1 , we consider it intrinsic to the star and not related to the properties of the data set as the spectral window function of the corresponding field B of the MOST 2006 observations does not show a signal at 2d −1 (see Fig. 3). The multi-sine fit to the data does not represent the observed variability perfectly, which leads us to suggest that additional frequencies might be buried in the noise and could be revealed with longer time series. With such a low number of detected frequencies the origin of the star's variability might be connected to pulsations or caused by rotational modulation due to spots. HD 261054 is very likely not a member of NGC 2264 because it is located furthest away from the cluster center and at a larger distance (see Section 2; Gaia Collaboration et al. 2016), has an offset position in the cluster CMD compared to the other B type stars and its proper motion values are significantly different from the overall cluster values as discussed earlier. As HD 261054 is listed as emission line star in the literature, its variations are likely to be connected to the Be phenomenon (e.g., Rivinius et al. 2013). For a more detailed study of HD 261054 and a possible interpretation as a Be star, longer photometric time series are required.

HD 47961
MOST observed HD 47961 in 2006 and in both fields in 2011/12. As field 1 and 2 were observed alternatingly in 2011/12, we merged the two data sets and ran the frequency analysis on the complete time series which consists of 25059 data points obtained during a total time base of 39.98 days. In the 2006 data set we find one significant frequency at ∼0.055d −1 and several peaks at similar amplitudes in the region between 0.77 and 1.35d −1 that are not formally significant (see Fig. 13). As the complete 2011/12 data set has a lower noise level and smaller frequency resolution, we can identify the low frequency peak in fact to consist of two frequencies F1 at 0.027d −1 and F2 at 0.053d −1 , where F2 is twice F1, and we find a third frequency F3 at 0.924d −1 . The origin of the photometric variability per se is inconclusive. Only in combination with the spectroscopic data are we able to suggest the source of the light variations as binarity (see Paper II).

Fundamental parameters
To conduct an evolutionary and asteroseismic inference, it is crucial to determine the stars' fundamental parameters. Therefore, we obtained high-resolution, high signal-to-noise spectroscopic and spectropolarimetric data for all eleven B stars in NGC 2264. A detailed analysis of all data will be discussed in Paper II. Here we only describe the determination of the atmospheric parameters -T eff , log g, υ sin i -for the four SPB type pulsators that are used as input for our asteroseismic analysis.
High-resolution spectroscopic data for HD 47469, HD 48012 and HD 261810 were obtained on 2012, January 31 with the ESPaDOnS spectropolarimeter of the Canada-France-Hawaii Telescope (CFHT). ESPaDOnS consists of a table-top, crossdispersed échelle spectrograph fed via a double optical fiber directly from a Cassegrain-mounted polarisation analysis module. Both Stokes I and V spectra were obtained throughout the 3700-10400 Å spectral range at a resolution of about 65 000. The spectra were reduced using the Libre-ESpRIT reduction pipeline (Donati et al. 1997). The spectra were then normalized by fitting a low-order polynomial to carefully selected continuum points. The signal-to-noise (S/N) ratios of the spectra for HD 47469, HD 48012 and HD 261810 were 275, 162 and 261 per 1.8km s −1 , respectively, measured at a reference wavelength of 5000 Å.
The spectrum for HD 261878 was observed on 2011, December 31 with the Robert G. Tull Coudé Spectrograph (TS), which is mounted on the 2.7-m telescope of Mc Donald Observatory and has a resolving power R ∼ 52 000. The spectra A&A proofs: manuscript no. NGC2246Bstars_phot Notes. Significant frequencies, amplitudes, signal-to-noise values and SigSpec significances identified from the CoRoT data set from 2008 (amp C , S/N C ,sig C ) and the MOST data set from 2011/12 (amp 2 , S/N 2 ,sig 2 ), and the identified multiples of F1 (multiples). The last-digit errors of the frequencies and periods are given in parentheses; frequency errors were computed according to Kallinger et al. (2008) and the errors in periods were derived using standard error propagation. cover a wavelength range of 3633-10849 Å with gaps between the échelle orders at wavelengths longer than 5880 Å. Bias and flat field frames were obtained at the beginning of each night, while several Th-Ar comparison lamp spectra were obtained each night for wavelength calibration purposes. The reduction was performed using the Image Reduction and Analysis Facil-  Notes. Same format as Table 5. ity 1 (IRAF). The spectra were normalised by fitting a low order polynomial to carefully selected continuum points. As Hα is not covered in TS spectra and Hβ is affected by a reflection in the imaging system of the spectrograph, we employed the Hγ line in the atmospheric parameter determination in HD 261878. We normalized the Hγ line continuum using the artificial flat-fielding technique described in Barklem et al. (2002), which has proven to be successful with TS data (e.g., Zwintz et al. 2013). The S/N of HD 261878' s spectrum is 101 per pixel at λ ∼5000Å.
For the determination of the atmospheric parameters, we computed synthetic spectra with a hybrid non-LTE approach in analogy to Nieva & Przybilla (2012). This is based on atmospheric structures calculated with the Atlas9 code (Kurucz  Fig. 11. Light curve of HD 261054 (grey points) and multi-sine fit with the three intrinsic frequencies (solid black line). 1993). Line-formation calculations in non-LTE were performed with recent versions of Detail and Surface (both updated by K. Giddings 1981;Butler & Giddings 1985). The validity of the approach was verified by direct comparison with hydrostatic and hydrodynamic line-blanketed non-LTE model atmospheres (Nieva & Przybilla 2007;Nieva & Simón-Díaz 2011;Przybilla et al. 2011). Non-LTE level populations and synthetic spectra of all elements were computed using updated model atoms as summarized in Table 2 of Przybilla et al. (2016). An additional challenge for the analysis of mid-/late B-type stars as discussed here are the overall reduced number of metal lines, which are also typically much weaker than in the earlier B-type stars. Moreover, many of the mid-/late B-type stars show intermediate to large υ sin i values, such that often only a small number of metal lines remains for the quantitative analysis besides the H i and He i lines. The only metal ionization equilibria available for the atmospheric parameter determination (showing sensitivity to both T eff and log g variations) are those of O i/ii and Si ii/iii. This renders the atmospheric parameter determination less accurate than that of sharp-lined and hotter stars (like those analyzed by Nieva & Przybilla 2012), for which the redundancy from additional ionization equilibria provides tighter constraints. The υ sin i values were derived from fits of the model spectra to observed metal lines.
The resulting atmospheric parameters of the four SPB type pulsators of our star sample are summarized in Table 12. Examples for the quality of the fits to the observed spectra of the four stars in the region of Hγ are shown in Figure 14. In the course of our analysis, we identified HD 47469 to be a binary system. We therefore list the atmospheric parameters of both components in Table 12. From the present data it is impossible to identify if one of the two components is pulsating or if both are SPB stars. The detailed analysis of the non-LTE quantitative spectral analysis will be discussed in Paper II.  HD 261810 17200 ± 400 4.25 ± 0.10 180 ± 15 5 HD 261878 15200 ± 400 4.30 ± 0.10 40 ± 5

MESA Evolutionary Tracks
We employed the state-of-the-art one-dimensional stellar structure and evolution code MESA (Paxton et al. 2011(Paxton et al. , 2013(Paxton et al. , 2015 with the choice of physical input based on thorough asteroseismic inferences from the modeling of two Kepler B8-type SPB stars KIC 10526294  and KIC 7760680 ). We computed a coarse grid of non-rotating models by varying the initial mass M ini from 4.0 to 8.0 M , in steps of 0.5 M . MESA incorporates the effect of rotation on the equilibrium structure of stars by solving the modified structure equations presented by Endal & Sofia (1976, 1978, including the effects of rotational mixing by several instabilities proposed by Heger et al. (2000). The latter formalism contains a handful of unconstrained scaling coefficients which are non-trivial to constrain (see e.g. Brott et al. 2011). Considering that NGC 2264 is 3 to 8 Myr old, the effect of rotation on redistributing the internal angular momentum in the target stars is negligible. Given the limited data at our disposal for each star, and the fact that the pulsation modes of these stars are unidentified, it is inconclusive to incorporate the effect of rotation on the structure; instead, we use non-rotating models, and assume the stars are uniformly rotating. The effect of rotation is only incorporated on the pulsations. For KIC 7760680, Moravveji (2016) demonstrated that differentially rotating models fail to reproduce the observed period spacing. The evolution of each model began from the pre-main sequence (pre-MS) phase, and was terminated as soon as it reached the terminal age main sequence (TAMS, X core = 10 −6 ). We adopted the Asplund et al. (2009) initial mixture with the initial composition (X ini , Y ini , Z ini ) = (0.710, 0.276, 0.014), consistent with the Cosmic Galactic Standard composition of Nieva & Przybilla (2012) for early-type stars in the solar neighbourhood. We included convective core overshooting during the main sequence (MS) phase with the exponential overshoot prescription of Freytag et al. (1996) and Herwig (2000), which outperforms the classical step function prescription in reproducing the high precision g-mode frequencies of the above two Kepler stars (see Fig. 15. Location of the four SPB type cluster members in the Kiel diagram. Each star is given the number according to Table 12. The pre-MS contraction phase of each track is shown with solid lines, and the main sequence is marked by dotted lines for the range between 4 and 8M in steps of 0.5M . The positions of the four SPB stars are so close to the ZAMS that even our asteroseismic modeling cannot discriminate whether or not they are pre-MS stars. e.g., Moravveji et al. 2015Moravveji et al. , 2016 . The free overshooting parameter was kept fixed at f ov = 0.025. For the same reason, we applied a minimum extra diffusive mixing in the radiative interior (beyond the overshoot region), with the value D ext = 100 (cm 2 s −1 ), which is close to the two values found for KIC 10526294 and KIC 7760680. Along each evolutionary track, we stored an equilibrium model every ∼ 100 K drop in effective temperature T eff , provided T eff of each model was in the range 11 000 to 18 000 K. The Kiel diagram (i.e., T eff versus log g; see Fig. 15) shows the pre-MS (solid) and main sequence (MS, dotted) evolutionary tracks where the positions of the four NGC 2264 SPB candidates derived from our non-LTE analysis are marked with the numbers given in Table 12.
Clearly, all candidate stars lie close to the ZAMS. The pre-MS evolution occurs on the Kelvin-Helmholtz contraction timescale, and the MS evolution occurs on the nuclear timescale. Given the large difference between these two timescales, the likelihood for these stars being on the pre-MS phase is slim. Consequently, in the following seismic modeling, we adopt MS models to represent HD 261810 and HD 47469.
using the spectroscopically determined fundamental parameters (see Table 12; for more details see Paper II), with Based on this grid, we found the closest model to the HRdiagram position of each star, using the MS part of evolutionary tracks to assess the model ages, accordingly. Table 13 gives a summary of the inferred initial masses, ages and core hydrogen mass fraction X c depending on the star being in the pre-MS or MS phase. The inferred age spread is narrow, even with the ambiguity in the stars' burning status. The two components of the HD 47469 system are expected to be coeval. From the 1σ uncertainty of T eff and log g of both components, the maximum possible age of the primary and secondary are 94 and 71 Myr, respectively. Therefore, the reported ages in Table 13 are agreeably close to one another.
The basic analysis presented above confirms previous assumptions that SPB pulsators that are clearly in their pre-MS phase are statistically hard to find as they evolve quite rapidly within a few million years through their pre-MS evolutionary stage. Even in a cluster as young as NGC 2264 the late B type objects that would be candidates for pre-MS SPB pulsation seem to already have started hydrogen core burning and have to be considered early ZAMS objects. The available data for the young SPBs in NGC 2264 suffer from short duration (of < 40 days), and high noise level. Such data sets only allow a search for indications of g-mode period spacings. However, these time bases are insufficient for a more thorough asteroseismic analysis that would allow a study of possible structural differences (e.g., chemical gradients, mixing, convective core overshooting etc.) between young ZAMS and more evolved MS SPB stars. For the latter purpose, we would need period spacing series of at least a factor 2 longer that might be provided by future space missions such as TESS (Ricker et al. 2015) or PLATO (Rauer et al. 2014).

Period spacings in the SPB stars
In rotating stars, the pulsation eigenfrequeincies in the corotating frame are sensitive to rotation (Unno et al. 1989). The period spacing of zonal and prograde modes decreases monotonically with increasing rotation velocity (Townsend 2003;Bouabid et al. 2013). Using the results of our frequency analysis for the four SPB stars HD 47469, HD 48012, HD 261810 and HD 261878, we investigated the presence of regular gravity (g) mode period spacings and the deviations from their averages. We inspected the period content for each of the stars separately and computed the corresponding period spacings. Because the time bases of the MOST observations are rather short (∼ 24 and 39 days), our search for nearly equidistant period spacings is limited, and does not allow us to investigate rotational splittings of g-modes. From the four SPB stars in NGC 2264, we find indications for period spacing patterns in two stars: HD 261810 and HD 48012.
The star HD 261810 exhibits a regular period spacing structure (see Fig. 16), assuming that all modes belong to a series with identical spherical degree . The longest series we find includes five consecutive peaks, but the similar spacing can also be found between two smaller periods (around 0.2 days). The average spacing is 1495 seconds (see Fig. 16) with a maximum absolute deviation of 230 seconds.
The pulsational periods found in HD 48012 also show an average spacing of 669 seconds with a maximum absolute deviation of 120 seconds which is illustrated in Fig. 17. Figure 18 shows a comparison between the short period spacing patterns identified from the MOST data of HD 48012 (yellow squares) and HD 261810 (yellow circles) and seven stars observed by the Kepler mission (grey squares, Pápics et al. 2016). The authors discuss that stars with higher rotational velocity show their strongest excited modes and their main period series towards shorter periods, and their average spacing values are smaller. Our two young SPB stars, HD 48012 and  HD 261810, with υ sin i values of 225±20 km s −1 and 180±15 km s −1 follow that relation nicely as they belong to the fastest rotators investigated, they show the shortest periods and smallest spacing values.
No clear spacing pattern can be identified from the 12 pulsation periods found in the HD 47469 data. The reasons might also be connected to the fact that HD 47469 was discovered to be a spectroscopic binary (see Paper II) that might complicate the pulsational analysis to the extent that both components might be pulsating.
The present data for HD 261878 are too sparse to calculate an average observed spacing and its deviations. It seems that periods at lower amplitudes might be buried in the noise and could fill the gaps between the six periods identified by us. It is therefore obvious that longer photometric time series with lower noise levels are necessary to recover the possibly missing periods and reveal the complete period spacing pattern.

Asteroseismic Inference for HD 261810
As explained in Sect. 3.4.4, HD 261810 exhibits five g-modes with nearly equidistant period spacing (F1, F4, F5, F6, and F10 in Table 6). From Fig. 15, it can be seen that HD 261810 is very close to the zero-age main sequence (ZAMS), with υ sin i = 180±15 km s −1 .From the light curve analysis, it is impossible to identify the angular quantum numbers ( , m), since there are no reconcilable frequency splittings around each peak (see Fig. 6). For high rotation velocities, rotational splittings will be large and asymmetric.
To remedy this, we benefit from the available seismic information, i.e., the periods of the observed five g-modes explained above, and examine if a seismic analysis can help iden- tifying the observed modes. We computed mode stability analysis with the state-of-the-art linear nonradial nonadiabatic onedimensional oscillation code GYRE (v. 4.4, Townsend & Teitler 2013). GYRE accounts for the effect of rotation on pulsation eigensolutions by employing the traditional approximation of rotation (Eckart 1960), which is accomplished by ignoring the θcomponent of angular velocity in the governing oscillation equations (Lee & Saio 1986Townsend 2003). This is a valid approach, because Townsend (2005a,b) has already shown that the coupling of the κ-mechanism and the Coriolis force provides a sufficient restoring force for the excitation of gravito-inertial waves in rotating SPB stars. Two confirmed examples of gravitoinertial waves in B stars are HD 43317 (Pápics et al. 2012) andKIC 7760680 (Moravveji et al. 2016).
As an input to GYRE, we use a MESA model that is closest to the position of HD 261810 on the Kiel diagram (Fig. 15), and has just started burning hydrogen in the core. Given the projected rotational velocity υ sin i = 180±15 km s −1 , and the five regularly spaced periods, we attempt to reproduce the observed modes. To that end, we vary the mode spherical degrees ( , m) and the uniform rotation Ω rot . At the same time, we try to satisfy the following five constraints regarding the observed frequencies: (a) they must remain positive-valued in the co-rotating frame, (b) they must be predominantly excited, (c) they must be consecutive in radial order, as inferred from observations, (d) the equatorial rotation velocity must exceed the projected rotation velocity, i.e. v rot ≥ v sin i to ensure the projection factor of the inclination angle stays between zero and unity | sin i| ≤ 1, and (e) the observed g-modes are zonal or prograde m ≥ 0, given the negative slope of the period spacing (bottom right panel of Fig. 18). The latter constraint stems from the theoretical predictions of Townsend (2003) and Bouabid et al. (2013), K. Zwintz et al.: Young B stars in NGC 2264 Fig. 18. Comparison of the period spacing patterns of HD 48012 (yellow squares) and HD 261810 (green circles) to the seven SPB stars observed by the Kepler space mission and discussed in Pápics et al. (2016). . From flux cancellation of unresolved stellar disks, we know that the dipole ( = 1) and quadrupole ( = 2) modes have greater chance to be observed in white-light photometry than higher-degree modes do (Aerts et al. 2010). For these reasons, we restrict our further seismic analysis to dipole ( = 1) and quadrupole ( = 2) zonal and prograde g-modes. Townsend (2003) and Bouabid et al. (2013) have already demonstrated that for a fixed harmonic degree and uniform rotation velocity, the g-mode period spacing for tesseral prograde modes m > 0 has a steeper slope than that of zonal m = 0 modes. In other words, to reproduce the observed period spacing slope of HD 261810, we can either adopt much higher rotational velocity for zonal modes, or adopt a lower rotational velocity for tesseral modes. In the following, we adopt ( , m) = (1, 0), (1, +1), (2, 0), (2, +1) and (2, +2), and try to simultaneously fulfill the five constraints (a) to (e) given above.
For each choice of and m, we search for a rotation rate for which the five observed periods coincide with five consecutive g-modes. Our findings are summarized in Table 14. The first column gives the choice of spherical harmonics, and the second column gives the optimal rotation velocity v rot from a model that sufficiently matches the observed periods. The third column gives the corresponding rotation rate Ω rot /Ω K in percent. The Keplerian break up velocity is Ω K = GM /R 3 1/2 . The fourth column gives the ratio of the rotation velocity to the projected rotation velocity v rot /v sin i, where i is the inclination angle in degrees. The last column concludes if the corresponding harmonic degrees can provide a possible explanation to the observed modes. Figure 19 shows the evolution of predicted excited (blue) versus damped (grey) mode periods in the inertial (observer's) frame versus the equatorial rotation velocity on the lower abscissa. The upper abscissa shows Ω rot /Ω K . The horizontal red lines mark the five observed periods of HD 261810. The dipole prograde ( , m) = (1, 1) and quadrupole zonal ( , m) = (2, 0) modes are shown on the left (and right) panels. The best matching modes are shown with larger blue squares. Similar figures for other spherical degrees are not presented here for brevity.
The dipole zonal modes cannot explain the observed modes, unless the star rotates 85% its break up rotation frequency. At such rotation rates, B stars are expected to form a decretion disk, turn into B[e]-type stars, and exhibit emission spectral lines (e.g., Rivinius et al. 2013). However, this contradicts the observed spectrum of this star (Paper II), which does not show emission lines. Thus, we interpret this as an unlikely circumstance, and exclude dipole zonal modes. The quadrupole tesseral modes provide a possible match to the observed periods, but only for rotational velocities significantly smaller than the projected rotational velocity; as the latter contradicts with the constraint (d), we exclude quadrupole tesseral modes, too.
The dipole prograde and quadrupole zonal modes provide a degenerate solution to the observations. They both equally ful-  fill all modeling constraints (a) to (e). The shortest-period dipole prograde mode that fits the observations is found to be stable. This is not a major concern since a slight increase in metallicity can increase the height of the iron opacity bump, and destabilize the additional mode. For young B-stars, a slight metallicity excess is justifiable. Despite the fact that the two above models are degenerate, the valid range of the rotational velocity and the inclination angle are narrow. The two plausible rotational velocities, 236 or 267 km s −1 , are only ∼12% different. The two inclination angles of 41 or 48 degrees imply that HD 261810 is observed at an intermediate inclination angle which keeps both sets of modes visible.

Asteroseismic Inference for HD 48012
The star HD 48012 has the highest υ sin i of the sample with υ sin i= 225±25 km s −1 . It exhibits six significant pulsation frequencies, five of which form a period spacing series with a negative slope (Fig 17). In this regard, the same analysis that was already presented in the previous section can be repeated for this star. Table 15 and Fig. 20 present the modeling results. The dipole zonal modes cannot represent the observed periods because they exhibit no matching patterns even up to 80% break up velocity.
The same is true for the quadrupole tesseral modes. The dipole prograde modes can reproduce the observed periods ( Fig. 20 left panel), provided the star rotates uniformly at 45.75% break up velocity, or with an equatorial velocity of 258 km s −1 . In this case, the inclination angle is 61 degrees. Alternatively, the quadrupole zonal modes can also provide a fit to the observed periods ( Fig. 20 right panel), provided the star rotates uniformly at 74% break up velocity. This corresponds to an equatorial rotational velocity of 416 km s −1 with an inclination angle of 33 degrees. Similar to the previous star, the fast rotation rate deduced from quadrupole zonal modes is unlikely, because there is no signature of B[e] phenomenon in the spectrum of this star. Consequently, we reject this identification for the modes.
The validity domain of the traditional approximation, compared to the results of two-dimensional frequency computations is limited to ∼ 40% (Ballot et al. 2012;Prat et al. 2016;Ouazzani et al. 2016). The proposed solution for HD 48012 (Table 15) exceeds this validity domain. To a much lesser extent, the same is true for HD 261810. Consequently, fitting the observed periods of HD 48012 with quadrupole zonal modes -which requires Ω rot /Ω K = 74% -is unreliable. Based on this, the corresponding rotational velocity and inclination angles are subject to large model uncertainties.  Fig. 19 for HD 48012. Note that the projected rotation velocity is v sin i = 220 ± 20 km s −1 . The ( , m) = (2, 0) is rejected because it requires high rotation velocities, giving rise to B[e] phenomenon, which is not observed in spectroscopy.

Summary and conclusions
We have conducted a first comprehensive investigation of the variability of B stars in the young open cluster NGC 2264 using space photometry obtained by the MOST, CoRoT and Spitzer satellites. This will be complemented by a detailed study of the fundamental parameters, chemical abundances, presence of magnetic fields, spectral energy distributions and distances in the second part of this work (i.e., Paper II by Przybilla et al. 2017, in prep.). Our sample comprises in total eleven stars in the field of the cluster where one object, HD 261054, was identified as non-member with Be type variability and the remaining ten stars are likely cluster members according to their proper motions, radial velocities, distances and positions in the cluster colour-magnitude diagram. These ten members of NGC 2264 show variability due to SPB type pulsations (HD 47469, HD 48012, HD 261810, HD 261878), due to rotational modulation caused by spots on their surfaces (HD 47777, HD 47887, HD 261903, NGC 2264 137, and possibly HD 261938) or are connected to binarity (HD 47961 and possibly HD 261938) which will be discussed in detail in Paper II. The WEBDA database of open clusters 2 lists 32 B type stars in the field of NGC 2264; twenty-six of these are cluster members according to their proper motions (Vasilevskis et al. 1965), their radial velocities (Fossati et al. 2014) or identified to be likely members in our study. The ten variable stars analyzed here are included in this list, and therefore correspond to a percentage of 38% variable and 15% pulsating stars among the spectral type B objects in NGC 2264.
Originally, we intended to search for pre-MS SPB type stars in NGC 2264. Using our theoretical models computed with the MESA stellar structure and evolution code (Paxton et al. 2011(Paxton et al. , 2013(Paxton et al. , 2015 in combination with the stars' fundamental parameters derived from high-resolution spectroscopy (Paper II), we find them to be located near the ZAMS in the HR-diagram (Fig.  15). Although their ages computed from these models lie between only ∼ 1 and 6 million years, all stars seem to be early ZAMS objects that are not in their pre-MS phase any more. This illustrates how fast early stellar evolution progresses for (late) B type stars and that bona-fide pre-MS SPB stars are statistically hard to discover.
Two of the four SPB type pulsators show regularly spaced gmode period series which we used for a first asteroseismic analysis. For HD 261810 two possible mode identifications are either ( ,m) = (1,+1) or ( ,m) = (2,0). For HD 48012, the only plausible solution is ( ,m) = (1,+1). We also compare the period spacings of our two stars (HD 48012 and HD 261810) to those of seven SPB stars observed by the Kepler main mission (Borucki et al. 2010) and analyzed by Pápics et al. (2016). With projected rotational velocities, υ sin i, of 225±25 and 180±15km s −1 for HD 48012 and HD 261810, respectively, they are among the fastest rotators in that sequence, and, hence, show the shortest periods and period spacings.
Both HD 261810 and HD 48012 are fast rotators, with their Ω rot exceeding ∼ 40% of their breakup rotation frequency (see Tables 14 and 15). For such high rotation rates, the validity of traditional approximation of rotation is questionable (Mirouh et al. 2016;Prat et al. 2016). However, those formalisms were developed based on polytropic models, and the effect of chemical stratification on the pulsational behavior was ignored. As a result, our findings regarding the rotational velocity and inferred inclination angles might change if one goes beyond the traditional approximation of rotation, or uses 2D pulsation codes (Ouazzani et al. 2016). Consequently, these two stars are ideal test cases for more sophisticated models of rotating pulsating B stars, if higher-precision data with longer time baseline is collected in the future with e.g. the PLATO 2.0 mission.
The currently available photometric time series for the young B stars have time bases of less than 40 days. Hence, the noise level is too high to detect longer period spacing series for HD 48012 and HD 261810 and to discover possible period spacing patterns for HD 47469 and HD 261878. With longer photometric time series, e.g., from future space missions like TESS (Ricker et al. 2015) and PLATO (Rauer et al. 2014), we will also be able to conduct a more in-depth asteroseismic analysis of the internal structures of young ZAMS B type pulsators and compare them to their more evolved counterparts. With such data we will also be able to continue our search for bona fide pre-MS SPB stars.