\section{Introduction} \label{intro}


The {\it Kepler} spacecraft \citep{2010ApJ...713L..79K, 2010Sci...327..977B} was originally launched to detect Earth-like transiting planets, but 
it also shed new light on many other aspects of stellar astrophysics \citep[e.g.][]{2019arXiv190605587P, 2019MNRAS.485..560C, 2019MNRAS.484..451H}.
In particular, new discoveries about pulsating stars
provided important input for asteroseismology, opening new frontiers and posing new questions.
Some of those concern the classes of the A- and F-type
pulsating variables located on or close to the main sequence, i.e., 
$\delta$\,Scuti (\dst) and $\gamma$\,Doradus (\gd) stars.


The \dst\, stars generally exhibit pressure ($p$)
modes which are excited by the $\kappa$ mechanism in the \ion{He}{ii} ionization zone. 
They typically oscillate with frequencies higher than $\sim$5~d$^{-1}$. 
The \dst\, stars are mostly placed in the lower part of the classical instability strip. 
Recently, \citet{2019MNRAS.485.2380M} classified a sample of over 15\,000 {\it Kepler} A-type  and F-type targets into
\dst\, and non-\dst\, stars, also providing a subdivision in groups on the basis of the observed pulsation properties. They also
defined a new empirical instability strip for \dst\, stars.

The \gd\, stars are a little cooler than the \dst\, variables. The gravity ($g$) mode oscillations of 
\gd\, stars are believed to be excited by the convective blocking mechanism \citep{2000ApJ...542L..57G} with frequencies typically 
lower than 5~\cd. 
In a recent study, it was suggested that the \gd-type pulsations are caused by the combination of $\kappa$ mechanism and 
the convection-oscillation coupling \citep{2016MNRAS.457.3163X}. 
The \gd\, domain
partially overlaps the cool border of the theoretical instability strip of \dst\, stars. In this part, new variables 
called A-F type hybrid stars were predicted \citep{1999MNRAS.309L..19H, 2004A&A...414L..17D}. These stars display both \dst\, and \gd\, type pulsations,
i.e., $p$ and $g$ modes. 

Before {\it Kepler}, only a few
hybrid stars had been discovered by means of ground-based observations \citep{2005AJ....129.2026H,2008A&A...489.1213U,2009MNRAS.398.1339H} and
their positions in the Hertzsprung-Russell (H-R) diagram  matched very well with
the theoretically predicted region. The precision and near continuous nature of {\it Kepler} photometry revolutionized the field, showing that apparently there are many
A-F type hybrid stars located beyond the theoretical region \citep{2010ApJ...713L.192G, 2011A&A...534A.125U}.
On the other hand, it was also found that some \dst\, and \gd\, stars are also placed outside the respective instability strips.
These results give rise to conflicts with theory. Therefore, it was timely to fix the exact positions
of these variables  in the H-R diagram and to revise the borders of their instability strips. For these reasons, spectroscopic
studies were carried out to determine accurate atmospheric parameters of \dst, \gd, and hybrid
stars (e.g. \citealt{2012MNRAS.422.2960T, 2015MNRAS.450.2764N, 2017MNRAS.470.2870N}).
These investigations  confirmed that some stars are actually outside the theoretical instability strips.

In particular, \citet{2014MNRAS.437.1476B} and \citet{2016MNRAS.460.1318B} 
noticed how 
 some gravity-mode pulsators are located close to, and even beyond,   
 the hot border of \dst\, instability strip, defining them as hot \gd\, stars. These variables are 
remarkable objects because 
if the theory of convective driving would apply to stars 
located in this part of the H-R diagram then they should not have a  
sufficient convective envelope to drive the \gd-type pulsations. 
It has been suggested that the hot \gd\, stars could be A- or B-type stars with a cooler
\gd\, companion,  or simply stars with wrong determinations of the effective temperature (\teff).
A spectroscopic study of hot \gd\, stars was undertaken \citep{2016MNRAS.460.1318B} and
it turned out that the resulting \teff\, values were mostly consistent with those given in the 
{\sl Kepler} input catalog \citep[KIC; ][]{2011AJ....142..112B}. However, the binary nature of the stars 
was not assessed and also the chemical composition of these variables could not be probed due to the low resolution 
of the spectra.
 
Another explanation is that the hot \gd\, stars are actually rapidly-rotating slowly pulsating B (SPB) stars 
\citep{2016MNRAS.460.1318B}. 
Due to gravity darkening, their equatorial zones appear cooler than the rest of the surface, then they are classified as hot A-stars and, hence, 
as hot \gd\, variables.
If this hypothesis is true, all hot \gd\, stars should rotate with high rotational velocities. 
There are also some hot hybrid stars which are located in the same area of the H-R diagram. 
In both cases the \gd-type pulsation conflicts with the theory, 
since would require that the convective blocking mechanism is continuing to be active in the hottest A stars \citep{2014MNRAS.437.1476B}. 
The $g$ mode pulsation in those \gd\, stars can also be explained by the radiative $\kappa$ mechanism and the 
coupling between oscillation and convection \citep{2016MNRAS.457.3163X}.


\begin{figure}
\includegraphics[width=8cm,angle=0]{Fig1_NOT_firts_positions.ps}
\caption{The positions of the selected candidate $\gamma$ Dor (star symbols) and hybrid (diamonds) stars.
The parameters of the stars were taken from \citet{2014ApJS..211....2H}. The theoretical instability strips of the 
$\gamma$ Dor (dashed lines) and $\delta$ Sct (solid lines) stars were taken from \citet{2005A&A...435..927D}. 
The evolutionary tracks (Z\,=\,0.02) were 
adopted from \citet{2016MNRAS.458.2307K}.}
\label{NOT_firts_positions}
\end{figure}


In this study, a new spectroscopic survey of a selected sample of hot \gd\, and hybrid stars is presented. We focused on the candidate \gd\,
and hybrid stars located close to the hot border of the \dst\, instability strip.
We aim to answer the following questions. First, are the \teff\, values of these hot objects correct? Second,
do the hot \gd\, stars have hotter A- or B-type companions? Third, are the hybrid stars members of binary or
multiple systems? Can those hot variables be rapidly-rotating stars?
Additionally, are  hot \gd\, and hybrid stars chemically peculiar?
Therefore, we present a detailed spectroscopic analysis to determine the atmospheric parameters
(\teff, surface gravity \logg, and microturbulence velocity $\xi$), projected rotatinal velocities, binarity effects, and surface chemical abundances.



The target selection is described in Sect.~\ref{tarsel}. 
The spectroscopic observations, data reduction and normalisation in Sect.~\ref{spobs}. The determinations of the atmospheric parameters from photometric indices, spectral energy 
distribution (SED), the analysis of Balmer and metal lines are discussed in Sect.~\ref{atpar}. The chemical abundance analysis is presented in Sect.~\ref{anche}. The frequency analysis of the targets is performed in Sect.~\ref{freana}.
The discussion and conclusions are given in Sects.~\ref{disc} and \ref{conclu}, respectively. 


\section{Target Selection}\label{tarsel}

We started the selection process from the  \gd\, and hybrid star candidates listed in the study 
of \citet{2011A&A...534A.125U}. They were proposed on the basis of their {\sl Kepler} light curves and preliminary  atmospheric parameters listed in the KIC. 
To better characterize them, 
we used the improved parameters reported in the revised 
{\sl Kepler} catalog of 
\citet[][hereafter H14]{2014ApJS..211....2H}, the latest available at the time of our observations (October 2016). 
This catalog is based on a compilation of literature values for atmospheric properties  derived from 
a variety of observational techniques.

