\section{Introduction}
\label{sec:sec1}

A movie review is a subjective evaluation written by a moviegoer, usually consisting of two kinds of data: 1) textual reviews and 2) ratings. Textual reviews are qualitative evaluations; ratings are quantitative evaluations expressed by scores. Fig.~\ref{fig:fig1} shows a sample movie review from IMDb\footnote{https://www.imdb.com/}, a leading website managing movie reviews along with other information. The example shows the rating and the textual review for the movie. We note that these two types of information are closely related: the rating may be considered a numerical summarization reflecting the textual review.

\begin{figure}[!h]
\centering {
  \includegraphics[width=13cm]{img/fig1.png}
}
\caption{Sample movie review}
\label{fig:fig1}
\end{figure}

Movie reviews have a significant impact on the decision making of potential moviegoers~\cite{chintagunta2010effects}. However, a certain group of movie reviewers may generate indistinguishable and/or intentionally malicious reviews so as to bias the overall movie rating. While such groups ought to lack credibility, many well-known existing movie review management systems including Metacritic~\footnote{https://www.metacritic.com/} and Rotten Tomatoes~\footnote{https://www.rottentomatoes.com} do not consider the credibility of movie reviews in calculating overall movie ratings.

To enhance the evaluation of ratings in movie review systems, it is necessary to make the credibility of movie reviews a factor in the overall movie rating. For this, we need a method to classify movie reviews according to their credibility. However, this is challenging because even experts cannot clearly judge movie review credibility and the number of movie reviews is very large, which requires much time and effort for manual annotation. Table~\ref{tab:tab6} shows an example of different annotation results for a real review from Rotten Tomatoes\footnote{https://www.rottentomatoes.com/}.

\begin{table}
\center{
\caption{An example of varying annotation results}
\setlength{\tabcolsep}{3pt}
\begin{tabular}{|c|c|c|c|}
\hline
\textbf{Textual reviews} & \textbf{Ratings}  & \textbf{Annotation1}  & \textbf{Annotation2}  \\ \hline\hline
\makecell{Really bad. I'm a huge fan of the Marvel \\movies in general, but this one stands out \\ as a big disappointment.} & 4 & Distrusted & Trusted \\ \hline

\end{tabular}

\label{tab:tab6}
}
\end{table}

In this paper, we propose a weakly supervised learning method for judging the credibility of movie reviews. It annotates entire data sets based on a rule-based function instead of manual annotation by human beings, thus annotating at a very fast speed. The most influential factor in the accuracy of weakly supervised learning is a criterion for classification. In this paper, we use a simple and clear criterion based on the historical movie ratings of movie reviewers. Fig.~\ref{fig:fig2} shows actual examples of movie reviewers submitting indistinguishable reviews for multiple movies, which are consequently considered untrustworthy. Fig.~\ref{fig:fig2-a} is an example where a movie reviewer has given the maximum 10 points to all the movies they have reviewed. We indicate that indistinguishable overall positive reviews with similar expressions are written regardless of the specific movies reviewed. Fig.~\ref{fig:fig2-b} shows another example where a movie reviewer has given the minimum 0 points and overall negative reviews to all the movies reviewed. Hence, we conclude by evaluating their historical reviews that reviews that will be written by those movie reviewers also have no discrimination.

\begin{figure}[!h]
\begin{subfigure}[t]{0.47\linewidth}
  \centering

  \includegraphics[width=.99\textwidth]{img/fig2-a.png}
  \caption{A reviewer giving the maximum points}
  \label{fig:fig2-a}
\end{subfigure}
\begin{subfigure}[t]{0.47\linewidth}
  \centering

  \includegraphics[width=.99\textwidth]{img/fig2-b.png}
  \caption{A reviewer giving the minimum points}
  \label{fig:fig2-b}
\end{subfigure}%
\caption{Examples of movie reviewers without discrimination}
\label{fig:fig2}
\end{figure}

The contributions of the paper is summarized as follows.

\begin{itemize}
    \item We propose a weakly learning method for judging the credibility of movie reviews. The proposed method enables fast annotation according to predefined rules. To reduce accuracy limitations in weakly supervised learning, we use a simple but effective criterion to classify movie reviews clearly according to their credibility, which we call \emph{historical credibility}.
    \item The proposed learning method is \emph{efficient}. Specifically, we show that the proposed method can annotate 8,000 movie reviews in only 0.712 seconds, which occupies about 0.145\% of the entire training process when we use Word2Vec as the text representation model. This is a notable result if the annotation results are effective because the usual (human) annotation approach to classifying movie reviews based on credibility takes much time and effort even for experts.
    \item The proposed learning model is \emph{effective}. Through extensive experiments, we measure the accuracy of the proposed learning model and that of the chosen comparison method on 40,000 movie reviews collected for five different movies. The comparison method uses the helpfulness votes of other reviewers as the criterion to judge the credibility of movie reviews, a criterion that has been used in much research. Our finding is that the proposed learning model shows an improvement in accuracy of the learning model over the comparison method by 1.57\% $\sim$ 4.54\% depending on the text representation model, machine learning technique, and movie genres chosen.
\end{itemize}

The paper is organized as follows. In Section~\ref{sec:sec2}, we describe related work. In Section~\ref{sec:sec3}, we describe text representation models and machine learning techniques that will be used in the proposed learning method. In Section~\ref{sec:sec4}, we propose a weakly supervised learning method using historical credibility as a criterion for judging the credibility of movie reviews. In Section~\ref{sec:sec5}, we explain how the experimental results show the efficiency and accuracy of the proposed method. In Section~\ref{sec:sec6}, we further discuss the proposed learning method. Section~\ref{sec:sec7} concludes the paper.

\section{Related Work}
\label{sec:sec2}

\noindent {\bf Research on weakly supervised learning:} Weakly supervised learning allows us to annotate large-scale data at a very fast rate, at the cost of potentially incurring some inaccurate annotation. Typical annotation methods for weakly supervised learning can be classified into three categories~\cite{dong2018data}: 1) crowdsourcing, 2) distant supervision, and 3) labeling functions. Crowdsourcing is a form of supervised learning, but annotation of data is performed by several annotators who are not experts~\cite{raykar2010learning}. Distant supervision extracts the structure from a data set using large-scale existing databases; the extracted structure is used for annotation~\cite{mintz2009distant}. Labeling functions are constructed from defined rules for annotating data~\cite{ratner2017snorkel, ratner2016data}. This allows annotation to be done quickly, but the results include some noise. The weakly supervised learning used in this paper employs labeling functions because we annotate the entire set of movie reviews based on the predefined criterion, which is obtained from analysis of historical movie ratings.

Recently, weakly supervised learning has been adopted in many studies. For instance, Lee et al. have applied weakly supervised learning~\cite{lee2018sentiment} to sentiment classification, using weakly supervised learning based on a convolutional neural network to identify keywords for classifying positive and negative sentences. Lin et al. have used weakly supervised learning for understanding the sentiments of users from social media content containing different data models such as texts and images~\cite{lin2018multi}.

\noindent {\bf Research on movie reviews:} For movie reviews, most research has focused on sentiment analysis of textual reviews, because movie reviews depend on the sentiment experienced after watching the movie~\cite{serrano2015sentiment}. Topal et al. have constructed an emotion map according to emotions expressed after watching the movie and have recommended movies to moviegoers based on this emotion map~\cite{topal2016movie}. Manek et al. have performed sentiment classification for movie reviews using a Gini index based on feature selection and a support vector machine (SVM) classifier~\cite{manek2017aspect}. Chakraborty  et al. have compared the performance of two clustering algorithms, namely k-means and k-means++, for emotion classification of movie reviews by using Word2Vec as a classifier of the textual review~\cite{chakraborty2018comparative}. Alsaqer et al. have improved the accuracy of movie review summarization based on sentiment analysis~\cite{alsaqer2017movie}. He et al. have proposed a self-training method using some labeled features extracted from an existing sentiment lexicon for sentiment analysis~\cite{he2011self}. Elmurngi et al. have detected fake movie reviews based on sentiment analysis using supervised learning methods such as SVM and naïve Bayes~\cite{elmurngi2018fake}.

\noindent {\bf Research on credibility in online review systems:} To the best of our knowledge, there are no public methods that consider the credibility of movie reviews. Many well-known existing movie scoring systems including Metacritic and Rotten Tomatoes do not consider movie review credibility in calculating the overall movie rating. IMDb is known to consider movie reviewers’ credibility; however, since the specific methods are not disclosed, we cannot evaluate their effectiveness.

On the other hand, much effort has been invested in research on review credibility in other online review systems. It is well known that judging the credibility of online reviews is important because they affect consumers’ decision making~\cite{lee2011helpful, mudambi2010makes, ham2019subjective}. Liu et al. have evaluated the effectiveness of reviews of travel products based on two aspects~\cite{liu2015makes}: 1) reviewer information such as profiles and 2) reviews including some quantitative information (i.e., ratings and length of reviews). As a result, they have shown that qualitative judgments are the most influential factors for review credibility. Ghose et al. have explored multiple aspects of textual reviews, namely subjectivity levels, readability, and spelling errors, to understand the helpfulness of reviews in terms of economic and social outcomes~\cite{ghose2010estimating}. Hochmeister et al. have compared destination experts (i.e., the most active members) in TripAdvisor with general reviewers and have shown that destination experts have more influence in the TripAdvisor community~\cite{hochmeister2013destination}. Barbado et al. have proposed a method for detecting fake online reviews in the domain of consumer electronics businesses (i.e., studying Yelp businesses)~\cite{barbado2019framework}, labeling reviews in the data set used as trustworthy or fake. Reyes et al. have studied five chosen factors in online reviews that can affect the decision making of customers by exploring TripAdvisor reviews~\cite{reyes2019impact}.

Some research results have shown that helpfulness votes from other reviewers are a critical measure for evaluating review credibility. Fang et al. have conducted empirical analysis to explore two factors affecting the value of reviews: 1) the reviews themselves and 2) reviewer characteristics~\cite{fang2016analysis}, using the helpfulness vote as the criterion for determining review value. Shan et al. have shown that the reviewer reputation generated by peer ratings has great influence on credibility evaluation~\cite{shan2016credible}. Yang et al. have examined six factors (specifically, reviewer location, reviewer level, reviewer helpfulness vote, review rating, review length, and review photo) and have shown that the review rating and reviewer helpfulness vote are the most influential factors~\cite{yang2017exploring}. Gang et al. have examined the relationship between emotions and review helpfulness~\cite{ren2019examining}. For this, they have defined three kinds of emotion embedded in online reviews (specifically, anger, fear, and sadness) and have analyzed the different effects of each emotion on perceived review helpfulness.

\section{Preliminaries}
\label{sec:sec3}

\subsection{Text representation models}
\label{sec:sec3-1}

