\section{#1}\setcounter{equation}{0}}
\renewcommand{\theequation}{\thesection.\arabic{equation}}

\newcommand{\ket}[1]{|#1\rangle}

\begin{document}

\baselineskip=17pt

\begin{titlepage}
\rightline{\tt arXiv:1207.3335}
\rightline{\tt  UT-Komaba/12-6}
\begin{center}
\vskip 0.7cm
{\Large \bf {Gauge-invariant observables
and marginal deformations}}\\
\vskip 0.4cm {\Large \bf {in open string field theory}} \vskip 1.0cm {\large {Mat\v{e}j
Kudrna,${}^{1}$ Toru Masuda,${}^2$ Yuji Okawa,${}^2$}} \vskip 0.3cm {\large {Martin Schnabl${}^1$
and Kenichiro Yoshida${}^3$}} \vskip 0.8cm
${}^1${\it {Institute of Physics AS CR}}\\
{\it {Na Slovance 2, Prague 8, Czech Republic}}\\
\vskip 0.3cm
${}^2${\it {Institute of Physics, The University of Tokyo}}\\
{\it {3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan}}\\
\vskip 0.3cm
${}^3${\it {Ministry of Internal Affairs and Communications}}\\
{\it {2-1-2 Kasumigaseki, Chiyoda-ku, Tokyo 100-8926, Japan}}\\
\vskip 0.3cm
matej.kudrna@email.cz, masudatoru@gmail.com, okawa@hep1.c.u-tokyo.ac.jp,\\
schnabl.martin@gmail.com, ken16.yoshida@gmail.com

\vskip 1.0cm

{\bf Abstract}
\end{center}

\noindent The level-truncation analysis of open string field theory for a class of periodic
marginal deformations indicates that a branch of solutions in Siegel gauge exists only for a finite
range of values of the marginal field. The periodicity in the deformation parameter is thus
obscure. We use the relation between gauge-invariant observables and the closed string tadpole on a
disk conjectured by Ellwood to construct a map between the deformation parameter of the boundary
conformal field theory and the parameter labeling classical solutions of open string field theory.
We evaluate the gauge-invariant observables for the numerical solutions in Siegel gauge up to level 12
and find that our results qualitatively agree with the analysis by Sen using the energy-momentum
tensor and are consistent with the picture that the finite range of the branch covers one
fundamental domain of the periodic moduli space.
\end{titlepage}

\newpage

\section{Introduction}
\setcounter{equation}{0}

String field theory can potentially be
a background-independent formulation of string theory.
The current construction of string field theory, however,
requires a choice of a consistent background,
and different backgrounds are expected to be described
by classical solutions of the theory
based on the original background we chose.

Thus the first step to address the problem of background independence in string field theory is the
construction of classical solutions. There has been remarkable progress in constructing analytic
solutions in open string field theory~\cite{Witten:1985cc} since the first construction of an
analytic solution for tachyon condensation by one of the authors~\cite{Schnabl:2005gv}. In
particular, analytic solutions for marginal deformations of boundary conformal field theory have
been constructed and intensively studied~\cite{Schnabl:2007az}--\cite{Noumi:2011kn}. A systematic
procedure to construct analytic solutions for marginal deformations to all orders in the
deformation parameter has been presented in~\cite{Kiermaier:2007vu,Kiermaier:2007ki}. On the other
hand, there has been a puzzle associated with marginal deformations in the level-truncation
analysis in Siegel gauge carried out in~\cite{Sen:2000hx}, where the deformation parameter was
treated nonperturbatively. This is the problem we discuss in this paper. Let us first explain the
setup and describe the puzzle.

When open bosonic string theory on a D25-brane is compactified on a circle of the self-dual radius
$\sqrt{\alpha'}$, there are three marginal operators in the boundary conformal field theory (CFT)
describing the compactified direction. The three operators are given by
\begin{equation}
V_1 (t) \equiv
\sqrt{2} : \cos \frac{X (t)}{\sqrt{\alpha'}} : \,, \quad
V_2 (t) \equiv
\sqrt{2} : \sin \frac{X (t)}{\sqrt{\alpha'}} : \,, \quad
V_3 (t) \equiv
\frac{i}{\sqrt{2 \alpha'}} \, \partial_t X (t) \,,
\end{equation}
where $X (t)$ is the coordinate of the compactified direction. The marginal deformation by the
operator $V_3$ corresponds to turning on a constant mode of the gauge field on the D-brane and is
exactly marginal. Since the direction of the coordinate $X$ is compactified, the one-dimensional
moduli space of this deformation is periodic. These three operators are related by the enhanced
$SU(2)$ symmetry at the self-dual radius, and so the marginal deformations by the other two
operators $V_1$ and $V_2$ are also exactly marginal and the moduli space for each deformation is
periodic. It is known that the deformation by the operator $V_1$ changes the original Neumann
boundary condition to a Dirichlet boundary condition at a special point of the moduli space. In
other words, the original D25-brane is deformed to a D24-brane by this marginal deformation.

We thus expect that the equations of motion of open string field theory formulated around a
D25-brane compactified on the self-dual radius have a one-parameter family of solutions associated
with each of the three marginal deformations, and the moduli space is periodic for each case.
In~\cite{Sen:2000hx} Sen and Zwiebach constructed solutions using level truncation in Siegel gauge.
They computed the effective potential of the massless mode associated with the marginal deformation
by solving the equations of motion for other fields. They found that the effective potential is
approximately flat and becomes flatter as the level of the approximation is increased, which is in
accord with the expectation that we have a one-parameter family of solutions. However, they also
found that the branch of the effective potential is truncated at a finite distance from the origin.
The effective potential does not exist beyond that point, and this result seems stable as the
truncation level is increased. It has been a long-standing problem to understand the nature of this
phenomenon. First of all, the periodicity in the moduli space is obscure. In particular, it is
important to know whether or not the point of the moduli space corresponding to a Dirichlet
boundary condition in the case of the deformation by $V_1$ is within the branch.

In order to investigate this problem it is helpful to construct a map between the deformation
parameter $\lambda_{\rm BCFT}$ of boundary CFT and the parameter $\lambda_{\rm SFT}$ labeling
solutions of open string field theory. However, this has in general been a difficult problem.
In~\cite{Sen:2000hx} Sen and Zwiebach attempted to obtain information on the map in the case of the
deformation associated with $V_1$ by slightly increasing the radius $R$ of the compactification.
The marginal deformation becomes a relevant deformation and the effective potential develops a
local minimum corresponding to tachyon condensation. They used the location of the local minimum to
identify the point of the moduli space corresponding to a Dirichlet boundary condition. However,
the extrapolation to the self-dual radius $R \to \sqrt{\alpha'}$ is not smooth, and they were not
able to obtain a definite conclusion.

Later Sen developed a different method to construct a map between $\lambda_{\rm BCFT}$ and
$\lambda_{\rm SFT}$ in~\cite{Sen:2004cq}. It is based on the energy-momentum tensor in spacetime.
Its dependence on $\lambda_{\rm BCFT}$ can be calculated from the boundary state. Its dependence on
$\lambda_{\rm SFT}$ was calculated from the dependence of the effective potential on the
compactification radius $R$. The two results were combined and a map between $\lambda_{\rm BCFT}$
and $\lambda_{\rm SFT}$ was constructed based on numerical results by level truncation up to
level~4.\footnote{
The analysis was extended to higher levels by Andr\'{e} Kurs
in his Senior Thesis at Princeton University~\cite{Kurs}.}
In order to use this method, however, it is necessary to calculate the dependence of the
effective potential on the compactification radius $R$. Furthermore, this method cannot be used for
arbitrary marginal deformations. For instance, it cannot be used directly for the marginal
deformation by $V_3$ because the relevant component of the energy-momentum tensor does not depend
on the deformation parameter, although we can use the $SU(2)$ symmetry to convert the result for
$V_1$ into that for $V_3$.

In this paper we present a new approach to the construction of a map between the deformation
parameter $\lambda_{\rm BCFT}$ of boundary CFT and the parameter $\lambda_{\rm SFT}$ labeling
solutions of open string field theory for marginal deformations. It is based on a relation between
gauge-invariant observables discovered in~\cite{Hashimoto:2001sm, Gaiotto:2001ji} and the closed
string tadpole on a disk conjectured by Ellwood~\cite{Ellwood:2008jh}. If the closed string tadpole
depends on the marginal deformation of boundary CFT, we can construct a map between $\lambda_{\rm
BCFT}$ and $\lambda_{\rm SFT}$ by calculating the gauge-invariant observable for the solutions of
open string field theory. This method can in principle be used for any marginal deformation if
there is a closed string tadpole which depends on the deformation parameter. Furthermore, the
calculation on the string field theory side is much simpler than that of the method
in~\cite{Sen:2004cq}. The gauge-invariant observables depend linearly on the open string field and
can be calculated by contracting the open string field and the identity state with an additional
insertion of an on-shell closed string vertex operator.

The paper is organized as follows. In section~\ref{s-Marginal} we briefly review the construction
of the solutions for marginal deformations by Sen and Zwiebach in level truncation. In
section~\ref{s-Ellwood} we apply the conjecture by Ellwood to the solutions for marginal
deformations and explain how one can relate the parameter of the solutions in open string field
theory with that of the corresponding boundary CFT. In section~\ref{s-level4} we illustrate our
computation using level truncation up to level~4. In section~\ref{s-level12} we summarize our main
numerical results for solutions with various values of $\lambda_{\mathrm{SFT}}$ obtained up to
level~12. Our results are consistent with the picture that the finite range of the branch covers
just one fundamental domain of the periodic moduli space. Further supporting tables and plots are
presented in appendix~\ref{a-tables}.

\section{Marginal deformations in Siegel gauge}
\label{s-Marginal}
\setcounter{equation}{0}

In this section we briefly review the result by Sen and Zwiebach in~\cite{Sen:2000hx} for marginal
deformations in open string field theory using level truncation. We consider open string field
theory for a D25-brane in a 26-dimensional flat spacetime with one of the spatial directions
compactified on a circle of the self-dual radius and the marginal deformation generated by $V_1
(t)$.

The important features can already be seen at level~1.
The string field truncated to level~1 is given by
\begin{equation}
t_0 \, c_1 \ket{0}
+ t_1 : \cos \frac{X(0)}{\sqrt{\alpha'}} : c_1 \ket{0} \,,
\end{equation}
and the potential for the modes $t_0$ and $t_1$
with the normalization used in~\cite{Sen:2000hx} is
\begin{equation}
V (t_0, t_1) = - \frac{1}{2} \, t_0^2
+ \frac{27 \sqrt{3}}{64} \, t_0^3
+ \frac{3 \sqrt{3}}{8} t_0 \, t_1^2 \,.
\end{equation}
Let us derive the effective potential for $t_1$
by solving the equation of motion for $t_0$.
Since the equation of motion for $t_0$ is a quadratic equation,
there are two solutions:
\begin{equation}
t_0^M = \frac{4}{81 \sqrt{3}}
\left( - \sqrt{64 - \frac{729}{2} \, t_1^2} + 8 \right) \,, \qquad
t_0^V = \frac{4}{81 \sqrt{3}}
\left( \sqrt{64 - \frac{729}{2} \, t_1^2} + 8 \right) \,.
\end{equation}
The superscript $M$ denotes the {\it marginal} branch.
The solution $t_0^M$ vanishes as $t_1 \to 0$
so that this branch is associated with the marginal deformation.
The superscript $V$ denotes the {\it vacuum} branch.
The solution $t_0^V$ is associated
with the branch connected to the tachyon vacuum solution.
We can see from these expressions that
we do not have real solutions when
\begin{equation}
64 - \frac{729}{2} \, t_1^2 < 0 \,.
\end{equation}
The critical values $\pm \bar{t}_1$ for $t_1$ are thus given by
\begin{equation}
\label{lambda_crit}
\bar{t}_1 = \frac{8 \sqrt{2}}{27}
\simeq 0.419 \,.
\end{equation}
The two branches meet at these critical values.
Since the effective potential is even with respect to $t_1$,
from now on we focus our attention on the region $t_1 \ge 0$.

Corresponding to the two solutions $t_0^M$ and $t_0^V$,
there are two branches for the effective potential.
By substituting $t_0^M$ and $t_0^V$ to $V(t_0, t_1)$, we obtain
\begin{equation}
\begin{split}
V^M (t_1) & = \frac{2}{59049} \,
\Bigl( \, -512 + 4374 \, t_1^2
+ ( \, 64 -\frac{729}{2} \, t_1^2 \, )^{\frac{3}{2}} \, \Bigr)
= \frac{27}{128} \, t_1^4 + \frac{6561}{32768} \, t_1^6
+ {\cal O}(t_1^8) \,, \\
V^V (t_1) & = \frac{2}{59049} \,
\Bigl( \, -512 + 4374 \, t_1^2
- ( \, 64 -\frac{729}{2} \, t_1^2 \, )^{\frac{3}{2}} \, \Bigr)
= - \frac{2048}{59049} + \frac{8}{27} \, t_1^2
+ {\cal O}(t_1^4) \,.
\end{split}
\end{equation}
The two branches $V^M (t_1)$ and $V^V (t_1)$
of the effective potential are shown in figure~\ref{figure-1}.

\begin{figure}[h]
\centerline{\hbox{\epsfig{figure=lowest.eps, width=12cm}}}
\caption{The two branches of the effective potential for $t_1$
at level~1.
The solid line is for $2 \pi^2 V^M (t_1)$
and the dashed line is for $2 \pi^2 V^V (t_1)$.}
\label{figure-1}
\end{figure}

We expect an exactly flat potential for the marginal branch because the deformation by $V_1 (t)$ is
known to be exactly marginal. The potential $V^M (t_1)$ at level~1 is not exactly flat, but this is
considered to be an artifact of level truncation. The higher-order analysis in level truncation
shows that the effective potential for $t_1$ on the marginal branch becomes flatter as the
truncation level is increased, as expected. See figure~\ref{fig:energy}.\footnote{The difference
between the energy density $E_{tot}$ computed from the full action and the energy density $E_{kin}$
computed from the kinetic term is proportional to
$t_1$ times the equation of motion for $t_1$. }
It has also been shown
analytically that the coefficient in front of $t_1^4$ in $V^M (t_1)$ vanishes in the limit where
the truncation level becomes infinite~\cite{Berkovits:2003ny}. In this deformation, we can
therefore use $t_1$ as the label $\lambda_{\rm SFT}$ of approximate numerical solutions in level
truncation.

On the other hand, the existence of the critical value of $t_1$ does not seem to be an artifact of
level truncation. The results at higher orders in figure~\ref{fig:energy} show that the critical
value persists and the position of the critical value does not move significantly as the truncation
level is increased. The analysis using level truncation thus indicates that the effective potential
$V^M (t_1)$ on the marginal branch in Siegel gauge becomes exactly flat and terminates at a finite
critical value in the limit where the truncation level becomes infinite.

\begin{figure}[ht]
 \begin{minipage}[b]{0.5\linewidth}
 \centering
 \resizebox{3.2in}{1.7in}{\includegraphics[scale=1]{Energy_tot_MT.eps}}
 \end{minipage}
 \hspace{-0.2cm}
 \begin{minipage}[b]{0.5\linewidth}
 \centering
 \resizebox{3.2in}{1.7in}{\includegraphics[scale=1]{Energy_kin_MT.eps}}
 \end{minipage}
 \caption{\label{fig:energy}{The energy density
 of the marginal and vacuum branches computed from the full action (the left graph)
 and from the kinetic term (the right graph).
 We have chosen $t_1$ as the parameter $\lambda$ of the branches.
 The energy density is measured from the tachyon vacuum
and normalized by the D25-brane tension.
 In this and other figures, the colors follow the spectrum: higher level results are depicted using shorter wavelength colors.
 For explicit color coding, see figure~\ref{Fig:legend}.}}
 \end{figure}

\section{Gauge-invariant observables and the closed string tadpole}
\setcounter{equation}{0}
\label{s-Ellwood}

In~\cite{Ellwood:2008jh} Ellwood conjectured a relation between gauge-invariant observables of open
string field theory discovered in~\cite{Hashimoto:2001sm, Gaiotto:2001ji} and the closed string
tadpole on a disk. The gauge-invariant observable $W ( \phi \,, \mathcal{V} )$ for an open string
state $\phi$ and an on-shell closed string vertex operator $\mathcal{V}$ is defined by the
following correlation function on the upper half-plane:
\begin{equation}
W ( \phi \,, \mathcal{V} )
= \langle \, \mathcal{V} (i) \, f_I \circ \phi (0) \,
\rangle_{\rm UHP} \,.
\end{equation}
Here $\phi (0)$ is the operator corresponding to the state $\phi$
in the state-operator correspondence
and we denoted by $f_I \circ \phi (0)$
the conformal transformation of $\phi (0)$ under the map $f_I (\xi)$ associated with the
identity state:
\begin{equation}
f_I (\xi) = \frac{2 \, \xi}{1 - \xi^2} \,.
\end{equation}
The closed string tadpole on a unit disk
defined in~\cite{Ellwood:2008jh} is
\begin{equation}
\label{tadpole-disk}
\mathcal{A}^{\rm disk} ( \mathcal{V} )
= -\frac{e^{-i \theta}}{2 \pi i} \,
\langle \, \mathcal{V} (0) \, c (e^{i \theta}) \,
\rangle_{\rm disk} \,.
\end{equation}
This should be independent of $\theta$.
The unit disk with a complex coordinate $w$ can be mapped
to the upper half-plane of $z$
by the following conformal transformation:
\begin{equation}
z = i \, \frac{1-w}{1+w} \,.
\end{equation}
The correlation function on a unit disk in~(\ref{tadpole-disk})
can be mapped under this conformal transformation to
\begin{equation}
\mathcal{A}^{\rm disk} ( \mathcal{V} )
= - \frac{1}{\pi} \,
\frac{1}{1 + t^2} \,
\langle \, \mathcal{V} (i) \,
c (t) \, \rangle_{\rm UHP}
\quad \text{with} \quad
t = \tan \frac{\theta}{2} \,.
\end{equation}
It is easy to confirm that
$\mathcal{A}^{\rm disk} ( \mathcal{V} )$
is independent of $t$
when the ghost part of $\mathcal{V}$ is $c \tilde{c}$.
Ellwood conjectured the following relation:
\begin{equation}
W ( \Psi, \mathcal{V} )
= \mathcal{A}^{\rm disk}_\Psi ( \mathcal{V} )
- \mathcal{A}^{\rm disk}_0 ( \mathcal{V} ) \,,
\label{Ellwood-conjecture}
\end{equation}
where $\mathcal{A}^{\rm disk}_0$ is the closed string tadpole
with the original boundary condition
and $\mathcal{A}^{\rm disk}_\Psi$ is the closed string tadpole
with the boundary condition
corresponding to the classical solution $\Psi$.
If the boundary CFT has marginal deformations
labeled by $\lambda_{\rm BCFT}$,
we expect that the equation of motion of open string field theory
has a one-parameter family of solutions
labeled by $\lambda_{\rm SFT}$.
The left-hand side of~(\ref{Ellwood-conjecture})
is a function of $\lambda_{\rm SFT}$.
The closed string tadpole appearing
on the right-hand side of~(\ref{Ellwood-conjecture})
is a function of $\lambda_{\rm BCFT}$.
We can thus obtain a map
between $\lambda_{\rm BCFT}$ and $\lambda_{\rm SFT}$ from
the conjectured relation~(\ref{Ellwood-conjecture}).

We consider the marginal deformation by the cosine potential $V_1 (t)$. We have to choose a closed
string vertex operator $\mathcal{V}$ such that the one-point function $\mathcal{A}^{\rm disk} (
\mathcal{V} )$ has a nontrivial dependence on $\lambda_{\rm BCFT}$. We choose
\begin{equation}
\label{cal-V}
\mathcal{V} = - \frac{2i}{\alpha'} \,
c \tilde{c} \partial X \bar{\partial} X \,.
\end{equation}
The normalization of $\mathcal{V}$ is the same as
that in~\cite{Kawano:2008ry}.
The dependence of $\mathcal{A}^{\rm disk} ( \mathcal{V} )$ on $\lambda_{\rm BCFT}$
can be easily determined using the boundary state~\cite{Recknagel:1998ih},
as studied, for example, in~\cite{Sen:2002nu}.
We have
\begin{equation}
\mathcal{A}^{\rm disk} ( \mathcal{V} )
\propto \cos \, ( \, 2 \pi \lambda_{\rm BCFT} \, ) \,.
\end{equation}
It is a periodic function of $\lambda_{\rm BCFT}$,
and the Dirichlet boundary condition corresponds to
the point $\lambda_{\rm BCFT} = 1/2$.
The overall normalization of $\mathcal{A}^{\rm disk} ( \mathcal{V} )$
can be determined by evaluating the one-point function
of $\mathcal{V}$ on a disk with the Neumann boundary condition,
which corresponds to $\lambda_{\rm BCFT} = 0$.
Our normalization of correlation functions on the upper half-plane is
\begin{equation}
\langle \, c(z_1) \, c(z_2) \, c(z_3) \, \rangle_{\rm UHP}
= (z_1 - z_2) (z_1 - z_3) (z_2 - z_3) \,.
\end{equation}
Here and in what follows
we divide correlation functions
by the spacetime volume factor.
The one-point function of $\mathcal{V}$
on the upper half-plane with the Neumann boundary condition is given by
\begin{equation}
-\frac{1}{\pi} \,
\frac{1}{1 + t^2} \,
\Bigl( \, -\frac{2i}{\alpha'} \, \Bigr)
\langle \, c \tilde{c} \partial X \bar{\partial} X (i) \,
c (t) \, \rangle_{\rm UHP}
= -\frac{1}{2 \pi} \,.
\end{equation}
This fixes the overall normalization of $\mathcal{A}^{\rm disk} ( \mathcal{V} )$
to give
\begin{equation}
\mathcal{A}^{\rm disk} ( \mathcal{V} )
= -\frac{1}{2 \pi} \,
\cos \, ( \, 2 \pi \lambda_{\rm BCFT} \, ) \,.
\end{equation}
Using the relation~(\ref{Ellwood-conjecture}) we have
\begin{equation}
\label{W-lambda_BCFT}
W ( \Psi, \mathcal{V} )
= -\frac{1}{2 \pi} \, \Bigl[ \,
\cos \, ( \, 2 \pi \lambda_{\rm BCFT} \, ) - 1 \,
\Bigr] \,.
\end{equation}
It is convenient to introduce
\begin{equation}
\label{W_XX}
W_{XX}
\equiv 1 - 2 \pi \, W ( \Psi, \mathcal{V} )
\end{equation}
for $\mathcal{V}$ given by (\ref{cal-V}).
Then we have
\begin{equation}
W_{XX}
= \cos \, ( \, 2 \pi \lambda_{\rm BCFT} \, ) \,.
\label{map}
\end{equation}
As we mentioned before, the left-hand side is
a function of $\lambda_{\rm SFT}$ labeling the solutions so that we can
derive a map between $\lambda_{\rm SFT}$ and $\lambda_{\rm BCFT}$ from this relation. In
particular, the value of the gauge-invariant observable for the solution corresponding to the
Dirichlet boundary condition $\lambda_{\rm BCFT} = 1/2$ is given by $W_{XX} = -1$.

\section{Evaluation of the gauge-invariant observables}
\setcounter{equation}{0}
\label{s-level4}

In this section we illustrate our computation of the gauge-invariant observables
for the solutions corresponding to the cosine deformations reviewed in section~\ref{s-Marginal}.
The solutions in Siegel gauge
were constructed by Sen and Zwiebach in~\cite{Sen:2000hx}
up to level~4.
We expand the string field in level $\ell$ as follows:
\begin{equation}
| \, \Psi \, \rangle
= \sum_{\ell = 0}^\infty | \, \Psi^{(\ell)} \, \rangle \,.
\end{equation}
The expressions up to level~4
in the notation used in section~3 of~\cite{Sen:2000hx}
are given by
\begin{equation}
\begin{split}
\label{psi}
| \, \Psi^{(0)} \, \rangle
& = t_0 \, c_1 \ket{0} \,, \\
| \, \Psi^{(1)} \, \rangle
& = t_1  : \cos \frac{X(0)}{\sqrt{\alpha'}} : c_1 \ket{0} \,, \\
| \, \Psi^{(2)} \, \rangle
& = \bigl( \, u_0\,c_{-1}\, b_{-1} + v_0\, L_{-2}^X
+w_0\,L_{-2}^{\prime}  \, \bigr) \, c_1\, \ket{0} \,, \\
| \, \Psi^{(3)} \, \rangle
& = \bigl( \, u_1\, c_{-1}\, b_{-1} + v_1\, L_{-2}^X +w_1\,L_{-2}^{\prime}
+ z_1 L_{-1}^X \,L_{-1}^X \, \bigr) \,\ket{\varphi_t}  \quad
\text{with} \quad
\ket{\varphi_t} = {}: \cos \frac{X(0)}{\sqrt{\alpha'}} : c_1 \ket{0}\,, \\
| \, \Psi^{(4)} \, \rangle
& = \tilde g \, \ket{p_4} + t_2 \, \ket{\chi}
+\Bigl[ \, a\,L_{-4}^X + \bar a\, L_{-4}^{\prime}
+b\, L_{-2}^X\, L_{-2}^X +\bar b\, L_{-2}^{\prime}\, L_{-2}^{\prime}
+\hat b\, L_{-2}^{\prime}\, L_{-2}^X \\
& \qquad \quad
+ c \, c_{-3}\, b_{-1}
+ d \, b_{-3}\, c_{-1}+ e\, b_{-2}\, c_{-2}
+ ( \, f \, L^X_{-2} +\bar f\, L^{\prime}_{-2} \, ) \, c_{-1}\, b_{-1} \, \Bigr] \,
c_1 \ket{0} \,,
\end{split}
\end{equation}
where
$L_n^{X}$ are the Virasoro generators associated with the field $X (t)$ describing the compactified
direction and $L_{n}^{\prime}$ are the Virasoro generators associated with the rest of the matter
fields. The primary field $\ket{p_4}$ at level~4 is given by
\begin{equation}
\ket{p_4}
=\bigl( \, \alpha_{-3}^X\, \alpha_{-1}^X
-\frac{3}{4}\,(\alpha_{-2}^X)^2
-\frac{1}{2}\,(\alpha_{-1}^X)^4 \, \bigr) \, c_1\, \ket{0} \,,
\end{equation}
and $\ket{\chi}$ is
\begin{equation}
\ket{\chi}= {}: \cos \frac{2\, X(0)}{\sqrt{\alpha'}}  : \ket{0} \,.
\end{equation}

The gauge-invariant observable $W ( \phi, \mathcal{V} )$
can be expressed as a BPZ inner product of $\phi$ with an open string state.
For $\mathcal{V}$ given by (\ref{cal-V}),
the gauge-invariant observable is
\begin{equation}
W ( \phi, \mathcal{V} )
= \langle \, \Phi_{XX} \,, \phi \, \rangle \,,
\end{equation}
where the open string state $\Phi_{XX}$ has been constructed in~\cite{Kawano:2008ry} and is given
by
\begin{equation}
\ket{\Phi_{XX}}
= \left (\frac{1}{4}
-4\sum_{n,\,m=1}^{\infty}
i^{m-n}\,m\,n\,\alpha^X_{-m}\,\alpha _{-n}^X \right )
e^{E}c_0\,c_1\,\ket{0}
\end{equation}
with
\begin{equation}
E=-\frac{1}{2} \sum_{n=1}^\infty \,
\frac{(-1)^n}{n}\alpha_{-n}\cdot \alpha_{-n}
+ \sum _{n=1}^\infty \, (-1)^n\,c_{-n}\,b_{-n}\,.
\end{equation}
We denote the gauge-invariant observable
$W ( \Psi, \mathcal{V} )$ truncated to level~$\ell$
by $W^{(\ell)} ( \Psi, \mathcal{V} )$.
The quantity $W^{(4)} ( \Psi, \mathcal{V} )$
is given by
\begin{equation}
\begin{split}
\label{W^(4)}
W^{(4)} ( \Psi, \mathcal{V} )
&= \sum_{\ell = 0}^4 \,
\langle \, \Phi_{XX}, \Psi^{(\ell)} \, \rangle\\
&= \frac{1}{4} \, t_0
+ \frac{1}{4} \, u_0 - \frac{31}{8} \, v_0 + \frac{25}{8} \, w_0
+ \frac{31}{4} \, a -\frac{25}{4} \, \bar a
+ \frac{963}{16} \, b + \frac{675}{16} \, \bar b
-\frac{775}{16}\, \hat b \\
& \quad~ + \frac{1}{4}\, e
-\frac{31}{8}\, f +\frac{25}{8}\, \bar f +192 \, \tilde g \,.
\end{split}
\end{equation}
Note that there are no contributions from
$\Psi^{(1)}$ and $\Psi^{(3)}$
because of momentum conservation.

In the numerical solution constructed by Sen and Zwiebach in~\cite{Sen:2000hx}, the component
fields $t_0$, $t_1$, $u_0$, $\ldots$ are given as functions of $t_1$. We therefore obtain $W^{(4)}
( \Psi, \mathcal{V} )$ for the solution as a function of $t_1$, which we are using as the label
$\lambda_{\rm SFT}$. This is how we compute the gauge-invariant observable as a function of
$\lambda_{\rm SFT}$ in level truncation. We then compare the result with the
expression~(\ref{W-lambda_BCFT}) of the gauge-invariant observable as a function of $\lambda_{\rm
BCFT}$ to find a relation between $\lambda_{\rm SFT}$ and $\lambda_{\rm BCFT}$ numerically.

\section{Numerical evaluation to level 12}
\setcounter{equation}{0}
\label{s-level12}

The method illustrated in the preceding section at level~4 can be extended to higher levels. We had
generated all required vertices using the conservation laws of \cite{Rastelli:2000iu} and then
solved the equations of motion
by Newton's iterative method.\footnote{
Newton's method requires
a choice of an initial approximate solution, which is then iteratively improved to any desired
accuracy. We define it as $\frac{\|\Psi^{(i)}-\psi^{(i-1)}\|}{\|\psi^{(i)}\|}$ and stop the
iteration when we reach $10^{-8}$.
The  norm is defined, as usual, by a square root of the sum of coefficients squared
in the basis formed by $\alpha$ oscillators of the compact direction,
$bc$ oscillators of the ghost sector,
and the Virasoro generators in the rest of the matter sector.
As the starting point we either used results for the solution
with the same $\lambda$ found at lower level or, sometimes more efficiently, used a different
solution with a neighboring value of $\lambda$ at the same level.
 } With
the help of a computer cluster we were able to perform all our computations to level~$(12, 36)$.

To avoid possible confusion between the marginal and vacuum branches, but also out of curiosity, we studied numerically both branches for various values of $\lambda$, up to the critical
level-dependent $\lambda_{crit}$, where both branches meet. The value of $\lambda_{crit}$ at level
1 is given by the simple analytic formula (\ref{lambda_crit}). At higher levels we employed the
interval bisection method and studied
whether Newton's method
converges for given $\lambda$.
The results are given in table \ref{t-lambda_crit}.
\begin{table}[h]
\begin{equation}
\begin{array}{|c|c|r|c|}\hline
L & \lambda_{crit} &L & \lambda_{crit}\\\hline
1  & 0.419026 & 7  & 0.469239\\
2  & 0.454866 & 8  & 0.468645\\
3  & 0.467900 & 9  & 0.468623\\
4  & 0.468614 & 10 & 0.468160\\
5  & 0.469761 & 11 & 0.468109\\
6  & 0.469100 & 12 & 0.467748\\\hline
\end{array}\nonumber
\end{equation}
\caption{Results for $\lambda_{crit}$ at levels $L=1,\ldots,12$. The values are rounded down to six
digits, so that the solution still exists for these values. It can be found within about 5
iterations if we use as the starting point
for Newton's method
the highest $\lambda$ solution
found in the previous steps at the same level. Increasing the last digit by one,
Newton's method
would not anymore converge within~50 iterations at lower levels and at least~10 iterations
at level~12.} \label{t-lambda_crit}
\end{table}
Interestingly, the critical value grows at low levels and is largest at level 5, and after one
oscillation it monotonically decreases. Simple linear extrapolation in $1/L$ to the infinite level
gives the critical value of about $0.466$.
Values of the energy and the gauge-invariant observables
at $\lambda_{crit}$ for each given level $L=1,\ldots,12$ are presented
in table \ref{t-critical} of appendix~\ref{a-tables}.

We have constructed the marginal and vacuum branches for about 15 different values of $\lambda$.
Figure \ref{fig:energy}, presented in section~\ref{s-Marginal}, shows the energy density of the
marginal and vacuum branches meeting at $\lambda_{crit}$. The quantity $W_{XX}$ defined
in~(\ref{W_XX}) on the marginal branch for various levels is plotted in the left graph of
figure~\ref{fig:WXX_M}. Its dependence on the level is rather erratic, so we tried to apply
Pad\'{e}-Borel resummation of contributions of different levels for each solution. The right graph
of figure~\ref{fig:WXX_M} shows the resulting improved values of $W_{XX}$, together with a simple
extrapolation to the infinite level. More data are presented in appendix \ref{a-tables}.

\begin{figure}[ht]
 \begin{minipage}[b]{0.5\linewidth}
 \centering
 \resizebox{3.2in}{1.7in}{\includegraphics[scale=1]{WXX_M.eps}}
 \end{minipage}
 \hspace{-0.2cm}
 \begin{minipage}[b]{0.5\linewidth}
 \centering
 \resizebox{3.2in}{1.7in}{\includegraphics[scale=1]{WXX_PB_M.eps}}
 \end{minipage}
 \caption{\label{fig:WXX_M}{
 The values of $W_{XX}$ at levels $L=2,\ldots,12$ (left) and their Pad\'{e}-Borel improvements at levels $L=4,\ldots,12$, together with a fit to the infinite level in black (right).}}
 \end{figure}

One may notice that starting at level 3 the results at odd levels are close to the values found at
the previous even level. The reason is most likely that at even levels new fields with zero
momentum appear, and they seem to have larger influence than the fields with nonzero momentum. The
similarity of odd and even levels is even greater for $W_{XX}$ because it receives contributions
only from fields at even levels, meaning that any change at odd levels can arise only from the
change of the solution and not from
adding new fields.

Table \ref{t-observables} summarizes linear extrapolations in $1/L$ of various quantities of
interest to the infinite level
together with the estimated statistical error.
We estimate the error by considering extrapolations
using different numbers of data points
and computing the sample standard
deviation for such results. For $W_{XX}$ we use the Pad\'{e}-Borel resummed version. We also
computed $W_{00}$ defined by
\begin{equation}
W_{00} = 1 +2 \pi \, W ( \Psi, \mathcal{V} )
\end{equation}
with
\begin{equation}
\mathcal{V} = - \frac{2i}{\alpha'} \,
c \tilde{c} \partial X^0 \bar{\partial} X^0 \,.
\end{equation}
Up to approximately $\lambda=0.325$ the extrapolated values of the energy density $E_{tot}$
computed from the full action, the energy density $E_{kin}$ computed from the kinetic term, and
$W_{00}$ are consistent with the expected value 1. Apparently, for higher values of $\lambda$ our
method underestimates the error of the extrapolated values.

\begin{table}
\begin{equation}
\begin{array}{|l|ll|ll|ll|ll|}\hline
 \lambda & E_{tot} & \sigma_{E_{tot}} & E_{kin} & \sigma_{E_{kin}} & W_{00} &\sigma_{W00} & W_{XX} & \sigma_{W_{XX}}\\ \hline
 0.05 & 1. & 6.3 \times 10^{-7} & 1. & 2.1 \times 10^{-7} & 1. & 0.00007 & 0.9493 & 0.0018 \\
 0.1 & 1. & 0.00001 & 1. & 3.4 \times 10^{-6} & 1. & 0.00028 & 0.8051 & 0.0025 \\
 0.15 & 1. & 0.000051 & 1. & 0.000018 & 1.0001 & 0.00066 & 0.5686 & 0.0046 \\
 0.2 & 1. & 0.00016 & 1. & 0.000057 & 1. & 0.0012 & 0.2572 & 0.014 \\
 0.25 & 0.9999 & 0.0004 & 1. & 0.00014 & 0.9996 & 0.002 & -0.1338 & 0.012 \\
 0.3 & 0.9999 & 0.00083 & 0.9997 & 0.00026 & 0.9978 & 0.0029 & -0.5499 & 0.0022 \\
 0.325 & 1. & 0.0011 & 0.999 & 0.00033 & 0.9954 & 0.0034 & -0.7449 & 0.01 \\
 0.35 & 1.0004 & 0.0015 & 0.9972 & 0.00041 & 0.9905 & 0.0039 & -0.9096 & 0.013 \\
 0.375 & 1.0016 & 0.0019 & 0.992 & 0.00052 & 0.9798 & 0.0043 & -1.0251 & 0.027 \\
 0.4 & 1.0049 & 0.0023 & 0.9781 & 0.00035 & 0.9565 & 0.0046 & -1.076 & 0.044 \\
 0.425 & 1.013 & 0.0028 & 0.9433 & 0.0013 & 0.9075 & 0.005 & -1.0649 & 0.059 \\
 0.45 & 1.0321 & 0.0035 & 0.8668 & 0.0046 & 0.8116 & 0.0064 & -1.0562 & 0.076 \\
 0.46 & 1.0451 & 0.0037 & 0.8112 & 0.0067 & 0.7453 & 0.0077 & -1.0778 & 0.085 \\
 0.4675 & 1.0585 & 0.0039 & 0.7124 & 0.019 & 0.6326 & 0.02 & -1.1149 & 0.092 \\\hline
\end{array}\nonumber
\end{equation}
\caption{Extrapolation to the infinite level
and its possible statistical error for the four basic
quantities computed at several values of $\lambda$ on the marginal branch. Up to about
$\lambda=0.325$ the first three observables are within the $3\sigma$ tolerance away from the
expected value of $1$.} \label{t-observables}
\end{table}

\begin{figure}[ht]
 \centering
 \includegraphics[scale=1.4]{WXX_conjecture.eps}
 \caption{\label{fig:conjecture}{Extrapolation to the infinite level
 of the Pad\'{e}-Borel resummed version of $W_{XX}$.
 The black solid line is the infinite-level fit itself, while the
 shaded regions around it indicate $\sigma$ and $3\sigma$
 uncertainty ranges
 in the extrapolation.
 The magenta line shows the (unimproved) data at level~12 for comparison.
 The black dot-dashed line is
the one-parameter fit of the form $\cos\left(2\pi(\lambda + a \lambda^3)\right)$.}}
 \end{figure}

The results for $W_{XX}$ together with its estimated error are presented in~figure \ref{fig:conjecture}.
The one-parameter fit of the form $\cos\left(2\pi(\lambda + a \lambda^3)\right)$ is also shown,
where the constant $a$ has been fitted
so that  the extrapolation to the infinite level
coincides with the
boundary CFT result at lower values of $\lambda$. The best fit value for $a$ is about $0.99$ at
level 4 and grows at higher levels.
For the extrapolation to the infinite level
we find $a \approx 1.48$.
It would be interesting to calculate
this value analytically. The results of $W_{XX}$ qualitatively
agree with the analysis by Sen~\cite{Sen:2004cq} using the energy-momentum tensor and are
consistent with the picture that the finite range of the branch covers one fundamental domain of
the periodic moduli space. While our results do not provide sufficient evidence that the branch
covers {\it exactly} one fundamental domain, we can safely exclude the possibility that the branch
covers the fundamental domain many times. We can also exclude the opposite possibility that the
branch covers, say, less than 75\% of the fundamental domain. If the branch indeed covers exactly
one fundamental domain, it would be interesting to understand why it is the case.

\section{Discussion}
\setcounter{equation}{0}

In this paper we presented a method of constructing a map between the deformation parameter
$\lambda_{\rm BCFT}$ of boundary CFT and the parameter $\lambda_{\rm SFT}$ labeling solutions of
open string field theory for marginal deformations. While we applied our method to the specific
problem of covering moduli for solutions in Siegel gauge, the basic idea is universal and we can
use it for other problems. For example, a similar problem for the system of separated D-branes has
recently been discussed in~\cite{Karczmarek:2012pn, Longton:2012ei} and it would be interesting to
compute the gauge-invariant observables in the system.

If we construct the boundary state
from numerical solutions
following the proposal in~\cite{Kiermaier:2008qu},
we will be able to obtain more information
on the boundary CFT.
However, construction of the boundary state
for numerical solutions seems to be challenging.
In particular, the closed string state
proposed in~\cite{Kiermaier:2008qu}
is conjectured to reproduce the boundary state
up to a BRST-exact term,
and such a BRST-exact term would obscure
the boundary state especially for approximate solutions
constructed numerically.\footnote{
After this paper was submitted to arXiv,
another approach to the construction of the boundary state
was proposed in~\cite{Kudrna:2012re},
which is more suitable for the application to numerical solutions.}

Our results are qualitatively consistent
with the analysis by Sen~\cite{Sen:2004cq}
using the energy-momentum tensor.
This is not unexpected because both methods
are based on coupling to an infinitesimal closed string field.
However, the closed string is described
by an unintegrated vertex operator
in the gauge-invariant observables
and also in the boundary state construction in~\cite{Kiermaier:2008qu},
while the infinitesimal change in the closed string background
is related more closely to an integrated vertex operator.
It would be an interesting problem
to convert the unintegrated vertex operator
to the integrated vertex operator
in the framework of the gauge-invariant observables
and the boundary state construction~\cite{Kiermaier:2008qu}.

\bigskip

\noindent
\section*{Acknowledgments}

\medskip

We would like to thank Ted Erler, Mitsuhiro Kato, Michael Kiermaier, Carlo Maccaferri and Barton Zwiebach for helpful discussions and to Masaki Murata for careful reading of the manuscript.
The access to computing and storage facilities owned by parties and
projects contributing to the Czech National Grid Infrastructure
MetaCentrum, provided under the programme ``Projects of Large
Infrastructure for Research, Development, and Innovations"
(LM2010005), and to CERIT-SC computing facilities provided under the
programme Center CERIT Scientific Cloud, part of the Operational
Program Research and Development for Innovations, reg. no. CZ.
1.05/3.2.00/08.0144 is highly appreciated.
The work of Y.O. was supported in part
by Grant-in-Aid for Scientific Research~(B) No.~20340048
and Grant-in-Aid for Scientific Research~(C) No.~24540254
from the Japan Society for the Promotion of Science (JSPS).
The work of M.S. was supported by the EURYI grant GACR EYI/07/E010 from EUROHORC and ESF.
The work of Y.O. and M.S. was also supported in part
by the M\v{S}MT contract No. LH11106 and
by JSPS and the
Academy of Sciences of the Czech Republic (ASCR)
under the Research Cooperative Program between Japan and the Czech Republic.


\section{Introduction}

Manipulating the molecular rotational degree of freedom remains a
goal of primary interest in photochemistry with applications
extending from chemical reactivity to nanoscale design
\cite{newrabitz,shapiro,rice}. In this framework, molecular
alignment and orientation constitute a well-established topic both
from the experimental and theoretical points of view
\cite{stapelfeldt,seideman}. While the alignment process is now
well understood in the adiabatic \cite{friedrich} or sudden regime
\cite{averbukh2} with recent extensions such as the molecular
classical rotation \cite{averbukh1}, the deflection of aligned
molecules \cite{averbukh3}, the planar alignment
\cite{sugnyplanar1,sugnyplanar2,hcn} or the analysis of the
dissipation effects \cite{seidemandiss,sugny1}, work remains to be
done in order to control and produce molecular orientation with a
sufficient high efficiency. On the theoretical side, several basic
mechanisms have been proposed, built on intuitive or optimal
control strategies \cite{salomon,lapert}. Among others, we can
cite the kick mechanism which consists in a sudden impact to the
molecule by a half-cycle pulse (HCP)
\cite{HCPexp,kick1,kick2,kick3,kick4,kick5}, its combination with
a laser field \cite{averbukhHCP1,averbukhHCP2,daems} or the
$(\omega-2\omega)$ scheme
\cite{KlingCO,tehini,zhang,tehini2,2colorchinois,multicolorchinois,twocolorchinois,coreehyper}.
Due to the efficiency of the first process based on its asymmetric
temporal shape, most of the theoretical works have focused on its
application \cite{arvieu,sugny7,sugny4,sugny3,HCPcoree}. However,
recent studies \cite{ortigoso,sugny6} have pointed out the
inherent experimental difficulties associated to the use of such
pulses, which are distorted when they propagate in free space as
well as through focusing optics \cite{HCPpropa}. This phenomenon,
due to the DC component of the field, makes thus problematic the
experimental implementation of these techniques in the control of
molecular rotation. In this framework, a fundamental question is
whether it is possible to orient linear molecules in the THz
regime by using only zero area laser pulses. These latter do not
contain DC field and are therefore free of these propagation
distortions. This problem has been recently addressed
theoretically
\cite{sugny6,THzchinois,THzchinois2,THzchinois3,chinese} and
experimentally \cite{orientationHCP}, but no systematic study of
the efficiency of this process has been done. Note that this
question is not trivial since the sudden impact approximation
predicts no post-pulse orientation in this regime
\cite{kick2,kick3,sugny6}. We present in this paper a complete
analysis of this control strategy by considering different linear
molecules. We establish under which conditions this process is
efficient and we analyze its robustness with respect to
temperature effects. Two mechanisms leading to molecular
orientation are identified. The first one is valid at low
temperature, while the second process is only efficient for higher
temperatures. In this second non-intuitive control scheme, we show
the positive role of temperature effects in the orientation
mechanism.

The paper is organized as follows. The model system is presented
in Sec. \ref{sec2}. The different numerical results are discussed
in Sec. \ref{sec3}. Conclusions and prospective views are given in
Sec. \ref{sec4}.
\section{The model system}\label{sec2}
We consider the control of a linear polar molecule by a linearly
polarized THz laser field $E(t)$ of zero area. The molecule is
assumed to be in its ground vibronic state. Within the rigid rotor
approximation, the Hamiltonian of the system can be written as
\begin{equation}
H(t)=BJ^2-\mu_0 E(t)\cos\theta,
\end{equation}
where $B$ is the rotational constant, $J^2$ the angular momentum
operator and $\mu_0$ the permanent dipolar moment. We neglect in
this paper the effect of the polarizability components since the
maximum intensity of the electric field remains moderate. The
units used throughout the paper are atomic units unless otherwise
specified. At non zero temperature, the system can be either
described by a density matrix $\rho(t)$ with a dynamics governed
by the von Neumann equation
\begin{equation}\label{eq1}
i\frac{\partial \rho(t)}{\partial t}=[H(t),\rho(t)],
\end{equation}
where $\rho(0)$ is the canonical density operator at thermal
equilibrium, or by a set of wave functions
$|\psi_{J_0,M_0}(t)\rangle$ satisfying each the Schr\"odinger
equation:
\begin{equation}\label{eq1b}
i\frac{\partial |\psi_{J_0,M_0}\rangle}{\partial
t}=H(t)|\psi_{J_0,M_0}\rangle,
\end{equation}
with as initial state $|\psi_{J_0,M_0}(t=0)\rangle
=|J_0,M_0\rangle$. The second representation given by Eq.
(\ref{eq1b}), which is more suited to the control mechanisms
interpretation, will be used in this work. The expectation value
$\langle \cos\theta\rangle$ defined by
\begin{equation}\label{eq2}
\langle \cos\theta\rangle
(t)=\frac{1}{Z}\sum_{J_0=0}^{+\infty}c_{J_0}\sum_{M_0=-J_0}^{M_0=J_0}\langle
\psi_{J_0,M_0}|\cos\theta |\psi_{J_0,M_0}\rangle,
\end{equation}
is taken as a quantitative measure of orientation, with the
weights $c_{J}=e^{-BJ(J+1)/(k_BT)}$ and the partition function
$Z=\sum_{J=0}^\infty\sum_{M=-J}^J c_J$, where $T$ is the
temperature and $k_B$ the Boltzman constant. In the rest of the
paper, we will need to separate the zero temperature contribution
to the thermal one, which are respectively denoted by $\langle
\cos\theta \rangle _0$ and $\langle \cos\theta \rangle _T$. In the
zero temperature response, the sum of Eq. (\ref{eq2}) is carried
out only for $J_0=0$, while for the thermal contribution, the
expectation value is computed over the other values of $J_0$.

The electric field is assumed to be of the form
\begin{eqnarray*}
E(t) &=& E_0f(t)=E_0\cos^2(\pi\frac{t}{\delta})\sin(2\pi ft),\quad t\in [ -\delta/2,\delta/2],\\
E(t) &=& 0\quad\textrm{otherwise},
\end{eqnarray*}
where $E_0$ is the amplitude of the pulse, $\delta$ its duration
and $f$ its central frequency. By symmetry, this field has a zero
area for any values of $\delta$ and $f$. Note that the sudden
impact approximation can be applied if the pulse duration $\delta$
is small with respect to the rotational period
$T_{\textrm{per}}=\pi/B$. This also means that noticeable
orientation can be produced only for sufficient large values of
$\delta$ \cite{kick2,sugny6}.

The Schr\"odinger equation (\ref{eq1b}) can be written as
\begin{equation}
i\frac{\partial |\psi_{J_0,M_0}\rangle}{\partial
\tau}=[J^2-A\cos\theta f(\tau)]|\psi_{J_0,M_0}\rangle,
\end{equation}
where the new dimensionless parameters are defined by:
\begin{eqnarray*}
\tau = Bt,~A = \frac{\mu_0E_0}{B},~F = \frac{f}{B},~D = B\delta,~
\tilde{T} = \frac{Tk_B}{B},
\end{eqnarray*}
with $f(\tau)=\cos^2(\pi\frac{\tau}{D})\sin(2\pi F\tau)$. In these
coordinates, the rotational period becomes $T_{\textrm{per}}=\pi$
and the pulse duration is $D$. These effective parameters
completely describe the dynamical evolution of the system and the
degree of molecular orientation produced. They will give a general
understanding of the orientation mechanism free of any particular
molecule. Note that the rotational constant $B$ is a crucial
parameter in these new coordinates since the effective field,
frequency and duration $(A,F,D)$ depend on $B$. For numerical
applications, we will consider the parameters listed in Table
\ref{tab1}.
\begin{table}\caption{\label{tab1}  Molecular parameters of different molecules used in the numerical computations. Numerical values are taken to be
$E_0=2.19$ MV/cm, $\delta=5$ ps and $f=0.5$ THz.}
\begin{center}
\begin{tabular}{c|c|c|c|c|c}
\hline
  Molecule& OCS &  HF  & LiH & CO & LiCl\\
  \hline \hline
  $B~(\textrm{cm}^{-1})$ &0.2029&20.956&7.513&1.931&1.345\\
\hline
 $ \mu_0~(\textrm{debye})$ &  0.712&1.820&5.88&0.112&6.33\\
\hline
 $A$&117.8497&2.9167&26.2842&1.9479&158.0563\\
\hline
$F$  &13.0823&0.1267&0.3533&1.3746&1.9735 \\
\hline $D$ & 0.1911&19.7371&7.0760&1.8187&1.2668   \\
\hline
 \end{tabular}
\end{center}
\end{table}
\section{Numerical Results}\label{sec3}
We begin our study by a general analysis of the maximum degree of
orientation that can be reached as a function of the rotational
constant $B$ and the temperature $T$. A fictive molecule with a
permanent dipolar moment of 1 debye has been considered in the
computation. The field amplitude is assumed to be $E_0=2.19$ MV/cm
(i.e. with a peak amplitude of 2 MV/cm), which corresponds roughly
to the maximum amplitude of THz pulses actually available
experimentally \cite{wolf}. This dependance is shown in Fig.
\ref{fig1} where two zones of high orientation can be clearly
distinguished: an upper zone, denoted (I), associated to high
values of $B$ and low temperature and a flat zone, denoted (II)
requiring small rotational constants and larger temperatures up to
250 K.
\begin{figure}[htbp]
\centering
\includegraphics[scale=0.6]{fig1.eps}
\caption{(Color online) Maximum orientation as a function of the
rotational constant $B$ and the temperature $T$. The field
parameters are $\delta=5$ ps, $f=0.5$ THz and $E_0=2.19$ MV/cm.
Two regions (I) and (II) of the diagram are delimited by arcs of
ellipses. \label{fig1}}
\end{figure}
If the rotational constant is lower than 2 cm$^{-1}$, we first
notice that a very weak orientation is produced at $T=0$ K. The
rotational period varying as the inverse of $B$, this result can
be explained through the sudden impact approximation
\cite{kick2,kick3,sugny6}. Another standard feature shown in
previous studies is the detrimental effect of temperature on
molecular orientation \cite{seideman,kick2,kick3}. This behavior
can be recovered in the case of region (I) where almost no
orientation is obtained for a temperature larger than 100 K. By
comparison, the orientation observed in zone (II) is rather
unexpected, since the temperature effect becomes positive, no
orientation being produced at very low temperature. One of the
goal of this work will be to explore the basic mechanism at the
origin of this non trivial phenomenon.\\
\textbf{The orientation mechanism.}\\
We first present a general spectral analysis of the control
problem. We use in the following the dimensionless coordinates
introduced in Sec. \ref{sec2}, for which the molecular spectrum is
a discrete spectrum composed of the frequencies $\omega_J=2(J+1)$.
To simplify the discussion, we associate to each frequency a
weight $P$ which is defined as the average of the initial thermal
population of the levels $J$ and $J+1$, i.e. $P=(c_J+c_{J+1})/2$.
The spectrum of the control field is proportional to the one of
the function $f(\tau)$ and can roughly be approximated by a
gaussian spectrum centered at $\nu=F$ with a bandwidth
proportional to $1/D$. Different spectra are represented in Fig.
\ref{fig2} and \ref{fig3} for different values of $D$ and $F$,
respectively. In the case of Fig. \ref{fig2}, only the width of
the field spectrum is changed, while only the central frequency is
modified in Fig. \ref{fig3}. From a qualitative point of view, the
overlap between the molecular and the field spectra is a necessary
condition to produce molecular orientation. Figure \ref{fig2}
displays two cases with different values of $D$, i.e. a pulse with
a large spectrum overlapping many transitions including the first
one, and a pulse with a larger width $D$ corresponding to a
narrower spectrum. The first pulse can produce orientation at zero
and non-zero temperatures, while the second one is only efficient
at non zero temperature. This behavior is confirmed in the top
panel of Fig. \ref{fig4}. Figure \ref{fig3} depicts two different
situations, a first one with a low main frequency which overlaps
with the first frequency spectrum and a second case where the main
frequency is far from the first molecular transition. The first
pulse should work at zero and nonzero temperatures due to the
overlap with other transitions while in the second example, we
cannot expect orientation at zero temperature but only for
$\tilde{T}\neq 0$. The corresponding orientation responses are
given in the bottom panel of Fig. \ref{fig4}. Comparing in Fig.
\ref{fig4} the efficiency of the pulses of Fig. \ref{fig2} and
\ref{fig3}, one observes that at $\tilde{T}\neq 0$, a stronger
orientation is obtained in the first case of Fig. \ref{fig2} due
to a larger overlap of the field spectrum with the population
distribution. In other words, a noticeable orientation is produced
only for a large bandwidth of the laser field.

\begin{figure}[htbp]
{\centering\includegraphics[scale=0.6]{fig2.eps}} \caption{(Color
online) Plot of the spectral distribution of the field at
$\tilde{T}=0$ and $\tilde{T}=50$ ($T=143.9$ K for $B=2$ cm$^{-1}$)
on the top and bottom panels, respectively. The amplitude of the
field spectrum is arbitrary. Control fields parameters are taken
to be $D=1$ ($\delta=2.654$ ps for $B=2$ cm$^{-1}$) for the blue
(black) line and $D=3$ ($\delta=7.963$ ps for $B=2$ cm$^{-1}$) for
the green (gray) line. Other parameters are fixed to $F=2$
($f=0.753$ THz for $B=2$ cm$^{-1}$) and $A=4$. The columns
represent the population distribution $P=(c_J+c_{J+1})/2$ at the
corresponding frequency $\omega=2(J+1)$. The quantities $P$ and
$\omega$ are unitless.\label{fig2}}
\end{figure}
\begin{figure}[htbp]
{\centering\includegraphics[scale=0.6]{fig3.eps}} \caption{(Color
online) Same as Fig. \ref{fig2} but for $F=0.5$ ($f_1=0.188$ THz
for $B=2$ cm$^{-1}$) in blue or dark and $F=2.5$ ($f=0.942$ THz
for $B=2$ cm$^{-1}$) in green or gray. The other parameters are
given by $D=3$ ($\delta=7.96$ ps for $B=2$ cm$^{-1}$) and
$A=4$.\label{fig3}}
\end{figure}
\begin{figure}[htbp]
\centering\includegraphics[scale=0.5]{fig4.eps} \caption{(Color
online) Time evolution of $\langle \cos\theta\rangle$ for the
cases of Fig. \ref{fig2} (top panel) and \ref{fig3} (bottom
panel). The same color code as in Fig. \ref{fig2} and \ref{fig3}
has been used. Solid and dashed lines depict respectively the
orientation at $\tilde{T}=0$ and $\tilde{T}=50$. The field is on
for negative times and switches off at $\tau=0$. The variable
$\tau$ is unitless.\label{fig4}}
\end{figure}

The interpretation of the previous results leads to two basic
mechanisms governing the orientation response. The orientation
created at $T\simeq 0$ K can be associated to a rotational ladder
climbing mechanism from the ground energy level, which consists
here in the successive excitation of neighboring rotational
levels. This process is efficient only if the field spectrum
allows to excite the first rotational frequencies with a
sufficient high intensity. Note that this control scheme has
already been identified in the literature to produce molecular
orientation \cite{salomon}. The rotational population distribution
being displaced to high $J$ levels with increasing temperatures,
this mechanism then loses its efficiency very fast. This
orientation will be called a \emph{zero-temperature orientation}
and can be quantitatively measured by the partial expectation
value $\langle \cos\theta\rangle _0$. In the case of an initial
thermal distribution of rotational states, a new mechanism occurs.
This situation corresponds to molecules with small values of $B$,
i.e. with a quite large frequency $F$. As can be seen in Fig.
\ref{fig3} for $F=2.5$, there is no overlap with the population
distribution at low temperature and thus no orientation. In this
second mechanism at high temperature, the control field excites
simultaneously many rotational frequencies (i.e. with a large
laser bandwidth), creating several rotational waves packets
$|\psi_{J_0,M_0}(t)\rangle$ which interfere constructively to
produce a noticeable orientation, denoted \emph{thermal
orientation}. This latter can be computed through the expectation
value $\langle \cos\theta\rangle _T$. In this scenario, the
orientation reaches a maximum at a temperature different from zero
and presents a slow decrease with increasing temperatures.

These different conclusions can be checked in Fig. \ref{fig5}
where the zero-temperature and the thermal orientation responses
have been plotted. One clearly sees in this figure that the zone
(I) is associated to $\langle \cos\theta\rangle _0$, where only
the initial state $|0,0\rangle$ is considered. The efficiency of
THz pulses in the region (II) can be interpreted as a thermal
orientation such that $\langle \cos\theta\rangle \simeq \langle
\cos\theta\rangle _T$. Note that the sum of the two figures
\ref{fig5} does not give exactly the result of Fig. \ref{fig1} due
to destructive interferences between the two orientation
responses.

The transition from a regime with thermal orientation to a regime
with a zero-temperature orientation can be understood from the
definition of the parameters $(F,D)$. Decreasing the value of $B$
is equivalent to increase the effective frequency $F$ and to
decrease the temporal width $D$ of the pulse. Starting from the
zone (I) where the spectrum of the field overlaps the first
molecular transition, a decrease of the rotational constant $B$
will shift the field spectrum to higher frequencies $F$. The
overlap with the first frequency is removed and no zero
temperature orientation is possible. At the same time, as a
consequence of the increase of the width of the field spectrum,
more and more molecular transitions can be excited by the pulse,
leading thus to a thermal orientation.

\begin{figure}[htbp]
\centering\includegraphics[scale=0.6]{fig5a.eps}
\centering\includegraphics[scale=0.6]{fig5b.eps} \caption{(Color
online) Same as Fig. \ref{fig1} but for the zero temperature
contribution $\langle \cos\theta\rangle_0$ to the molecular
orientation (top) and for the thermal one $\langle
\cos\theta\rangle_T$ (bottom).\label{fig5}}
\end{figure}

\begin{figure}[htbp]
\centering\includegraphics[scale=0.5]{fig6.eps} \caption{(Color
online) The top panel represents the evolution of the maximum of
orientation as a function of the temperature for the molecules of
Table \ref{tab1}. The field parameters of Fig. \ref{fig1} have
been used. The contribution of the zero temperature (red or dark
gray) and thermal (green of light gray) wave packets have been
plotted for the LiCl (solid lines) and HF (dashed lines) molecules
in the bottom panel. The blue (black) lines corresponds to the
total orientation response. The temperature $T$ is expressed in
Kelvin.\label{fig6}}
\end{figure}
The evolution of the orientation for the molecules of Table
\ref{tab1} is displayed in Fig. \ref{fig6}. All the molecules
considered belong to the zone (I), except for the CO and OCS
molecules which present an orientation close to 0 at $T\simeq 0 K$
and an orientation of the order of 0.15  and 0.1 at $T\simeq 10 K$
and $T\simeq 150 K$, respectively. The efficiency of the control
scheme at high temperature becomes even better than the one
obtained for other molecules such as HF or LiH, which is very good
only at low temperature. Note that the CO and LiCl molecules have
quite the same parameters $F$ and $D$. The difference of
orientation response comes from the difference of two order of
magnitude in their effective amplitude $A$. A ladder climbing
process at $T=0$ K is possible for LiCl since the tails of the
spectrum are sufficient, while it is not the case for CO. The two
molecules present however a thermal orientation. We also point out
that the orientation of OCS at high temperature is due to its low
value of $D$ (very broad spectrum) combined with a high value of
$F$ and a large effective amplitude of the field. In Fig.
\ref{fig6}, note the smooth evolution of the orientation, except
in the cases of the LiH and LiCl molecules where a slope change
occurs. This feature can be explained by a transition from a
zero-temperature orientation to a thermal one as shown in Fig.
\ref{fig6} (bottom). It can be seen that the non smooth point of
the curve for the LiCl molecule can be viewed as a limit point for
which the thermal orientation becomes predominant with respect to
the zero-temperature one.

\begin{figure}[htbp]
\centering\includegraphics[scale=0.5]{fig7.eps} \caption{(Color
online) Maximum orientation as a function of the field amplitude
$E_0$ and the temperature $T$. The other field parameters are
$\delta=5$ ps and $f=0.5$ THz. The LiCl molecule is considered for
the computations. The temperature $T$ is expressed in
Kelvin.\label{fig7}}
\end{figure}
We finally explore the dependance of the orientation with respect
to the field strength. A global description of this sensitivity as
a function of the temperature is displayed in Fig. \ref{fig7}. At
zero temperature, we observe a quadratic increase of the
orientation up to $E_0\simeq 0.6$ MV/cm. Higher amplitudes lead to
a more chaotic evolution of the orientation characterized by the
occurrence of maxima and minima. At high temperature, the same
erratic distribution of the orientation response can be seen with,
however, a larger periodicity. We thus conclude that the molecular
orientation is less sensitive to amplitude field changes at
non-zero than at zero temperature. We also observe that the
thermal orientation can be produced at a lower intensity than the
zero temperature one.

\section{Conclusion}\label{sec4}
This paper has focused on the use of THz laser pulses for
controlling the orientation dynamics of linear molecules.
Numerical tests have shown the efficiency of the proposed control
scheme, even at high temperature for some molecules. Until now,
most of the works have envisaged the production of molecular
orientation with short laser pulses characterized by a non zero
time average of the electric field as, e.g., with the use of HCPs.
Indeed, the sudden impact approximation predicts no post-pulse
orientation when this average vanishes. Opposite to this accepted
fact, we have shown that a significant orientation can be obtained
in two different regimes. The first one corresponds to a standard
situation with low temperature and high rotational constants. In
the second case, for an adequate choice of the pulse parameters,
we have established that the temperature plays an active role in
the production of molecular orientation. This study calls for
further experimental investigation of the use of such laser pulses
in order to complete the initial work of Ref.
\cite{orientationHCP}. In particular, one objective could be to
demonstrate experimentally the existence of thermal orientation in
molecules such that OCS or CO.

\noindent \textbf{Acknowledgment}\\
We are grateful to E. Hertz for discussions.