There are different classes of pulsating stars in the HR diagram. While many of those classes are undisputed some remain a mystery such as the objects historically called Maia variables. Whereas the presence of such a class was suggested seven decades ago no pulsational driving mechanism is known that could excite shortperiod oscillations in these late B to early Atype stars. Alternative hypotheses that would render the reports of variability of those stars erroneous have been proposed such as incorrect effective temperatures binarity or rapid rotation but no certain conclusions have been reached yet. Therefore the existence of these variables as a homogeneous class of pulsating star is still under discussion. Meanwhile many new candidates of these variables have been claimed especially by using photometric observations of space telescopes. In this study we examined 31 objects that are alleged members of this hypothetical group and carried out detailed spectroscopic and photometric analyses to test the proposed hypotheses for their cause of variability. The Ttextrmeff  log g  v sin i  and chemical abundances of the targets were determined and the TESS photometric data were examined. As a result we found that most of these targets are located inside the  Scuti  Cephei or SPB star instability strips a few show evidence for binarity and others for rapid rotation. We give arguments that none of the apparently rapid pulsations in our targets is caused by a star outside any known instability strip. By extrapolation we argue that most stars proposed as pulsators outside wellestablished instability domains are misclassified. Hence there is no sufficient evidence justifying the existence of a class of pulsating stars formerly known as the Maia variables.
['Stars: oscillations', 'stars: variables: general', 'techniques: photometric', 'techniques: spectroscopic', 'Astrophysics - Solar and Stellar Astrophysics']
On the existence of Maia variables
{"Onthe Existence of 'Maia variables'": 'F. Kahraman Aliçavu/uni015F, 1,2 G. Handler, 3 S. Chowdhury, 3 E. Niemczura, 4 R. Jayaraman, 5 P. De Cat, 6 D. Ozuyar, 7 and F. Aliçavu/uni015F 1,2  \n- 1 Çanakkale Onsekiz Mart University, Science Faculty, Physics Department, 17100, Canakkale, Türkiye\n- 2 Çanakkale Onsekiz Mart University, Astrophysics Research Center and Ulupınar Observatory, TR-17100, Çanakkale, Türkiye\n- 3 Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warsaw, Poland\n- 4 University of Wrocław, Astronomical Institute, Kopernika 11, 51-622, Wrocław, Poland\n- 5 Department of Physics, and Kavli Institute for Astrophysics and Space Reserch, M.I.T., Cambridge, MA 02139, USA\n- 6 Royal Observatory of Belgium, Ringlaan 3, B-1180 Brussel, Belgium\n- 7 Ankara University, Faculty of Science, Dept. of Astronomy and Space Sciences, 06100, Tandogan, Ankara, Türkiye\n- Author for correspondence: F. Kahraman Aliçavu/uni015F & G. Handler, Email: filizkahraman01@gmail.com, gerald@camk.edu.pl .  \n(Received dd Mmm YYYY; revised dd Mmm YYYY; accepted dd Mmm YYYY; first published online dd Mmm YYYY)', 'Abstract': 'There are different classes of pulsating stars in the H-R diagram. While many of those classes are undisputed, some remain a mystery such as the objects historically called "Maia variables". Whereas the presence of such a class was suggested seven decades ago, no pulsational driving mechanism is known that could excite short-period oscillations in these late B to early A-type stars. Alternative hypotheses that would render the reports of variability of those stars erroneous have been proposed such as incorrect effective temperatures, binarity or rapid rotation, but no certain conclusions have been reached yet. Therefore the existence of these variables as a homogeneous class of pulsating star is still under discussion. Meanwhile, many new candidates of these variables have been claimed especially by using photometric observations of space telescopes. In this study, we examined 31 objects that are alleged members of this hypothetical group and carried out detailed spectroscopic and photometric analyses to test the proposed hypotheses for their cause of variability. The T eff , log g , v sin i , and chemical abundances of the targets were determined and the TESS photometric data were examined. As a result, we found that most of these targets are located inside the δ Scuti, β Cephei, or SPB star instability strips, a few show evidence for binarity and others for rapid rotation. We give arguments that none of the apparently rapid pulsations in our targets is caused by a star outside any known instability strip. By extrapolation, we argue that most stars proposed as pulsators outside well-established instability domains are misclassified. Hence there is no sufficient evidence justifying the existence of a class of pulsating stars formerly known as the "Maia variables".  \nKeywords: stars: oscillations, stars: variables: general, techniques: photometric, techniques: spectroscopic', '1. INTRODUCTION': 'Perhaps the most efficient and accurate way to investigate stellar interior structure is asteroseismology which uses the eigenfrequencies of stellar oscillations to probe the interiors of stars. Using asteroseismology phenomena such as internal rotation, angular momentum, and convective core overshooting have been successfully examined (e.g., Christophe et al., 2018; Townsend et al., 2018; Walczak et al., 2019; Kurtz, 2022).  \nThe loci of several different classes of pulsating stars can be delineated in the Hertzsprung-Russell (H-R) diagram. These loci are called "instability strips" and comprise groups of different evolutionary statuses. Among those on the main sequence, six types of pulsating stars are commonly known. Moving from hotter to cooler temperatures, these variables are the β Cephei stars, the slowly pulsating B-type (SPB) stars, the δ Scuti stars, the rapidly oscillating Ap (roAp) stars, the γ Doradus stars and solar-like oscillators. The β Cephei stars are objects of spectral types B0 - B3 showing low-order pressure (p) and gravity (g) modes with frequencies > ∼ 3.5 d -1 (Stankov & Handler, 2005). On the other hand, the SPB variables are B3 - B9 stars which exhibit high-order g-mode oscillations with frequencies < ∼ 3.5 d -1 (Aerts et al., 2010). The δ Scuti stars exhibit highfrequencies ( > ∼ 5d -1 ) similar to β Cephei variables, however, their spectral types range from A2 to F5 (Chang et al., 2013;  \nUytterhoeven et al., 2011). The roAp stars comprise similar spectral types, but are chemically peculiar due to the presence of strong magnetic fields and pulsate in high-order pressure modes with even higher frequencies ( > ∼ 50d -1 ) (Holdsworth et al., 2024). The γ Doradus stars are low-frequency ( < ∼ 5 d -1 ) g mode pulsators with spectral types of F2 - F6 (Uytterhoeven et al., 2011). Solar-like oscillators are of even later spectral type on the main sequence and show stochastic high-order p mode oscillations with frequencies around 300 d -1 (e.g., Jackiewicz, 2021). All these pulsators have their own instability strips where theoretical models can explain the observed types of oscillations. Consequently, almost the whole main sequence is populated with pulsating stars with spectral types ranging from early B to late G, but with a marked gap between spectral types B9 and A2. The reason for this gap is that both B-type variable star classes show pulsations triggered by the κ mechanism operating in the partial ionization zone of iron group elements (Dziembowski et al., 1993), but the δ Scuti stars are triggered by the κ and/or the turbulent mechanism acting in the partial ionization zones of Helium (Baker & Kippenhahn, 1962; Antoci et al., 2014) and hydrogen (e.g., Cunha et al., 2013). In the gap between those groups of pulsators at late B/early A spectral type, none of the currently known pulsation mechanisms are able to drive pulsations in models, consistent  \n<!-- image -->  \nwith observations. In particular, no pulsations with intrinsic frequencies in excess of 3 d -1 should exist in stars with spectral types between B3 and A2.  \nHowever, a group of high-frequency pulsating stars in this gap and in the cool part of the SPB instability strip has been proposed by several authors (e.g., Struve, 1955; Lehmann et al., 1995). These variables seem to have high frequencies like β Cephei and δ Scuti stars but they are too cool to be a β Cephei and too hot to be a δ Scuti pulsator. Hence they appear to populate the gap where pressure mode pulsations should not exist. These variables have historically been known as the "Maia stars" following the apparent prototype of these variables, Maia in the Pleiades cluster. Maia was first thought to be a high-frequency pulsator (Struve, 1955), but this claim was refuted several times (e.g., McNamara, 1985, and references therein) and in the end it was shown that Maia is actually a rotational variable (HgMn star, White et al., 2017) without any detectable pulsational variability. However, recent studies, especially those based on space-based data, revived interest in those objects (Balona et al., 2016; Balona & Ozuyar, 2020; Balona, 2023).  \nFurthermore, examples of objects with slower apparently pulsational variations have been reported in the gap between the main sequence instability strips (e.g., Aerts & Kolenberg, 2005; Degroote et al., 2009; Mowlavi et al., 2013; Sharma et al., 2022), which are not to be confused with the short period variables mentioned above. Nevertheless, the hypotheses to explain all these variable stars are essentially the same. We call the first one the "rapid rotation hypothesis". This idea would explain the apparent pulsators as rapidly rotating SPB or β Cephei stars seen equator-on. Because of this fast rotation, the star would have an ellipsoidal-like shape, and its effective temperature ( T eff ) would be lower near the equator and higher near the poles. This would have several consequences. First, the observationally determined T eff would be different from the overall value due to projection effects, modifying the star\'s position in the HR Diagram. Second, the rapid rotation would shift the observed pulsational mode frequencies (Townsend, 2005; Salmon et al., 2014; Szewczuk & Daszyńska-Daszkiewicz, 2017; Mowlavi et al., 2016) to higher values. Third, Rossby waves (r modes) could be excited (Szewczuk & Daszyńska-Daszkiewicz, 2017; Saio et al., 2018). Fourth, the non-uniform temperature distribution in the star could facilitate pulsational excitation near the stellar poles or the equator (cf. Kurtz et al., 2023). Fifth, rapid rotation would extend the SPB instability strip, also towards lower temperatures (Townsend, 2005). Another hypothesis to explain pulsation outside of any known instability strip is associated with binarity (hereinafter the "binary hypothesis"). The pulsating stars could be components of binary systems consisting of a non-pulsating late B-type star, which dominates the spectrum, and a δ Scuti pulsator. Given that most B-type stars have companions, and often more than one (Moe &Di Stefano, 2017), such a scenario is not unlikely. Therefore, apparent p-mode pulsations would be associated with a star too hot to have them excited. The third hypothesis involves simple observational uncertainties that lead to incorrect T eff  \nvalues (hereinafter the "wrong T eff hypothesis"). Most of the variables under question were observed by the Transiting Exoplanet Survey Satellite ( TESS , Ricker et al., 2014). The T eff of these B type systems given in the TESS input catalogue (Brown et al., 2011; Stassun et al., 2019) may not always be reliable. These T eff parameters are in most cases photometric estimates based on a system lacking a U-band filter, and that carry some uncertainties regarding dereddening for stars located close to the Galactic plane (Stassun et al., 2019). Therefore the T eff values for hot stars from the TIC catalogue have to be treated with caution. Finally, regarding TESS data, a "contamination hypothesis" must be considered. Due to the large (21") size of TESS \'s CCD pixels on the sky, the light of nearby stars may fall into the photometric aperture, hence, for the purpose of example, an unrelated δ Scuti star may introduce short-period pulsational variations into the light curve of a star outside of any pulsational instability strip.  \nBalona et al. (2016); Balona & Ozuyar (2020) tested some of these hypotheses, but without an adequate amount of good quality spectroscopic data. In a recent study Balona (2023) even suggested that Maia variables are an extension of δ Scuti stars and there is an interplay between available pulsational driving mechanisms, but this conclusion was partly based on inhomogeneous literature data. To improve on this situation one needs to study these variable stars with medium to highresolution spectroscopic data to determine some important parameters (e.g. T eff , v sin i ) accurately.  \nIn a recent study, Kurtz et al. (2023) analyzed an apparently genuine, very rapidly rotating p-mode pulsator of spectral type A0 (HD42477), which is located between the β Cephei and δ Scuti instability strips, by taking into account the incorrect T eff , contamination, and binarity hypotheses. They ruled out or placed strong constraints against all these hypotheses and concluded that p-mode pulsations can be present in the mentioned region of the H-R diagram. Their excitation would be explained by coupling with several g and r modes or by the classical κ -γ mechanism operating in the He/I.sc/I.sc ionization zone in the stellar equatorial regions, hence is understood based on current astrophysical knowledge.  \nRecapitulating, there is one well-documented case of pressure mode pulsation between the instability strips of the SPB and δ Scuti stars. In addition, there are claims of hundreds of such stars based on space photometry data, but without proper assessments of whether the fast rotation, the binary, the contamination, or the wrong T eff hypotheses would apply. Given the sheer number of these objects, is it really possible that most of them were claimed by misfortune?  \nIn this study, we present detailed analyses of some candidate variables in the gap between the main sequence instability strips to test the proposed hypotheses and to establish whether a class of "Maia variables" a is justified. We use TESS photometric data and medium and/or high-resolution spectra and the paper is organized as follows. In Sect. 2, we introduce the target selection and observations. In Sect. 3, the spectroscopic  \na Given that Maia itself is not a pulsating star, a different name for this prospective class would be required.  \nanalyses including radial velocity measurements, determination of atmospheric parameters, and chemical abundances are presented. A contamination analysis on TESS data is introduced in Sect. 4. Discussions and conclusions are given in Sect. 5, and 6, respectively.', '2. Target selection and observations': 'Our variable star candidates were selected from the list of Balona & Ozuyar (2020). These authors defined priority classes considering the methods used to estimate T eff values. The highest priority was given for spectroscopic T eff while the lowest priority was defined for T eff estimated from the spectral types. In between, there are T eff estimations from some photometric colours. In this study, we preferentially chose stars of lower priority, as this adds more value to our spectral analyses for accurate T eff determinations.  \nFigure 1. The initial positions of the selected candidates in the H-R diagram as derived by Balona & Ozuyar (2020). The theoretical instability strips of the β Cephei and SPB stars were taken from Pamyatnykh (1999). The domain of the δ Scuti stars was taken from Murphy et al. (2019).  \n<!-- image -->  \nWe reanalyzed the TESS Simple Aperture Photometry (SAP) 2-min cadence light curves of our selected candidates to confirm the presence of pressure mode oscillations in these variables. We made use of the Period04 (Lenz & Breger, 2005) package, which utilizes simultaneous multi-frequency sinewave fitting with input parameter optimization. After removing the sine-wave fits from the original data, the residuals were then re-analyzed in the same manner to detect further periodicities. The signals were extracted until we reached a Signal-to-Noise ratio (SNR) of 4.5 (Baran & Koen, 2021). In addition, the initial positions of the selected candidates in the H-R diagram were checked using the luminosity ( L ) values calculated from the Gaia parallaxes (Gaia Collaboration, 2022), interstellar reddening, E(B-V), calculated from the interstellar extinction map (Amôres & Lépine, 2005), the bolometric corrections (Flower, 1996) and the T eff values given by Balona & Ozuyar (2020) (see, Fig. 1). As a result, 31 candidates were selected for this study; they are listed in Table 1.  \nFor these candidates, we carried out spectroscopic surveys. Some spectra were taken from the Fibre-fed Échelle Spectrograph (FIES, Telting et al., 2014) that is a cross-dispersed high-resolution échelle spectrograph attached to the 2.56-m  \nTable 1. The list of the targets. The number of the taken spectra from FIES and HERMES and their signal-to-noise ratios (SNR) are given in the last two columns, respectively. These spectra were taken between 2020 and 2021.  \nNordic Optical Telescope of the Roque de los Muchachos Observatory (ORM, La Palma, Spain). It has three different resolution options; low ( R = 25000), medium ( R = 46000) and high ( R = 67000). Considering the T eff values and the V-band magnitudes of the target stars, the low-resolution mode was used during the observations to achieve a higher signal-to-noise ratio (SNR). In addition, this resolving power is good enough for spectral analysis of systems having medium to high projected rotational velocity ( v sin i ≳ 30 kms -1 ). We also gathered some spectra from the High Efficiency and Resolution Mercator Échelle Spectrograph (HERMES, Raskin et al., 2011) which is a high-resolution ( R = 85000) fibre-fed échelle spectrograph attached to the 1.2-m Mercator telescope at the Roque de los Muchachos Observatory (ORM, La Palma, Spain). Some stars have spectra from both spectrometers, while some have only from one. For each target, we tried to obtain at least two spectra taken at different times to check for possible binary-induced radial velocity shifts. The information about the spectroscopic surveys is given in Table 1.', '3.1 Radial velocity measurements': 'Todiscover possible variability of the targets induced by orbital motion in a binary system, we used spectra taken at different epochs. There are at least a few hours (mostly a day) or a maximum of around two months between two consecutive spectra of the same star for the same spectrometer.  \nThe radial velocities of the systems were measured using the cross-correlation technique. The IRAF b (Tody, 1986) FXCOR task was used in the investigation. As templates, some synthetic spectra were generated with different T eff values of 8000 K, 10000 K and 15000 K.These atmospheric models were generated using local thermodynamic equilibrium (LTE) ATLAS9 models (Kurucz, 1993) and the SYNTHE code (Kurucz & Avrett, 1981). Each spectrum of the individual stars was analyzed. However, for most of the stars, no trace of another component and no significant difference between the radial velocity measurements were found within error bars. There are two systems, TIC 202567458 and TIC 251196433, whose line profile shapes differ from those of single stars. However, because of high rotation and low SNR values, we could not detect double peaks in cross-correlation analyses of these systems. The line profile shapes of the two atypical systems could be caused by pulsations as well. We classified these systems as suspected binaries. Table A1 gives a list of the radial velocity measurements.', '3.2 Atmospheric parameters': 'For the analysis of stars having T eff lower than ∼ 15000K LTE model atmospheres are convenient (Przybilla et al., 2011). Therefore at the beginning of our study, to estimate the T eff , surface gravity (log g ), and v sin i parameters we used LTE ATLAS9models (Kurucz, 1993) and SYNTHE code (Kurucz & Avrett, 1981) to generate theoretical spectra. The v sin i values  \nwere determined by the profile fitting method using unblended metal or helium (He) lines (Gray, 2008). The T eff and log g parameters were derived utilizing the hydrogen lines. During the analysis, the minimization method was used to obtain the atmospheric parameters (Catanzaro et al., 2004). Since most stars in our list seem to have moderate to high rotational velocity, rather than taking into account the resolving power of the spectrometers we preferred the spectra which have the highest SNR, if a star was observed by both spectrometers.  \nAccording to initially derived T eff parameters from the LTE models, some stars were found hotter than 15000 K. For these stars, the same analysis was carried out by using the nonLTE Tlusty BSTAR2006 grids (Lanz & Hubeny, 2007). These non-LTE models were generated using the Synspec (Hubeny & Lanz, 2011) program. As a result of LTE and non-LTE analyses, the estimated parameters are listed in Table 2. The uncertainties of the determined parameters were estimated by the 1σ tolerance in the goodness-of-fit parameter. The consistency between the theoretical and observed spectra is shown in Fig. 2 for some samples.', '3.3 Chemical abundances': 'After we determined the parameters given in Table 2, in the next step we singled out stars that are suitable for chemical abundance analysis. The stars having v sin i lower than ∼ 200kms -1 , and spectra with a good SNR were chosen for that purpose. During this examination, previously derived atmospheric parameters were taken as input and they were updated utilizing the excitation and ionization equilibrium as described by Kahraman Aliçavuş et al. (2016).  \nIn the analysis, the Kurucz line list c was used for the line identification and we used the LTE ATLAS9 models. The spectral synthesis method was considered in the analysis. As a result of the investigation, updated atmospheric parameters and chemical abundances were derived. The list of the updated parameters is given in Table 2, while the element abundances are presented in Table A2. The distribution of the final T eff , log g , and v sin i parameters for all targets are shown in Fig. 3. In addition, the Fe abundance distribution of stars is shown in Fig. 4.', '4. Examination of contamination': "To complement our spectroscopic analyses, we used data from TESS . However, the TESS pixel scale is 21', which may induce contamination of the photometric signal from stars located nearby in the sky. To characterize the extent of the contamination in the TESS data of our candidates, we worked directly with target pixel files (TPFs).  \nFor any star in which the contamination percentage was greater than 1%, we queried Gaia DR3 for all the stars brighter than a G mag of 17 within a 200' radius of the TESS source. Any targets that had a star within the SPOC photometric aperture that had ∆ G mag < 5 were thoroughly investigated through the use of custom apertures for photometry. We  \nTable 2. Atmospheric parameters of the stars determined from the hydrogen Balmer lines and the updated parameters. An asterisk represents a fixed log g value for TIC 352781227 that we could not determine because of rapid rotation and low SNR. Since some stars are not suitable for further analysis because of their large rotational velocities and low SNR values of their spectra, only some of the parameters could be updated (see, Sect 3.3). E ( B -V ) was determined from the sodium lines and the final L parameters are given. The Gaia RUWE parameters are also listed in the last column.  \nFigure 2. Upper panels: Theoretical spectral fits (red dashed lines) to the observed spectra of target stars. Lower panel: Residuals.  \n<!-- image -->  \nFigure 3. Distributions of the T eff , log g and v sin i parameters.  \n<!-- image -->  \nFigure 4. The di/fferences between the obtained and solar (Asplund et al., 2009) Fe abundances for each star. The x-axis represents the TIC numbers.  \n<!-- image -->  \nchecked to see whether we recover the same periodicities for the Maia stars from the SPOC light curve and the light curve calculated using the custom aperture, which will allow us to have confidence in our identification of the pulsation modes.  \nAs part of this analysis, 22 of the 31 stars were found to have contamination ratios greater than 1%. Several individual stars were selected for individual discussion in Appendix 1 if one of the following criteria was true: (i) There are changes in the number or location of significant frequencies when the periodogram of the custom light curve is calculated, (ii) there exist potentially contaminating stars that are within the SPOC photometric aperture and have a ∆ G mag ≤ 3 compared to the target of interest, or (iii) there is a star with G mag ≤ 14 whose signal falls into the same pixel as the star of interest. Some of the custom-generated light curves were corrected with cotrending basis vectors to remove long-term trends. The stars that are not discussed in Appendix 1 can be considered to have intrinsic variability that is appropriately reflected in the SPOC-generated light curves.  \nNote that there is a significant tradeoff between the custom light curves and the light curves generated by the SPOC pipeline-while the former will isolate signals to the star of interest, the white noise level will be somewhat higher because the detrending and photometry for this light curve are more rudimentary than the tools SPOC uses. This may mean that some of the lower S/N peaks in the periodogram for the SPOC light curve are washed out by the higher white noise level. Cases where the white noise level of the custom light curve's periodogram is comparable to the height of the peak in the periodogram of the SPOC light curve, are noted below.  \nAdditionally, in some cases, we used the TESS -Localize package (Higgins & Bell, 2022) to identify the most likely source of some of the frequency peaks in the periodogram. All calls to TESS -Localize were used with the principal\\_components parameter set to 2, in order to de-trend the light curve generated from the TPF.  \nIn Fig. 5, a plot of the TESS pixel files of TIC 64774351, TIC352781227, TIC467219052 and TIC450085523 is presented. Custom apertures were created containing just the target stars; for some sources, multiple apertures were created to verify the results, as discussed in Appendix 1. The periodograms for the custom apertures are shown in blue in the  \nsame figure. According to periodograms and the discussion given in Appendix 1, for TIC 352781227, TIC 467219052, and TIC450085523 we do not find any significant contamination that would cause frequencies higher than 3 d -1 . However, for TIC64774351 the result of contamination analysis could not give a reliable result for high-frequencies seen in the original periodogram.  \nIn Fig. 6, the heavily contaminated stars are shown. For TIC255996252 according to our investigation, we found that the target of interest is the original source of high-frequencies (for detail see Appendix 1). However, for other systems, we were unable to distinguish whether or not the high-frequency peaks come from our targets.", '5. Discussion': 'We now proceed to the overall discussion of the variability of our sample of target stars. Each individual star is considered in more detail in Appendix 2. We present the positions of our targets in the theoretical H-R Diagram in Fig. 7. The frequency spectra of all systems are shown in Fig A1 in order from the hottest star to the coldest star.  \nThe positions of the candidates in this diagram take into account the final atmospheric parameters, the structure of frequency spectra, and the v sin i values. To plot the final H-R diagram the L parameters were recalculated taking into account the final T eff values. Additionally, the E(B-V) parameters were re-determined from the interstellar sodium lines as described in the study of Kahraman Aliçavuş et al. (2016) (cf. Table 2).', '5.1 Incorrect T eff hypothesis': 'In stark contrast to Fig. 1 most of the target stars now lie in known pulsational instability regions, suggesting that most literature T eff values used to select them are inaccurate. We compared our spectral T eff values with the T eff of the TIC and Balona & Ozuyar (2020) in Fig. 8. As can be seen from the figure there are significant differences between the T eff values. This result is an expected conclusion for the TIC T eff values for hot stars because of the uncertain reddening corrections for stars in the Galactic plane. A similar explanation may apply to the T eff values used by Balona & Ozuyar (2020), as those were derived mostly from photometric colours. According to this comparison, we conclude that for most targets previously derived T eff values show significant discrepancies from the spectral ones. Even so, in Fig. 7 there are still some stars located outside of any known pulsational instability strip.  \nAt this point it is important to stress that the instability strip boundaries depicted in Figs. 1 and 7 are not immutable. The theoretical boundaries ( β Cep stars, SPB stars) depend on the input physics used, as discussed by (e.g., Aerts & Tkachenko, 2023), and the observational boundaries for the δ Scuti stars indicate where pulsations are less likely (Murphy et al., 2019). Nevertheless, these Figures give an idea in which stars pulsations could be difficult to explain.  \nFigure 5. A plot of the TESS target pixel files of 4 of the 7 stars discussed in Appendix 1. These stars were the ones for which custom light curves were made. The top panel shows the target pixel files, with the star of interest shown in green (and indicated in the title) and the potentially contaminating stars shown in tan. The middle panel shows the light curve created using a custom aperture (shown in blue in the top panel); the bottom panel shows the light curve generated via the SPOC pipeline. Note that the y-axis scales di/ffer between the middle and bottom panels as a result of the fact that the custom aperture imperfectly captures the flux from the star of interest. This, however, provides the advantage of allowing us to highlight the peak locations.  \n<!-- image -->  \nFigure 6. Eachpanelcontainsthetargetpixelfileforastarthatis so heavily contaminated that we are unable to determine where the true source of the frequency peaks in the periodogram lies. As in Figure 5, the green points represent the source of interest, while the tan points represent potentially contaminating stars.  \n<!-- image -->  \nFigure 7. The final positions of the selected Maia candidates in the H-R diagram. The theoretical instability strips of β Cephei , SPB and δ Scuti stars were taken the same as in Fig. 1. A detailed explanation of the brief expressions of frequencies in the figure can be found in Sect. 5. Some targets are not included in this figure, since they were excluded from the Maia candidate list according to the results of our contamination analysis.  \n<!-- image -->  \nFigure8. Comparisonofthe T eff of the TIC (Stassun et al., 2019) (upper panel) and Balona & Ozuyar (2020) (lower panel) with the spectral T eff determined in this study. The x-axis shows the spectral T eff values while the y-axis illustrates the T eff di/fferences between the spectral and TIC (upper panel) and Balona & Ozuyar (2020) (lower panel).  \n<!-- image -->', '5.2 Binary hypothesis': 'According to the binary hypothesis, it was suggested that late B/early A type pressure mode pulsators could be binary systems including a B-type star with a δ Scuti companion. Hence, the system seems to be an oscillating star but it is too cool to be a β Cephei and too hot to be a δ Scuti star. As a result of our investigation of the spectral line profiles and radial velocities of the stars in our sample, we did not reveal any obvious binary systems. Only two systems were classified as suspected binary systems because of their line profile changes, which however could be well caused by pulsations. Additionally, our method is not sensitive to binary systems with long orbital periods (more than a few months) given the temporal sampling of our spectra. Also, many stars in our sample show high rotational velocity (>100 km s -1 ) and some spectra have low SNR, which affects the radial velocity measurements and the detection of a possible binary component in the systems. Consequently, our spectroscopic investigation of possible binarity remained mostly inconclusive.  \nAnother possibility to detect binaries would be to search for phase modulations of the pulsational signals (Murphy et al., 2014). Unfortunately, we were not able to find any evidence thereof, presumably owing to the small amplitudes of the pulsations detected and the relatively short duration of the individual TESS light curves.  \nWe also examined the Gaia Renormalized Unit Weight Error (RUWE Lindegren et al. 2021). RUWE is a parameter that shows the degeneration in the astrometric parameters determined with Gaia (Gaia Collaboration et al., 2021). High values of RUWE (higher than 1.4) indicate that the system of interest could be a binary (Lindegren et al., 2018, 2021). Stassun & Torres (2021) also point out that any RUWE value above 1.0 may already be indicative of binarity. Therefore, we collected the RUWE parameters of our targets from Gaia DR3 (Gaia Collaboration et al., 2021) which are given in Table 2. As can be seen from the table, there are nine systems that have RUWE>1.4. The largest values, 11.59 and 20.59, relate to TIC 36557487, and TIC 367351021, respectively. One suspected binary system from our radial velocity analysis (TIC 251196433) has a RUWE value of 2.01. The systems with high RUWE values (> 1.4) are likely astrometric binaries, even if none of our targets is in the list of non-single stars according to Gaia DR3. Catalogues of double and/or multiple stars were also searched (e.g. Balega et al., 2012; Eggleton & Tokovinin, 2008; Fabricius et al., 2002; Mason et al., 2001; 199, 1997; Dommanget & Nys, 1995). Indeed some of our targets were found in these catalogs and are classified as double or multiple systems. TIC 36557487, TIC 64774351, and TIC354793407 are members of visual double systems.', '5.3 Rapid rotation hypothesis': "As discussed in the introduction, rapid rotation could a) modify a star's position in the H-R Diagram, b) shift the observed pulsation frequencies, c) excite Rossby waves, and d) provide means of κ -γ -driven excitation of p and g modes. To evaluate our candidate stars with respect to rotation we  \nfirst took into account the stars which show frequencies between 1-10 d -1 . There are four stars that show frequencies in that range (TIC 251196433, TIC 352781227, TIC 377099704, TIC467219052) that are shown with black dots in Fig. 7 and two of them exhibit high rotation; TIC 251196433 and TIC352781227. TIC251196433 is positioned slightly below the instability strip, situated between SPB and β Cephei stars. This fast-rotating star ( ≳ 255kms -1 ) could either be a hybrid pulsator or an SPB star displaying shifted pulsation frequencies. TIC352781227 is also a fast-rotating star ( ≳ 270kms -1 ) and placed between the δ Scuti and SPB instability strip. For this star, high rotational velocity could cause low frequencies to be shifted towards high frequencies and since we likely see this star close to equator-on, our T eff determination would be too low. Therefore, this variable is likely an SPB star. A more detailed discussion on this object is also given in Appendix 2.  \nThe other stars showing frequencies between 1-10 d -1 have low to moderate v sin i . One (TIC467219052) is simply a somewhat evolved δ Scuti star, while the other is located in the SPB domain, very close to the β Cephei territory (TIC377099704). These stars are hence well placed in known instability strips. Both show both low and high frequencies, these stars could be SPB β Cephei and a δ Scuti γ Doradus hybrids.  \nThere is a group of stars that are distinguished by two close frequencies dominant in the Fourier spectra (TIC 26369887, TIC36557487, TIC312637783, TIC 374038418, TIC 434340780). As depicted in Fig. 7 most of these stars exhibit rapid rotation with some falling between the SPB and δ Scuti domains. Because of high rotational velocity, the rapidly rotating ones could appear cooler than their actual temperature, leading to a shift in their pulsation frequencies towards higher values. We therefore classify them as cool SPB stars. Further discussions on these targets are provided in Appendix 2.  \nWhen examining other stars exhibiting oscillations in a wide frequency range (up to 50 d -1 ), we noticed that most of them are located inside the δ Scuti instability strip. However, six of them are placed beyond the hot border of δ Scuti instability strip. Among these, TIC 20534584 lies between the δ Scuti and SPB instability strips. Notably, it possesses a moderate to high v sin i (150 km s -1 ), and its relatively high RUWE parameter (2.26) suggests a potential binary nature with a hotter companion, placing it beyond the hot border of δ Scuti instability strip. There are also four stars (TIC 83803744, TIC 192780911, TIC236785664, TIC354793407) placed inside the SPB domain and showing frequencies up to 40 d -1 . While two of them have low v sin i values, TIC 83803744 exhibits a low RUWE parameter, and TIC 354793407 has a RUWE of 1.81. Regarding TIC 83803744 its high frequencies lack a clear explanation. However, it is classified as HgMn star (see Appendix 2). As HgMn stars are often found in binaries (e.g., Gerbaldi et al. 1985), which suggests a binary nature. The remaining two systems (TIC 192780911 and TIC 236785664) exhibit high rotational velocities with high values of the RUWE parameter indicative of a possible binary nature. TIC 192780911 displays a weak signal at around 3.9 -1 , possibly indicating an SPBtype pulsation modified by the rapid rotation. TIC 236785664  \nshows both low and high frequencies typically of δ Scuti stars, yet it's T eff exceeds the typical range for δ Scuti stars. Finally, TIC202567458, a rapid rotator, manifests a pulsation spectrum comprising apparent g-mode pulsations, their linear combinations, and a stochastic component implying the presence of Internal Gravity Waves (IGWs, e.g., see Ratnasingam et al. 2023; Bowman et al. 2019). The higher frequency signals are compatible with combination frequencies.", '5.4 Variability classification of the target stars': 'We summarize the discussions above and the results on the individual stars presented in Appendix 2 in Table 3, respectively. Out of the 31 stars considered in this paper and previously claimed as variables outside of any known instability strip, the largest group (14 objects) are normal δ Scuti stars that simply had overestimated T eff . Six stars are gravity mode pulsators with combination frequencies that protrude into the pressure mode domain, but these higher frequencies are not due to independent pulsation modes (cf. Kurtz et al., 2015). Among the hottest stars, we found one star with gravity modes and Internal Gravity waves, two β Cephei pulsators, and two β Cephei/SPB "hybrids".  \nRegarding the remaining six objects, the TESS light curve of one target was revealed to be contaminated by a background ellipsoidal variable. Two stars were identified as SPB stars with a very fast rotation that presumably modifies the pulsation frequencies in the observer\'s reference frame to values that are comparable with those of pressure modes. That leaves three stars, two of which are in the SPB star instability strip and one is in between the SPB and the δ Scuti star domains. Two of these stars are rotational variables with additional frequencies in the δ Scuti range. For both of those, there is some indirect evidence of binarity, one TIC 354793407 = HD 21699 has a RUWE of 1.81, and the other TIC 83803744 = HD 1009 is a HgMn star which in most cases are binaries. Consequently, the short-period signals in these star\'s light curves might well arise from a δ Scuti companion. The last object is a rapidly rotating Be star with a somewhat high RUWE and several frequencies typical for δ Scuti stars. This latter object, TIC 236785664 = HD 175511, deserves some further investigation.', '6. Conclusions': 'In this study, a detailed analysis of 31 stars that have been claimed as late B/early A type short period pulsators is presented. In this part of the H-R diagram no such pulsations are expected, so these objects need to be explained. There are four different hypotheses that are to be considered: incorrect T eff , binarity, rapid rotation and possible contamination of the light curves. To this end, a detailed spectroscopic investigation was carried out for our candidate variables. For all stars accurate stellar parameters such as T eff , log g , v sin i were determined. Furthermore, individual element abundances were derived for all stars where this was possible. In addition, a contamination analysis was carried out to be sure that the determined high pulsation frequencies arise from the stars of interest.  \nTable 3. Variability classification of our target stars.  \nThere are several stars in our sample whose frequency spectra appear phenomenologically similar, hence we analysed them in groups. The first group shows a wide frequency range extending to ∼ 40d -1 . Another group exhibits two dominant frequencies close to each other and the final group possesses frequencies between 1 - 10 d -1 . Taking into account these groups, the positions of the stars in the H-R diagram, v sin i values, the radial velocity measurements, and the RUWE parameter we discussed the proposed hypotheses of tthe existence of the so-called "Maia variables" located on the main-sequence in between the instability strips of the SPB and δ Scuti stars. First, we conclude that incorrect T eff values by themselves can not explain all of these stars even if for some systems there are significant differences between the spectroscopic and photometric T eff values (see Fig. 8). The binary hypothesis alone is also insufficient. When we examined the high rotation hypothesis by considering the systems\' positions in the H-R diagram and the structure of frequency spectra. In this way, we identified some stars whose frequencies were likely shifted towards higher values due to rapid rotation.  \nHowever, taking into account a combination of all hypotheses, most of the 31 stars we studied can be explained easily. 27 are members of established variability classes, ranging from Internal Gravity Wave pulsators, β Cephei stars, SPB stars to δ Scuti and γ Dor stars. The variability of one object is due to a contaminating neighbouring star. There are three stars that are not that easy to explain. Two of them are hot rotational variables which also exhibit short-period pulsations, but for both of them, we found indirect evidence for binarity. The last star is a rapidly rotating Be star in the SPB instability strip that also shows multiperiodic short-period pulsations, but also, in this case, the Gaia RUWE parameter is suspicious in terms of binarity.  \nExtrapolating our results to the hundreds of main sequence pulsators of B/A spectral type outside of any known instability region claimed recently (Balona, 2023), we expect that most if not all of these objects can be explained in different ways as well. From today\'s point of view, there is thus no need to claim the existence of an unexplained class of variable stars in this domain of the H-R diagram, let alone call them "Maia variables" as Maia itself is not a pulsating star (White et al., 2017).', 'Acknowledgement': "The authors would like to thank the reviewer for useful comments and suggestions that helped to improve the publication. This study has been supported by the Scientific and Technological Research Council (TUBITAK) project 121F474. GH thanks the Polish National Center for Science (NCN) for supporting the study through grants 2015/18/A/ST9/00578 and 2021/43/B/ST9/02972. Horizon 2020: (OPTICON) This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 730890. The calculations have been partly carried out using resources provided by the Wrocław Centre for Networking and Supercomputing (http://www.wcss.pl),  \nGrant No 214. This material reflects only the authors views and the Commission is not liable for any use that may be made of the information contained therein. Based on observations made with the Mercator Telescope, operated on the island of La Palma by the Flemish Community, at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofìsica de Canarias. The TESS data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). Funding for the TESS mission is provided by the NASA Explorer Program. This work has made use of data from the European Space Agency (ESA) mission Gaia (http://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, http://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular, the institutions participating in the Gaia Multilateral Agreement. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. We are grateful to Vichi Antoci, Tim Bedding, Dominic Bowman, Don Kurtz, Nami Mowlavi and Simon Murphy for their comments on a draft version of this paper.", 'References': '1997, ESA Special Publication, Vol. 1200, The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission  \n- Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology Aerts, C., & Kolenberg, K. 2005, A&A, 431, 615\n- Aerts, C., & Tkachenko, A. 2023, arXiv e-prints, arXiv:2311.08453\n- Amôres, E. B., & Lépine, J. R. D. 2005, AJ, 130, 659\n- Antoci, V., Cunha, M., Houdek, G., et al. 2014, ApJ, 796, 118\n- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Baker, N., & Kippenhahn, R. 1962, ZAp, 54, 114\n- Balega, Y. Y., Dyachenko, V. V., Maksimov, A. F., et al. 2012, Astrophysical Bulletin, 67, 44\n- Balona, L. A. 2023, Frontiers in Astronomy and Space Sciences, 10, 1266750\n- Balona, L. A., & Ozuyar, D. 2020, MNRAS, 493, 5871\n- Balona, L. A., Engelbrecht, C. A., Joshi, Y. C., et al. 2016, MNRAS, 460, 1318 Baran, A. S., & Koen, C. 2021, Acta Astron., 71, 113\n- Bowman, D. M., Burssens, S., Pedersen, M. G., et al. 2019, Nature Astronomy, 3, 760\n- Brown, T. M., Latham, D. W., Everett, M. E., & Esquerdo, G. A. 2011, AJ, 142, 112\n- Catanzaro, G., Leone, F., & Dall, T. H. 2004, A&A, 425, 641\n- Chang, S. W., Protopapas, P., Kim, D. W., & Byun, Y. I. 2013, AJ, 145, 132\n- Charpinet, S., Brassard, P., Fontaine, G., et al. 2019, A&A, 632, A90\n- Christophe, S., Ballot, J., Ouazzani, R. M., Antoci, V., & Salmon, S. J. A. J. 2018, A&A, 618, A47\n- Cunha, M. S., Alentiev, D., Brandão, I. M., & Perraut, K. 2013, MNRAS, 436, 1639\n- Degroote, P., Aerts, C., Ollivier, M., et al. 2009, A&A, 506, 471\n- Dommanget, J., & Nys, O. 1995, in Astronomical and Astrophysical Objectives of Sub-Milliarcsecond Optical Astrometry, ed. E. Hog & P. K. Seidelmann, Vol. 166, 395\n- Dziembowski, W. A., Moskalik, P., & Pamyatnykh, A. A. 1993, MNRAS, 265, 588\n- Eggleton, P. P., & Tokovinin, A. A. 2008, MNRAS, 389, 869\n- Fabricius, C., Høg, E., Makarov, V. V., et al. 2002, A&A, 384, 180\n- Flower, P. J. 1996, ApJ, 469, 355\n- Gaia Collaboration. 2022, VizieR Online Data Catalog, I/355\n- Gaia Collaboration, Smart, R. L., Sarro, L. M., et al. 2021, A&A, 649, A6 Gerbaldi, M., Floquet, M., & Hauck, B. 1985, A&A, 146, 341\n- Gray, D. F. 2008, The Observation and Analysis of Stellar Photospheres Higgins, M. E., & Bell, K. J. 2022, TESS-Localize: Localize variable star signatures in TESS Photometry, Astrophysics Source Code Library, record ascl:2204.005, ascl:2204.005\n- Holdsworth, D. L., Cunha, M. S., Lares-Martiz, M., et al. 2024, MNRAS, 527, 9548\n- Hubeny, I., & Lanz, T. 2011, Synspec: General Spectrum Synthesis Program, Astrophysics Source Code Library, record ascl:1109.022, ascl:1109.022 Jackiewicz, J. 2021, Frontiers in Astronomy and Space Sciences, 7, 102 Kahraman Aliçavuş, F., Niemczura, E., De Cat, P., et al. 2016, MNRAS, 458, 2307\n- Kochukhov, O., Khalack, V., Kobzar, O., et al. 2021, MNRAS, 506, 5328 Kurtz, D. W. 2022, ARA&A, 60, 31\n- Kurtz, D. W., Jayaraman, R., Sowicka, P., et al. 2023, MNRAS, 521, 4765\n- Kurtz, D. W., Shibahashi, H., Murphy, S. J., Bedding, T. R., & Bowman, D. M. 2015, MNRAS, 450, 3015\n- Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CD-ROM No. 13. Cambridge, 13\n- Kurucz, R. L., & Avrett, E. H. 1981, SAO Special Report, 391\n- Labadie-Bartz, J., Carciofi, A. C., Henrique de Amorim, T., et al. 2022, AJ, 163, 226\n- Labadie-Bartz, J., Handler, G., Pepper, J., et al. 2020, AJ, 160, 32 Lanz, T., & Hubeny, I. 2007, ApJS, 169, 83  \nLehmann, H., Scholz, G., Hildebrandt, G., et al. 1995, A&A, 300, 783 Lenz, P., & Breger, M. 2005, Communications in Asteroseismology, 146, 53 Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & Worley,  \n- C. E. 2001, AJ, 122, 3466\n- McNamara, B. J. 1985, ApJ, 289, 213\n- Moe, M., & Di Stefano, R. 2017, ApJS, 230, 15\n- Mowlavi, N., Barblan, F., Saesen, S., & Eyer, L. 2013, A&A, 554, A108\n- Mowlavi, N., Saesen, S., Semaan, T., et al. 2016, A&A, 595, L1\n- Murphy, S. J., Bedding, T. R., Shibahashi, H., Kurtz, D. W., & Kjeldsen, H. 2014, MNRAS, 441, 2515\n- Murphy, S. J., Hey, D., Van Reeth, T., & Bedding, T. R. 2019, MNRAS, 485, 2380\n- Pamyatnykh, A. A. 1999, Acta Astron., 49, 119\n- Przybilla, N., Nieva, M.-F., & Butler, K. 2011, in Journal of Physics Conference Series, Vol. 328, Journal of Physics Conference Series, 012015\n- Raskin, G., van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 Ratnasingam, R. P., Rogers, T. M., Chowdhury, S., et al. 2023, A&A, 674,\n- A134\n- Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, ed. J. Oschmann, Jacobus M., M. Clampin, G. G. Fazio, & H. A. MacEwen, 914320\n- Saio, H., Kurtz, D. W., Murphy, S. J., Antoci, V. L., & Lee, U. 2018, MNRAS, 474, 2774\n- Salmon, S. J. A. J., Montalbán, J., Reese, D. R., Dupret, M. A., & Eggenberger, P. 2014, A&A, 569, A18\n- Sharma, A. N., Bedding, T. R., Saio, H., & White, T. R. 2022, MNRAS, 515, 828\n- Shore, S. N., Brown, D. N., & Sonneborn, G. 1987, AJ, 94, 737\n- Stankov, A., & Handler, G. 2005, ApJS, 158, 193\n- Stassun, K. G., & Torres, G. 2021, ApJ, 907, L33\n- Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138\n- Struve, O. 1955, S&T, 14, 461\n- Szewczuk, W., & Daszyńska-Daszkiewicz, J. 2017, MNRAS, 469, 13 Telting, J. H., Avila, G., Buchhave, L., et al. 2014, Astronomische Nachrichten,\n- 335, 41  \nTody, D. 1986, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 627, Instrumentation in astronomy VI, ed. D. L. Crawford, 733  \nTownsend, R. H. D. 2005, MNRAS, 360, 465  \nTownsend, R. H. D., Goldstein, J., & Zweibel, E. G. 2018, MNRAS, 475, 879 Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125 Walczak, P., Daszyńska-Daszkiewicz, J., Pigulski, A., et al. 2019, MNRAS, 485, 3544  \nWhite, T. R., Pope, B. J. S., Antoci, V., et al. 2017, MNRAS, 471, 2882', 'Appendix 1. Contamination analysis of individual stars Appendix 1.1 TIC 64774351': "The SPOC-derived light curve for TIC 64774351 has only 90.3% of its flux arising from the star of interest. There are two stars with 9 ≤ G mag ≤ 10 within 100' of TIC 64774351 (BD+52 2951 and BD+52 2952); the former is especially a key potential source of contamination. A light curve created from an aperture that does not include these two stars still shows a highly significant peak at 0.203 d, as shown in the left panel of Figure 5. Another choice of aperture, centered on the pixels to the right of the star of interest, provides some evidence for the peak at ∼ 5d -1 being associated with this star. However, it is unable to resolve the forest of peaks at higher frequencies (e.g., those around 10-20 d -1 ) neither can the periodogram generated from the custom light curve for the aperture shown in Figure 5. However, the white noise level is higher in the periodogram created from the custom aperture, which could hide these peaks, which are already relatively low signal-tonoise in the SPOC light curve. The nearby star, BD+52 2951 (indicated in tan in the top left panel of Figure 5), has a listed spectral type of F8, suggesting that it is likely not in the classical instability strip and should thus not pulsate.  \nAn analysis of this star using the TESS -Localize tool suggests that the lower frequencies (those less than ∼ 5d -1 ) are highly likely to be associated with BD+52 2952. However, the analysis was unreliable for the higher frequencies (i.e., those greater than ∼ 5d -1 ), so in this case, we are unable to use TESS -Localize to draw any conclusions.", 'Appendix 1.2 TIC 251196433': "The SPOC-derived light curve for this star has only 94.1% of its flux arising from the star of interest. There exists a star with G mag = 14 within 4' (Gaia DR3 460981833479443456; henceforth referred to as GDR3-456); their locations in the TPF are shown in Figure 6. These two stars fall within the same pixel, making it nearly impossible to disentangle their respective variability via the use of custom apertures. LabadieBartz et al. (2020) identified TIC 251196433 as a β Cep star based on its low-frequency variability (most of the peaks are at frequencies less than 10 d -1 ). However, this is also the regime of certain compact g mode pulsators.  \nGiven that GDR3-456 has a similar parallax to TIC 251196433, they likely lie at similar distances, and so we can assume that the extinction coefficients (and thus, color excess) for both stars are roughly identical. Using the Gaia value of E( G BP -G RP ) for TIC251196433, we see that the color of GDR3-456 is around -0.17, which could make it a cool sdB star. The pulsation  \namplitude is 0.7% for the highest amplitude frequency in the periodogram, which is too high for an sdB g mode pulsator and is thus likely associated with TIC 251196433. However, the other peaks in the periodogram represent pulsation amplitudes of 0.15% or less; these modes could arise from the contaminating star GDR3-456. d Further analyses, including speckle imaging and ground-based spectroscopy, may be needed to tease out the properties of this star and separate its modes from TIC251196433.", 'Appendix 1.3 TIC 255996252': "The SPOC light curve for this star has only 90.8% of its flux arising from the star of interest. As for TIC 251196433, and as shown in Figure 6, this star is nearly coincident in the sky with another star, Gaia DR3 2274443263522286208 ( G mag = 12.49; henceforth GDR3-208). This star is located at a similar distance as TIC 255996252, at around 210 pc based on parallax. GDR3208, however, does not have any information about its color from Gaia DR3 or prior survey missions such as Hipparcos. Assuming similar extinctions A G (a value is provided for this for TIC255996252 in Gaia DR3), we get an absolute magnitude in G -band of 5.4, suggesting that it is likely a G-type star. Given that such stars' pulsation modes are convectively driven and exhibit clear frequency spacings, this star likely does not contaminate the light curve of TIC 255996252, and all the peaks seen in the star's Fourier spectrum can be taken to be intrinsic to this star.", 'Appendix 1.4 TIC 352781227': "The SPOC light curve for this star has only 84.5% of the flux arising from TIC 352781227. This star has G mag = 11.28, and there is a nearby star (within 30') with G mag = 12.82. We created a custom aperture for the target pixel file containing just this star (shown in light blue shading in Figure 5, compared to the SPOC aperture in red shading). While this of course does not contain the entirety of the flux from this star (and so the peak amplitudes in the periodogram will be correspondingly lower), the periodogram should capture the frequencies of the peaks intrinsic to this star's variability. We find that the noise level of the periodogram is too high in the low-frequency regime ( ν ≲ 3d -1 ), so no conclusions can be drawn about these frequencies; however, the higher-frequency ones ( ν ≥ 3d -1 ) are very likely intrinsic to the star.  \nAnother aperture that includes only the two topmost pixels in the SPOC aperture exhibits a peak at 6.5 d -1 , suggesting that it is intrinsic to the star. No evidence for higher frequencies is present, and the noise level at low frequencies is too low to draw any meaningful conclusions.", 'Appendix 1.5 TIC 424032547': 'Of all the stars in our sample, this star has the highest likelihood for contamination, as only 72.7% of the flux in the  \nd For an example of sdB stars pulsating in this regime, see Charpinet et al. (2019).  \nFigure A1. Amplitude spectra of the targets ordered by decreasing T eff . The TESS data sectors used for those plots are indicated.  \n<!-- image -->  \naperture comes from the star of interest. There is one significant peak at 5.03 d -1 , with at least 5 significant harmonics in the periodogram. An analysis of this TPF with TESS -Localize suggests that the probability that this frequency (and, by extension, its harmonics) is associated with that particular star is quite low. Indeed, as noted in Labadie-Bartz et al. (2020), this star is located in a dense region of an open cluster, such that the light curve and target pixel file suffer from significant contamination and saturation that could affect the SPOC light curve.  \nFurther analysis using TESS -Localize provides some support to our conclusion that the lower frequencies arise from this star. The fits are not significant, but have p < 0.2, suggesting that this conclusion may be somewhat likely.', 'Appendix 2. Discussion of individual stars': 'Here we discuss each star according to their positions in the final H-R diagram, starting from the hottest objects and progressing towards lower T eff . We give arguments about whether the given objects can be considered members of an established class of variable stars, or whether they cannot be explained this way.', '). Appendix 2.1 Stars in the β Cephei domain': 'TIC 202567458 = HD 332044: most of the Fourier spectrum of this star\'s pulsations can be described by a basic multitude of signals between 2.91-3.24 d -1 and a single frequency of 1.46 d -1 as well as a surprisingly long series of peaks at combination multiples thereof. Pressure modes are hard to discern. In addition, there is a 1/f "noise floor" reminiscing the signals usually connected with Internal Gravity Waves (IGWs, e.g., see Ratnasingam et al. 2023). Given that the star lies close to the zero-age main sequence in Fig. 7, we conclude that the apparently coherent signals are due to gravity mode pulsations and their combination frequencies (Kurtz et al., 2015). Thus there is no reason for calling this object a "Maia" variable.  \nTIC 312637783 = HD 228456: this star shows multiple pulsation frequencies between 8-9 d -1 and 11 d -1 as well as some combination frequencies thereof. It has a high RUWE parameter of 3.53 (cf. Table 2). With this and the position of the star in Fig. 7, we classify this star as a suspected binary containing at least one β Cephei pulsator. There is nothing unusual about this object.  \nTIC 374038418 = BD+60 416: the Fourier spectrum is dominated by three coherent frequencies typical of β Cephei pulsators. Superposed on that is some non-coherent variability which could be due to IGWs or mass loss. All these phenomena are consistent with this object\'s temperature and luminosity put it in the β Cephei instability strip near the end of the main sequence phase. There is no need to invoke a possible "Maia" phenomenon.  \nTIC 377099704 = HD 177195: this star is located at the red edge of the β Cephei instability strip and inside the instability strip of the Slowly Pulsating B (SPB stars). Its pulsation spectrum is very complicated, with many modes in both the expected p and g mode frequency domains that do not seem to be harmonically related. We classify it as a "hybrid" β Cep/SPB pulsator.  \nTIC 251196433 = BD+57 655: this is a very rapidly rotating star with a large RUWE value suggestive of binarity. In Sect.3.1 we also commented on its line-profile shape atypical of a single star. This object is located near the edges of both the β Cep and SPB instability strips. Its pulsation spectrum is dominated by an apparent p mode, and has many frequencies in  \nAdditional verification that the observed signal does not arise from this particular star lies in the VARIABLE flag assigned to stars in the Gaia catalog. There are two nearby stars (designated Gaia DR3 1823376830245623296 and DR3 1823376825919999232), which are somewhat fainter; however, both are designated with the VARIABLE flag in Gaia, while TIC 424032547 is not (its value for this flag is NOT\\_AVAILABLE While this is not a sufficient condition to say that the observed variability does not arise from TIC 424032547, it is further strong evidence that suggests that the SPOC light curve does not capture the intrinsic variability for this star. Additionally, an aperture with the rightmost 4 pixels from the SPOC aperture does not show any evidence of variability at the observed frequencies, providing further evidence that the observed signals in the periodogram do not arise from the star of interest.  \nThe nearby contaminants are shown in the left panel of Figure 6.', 'Appendix 1.6 TIC 450085823': "There is a nearby star with G mag = 14.04, as shown in the rightmost panel of Figure 5. This star contributes 2.73% of the flux to the aperture used to construct the SPOC light curve, so it will likely have a minimal impact on the frequencies recovered from the light curve. We constructed a custom light curve using a two-pixel aperture and recovered most of the large amplitude peaks visible in the SPOC light curve, albeit at a lower S/N ratio. The peak at ν ∼ 1.8 d -1 is not visible in the custom light curve, perhaps due to the fact that it arises from the contaminating star. However, it could also simply be hidden below the noise. The high-frequency peaks are still visible in the custom light curve's analysis.", 'Appendix 1.7 TIC 467219052': 'This star has two potential nearby contaminants with Gaia magnitudes 12.06 and 12.94 (as seen in the third panel from left in Figure 5); the SPOC light curve captures 88% of the flux from the star of interest. A custom light curve from the single-pixel aperture shows that the largest peaks below 7d -1 are likely intrinsic to this star. The peak at ∼ 11.8 d -1 is significant in the custom light curve but is not present in the SPOC light curve. This may arise from some systematic variability in the custom light curve that was unable to be removed using the CBVs. Selecting the 6 leftmost pixels in the SPOC aperture and constructing a light curve yields a periodogram with significant peaks that agree with the SPOCgenerated light curve.  \nboth the β Cep and SPB domains plus some combination frequencies. According to us this is another "hybrid" β Cep/SPB pulsator.', 'Appendix 2.2 Stars in the SPB domain': 'TIC 354793407 = HD 21699: this star\'s light curve is dominated by a 2.5-d variation that is linked to rotation (Shore et al., 1987). Superposed on that is low-amplitude multiperiodic variability with periods around 40 minutes, typical for unevolved δ Scuti stars. At T eff = 16000K, this star is too hot to be a δ Scuti star (cf. Table 2). However, it also has a large RUWE of 1.81, suggestive of binarity. A hypothesized main sequence δ Scuti companion would be expected to have about 1/60 of the luminosity of the chemically peculiar primary. Therefore the small (< 0.23 mmag) short-period pulsation amplitudes in the TESS light curve would have intrinsic values of about 14 mmag and less, quite plausible for such an object. We therefore suspect that a binary companion is responsible for the short-period pulsations.  \nTIC 36557487 = HD 49643: similar to TIC 202567458, this star shows groups of closely spaced frequencies between 1.98 and 3.13 d -1 and harmonics thereof, that can be discerned out to 37 d -1 in the Sector 6 data and that we analogically classify as g mode pulsations. Apart from that, in TESS Sector 33 the light curve shows a sudden spectacular increase in total pulsation amplitude by a factor of about 15 in connection with a brightening event. We, therefore, suspect this is a hitherto unrecognized Be star as discussed e.g. by (Labadie-Bartz et al., 2022). In any case, we find no reason why this object should be a "Maia" variable.  \nTIC 424032547 = HD 350990: the light curve of this object is dominated by a frequency of 5.025 d -1 whose amplitude is different within the three sectors it was observed. Interestingly, the amplitude of this variation is inversely correlated with the CROWDSAP parameter, i.e. the larger the flux contribution of the target in the photometric aperture, the lower the amplitude. This suggests that this object is not the source of the variability (cf. Sect.Appendix 1.5). Furthermore, careful prewhitening of the dominant frequency and its harmonics reveals the presence of a weak subharmonic. A phase fold of the light curve around the subharmonic frequency yields a light curve shape indicative of an ellipsoidal variable.  \nTIC 236785664 = HD 175511 is classified as a late-type Be star. The Fourier spectrum of its TESS light curve shows a single frequency near 1.24 d -1 and a multitude of peaks between 12 and 43 d -1 that are not harmonically related. Such pulsation frequencies would be typical for a δ Scuti star. At T eff = 12500K this star is much hotter than such pulsators. Even its rapid rotation v sin i = 260 km s -1 should not create sufficient oblateness that its equatorial regions would be cool enough to excite pulsations. Only its Gaia RUWE parameter (1.20) gives a hint of possible binarity, which by itself we do not find sufficient to rule this object out as a "Maia" candidate. One may suspect that this star could be a case related to HD 42477. However, an examination for possible coupling between the g and p modes yielded a null result.  \nTIC 192780911 = HD 19921: the TESS data show little intrinsic variability except for a weak signal at 3.9 d -1 , which could be an SPB-type pulsation modified by the rapid rotation of the star ( v sin i = 273 km s -1 ). The RUWE value of 1.77 is suggestive of binarity.  \nTIC 16781787 = HD 207674: this object\'s variability is predominantly singly periodic with a frequency of 6.0 d -1 . Two very weak signals near 2.35 d -1 are present in the S56 data, but are not subharmonics. As the strongest frequency is too high to be caused by rotation or binarity we speculate that this star\'s very rapid rotation shifts the dominant pulsation frequency with an estimated amount of > 2 d -1 to such a large value for an SPB star.', 'Appendix 2.3 Stars between the SPB and δ Scuti strips': 'TIC 64774351 = HD 204905: the frequency spectrum of this star changed between the Sector 15/16 and 56/57 observations, but the range of dominating frequencies remained around 4 - 5 d -1 . The Fourier Transform contains peaks between 1 31d -1 that do not appear to cluster in groups. However, a closer look reveals that the higher frequencies can in all cases be explained by combinations of lower frequencies, with these few base frequencies located in the g mode domain. We therefore surmise that this star is related to those with frequency groupings. The high-frequency signals are then not intrinsic normal modes of pulsation but combination frequencies.  \nTIC 26369887 = HD 222142: there is a frequency group between 3.05-3.5 d -1 that, together with a single frequency of 2.47 d -1 , appears to form all other observed frequencies by linear combinations. This object then is a very cool SPB star. Its rapid rotation could have an influence on pulsational excitation. There are no high frequencies that would make this object a "Maia" candidate.  \nTIC 83803744 = HD 1009 is a HgMn star that has two periods consistently present in the TESS data sets, one at 0.89 d -1 that is apparently the first harmonic of its rotation period (Kochukhov et al., 2021) and a low amplitude (0.04 - 0.05 mmag) apparent pulsation at 12.36 d -1 . Even though our two radial velocity measurements are consistent within the errors, it is known that HgMn stars often occur in binary systems (e.g., Gerbaldi et al. 1985) which is why we are reluctant to classify the apparent pulsation as "Maia" variability: the pulsations could arise from a δ Scuti companion.  \nTIC 20534584 = HD 176052: the low-frequency part of the Fourier Transform of the TESS light curve of this fairly rapid rotator ( v sin i = 150 km s -1 ) is constituted by harmonically related frequencies based on a group around 3.35-3.40 d -1 . A single significant higher frequency peak occurs at 22.64 d -1 with an amplitude of only 0.16 mmag (the peak at 27.69 d -1 is not significant). Given the high RUWE value of this object (2.26), we suspect that the high-frequency variation comes from a companion of the δ Scuti type.  \nTIC 352781227 = HD 351435: this is a very rapid rotator with a consequently large uncertainty in its effective temperature (cf. Table 2). The pulsation spectra are clearly variable in amplitude from sector to sector of data and even within single  \nsectors, but the underlying frequencies are essentially the same, suggesting these are unresolved groups of frequencies. The signals below 10 d -1 can be fairly easily explained with combinations of two basic frequencies near 1.43 and 3.84 d -1 , and the higher amplitudes of the combination frequencies with heavy geometrical cancellation of the base frequencies (see Kurtz et al., 2015). The single higher frequency signal at 13.66 d -1 can be matched with a third-order combination frequency and may also be an unresolved group. Thus this object could be a pure g mode pulsator in spite of the observed high frequencies or a p mode that is driven in the equatorial regions of this oblate star(Kurtz et al., 2023).  \nTIC 434340780 = BD+61 2380 is a rapidly rotating star of late B spectral type. Its two pulsation frequencies are too low for p modes at this star\'s position in the HR Diagram, which is why we interpret them in terms of g modes shifted to high frequency by rapid rotation ( f rot ≈ 2d -1 ).', 'Appendix 2.4 Stars in the δ Scuti strip': 'The behaviour of all these objects is consistent with that of δ Scuti pulsators, which is also consistent with their positions in the HR Diagram. Nevertheless, they will be discussed briefly in what follows.  \nTIC 201739321 = HD 356 is located at the blue edge of the δ Scuti instability strip. The short-period variations can therefore be classified as δ Scuti pulsations. The two signals below 10 d -1 are harmonically related.  \nTIC 317512446 = HD 220855 is located just inside the δ Scuti instability strip. This is consistent with its two shortperiod pulsation frequencies.  \nTIC 27804376 = HD 210661 and TIC 26368855 = HD 222017 are two multiperiodic δ Scuti stars near the hot border of their pulsational instability strip with intermediate projected rotational velocities.  \nTIC 255996252 = HD 196358 has similar characteristics to the two stars just discussed, but its pulsation spectrum is not as rich.  \nTIC 450085823 = HD 287149 shows multiple pulsation frequencies of the δ Scuti type below 20 d -1 . Its pulsation spectrum exhibits considerable amplitude changes between its two TESS light curves taken two years apart.  \nTIC 402910080 = HD 193989 is a very rapidly rotating star with a wide range of δ Scuti type pulsation frequencies observed. Amplitude changes between the different epochs of TESS observations are easily discernable.  \nTIC 387757610 = HD 200506 is another rapidly rotating δ Scuti star. Its pulsation spectrum appears fairly stable over time and it also exhibits g mode pulsations of the γ Doradus type plus combination frequencies thereof.  \nTIC 377443211 = HD 185757 also rotates rapidly. It shows a wide range of δ Scuti type pulsation frequencies whose amplitudes somewhat change over time.  \nTIC 90647390 = HD 278751 shows a few δ Scuti type pulsations and some combination frequencies thereof. Its Gaia RUWE value strongly suggests binarity.  \nTIC 467219052 = HD 239419 is the slowest pulsating  \nδ Scuti star in our sample, and it is also farthest away from the zero-age main sequence. Its pulsation amplitudes vary somewhat over time.  \nTIC 367351021 = HD 216486 has the highest RUWE value among our targets and a set of almost equally spaced δ Scuti type pulsation frequencies.  \nTIC 50850670 = BD+65 291 belongs to the δ Scuti stars of intermediate rotational velocities. It has three distinct groups of pulsation frequencies, the lowest frequency group likely due to γ Doradus type pulsations. The two other groups are separated by a factor of two in frequency but the individual modes are not harmonically related.  \nTIC 335484421 = BD+62 2044 is the coolest δ Scuti star in our sample, with a wide range and a large number of pulsation frequencies.  \nTable A1. Radial velocity measurements. * represents the cannot measure radial velocity value due to low signal-to-noise level.  \nTable A1. Continuation.  \nTable A2. The result of the analysis of chemical abundances. The numbers given in the brackets represent the number of used lines in the analysis.Table A2. Continuation.  \nTable A2. Continuation.'}