In this paper, we build a learning model to classify movie reviews based on their credibility. To annotate movie reviews based on textual content, we need text representation models. We use two kinds of bag-of-words models: 1) TF-IDF and 2) Word2Vec. TF-IDF is a representative model for representing the text~\cite{jones2004statistical}. Word2Vec is a vector space model that embeds a corpus of text into a vector space~\cite{mikolov2013distributed}. The purpose of this paper is not to propose a specific novel text representation model. Instead, we choose some representative text representation models to show the effectiveness of the proposed learning method; we could equally apply any other models in place of TF-IDF and Word2Vec to the proposed method.

\subsubsection{TF-IDF}

TF-IDF represents the relative importance of words in a document based on the entire set of documents~\cite{wu2008interpreting}. Eq.~(\ref{equ1-1}) is the formula for obtaining TF-IDF, which is calculated by multiplying TF (text frequency) and IDF (inverse document frequency). Essentially, the result of TF-IDF increases as the frequency of given words in a document increases and the number of documents containing given words decreases. Eq.~(\ref{equ1-2}) and Eq.~(\ref{equ1-3}) are formulas for obtaining TF and IDF, respectively. TF is the frequency of a particular word in the document. Specifically, $f(i,j)$ is the frequency of a word $i$ in a document $j$ while $\max  \left\{ f(w,d):w\in d \right\}$ is the maximum of $f(i,j)$ of the words in the document. IDF is the total number of documents $|D|$ divided by the number of documents containing a word $t$. IDF uses a logarithm to prevent the result from increasing too much as $|D|$ increases~\cite{chowdhury2010introduction}.

\begin{eqnarray}
\mathrm{TF-IDF_{i,j}} = TF_{i,j} \times IDF_{i}
\label{equ1-1}
\end{eqnarray}

\begin{eqnarray}
\mathrm{TF_{i, j}} = 0.5 + \frac{0.5 \times f(i, j)}{\max  \left\{ f(w,d):w\in d \right\} }
\label{equ1-2}
\end{eqnarray}

\begin{eqnarray}
\mathrm{IDF_{i}} =  \log\left( \frac{\mid D\mid}{count \left\{ d:t\in d  \right\}   }\right)
\label{equ1-3}
\end{eqnarray}

\subsubsection{Word2Vec}

Word2Vec is a word embedding model that maps words contained in the document into a vector space, where it tries to preserve contextual meanings~\cite{mikolov2013distributed}. Word2Vec supports two kinds of architectures~\cite{ling2015two}: 1) continuous bag-of-words~(simply, CBOW) and 2) continuous skip-gram, as shown in Fig.~\ref{fig:word2vec}. The continuous CBOW architecture tries to predict the center word from the context words. For example, in a given sentence, ``This movie \_ is wonderful,'' this architecture tries to find the most likely word to fill in the blank from the context. As a result, it shows better performance when we predict frequently used words. Conversely, the continuous skip-gram architecture tries to predict the context from a given center word. For example, for the given center word ``soundtrack,'' this architecture tries to predict the surrounding context.

\begin{figure}[!h]
\begin{subfigure}[b]{.47\linewidth}
  \centering
  \includegraphics[width=0.95\linewidth]{img/cbow.png}
  \caption{Continuous bag-of-words~(CBOW)}
  \label{fig:cbox}
\end{subfigure}
\begin{subfigure}[b]{.48\linewidth}
  \centering
  \includegraphics[width=0.97\linewidth]{img/skipgram.png}
  \caption{Continuous skip-gram}
  \label{fig:skipgram}
\end{subfigure}
\caption{Two architectures for Word2Vec~\cite{ling2015two}}
\label{fig:word2vec}
\end{figure}

\subsection{Machine learning techniques}
\label{sec:sec3-2}

For the classification of movie reviews, we apply machine learning techniques to the textual reviews. A variety of text classifiers using machine learning techniques have been proposed~\cite{ikonomakis2005text}: decision trees, naïve Bayes, rule induction, neural networks, nearest neighbors, and support vector machines~(SVMs). In this paper, we apply two techniques that have been studied for movie reviews~\cite{pang2002thumbs}: 1) naïve Bayes and 2) SVM. (As previously noted, we could apply any other machine learning techniques instead of naïve Bayes and SVM for the proposed method.)

\subsubsection{Naïve Bayes classifier}

The naïve Bayes classifier has been widely used for document classification because it is quite simple but shows high accuracy~\cite{ch1997bayesian}. The name ``naïve Bayes’’ refers to a classifier based on Bayes' theorem. The important characteristic of the naïve Bayes classifier is that the features used for the classifier are stochastically independent of each other.

When a set of classes $C$ = \{${\mathbf c}_1$, ${\mathbf c}_2$, ..., ${\mathbf c}_n$\} is given, the naïve Bayes classifier predicts the ${\mathbf c}_i$ that is most likely to include a given document $d$ based on Bayes' theorem. Eq.~(\ref{equl-4}) shows Bayes' theorem, which gives the probability that a given document $d$ is included in a category ${\mathbf c}_i$. $P(d)$ is the probability of choosing a document $d$ from all the documents. $P$(${\mathbf c}_i$) is the number of documents in a category ${\mathbf c}_i$ divided by the total number of documents. $P$($d$ $\mid$ ${\mathbf c}_i$) is, for a given category ${\mathbf c}_i$, the probability of choosing document $d$ in that category. When a set of classes \{${\mathbf c}_1$, ${\mathbf c}_2$\} is given, we can predict the category in which a document $d$ is included by comparing $P$(${\mathbf c}_1\mid d$) and $P$(${\mathbf c}_2\mid d$).

\begin{eqnarray}
\mathrm{P(c_{i}\mid d)}=\frac{P(c_{i} ,d)}{P(d)} =  \frac{P(c_{i})P(d\mid c_{i})}{P(d)}
\label{equl-4}
\end{eqnarray}

\subsubsection{Support vector machines}

A support vector machine~(SVM) finds the optimal hyperplane for a linear classification on data on a multidimensional space~\cite{vapnik2013nature}. Here, we can have multiple possible hyperplanes that can classify the data. The SVM finds a hyperplane having the maximal margin. Fig.~\ref{fig:fig3} shows an example of SVM-based classification. Here, when we compare candidate hyperplanes ${\mathbf B}_1$ and ${\mathbf B}_2$ in terms of the margin between the classified categories, ${\mathbf B}_1$ is chosen as the result hyperplane because its margin is greater than that of ${\mathbf B}_2$. The maximal margins concept aims to minimize incorrect classifications.

\begin{figure}[!h]
\centering{
\includegraphics[width=0.5\linewidth]{img/fig3.png}
}
\caption{An example of the SVM classifier}
\label{fig:fig3}
\end{figure}

We can extend SVM to the case of nonlinear classification. For this, the kernel function has been introduced~\cite{wilk2016comparative}, defining additional functions for describing new dimensions required in nonlinear classification. Fig.~\ref{fig:fig4} shows the effect of nonlinear classification using the SVM kernel function. Fig.~\ref{fig:fig4-a} shows objects that are not classified by linear classification in SVM; Fig.~\ref{fig:fig4-b} shows that they can be classified by nonlinear classification using a kernel function. Kernel functions commonly used for the SVM classifier are the polynomial kernel, the Gaussian kernel, and the radial basis function~\cite{wilk2016comparative}.

\begin{figure}[!h]
\begin{subfigure}[b]{.48\linewidth}
  \centering
  \includegraphics[width=0.95\linewidth]{img/fig4_a.png}
  \caption{Linear classification}
  \label{fig:fig4-a}
\end{subfigure}
\begin{subfigure}[b]{.47\linewidth}
  \centering
  \includegraphics[width=0.95\linewidth]{img/fig4_b.png}
  \caption{Nonlinear classification}
  \label{fig:fig4-b}
\end{subfigure}
\caption{The effect of nonlinear classification using a kernel function in SVM}
\label{fig:fig4}
\end{figure}

\section{The Proposed Method}
\label{sec:sec4}

In this section, we propose a new learning method for judging the credibility of movie reviews. Our research objectives are summarized as follows: the proposed learning model 1) should be efficient at dealing with large-scale movie review data sets and 2) should provide accurate results based on a clear criterion for judging the credibility of movie reviews. To satisfy the first objective, we propose a method based on weakly supervised learning, which allows fast annotation by a predefined rule. To satisfy the second objective, we propose a new criterion for judging the credibility of movie reviews, based essentially on historical rating records. In Section~\ref{sec:sec4-1}, we will present the overall framework based on weakly supervised learning. In Section~\ref{sec:sec4-2}, we will formally define a new criterion for judging credibility.

\subsection{The overall framework}
\label{sec:sec4-1}

Fig.~\ref{fig:fig5} shows the overall framework for classifying movie reviews based on credibility. It consists of 1) a learning phase and 2) an evaluation phase. In the learning phase, we first analyze historical review ratings for each reviewer so as to determine movie review credibility. In Section~\ref{sec:sec4-2}, we present a new criterion for determining movie review credibility. Next, we annotate the entire movie review as ``trusted'' or ``distrusted'' based on the proposed criterion. Here, we note that the annotation of the entire data set is performed at a very fast speed because we can annotate movie reviews based on the defined rule using a clear criterion without manual effort. Finally, we build a learning model by obtaining a vector of keywords from textual reviews using two text representation models, namely Word2Vec and TF-IDF, and by applying two machine learning techniques, namely naïve Bayes classifier and SVM, to the annotated data sets. In the evaluation phase, we can use the learning model to judge the credibility of new movie reviews based on the textual review information. Here, thanks to the learning model, we can even judge the credibility of movie reviews that are written by first-time reviewers.

\begin{figure}[!h]
\begin{subfigure}[b]{1\linewidth}
  \centering
  \includegraphics[width=0.97\linewidth]{img/fig5-a.png}
  \caption{Learning phase}
  \label{fig:fig5-a}
\end{subfigure}
\vspace*{0.3cm}
\begin{subfigure}[b]{1\linewidth}
  \centering
  \includegraphics[width=0.97\linewidth]{img/fig5-b.png}
  \caption{Evaluation phase}
  \label{fig:fig5-b}
\end{subfigure}
\caption{A weakly supervised learning method for classifying movie reviews based on credibility}
\label{fig:fig5}
\end{figure}

\subsection{Historical credibility for weakly supervised learning}
\label{sec:sec4-2}

In this section we present a simple criterion, which we call \emph{historical credibility}, for weakly supervised learning. We need a criterion that can clearly classify movie reviews according to their credibility, utilizing the historical movie ratings of each movie reviewer. The basic idea is to identify movie reviewers who have given indistinguishable ratings and to define all movie reviews written by them as distrusted reviews. We define historical credibility as in Definition~1.

