A Data-Analytic Evaluation of Smoking Cessation via Online Hypnosis

Abstract: In this article we report the observations of our field study which investigates the success rates of online hypnosis for smoking cessation. Following an international call in Germany, Austria and Switzerland for subjects willing to stop smoking $83$ participants contributed to this study, which took place during the lockdowns of the COVID-19 pandemic. All subjects received two online hypnosis sessions in small groups and additionally got two MP3 audio files for invididual intensification. The participants sent back several questionnaires, from which we present the evaluation results, and we discuss the problems observed during the study. It turned out that $86\%$ of the participants are smoke-free after $6$ months as long as they participated according to the intervention protocol. Among all smokers after six months we have observed a reduction of $32\%$ in cigarette consumption. This report shall be seen as a statistical evaluation of the collected data and not as the report of a fully pre-designed study in a narrower sense.

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

According to the German Federal Ministry of Health \cite{bmg-tote} smoking is the most avoidable health risk in Germany: each year in average $127.000$ deaths in Germany are counted which are related to tobacco consumption. In 2019 about $458.000$ patients have received hospital treatments due to smoking-specific diseases like lung cancer or COPD, which is an increase of $18\%$ compared to the year 2010, see German Federal Office of Statistics \cite{destatis1}. Governmental programs have been established in order to decrease the rate of smokers. Morever,
during the last decade anti-smoking laws in Germany and across Europe have been tightened in order to protect non-smokers in public areas like train stations, airports or the  interior of restaurants. This seems to imply an increasing demand of smokers in supporting them to quit smoking, compare with German Federal Ministry of Health \cite{bmg-rauchfrei}.
\par
Many smokers seem to think about quitting smoking, but most of them are not successful without any additional help or intervention and there are high fall-back rates after some time without smoking. According to Hughes, Keely and Naud \cite{hughes-keely-naud} and Giovino et al. \cite{giovino} only $3\%-5\%$ of smokers, who try to quit on their own, are successful after $6$ to $12$ months.
There is a big variety of methods to support smokers in their wish to become abstinent of cigarettes, e.g., single and group consultations, pharmalogical interventions, behaviour and/or psychological therapies and also intervention by hypnosis. In Germany the federal government proposed the program \textit{Rauchfrei} (smoke-free) which should decrease the rate of smokers, see also \cite{bmg-rauchfrei}. Amongst others, an evaluation and good summary of the intermediate results can be found, e.g., in Wenig et al. \cite{wenig}.
\par
In this article we are interested in the effectiveness of smoking cessation with the help of hypnosis, a classical  method for quitting smoking. Typically, smoking intervention by hypnosis uses a concept where people meet a hypnotherapist in person for one or more sessions. During the COVID-19 pandemic meetings in person were restricted in a very extended way laws. In many areas of daily life meetings went \textit{online}. While online hypnosis was offered by a rather limited number of hypnotists before the pandemic, it seems that nowadays many hypnotists started to offer online hypnosis sessions  for smoking cessation or weight loss. During an online session via online conference tools like Skype, Zoom, etc., the patients stay at home and turn on their webcams, and the hypnotist may perform his session almost in the same way as in a real-life session. Nonetheless, there are obviously still differences between online hypnosis sessions and sessions in real life. On the one hand side, the direct contact between hypnotist and subject is lost, which may reduce the chance of success; on the other hand side, subjects are staying in their comfort zone at home, where they may relax even easier than in an external room provided by the hypnotist. This leads to the question of effectiveness of online hypnosis, and
this question was the starting point for our own field study, in which we want to evaluate the \textit{effectiveness of online hypnosis for smoking cessation}.
For this purpose, we made a call for smokers wishing to stop smoking by online hypnosis intervention. Our hypnosis intervention protocol consists of two online hypnosis sessions and individual post-session treatment by audio files. Effectiveness is examined by questionnaires after each session, and additionally after three and six months. The protocol will be described in detail in Section \ref{sec:method}.
At this point we remark that we do \textit{not} aim on comparing the success rate of hypnotic intervention with other intervention types, but we have the following goals:
\begin{itemize}
\item estimate the \textit{success rate by intervention via online hypnosis},
\item determine the \textit{main criteria which influence the success rate},
\item and detect \textit{possible problems} which may lead more likely to flops.
\end{itemize}
Let us provide some results on effectiveness of smoking intervention by hypnosis. Barrios \cite{barrios} summarized the results of several studies concerning effectiveness of different intervention methods for smoking cessation. Averaging over several studies, where participants received hypnosis intervention for smoking cessation (six sessions in average), he came to a success rate of $93\%$. Other therapies seem to be much less effective: $72\%$ of smokers, who received a behavioural therapy, were successful after 22 sessions in average, while only $38\%$ of smokers, who received a psychoanalytic treatment over several years, were successful after 600 sessions in average. Amongst others, we want to mention the following articles which evaluated the success rate of hypnotic intervention for smoking cessation:
Johnson and Karkut \cite{johnson-karkut} came to a $86.5\%$ success rate after three months when combining hypnosis and aversion techniques. In particular, there was no gender-specific difference in the success rate.  Barber \cite{barber} calculated a success rate of $90.6\%$ when using an integrated approach combining hypnotic methods and rapid smoking treatment protocol. The article of Elkins and Rajab \cite{elkins2004} presents preliminary data of a three-session hypnosis intervention for smoking cessation, which comes close to the method of hypnotic intervention used in our article; they came to a success rate of $81\%$. However, the picture of effectiveness of hypnotic intervention for smoking cessation is still mixed: e.g.,
the article of Ahijevych, Yerardi and Nedilsky \cite{ahijevych}  ended up with a success rate of less than $25\%$.  Nonetheless, the critical article of Green and Lynn \cite{green-lynn} compared $59$ studies and they came to the conclusion to classify hypnosis as a ``possibly efficacious'' treatment for smoking cessation.
\par
This paper is organized as follows: in Section \ref{sec:subjects} we explain how we found participants and we give a short statistical overview on the composition of the participants. In Section \ref{sec:method} we present the protocol of our study, including a detailed description of the methods used for the hypnotic intervention for smoking cessation. Finally, in Section \ref{sec:results} the main results of our study are presented, and the results and problems observed are discussed in Section \ref{sec:discussion}.

\section{Subjects}
\label{sec:subjects}