High-resolution spectroscopy of a well-constrained sample of stars should provide reliable answers to the questions posed in 
Sect.\,\ref{intro}, especially when combining this technique 
with a homogeneous reduction and analysis of the spectra taken with a single instrument.

To refine our selection, we considered that the \gd\, stars have \teff\, values in the  range  6900\,$-$\,7300\,K \citep{2015ApJS..218...27V} and that 
the hybrid stars are expected to be found at the intersection  of the \gd\, and \dst\, instability strips,
 where \teff\, changes approximately from 6600 to 7300\,K. 
Given that the typical \teff\, uncertainty is $\sim$300~K \citep[][]{2017MNRAS.470.2870N} and hot \gd\, stars 
definition by \citet{2016MNRAS.460.1318B}, we selected 
 \gd\, and hybrid candidates with \teff\,$\geq$\,7500\,K. 


In the end, our selection contained nine hot \gd\, and fifteen hot hybrid candidate stars 
(Table\,\ref{infotable}), for which high-resolution spectroscopy was not available. 
Some of them were observed with  low-resolution spectroscopy 
\citep{2016A&A...594A..39F, 2016MNRAS.460.1318B}, but  
it is not possible to obtain reliable \logg\, and metallicity values from this technique. 



The positions of the targets in the log\,\teff\,$-$\logg\, diagram are shown in Fig.\,\ref{NOT_firts_positions}. As can be 
seen from the figure, all stars are located beyond the blue edge of \gd\, instability strip. 
The \teff\, values of the selected targets were also checked by using the updated parameters of the stars \citep[][hereafter M17]{2017ApJS..229...30M}. 
It turned out that  all targets except KIC\,11508397 ($\sim$400\,K cooler) still have \teff\, values in the same range of the selection criteria. 
We used the M17 parameters in the following investigations. 



\section[]{Spectroscopic Observations}\label{spobs}

The stars were observed with the Fibre-fed \'{E}chelle Spectrograph (FIES),  a cross-dispersed 
spectrograph mounted on the 2.56-m Nordic Optical Telescope of the Roque de los Muchachos Observatory in La Palma \citep{2014AN....335...41T}. 
The spectrograph offers three resolving power options. The maximum resolving power is $R$\,=\,67\,000, while the medium and 
low resolving powers are $R$\,=\,46\,000 and $R$\,=\,25\,000, respectively, covering the wavelength range of $3700-8300$\,{\AA}. 

Since we aimed to derive atmospheric parameters and chemical abundances of the targets, 
we opted for the medium-resolution ($R$\,=\,46\,000) configuration,  
taking into account the average brightness of the sample ($V$\,=\,10.6 mag).
Observations were performed in the  first halves of the nights from  October 13 to 19, 2016.
To examine the binarity nature of the targets, we tried to take at least two spectra per star on different nights. 
For some stars, only one spectrum could be taken due to the weather conditions and the limited observing time. 
The number of the spectra for each star and the signal-to-noise (S/N) ratios at around 5500\,{\AA} for combined spectra
are given in Table\,\ref{infotable}.