\begin{definition}
\emph{Historical credibility} is a criterion that evaluates review credibility based on the historical review ratings for each reviewer, defined as follows: all reviews are classified as either trusted or distrusted. The \emph{distrusted reviews} are defined as those written by a reviewer $R$ satisfying the following two conditions: 1) the standard deviation of all past movie ratings $X_R$ is 0 and 2) the average over $X_R$ is either the lowest score (e.g., 1 point) or the highest score (e.g., 10 points). All other reviews are defined as the \emph{trusted reviews}. Here, we exclude a single review written by a reviewer because it is classified as distrusted by the definition but cannot be judged.\hfill
\end{definition}

As a comparison criterion for the proposed historical credibility, we choose the helpfulness vote. In general review systems, the helpfulness vote has been used as a meaningful measure for determining the value of the review, as explained in Section~\ref{sec:sec2}. Thus, Shan et al. and Yang et al. have shown that the helpfulness vote has a great influence on credibility evaluation~\cite{shan2016credible, yang2017exploring}. Further, Fang et al. have used the helpfulness vote as a criterion to evaluate the influence of various review factors on review credibility~\cite{fang2016analysis}. Similarly, with movie reviews, the helpfulness vote can be an effective criterion in determining review value.

Table~\ref{tab:tab1} shows annotation methods for historical credibility as the proposed criterion and the helpfulness vote as the comparison criterion. For historical credibility, we define distrusted reviews by Definition~1 and all other reviews as trusted. For the helpfulness vote, we define each movie review as trusted when the number of helpful votes for the review exceeds the number of unhelpful votes.

\begin{table}
\fontsize{9}{11}{11}\selectfont
\centering
\caption{Data annotation based on credibility criteria
}
\label{tab:tab1}
\begin{tabular}{|c|p{4cm}|p{4cm}|}
\hline
\textbf{Movie reviews}
& \textbf{Historical credibility}
& \textbf{Helpfulness vote~\cite{shan2016credible, yang2017exploring, fang2016analysis}}

\\ \hline \hline
Distrusted review
& Review that is defined as distrusted by Definition~1
& Review with fewer helpful than unhelpful votes, or equal numbers of each kind of vote
\\ \hline
Trusted review
& Review that is defined as trusted by Definition~1
& Review with more helpful than unhelpful votes
\\ \hline
\end{tabular}
\end{table}

\section{Performance Evaluation}
\label{sec:sec5}

\subsection{Data collection}
\label{sec:sec5-1}

In this section, we describe the method we used to collect movie reviews and present the results of the data collection. We use BeautifulSoup and Selenium to implement a crawler to collect movie reviews from a movie website. The target movie website is Naver Movie\footnote{https://movie.naver.com/.}, the largest movie website in Korea\footnote{As of October 17, 2019, Naver was managing 12,539,719 ratings and 2,086,503 textual reviews.}. From this website, we can collect all the historical movie reviews for a specific user; this is an essential condition for verifying the concept proposed in this paper. In addition, we collect the number of helpful or unhelpful votes for each movie review, as required for the comparison method. Only if these conditions are satisfied can we apply our concept into data sets crawled from any other website carrying movie reviews.

We consider various movie genres, comparing learning for movie reviews within the same genre and across different genres. For this purpose, reviews of a total of five movies from three different genres were collected. We collect movie reviews having both movie ratings and textual reviews. The movie ratings range from 1 through 10 points; the textual reviews are described in up to 140 characters. We also collect all the historical movie ratings for the movie reviewers included for the five movies chosen. The total number of historical movie ratings for these movie reviewers was 87,125. Table~\ref{tab:table2} shows the characteristics of the movie reviews collected. We collected 40,000 movie reviews, 8,000 for each movie. Here, 32,000 reviews are used for training and 8,000 are used for testing.

\begin{table}[]
\centering
\caption{Characteristics of collected movie reviews}
\label{tab:table2}
\begin{tabular}{|c|c|c|c|c|}
\hline
\textbf{Genres}        & \textbf{Movies}                                               & \textbf{Release date} & \textbf{\begin{tabular}[c]{@{}c@{}}Training \\ data\end{tabular}} & \textbf{\begin{tabular}[c]{@{}c@{}}Testing\\ data\end{tabular}} \\ \hline \hline
\multirow{3}{*}{Drama} & Assassination                                                          & 22 Jul. 2015          & 6,400                                                                      & 1,600                                                                     \\ \cline{2-5}
                       & \begin{tabular}[c]{@{}c@{}}I can speak\end{tabular} & 21 Sep. 2017          & 6,400                                                                      & 1,600                                                                     \\ \cline{2-5}
                       & The Spy Gone North                                                 & 08 Aug. 2018          & 6,400                                                                      & 1,600                                                                     \\ \hline
Comedy                 & \begin{tabular}[c]{@{}c@{}}Intimate Strangers\end{tabular}                                           & 31 Oct. 2018          & 6,400                                                                      & 1,600                                                                     \\ \hline
Action                 & Confidential Assignment                                                   & 18 Jan. 2017          & 6,400                                                                      & 1,600                                                                     \\ \hline
\hline
\multicolumn{3}{|c|}{Total}                                                                                & 32,000                                                                      & 8,000                                                                    \\ \hline
\end{tabular}
\end{table}

\subsection{Experimental methods}
\label{sec:sec5-2}

In the experiments, we measure the elapsed learning time for movie reviews to check the efficiency of annotation in the proposed learning method. In addition, we measure the accuracy of the proposed learning method using historical credibility and that of the comparison method using the helpfulness vote as the classification criterion. Here, we measure the accuracy of combinations of the text representation models (TF-IDF and Word2Vec) and machine learning techniques (naïve Bayes and SVM). For TF-IDF and Word2Vec, we first need to extract keywords from the movie reviews. For this, we use a Korean keyword extractor \footnote{https://pypi.org/project/python-mecab-ko/} and select the top 20 keywords based on frequency for each movie review. In TF-IDF, we apply a logarithm to TF to smooth its significance because the same proper nouns, such as actor/actress and movie names, tend to be repeated in movie reviews. In Word2Vec, we use a skip-gram to consider the relevance of given keywords such as actor/actress and movie names in the movie reviews. We use a Gaussian kernel function for SVM, which has been widely used to classify textual data on high dimensions~\cite{keerthi2003asymptotic}. For the experiment, the entire data set is annotated according to the annotation methods presented in Section~4.2.

The experimental results consist of three parts:

\begin{enumerate}
    \item Elapsed time for the proposed method
    \item Accuracy comparison between the proposed method and the helpfulness vote
    \item Variation in accuracy of the proposed method
\end{enumerate}

In the first part, we show the efficiency of the proposed learning method. In the second part, we use the following experimental variables: 1) textual representation models, 2) machine learning techniques, and 3) variety of movie genres. In the third part, we use the following experimental variables: 1) number of movies, 2) machine learning techniques, and 3) textual representation models.

\subsection{Experimental results}
\label{sec:sec5-3}

\subsubsection*{Elapsed time for the proposed method}

Table~\ref{tab:tab4} shows the average elapsed learning time for the five chosen movies using the proposed learning method using a combination of text representation models and machine learning techniques. Here, to measure the elapsed time, we use 8,000 movie reviews for each movie. We note that annotation of 8,000 movie reviews requires only 0.712 seconds on average. In particular, we note that when we use Word2Vec as the text representation model, which shows better performance than TF-IDF as will be shown in Section 5.3.3, the portion of the entire training process occupied by annotation is only about 0.145\% $\sim$ 0.280\%. By comparison, manual annotation to judge the credibility of 8,000 movie reviews requires much more time in absolute terms; the helpfulness vote approach also requires much time to collect enough votes after the movie reviews are written.

\begin{table}
\fontsize{9}{10}\selectfont
\center{
\caption{ Elapsed learning time for movie reviews (seconds)}
\setlength{\tabcolsep}{3pt}
\begin{tabular}{|c|c|c|c|c|}

\hline
Combination & \textbf{TF-IDF + NB} & \textbf{W2V + NB}  & \textbf{TF-IDF + SVM} & \textbf{W2V + SVM} \\
\hline
\hline Annotation & \multicolumn{4}{|c|}{0.712}  \\
\hline Training & 0.122 & 195.822 & 3.533 & 288.114 \\
\hline Testing & 0.040 & 57.615 & 0.027 & 200.070 \\
\hline
\hline  \textbf{Total} & 0.874 & 254.149 & 4.273 & 488.896 \\

\hline
\end{tabular}
\label{tab:tab4}
}
\end{table}

\subsubsection{Accuracy comparison between the proposed method and the helpfulness vote}

Fig.~\ref{fig:fig6} compares the accuracy of the methods by credibility criteria, namely the proposed historical credibility and the helpfulness vote. Here, we measure the accuracy for each combination of text representation models, namely TF-IDF or Word2Vec, and machine learning techniques, namely naïve Bayes or SVM. The results indicate that the accuracy of the proposed historical credibility is comparable to or even better than that of the helpfulness vote for all combinations of text representation models and machine learning techniques.

Fig.~\ref{fig:fig6-a} shows the accuracy when we use TF-IDF as the text representation model and naïve Bayes as the machine learning technique. The result shows that mean accuracy for the helpfulness vote is 0.506, while that for historical credibility is 0.526. That is, the proposed historical credibility shows better accuracy than the helpfulness vote by a margin of 3.92\%. Fig.~\ref{fig:fig6-b} shows the accuracy when we use TF-IDF with SVM as the machine learning technique. The result shows that the proposed historical credibility has an accuracy improved over the helpfulness vote by 4.41\%. Fig.~\ref{fig:fig6-c} shows the accuracy when we use Word2Vec as the text representation model and naïve Bayes as the machine learning technique. The result shows that the proposed historical credibility has an accuracy improved over the helpfulness vote by 1.55\%. Fig.~\ref{fig:fig6-d} shows that the accuracy of the proposed historical credibility when we use Word2Vec and SVM is improved over the accuracy of the helpfulness vote by 2.49\%.

\begin{figure}[!h]
  \begin{subfigure}[b]{0.5\linewidth}
    \centering
    \includegraphics[width=0.95\linewidth]{img/fig6_a.png}
    \caption{Naïve Bayes with TF-IDF}
    \label{fig:fig6-a}
    \vspace{4ex}
  \end{subfigure
  \begin{subfigure}[b]{0.5\linewidth}
    \centering
    \includegraphics[width=0.95\linewidth]{img/fig6_b.png}
    \caption{SVM with TF-IDF}
    \label{fig:fig6-b}
    \vspace{4ex}
  \end{subfigure}
  \begin{subfigure}[b]{0.5\linewidth}
    \centering
    \includegraphics[width=0.95\linewidth]{img/fig6_c.png}
    \caption{Naïve Bayes with Word2Vec}
    \label{fig:fig6-c}
  \end{subfigure
  \begin{subfigure}[b]{0.5\linewidth}
    \centering
    \includegraphics[width=0.95\linewidth]{img/fig6_d.png}
    \caption{SVM with Word2Vec}
    \label{fig:fig6-d}
  \end{subfigure}
  \caption{Accuracy comparison by combination of text representation models and machine learning techniques}
  \label{fig:fig6}
\end{figure}

Fig.~\ref{fig:fig9} compares the accuracy of the methods by credibility criteria, namely historical credibility and helpfulness vote, and by variety of movie genres. Here, we use three movies having the same genre, drama, and then having different genres: drama, comedy, and action. The result shows that in improvement in accuracy of the proposed historical credibility over that of the helpfulness vote in both cases. Specifically, the proposed historical credibility shows accuracy improvements of 1.57\% in the case of the same genre and 4.54\% in the case of the three different genres.

\begin{figure}[!h]
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.99\linewidth]{img/fig9_a.png}
  \caption{The same movie genre: drama}
  \label{fig:fig9-a}
\end{subfigure}
\quad
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.99\linewidth]{img/fig9_b.png}
  \caption{ Different movie genres: drama, comedy, and action}
  \label{fig:fig9-b}