The participants for this field study were acquired by calls in social media (Facebook, Twitter, Instagram), regional newspapers (Dorstener Zeitung, Passauer Neue Presse, Grazer Woche), and via radio channels (UnserRadio, Radio Galaxy). Finally, $99$ subjects from Germany, Austria and Switzerland were taken into account for this study, consisting of $56$ ($57\%$)  female and $43$ ($43\%$) male subjects. All interested subjects were accepted for participation unless severe psychological problems in the past were given, which led to an exclusion of participation.
The participants had to pay a small nominal fee in order to cover the study's expenses (for this study, \textit{no} fundings were  requested)  and to prevent participation of people, who just want to ``try'' without any serious intention.
\par
The average age of the participants was $42$ (median: $40$), while the minimum age was $23$ and the maximum age was $70$ with a standard deviation of $11.4$. Furthermore, $52\%$ of the participants were married, while $16\%$ were divorced. We also remark that $36\%$ of the participants were suffering from overweight.
\par
The average daily amount of cigarettes smoked by the participants \textit{before} the intervention was $19.4$ cigarettes per day (median: $20$), while the minimum was $4$ cigarettes per day and the maximum was $50$ cigarettes per day. The average number of years of smoking was $24.9$ years (median: $23$), with a range between $5$  and $50$ years. The participants' strength of wish to stop smoking and their self-confidence in the intervention is sketched in Figure \ref{fig:wish}.
\begin{figure}[H]
\includegraphics[scale=.25]{Wunschaufzuhoeren.jpeg}
\hspace{0.5cm}
\includegraphics[scale=.25]{Sicherheit0.jpeg}
\caption{Participants' strength of wish and self-confidence to quit smoking.}
\label{fig:wish}
\end{figure}
As one can see, the majority has a strong or very strong wish to quit smoking; nonetheless, around $22$ subjects (strength of wish at level $5$ or $6$ or self-confidence level $1$ or $2$) seem to participate without having a strong desire to stop smoking and/or have rather less self-confidence in the intervention. One third of the particpants have an intermediate confidence in the hypnosis intervention. It is clear that the subject's collaboration in some way is indispensable for the success, but we want to study the main connections between sufficient collaboration/motivation and success.
\par
We also remark that  $91\%$ among the $99$ participants reported that they have failed in previous unassisted attempts to quit smoking permanently. Moreover, $23\%$ of the participants already had general experience with hypnosis.
\par
Finally, we want to give insights in which situations the subjects were smoking before the intervention. Each participant had to give a grade (on a scale from $1$ to $5$) for a given situation, in which he smokes very often (grade $1$) or very rarely (grade $5$).
Figure \ref{fig:smoking-situations} lists the favourite reasons, where participants gave grade $1$:
\begin{figure}[H]
\includegraphics[scale=.3]{Rauchgruende.jpeg}
\caption{Participants' main situations for smoking.}
\label{fig:smoking-situations}
\end{figure}
E.g., for $46.5\%$ of the participants it is quite hard \textit{not} to smoke when consuming alcohol. Furthermore,  $63\%$ of the participants tend to smoke when drinking alcohol, after standing up in the morning or when having stress.

The participants were splitted up into different groups according to their stress level, psychological problems in the past (these persons received an induction leading into a lighter state of trance in order to avoid any problems occuring from the past) and other secondary factors like the ones mentioned above. The distribution of the participants' stress level and previous experience with hypnosis and/or psychotherapies is sketched in \mbox{Figure \ref{fig:stresslevel}.}
\begin{figure}[H]
\includegraphics[scale=.24]{Stresslevel.jpeg}
\hspace{0.1cm}
\includegraphics[scale=.24]{Erfahrungen.jpeg}
\caption{Participants' stress level (0=very low, 10=very strong) and experiences with hypnosis or psychotherapies.}
\label{fig:stresslevel}
\end{figure}

Although we tried to de-randomize the group structure, we mention that it was only possible in a very limited way to adapt the suggestions given to the different groups during  the hypnosis sessions.
\par
Finally, we remark that we even had $15$ participants more, who attended the first hypnosis sessions but did \textit{not} send in the questionnaire after this session or who did \textit{not} attend the second session; those particpants were excluded from the evaluation.

\section{Method}
\label{sec:method}

In this section we describe the intervention protocol, which consists of the following four main points: introduction to hypnosis, a first hypnosis session, individual treatment with MP3 audio files and a second hypnosis session. We will explain these stages in detail in the following and we also give a comment on the data collection. Due to the COVID-19 pandemic, all hypnosis sessions took place online via \textit{Zoom}. The participants were splitted up into different rather homogeneous groups (minimal group size was $5$, maximum group size was $20$) by taking into account different stress level, psychological problems in the past and other secondary factors.
\par
All online meetings and hypnosis sessions were held by the authors of this article, who are both certified hypnocoaches.

\subsection{Introduction to Hypnosis}

All participants were invited to join an online meeting via Zoom, where we explained our study, formulated the goals and gave an explanatory introduction to hypnosis. It was our main goal to reduce or even eliminate false opinions on hypnosis (e.g., originating from TV shows etc.) and to explain what happens during a hypnosis session. This education was very important in order to prepare all participants as good as possible. Moreover, all participants received guidelines for the rules of conduct for the hours before and after the hypnosis sessions (e.g., we remarked that people should not drink coffee just before the sessions or have to strictly avoid testing after the session if smoking still ``works'').

\subsection{Session 1}

In the first session we used a classical hypnosis approach for anti-smoking intervention consisting of the following parts:
\begin{enumerate}
\item Classical induction with encouraging people to quit smoking.
\item $7+2$ deepening: this deepening uses the approach that people, in general, can only think of $7$ plus or minus $2$ things at the same time. Subjects are requested to think at different parts of the body step by step. This should overload the consciousness and help to improve the relaxation.
\item Dissociation of mental blockades: this part aims on preparing the subjects to get rid of mental problems  and to become open for the following part of suggestions. The subjects went on a virtual  journey to a waterfall, under which all problems and troubles in their mind and soul could be washed away.
\item Suggestion booster: this part shall boost the upcoming suggestions.
\item Commitment for stopping to smoke: in this part the subjects are prepared to stop smoking; they are encouraged to get rid of smoking.
\item Posthypnotic suggestions: in this part the main suggestions are given which aim on the non-necessity of smoking, better feelings and emotions without cigarettes, and so on. Suggestions were kept quite general in order to reach all participants.
\item Integration of the goal: after having implemented the anti-smoking suggestions some further suggestions were given in order to deepen the suggestions and to intensify them.
\item Closing the session: people were led out of trance.
\end{enumerate}
The presented structure of this session follows a well-known scheme which is widely used by hypnotherapists for smoking cessation.

\subsection{MP3 Audio Files \& Post-Session Treatment}

Three days after the first session all participants received two MP3 files for individual treatment and intensification. The first audio file is a classical $15$-minutes short hypnosis with suggestions for deepening the hypnotic anti-smoking suggestions. Furthermore, the participants got a subliminal MP3 file with anti-smoking suggestions. The participants were encouraged to listen to these MP3's whenever  the impulse to smoke increases and as often as they wanted (e.g., during lunch breaks, in the evening, etc.), but \textit{at least once per day during the first two weeks} unless the impulse to smoke is already gone. This task is an important part of the intervention protocol, which is requested to be followed by the participants.

\subsection{Session 2}

The second hypnosis session took place \textit{one week after the first session}. It consisted of the following parts:
\begin{enumerate}
\item Classical induction with eye fixation.
\item Brain twister deepening: this deepening sequence aims on creating some confusion, which overloads the subject's consciousness resulting in an even deeper state of relaxation.
\item Dissociation of mental blockades:  analogously to the first session the subjects were asked to  imagine an ice block consisting of all problems and troubles in their mind, which is then melted away.
\item Suggestion booster: this part shall boost the upcoming suggestions.
\item Commitment for stopping to smoke and reviewing the personal results so far: subjects were encouraged to look back on their (maybe partial) success and to improve their trust in a successful smoking cessation.
\item Posthypnotic suggestions: these suggestion aim once again on better feelings and emotions without cigarettes, and so on.
\item Integration of the goal: some further suggestions were given in order to deepen the anti-smoking suggestions and to intensify them.
\item Closing the session: people were led out of trance.
\end{enumerate}
The structure of this session follows also a well-known scheme which is widely used by hypnotherapists for a second hypnosis session in the context of smoking cessation. We note that the second session should \textit{not} be seen as a ``second chance'', but as an intensification of the first session which should lead to a sustainable smoke-free life. The participants were informed about this fact during the introduction event.

\subsection{Data Collection}

The participants were asked to submit questionnaires at the follwing stages:
\begin{enumerate}
\item Before the first session for getting information about the smoking behaviour and to exclude psychological problems.
\item Three days after Session $1$.
\item One week after Session $2$, that is, two weeks after Session $1$.
\item Approx. six weeks after Session $2$.
\item Approx. six months after Session $2$.
\end{enumerate}

\section{Results}
\label{sec:results}

We present the main results of our study in this section which is structured as follows. First, we identify three different motivation-related groups of participants identified according to the observations we have made during the intervention process (in particular, during the hypnosis sessions). Afterwards we present the success rates of the anti-smoking intervention considered from different points of view. Finally, we compare the groups and give insights into each group, where we detect some factors which are in correlation with success or failure of the intervention.

\subsection{Individual Participation Related Group Assignment}
\label{subsec:groups}

One of the main observations during the intervention process was that the subjects attended the hypnosis sessions quite differently. While there was a group of subjects who really wanted to quit smoking and put effort into the anti-smoking process, there was another group of subjects who did \textit{not} seem to participate with the required seriousness or with sufficient effort. The later group contains a subgroup of subjects who did \textit{not} attend the hypnosis sessions in a serious manner observed from an objective point of view;  e.g., people in this subgroup  stood up during the hypnosis session or pet their dogs/cats, and so on.
These people definitely were \textit{not} in any state of trance and we put them into a \textit{Placebo effect group}, where we do \textit{not} expect any success due to the effects of hypnosis. This group shall serve as a control group, where people possibly quit smoking due to solidarity or by their own effort.
\par
The composition of the remaining  subjects (that is, all participants not belonging to the \textit{Placebo effect group}) was investigated closer according to their degree of following our intervention protocol concerning post-session treatment. It turned out that a good indicator for the degree of following our guidelines is given by the \textit{intensity of the individual post-treatment} starting after Session $1$. The participants were asked to listen to the audio files \textit{at least once per day} for the first two weeks after Session $1$ (unless there is no impulse to smoke any more), depending on how difficult it is to resist smoking. While there were many subjects, who listened to the audio files almost daily, also many subjects did listen to the audio files only very rarely. We have identified this difference as a good indicator whether subjects participated in a reasonable motivated manner, and therefore we classify the participants accordingly: since there were $12$ days between the receipt of the audio files and day $14$ after Session $1$, we classify those participants as \textit{motivated} who finally quitted smoking after Session $2$ and are still non-smoker after six months or who have listened to the audio files at least $12$ times during the first two weeks. In other words, those participants were immediately successful after Session $2$ and/or they have at least tried to put enough own effort into the intervention process.
\par
As a summary, we have split the participants up into one of the following (disjoint) groups:
\begin{itemize}
\item \textit{Placebo effect group:} participants, who did \textit{not} attend the hypnosis session in a serious manner.
\item \textit{Lazy group:}  participants, who did attend the hypnosis session in a mostly serious manner but did \textit{not} follow the intervention protocol as requested.
\item \textit{Motivated group:}  participants, who did attend the hypnosis session in a serious manner and did follow the intervention protocol as requested.
\end{itemize}
The group sizes are as follows:
\begin{center}
\tablinesep=2ex\tabcolsep=10pt
\begin{tabular}{|c|c|c|}
\hline
Group & Number of participants & Percentage\\
\hline
\hline
\textit{Motivated group} & $28$ & $28.3\%$ \\
\textit{Lazy group} & $55$ & $55.6\%$\\
\textit{Placebo effect group} & $16$ & $16.1\%$\\
\hline
Total & $99$ & $100\%$\\
\hline
\end{tabular}
\end{center}

For the evaluation of the effectiveness of online hypnosis for smoking cessation the \textit{Motivated group} is essentially most important.
In Section \ref{subsec:motivated-vs-lazy} we demonstrate the impact when separating   the \textit{Motivated group} and the \textit{Lazy group} by a lower number of minimal audio files listenings, and in Sections \ref{subsec:discussion-decomposition} and \ref{subsec:discussion-M-L} we discuss the question of sufficient motivation regarding the separation between the \textit{Motivated group} and the \textit{Lazy group}.

\subsection{Success Rates}
\label{subsec:success-rates}

All participants reported three days after Session $1$ whether they have smoked since Session $1$. One week after Session $2$ they reported whether they have smoked since Session $2$. Approx. six weeks and approx. three months after Session $2$ the participants reported whether they are still/again smoking or not. The evolution of non-smoker rates for all three groups are sketched in Figure \ref{Nonsmoker-lazy-main}.
\begin{figure}[H]
\includegraphics[scale=.33]{Success-Rates4.jpeg}
\caption{Evolution of non-smoker rates by group.}
\label{Nonsmoker-lazy-main}
\end{figure}

We remark that one person in the \textit{Placebo effect group} does \textit{not} smoke after six months, while two participants in the \textit{Lazy group} are \textit{not} smoking after six months.

\subsection{Amount of Cigarettes Consumption}
\label{subsec:amount-cigarettes}

The main goal of the anti-smoking hypnosis intervention was, of course, to eliminate tobacco consumption. Nonetheless, we take a closer look on the smoking behaviour after six months in the different groups \textit{among smokers}. The average amount of daily cigarettes before the intervention process and after six months \textit{among smokers}  is summarized by group  in the following table, where the corresponding standard deviation is noted in brackets:

\begin{center}
\tablinesep=2ex\tabcolsep=10pt
\begin{tabular}{|c|c|c|}
\hline
Group & Mean before intervention  & Mean after $6$ months  \\
\hline
\hline
\textit{Motivated group} & $21.8$ ($9.6$) & $14.8$ ($17.3$)\\
\textit{Lazy group} & $18.4$  ($6.5$) &  $12.7$ ($6.1$)\\
\textit{Placebo effect group} & $24.7$ ($7.8$) & $17.9$ ($4.4$) \\
\hline
All groups & $19.9$ ($7.3$) & $13.9$ ($7.0$)\\
\hline
\end{tabular}
\end{center}

The differences of the individual  average number of daily cigarettes \textit{before the intervention} with the individual average number of daily cigarettes \textit{after six months} are sketched in Figure \ref{cigarettes-decrease}.

\begin{figure}[h]
\includegraphics[scale=.33]{Rauchrueckgang-absolut3neu.jpeg}
\caption{Decrease in the average number of daily cigarettes among smokers after six months.}
\label{cigarettes-decrease}
\end{figure}

Among all smokers we observe an average decrease of $6.3$ cigarettes per day after six months .
\par
We also present the estimated densities of the groupwise differences after two weeks and after six months which are drawn in Figure \ref{fig:densities}.
\begin{figure}[H]
\includegraphics[scale=.23]{densities2-neu.jpeg}
\hspace{.5cm}
\includegraphics[scale=.23]{densities-neu.jpeg}
\caption{Estimated densities of the decrease in the average number of daily cigarettes among smokers.}
\label{fig:densities}
\end{figure}

The \textit{relative decrease} in daily cigarettes consumed before the intervention and after six months \textit{among smokers} is sketched in Figure \ref{cigarettes-decrease-prozentual}.

\begin{figure}[h]
\includegraphics[scale=.35]{Rueckgang-prozentual-neu.jpeg}
\caption{Relative decrease in daily cigarettes among smokers after six months.}
\label{cigarettes-decrease-prozentual}
\end{figure}

The evolution of the average relative decreases in daily cigarettes \textit{among smokers} is sketched in Figure \ref{cigarettes-average-evolution}.
\begin{center}
\begin{figure}[H]
\includegraphics[scale=.38]{Rauchrueckgang-Entwicklung.jpeg}
\caption{Average decrease of daily cigarettes  \textit{among} smokers.}
\label{cigarettes-average-evolution}
\end{figure}
\end{center}

The  boxplots in Figure \ref{cigarettes-decrease-prozentual-boxplot} compare the location of the individual relative decreases by group:

\begin{figure}[h]
\includegraphics[scale=.25]{Zigarettendifferenz3-Boxplot.jpeg}
\hspace{0.5cm}
\includegraphics[scale=.25]{Zigarettendifferenz-Boxplot.jpeg}
\caption{Relative decrease in daily cigarettes among smokers.}
\label{cigarettes-decrease-prozentual-boxplot}
\end{figure}

Among all participants who are still smoking after six months we have an average decrease of $31 \%$ in tobacco consumption.

\subsection{Group Comparision}

\subsubsection{Self-confidence and Strength of Wish}
\label{subsec:self-confidence}

We compare two important \textit{a-priori} indicators for the motivation of the participants: the self-confidence score and the strength of wish to become non-smoker; see Figure \ref{fig:wish}.
\par
Comparing the self-confidence score by group on a scale from $1$ (very low) to $5$ (very high)
\textit{before} the first hypnosis session  (see Figure \ref{fig:wish}) yields the following picture:

\begin{center}
\tablinesep=2ex\tabcolsep=10pt
\begin{tabular}{|c|c|c|}
\hline
Group & Mean & Standard deviation \\
\hline
\hline
\textit{Motivated group} & $3.46$ & $0.96$ \\
\textit{Lazy group} & $3.24$ & $0.94$ \\
\textit{Placebo effect group} & $3.56$ & $0.72$ \\
\hline
\end{tabular}
\end{center}

The strength of wish score on  a scale from $1$ (very low strength) to $10$ (very high strength) from Figure \ref{fig:wish} by group is summarised in the following table:

\begin{center}
\tablinesep=2ex\tabcolsep=10pt
\begin{tabular}{|c|c|c|}
\hline
Group & Mean & Standard deviation \\
\hline
\hline
\textit{Motivated group} & $8.18$ & $1.42$ \\
\textit{Lazy group} & $8.55$ & $1.54$ \\
\textit{Placebo effect group} & $8.19$ & $1.91$ \\
\hline
\end{tabular}
\end{center}

The self-confidence \textit{after} Session $1$ -- now on a scale from $1$ (very low) to $10$ (very high) -- is as follows:

\begin{center}
\tablinesep=2ex\tabcolsep=10pt
\begin{tabular}{|c|c|c|}
\hline
Group & Mean & Standard deviation \\
\hline
\hline
\textit{Motivated group} & $7.75$ & $2.10$ \\
\textit{Lazy group} & $6.36$ & $2.41$ \\
\textit{Placebo effect group} & $5.31$ & $2.27$ \\
\hline
\end{tabular}
\end{center}

We remark that the scale change in the self-confidence score was made in order to get better insights, although comparing the scores before and after Session $1$ is now a bit vague.

\subsubsection{Motivated Group vs. Lazy Group}
\label{subsec:motivated-vs-lazy}

We separated the subjects in the \textit{Motivated Group} and the \textit{Lazy Group}  by considering the number of audio listenings within the first two weeks of the intervention.
The participants were requested to listen to the audio files at least once per day  during this period, hence at least $12$ times. A legitimate question is now what happens if we replace this critical number by a lower number. If we put all participants not belonging to the \textit{Placebo effect group} into the \textit{Motivated group} if they have \textit{not} smoked since Session $2$ any more (maybe with a lower number of listenings because it was not hard for them to resist smoking) or have listened to the audio files at least $N$ times, then we obtain a non-smoker rate after six months depending on $N$ as sketched in Figure \ref{fig:success-listenings}.
\begin{figure}[H]
\includegraphics[scale=.28]{Success-Rate-dependent-hearings-neu2.jpeg}
\caption{Dependency of the success rate on the number of audio file listenings during the first two weeks.}
\label{fig:success-listenings}
\end{figure}
Observe that the hightest jump occurs at $N=12$.

\subsubsection{Lazy Group vs. Placebo Effect Group}
\label{subsub:compare}

As we have seen in Figure \ref{Nonsmoker-lazy-main}  there is some significant difference between the \textit{Lazy group} and the \textit{Placebo effect group}  at the beginning of the intervention which disappears towards the end of the observation period.
We illustrate this observation with the help of  two indicators: the decrease in daily cigarettes and the subjects' possible observation that cigarettes taste bad.
\begin{enumerate}
\item A Kolmogoroff-Smirnov test suggest at a $p$-value of $0.046$ that the average decrease in daily cigarettes after six weeks in the \textit{Lazy group} is bigger than in the \textit{Placebo effect group}. However, after six months this significance is disappearing, that is, we can \textit{not} reject the conjecture that the average number of daily cigarettes after six months  follow the same distribution in both groups.
\item
The participants were asked at different stages of the intervention whether cigarettes taste the same as before the intervention or not. Figure \ref{cigarettes-bad-taste} shows the results:
\begin{center}
\begin{figure}[H]
\includegraphics[scale=.30]{Zigarettengeschmack-neu2.jpeg}
\caption{Rate of participants who do \textit{not} smoke or for whom cigarettes are tasting bad.}
\label{cigarettes-bad-taste}
\end{figure}
\end{center}
We observe a different behaviour during the first two stages of the intervention process, but the gap between \textit{Lazy group}  and \textit{Placebo effect group} closes after six months. A $95\%$ confidence interval for the probability that a random member in the \textit{Lazy group} does not smoke or finds cigarettes non-tasteful after Session $1$ is given by the Clopper-Pearson interval $[0.755,  0.947]$. The observed value of $0.38$ in the \textit{Placebo effect group} lies far outside of this interval, which affirms the conjecture that both groups behave differently during the first stage. After six months the values are so close together such that we can \textit{not} reject the conjecture that both groups now follow the same law.
\end{enumerate}

\subsection{Group Specific Observations}

In this section we investigate group specific criteria for success or failure of the smoking intervention. We use Fisher's exact test in order to detect possible stochastic dependencies of different factors with the success.

\subsubsection{Motivated Group}

Consider the main situations for smoking listed in Figure \ref{fig:smoking-situations}. Among the main reasons only the combination of smoking when hungry seems to have a negative effect in this group: independence of success after six months and smoking habits when hungry can be rejected with a $p$-value of $0.011$. Furthermore, we detect some light
 negative dependency with being nervous or having fear before the hypnosis sessions with a $p$-value of $0.094$.
Since there are only four smokers after six months in this group, the statistical dependencies are a bit uncertain and should be taken with caution only.

\subsubsection{Lazy group}
\label{subsec:lazy}

In this group we consider the success of being abstinent \textit{after the first two weeks}. We can \textit{reject} stochastic independence of the success after two weeks with the following factors:
\begin{center}
\begin{tabular}{|c|c|}
\hline
Factor & $p$-value \\
\hline
\hline
Smoking when hungry & $0.012$\\
Smoking when having lack of concentration & $0.014$\\
Individual stress level & $0.015$ \\
After standing up in the morning & $0.040$ \\
Relationship & $0.037$\\
Influence from other persons &  $0.054$ \\
Smoking when being angry & $0.074$ \\
\hline
\end{tabular}
\end{center}

\subsubsection{Placebo Effect Group}

In the \textit{Placebo effect group} we did \textit{not} detect any specific factors/circumstances which let people tend to attend the hypnosis sessions inappropriately.

\subsection{General Summary}
\label{subsub:general-result}

From Subsections \ref{subsec:success-rates} and \ref{subsec:amount-cigarettes} we can draw the following summary: we observe  that $85\%$ among \textit{all} participants do \textit{not} smoke after six months or have at least reduced their daily tobacco consumption.

\section{Discussion}
\label{sec:discussion}

\subsection{Summary of the Study}
This study had the aim to investigate the effectiveness of smoking cessation via \textit{online hypnosis} in small groups, which has \textit{not} been studied as far as the authors know. For this purpose, participants were acquired by an international call over social media, regional newspapers and radio channels. All interested applicants were allowed to participate as long as there were no contraindications (e.g., due to psychological problems in the past).
The participants received two online hypnosis sessions in small groups and also audio files for individual post-treatment and intensification. Furthermore, they were briefed on the intervention process and their required collaboration in advance.
 At different stages the participants reported whether they are still abstinent of cigarettes and tobacco or not.
\par
In order to allow an adequate evaluation of the effectiveness we had to split all participants up into different groups according to their willingness/motivation to contribute to the smoking cessation process: a \textit{Placebo effect group}, a group with rather less or intermediate level of motivation which did \textit{not} satisfy the expected level of collaboration (called  \textit{Lazy group}), and a group where the motivation satisfied our intervention protocol (called  \textit{Motivated group}). The study resulted in a $86\%$ abstinency rate in the  \textit{Motivated group} after six months, while most of the participants in the other groups are  smoking again after $6$ months. Moreover, \textit{among all smokers} there is an average decrease of $31\%$ of tobacco consumption after six months since the beginning of the intervention.

\subsection{Decomposition of Participants}
\label{subsec:discussion-decomposition}

Let us discuss the decomposition of the participants into different groups regarding their motivation. While it was obvious whether a subject belongs to the \textit{Placebo effect group} or not due to his manner of attending the hypnosis sessions, the remaining subjects (beyond the \textit{Placebo effect group}) formed still a quite unclear picture. Our observation was that some subjects participated well-motivated while others didn't. This was the starting point to consider motivated subjects separately and to detect an objective indicator which separates subjects according to their motivation.
The assignment of these subjects either to the \textit{Motivated group} or to the \textit{Lazy group} was, of course, less obvious at first glance. Therefore, we demanded for members to be in the \textit{Motivated group} that our guidelines (in particular, to listen to the audio files at least once per day during the first two weeks after Session $1$) are followed  during the whole intervention process.
\par
For this purpose, we have found an objective indicator by considering the number of audio listenings within the first two weeks of the intervention. Continuous post-session treatment by listening to audio files is a good support to strengten the abstinence of tobacco and to deepen the suggestions received during the hypnosis sessions. From this point of view, participants, who listened to the audio files rarely during the first two weeks of intervention, obviously seemed  to have some lack of motivation.
Setting the critical number of audio listenings to $12$ is justified by different reasons: first, $12$ listenings correspond to an average number of listenings of once per day during the first two weeks; second, Figure \ref{fig:success-listenings}  illustrates how the number of listenings influences the success rate, where the highest jump occurs at $N=12$, a reasonable change point to choose.
\par
From our point of view, subjects in the \textit{Lazy group}  seem to be rather to lazy to follow the complete intervention protocol; instead, they maybe hoped that the hypnosis sessions does all the work for them and no own effort is required in order to strengthen the effects of hypnosis and to complete the intervention process, while most of the people in the \textit{Motivated group} seemed to have a stronger motivation.
We note that not all non-smokers after six months in the \textit{Motivated group} listened to the audio files at least $12$ times; these participants were already completely abstinent of tobacco immediately after Session $2$, since the hypnosis seemed to work very well for them. Of course, motivation cannot be classified by the number of audio listenings solely, but -- in absence of other reasonable, objective indicators -- this number is at least one obvious indicator for the degree of  the subjects' personal motivation and collaboration. One critical point, of course, is that we have an \textit{a-posteriori} group decomposition and there is the question which factors influence the effectiveness and motivation beyond the individual audio listenings. Further work has to be done in this direction. For a discussion regarding differences between the \textit{Motivated group} and the \textit{Lazy group}, we refer to Subsection \ref{subsec:discussion-M-L}.
\par
Recall that we have split up the participants into small groups for the single hypnosis sessions. We note that we did \textit{not} observe any different behaviours concerning which hypnosis session has been visited by the participants.
\par
Let us also remark that our big call for participants together with the very low  participation fee of $25$\,\euro{}  (we note that hypnosis intervention of our proposed form starts from approx. $200$\,\euro{} in Germany) presumably led to participation of many people, who rather wanted to get the experience of hypnosis without putting too much \textit{own} effort into the whole intervention process. This could be an explanation why the \textit{Lazy group} has $51$ members and the \textit{Motivated group} consists of $28$ members only.
\par
Overall, the presented reasoning underlines the necessity  of motivation and collaboration which led to the presented decomposition of the participants into  different groups.

\subsection{Success Rates}

As one can see in Figure \ref{Nonsmoker-lazy-main} the motivation and willingness to put some own effort into the intervention process is \textit{essential} for the success to become non-smoker. In particular, observe the higher rate of non-smokers in the \textit{Motivated group} one week after Session $2$. We interpret this higher rate after Session $2$ as the necessity of an additional hypnosis session for deepening the suggestions and improving the effects: some smokers were not able to stop smoking completely after Session $1$ but Session $2$ let them finally stop.
 However, the non-smoker rates in the \textit{Lazy group} are decreasing quickly and are already quite low after six weeks.
 \par
Our main conclusion is that the hypnosis sessions in small groups lay the base for a successful smoking intervention but  some individual post-session treatment is inevitable in most cases for a sustainable effect. In particular, sufficient motivation and putting some own effort into the intervention process is indispensable in most cases.

\subsection{Amount of Daily Cigarettes}

An astonishing evaluation result is the significant decrease in the average number of daily cigarettes as reported in Subsection \ref{subsec:amount-cigarettes}. Although the hypnosis intervention did in most cases \textit{not} lead to tobacco abstinence for the participants in the \textit{Lazy group}  and in the \textit{Placebo effect group}, we observe at least some ``partial success'' in a noteworthy decrease in the amount of daily cigarettes consumed. In average, the participants who still smoke after six months have reduced their tobacco consumption by $31\%$.
We remark that it is rather unclear whether the  decrease in daily cigarettes in the \textit{Lazy group} arises as some partial hypnotic effect or as ``solidarity'' like in the \textit{Placebo effect group}.
We compare and discuss the groupwise decrease in daily cigarettes in the upcoming sections in more detail.

\subsection{\textit{Motivated Group} vs. \textit{Lazy Group}}
\label{subsec:discussion-M-L}

Let us discuss the different evolution in both groups. First of all, the success rates in Figure \ref{Nonsmoker-lazy-main} demonstrate a significant different behaviour: while the number of non-smokers in the \textit{Motivated group}  is very high and increasing or constant on a high level, the number of non-smokers in the \textit{Lazy group} is monotonously decreasing  already after Session $2$. The boxplots in Figure \ref{cigarettes-decrease-prozentual-boxplot} suggest that the  \textit{Lazy group} has a behaviour similar rather to the \textit{Motivated group} after six weeks but becoming similar to the \textit{Placebo effect group}
after six months. Our interpretation is that the hypnosis sessions might have had some effects, but these effects were not sustainable for various reasons like an insufficient collaboration regarding the intervention protocol. In particular, we want to mention that several people in the \textit{Lazy group} consumed cigarettes within three hours after the first/second hypnosis sessions. Once again, this shall be regarded as an inadequate collaboration. In particular, those subjects did not listen to the audio files often, and therefore are assigned to the \textit{Lazy group}.
\par
The effect of the intensity of the individual post-session treatment (i.e., the number of audio listenings) is demonstrated in Figure \ref{fig:success-listenings} which suggests that participants shall be well-prepared in advance that they have to contribute to the success after the hypnosis sessions.
\par
Concerning the self-confidence and strength of wish in becoming non-smoker, the  \textit{Motivated group} had more self-confidence, but the \textit{Lazy group} had in average a higher strength of wish, see Subsection \ref{subsec:self-confidence}. The higher strength of wish in the \textit{Lazy group}  could possibly serve as an indicator for too high expectations. The lower self-confidence after Session $1$ in the \textit{Lazy group} can be interpreted as a hint that people have lost some trust in the success of  the intervention process, which could be the starting point for the lack of motivation afterwards.
\par
Our conclusion is that people shall be motivated after Session $1$ to follow the intervention protocol and not to expect miracles; instead, they shall be aware of an ongoing intervention process which affords collaboration.
\par
As a summary, if a subject attended the hypnosis sessions in a serious manner and followed the intervention protocol as requested, then there is  a $86\%$ chance that this subject is \textit{abstinent of tobacco consumption after six months}.

\subsection{\textit{Lazy Group} vs. \textit{Placebo Effect Group}}

We compare the evolution of both groups, where the \textit{Placebo effect group} serves as a control group.  When considering non-smoker rates in Figure \ref{Nonsmoker-lazy-main}  the rates of the \textit{Placebo effect group}   are less or equal than the corresponding ones in the  \textit{Lazy group} at each evaluation stages, and in the first two evaluation stages the numbers differ even significantly.
However, the gap of the non-smoker rates  closes towards the end of the evaluation period. In Section \ref{subsub:compare} we have verified via statistical tests which affirm the conjecture
that there is some significant different evolution in  both   groups at the beginning of the intervention process, which disappears after six months. A similar observation is obtained when considering the estimated densities of the reduction of daily cigarettes in Figure \ref{fig:densities}:  while the density of the \textit{Lazy group} has a heavy tail on the right hand side after six weeks, the density after six months looks like rather similar to the corresponding density of the \textit{Placebo effect group}. Figure \ref{cigarettes-decrease-prozentual-boxplot} is also a hint that the \textit{Lazy group} has a different behaviour during the first half of the intervention process compared to the second half.
\par
We remark that there are also some other general  aspects which are different in both groups: the \textit{Placebo effect group} has the biggest tobacco consumption among all three groups and the \textit{Lazy group} has the lowest one. Moreover, as pointed out in Subsection \ref{subsec:lazy} fall-back rates in the \textit{Lazy group} seem to be strongly dependent on specific situations, where the participants are \textit{not} strong enough in resisting to smoke. However, it seems to be quite hard or even impossible (due to insufficient data) to take these observations into account for our evaluation in a reasonable way.
\par
Overall,  from the observations made above we come to the conclusion that there were some presumably hypnotic effects in the \textit{Lazy group} during the first two stages which were \textit{not} sustainable.
The subjects in the \textit{Lazy group} probably lowered their chances to become smoke-free by not following our intervention protocol guidelines. This emphasises once again the necessity of post-treatments.

\subsection{General Discussion}

Let us discuss the general setting of online hypnosis itself. When people visit a hypnotherapist in his office rooms, people can become nervous which maybe omits a sufficient deep state of trance. Vice versa, during an online hypnosis session people can stay at their comfort zone at home and can feel more safety due to the physical distance allowing better relaxation. However, according to our experience this distance can lead to much less seriousness in attending the hypnosis sessions: people may move around and stay not concentrated. A behaviour like this seems to be much less probable when people visit a hypnotherapist in person and pay a multiple for his service. This seems to lead to a rather higher drop-out rate for online hypnosis sessions compared to sessions in person. We have observed these effects with people in the \textit{Placebo effect group}.
\par
We underline that this study is only a first evaluation of the efficiency of online hypnosis for smoking cessation. In particular, we didn't make any comparison to alternative methods for smoking cessation. Furthermore, it was not possible to adapt the post-hypnotic suggestions for each single participant, which could possibly reduce the fall back rates after Session $1$ (e.g., some participants started to smoke the morning after). Therefore, it is strongly recommended to adapt the post-hypnotic suggestions individually, since every smoker has specific preferred situation in which he smokes. Post-hypnotic suggestions targeting these situations will very likely have an even better impact on the success probability. It remains open at this point how specific suggestions will improve the success rates.
Further open questions consider possible existence of \textit{a-priori} measureable, individual factors which influence the success rate.
\par
As a summary, we have seen that online hypnosis works well for smoking cessation under some premises. People have to be aware how hypnosis works, they shall participate in an appropriate manner (which means, they shall allow themselves to relax) and, in particular, they shall understand that hypnosis still needs some own motivation and some own effort (i.e., individual post-session treatment) for being successful. In other words, we may conclude: if a subject shows this discipline then the failure rate is quite low. At this point we clearly recommend that subjects are getting well-prepared before the start of the intervention and  subjects are motivated even after Session $1$ to follow the intervention protocol.

\begin{comment}
\newpage
\section{old}

\begin{figure}[H]
\includegraphics[scale=.4]{Nonsmoker-rate-overall.jpeg}
\caption{Evolution of non-smoker rate}
\label{Nonsmoker-overall}
\end{figure}

\begin{figure}[H]
\includegraphics[scale=.4]{Nonsmoker-rate-gender.jpeg}
\caption{Gender-dependent evolution of non-smoker rate}
\label{Nonsmoker-gender}
\end{figure}

While in average the participants smoked $19$ cigarettes daily, after $6$ months, the average of daily cigarettes reduced to $9.3$ cigarettes per day, that is, in average the participantes smoke $51.1\%$ less after $6$ months compared to the time before starting the intervention. Moreover, $82.6\%$ of the participants smoke less after $6$ months than before the first hypnosis session.

\begin{figure}[H]
\includegraphics[scale=.4]{RelativeTageszigarettenVergleich.jpeg}
\caption{Relative decrease in daily cigarettes after 6 months compared to before the intervention process}
\label{fig:RelativeTageszigarettenVergeich}
\end{figure}

\begin{figure}[H]
\includegraphics[scale=.4]{TageszigarettenVergleich.jpeg}
\caption{Absolute decrease in daily cigarettes after 6 months compared to before the intervention process}
\label{fig:TageszigarettenVergeich}
\end{figure}

Another interesting fact is the time of the participants' first daily cigarette after standing up in the morning: almost half of the participants smoke their first cigarette almost immediately after waking-up, and $74\%$ smoke within half an hour after waking-up, see Figure \ref{fig:first-cigarette}.
\begin{figure}[H]
\includegraphics[scale=.29]{ErsteTageszigarette.jpeg}
\caption{Time of the participants' daily first cigarette after standing up in the morning.}
\label{fig:first-cigarette}
\end{figure}

\subsection{Still open questions}

\begin{itemize}
\item Zusammenhang Erfolg mit Wunsch, Sicherhheit aufzuh\"oren, Hypnoseerfahrung, fr\"uhere Versuchen, verschiedene Gruppen (Stress, Psycho etc)

\item Erste Tageszigarette noch reinnehmen?
\item Anpassung der suggestionen an Gruppen

\end{itemize}

$$
\begin{array}{r|c|c|c|c}
&\textrm{Main group} & \textrm{Lazy group} & \textrm{Placebo effect group} & \textrm{Overall}\\
\textrm{Number of participants} & 28 & 51 & 20 & 99 \\
\textrm{Mean wish} & 8.19 & 8.55 & 8.25 & 8.38 \\
\textrm{Mean stress level} & 5.36 & 5.1 & 5.7 & 5.30 \\
\textrm{Mean subliminal listenings FB2} & 6 & 6.37 & 5.45 & 6.08 \\
\textrm{Mean subliminal listenings FB3} & 2.93 & 6.16 & 4.85 & 4.98 \\
\textrm{Mean audio listenings FB2} & 7.04 & 6.31 & 6.30 & 6.51\\
\textrm{Mean audio listenings FB3} & 6.04 & 6.33 & 5.05 & 5.99\\
\textrm{Mean years of smoking} & 22.9 & 25.0 & 27.4 & 24.9\\
\textrm{Mean age of starting smoking} & 15.1& 16.5 & 15.1&15.8\\
\textrm{Selfconfidence in quitting smoking} & 3.46 & 3.22 & 3.55 &
\end{array}
$$

\newpage

\section*{Appendix}
\label{sec:appendix}

\begin{figure}[H]
\includegraphics[scale=.4]{Alter.jpeg}
\caption{Age of participants}
\label{fig:alter}
\end{figure}

\begin{figure}[H]
\includegraphics[scale=.4]{ZigarettenProTag.jpeg}
\caption{Average number of daily cigarettes}
\label{fig:zigaretten-pro-tag}
\end{figure}

\begin{figure}[H]
\includegraphics[scale=.4]{Rauchbeginn.jpeg}
\caption{Age of starting smoking}
\label{fig:rauchbeginn}
\end{figure}

\begin{figure}[H]
\includegraphics[scale=.4]{Raucherjahre.jpeg}
\caption{Years of smoking}
\label{fig:raucherjahre}
\end{figure}

\end{comment}

\bibliographystyle{abbrv}


HIP 56948: A Solar Twin With a Low Lithium Abundance

Abstract: For more than a decade, 18 Sco (HD 146233) has been considered the star that most closely resembles the Sun, even though significant differences such as its Li content, which is about three times solar, exist. Using high resolution, high S/N spectra obtained at McDonald Observatory, we show that the stars HIP 56948 and HIP 73815 are very similar to the Sun in both stellar parameters and chemical composition, including a low Li abundance, which was previously thought to be peculiar in the Sun. HIP 56948, in particular, has stellar parameters identical to solar within the observational uncertainties, being thus the best solar twin known to date. HIP 56948 is also similar to the Sun in its lack of hot Jupiters. Considering the age of this star (1+/-1 Gyr older than the Sun) and its location and orbit around the Galaxy, if terrestrial planets exist around it, they may have had enough time to develop complex life, making it a prime target for SETI.

\section{Introduction}
Even though it is the star we know best, the Sun can not be used as the primary calibrator in stellar astrophysics because it is exceedingly bright for the instrumentation utilized for distant stars. In many cases, one would like to observe a calibrating star whose fundamental properties (temperature, luminosity, metallicity, mass and age) are identical to solar but can be analyzed in a similar way as the objects of interest.

The study of ``solar twins'' (Cayrel de Strobel 1996), stars that are spectroscopically and photometrically identical to the Sun, is necessary to set the zero point of fundamental calibrations in astrophysics, such as the color-temperature relations and the absolute flux scale (e.g., Colina et~al. 1996, Holmberg et~al. 2006). Further applications extend to many other fields. For example, in a study of the physical properties and classification of asteroids within the GAIA mission (e.g., Perryman 2005), the solar reflectance spectrum must be first removed from the asteroid observations to allow for their accurate physical modeling (A.~Korn, private communication). This can be achieved with faint solar twins, which will be found and understood more easily once the bright ones have been reliably characterized.

Solar twins can also be used to validate stellar interior and evolution models. For example, the
depletion in the observed solar Li abundance by about a factor of 200 relative to that found in
meteorites is not explained by standard stellar evolution models. Since more physically
realistic models, calculated from basic principles, cannot be computed yet, non-standard stellar
evolution models have to introduce free parameters to reproduce the observed Li
depletion (e.g., D'Antona \& Mazzitelli 1984; Ventura et~al. 1998; Charbonnel \& Talon 2005).
Solar twins and analogs with well-determined physical and chemical properties can help validate models by providing a sample covering a wide range of Li depletion. The Sun appears to be lithium-poor by a factor of 10 compared to similar solar type disk stars (Lambert \& Reddy 2004), a fact that has led to the suggestion that the Sun is peculiar in Li and therefore of dubious value for calibrating non-standard models of Li depletion.

Only three solar twins are currently known: 18~Sco (Porto de Mello \& da Silva 1997), HD~98618 (Mel\'endez et~al. 2006)
and HIP~100963 (Takeda et~al. 2007). However, these three stars have a Li abundance a factor
of 3 to 6 higher than solar. Furthermore, according to our criteria, these stars are solar twins
only marginally and should be called ``quasi solar twins'' instead.\footnote{Our definition of solar twin requires
the star to have fundamental physical parameters (e.g. T$_{\rm eff}$, log {\it g}, [Fe/H], $v_t$) identical to solar within
a $1\sigma$ radius of observational and systematic errors. As we will see in $\S$2, HIP 56948 is the only star that can strictly
be called solar twin. Indeed, the average of six high-precision differential analyses
(Table~3 of Mel\'endez et al. 2006, Takeda et al. 2007, our Table~1)
shows that 18 Sco, the previous best solar-twin candidate, has [Fe/H] = +0.035 ($\sigma$ = 0.010) dex,
i.e., significantly different from solar.}

Here we report the discovery of two solar twins with a low Li abundance. This study is part of a survey of solar twins and analogs in the northern hemisphere that is being completed at McDonald Observatory and will be described thoroughly in a forthcoming publication (Ram\'{\i}rez et~al. 2008).

\section{Data and Analysis}

We selected candidate solar twins from a sample of more than 100,000 stars listed in the \textit{Hipparcos} catalog by employing precise blue-optical-red-infrared color:T$_{\rm eff}$  relations (Ram\'{\i}rez \& Mel\'endez 2005) with a zero-point correction based on our pilot solar-twin survey (Mel\'endez et~al. 2006), trigonometric parallaxes, and, when available, age indicators.

High resolution ($R\sim60,000$), high signal-to-noise ratio ($S/N\sim200-700$) spectra were obtained for 23 solar twin candidates and the asteroid Vesta (reflected
solar spectrum) with the 2dcoud\'e spectrograph (Tull et~al. 1995) on the 2.7-m Harlan J. Smith Telescope at McDonald Observatory in April 2007. Standard data reduction
procedures were applied to the spectra using IRAF. Our spectra cover a wide, albeit incomplete, spectral range ($380-1000$~nm), and includes iron lines in a broad range of excitation potentials ($0-5$~eV) and line strengths. Sample spectra of the best solar twins are shown in Fig.~1.

The first analysis of the data was model-independent. Following Mel\'endez et~al. (2006), we measured the relative difference in equivalent width, $\delta W_r =
(W_\lambda^* - W_\lambda^\odot)/W_\lambda^\odot$, for a large set of \ion{Fe}{1} lines, using their median value $<\delta W_r>$ as a proxy for metallicity, while the slope in the $\delta W_r$ vs. excitation potential relation was used as a temperature indicator. Stars for which these values were zero within the observational uncertainties were considered for further detailed analysis. The results were corroborated using line-depths instead of equivalent widths. In this way we found that HIP~56948 is spectroscopically indistinguishable from the Sun and significantly a better solar twin than the other two best previously known candidates (18~Sco and HD~98618). Another close twin was identified, HIP~73815, which is, however, not as good as HIP~56948.

A detailed model-atmosphere analysis of the best candidates was performed to verify the outcome of the empirical solar-twin test and estimate precise stellar parameters, as in Mel\'endez et~al. (2006). Kurucz model atmospheres and mostly laboratory $gf$-values were employed (see Mel\'endez et al. 2006). Additionally, we measured the chromospheric emission $S$ index (Wright et al. 2004) and relative $v\sin i$ values from our spectra. The model-atmosphere analysis confirmed the empirical results, reinforcing the solar twin status of 18~Sco and HD~98618, and adding two more stars, HIP~56918 and HIP~73815, to the list of known solar twins. The parameters derived for these four solar twins are given in Table~1. HIP~56948 has solar parameters within 1-$\sigma$, while the other three candidates are solar twins within 2-$\sigma$. Relative chemical compositions are shown in Fig.~2. Except for the case of lithium (see \S3), the solar twins studied here have solar abundances within the uncertainties.

The ages of the solar twins discussed here were determined using Y$^2$ isochrones (Demarque et~al. 2004) and the chromospheric emission index $S$ (e.g., Wright et~al. 2004). Both HIP~56948 and HIP~73815 are older than the Sun by 1.2 and 1.7 Gyrs, respectively. On the other hand, 18~Sco is about 1.2~Gyrs younger than the Sun, while HD~98618 has an age similar to solar (see Table 1).

A radial velocity monitoring of what we consider to be the best solar twin ever known, HIP~56948, has been started by the McDonald Observatory planet search program (e.g., Endl et~al. 2005). Preliminary results reveal that no hot Jupiter is present around the star. Giant planets around the other two solar twins (18~Sco and HD~98618) have not been reported either by the California and Carnegie Planet Search Project (Marcy et~al. 2005). Details on the radial velocity follow up will be given in a forthcoming paper (Ram\'{\i}rez et~al. 2008).

We have also determined the Galactic space velocities and orbital parameters of the four solar twins. Like the Sun, they are in low eccentricity orbits ($e<0.1$) close to the galactic plane at a galactocentric distance of about 8\,kpc and, therefore, have an almost 100\% probability of being members of the thin disk. Furthermore, the maximum distance from the galactic plane for the orbit of HIP~56948 is only about 50\,pc.

\section{The solar twins HIP 56948 and HIP 73815}

For many years, the Sun has been thought to be peculiar in its low lithium content. This work, as well as a recent study by Takeda et al. (2007), show that this is not the case. We have found that both HIP~56948 and HIP~73815 have a low lithium abundance, similar to solar (Table 1 and Figs. 1--2), while Takeda et al. have found that some stars closely resembling the Sun also have a low Li abundance. Curiously, the other two quasi-twins (HD~98168, 18~Sco) that we have analyzed have a Li abundance about three times higher.
HIP~100963, the other quasi-solar twin recently identified by Takeda et al. (2007), shows an even much higher Li abundance (about 6 times solar).

In order to understand the large spread in Li abundances among the five known solar twins, more high resolution, high S/N spectroscopic studies are needed to improve the precision of the
Li abundances of the Li-poor stars and explore other fragile elements such as beryllium.
As noted by several authors (e.g., D'Antona \& Mazzitelli 1984;
Deliyannis \& Pinsonneault 1997; Ventura et~al. 1998; Charbonnel \& Talon 2005), the amount of Li depletion
depends on several ingredients such as mass, age, convection treatment,
extra-mixing mechanisms (e.g., rotation), mass loss, metallicity, and magnetic fields.
Since Li and Be are destroyed at different temperatures
($\approx$ 2.5 and 3.5 $\times 10^6$ K; Deliyannis \& Pinsonneault 1997), their observed abundances can be used to reveal the dominant mechanism(s) causing only some stars to have such low Li abundances.

Two apparent trends in our small solar twin sample may help understand the large Li spread.
One is that the twins with the highest masses
(18 Sco and HD 98618, $\approx 3-4$\,\% more massive than solar) are less depleted in Li,
in qualitative, albeit not quantitative (apparently too large Li destruction for the small mass variation), agreement with the models of Li depletion cited above. Note, however, that the solar analogs
16 Cyg A and 16 Cyg B have masses that differ by only about 2 \% (as estimated
by our method and the relative stellar parameters determined by Deliyannis et al. 2000), yet they show very large differences (about a factor of 5) in their Li abundances (Deliyannis et al. 2000), which suggests a depletion that depends exponentially on mass. The other apparent trend is that the Li abundances in the solar twins and the Sun seem to be correlated with age, with a decrease
of about 0.1 dex in Li abundance for an increase of 0.35 Gyrs in age, in good
agreement with the predicted depletion of about 0.1 dex in Li per 0.3 Gyrs for a solar model including internal gravity waves (Fig.~2 in Charbonnel \& Talon 2005). This age correlation seems to apply
also to the 16~Cyg system, which according to our methods should be about 2--3~Gyr older than the Sun, i.e., they should have low Li abundances, and they indeed show Li abundances similar to solar (Deliyannis et~al. 2000) although 16~Cyg~B has a Li abundance about a factor of 5 lower than 16~Cyg~A.
Perhaps both (age and mass) correlations are correct, explaining thus
why 16~Cyg~A and B have low Li abundances (age effect) and at the same
time why 16~Cyg~B has much lower Li abundance than 16~Cyg~A (mass effect).
We stress here the preliminary nature of our results. Many more solar twins must be identified to confirm or reject the apparent trends with mass and/or age and help constrain stellar models of Li depletion reliably.

\section{Conclusions}

We have found two solar-twins (HIP~56948 and HIP~73815) that have a low Li abundance, similar to solar. This is remarkable given that this is a Li abundance considerably lower (by a factor of $3-6$) than that found in three previously identified solar twins (18~Sco, HD~98618, and HIP~100963). The determination of precise abundances of other fragile elements (e.g., Be) in these solar twins may lead to an understanding of the still unknown mechanism responsible for the large Li depletion in the Sun.

Neither HIP~56948, 18~Sco, nor HD~98618 have hot Jupiters. Long-term radial velocity
monitoring of these objects is encouraged, given that it could reveal, in about a
decade, whether gaseous giant planets reside within the planetary systems that these
stars may host.

HIP~56948 has physical properties and a chemical composition identical to that of the Sun, as well as a  similar age. Stellar ages are key in the search for complex life in the universe, especially if one considers that the formation timescale of complex life is similar to that for the Earth. SETI searches can be optimized by selecting stars with ages similar to solar. Furthermore, given its orbital properties, HIP 56948 is well within the Galactic Habitable Zone (e.g., Lineweaver et~al. 2004). Considering all these facts, HIP~56948 is an excellent candidate for life outside the solar system.

The two new solar twins presented here have not been mentioned specifically as prime candidates for habitable systems in the literature, probably because ours is the first detailed spectroscopic analysis of them. Although the twin candidates are included in a catalog of nearby habitable systems (Turnbull \& Tarter 2003) based on the \textit{Hipparcos} catalog, so are thousands of other FGKM stars with parallaxes. We suggest SETI programs to give our twins top priority in their quest for intelligent life in the universe.

Altogether, the four solar twins presented in Table~1 cover the immediate past, present and future of the Sun. Long-term studies of these stars can help us understand variations of the solar activity cycle,
which are very difficult to predict theoretically. Sunspot studies cover only a few
centuries, while indirect measurements of solar activity extend to only about ten millennia
(e.g., Usoskin et~al. 2007). Solar twins offer the opportunity to study solar activity over much longer timescales. In fact, chromospheric emission measurements of 18 Sco reveal an activity cycle of about seven years (Hall et~al. 2007),  significantly shorter than the 11-year solar cycle. In addition, these studies may help to understand anomalous variations such as the extreme low activity of the Maunder Minimum (Maunder 1890, Eddy 1976), which can be particularly interesting if they are related to climate changes on Earth (Eddy 1976, Shindell et~al. 2001). Apparently, we know no Maunder minimum star yet (Wright 2004).

We have started photometric measurements ($uvby$ and $UBVRI$) of our solar twin sample.
The photometry of the solar twins will be employed to derive precise solar colors (the average of the colors of solar twins) and determine accurately the zero-points of the color-T$_{\rm eff}$  relations, for example those of Ram\'{\i}rez \& Mel\'endez (2005), which are reliable in a relative sense: temperature differences between stars are known with good precision but the actual T$_{\rm eff}$  values may be inaccurate due to potential systematic errors in the absolute flux calibration. For accurate estimates of T$_{\rm eff}$  values, a zero point correction is required. Since the temperature of the Sun is known with extreme accuracy, a very well determined solar color (e.g., four solar twins can provide a solar $B-V$ color with an uncertainty of about 0.005~mag) will result in a zero point for the T$_{\rm eff}$  scale as accurate as 25~K. We encourage other groups to observe these stars and extend these fundamental calibrations to as many photometric systems as possible.

Besides the astrophysical significance of solar twins, there is another motivation to look for and study them in full detail: they help answering the question of whether the Sun is unique or not (Gustafsson 1998, Gonzalez 1999), a question that has important philosophical consequences. An anomalous Sun favors some forms of the anthropic principle, those that suggest that our location in the universe is privileged for our existence as observers, but the identification of HIP~56948 as a real solar twin, although not disproving it completely, can be used as an argument against it. Also, due to its potential for hosting habitable planets, HIP~56948 may support a Copernican view for habitability in the universe.

\acknowledgements

IR would like to thank the support of the Robert~A.~Welch Foundation of Houston, Texas to David~L.~Lambert. JM acknowledges support from the Australian Research Council to Martin Asplund. We thank W.~Cochran and M.~Endl for providing the first results on the search for planets around HIP~56948 and the anonymous referee for comments that helped improve the manuscript.


Quasars as Probes of Late Reionization and Early Structure Formation

Abstract: Observations of QSOs at z ~ 5.7 - 6.4 show the appearance of Gunn-Peterson troughs around z ~ 6, and a change in the slope of the IGM optical depth tau(z) near z ~ 5.5. These results are interpreted as a signature of the end of the reionization era, which probably started at considerably higher redshifts. However, there also appears to be a substantial cosmic variance in the transmission of the IGM, both along some lines of sight, and among different lines of sight, in this intriguing redshift regime. We suggest that this is indicative of a spatially uneven reionization, possibly caused by the bias-driven primordial clustering of the reionization sources. There is also some independent evidence for a strong clustering of QSOs at z ~ 4 - 5 and galaxies around them, supporting the idea of the strong biasing of the first luminous sources at these redshifts. Larger samples of high-z QSOs are needed in order to provide improved, statistically significant constraints for the models of these phenomena. We expect that the Palomar-Quest (PQ) survey will soon provide a new set of QSOs to be used as cosmological probes in this redshift regime.

\section{Introduction}

Understanding of the cosmic reionization era -- from the appearance of the
first luminous sources perhaps at $z \sim 20$, to the end of the IGM phase
transition from neutral to essentially fully ionized hydrogen at $z \sim 6$
-- is now perhaps the focal arena of cosmology.  This field connects a
number of fundamental astrophysical processes, formation of the first stars,
galaxies, AGN, and large-scale structure, and their effect on the early IGM,
more than any other redshift range studied so far.  For recent
reviews see, e.g., \cite{loe+01, bar+01, mir+03, djo+05}.

There are so far three kinds of observational constraints on the nature and
extent of the reionization era.  Historically the first was the evidence of
the approach to the reionization, or its final stages, at $z \sim 6$ from
the spectra of QSOs in this redshift range
\cite{bec+01, fan+01, fan+03, djo+01, whi+03}.

The primary evidence is the existence of Gunn-Peterson (hereafter GP)
\cite{gun+65} absorption troughs in the spectra of $z > 6$ QSOs.
Such features are now seen along all lines of sight at $z > 6$ where
adequate data exist.
There is also a growing evidence for a change in the slope of the optical
depth vs. redshift, $\tau(z)$~\cite{fan+02,cen+02,whi+03}
around $z \sim 5.5 - 6$, indicating that some qualitative change in the
physics and geometry of IGM occurs at these redshifts.
Transition from the thickening of the normal Ly$\alpha$ forest out to
$z \sim 5$ to the UV-opaque IGM at $z \sim 6$ is fairly dramatic (see Fig. 1).
As discoveries of more such objects continue \cite{fan+04},
spectroscopy of high-$z$ QSOs remains one of the principal empirical
approaches to mapping of the final stages of reionization.

\begin{figure}[!t]
\centerline{\psfig{file=dba_fig1.eps,width=5.0in}}
\caption{
From the forest to the fog: continuum-normalized QSO absorption spectra at
progressively higher redshift windows.  The forest thickens dramatically at
$z \sim 5.5$ and GP-type troughs are ubiquitous at $z \ge 6$.
Spectra of several high-$z$
SDSS QSOs, obtained at the Keck-II telescope with the ESI instrument were used.
From \cite{djo+02,djo+05}.
}
\label{}
\end{figure}

The picture was complicated by the evidence for an $early$ reionization at
$z \sim 10 - 20$ from the CMBR measurements by the WMAP satellite~\cite{kog+03}.
While the confirmation of this important result is still to come, we note
that it is not necessarily contradictory to the QSO results: the history of
the luminous sources and their effect on the IGM was probably highly complex,
and there was a finite time interval from the appearance of the first sources
of UV photons at the end of the cosmic dark ages, and the completion of the
reionization phase transition a few hundred million years later.

In parallel, a number of groups have studied star-forming galaxies at
$z \sim 6 - 7$, and measurements of the Ly$\alpha$ emission line luminosity
function evolution provide another useful observational constraint
\cite{mal+04,ste+05}.
While the QSO absorption
spectra probe the neutral hydrogen fraction regime $x_{H I} \leq 10^{-2}$,
this method is sensitive to the range $x_{H I} \sim 10^{-1} - 10^{0}$.

In the future, using H I hyperfine transition line at 21 cm promises to
be another powerful probe of the earlier stages of reionization, as
reviewed in a number of talks at this conference.  Detection and spectroscopy
of GRB afterglows from the first generations of massive stars will be another
spectacular method to study the astrophysics of reionization.

If, as both the data and theory now suggest, the distribution of the first
luminous sources was highly clustered, leading to a spatially uneven
reionization phase transition, direct measurements of the primordial
clustering would provide a complementary insight into the overall picture.
Luminous QSOs at these redshifts are potentially a powerful probe of such
bias-driven clustering; see, e.g.,~\cite{djo+05} and references therein.

Here we review some of the current evidence for a spatially inhomogeneous
reionization, and the possible uses of high-$z$ QSOs to study the late
stages of this fundamental cosmological era.  In some ways, QSOs provide
unique observational constraints, and are in many ways complementary to
other approaches.  We also outline some possibilities for the future work.

\section{Evidence for a Clumpy Reionization?}

The reionization era must be extended in redshift, from the time the first
luminous sources turn on, to the final percolation of the reionized bubbles
universe-wide.  At any given intermediate redshift, there will be some
spatially uneven distribution of the neutral and ionized IGM, whose
characteristic physical scales and topology will depend on the clustering
and luminosity histories of the local reionization sources.  In principle,
such variations in the geometry/topology of reionization should be
detectable as variations in the IGM transmission, both along and among
different lines of sight to sources still within the incomplete reionization
regime.

Such variations are indeed observed, and can in principle be used to place
quantitative constraints on theoretical models~\cite{djo+02,djo+05,bog+06};
see Figs. 2 and 3.  Measurements of the Ly$\alpha$ and Ly$\beta$ optical
depth to a number of SDSS high-$z$ QSOs suggest a comparable cosmic
variance in the IGM transmission in this redshift regime \cite{str+05}.

\begin{figure}[!t]
\centerline{\psfig{file=dba_fig2.eps,width=3.5in}}
\caption{
A substantial variety in the IGM absorption along 4 different QSO lines of
sight is apparent in these continuum-normalized spectra covering the same
redshift window, which just follows the reionization epoch.  This variance
may be a residual effect of an uneven end to reionization,
as expected from the biased structure formation picture.
From \cite{djo+02,djo+05}.
}
\label{}
\end{figure}

What is interesting here is that there may be an apparent increase in the
cosmic variance of the IGM transmission as one approaches the ostensible
end-of-the-reionization redshift, $z \sim 6$.  The variance is impossible
to measure as the flux disappears at higher redshifts, but of course a
substantial variation in the neutral hydrogen fraction could be present
in the optically thick regime; this calls for some alternative measurement
methodology.

\begin{figure}[!t]
\centerline{\psfig{file=dba_fig3.eps,width=5.0in}}
\caption{
The behavior of the transparency of the IGM at high redshifts, shown as the
transmission on the left, and as the optical depth on the right.  The bottom
panels indicate the number of lines of sight used at every redshift.
The solid lines indicate the mean, and the hashed areas the spread among the
different lines of sight.  The data have been smoothed with a Gaussian with
a $\sigma(z) = 0.02$; a wider smoothing or binning, preferred by some authors,
can wash out important features.
Note the change in the slope of $\tau(z)$ around $z \approx 5.8$; this
corresponds to the dramatic change in the ionizing flux and the fraction of
the neutral hydrogen at these redshifts \cite{fan+02}, which
is interpreted as the signature of the end of the reionization era.
There is also a hint that the spread among different lines of sight (i.e.,
the cosmic variance in the IGM absorption properties) increases just before
the reionization redshift is reached.
From~\cite{djo+02,djo+05}.
}
\label{}
\end{figure}

Interestingly, high quality QSO spectra at $z > 6$ can be used to
measure the sizes of the H II regions around these objects, and then to
constrain the neutral hydrogen fractions in the regime
$x_{H I} \sim 10^{-1} - 10^{0}$~\cite{wl04a,mes+04}
 From a single QSO, Wyithe \& Loeb conclude that $x_{H I} > 10^{-1}$
in its vicinity.  However, this is the $x_{H I}$ range also probed by the
census of the Ly$\alpha$ emitters~\cite{mal+04,ste+05},
where the preliminary results indicate $\langle x_{H I} \rangle < 10^{-1}$
in the general field.  There are both small number statistics and
interpretative uncertainties in both
methods, and clearly at least a few more QSO Str\"omgren sphere measurements
are needed.  If the two approaches continue to give systematically different
results as the data accumulate, that may be
indicative of systematic differences in the large-scale environments, e.g.,
with QSOs being situated in particularly dense regions.

A strong spatial nonuniformity in the final stages of reionization could
be caused by a strong clustering of the reionization sources.  Since there
was not enough time for a substantial large scale structure growth at these
redshifts, the clustering would have to be primordial, e.g., reflecting a
generic biased galaxy formation scenario, where the first luminous sources
form in and around the highest peaks of the primordial density field,
which are a priori strongly clustered almost regardless of any details
of the assumed structure formation models~\cite{kai+84}.  Thus, we expect
that there was a fundamental connection between the early structure formation
and the geometry and topology of the reionization, as reflected in the cosmic
variance of the IGM transmission at these redshifts.

\section{QSOs as Markers of the Primordial Large-Scale Structure}

There is now a substantial and growing evidence for a strong bias among
the field galaxies at redshifts approaching 6; see, e.g.,~\cite{djo+05}
for a review and references.  This is entirely as expected from models of
biased galaxy formation.

Many or all luminous QSOs at these redshifts are also likely situated in
massive host galaxy halos.  This is suggested by the high masses of their
central black holes ($M_\bullet > 10^9 M_\odot$);
see, e.g.,~\cite{wil+03,die+04,mcl+04};
as well as their highly evolved chemical abundances
\cite{ham+93,ham+99,mat+93,die+03a,die+03b}.
Direct observations of circum-QSO regions at
$z > 4$ show a considerably enhanced density of star forming galaxies in
their vicinity~\cite{djo+98,djo+99a,djo+99b,sti+05}
consistent with the idea that luminous high-$z$ QSOs may mark sites of
future rich clusters of galaxies.

QSOs themselves appear to be strongly clustered at these redshifts;
see~\cite{djo+03a,djo+05} and references therein.  The implied
clustering lengths are comparable to those of the rich clusters today,
tens of Mpc.
An even more intriguing possibility is that QSO clustering at $z > 4$
extends to physical scales of $\sim 100 - 200 h^{-1}$ comoving Mpc
\cite{djo+98,djo+99b}.
These highly preliminary, uncertain, but
suggestive results must be confirmed by careful analysis of statistical
samples of high-$z$ QSOs.

If QSO clustering at $z \sim 4 - 6$ is real,
it has to be bias-driven.  It would be a nearly unique probe of the
clustering of luminous sources at these redshifts, on physical scales
$\sim 10 - 100$ comoving Mpc, and directly relevant for the interpretation
of the cosmic variance of IGM transmission at the end of the reionization.

Also, if substantial spatial variations in the distribution of the ionized gas
existed at these redshifts, they might have also left an observable imprint
in the CMBR fluctuations, essentially the clustered S-Z effect.  In this
context, it is intriguing to recall the excess power at high angular
frequencies, corresponding to the comoving separations of $\sim 10 - 20$ Mpc
at the CMBR photosphere; see \cite{rea+04} and references therein.

\section{The Palomar-Quest Survey}

As the preceeding discussion indicates, substantially larger samples of
high-$z$ QSOs in the late- and post-recombination redshift regime are
needed in order to provide better observational constraints on the
physical parameters, geometry, and evolution of the IGM in this crucial
era, as well as to quantify better the early, possibly highly biased
large-scale structure.

The Palomar-Quest (PQ) Digital Synoptic Sky Survey
\cite{pqwebsite, bau+03,rab+03,mah+03,djo+03b,djo+04}
will be a major new provider of high-$z$ QSOs at $z \sim 4 - 6.5$.
The PQ survey is a collaborative
effort between Yale, Caltech, Indiana University and NCSA in the USA,
with collaborations with INAOE in Mexico, and other groups.  The data are
taken at the 48-inch Samuel Oschin Schmidt telescope at Palomar
Observatory, using a special 112-CCD camera built for this purpose.
The data rate is about 1 TB per month.  Data are reduced and archived
at multiple locations (Caltech, NCSA, Yale).
National Virtual Observatory (NVO)
standards, protocols, and connections are built in from the start, and
will facilitate a very broad data access and analysis.

Approximately 50\% of the time is used for PQ survey drift scans, and
the rest is used largely in the traditional point-and-stare mode by
groups at JPL and Caltech for exploration of the Solar System.
In the drift scan mode, the PQ survey covers strips $\sim 4.6^\circ$
wide at a constant Dec, with an area coverage of $\sim 500\,deg^{2}$
per clear night, in 4 filters.  The pixel scale is 0.878 arcsec.  The
scans are obtained in the range $-25^{\circ} \leq \delta
\leq+25^{\circ}$, giving the total useful survey area of
$\sim 15,000\,deg^{2}$.  About a third of the survey area overlaps with
that of SDSS, thus providing a valuable cross-check and cross-calibration.
The survey uses two sets of
filters, the traditional $UBRI$, and the SDSS $r'i'z'z'$, with about
an equal time share in each.  Typical limiting magnitudes in a single
pass in good conditions are: $r' \sim 21.5$ mag, $i' \sim 20.5$ mag,
$z' \sim 19.5$ mag, $R \sim 22.0$ mag, and $I \sim 21$ mag.  Coadding
of about 6 passes reaches the depths comparable to SDSS.  As the
survey unfolds, we anticipate that we will cover much of the area
to at least a magnitude deeper, and a subset of the total survey
area deeper yet.

One of the key scientific goals of the PQ survey is the discovery
of a significant number of high-$z$ QSOs, to be used as a statistical
sample for studies of early structure formation and reionization.  The
methodology is identical to that used by SDSS, DPOSS, and other
multi-color high-$z$ QSO surveys to date, using colors as
discriminants for objects classified as point sources in the images.
High-$z$ QSO candidates
are identified in the PQ data using $BRIr'i'z'$ colors, which in
appropriate combinations can be used to find QSOs at $z \sim 4 - 6.5$.
IR photometry is necessary to separate QSO candidates at $z > 5.5$
from late-M and brown dwarfs.  We use coadds of 2MASS images to eliminate
some of the brighter dwarf contaminants, but most candidates
require a deeper IR imaging.  This is currently obtained at the Yale
SMARTS telescopes at CTIO, and at the INAOE/OAGH 2.1-m telescope at
Cananea, Mexico.  The surviving color-selected candidates are then followed
spectroscopically at the Palomar 200-inch telescope, and QSOs at $z > 5$
are selected for deeper and more extensive follow-up studies at the Keck
and other facilities.

The first high-$z$ QSOs have been found, but the survey is still in an early
stage.  We anticipate that we will double the sample of the SDSS high-$z$
QSOs over the next few years.  Studies of the joint SDSS+PQ sample will be
used to provide new observational constraints for structure formation at
$z \sim 4 - 6.5$, and late stages of reionization at $z \sim 6$.

\section{Concluding Comments}

As the reionization becomes complete, the residual H I fraction rapidly drops
over large spatial scales.  QSO absorption spectra are the $only$ currently
known probe of this IGM regime, $\langle x_{H I} \rangle \leq 10^{-2}$.
If this happens at $z \sim 6$, as the data seem to indicate, then
spectroscopy of QSOs in this redshift regime will remain a unique probe
of the late stages of reionization for some time to come.

Moreover, the data indicate that a considerable cosmic variance is present
in the transmission properties of the IGM at these redshifts, along
different lines of sight.  If the characteristic physical scales of the
late reionization topology are in the range of tens to hundreds of comoving
Mpc, as the recent theoretical studies suggest \cite{bl04, wl04b, of05},
then statistically significant samples of QSOs will provide the first solid
observational constraints for theoretical models.

If there is indeed a substantial, bias-driven clustering of the luminous
sources responsible for the reionization (which would then naturally be
spatially clumpy and lead to the observed cosmic variance in the IGM
transmission properties), QSOs may again prove to be useful probes of such
clustering, especially if they are associated with some of the densest
peaks of the density field at the time.  Until galaxy redshift surveys at
$z \sim 6$ start spanning comparable comoving volumes and $>> 100$
comoving Mpc scales -- which may be a while in coming -- QSOs may be
again the only viable probe of such primordial clustering and biasing.

Therefore, we expect that studies of larger high-$z$ samples than currently
available, coming from SDSS, PQ, and other surveys, will continue to
generate ever more significant observational constraints for the models,
leading to a better physical understanding of the final stages of the
reionization era.

\noindent {\bf Acknowledgements.}~~
We are indebted to numerous collaborators, and especially members of the
PQ survey team, for sharing the work and the excitement of a number of
research projects, past, current, and future, on which this paper is
based.  We also thank the staff of Palomar and Keck observatories for
their expert help during numerous observing runs.  This work was
supported in part by the NSF grant AST-0407448, and by the Ajax Foundation.
Processing of the PQ data is supported in part by the NSF grant AST-0326524,
and by resources at NCSA.