The spectra were reduced by the dedicated pipeline FIEStool \citep{2014AN....335...41T}. The standard reduction procedure 
was applied. Bias subtraction, correction of flat-field, scattered light extraction, wavelength calibration, 
and merging of orders were performed for each spectrum. Normalisation of the reduced spectra was carried out 
manually by using the \textit{continuum} task of the NOAO/IRAF package\footnote{http://iraf.noao.edu/}. 



\begin{table}
\centering
  \caption{Log of the observations (from October 13 to 19, 2016):
number of the star (N$_i$) used in the paper, KIC ID, V-magnitude, signal-to-noise ratio (S/N), 
number of taken spectra (N$_S$), and input classification \citep {2011A&A...534A.125U}.} 
    \label{infotable}
  \small
  \begin{tabular*}{0.95\linewidth}{@{\extracolsep{\fill}}lrcccr}
  \hline
\noalign{\smallskip}
N$_i$  &KIC  	  & V     &  S/N  & N$_s$  & \multicolumn{1}{c}{Input} \\
   &           & [mag] &       &    &\multicolumn{1}{c}{classifcation}\\
\noalign{\smallskip}
 \hline
\noalign{\smallskip}
1	      &2168333  & 10.02 &   90 & 2 & Hybrid\\
2         &3119604  & 10.90 &   50 & 1 & Hybrid\\    
3         &3231406  & 10.38 &   85 & 2 & Hybrid\\
4         &3240556  & 10.10 &   75 & 2 & Hybrid\\
5         &3245420  & 10.52 &   75 & 2 & Hybrid\\
6         &3868032  & 10.44 &   70 & 2 & \gd\\
7         &4677684  & 10.19 &   85 & 2 & \gd\\
8         &4768677  & 10.91 &   60 & 2 & Hybrid\\    
9         &5180796  & 10.11 &   85 & 2 & \gd\\
10        &5630362  & 10.69 &   70 & 2 & \gd\\
11        &6199731  & 10.94 &   50 & 1 & Hybrid\\
12        &6500578  & 10.77 &   80 & 2 & \gd\\
13        &6776331  & 10.71 &   50 & 1 & Hybrid\\
14        &7694191  & 10.94 &   45 & 1 & \gd\\    
15        &7732458  & 10.85 &   75 & 2 & Hybrid\\
16        &9052363  & 10.64 &   60 & 2 & Hybrid\\
17        &9775385  & 11.05 &   45 & 1 & Hybrid\\
18        &10281360 & 11.06 &   50 & 1 & \gd\\
19        &11197934 & 10.81 &   60 & 2 & Hybrid\\
20        &11199412 & 10.90 &   50 & 1 & \gd\\    
21        &11508397 & 10.65 &   80 & 2 & Hybrid\\
22        &11612274 & 10.44 &   85 & 4 & \gd\\
23        &11718839 & 10.73 &   60 & 2 & Hybrid\\
24        &11822666 & 10.69 &   70 & 2 & Hybrid\\
\hline
\end{tabular*}
\end{table}



\section{Determination of the Atmospheric parameters}\label{atpar}



The spectroscopic atmospheric parameters (\teff, \logg, and $\xi$) were determined using the Balmer and iron lines, 
as done in previous several studies \citep[e.g., ][]{2010EAS....43..167N, 2012MNRAS.422.2960T, 2017MNRAS.470.2870N}. 
The approach we used in the analysis of the Balmer lines has been successfully applied in other papers \citep{2011MNRAS.411.1167C,2012MNRAS.421.1222C, 2013MNRAS.431.3258C}.
In practice, the procedure minimized the difference between observed and synthetic spectra, using the $\chi^2$ as
goodness-of-fit parameter. 
Since the rotational velocity affects the profile of the lines, we determined initial estimates of
projected rotational velocity (\vsini) values  by using the least-squares deconvolution (LSD) technique \citep{1997MNRAS.291..658D}.
The \vsini\, values were measured from the zero positions of the Fourier transform of the mean line profiles (Table\,\ref{initialvsini}).

\begin{table*}
\centering
\footnotesize
  \caption{The atmospheric parameters taken from M17 and our \teff\,values derived from the analysis of 
  Balmer lines. ${\it E(B-V)}$ values (Sect.~\ref{gaia}) are also listed.}
  \label{initialvalues}
  \begin{tabular*}{0.95\linewidth}{@{\extracolsep{\fill}}lrccccc@{}}
\toprule 
  Number   & \multicolumn{1}{c}{KIC} &  {$E(B-V)$}     & \teff$_{M17}$         & \logg$_{M17}$          &  \teff$_{\rm Hlines}$           \\
           &                         & [mag]           &     [K]               &    [dex]               &    [K]                       \\
           &                         &  $\pm$\,0.02    &                       &                        &                 		                  \\
\midrule
   1        &2168333 &  0.03           & 8363$^{+197}_{-395}$ &  3.80$^{+0.36}_{-0.15}$ &   8000\,$\pm$\,300       \\
   2        &3119604 &  0.03           & 8383$^{+233}_{-350}$ &  4.10$^{+0.14}_{-0.15}$ &   7900\,$\pm$\,360      \\
   3        &3231406 &  0.02           & 8005$^{+222}_{-333}$ &  3.77$^{+0.40}_{-0.07}$ &   7600\,$\pm$\,270       \\
   4        &3240556 &  0.03           & 8615$^{+238}_{-374}$ &  3.78$^{+0.41}_{-0.14}$ &   8000\,$\pm$\,300      \\
   5        &3245420 &  0.11           & 8079$^{+225}_{-338}$ &  3.72$^{+0.43}_{-0.11}$ &   7900\,$\pm$\,320       \\
   6        &3868032 &  0.10           & 8564$^{+234}_{-402}$ &  4.04$^{+0.16}_{-0.16}$ &   8300\,$\pm$\,610       \\
   7        &4677684 &  0.22           & 8602$^{+68}_{-94}$   &  3.78$^{+0.28}_{-0.11}$ &   8600\,$\pm$\,480       \\
   8        &4768677 &  0.20           & 8584$^{+77}_{-86}$   &  3.76$^{+0.26}_{-0.03}$ &   8800\,$\pm$\,830       \\
   9        &5180796 &  0.15           & 8076$^{+64}_{-96}$   &  3.82$^{+0.22}_{-0.09}$ &   8100\,$\pm$\,280       \\
    10      &5630362 & 0.11            & 7821$^{+78}_{-78}$   &  3.88$^{+0.17}_{-0.07}$ &   7700\,$\pm$\,170      \\
    11      &6199731 & 0.05            & 8040$^{+251}_{-306}$ &  3.64$^{+0.55}_{-0.09}$ &   7500\,$\pm$\,320       \\
    12      &6500578 & 0.28            & 8072$^{+223}_{-363}$ &  3.72$^{+0.42}_{-0.10}$ &   8200\,$\pm$\,450      \\
    13      &6776331 & 0.08            & 7870$^{+244}_{-325}$ &  3.70$^{+0.46}_{-0.08}$ &   7600\,$\pm$\,400       \\
    14      &7694191 & 0.19            & 8070$^{+251}_{-334}$ &  3.66$^{+0.50}_{-0.09}$ &   8300\,$\pm$\,580      \\
    15      &7732458 & 0.02            & 7766$^{+216}_{-324}$ &  3.64$^{+0.49}_{-0.09}$ &   7300\,$\pm$\,230      \\
    16      &9052363 & 0.11            & 7810$^{+216}_{-324}$ &  3.74$^{+0.42}_{-0.10}$ &   7800\,$\pm$\,360      \\
    17      &9775385 & 0.05            & 7675$^{+211}_{-316}$ &  4.05$^{+0.17}_{-0.15}$ &   7300\,$\pm$\,300      \\
    18      &10281360& 0.06	       & 7776$^{+216}_{-325}$ &  4.02$^{+0.22}_{-0.14}$ &   7200\,$\pm$\,280      \\
    19      &11197934& 0.06            & 7641$^{+68}_{-91}$   &  3.83$^{+0.22}_{-0.06}$ &   7600\,$\pm$\,360      \\
    20      &11199412& 0.02            & 7693$^{+239}_{-319}$ &  3.72$^{+0.46}_{-0.08}$ &   7200\,$\pm$\,280      \\
    21      &11508397& 0.00            & 7287$^{+76}_{-87}$   &  3.87$^{+0.18}_{-0.10}$ &   7200\,$\pm$\,240       \\
    22      &11612274& 0.00            & 7694$^{+214}_{-322}$ &  3.60$^{+0.54}_{-0.06}$ &   7000\,$\pm$\,190      \\
    23      &11718839& 0.03            & 8443$^{+233}_{-367}$ &  3.76$^{+0.42}_{-0.14}$ &   8100\,$\pm$\,500      \\
    24      &11822666& 0.03            & 8615$^{+238}_{-374}$ &  3.78$^{+0.42}_{-0.14}$ &   8200\,$\pm$\,540      \\
\bottomrule
\end{tabular*}
\end{table*}






Synthetic spectra were generated in three steps. First, we computed LTE atmospheric models
using the ATLAS9 code \citep{1993KurCD..13.....K,kur93}. Second, the stellar spectra were then synthesized using SYNTHE \citep{1981SAOSR.391.....K}. 
Third, the spectra were convolved with the instrumental and rotational broadenings.

\teff\, was estimated by computing the ATLAS9 model atmosphere which gave the best match
between the observed $H_{\alpha}$, $H_{\beta}$, $H_{\gamma}$, and $H_{\delta}$ lines profile and those computed with SYNTHE.
The models were computed using solar opacity distribution functions (ODF). It is also known 
that Balmer lines are not sensitive to \logg\, parameter for \teff\,$\lesssim$ 8000\,K \citep{2002A&A...395..601S}. 
Considering the \teff\, range of our targets and the error bars, we fixed the \logg\, and metallicity values to
4.0\,dex and solar, respectively.


The Balmer lines are located far from the edges of the \'{e}chelle orders.
The simultaneous fitting of four lines led to a final solution
at the intersection of the four $\chi^2$ iso-surfaces. An important source of uncertainty arose from 
difficulties in the normalization since it is always challenging for Balmer lines in \'{e}chelle spectra. 
The uncertainties on the \teff\, values were estimated by introducing a 1$\sigma$ change in the 
normalization level and a 0.2\,dex error on \logg\, and metallicity. 
We also considered the  errors on  the initial \vsini\, values (Table~\ref{initialvsini}).
These uncertainties were summed in quadrature with the errors obtained by the
fitting procedure. The final results for \teff\, values and their errors are reported in Table\,\ref{initialvalues}.

\begin{figure*}
\centering
\includegraphics[width=18cm,angle=0]{Fig_2_NOT_paramater_dist_vsini.ps}
\caption{Distributions of the final atmospheric parameters determined from the iron lines analysis and \vsini\, values.}
\label{NOT_paramater_dist}
\end{figure*}



Final atmospheric parameters were derived by using the 
excitation and the ionisation potentials of metal lines. For the correct atmospheric parameters of a star, 
all lines of the same element should give the same chemical abundance. The relationship between 
the chemical abundance and the excitation, ionisation potentials of the 
same element should be flat. 

We used Fe lines in the analysis, since they are the most numerous lines
in the spectra having the \teff\, range of our stars. The ATLAS9 \citep{1993KurCD..13.....K, kur93} model atmospheres were synthesized by using the SYNTHE code \citep{1981SAOSR.391.....K} 
in this and the following chemical abundance analysis.
 A synthetic spectrum is adjusted until it fits well the observed spectrum looking at $\chi^2$ parameter (for more details see
\citeauthor{2015MNRAS.450.2764N} \citeyear{2015MNRAS.450.2764N},
\citeauthor{2016MNRAS.458.2307K} \citeyear{2016MNRAS.458.2307K}).
The Fe lines of all stars were analysed for a range of \teff, \logg, and $\xi$ with a step of 
100\,K, 0.1\,dex and 0.1\,\kms, respectively. The range of the atmospheric parameters was 
selected taking into account the initial \teff\, values derived from Balmer lines 
and the values given by M17 (Table\,\ref{initialvalues}). After the analysis 
was performed, we determined \teff\, and \logg\, values considering the excitation potential$-$abundance and 
the ionisation potential$-$abundance relations, respectively. The $\xi$ values were also obtained by checking the dependence 
between abundance and line strength. The obtained  atmospheric parameters are given in Table\,\ref{Feresult}. The 
uncertainties in the parameters were estimated by checking how much the parameters change for $\sim$$5$\,\% differences in the excitation potential$-$abundance, 
ionisation potential$-$abundance, and the abundance$-$line strength relationships. 


The distributions of the derived atmospheric parameters are shown in Fig.\,\ref{NOT_paramater_dist}. 
The final \teff, \logg, and $\xi$ ranges were obtained to be $7200-8600$\,K, $3.5-4.4$\,dex, and $1.5-4.4$\,\kms, 
respectively. The final atmospheric parameters were compared with the M17 atmospheric 
parameters. In most cases \teff\, values are consistent with each other within errors. 
However, in Fig.~\ref{NOT_FeTeff-logKIC}, the final \logg\, values were also compared 
with the \logg\, values given by M17. The final \logg\, values are 
generally higher than the M17 \logg\, values for the \teff\, ranges of our stars. Additionally, 
it appears that the \logg\, difference between 
our spectroscopic values  and M17 ones increases with the growing \logg. 
The relationship between those parameters is shown in the lower panel of Fig.~\ref{NOT_FeTeff-logKIC}.

 


\begin{figure}
 \centering
 \begin{minipage}[b]{0.43\textwidth}
  \includegraphics[height=3.7cm, width=0.9\textwidth]{Fig_3_1_NOT_FeTeff-logKIC.ps}
  \end{minipage}
 \begin{minipage}[b]{0.43\textwidth}
  \includegraphics[height=3.7cm, width=0.9\textwidth]{Fig_3_2_NOT_logFe-logKIC_fit.ps}
 \end{minipage}
   \caption{Differences between the final spectroscopic \logg\, and the \logg\, taken from M17 as a function
of the final spectroscopic \teff\, (upper panel) and \logg\, (lower panel) obtained from the iron line analysis. Dashed lines represent 1-$\sigma$ levels. 
The Spearman's rank correlation coefficient (R) and probability (the number after comma) are given in the right corner of the right panel.}
 \label{NOT_FeTeff-logKIC}
\end{figure}

Relations between $\xi$ and \teff\, and \logg\, were checked as shown in Fig.\,\ref{NOT_mic-_Teff_log}. 
The $\xi$\,$-$\,\teff\, relationship has already been examined in several studies (e.g.,
 \citealp{2009A&A...503..973L, 2014psce.conf..193G, 2017MNRAS.470.2870N}). 
According to these studies, a decline in $\xi$ is expected for the \teff\, value higher than about 7000\,K. 
The $\xi$\,$-$\,\logg\, relationship was also examined by \citet{2001AJ....121.2159G}. They showed a relation between these 
two parameters for spectral types from  A5 to G2. According to this relation, the  $\xi$ values decrease with increasing \logg. 
However, as our sample spans a  \teff\, range narrower than the A5\,$-$\,G2 one, 
we could not verify these relationships (Fig.\,\ref{NOT_mic-_Teff_log}).

\begin{figure}
\includegraphics[width=8.5cm,angle=0]{Fig4_NOT_mic-_Teff_log.ps}
\caption{The final $\xi$ as a function of the final spectroscopic \teff\, (right panel) and \logg\, (left panel).}
\label{NOT_mic-_Teff_log}
\end{figure}

\begin{table*}
\centering
\footnotesize
  \caption{The final atmospheric parameters derived from the iron lines analysis, the \vsini\, values and the Fe abundances. The final pulsation type classification 
  is given in last column as a result of our frequency analysis. }
  \label{Feresult}
  \begin{tabular*}{0.95\linewidth}{@{\extracolsep{\fill}}lrccccrr@{}}
\toprule 
 Number & KIC   &  \teff            & \logg           & $\xi$            &$\log \epsilon$ (Fe)  & \multicolumn{1}{c}{\vsini}         & \multicolumn{1}{r}{Pulsation}     \\
        &       &      [K]          &  [dex]          &    [\kms]        &[dex]                 &  [\kms]        & \multicolumn{1}{r}{type}\\   
\midrule
    1    &2168333 & 8400\,$\pm$\,200  & 3.7\,$\pm$\,0.2 & 2.0\,$\pm$\,0.5  & 7.47\,$\pm$\,0.27    & 175\,$\pm$\,9  & \dst\,\\
    2    &3119604 & 7700\,$\pm$\,200  & 3.9\,$\pm$\,0.3 & 2.9\,$\pm$\,0.2  & 7.23\,$\pm$\,0.34    & 92\,$\pm$\,5   & \dst\,\\
    3    &3231406 & 7900\,$\pm$\,100  & 3.8\,$\pm$\,0.2 & 2.1\,$\pm$\,0.2  & 7.63\,$\pm$\,0.31    & 169\,$\pm$\,10 & Hybrid\\
    4    &3240556 & 7800\,$\pm$\,200  & 4.4\,$\pm$\,0.2 & 1.5\,$\pm$\,0.5  & 7.51\,$\pm$\,0.35    & 213\,$\pm$\,9  & \dst\,\\
    5    &3245420 & 8100\,$\pm$\,200  & 3.7\,$\pm$\,0.1 & 2.0\,$\pm$\,0.2  & 7.58\,$\pm$\,0.32    & 151\,$\pm$\,6  & \dst\,\\
    6    &3868032 & 8400\,$\pm$\,200  & 4.0\,$\pm$\,0.3 & 3.5\,$\pm$\,0.4  & 7.11\,$\pm$\,0.32    & 181\,$\pm$\,9  & Non-pulsator\\
    7    &4677684 & 8500\,$\pm$\,200  & 3.5\,$\pm$\,0.2 & 1.7\,$\pm$\,0.3  & 7.52\,$\pm$\,0.28    & 71\,$\pm$\,3   & \gd\,\\
    8    &4768677 & 8600\,$\pm$\,300  & 4.2\,$\pm$\,0.3 & 2.0\,$\pm$\,0.5  & 6.94\,$\pm$\,0.50    & 256\,$\pm$\,14 & \dst\,\\
    9    &5180796 & 8400\,$\pm$\,200  & 4.0\,$\pm$\,0.2 & 2.5\,$\pm$\,0.2  & 7.73\,$\pm$\,0.30    & 152\,$\pm$\,6  & \gd\,\\
   10    &5630362 & 7700\,$\pm$\,300  & 3.7\,$\pm$\,0.3 & 3.0\,$\pm$\,0.5  & 7.34\,$\pm$\,0.33    & 230\,$\pm$\,12 & \gd\,\\
   11    &6199731 & 7800\,$\pm$\,200  & 4.1\,$\pm$\,0.3 & 4.4\,$\pm$\,0.4  & 7.23\,$\pm$\,0.35    & 236\,$\pm$\,10 & \dst\,\\
   12    &6500578 & 7700\,$\pm$\,300  & 3.8\,$\pm$\,0.2 & 1.9\,$\pm$\,0.3  & 7.35\,$\pm$\,0.30    & 95\,$\pm$\,6   & \gd\, \\
   13    &6776331 & 7900\,$\pm$\,200  & 3.9\,$\pm$\,0.1 & 2.3\,$\pm$\,0.2  & 7.65\,$\pm$\,0.32    & 49\,$\pm$\,3   & Hybrid\\
   14    &7694191 & 8400\,$\pm$\,200  & 4.1\,$\pm$\,0.2 & 3.9\,$\pm$\,0.2  & 7.31\,$\pm$\,0.36    & 76\,$\pm$\,6   & \gd\,\\
   15    &7732458 & 7800\,$\pm$\,200  & 4.1\,$\pm$\,0.2 & 2.6\,$\pm$\,0.2  & 7.73\,$\pm$\,0.34    & 87\,$\pm$\,5   & \dst\,\\
   16    &9052363 & 8000\,$\pm$\,200  & 4.1\,$\pm$\,0.2 & 3.5\,$\pm$\,0.3  & 7.40\,$\pm$\,0.33    & 109\,$\pm$\,5  & Non-pulsator\\
   17    &9775385 & 7400\,$\pm$\,200  & 4.2\,$\pm$\,0.2 & 4.1\,$\pm$\,0.2  & 7.44\,$\pm$\,0.36    & 71\,$\pm$\,3   & Hybrid\\
   18    &10281360& 7200\,$\pm$\,200  & 4.2\,$\pm$\,0.2 & 4.1\,$\pm$\,0.3  & 7.29\,$\pm$\,0.36    & 108\,$\pm$\,5  & \gd\,\\
   19    &11197934& 7600\,$\pm$\,200  & 3.9\,$\pm$\,0.2 & 3.0\,$\pm$\,0.5  & 7.46\,$\pm$\,0.31    & 267\,$\pm$\,18 & Hybrid\\
   20    &11199412& 7200\,$\pm$\,200  & 4.1\,$\pm$\,0.2 & 3.9\,$\pm$\,0.3  & 6.68\,$\pm$\,0.35    & 77\,$\pm$\,6   & \gd\,\\
   21    &11508397& 7200\,$\pm$\,200  & 3.9\,$\pm$\,0.3 & 1.5\,$\pm$\,0.5  & 7.68\,$\pm$\,0.30    & 240\,$\pm$\,11 & Hybrid\\
   22    &11612274& 7400\,$\pm$\,200  & 3.8\,$\pm$\,0.1 & 3.2\,$\pm$\,0.2  & 7.38\,$\pm$\,0.27    & 130\,$\pm$\,7  & \gd\,\\
   23    &11718839& 8100\,$\pm$\,200  & 3.9\,$\pm$\,0.1 & 3.4\,$\pm$\,0.2  & 7.45\,$\pm$\,0.32    & 57\,$\pm$\,3   & \dst\,\\
   24    &11822666& 8200\,$\pm$\,200  & 4.0\,$\pm$\,0.1 & 2.5\,$\pm$\,0.2  & 7.37\,$\pm$\,0.31    &115\,$\pm$\,7   & Hybrid\\
\bottomrule
\end{tabular*}
\end{table*} 




\section{Analysis of chemical abundances}\label{anche}

Chemical abundances of the stars were derived by performing spectrum synthesis. The derived atmospheric parameters 
were taken as input during the analysis and the chemical abundances of individual elements were determined in addition to \vsini. 

In the first step of the analysis, all lines in the spectra were divided into subsets line by line. In the case of rapidly 
rotating stars, lines are mostly blended. Therefore, for rapidly rotating stars wider ranges in wavelength were selected 
considering the normalisation level. For each subset the line identifications were done by using the line 
list of Kurucz\footnote{kurucz.harvard.edu/linelists.html}. Then these spectral subsets were analysed separately. 
The identified elements in each spectral subset and \vsini\, values were 
adjusted during the analysis. Taking the minimum 
difference between the observed and theoretical spectra, chemical abundances and \vsini\, values 
were obtained from each spectral subset. 
The range of \vsini\,  was found to be from 49 to 267~\kms\, 
(Table~\ref{Feresult}, Fig.\,\ref{NOT_paramater_dist}). 

The average values of the abundances of individual elements for each star 
are given in Table\,\ref{abundances}
and the uncertainties given in this table 
are the standard deviations. The total uncertainties were estimated by considering the 
errors in the obtained atmospheric parameters, assumptions in the model atmospheres, the resolution and the S/N ratio of spectra
\citep{2016MNRAS.458.2307K}.
As a result, the uncertainty in the 
obtained abundances was found to be 0.28\,dex on average. The total uncertainties were estimated for the Fe abundances (Table\,3).

\begin{figure}
\includegraphics[width=8.5cm,angle=0]{Fig5_KIC11199412_NOT_abundance_dis.ps}
\caption{Differences between the derived chemical abundances and the Solar values \citep{2009ARA&A..47..481A} as a function of elements for KIC\,11199412. 
Filled circles 
show the differences that were obtained from at least 5 and more lines, while others represent the opposite.}
\label{KIC11199412_NOT_abundance_dis}
\end{figure}


\begin{figure}
\includegraphics[width=8.5cm,height=3.5cm, angle=0]{Fig6_NOT_FeHcomparison_2figures.ps}
\caption{Differences between the spectroscopic [Fe/H] and [Fe/H] taken from H14 as a function
of the spectroscopic [Fe/H] (left panel) and \teff\, (right panel). The Spearman's rank correlation coefficient (R) and probability (the number after comma) are given 
in the right corner of the right panel.}
\label{NOT_FeHcomparison_2figures}
\end{figure}

Consequently, we found that all stars have chemical 
abundances similar to solar \citep{2009ARA&A..47..481A}. There are three stars, KIC\,3868032, KIC\,4768677 and KIC\,11199412, 
which seem to have moderately underabundant Fe ([Fe/H]\,$\lesssim$\, $-$\,0.50\,dex). However, when the uncertainties in Fe abundances 
of these stars were considered, we can say that only KIC\,11199412 has slightly lower Fe abundance. Additionally, this star displays 
underabundance in almost all elements compared to the solar abundances
(Fig.\,\ref{KIC11199412_NOT_abundance_dis}). 
We compared our [Fe/H] values and the values given by M17 (Fig.\,\ref{NOT_FeHcomparison_2figures}). 
According to the right panel, the difference between the spectroscopic [Fe/H] and the [Fe/H] values given by M17 increases 
after [Fe/H]=0.0\,dex. However, the Spearman's rank coefficient R do not support the significance of such a trend:
the correlation coefficient of $R=0.22$ and a probability value 0.31 indicate that there is no significant correlation.


Finally, KIC\,11199412 clearly shows weak
metal and \ion{Ca}{ii}\,K lines.
KIC\,11718839 also  exhibits
slightly weaker metal and \ion{Ca}{ii}\,K lines comparing the hydrogen spectral type,
but not enough to classify it as a peculiar star.
For sake of completeness, we report our final supervised spectral classification in
 Table~\ref{initialvsini}.




\section{Frequency analysis of photometric data}\label{freana}


To precisely classify the pulsational behaviour of the selected targets, we
performed an independent  frequency analysis of the {\it Kepler} data.
The original time series  were retrieved from the {\it MAST} archive\footnote{http://archive.stsci.edu/}. Keeping the
subdivision  into long- and short-cadence acquisition mode, the original data were normalized to the mean
values of each quarter, thus correcting instrumental drifts. During this procedure isolated outliers were removed.
When reconstructing the pulsational content of our variables, we used the short-cadence time series to investigate
the region above the Nyquist frequency ($f$\,=\,24.5~\cd) of the more numerous long-cadence data.

The accurate time series  were then analysed for their constituent frequencies with the goal to establish their modal content.
We have used the iterative both sine-wave fitting method (Vani\^cek 1971) and the  software package Period04 \citep{2005CoAst.146...53L}.
The final goal is to determine the pulsational characteristics with respect to the spectroscopic properties
and the position in the H-R diagram. The initial variability classification of our targets (\gd, hybrid) was
taken from \citet {2011A&A...534A.125U}. When our analysis was almost finished,  \citet{2019MNRAS.485.2380M} used {\it Gaia}--derived
luminosities to propose another classification scheme based on the skewness of the amplitudes in the Fourier spectra.
However, these authors focused on the \dst\, stars, without investigating  the frequency region below 5~\cd, where \gd\, and
hybrid stars are expected to show their $g$-modes.


\subsection{Variability not induced by pulsation}\label{rotat}

The most obvious case was that of KIC~3868032. The frequency spectrum clearly shows peaks at $f$\,=\,0.40~\cd, 2$f$, 3$f$,
and 4$f$. The spectroscopic analysis pointed out a mean profile with two superimposed components, with clearly
different rotational regimes, i.e., \vsini\,=\,180~\kms\, and \vsini\,$<100$~\kms. 
\citet {2011A&A...534A.125U} classified it as a \gd\, variable.
However, these peaks are associated with the orbital motion and cannot be ascribed to pulsation.
We noticed that after considering $f$ and harmonics,
a peak at $f_1$\,=\,1.68~\cd\, appears, close but not equal to 4$f$. It could be still an artifact of the
orbital/rotational effects, but its pulsational origin  cannot be ruled out.

The case of KIC~9052363 is similar. The light curve is almost flat, without any clear feature that could be ascribed
to pulsation. The frequency analysis reveals a couple of low-frequency, very-low amplitude peaks. Due to
their incoherence, small rotational effects  are the most plausible reason. 
Both stars are also reported as non-\dst\, by \citet{2019MNRAS.485.2380M}.



\begin{figure}
\includegraphics[width=8.5cm]{Fig7_gdorall5.eps}
\caption{Amplitude  spectra obtained by combining short- and long-cadence data of the variables showing 
a regime of low-frequency ($f\le6$~\cd) modes. They are ordered for increasing frequency of the highest peak.}
\label{gdor}
\end{figure}



\subsection {Low-frequency regime}\label{lowf}

Eight stars in our sample show a set of low-frequency peaks only. We ordered them by increasing frequency
of the strongest peak (Fig.~\ref{gdor}). We note that the range of frequencies is small, not exceeding
6~\cd. Occasional high frequency peaks appear in the frequency spectra not enough to claim that a clear double regime of pulsation is present. 
It is noteworthy that all stars in Fig.~\ref{gdor}\, were classified as \gd\, by \citet{2011A&A...534A.125U}.
Actually, they are the only stars classified as \gd\, in our sample and therefore we are in full agreement.

KIC~5630362 is the star showing the highest frequency ($f$\,=\,4.856~\cd, $P$\,=\,0.206~d), with an amplitude much
larger than those of the others (bottom panel).
It is unlikely that this frequency is that of the
fundamental radial mode, since the long period would suggest  an evolved \dst\, star: the very fast rotation
(230~\kms) and the gravity (\logg\,=\,3.7~dex) do not support such an hypothesis. Therefore, pulsators with 
largely predominant modes also exist in \gd\, stars, not only in \dst\, ones (i.e., the
Group A proposed by  \citealt{2019MNRAS.485.2380M}). 






\begin{figure}
\includegraphics[width=8.5cm]{Fig8_dsct6.eps}
\caption{Amplitude spectra obtained by combining short- and long-cadence data of the variables showing a 
prevailing regime of high-frequency modes. They are ordered for increasing frequency of the highest peak.
}
\label{dsct}
\end{figure}

\subsection {High-frequency regime}

None of the stars in our sample was classified as a pure \dst\, star. However, our frequency analysis
showed that some stars are characterized by a large set of high-frequency peaks accompanied by no
(or a few) low-frequency ones (Fig.~\ref{dsct}).
The case of KIC~3240556 is noteworthy, since this star shows a largely predominant mode.
The frequency is too high ($f$\,=\,25.2~\cd) to be due to rotation. Note that  the short-cadence
time series allowed us to determine the exact value of this high frequency very close to the
Nyquist frequency of the long-cadence data.
KIC~3240556 surely belongs to Group A in the \citet{2019MNRAS.485.2380M} classification scheme.


All these variables were classified as hybrid stars. However, 
we know from the analysis
of line-profile variations that very high-degree modes (up to $\ell$=14) are excited in \dst\, stars and they
are spectroscopically detectable \citep{2009A&A...506...85P,2012A&A...542A..24M}. 
The rotational splitting can shift the frequencies of retrograde modes of multiplets
toward low values, thus producing the bunch of peaks observed there in hybrid stars.
They also could be combination terms between high-frequency modes 
or again effects of the rotational modulation induced by changing spots and/or faculae on the stellar 
surfaces.  We emphasize that in general
this multiplicity of causes does not make a few low-frequency
peaks on their own a sufficient criterion to claim for an hybrid pulsational regime.





\subsection{Low and high-frequency regime}

\begin{figure}
\includegraphics[width=8cm]{Fig9_hybrid5.eps}
\caption{Amplitude spectra obtained by combining short- and long-cadence data of the variables showing clear regimes of both high-frequency
and low-frequency peaks.}
\label{hyb}
\end{figure}

Only six stars show simultaneously well-defined low- and high-frequency regimes
(Fig.~\ref{hyb}), with clear peaks in both regions.
However, it has to be noted  how the amplitude spectra differ: there are 
low- and high-frequencies both confined in separated groups and continuously distributed.
All these stars
had an input classification as hybrid variables

\section{Discussion}\label{disc}

We can now use the results from the frequency and atmospheric analyses to better define
the properties of our targets and closely look at how much they are really
peculiar, as suggested by previous works. Our new classification and parameters 
are reported in Table~\ref{Feresult}.
In the left panel of Fig.\,\ref{NOT_hr_filiz_logg}, the log\,\teff\,$-$\logg\, positions are shown  on the basis of
the predicted atmospheric parameters by M17, while in the right panel 
our final atmospheric parameters are used. As can be noticed, the larger changes in the positions are
mostly due to  the differences in the \logg\, values. 

\begin{figure*}
\includegraphics[width=18cm,angle=0]{Fig10_NOT_hr_filiz_logg_2figures.ps}
\caption{ {\it a)} The positions of the target stars, as numbered in Table~\ref{Feresult}, 
according to the atmospheric parameters given by M17.
 {\it b)} The positions of the stars according to the final atmospheric parameters. 
The theoretical instability strips of the 
$\gamma$ Dor (dashed-lines) and $\delta$ Sct (solid lines) stars were taken from \citet{2005A&A...435..927D}. 
The evolutionary tracks (Z=0.02) were 
adopted from \citet{2016MNRAS.458.2307K}.}
\label{NOT_hr_filiz_logg}
\end{figure*}


\subsection{Binary nature} \label{binary}
Composite spectra suggesting double-lined spectroscopic binaries are not observed,
except for the ellipsoidal variable KIC\,3868032 (Sect.~\ref{rotat}).
Seventeen stars in our sample have at least two spectra taken on different nights.
The radial velocities  were examined: if there was a companion, the radial velocities should vary due to the orbital motion.
We did not find any large variation in the radial velocities. Nevertheless,
the very small number of observations and their very limited time coverage do not allow us to clearly detect long-period
or small-amplitude or low-inclination binary systems. 
 For example,
KIC\,2168333 and KIC\,11508397  are known
binaries with very long orbital periods \citep[$\gtrsim$~350~d; ][]{2018MNRAS.474.4322M}, 
but the radial velocities obtained from their two spectra differ by 2.6 and 1.5~\kms\,only, respectively.






\subsection{Pulsation characteristics vs new atmospheric parameters}\label{classif}

Our detailed frequency analysis was able to clean the physical scenario since 
10 targets (8 \dst\, and 2 non-pulsating stars; Sects.~\ref{rotat} and ~\ref{lowf}, respectively) 
out of 24 (42\%) can be retired from the initial sample of \gd\, and hybrid pulsators.
In particular the 8  \dst\, stars erroneously classified as hybrid variables 
turned out to be normal
$p$-modes  pulsators, well inside the \dst\, instability strip. 
Additionally, the two non-pulsators remain located close to the hot border of the \dst\, domain, where 
the fraction of \dst\, pulsating stars is estimated to be 
$\sim$40\% \citep{2019MNRAS.485.2380M}.

Among the six hybrid variables, only  KIC\,11822666 is really a hot star (8200\,K),
while the \teff\,  of KIC\,9775385, KIC\,11197934, and KIC\,11508397 are very close or below the 7300\,K limit when
taking into account the error bars. KIC~6776331 and KIC~3231406 have intermediate values (7900~K).

A similar count applies to the 8 targets re-classified as pure \gd\, stars: five 
(KIC\,5630362, KIC\,6500578, KIC\,10281360, KIC\,11199412, and KIC\,11612274) have
 \teff\, values in agreement with that of the hot border of the classical \gd\, strip (assuming 1$\sigma$ error bars),
 though  KIC\,5630362 and  KIC\,6500578 seem to be more luminous than usual for \gd\, stars (Fig.\,\ref{NOT_hr_filiz_logg_sm}).
Note that the normal \gd\, pulsator KIC\,11199412 is the only star showing hints of chemical peculiarities (Sect.~\ref{chemabund}).
 In three cases (KIC\,4677684, KIC\,7694191, and KIC\,5180796) the \teff\, values are well above 8000\,K.

In the end, we have 8 stars (3 hybrid and 5 normal \gd\, pulsators)  
that show $g$-modes even if they are beyond the hot border of the \gd\, instability strip: namely, 
KIC~3231406, KIC\,6776331, KIC\,11822666, KIC\,4677684, KIC\,5180796, KIC\,5630362, KIC\,6500578, and KIC\,7694191.
They constitute our final bona-fide sample. We recall that the whole initial sample (24 stars)
was expected to be composed of such unusual \gd\, pulsators.

 
\subsection{Other atmospheric parameters}

When the obtained atmospheric parameters are taken into account, most stars are found that they have \teff\, values in agreement with those
given by M17. Only two stars, KIC\,10281360 and KIC\,11199412,  show spectroscopic \teff\, values about 500\,K
cooler than the previously given values by M17.  This explains why the previous classification put them so close to the hot border.

Our final \logg\, values were compared with the \logg\, values given by M17 (Fig.\,\ref {NOT_FeTeff-logKIC}).
In some cases, the final \logg\ values differ more than 0.4\,dex. Therefore, the positions of these stars change
in the log\,\teff\,$-$\logg\, diagram (Fig.\,\ref {NOT_hr_filiz_logg}). When taking into account the final positions
 in the log\,\teff\,$-$\logg\, diagram
and the evolutionary tracks, we can say that $\sim$$70$\,\% of stars have masses from 2.0 to 2.5 $M_{\sun}$, while the mass range of
the others is between 1.6 and 2.0 $M_{\sun}$.

The microturbulent velocity $\xi$ changes with \teff\, 
\citep{2001AJ....121.2159G, 2004IAUS..224..131S, 2009A&A...503..973L} 
and the $\xi$ value varies
from $\sim$1.5 to 3 \kms\, for \teff\, range of about 7000\,$-$\,7300. However, three cool stars
(KIC\,9775385, KIC\,10281360, KIC\,11199412) have $\xi$ values between 3.9 and 4.1~\kms with a maximum uncertainty of 0.3 \kms. These $\xi$ values are high for this \teff\, range, but the three stars
seem to be normally located on the cool border of the \gd\, instability strip. For the other hotter stars
the $\xi$ range is between 1.5 and 4.5 ($\pm$\,0.3)~\kms,  as expected for A-type stars.

Some target stars have low-resolution LAMOST spectroscopy in the literature. The atmospheric parameters derived from these 
studies are given in Table\,\ref{litatm}. When we compared these results with our final atmospheric parameters, we found significant differences 
and trends between the parameters obtained from the high- and low-resolution spectroscopy (Fig.\,\ref{lamostcomp}).


\begin{figure*}
\includegraphics[width=17cm]{Fig11_NOT_lamost_comp.ps}
\caption{Comparison of the atmospheric parameters obtained from high- and low-resolution spectroscopy. 
The subscripts $H$ and $L$ defines the high- and low-resolution spectroscopy, respectively. The Spearman's rank correlation coefficient (R) and probability (the number after comma) are given 
in the right corner of the panels. Blue lines represent the 1$\sigma$ levels.}
\label{lamostcomp}
\end{figure*}


\subsection{Chemical abundances}\label{chemabund}
Anomalous abundances, like the Am phenomenon \citep{2011ApJ...743..153H},
have  been considered a possible physical
explanation of the  hybrid \gd-\dst\, pulsation.
Our analysis did not reveal any abundance peculiarity in 
our bona-fide sample of eight hot \gd\, and hybrid stars,
in agreement with recent studies showing that most hybrid stars are
chemically normal  \citep{2015MNRAS.450.2764N, 2017MNRAS.470.2870N}.
We also searched for He lines, without any significant detection.
Only KIC\,11199412 (a normal \gd\, star) exhibits a moderate underabundance 
in almost all the  elements (Fig.\,\ref{KIC11199412_NOT_abundance_dis}).





\subsection{The \vsini\, values and the SPB hypothesis} \label{spb}
It has been suggested that 
hot \gd\, stars are actually rapidly-rotating SPB stars
\citep{2014A&A...569A..18S,2015MNRAS.451.1445}. 
Due to gravity darkening, their equatorial zones appear cooler than the rest of the surface, 
then they are classified as hot A-stars and, hence,
as hot \gd\, variables. 


\citet{2016MNRAS.460.1318B} determined 
an average \vsini\,=\,114\,\kms\, for a group of six hot \gd\, stars. 
Our five bona-fide hot \gd\, stars  show \vsini\, values ranging from
71 to 230~\kms\, (Table~\ref{Feresult}), with an average value of 125\,\kms.
In the spectral range B5-B9, main-sequence stars, like SPB variables, show a mean \vsini\,=\,144\,\kms\,
\citep[][]{2005ESASP.560..571G, 2016MNRAS.460.1318B}.
Since  hot \gd\, stars do not seem to rotate faster than normal SPB stars, 
 they do not constitute  a special SPB subclass. 
Moreover, we could not find spectral lines indicative of the B-type (Sect.~\ref{chemabund}), which should be present at least for SPB stars seen at intermediate or pole-on orientations. Finally, gravity darkening does not
seem to be able to lower the \teff\, of SPB stars sufficiently to approach the hot border of the $\delta$ Sct instability strip
\citep[see Fig.~4 in ][]{2014A&A...569A..18S}.
We prudently note that  the classification as hot or classical 
hybrid/\gd\, stars could  also be affected by the gravity darkening, especially in case of
fast rotators seen almost equator-on, like  KIC\,5630362, KIC\,11197934, and KIC\,11508397.

\subsection { {\it Gaia} parallaxes} \label{gaia}
As a final check, we used the {\it Gaia} parallaxes \citep{2018A&A...616A...1G} to investigate the positions of the stars
in the H-R diagram. We adopted  the bolometric corrections computed taking into account the \teff\, values
\citep{1996ApJ...469..355F}. The extinction coefficients ($A$$_{V}$) were estimated using the interstellar reddening
$E(B-V)$ values (Table\,\ref{initialvalues}).
These were obtained by measuring the  equivalent widths of the Na\,D lines
in our spectra and then applying the relation given by \cite{1997A&A...318..269M}.






\begin{figure}
\includegraphics[width=8cm,angle=0]{Fig12_hr_logL-Te_NOT_correct_dupret_sm.ps}
\caption{The positions of the stars in the H-R diagram, as numbered in Table~\ref{infotable}.
The theoretical instability strips of the
$\gamma$ Dor (dashed-lines) and $\delta$ Sct (solid black lines) stars were taken from \citet{2005A&A...435..927D}. The recently suggested
\dst\, instability strip \citep{2019MNRAS.485.2380M} is shown by brown solid lines. The evolutionary tracks (Z=0.02) were
adopted from \citet{2016MNRAS.458.2307K}.}
\label{NOT_hr_filiz_logg_sm}
\end{figure}


An offset of $-0.03$\,mas in the {\it Gaia} parallaxes was found \citep{2018A&A...616A...2L, 2019ApJ...878..136Z}. However,
when this offset was applied to a large sample ($\sim$\,15000 stars), a huge amount of stars
goes below the zero-age main-sequence with unreasonably low luminosities \citep{2019MNRAS.485.2380M}.
Therefore, to calculate our final Gaia luminosities we did not apply this offset to the parallaxes, accordingly with recent studies
\citep{2018A&A...616A..17A,2019MNRAS.485.2380M}. The positions of the stars in the H-R diagram are illustrated in Fig.~\ref{NOT_hr_filiz_logg_sm}.
Note that the new empirical  instability strip of \dst\, stars \citep{2019MNRAS.485.2380M} now incorporates
some variables located beyond the hot border of the previous theoretical instability strip.

As can be noticed from Fig.\,\ref{NOT_hr_filiz_logg_sm},
stars  7, 9, 10, 12, and 14 (see Table~\ref{Feresult} for their KIC identification)
have very high luminosities.
When the radii of those hot \gd\, stars were calculated using the
Gaia luminosities, we obtained values ranging from 4.3 to 6.0\,$R_\odot$.
We also calculated the radii of the hot \gd\, stars classified
by \citet{2016MNRAS.460.1318B}. It turned out that $\sim$70\% of their sample have larger radii (3.2\,$<$\,$R_\odot$\,$<$\,7.2).
Main-sequence A-F type  stars should have radii in
the range $1.5-2.7$\,$R_\odot$ \citep{2000asqu.book.....C}.
Therefore, the larger size of the hot \gd\, pulsators appears to  support the hypothesis that they are different from main-sequence A-F pulsators.
They could be evolved A-F stars or more massive stars (such as B-type stars) entering the classical
instability strip in their redward evolution.  
The impact of enhanced iron opacity on stellar pulsations can lead to a substantial
revision of the instability strips \citep{2016MNRAS.455L..67M} and then play a key role to
explain the location of our hot hybrid  and  \gd\, variables. Additionally, the higher luminosity values of these hot stars may be the result of binarity. 
These systems could be a member of
a long orbital period binary system which could not detect with our present
spectroscopic data. 



\section{Conclusions}\label{conclu}

We performed a detailed spectroscopic and photometric study of a group of 24 {\it Kepler} targets,
all claimed to be hot \gd\, or \dst-\gd\, hybrid stars, located well beyond the theoretical 
\gd\, instability strip \citep{2011A&A...534A.125U}.
The detailed frequency analysis of the {\it Kepler} time series and the 
determination of \teff\, and other atmospheric parameters  performed with the new FIES spectra allowed us to set
the bona-fide sample of such peculiar pulsators to five hot \gd\, stars and three hot hybrid stars. 
If on one hand we provide a strong confirmation that these peculiar pulsators exist,
on the other we reduce their recurrence to 1/3 to what originally evaluated. 
Therefore, we can assume that the physical explanation should reside in a mechanism or a cause applicable to 
a limited number of stars, not  to the vast majority of the \dst\, or \gd\, variables.


Searching for this still elusive explanation, 
we did not find any peculiarity in the atmospheric parameters and chemical abundances. 
The lack of composite spectra  does not  support 
the possibility that the hot hybrid stars are binaries in which one or both components are normal \dst\, and/or \gd\, pulsators. 
We investigated the rotation rate and we cannot support the hypothesis that hot \gd\, pulsators 
are actually fast-rotating SPB stars. 
On the other hand, {\it Gaia} luminosities suggested us
that hot \gd\, pulsators have larger radii and higher luminosities than normal main sequence A-F stars.
New  efforts, both theoretical and observational,  have to be
made to well constrain all the features of this new scenario.  



\section*{Acknowledgments}
The authors thank the anonymous referee for useful comments that have improved the presentation of our
results. 
The authors also wish to thank Simon Murphy for his useful comments on a first draft of the manuscript.
Based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical 
Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, under the 
proposal 54-003 (P.I. E.\,Poretti). The authors thank the whole NOT staff for the help in the observations, 
performed in visitor mode.
FKA and GH thank the Polish National Center for Science (NCN) for 
supporting the study through grant 2015/18/A/ST9/00578.
EN acknowledges support provided by the Polish National Science Center
through grant no. 2014/13/B/ST9/00902. 
This work has been partly supported by the Scientific and Technological Research Council of Turkey (TUBITAK) grant numbers 2214-A and 2211-C.
The calculations have been carried out in Wroc{\l}aw Centre for Networking and Supercomputing (http://www.wcss.pl), grant No.\,214. 
This research has made use of the SIMBAD data base, operated at CDS, 
Strasbourq, France.