\end{subfigure}

\caption{Accuracy comparison by variety of movie genres}
\label{fig:fig9}
\end{figure}

\subsubsection*{Variation in accuracy of the proposed method}

Fig.~\ref{fig:fig7} shows the accuracy of the proposed historical credibility by text representation model, namely TF-IDF and Word2Vec. Here, we select three movies from different genres: Assassination, Intimate Strangers, and Confidential Assignment. Fig.~\ref{fig:fig7-a} shows the accuracy of TF-IDF and that of Word2Vec when we keep the machine learning technique fixed as naïve Bayes; Fig.~\ref{fig:fig7-b} shows the respective accuracy with the machine learning technique fixed as SVM. Both results show better accuracy when using Word2Vec than when using TF-IDF. Specifically, for the naïve Bayes learning technique, mean accuracy using TF-IDF is 0.526 compared with 0.552 using Word2Vec; in the SVM technique, the mean accuracy using TF-IDF is 0.529 compared with 0.551 using Word2Vec.

\begin{figure}[!h]
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.98\linewidth]{img/fig7_a.png}
  \caption{Naïve Bayes}
  \label{fig:fig7-a}
\end{subfigure}
\quad
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.98\linewidth]{img/fig7_b.png}
  \caption{SVM}
  \label{fig:fig7-b}
\end{subfigure}

\caption{Accuracy variation by text representation model}
\label{fig:fig7}

\end{figure}

Fig.~\ref{fig:fig8} shows the accuracy by machine learning technique, namely naïve Bayes and SVM. Here, we select the same three movies from different genres as before. Fig.~\ref{fig:fig8-a} shows the accuracy of naïve Bayes and that of SVM when we keep the text representation model fixed as TF-IDF. Fig.~\ref{fig:fig8-b} shows the respective accuracy with the text representation model fixed as Word2Vec. The results show slightly better overall accuracy when using SVM than when using naïve Bayes. Specifically, in TF-IDF, mean accuracy using naïve Bayes is 0.526 compared with 0.529 using SVM. In Word2Vec, mean accuracy using naïve Bayes is 0.551 compared with 0.551 using SVM.

\begin{figure}[!h]
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.99\linewidth]{img/fig8_a.png}
  \caption{TF-IDF}
  \label{fig:fig8-a}
\end{subfigure}
\quad
\begin{subfigure}{0.47\linewidth}
  \centering
  \includegraphics[width=.99\linewidth]{img/fig8_b.png}
  \caption{Word2Vec}
  \label{fig:fig8-b}
\end{subfigure}

\caption{Accuracy variation by machine learning technique}
\label{fig:fig8}
\end{figure}

\section{Discussions}
\label{sec:sec6}

\noindent {\bf Language-independent method: } In this paper, we use movie reviews collected from the largest Korean movie review website, Naver Movie, because it satisfies two requirements to show the effectiveness of the proposed method: 1) it maintains all the historical movie reviews written by each user, as required for the proposed historical credibility method; 2) it maintains helpful and unhelpful votes for each review, as required for the comparison method based on the helpfulness vote. We can apply and test the proposed method to movie reviews collected from other websites only if they satisfy the two requirements above. We note that the proposed method is language-independent. The accuracy of the learning model might vary according to the language used in movie reviews. In this paper, we have focused on the relative accuracy difference between methods, not on the absolute value of the accuracy for the target language. This relative difference should be maintained in different languages.

\noindent {\bf Variation in accuracy according to text representation model and machine learning technique: } In this paper, we aim to show the effectiveness of the proposed learning method, rather than the superiority of a specific text representation model or a specific machine learning technique. For this reason, we choose two representative text representation models and two representative machine learning techniques. However, we could apply any other text representation model and/or machine learning technique to the proposed method.

\section{Conclusions}
\label{sec:sec7}

In this paper, we have proposed a weakly supervised learning method for judging the credibility of movie reviews. The proposed learning method has two advantages. First, it is efficient because it allows fast annotation according to predefined rules. Second, it is effective because it uses a simple but effective criterion to judge movie review credibility based on historical ratings.

For the experiments, we collected 40,000 actual movie reviews. We have shown that the proposed learning method is efficient because it can annotate 8,000 movie reviews in only 0.712 seconds, which occupies only about 0.145\% of the entire training process when we use Word2Vec as the text representation model. We have also measured the accuracy of the proposed learning method. By comparing accuracy between the proposed learning method and the comparison method based on the helpfulness vote, which has been used as an effective criterion in many studies, we have shown that the accuracy of the proposed method is better than that of the comparison method by 1.57\% $\sim$ 4.54\% depending on the text representation model, machine learning technique, and movie genres used.


\section{Introduction}

	In the past few years, because of the great achievements in observations, especially the first detection of the gravitational wave(GW) \cite{Abbott:2016blz,Abbott:2016nmj,Abbott:2017gyy} and  the first-ever picture captured for the black hole  \cite{Akiyama:2019cqa,Akiyama:2019eap}, the study of black hole physics has entered the golden era. These observations enable us to test general relativity (GR) in  the strong gravity regime \cite{Carter:1971zc,Robinson:1975bv,Chrusciel:2012jk}. In GR,  the Kerr metric describes the most general black hole spacetime (without charge). Besides mass and rotation, the Kerr black hole does not have additional hairs, which is the famous no-hair theorem. However this no-hair theorem is challenged in the modified gravity theories (MOGs), where the general hairy rotating black hole solutions were found \cite{Sotiriou:2015pka,Herdeiro:2015waa}. More interestingly it was found that the  Kerr  solution can also be allowed in MOGs, although in MOGs the perturbation on the same Kerr
	background is in general different from that in GR. This actually gives a  possible way to distinguish MOGs from GR though the study of perturbation dynamics.  Employing the current  available observational constraints, some MOGs, such as scalar-tensor theory, scalar-Einstein-Gauss-Bonnet theory (sEGB), dynamical Chern-Simon gravity (dCSG) and Lorentz-violating gravity (see, e.g., \cite{Berti:2018cxi} for a review) were found viable. It is intriguing to further disclose the signature of MOGs by uncovering the objective picture of wave dynamics. This in turn can check the no-hair theorem and unveil the hairy black holes imprints. Precise pictures in wave dynamics around hairy black holes can help understand better physics in MOGs through future GW detections.

	Among MOGs, the sEGB theory has recently attracted a lot of attentions. In this theory, an additional scalar field coupling to the Gauss-Bonnet (GB) term is added to the usual Einstein-Hilbert action. Spherical and axial symmetric black hole solutions in GR were found in the sEGB theory \cite{Antoniou:2017acq,Doneva:2017bvd,Silva:2017uqg,Cunha:2019dwb,Herdeiro:2020wei,Berti:2020kgk}. However, when the coupling between the scalar field and GB term becomes strong enough, the effective mass square can become negative in the vicinity of the black hole horizon, resulting in the tachyonic instability which leads to the so-called spontaneous scalarization and the formation of the hairy black holes.  The phenomenon of spontaneous scalarization was observed a long time ago in neutron stars, but there the instability is induced by the surrounding matter instead of the curvature \cite{Damour:1993hw}. Hairy black holes show the signature of MOGs. A lot of efforts have been devoted to study the physical implications of these hairy black holes, including black hole shadows \cite{Cunha:2019dwb}, the classical stability \cite{Blazquez-Salcedo:2018jnn,Silva:2018qhn},  etc. Results that show deviations from that of the Kerr solution indicate possible observational signatures to test sEGB theory in the future.

	Perturbations around black holes can reflect characteristic sounds of black holes which can be recorded by GW detectors. Studying matter perturbations such as the scalar field perturbations in the black hole background has both theoretical and observational interests. Theoretically, it has been known for a long time that the time evolution of the perturbations in the black hole background will generally experience three stages \cite{Vishveshwara:1970zz}. After the initial pulse, the perturbation field undergoes the damped oscillations,  called  the quasinormal ringing with frequencies and damping times  entirely fixed by the black hole parameters.  At late times, quasinormal oscillations are swamped by the relaxation process, which is the requirement of the black hole no-hair theorem. For more details, please refer to the reviews \cite{Kokkotas:1999bd,Nollert:1999ji,Berti:2009kk,Konoplya:2011qq}. The coupling between  the scalar field and GB term must influence the quasinormal modes (QNMs) and  their late time wave dynamics.  When the coupling is strong enough, the QNMs can show its effect and the late time wave behavior can reflect the destroy of the GR black hole,   because of the tachyonic instability.  The blow up of the perturbation tail illustrates the breakdown of the no-hair theorem and the birth of the hairy black hole. In observations, through GW detections, the MOGs can be confirmed if the signatures of the wave dynamics in sEGB theory are detected. It is expected that current and further GW detectors will detect the signatures of the ringdown and late-time stages \cite{Berti:2018vdi}.

	In GR, even for the most general Kerr background, thanks to the separability of the scalar perturbation equation, the QNMs have been extensively studied either in the frequency domain or time domain, for references on this topic please see \cite{Kokkotas:1999bd,Nollert:1999ji,Berti:2009kk,Konoplya:2011qq}. However, the situation becomes more complicated in sEGB, since the scalar perturbation equation cannot be separated with  usual methods. Because of this reason, in the study of the spontaneously scalarized Kerr black holes, the objective picture of perturbation wave dynamics was not presented explicitly in \cite{Cunha:2019dwb} when the coupling between the scalar field and GB curvature correction term $\alpha>0$.
	This technical obstacle is common in MOGs. Fortunately, several methods, mostly numerical, have been developed in the past few years to deal with such problems  \cite{Krivan:1996da,PazosAvalos:2004rp,Dolan:2011dx,Thuestad:2017ngu,Dolan:2012yt,Brito:2014nja,Zenginoglu:2007jw,Zenginoglu:2008pw,Zenginoglu:2008uc,Zenginoglu:2008wc,Zenginoglu:2009ey,Zenginoglu:2009hd,Zenginoglu:2010cq,Zenginoglu:2011jz,Racz:2011qu,Zenginoglu:2011zz,Cano:2020cao}. Among them, a general method is to study the full time evolution of the scalar perturbation by a $(1+1)$ or $(2+1)$-dimensional simulation. In  \cite{Dima:2020yac}, by projecting the scalar field onto a basis set of spherical harmonics to separate the angular dependencies, the authors considered the case with $\alpha<0$ and studied the time evolution of the scalar perturbation using  $(1+1)$-dimensional simulations. The results show that instability only occurs at sufficiently high spins $(a/M > \frac{1}{2})$ (see also Ref. \cite{Hod:2020jjy} for an analytical argument). For a  given big enough spin $a$, there exists a critical value $|\alpha_c|$, above which instability occurs. However,  $|\alpha_c|$ depends on the parameter $a$, and in particular
	it decreases as $a$ is increasing.  Further, these results were confirmed by using the $(2+1)$-dimensional simulation with a numerical method introduced in \cite{Doneva:2020nbb}. It was shown that the instability with $\alpha<0$ leads to the dubbed spin-induced scalarization and the final hairy black holes can be constructed \cite{Herdeiro:2020wei,Berti:2020kgk}. See also a related work Ref. \cite{Konoplya:2019hml} concerning the time evolution of a test and minimally coupled scalar field in this theory.

	In available numerical treatments, the tortoise coordinate was utilized to map the radial computational domain $(r_+, +\infty)$ to $(-\infty, +\infty)$ (here $r_+$ is the horizon). However, Price {\it et al}. pointed out that traditional numerical methods suffer the so-called ``outer boundary problem" \cite{Thuestad:2017ngu}: namely in practical calculations, one has to truncate the infinite radial computational domain to a finite range and put boundary conditions (ingoing and outgoing waves at the horizon and infinity,
	 respectively) at the outer edges, thus inevitably resulting spurious wave reflections from the edge which can spoil the evolution of the scalar perturbation at the late time. To overcome this problem, a technique called hyperboloidal compactification was developed in \cite{Zenginoglu:2007jw,Zenginoglu:2008pw,Zenginoglu:2008uc,Zenginoglu:2008wc,Zenginoglu:2009ey,Zenginoglu:2009hd,Zenginoglu:2010cq,Zenginoglu:2011jz,Zenginoglu:2011zz}, with which the radial computation domain is mapped to a finite range and the outer boundary conditions are satisfied automatically. Interestingly, armed with this numerical strategy, only until recently we obtain a well understanding on the late-time behaviors of  the scalar field  perturbation even for
	 the Kerr black hole \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}. This technique has also been successfully applied to study gravitational waveforms from large- and extreme-mass-ratio inspirals in
	  the  Kerr black hole extracted at null infinity for the first time \cite{Zenginoglu:2011zz}, and also the full-time evolution of scalar perturbations in dCSG \cite{Gao:2018acg}.

	It is interesting to further apply the hyperboloidal compactification technique to study the scalar perturbation in the sEGB theory when there is coupling between the scalar field and GB curvature correction.  We expect to get a complete and objective picture on the wave dynamics to describe precisely the time evolution of the scalar perturbation in the sEGB theory. Such a picture can not only provide us information about the stability of the black hole, serving as an independent confirmation of the scalarization discussions \cite{Cunha:2019dwb, Dima:2020yac, Doneva:2020nbb,Hod:2020jjy}, but also it can present further imprints on the scalar dynamics, including QNMs and the late-time behavior influenced by the coupling constant between the scalar field and GB term. Precise waveforms can reflect the effect of the coupling $\alpha$ in QNMs and suggest the potential verification of the sEGB theory.  We will examine the objective picture of the scalar evolution with $\alpha>0$ as well as $\alpha<0$ in the time domain,  and obtain a complete understanding on the effects of the coupling constant on the scalar wave dynamics.

	This paper is organized as follows. In Sec. II, we give a brief introduction of the sEGB theory and write out the scalar perturbation equation. In Sec. III, applying the hyperboloidal compactification technique, we cast the scalar perturbation equation into a form suitable for numerical calculations. In Sec. IV, we first check the validity of the hyperboloidal compactification technique for the case $\alpha=0$, and then we report numerical results for $\alpha>0$ and $\alpha<0$, respectively, in sEGB theory. The last section is devoted to summary and discussions.

	\section{Scalar-Einstein-Gauss-Bonnet theory and scalar perturbation equation}

	The action of sEGB theory is \cite{Antoniou:2017acq,Doneva:2017bvd,Silva:2017uqg,Cunha:2019dwb}
	\begin{eqnarray}
		S&=& \frac{1}{2\kappa}\int dx^4\sqrt{-g}\left(R-2\Lambda +\alpha{\cal L}_{GB}+{\cal L}_{\Phi}\right),\nonumber\\
		{\cal L}_{GB}&=& f(\Phi){\cal G},\nonumber\\
		{\cal L}_\Phi&=&-\frac{1}{2}\nabla^\mu\Phi\nabla_\mu\Phi-V(\Phi),\nonumber\\
		{\cal G}&=&R^2+R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}-4R_{\mu\nu}R^{\mu\nu},
	\end{eqnarray}
	where the scalar field $\Phi$ is non-minimally coupled to the Gauss-Bonnet term ${\cal G}$ with the coupling constant $\alpha$. $f(\Phi)$ is a function of the scalar field and $\Lambda$ is the cosmological constant. From the action, one can derive the equations of motion
	\begin{eqnarray}
		\nabla^2\Phi&=&\frac{dV}{d\Phi}-\alpha \frac{d f}{d\Phi}{\cal G},\label{ScalarEq}\\
		R_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R + \Lambda g_{\mu\nu}&=&\alpha T^{GB}_{\mu\nu}+T^\Phi_{\mu\nu},\label{MetricEq}\\
		T^{GB}_{\mu\nu}&=&2(\nabla_\mu\nabla_\nu f)R-2g_{\mu\nu}(\nabla_\rho\nabla^\rho f)R \nonumber\\
		&&-4(\nabla^\rho\nabla_\nu f)R_{\mu\rho}-4(\nabla^\rho\nabla_\mu f)R_{\nu\rho}\nonumber\\
		&&+4(\nabla^\rho\nabla_\rho f)R_{\mu\nu}+4g_{\mu\nu}(\nabla^\rho\nabla^\sigma f)R_{\rho\sigma}\nonumber\\
		&&-4(\nabla^\rho\nabla^\sigma f)R_{\mu\rho\nu\sigma},\nonumber\\
		T^\Phi_{\mu\nu}&=&\frac{1}{2}\nabla_\mu\Phi\nabla_\nu\Phi-\frac{1}{2} g_{\mu\nu}V(\Phi)-\frac{1}{4} g_{\mu\nu}\nabla^\rho\Phi\nabla_\rho \Phi.\nonumber
	\end{eqnarray}
	The theory admits GR black hole solutions with constant scalar profile $\Phi=\Phi_0$,  if
	\begin{eqnarray}
		V(\Phi_0)=0,\quad \frac{dV}{d\Phi}\bigg|_{\Phi_0}=0, \quad \frac{df}{d\Phi}\bigg|_{\Phi_0}=0.
	\end{eqnarray}
	In the following, we will consider a simple case as  \cite{Cunha:2019dwb} by choosing $\Lambda=0, V(\Phi)=0$ and
	\begin{eqnarray}
		f(\Phi)=\frac{1}{2\beta} \left(1-e^{-\beta \Phi^2}\right),
	\end{eqnarray}
	where $\beta>0$ is a constant.

	We are going to study the wave dynamics of scalar field perturbations on the background of GR black holes in the linear regime. More precisely, we consider  the Kerr black hole with $\Phi=0$. The metric in the
	Boyer-Lindquist coordinates is
	\begin{eqnarray}
		ds^2 = -\frac{\Delta}{\Sigma} (dt - a \sin^2 \theta d\phi)^2 + \frac{\Sigma}{\Delta} dr^2 + \Sigma d\theta^2 + \frac{\sin^2\theta}{\Sigma} (a dt - (r^2+a^2)d\varphi)^2,
	\end{eqnarray}
	where $\Delta\equiv r^2 - 2Mr+a^2$ and $\Sigma \equiv r^2 + a^2 \cos^2\theta$. In this case, the scalar perturbation equation (\ref{ScalarEq}) in the Kerr background reduces to
	\begin{eqnarray}
		&&\nabla^2\Phi=-\alpha {\cal G} \Phi, \label{ScalarEq2}\\
		&&{\cal G} = -\frac{48 M^2 \left(-32 r^6+48 r^4 \Sigma -18 r^2 \Sigma^2+\Sigma^3\right)}{\Sigma^{6}},\nonumber
	\end{eqnarray}
	where ${\cal G}$ is valued in the background. As one can see,  the scalar field  has an effective mass $m^2_{\rm eff}= -\alpha {\cal G}$, which is position dependent and approaches zero at infinity, with the sign depending on the coupling constant $\alpha$.
	When $\alpha=0$, the above equation describes wave propagations of a free scalar field in the Kerr background which has been studied thoroughly  \cite{Krivan:1996da,PazosAvalos:2004rp,Thuestad:2017ngu}. For  $\alpha>0$,  in the Schwarzschild background, one has ${\cal G} = 48 M^2/r^6$ so that the effective mass square is always negative, which leads to the tachyonic instability and the subsequent spontaneous scalarization. In \cite{Doneva:2017bvd,Silva:2017uqg} it was found that the instability occurs when $M/(2 \sqrt{\alpha}) \lesssim 0.587$ (note that in \cite{Doneva:2017bvd,Silva:2017uqg}  $\lambda^2$ was used to represent the coupling constant which is related to $\alpha$ by $\lambda^2 = 4 \alpha$), and becomes more violent for larger $\alpha$ or smaller $M$. However, for a Kerr background, the situation becomes more complicated, because  ${\cal G}$ is not a  monotonic function so that the effective mass square is not always negative. How does this effective mass affects the stability of the Kerr black hole configuration and triggers spontaneous scalarization  requires careful investigations.

	In the following sections, we will study carefully the time evolution of the scalar perturbation and obtain  object pictures on the influences of the coupling constant $\alpha$ on wave dynamics in the sEGB theory.

	\section{Setup of the numerical method}

	As mentioned in the introduction, traditional numerical methods usually encounter the so-called ``outer boundary problem", which can spoil the numerical precision, especially for the evolution of the scalar perturbation at the late time \cite{Thuestad:2017ngu}. To overcome this problem, we  apply the numerical strategy proposed in \cite{Zenginoglu:2007jw,Zenginoglu:2008pw,Zenginoglu:2008uc,Zenginoglu:2008wc,Zenginoglu:2009ey,Zenginoglu:2009hd,Zenginoglu:2010cq,Zenginoglu:2011jz,Zenginoglu:2011zz} to solve the scalar perturbation equation (\ref{ScalarEq2}), which contains two steps of coordinates transformations mainly. First, we define the ingoing coordinates $\{\tilde{t}, r, \theta, \tilde{\varphi}\}$ through the following transformation
	\begin{eqnarray}
		d\tilde{t} = dt + \frac{2 M r}{\Delta} dr,\qquad d\tilde{\varphi} = d\varphi + \frac{a}{\Delta} dr,
	\end{eqnarray}
	and the metric becomes
	\begin{eqnarray}
		ds^2 =&& -\left(1-\frac{2 M r}{\Sigma}\right) d\tilde{t}^2 - \frac{4 a M r}{\Sigma} \sin^2\theta d\tilde{t} d\tilde{\varphi} + \frac{4 M r}{\Sigma} d\tilde{t} dr\nonumber\\
		&& + \left(1+\frac{2 M r}{\Sigma}\right) dr^2 - 2 a \sin^2\theta \left(1+\frac{2 M r}{\Sigma}\right) dr d\tilde{\varphi} + \Sigma d\theta^2\nonumber\\
		&& + \left(r^2 + a^2 + \frac{2 M a^2 r \sin^2\theta}{\Sigma}\right) \sin^2\theta d\tilde{\varphi}. \label{KerrSchildMetric}
	\end{eqnarray}
	Considering the axial symmetry of the Kerr spacetime, the scalar perturbation can be decomposed as
	\begin{eqnarray}
		\Phi (\tilde{t},r,\theta,\tilde{\varphi}) = \frac{1}{r} \sum_m \Psi(\tilde{t},r,\theta) e^{i m \tilde{\varphi}}.
	\end{eqnarray}
	Substituting the above ansatz into  Eq.~(\ref{ScalarEq2}), the scalar perturbation equation in the ingoing  coordinates becomes
	\begin{eqnarray}
		A^{\tilde{t}\tilde{t}} \partial_{\tilde{t}}^2 \Psi + A^{\tilde{t} r} \partial_{\tilde{t}} \partial_{r} \Psi + A^{r r} \partial_{r}^2 \Psi + A^{\theta\theta} \partial_{\theta}^2 \Psi + B^{\tilde{t}} \partial_{\tilde{t}} \Psi + B^r \partial_r \Psi + B^\theta \partial_{\theta}\theta \Psi + C \Psi =0,
	\end{eqnarray}
	where
	\begin{eqnarray}
		A^{\tilde{t}\tilde{t}} &=& \Sigma + 2 M r,\nonumber\\
		A^{\tilde{t} r} &=& -4 M r,\nonumber\\
		A^{r r} &=& -\Delta,\nonumber\\
		A^{\theta\theta} &=& -1,\nonumber\\
		B^{\tilde{t}} &=& 2 M,\nonumber\\
		B^r &=& \frac{2}{r} (a^2 - M r) - 2 i m a,\nonumber\\
		B^\theta &=& -\cot\theta,\nonumber\\
		C &=& \frac{m^2}{\sin^2\theta}- \frac{2(a^2-M r)}{r^2} + \frac{2 i m a}{r} - \alpha \Sigma \cal{G}.
	\end{eqnarray}

	The second step is to use the technique of hyperboloidal compactifications developed in \cite{Harms:2014dqa},  which defines a compactified horizon-penetrating, hyperboloidal coordinates $(\tau,\rho,\theta,\tilde{\varphi})$ (HH coordinates) through the transformation \footnote{Of course, in this step, one can have other choices, such as the RT coordinates proposed in  \cite{Racz:2011qu}.}
	\begin{eqnarray}
		\tilde{t} = \tau + h(\rho),\qquad r = \frac{\rho}{\Omega(\rho)},\label{HHCoordinates}
	\end{eqnarray}
	where
	\begin{eqnarray}
		h(\rho) = \frac{\rho}{\Omega} - \rho -4 M \ln\Omega,\qquad \Omega (\rho) = 1- \frac{\rho}{S}.
	\end{eqnarray}
	Here $S$ is a free parameter controlling both the domain and the foliation. With this compactification, the semi-infinite radial domain outside the horizon $r\in [r_+, \infty)$ is mapped to a finite range $\rho \in [\rho_+, S)$ with the event horizon $r=r_+$ now locates at
	\begin{eqnarray}
		\rho_+ = \frac{a^2 S + M S^2 + \sqrt{M^2 S^4 - a^2 S^4}}{a^2 + 2 M S +S^2}.
	\end{eqnarray}
	Then we have the relations
	\begin{eqnarray}
		\partial_{\tilde{t}} = \partial_\tau,\qquad \partial_r = - H \partial_\tau + K \partial_\rho,
	\end{eqnarray}
	where $H \equiv \frac{d h}{dr} (\rho)$ and $K \equiv \frac{d\rho}{dr} (\rho)$. With these relations, the scalar perturbation equation in the HH coordinates can be written as
	\begin{eqnarray}
		\partial_\tau^2 \Psi = \tilde{A}^{\tau\rho} \partial_\tau \partial_\rho \Psi + \tilde{A}^{\rho\rho} \partial_\rho^2 \Psi + \tilde{A}^{\theta\theta} \partial_{\theta}^2 \Psi + \tilde{B}^\tau \partial_\tau \Psi + \tilde{B}^\rho \partial_\rho \Psi + \tilde{B}^\theta \partial_{\theta} \Psi + \tilde{C} \Psi,
	\end{eqnarray}
	where
	\begin{eqnarray}
		\{\tilde{A}^{\tau\rho}, \tilde{A}^{\rho\rho}, \tilde{A}^{\theta\theta}, \tilde{A}^\tau, \tilde{B}^\tau, \tilde{B}^\rho, \tilde{B}^\theta, \tilde{C}\} = -\frac{1}{A^{\tau\tau}} \{A^{\tau\rho}, A^{\rho\rho}, A^{\theta\theta}, A^\tau, B^\tau, B^\rho, B^\theta, C\},
	\end{eqnarray}
	and
	\begin{eqnarray}
		A^{\tau\tau} &=& A^{\tilde{t}\tilde{t}} - H A^{\tilde{t} r} + H^2 A^{rr},\nonumber\\
		A^{\tau\rho} &=& K A^{\tilde{t} r} - 2 K H A^{rr},\nonumber\\
		A^{\rho\rho} &=& K^2 A^{rr},\nonumber\\
		B^\tau &=& B^{\tilde{t}} - H B^r - \frac{d H}{d\rho} K A^{rr},\nonumber\\
		B^\rho &=& K \left(B^r + \frac{d K}{d\rho} A^{rr}\right),
	\end{eqnarray}
	To solve the above equation numerically in the time domain, it is convenient to rewrite it as two coupled first-order partial differential equations by introducing a new variable $\Pi \equiv \partial_\tau \Psi$,
	\begin{eqnarray}
		\partial_\tau \Psi &=& \Pi,\nonumber\\
		\partial_\tau \Pi &=& \tilde{B}^\tau \Pi + \tilde{A}^{\tau\rho} \partial_\rho \Pi + \tilde{A}^{\rho\rho} \partial_\rho^2 \Psi + \tilde{B}^\rho \partial_\rho \Psi+ \tilde{A}^{\theta\theta} \partial_{\theta}^2 \Psi   + \tilde{B}^\theta \partial_{\theta} \Psi + \tilde{C} \Psi.
	\end{eqnarray}

	Numerically, this form is suitable for the method mentioned above \cite{Schiesser}. Precisely, the derivatives in $\rho$ and $\theta$ directions are approximated by finite differences, and the evolution in
	the time direction is implemented with the fourth-order Runge-Kutta integrator. One of the main advantages of working in the HH coordinates is that the ingoing (outgoing) boundary condition
	at the horizon (infinity) is satisfied automatically, so we do not need to worry about the sophisticated ``outer boundary problem" which is hard to handle and influences the precision. At the poles in
	the angular direction $\theta=0$ and $\pi$, physical boundary conditions are needed
	\begin{eqnarray}
		\Psi|_{\theta=0, \pi} &=& 0 \quad {\rm for \ odd} \ m,\\
		\partial_\theta \Psi|_{\theta=0, \pi} &=& 0 \quad {\rm for \ even}\  m.
	\end{eqnarray}
	Following \cite{PazosAvalos:2004rp}, we use a staggered grid and add ghost points to implement these conditions.

	\section{Results}

	We consider the initial data of the scalar perturbation to be a Gaussian distribution localized outside the horizon at $\rho=\rho_c$,
	\begin{eqnarray}
		\Psi (\tau=0,\rho,\theta) \sim Y_{\ell m} e^{-\frac{(\rho-\rho_c)^2}{2 \sigma^2}},
	\end{eqnarray}
	where $Y_{\ell m}$ is the $\theta$-dependent part of the spherical harmonic function and $\sigma$ is the width of the Guassian distribution. Without loss of generality, we consider the perturbation to have time symmetry so that
	\begin{eqnarray}
		\Pi (\tau=0,\rho,\theta)=0.
	\end{eqnarray}
	In the following, we take $\rho_c=6 M$. Also, we set $M=1$ so that all quantities are measured in units of $M$. Observers are assumed to locate at $\rho =6 M$ and $\theta=\frac{\pi}{4}$.  We choose the free parameter $S=10$ as in \cite{Harms:2014dqa}.

	The Kerr spacetime is not spherically symmetric, except when $a=0$, so the mode-mixing phenomenon occurs \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}: a pure even (odd) initial $\ell$ multipole will excite other  even (odd) multipoles (denoted by $\ell'$) with the same $m$ as it  evolves. So, in this work we only consider $\ell=0$ and $\ell=1$ as representative even and odd multipoles, respectively. Also, we only consider axisymmetric perturbations with $m=0$ for simplicity.

	\subsection{$\alpha=0$}

	To check the reliability of our numerical strategy so as to further disclose  precisely the influence of $\alpha$ on the wave dynamics, we first consider the simple case where the scalar field does not couple to the GB term, i.e., $\alpha=0$. We compare our results with those in the ringdown phase and the late-time tail with available results either in the time domain or in the frequency domain. In this case with zero $\alpha$, the scalar field equation is separable for the $r$ and $\theta$ variables and can be cast into a Schr$\ddot{o}$dinger-like form. Then, QNMs can be computed directly in the frequency domain by using Leaver's continued fraction method \cite{Kokkotas:1999bd,Nollert:1999ji,Berti:2009kk,Konoplya:2011qq}.

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0$ ]{\includegraphics[width=0.47\textwidth]{Phi-alpha0a0m0.pdf}} \quad
		\subfigure[$~~a=0.5$ ]{\includegraphics[width=0.47\textwidth]{Phi-alpha0a05m0.pdf}}
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{Phi-alpha0a09m0.pdf}}
		\caption{(color online) Time evolution of the axisymmetric scalar perturbation for $a=0, 0.5$ and $0.9$. The initial multipoles we considered are $\ell=0, 1$ and $2$. Time is in units of $M$.}
		\label{alpha0Figs}
	\end{figure}

	In Fig. \ref{alpha0Figs}, the time evolutions of the scalar perturbation in the Kerr background are plotted for various spins. We only consider axisymmetric initial multipoles ($m=0$) with various $\ell$'s. From the figure, we can read off the quasinormal frequencies by using the Prony's method \cite{Berti:2007dg}. The results are listed in Table \ref{QNMsIndex}  denoted as $\omega_{\rm Prony}$. To make a comparison, the results derived by Leaver's method in the frequency domain \cite{Berti:2009kk,Dolan:2007mj} are also listed denoted as $\omega_{\rm Leaver}$. From the table, we can see that using the numerical method described in the last section, our obtained results  agree well with previous ones. The relatively large discrepancy for the $\ell=0$ case is due to the short duration of the ringdown phase so that Prony's method can not give quasinormal frequencies with sufficient precisions.

	\begin{table}[!htbp]
		\centering
		\begin{tabular}{p{1cm}<{\centering} p{1cm}<{\centering} p{3cm}<{\centering} p{3cm}<{\centering} p{1cm}<{\centering} p{2cm}<{\centering}}
			\hline
			\hline
			$a$ & $\ell$ & $\omega_{\rm Prony}$ & $\omega_{\rm Leaver}$ & $n_{\rm fitting}$ & $n_{\rm predicted}$\\
			\hline
			\multirow{3}*{$0$} & $0$ & $0.1084-0.1062i$ & $0.1105-0.1049i$ & $3.16$ & $3$\\
			&$1$ & $0.2929-0.0976i$ & $0.2929-0.0977i$ & $5.10$ & $5$\\
			&$2$ & $0.4836-0.0968i$ & $0.4836-0.0968i$ & $7.32$ & $7$\\
			\hline
			\multirow{3}*{$0.5$} & $0$ & $0.1124-0.0935i$ & $0.1123-0.1022i$ & $3.16$ & $3$\\
			& $1$ & $0.2979-0.0954i$ & $0.2979-0.0954i$ & $5.18$ & $5$\\
			& $2$ & $0.4920-0.0946i$ & $0.4920-0.0946i$ & $3.38$ & $3$\\
			\hline
			\multirow{3}*{$0.9$} & $0$ & $0.1186-0.0901i$ & $0.1138-0.0916i$ & $3.18$ & $3$\\
			& $1$ & $0.3108-0.0866i$ & $0.3108-0.0866i$ & $5.18$ & $5$\\
			& $2$ & $0.5148-0.0864i$ & $0.5148-0.0864i$ & $3.35$ & $3$\\
			\hline\hline
		\end{tabular}
		\caption{Dominant QNMs and the power-law index of the late-time tail for axisymmetric scalar perturbations with $a=0, 0.5$ and $0.9$. $n_{\rm predicted}$ refers to the power-law index according to Eq. (\ref{PowerLaw}) and $n_{\rm fitting}$ is our fitting result.}
		\label{QNMsIndex}
	\end{table}

	Now let us examine the late-time tail. Interestingly, it is only recently and with the hyperboloidal compactification technique that we have a precise and well understanding of the late-time behaviors of the scalar field perturbation in the Kerr black hole background \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}. For axisymmetric perturbations ($m=0$), the late-time behaviors of these multipoles exhibit a power-law decay $\Psi \sim \tau^{-n}$ at time infinities, where
	\begin{eqnarray}
		n = \left\{
		\begin{array}{ll}
			\ell+\ell'+3 \quad {\rm for} \ \ell=0,1\\
			\ell+\ell'+1 \quad {\rm otherwise}
		\end{array}
		\right. \label{PowerLaw}
	\end{eqnarray}
	In the full late-time, the dominant mode is $\ell'=0$ for the even $\ell$ case and $\ell'=1$ for the odd $\ell$  case (for more details we refer readers to Refs. \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}).

	From Fig. \ref{alpha0Figs}, we can see that in the late-time, the perturbation exhibits a power-law decay as expected. By fitting the late-time data with function $\tau^{-n}$, we can find the index that is listed in Table \ref{QNMsIndex} denoted as $n_{\rm fitting}$. It should be noted that the power-law decay (\ref{PowerLaw}) only holds exactly at very late-time, so to make a precise comparison we need to evolve the perturbation to a very late time which will be extremely time consuming. However, in the time range we show that the index $n_{\rm fitting}$ and its trend with time derived by fitting agrees well with the prediction (\ref{PowerLaw}).

	\subsection{$\alpha>0$}

	Now let us consider the case with $\alpha>0$. In Fig. \ref{alphaPositiveL0m0Fig}, the time evolution behaviors of the scalar perturbation for various positive $\alpha$'s are plotted. The initial multipole is taken as $\ell=0$. We first look at the wave propagation after the ringing. From the figure, one can see that when $a=0$, there exists a critical value  $\alpha_c \sim 0.726$, above which the perturbation grows,  indicating the occurrence of the instability. The instability happens earlier and becomes more violent as $\alpha$ further increases. Similar behaviors of the stability also appear for $a \neq 0$.  We learn that the instability starts at $\alpha_c=  0.75, 0.726, 0.71$ and $0.6$  for different rotations $a=0.3, 0.5, 0.7$ and $0.9$,  respectively. The dependence of $\alpha_c$ on $a$  is not simply monotonic as the increase of the spin from zero to the extreme value, $\alpha_c$ first increases when the hole starts to rotate and then changes to decrease as the hole rotates faster until the extreme value. Using a different and more precise numerical strategy described in the last section, we have shown  that our results agree quantitatively well with those computed in  \cite{Doneva:2017bvd,Silva:2017uqg}.

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0$ ]{\includegraphics[width=0.47\textwidth]{Phi-a0L0m0alphaPositive.pdf}} \quad
		\subfigure[$~~a=0.3$ ]{\includegraphics[width=0.47\textwidth]{Phi-a03L0m0alphaPositive.pdf}}
		\subfigure[$~~a=0.5$ ]{\includegraphics[width=0.47\textwidth]{Phi-a05L0m0alphaPositive.pdf}}\quad
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{Phi-a09L0m0alphaPositive.pdf}}
		\caption{(color online) Time evolution of the scalar perturbation for $a=0, 0.3, 0.5$ and $0.9$ with $\alpha>0$.  Initial multipole is fixed as $\ell=0$. Time is in units of $M$.}
		\label{alphaPositiveL0m0Fig}
	\end{figure}

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0$ ]{\includegraphics[width=0.47\textwidth]{Phi-a0L1m0alphaPositive.pdf}} \quad
		\subfigure[$~~a=0.5$ ]{\includegraphics[width=0.47\textwidth]{Phi-a05L1m0alphaPositive.pdf}}
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{Phi-a09L1m0alphaPositive.pdf}}
		\caption{(color online) Time evolution of the scalar perturbation for $a=0, 0.5$ and $0.9$ with $\alpha>0$. Initial multipole is fixed as $\ell=1$. Time is in units of $M$.}
		\label{alphaPositiveL1m0Fig}
	\end{figure}

	\begin{table}[!htbp]
		\centering
		\begin{tabular}{p{1cm}<{\centering} p{3cm}<{\centering} p{3cm}<{\centering} p{3cm}<{\centering}}
			\hline
			\hline
			\multirow{2}*{$\alpha$} & \multicolumn{3}{c}{$a$}\\
			\cline{2-4}
			&  $0$ & $0.5$ & $0.9$ \\
			\hline
			$-2.0$ & $0.3870-0.1105i$  & $0.3613-0.1067i$  & $-$\\
			$-1.5$ & $0.3633-0.1062i$  & $0.3456-0.1034i$  & $0.3387-0.1150i$\\
			$-1.0$ & $0.3397-0.1020i$ & $0.3300-0.1001i$  & $0.3118-0.1122i$\\
			$-0.5$ & $0.3163-0.0988i$ & $0.3139-0.0975i$ & $0.3088-0.1019i$ \\
			$0$ & $0.2929-0.0976i$ &  $0.2979-0.0954i$ & $0.3108-0.0866i$\\
			$0.5$ & $0.2703-0.1002i$ & $0.2817-0.0945i$ &  $0.3210-0.0747i$\\
			$0.8$ & $0.2580-0.1046i$ &  $0.2721-0.0948i$ & $0.3281-0.0711i$\\
			$1.0$ & $0.2512-0.1088i$ &  $0.2657-0.0955i$ & $0.3325-0.0700i$\\

			$2.0$ & $0.2394-0.1244i$ &  $0.2389-0.1067i$ & $0.3467-0.0729i$\\

			\hline\hline
		\end{tabular}
		\caption{Dominant quasinormal modes for the scalar perturbation with initial $\ell=1$ multipole. "-" means the duration of the ringdown phase is too short to give precise enough QNMs with Prony's method. }
		\label{QNMs}
	\end{table}

	We notice that instability also occurs for other multipoles. In Fig. \ref{alphaPositiveL1m0Fig}, the time evolution of the multipole $\ell=1$  for various values of $\alpha>0$ are shown. We observe similar wave dynamics
	behaviors due to the influence of $\alpha$ as the fundamental multipole mode. Besides, we learn that  $\alpha_c$ for $\ell=1$  is larger than its corresponding value for $\ell=0$.

	To gain a better understanding of the physics behind the above phenomena, we examine the influence of the coupling constant and the spin on the effective mass square $m^2_{\rm eff}= -\alpha {\cal G}$. When $a=0$, the background reduces to the Schwarzschild spacetime in which ${\cal G} = 48 M^2/r^6$, so that $m^2_{\rm eff}$ is always negative everywhere for any $\alpha>0$ and   becomes more negative when  $\alpha$ is further increased. Small negative $m^2_{\rm eff}<0$ is not sufficient to trigger the tachyonic instability. Only when $m^2_{\rm eff}$ is sufficiently negative ($\alpha>\alpha_c$), can the instability be developed and becomes more violent with the further increase of $\alpha$, for which
	  $m^2_{\rm eff}$ will become more negative. In Fig. \ref{mSquarePositiveAlpha}, $m^2_{\rm eff}$   is plotted for various $a \neq 0$ and $\alpha$. From the figure, one can clearly see that $m^2_{\rm eff}$ becomes negative somewhere even
	  when $\alpha$ is still smaller than the critical value. But this small negative effective mass square is not strong enough to change the bulk stability. As $\alpha$ continuously increases, at some fixed position close to the horizon, we see that $m^2_{\rm eff}$ becomes more and more negative.  The  instability for big enough $\alpha$ is triggered  by the negative enough $m^2_{\rm eff}$. Moreover, from the figure, one can see that with the increase of the spin, $m^2_{\rm eff}$ will be lifted up near the horizon in some $\theta$ direction. This cannot save the stability, since at other $\theta$ directions the effective mass square can become more negative, so that even smaller $\alpha_c$ can trigger the instability for fast rotating holes.

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0.3$ ]{\includegraphics[width=0.47\textwidth]{a03mSquareAlphaPositive.pdf}}\quad
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{a09mSquareAlphaPositive.pdf}}
		\caption{(color online) The effective mass square $m^2_{\rm eff}$  is plotted for $a=0.3$ and $0.9$. $\alpha$ is chosen as in Fig. \ref{alphaPositiveL0m0Fig}.}
		\label{mSquarePositiveAlpha}
	\end{figure}

	Now we examine the quasinormal ringing behavior. The influence of $\alpha$ on the ringdown phase is also explored by using the Prony's method, as shown in Table. \ref{QNMs}. Since the duration of the ringdown phase for the multipole  $\ell=0$ is too short to give precise enough QNMs by using Prony's method, we only present QNMs for  the $\ell=1$ multipole. From the table and the figure, one can see that when $a=0$, as $\alpha$ achieves  bigger positive values, the real part of QNMs $\Re \omega$ decreases,  while the imaginary part $|\Im \omega|$ increases. However, when $a \neq 0$, the situation becomes more complicated: as $\alpha$ increases from $0$,  $\Re \omega$ keeps decreasing for low spins but changes to increase for high spins, while  $|\Im \omega|$ first decreases and then increases. The non-monotonic change in $|\Im \omega|$ with the increase of the coupling constant becomes more obvious for fast rotation cases in QNMs. Since the QNMs  can be used  uniquely   to  identify  black
	holes (in GR) and shall be detected through GW observations (specially for the third generation detectors), it is clear that the fine structure in the QNMs due to the influence of the coupling between the scalar field and GB term  can be used to explore the signature of the sEGB theory.  Due to the  ``mode-mixing mechanism" \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}, one can expect that similar behaviors described above will also occur for higher multipoles.

	\subsection{$\alpha<0$}

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0$ ]{\includegraphics[width=0.47\textwidth]{Phi-a0L0m0alphaNegative.pdf}}\quad
		\subfigure[$~~a=0.5$ ]{\includegraphics[width=0.47\textwidth]{Phi-a05L0m0alphaNegative.pdf}}
		\caption{(color online) Time evolution of the scalar perturbation for $a=0$ and $0.5$ with $\alpha<0$. Initial multipole is fixed as $\ell=0$. Time is in units of $M$.}
		\label{alphaNegativeL0m0Fig}
	\end{figure}

	We now  consider the case with $\alpha<0$. In Fig. \ref{alphaNegativeL0m0Fig}, by fixing $a=0$ and $0.5$, the time evolution of the scalar perturbation for  multipole $\ell=0$ is plotted, where no instability is observed. This objective picture actually confirms  the argument that  there exists a minimum spin $a_{\rm min}=\frac{1}{2}$ and below which no instability can be triggered, which is irrespective of the value of $\alpha$  \cite{Dima:2020yac,Hod:2020jjy}.  To see more clearly  the influences of the spin and coupling constant on the onset of the instability, in Fig. \ref{alphaNegativeL0m0Fig1}, we fix the spin to some representative values above $a_{\rm min}$ and show the time evolution of the scalar perturbation for various $\alpha$.  From the figure, similar to the case $\alpha>0$, one can see that the instability only occurs when $\alpha$ exceeds a critical value $|\alpha_c|$ for a fixed spin. For  $a=0.7, 0,8$ and $0.9$, $\alpha_c$ is $-5.8,-2.0$ and $-0.92$,  respectively. This tells us that with the increase of the spin, $|\alpha_c|$ decreases quickly. For the extreme hole with the highest spin,  it is very simple to trigger instability because we only require the smallest $|\alpha_c|$. Moreover, we observe that above $|\alpha_c|$ for each chosen $a$, larger $|\alpha|$ makes the instability to appear earlier and behave more violently. Similar properties are also found for the multipole $\ell=1$.  These results agree well even quantitatively  with those found in \cite{Dima:2020yac,Hod:2020jjy,Doneva:2020nbb}.

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0.7$ ]{\includegraphics[width=0.47\textwidth]{Phi-a07L0m0alphaNegative.pdf}}\quad
		\subfigure[$~~a=0.8$ ]{\includegraphics[width=0.47\textwidth]{Phi-a08L0m0alphaNegative.pdf}}\quad
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{Phi-a09L0m0alphaNegative.pdf}}
		\caption{(color online) Time evolution of the scalar perturbation for various spins. $a$ is fixed as $0.7, 0.8$ and $0.9$ respectively. Initial multipole is fixed as $\ell=0$. Time is in units of $M$.}
		\label{alphaNegativeL0m0Fig1}
	\end{figure}

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0.7$ ]{\includegraphics[width=0.47\textwidth]{Phi-a07L1m0alphaNegative.pdf}}\quad
		\subfigure[$~~a=0.8$ ]{\includegraphics[width=0.47\textwidth]{Phi-a08L1m0alphaNegative.pdf}}\quad
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{Phi-a09L1m0alphaNegative.pdf}}
		\caption{(color online)  Time evolution of the scalar perturbation for various spins. $a$ is fixed as $0.7, 0.8$ and $0.9$ respectively. Initial multipole is fixed as $\ell=1$. Time is in units of $M$.}
		\label{alphaNegativeL1m0Fig}
	\end{figure}

	\begin{figure}[!htbp]
		\centering
		\subfigure[$~~a=0.5$ ]{\includegraphics[width=0.47\textwidth]{a05mSquareAlphaNegative.pdf}}\quad
		\subfigure[$~~a=0.9$ ]{\includegraphics[width=0.47\textwidth]{a09mSquareAlphaNegative.pdf}}
		\caption{(color online) The effective mass square $m^2_{\rm eff}$   is plotted for $a=0.5$ and $0.9$. $\alpha$ is chosen as in Figs. \ref{alphaNegativeL0m0Fig} and \ref{alphaNegativeL0m0Fig1}. In the right panel, the surface with $\alpha=-0.92$ overlaps largely with that with $\alpha=-1.0$.}
		\label{mSquareNegativeAlpha}
	\end{figure}

	Similar to the discussion for $\alpha>0$, we can also understand the instability from the behavior of the effective mass square $m^2_{\rm eff}$ and have a qualitative explanation. In Fig. \ref{mSquareNegativeAlpha}, $m^2_{\rm eff}$ as a function of the coordinates is plotted for various $\alpha<0$. From the left panel of the figure, one can see that when $a \leq a_{\rm min}=1/2$, $m^2_{\rm eff}$ is always positive everywhere for  any value of $\alpha<0$. This presents physically reason why there is no tachyonic instability for low spin $a \leq a_{\rm min}$. However, when $a>a_{\rm min}$, $a=0.9$ for example as shown in the right panel of the figure, we find that $m^2_{\rm eff}$ will become negative near the horizon, which explains the emergence of the tachyonic instability for sufficiently negative $\alpha$. Moreover, similar to the case with $\alpha>0$, as $|\alpha|$ is increased for a fixed spin,  $m^2_{\rm eff}$ becomes more and more negative at some fixed position, which explains the appearance of more and more violent instabilities.

	Besides the stability investigation in the final wave dynamics, we have used Prony's method and read off the QNMs in the time evolution of the perturbation in the ringdown phase for  $\ell=1$, see the results in Table \ref{QNMs}. We find that the result is simpler than the case with $\alpha>0$ mentioned above.  For a chosen $a$, with the decrease of $|\alpha|$, both $\Re \omega$ and $|\Im \omega|$ decrease. Differences in QNM frequencies from the GR situation show that in sEGB theory the characteristic sound is special. We expect that future precise detections of  GWs can observe the signature of the sEGB theory from its special sound frequency. Once again, due to the ``mode-mixing mechanism" \cite{Zenginoglu:2012us,Burko:2013bra,Thuestad:2017ngu}, one can expect that instability will also occur for higher  multipoles, and moreover even and odd modes will give roughly comparable contributions to the instability \cite{Dima:2020yac}.

	\section{Summary and Discussions}

	In this work, we have studied the massless scalar perturbation on the background of the Kerr black hole in sEGB theory. We have overcome the nonseparability problem and the ``outer boundary problem" by employing the hyperboloidal compactification technique to perform carefully the numerical computations and have got object pictures on the time evolution of the scalar field  perturbations. In addition to the information about the stability, we have also obtained the sound characteristics of the quasinormal ringing in the sEGB theory. For negative couplings between the scalar field and the GB term, we have presented objective support on the existence of the minimum rotation to break the Kerr stability, no matter how strong the coupling $\alpha$ is. This minimum $a$ requirement does not exist for positive $\alpha$s. For  fast rotating black holes, it requires smaller coupling $|\alpha|$ to trigger
	the instability. This result holds for negative couplings when the rotation is above a minimum value.  The stronger coupling $|\alpha|$ will bring the break down of the original Kerr configuration to happen earlier and the resulting instability becomes more violent.  We have also examined carefully the characteristic sound of the quasinormal ringing of the scalar perturbation outside the original Kerr black hole in the sEGB theory. We have unveiled the special sound once there is a coupling between the scalar field and the GB curvature and when it becomes stronger. The fine frequency difference uncovered can be used  as imprints to confirm the sEGB theory in the future precise GW detections.

	In addition to presenting objective pictures of the wave propagation in the sEGB theory, we have also examined the physical reason about the occurrence of the instability of the Kerr black hole once the coupling between the scalar field and GB term is strong enough. We have found that a big enough $|\alpha|$ always leads to a negative  effective   mass square, which triggers the  tachyonic instability and destroys the original Kerr background configuration in the sEGB theory.  With the increase of $|\alpha|$, the effective mass square  becomes more negative, which results in a more violent instability. For negative couplings, we have found the reason behind the minimum rotation in order to accommodate the instability, and the effective mass square is always positive no matter how strong  the coupling strength is.

    In this paper, we have concentrated ourselves on the massless scalar field perturbation. It is of great interest to generalize the discussion to the massive scalar field. For the massive scalar field, in addition to the tachyonic instability, there is the possibility for the superradiant instability to happen. There will be richer physics involved in studying the dynamics and the fate of the Kerr black hole. Furthermore, in this study we have only paid attention to the probe limit, the linear perturbation in the  Kerr black hole background. It is intriguing to extend our study to include the backreaction of the scalar field, which can present us the process of the formation of the scalar hair and the final spacetime of the hairy black hole.  Another interesting direction is to extend our study to consider the coupling of the scalar field to other topological curvature invariants beyond the Gauss-Bonnet term we considered here, such as Chern-Simons term \cite{Gao:2018acg}, Lovelock term, etc.

	\begin{acknowledgments}

		This work is partially supported by the National Natural Science Foundation of China (NNSFC) under Grant Nos. 12075202, 11675145,  11975203, and 12075207. J.S. is supported by ANID Chile through FONDECYT Grant  No.1170279.
	\end{acknowledgments}