\section{Introduction}

Along with the dramatic achievement in nanotechnology,
such as molecular beam epitaxy, atomic scale fabrication or advanced
lithography, mesoscopic physics and nanoelectronics are
undergoing a rapid development\cite{IBM,Buot}.
It has been a strong and definite trend in the miniaturization
of integrated circuits and components towards atomic-scale
dimensions\cite{Garcia} for the electronic device community.
When transport dimension reaches a characteristic dimension,
namely, the charge carrier inelastic coherence length, one must address
not only  quantum mechanical property but also
the discreteness of electron charge. Thus a correct quantum
theory is indispensable for the device physics in integrated circuits of
nanoelectronics.
The quantization of the circuit was carried out
\cite{Louisell} just likes that of a harmonic oscillator.
This only results energy quantization.
In fact, a different kind of fluctuation in
mesoscopic system
is due to the quantization of electronic charge.
We studied the quantization of electric circuit of LC-design
under the consideration of the discreteness of electric charge\cite{LiCh1}.

In this talk, we demonstrate the quantum theory for mesoscopic
electric circuits and it applications.
After recalling the Kirchoff's law in classical LC-design in next
section, we quantize the LC-design by considering the fact of
charge discreteness in section \ref{sec:quantization}, where
a finite-difference Schr\"{o}dinger equation for the mesoscopic electric
circuit is proposed.
In section \ref{sec:uncertainty},
the uncertainty relation for charge and current
is discussed.
In section \ref{sec:LC},
the finite-difference Schr\"{o}dinger equation for a mesoscopic
circuit of LC-design is turned to  Mathieu equation in
`p-representation' and solved exactly.
In section \ref{sec:L},
A gauge field
is introduced and a formula for persistent current which is a periodic
function of the magnetic flux is obtained.
In section \ref{sec:C}, The pure C-design is studied and the Coulomb
blockade is formulated.
In section \ref{sec:time}, we study the time evolution of the system,
particularly
the L-design in the presence of time-dependent source.
In section \ref{sec:dissipative},
we suggest an approach to include the dissipative effects in the circuit.
In section \ref{sec:discussion}, we summarize the main results and give some
discussions.

\section{The classical LC-design}\label{sec:classical}

A classical  non-dissipative electric
circuit of LC-design fulfills the
Kirchoff's law, {\it i.e.}

\[
\dot{q} = \frac{\partial H }{\partial p }, \,\,
\dot{p} = - \frac{\partial H }{ \partial q }
\]
with $H(t) = {\displaystyle
\frac{ 1 }{ 2L } p^2 + \frac{ 1 }{2C}q^2 + \varepsilon(t) q
} $ as the Hamiltonian.
The variable $q$ stands for the electric charge instead of
the conventional `coordinate', while its conjugation
$
{\displaystyle
p(t) = L dq /d t
} $
represents (apart from a factor L ) the  electric current instead of
the conventional `momentum'.

\section{Quantization}
\label{sec:quantization}

In order to take account of the discreteness of
electronic charge.
we must impose that
the eigenvalues of the self-adjoint operator
$\hat{q}$ take discrete
values\cite{LiCh1}, i.e.
\begin{equation}
\hat{q} | q > = n q_e | q >
\label{eq:c}\end{equation}
where $n \in {\sf \, Z \!\!\! Z \, } $ (set of integers) and
$q_e  = 1.602 \times 10^{-19}$ coulomb,
the elementary electric charge.  Obviously, any eigenstate of
$\hat{q}$ can be specified by an integer. This allows us to
introduce a minimum `shift operator'
$
\hat{Q} := \exp [iq_e \hat{p} /\hbar],
$
which  satisfies
\begin{eqnarray}
[ \hat{q}, \hat{Q} ] = - q_e \hat{Q}, \;\;
[ \hat{q}, \hat{Q}^{+} ] = q_e \hat{Q}^{+}, \nonumber \\[2mm]
\hat{Q}^{+} \hat{Q} = \hat{Q} \hat{Q}^{+} = 1.\hspace{10mm}
\label{eq:d}\end{eqnarray}
In ref.\cite{LiCh1} we
obtained the following
finite-difference Schr\"{o}dinger equation,
\begin{equation}
\left[
-\frac{\hbar^2}{2q_e L}(\nabla_{q_e }-\bar{\nabla}_{q_e })
+V(\hat{q})
\right]|\psi>=E|\psi>.
\label{eq:h}\end{equation}
This is the stationary equation, the time evolution will be studied in
section \ref{sec:time}.

\section{Uncertainty relation}
\label{sec:uncertainty}

We have shown \cite{LiCh1} that
the uncertainty relation for electric charge and electric current,
becomes
\begin{equation}
\Delta\hat{q}\cdot\Delta\hat{P} \geq \frac{\hbar}{ 2 }
( 1 + \frac{q^2_e }{\hbar^2} < \hat{H}_0 > ).
\label{eq:de}
\end{equation}
Obviously, this uncertainty relation recovers the usual Heisenberg
uncertainty relation if $ q_e $ goes to zero, i.e. the case that the
discreteness of electric charge vanishes. Moreover, the uncertainty relation
(\ref{eq:de}) has shown us further knowledge than
the traditional
Heisenberg uncertainty relation.

\section{The quantum LC-design}\label{sec:LC}

We only consider the adiabatic approximation here so that
$\varepsilon (t)$ is considered as a constant $\varepsilon$,
{\it i.e.}, $V(\hat{q})=C^{-1}\hat{q}^2/2 + \varepsilon\hat{q}$.
Let us consider a representation in which the operator $\hat{p}$ is
diagonal and called it as p-representation.
The transformation of wave functions between charge representation
and p-representation is given by
\begin{equation}
< n | \psi > = (\frac{q_e}{2\pi\hbar} )
\int^{\hbar(\frac{\pi}{q_e}) }_{-\hbar(\frac{\pi}{q_e} ) }
dp <p |\psi >
e^{-inq_e p /\hbar}
\label{eq:pton}\end{equation}

In the `p-representation', the finite-difference  Schr\"{o}dinger
equation becomes a differential equation for
$\tilde{\psi}(p):=<p|\psi>$,
\begin{equation}
\left[
- \frac{\hbar^2 }{2C} \frac{\partial^2 }{\partial p^2 }
- \frac{\hbar^2}{q_{e}^2 L }( \cos(\frac{q_e}{\hbar} p ) - 1 )
\right] \tilde{\psi } (p) = E\tilde{\psi } (p).
\label{eq:l}\end{equation}
which is the well known Mathieu equation \cite{Wang,Grad}. This equation
was ever appeared in \cite{Jurk} on the discussion of Pad\'{e} approximates.

In terms of the conventional notations \cite{Wang,Grad}, the wave
functions in p-representation can be solved as follows
$$
\tilde{\psi}^{+}_{l}(p) =
{\rm ce }_l (\frac{\pi}{2}- \frac{q_e }{ 2\hbar}p , \,\xi )
$$
or
\begin{equation}
\tilde{\psi}^{-}_{l+1}(p) =
{\rm se}_{l+1}(\frac{\pi}{2}-\frac{q_e}{2\hbar}p,\,\xi)
\label{eq:m}
\end{equation}
where the superscripts `+' and `-' specify  the even and odd parity
solutions respectively;
$l = 0, 1, 2, \cdots $;
$\xi = (2\hbar/q^{2}_e )^2
C/L$; ce$(z, \xi)$ and se$(z, \xi)$ are periodic Mathieu
functions.  In this case, there exist infinitely many eigenvalues
$\{ a_l \}$ and $\{ b_{l+1} \} $ which are not identically equal to zero.
Then the energy spectrum is expressed in terms of the
eigenvalues $a_l $, $b_{l}$ of Mathieu equation

\begin{eqnarray}
E^{+}_l =  \frac{q_{e}^2 }{8C}a_l(\xi)
+ \frac{\hbar^2 }{ q_{e}^2 L},
\nonumber \\
E^{-}_{l+1} =  \frac{q_{e}^2}{8C}b_{l+1}(\xi)
+ \frac{\hbar^2 }{ q_{e}^2 L}.
\label{eq:n}\end{eqnarray}

\section{Quantum L-design and persistent currents }
\label{sec:L}

Introducing an operator
$
\hat{G} :=\exp(-i\beta\hat{q}/\hbar)
$,
we can find that
$\hat{G} | p > = | p - \beta >$
and
$\hat{G}^+ | p > = | p + \beta >.$
For a unitary transformation to the eigenstates of
Schr\"{o}dinger operator given by
\[
| \psi > \rightarrow  | \psi' > = \hat{G} | \psi >,
\]
we find that the Schr\"{o}dinger equation (\ref{eq:h}) is not covariant.
This requests that we introduce a gauge field and define a reasonable
covariant discrete derivative. By making the definitions:
\begin{eqnarray}
D_{q_e} := e^{-i(q_e/\hbar)\phi }
     \frac{ \hat{Q} - e^{ i(q_e/\hbar)\phi } }
          { q_e } , \nonumber \\
\bar{D}_{q_e} := e^{i(q_e/\hbar)\phi }
\frac{ e^{-i(q_e/\hbar)\phi} - \hat{Q^+}}
     { q_e } ,
\label{eq:eg}
\end{eqnarray}
we can verify that they are covariant under a gauge transformation. The
gauge transformations are expressed as
\begin{eqnarray}
\hat{G} D_{q_e} \hat{G}^{-1} = D'_{q_e}, \nonumber \\
\hat{G} \bar{D}_{q_e} \hat{G}^{-1} = \bar{D}'_{q_e},
\label{eq:eh}
\end{eqnarray}
as long as the gauge field $\phi$ transforms in the following way
\[
\phi\rightarrow\phi'=\phi-\beta.
\]
Either the transformation law or the dimension of the field $\phi$
indicates that $\phi$ plays the role of the magnetic flux threading
the circuit.
Then the Schr\"{o}dinger equation for
a pure L-design in the presence of magnetic flux is given by,
\begin{equation}
-\frac{\hbar^2 }{2q_e L } (D_{q_e}\, - \bar{D}_{q_e} )
| \psi > = E | \psi >.
\label{eq:gh}
\end{equation}
The energy spectrum is easily solved
\begin{equation}
E(p,\phi) = \frac{2\hbar}{q^2_e} \sin^2
\left[\frac{q_e}{2\hbar}(p-\phi)
\right],
\end{equation}
which has oscillatory property with respect to
$\phi$ or $p$. Differing from the usual classical pure L-design, the
energy of a mesoscopic quantum pure L-design can not be large than
$2\hbar / q_e^2 $.
Clearly, the lowest energy states are such states that
$p = \phi + n h/q_e$.
Thus the eigenvalues of the electric current
of ground state are calculated
\begin{equation}
I(\phi) = \frac{\hbar}{q_e L} \sin( \frac{q_e}{\hbar} \phi ).
\label{eq:pc}
\end{equation}
Obviously, the electric current on a mesoscopic  circuit of pure L-design
is not null in the presence of a magnetic flux except
$\phi = n(h/q_e)$.
Clearly, this is a pure quantum characteristic.
(\ref{eq:pc}) exhibits that the persistent current in a mesoscopic L-design
is an observable quantity periodically depending on the flux $\phi$.
Because a mesoscopic metal ring is a natural pure L-design, the formula
(\ref{eq:pc}) is valid for persistent current on a single mesoscopic
ring\cite{Chand}. One can easily calculate the inductance of mesoscopic
metal ring and obtain
the formula for persistent currents
\begin{equation}
I(\phi) = \frac{\hbar}{8\pi r
   (\frac{1}{2}\ln(8r/a) - 1)q_e}
     \sin(\frac{q_e}{\hbar}\phi),
\end{equation}
where $r$ is the radius of the ring and $a$ is the radius of the metal wire.
Differing from the conventional formulation of the persistent
current on the basis of quantum dynamics for electrons, our formulation
presented a  method from a new point of view. Formally, the $I(\phi)$
we obtained is a sine function with periodicity of
$\phi_0 =h/q_e$,
But either the model that the electrons move freely
in an ideal ring\cite{Ch1}, or the model that the electrons have hard-core
interactions between them\cite{LiMa2} can only give the sawtooth-type
periodicity.

\section{Quantum C-design and Coulomb blockade}\label{sec:C}

We observe the Schr\"{o}dinger equation for a LC-design.
The mesoscopic capacity may be relatively very small
(about $10^{-8}F$) but the inductance of a macroscopic circuit
connecting to an adiabatic voltage source is relatively large
because the inductance of a circuit is proportional to the
area which the circuit span. Thus we may neglect the term
reversely proportional  to $L$ in the Schr\"{o}dinger equation
and study the equation for a pure C-design:
\begin{equation}
\left(\frac{1}{2C}\hat{q}^2 -\varepsilon\hat{q}
  \right)|\psi>=E|\psi>.
\end{equation}
Apparently, the Hamiltonian operator commutes with
the charge operator,
so they have simultaneous eigenstates.
The energy for the eigenstate $|n>$ is
\begin{equation}
E=\frac{1}{2C}(n q_e - C\varepsilon)^2
  -\frac{C}{2}\varepsilon^2,
\end{equation}
which involves both the charge quantum number and the voltage source.
After some analysis, we can find the relations between charge $q$
and the voltage $\varepsilon$ for the ground state,
\begin{eqnarray}
q=\sum_{m=0}^{\infty}
   \left\{ \theta[\varepsilon -(m+\frac{1}{2})\frac{q_e}{C}]\right.
              \nonumber\\
            \left. -\theta[-\varepsilon -(m+\frac{1}{2})\frac{q_e}{C}]
              \right\}q_e.
\label{eq:chargestep}
\end{eqnarray}
where $\theta(x)$ denotes the step function.
The corresponding eigenstates are
\begin{eqnarray}
|\psi(\varepsilon)>_{ground}=
   \sum_{m=-\infty}^{\infty}
     \left\{\theta[\varepsilon -(m-\frac{1}{2})\frac{q_e}{C}]
            \right. \nonumber \\
            -\left.\theta[-\varepsilon -(m+\frac{1}{2})\frac{q_e}{C}]
             \right\}|m>
\end{eqnarray}
The dependence of the currents on the time is obtained by taking
derivative of (\ref{eq:chargestep}),
\begin{eqnarray}
\frac{dq}{dt}=\sum_{m=0}^{\infty}q_e
   \left\{\delta[\varepsilon -(m+\frac{1}{2})\frac{q_e}{C}]
             \right.    \nonumber\\
            \left. +\delta[\varepsilon + (m+\frac{1}{2})\frac{q_e}{C}]
              \right\}\frac{d\varepsilon}{dt}.
\end{eqnarray}
Clearly, the currents are of the form of sharp pulses which occurs
periodically according to the changes of voltage. The voltage
difference between two pulses are $q_e/C$. This is the called Coulomb
blockade phenomena caused by the charge discreteness.

\section{Time evolution
with time dependent source}\label{sec:time}

We consider a L-design in the presence of time dependent source. The time
dependent Schr\"odinger equation for the circuit in charge representation is
given by
\begin{eqnarray}
i\hbar\frac{d}{dt}|\psi(t)>=\hspace{38mm}\nonumber\\
\left[
   -\frac{\hbar^2}{2q_eL}(\nabla_{q_e}-\bar\nabla_{q_e})
     +\varepsilon(t)\hat{q} \right]|\psi(t)>.
\label{eq:timeSch}
\end{eqnarray}
Expanding the state vector in terms of the orthonormal set of charge eigenstates
\begin{equation}
|\psi(t)>=\sum_{n=-\infty}^{\infty}u_n(t)|n>,
\end{equation}
we obtain the evolution equation for the amplitudes $u_n(t)$, namely
\begin{eqnarray}
i\frac{d}{dt}u_n(t)=[n\frac{q_e}{\hbar}\epsilon(t) +\frac{\hbar}{q_eL}]
   u_n(t)  \nonumber\\
-\frac{\hbar}{2q_eL}(u_{n+1}(t)+u_{n-1}).
\label{eq:timeevo}
\end{eqnarray}
Obviously, $du_n/dt=-i\omega u_n$ corresponds to the stationary case which
was discussed in section \ref{sec:L}.

Multiplying eq.(\ref{eq:pton}) with $\exp(inq_ep'/\hbar)$ and
summing over the index $n$, we obtain
\begin{equation}
<p|\psi>=\sum_{n=-\infty}^\infty <n|\psi>e^{inq_e p/\hbar}.
\end{equation}
Then eq. (\ref{eq:timeSch}) can be written in a convenient form
in $p$-representation, namely,
\begin{equation}
\left\{-\frac{\hbar}{q_eL}[\cos(\frac{q_e}{\hbar}p)-1]
 -i\frac{\hbar}{q_e}\varepsilon(t)\frac{\partial}{\partial p}\right\}
   \tilde\psi
   =i\frac{d}{dt}\tilde\psi,
\label{eq:partialeq}
\end{equation}
where
\begin{equation}
\tilde\psi(p)=<p|\psi>
=\sum_{n=-\infty}^\infty u_n(t)e^{inq_ep/\hbar}.
\end{equation}
The partial differential equation (\ref{eq:partialeq}) can be reduced to a first
order ordinary differential equation by means of the method of characteristics
\cite{character}. Although this reduction involves complicated mathematical
formulations, the solution of the above equation can be given in terms of
Bessel functions\cite{Dunlap}.
Using Graf's addition theorem for Bessel functions
\cite{Grad}, one is able to obtain
\begin{equation}
|u_n(t)|^2 =J^2_n(-\frac{\hbar}{q_e L}\sqrt{a^2(t)+b^2(t)}\,),
\end{equation}
where
\begin{eqnarray*}
a(t)&=&\int^{t}_0 \cos[f(t')]dt', \\
b(t)&=&-\int^{t}_0 \sin[f(t')]dt' \\
f(t)&=&\int^{t}_0 \varepsilon(t')dt'.
\end{eqnarray*}
The mean-square charge is obtained with the help of
the identity $\sum_n J_n^2(z) = z^2/2$,
\begin{equation}
<\hat{q}^2>=\frac{\hbar^2}{2L}[a^2(t)+b^2(t)].
\end{equation}
Some interesting physics phenomena, such as dynamic
localization etc are expected to be studied furthermore.

\section{About the dissipatives}
\label{sec:dissipative}

In order to take account of the effects of dissipatives in the circuit, we
introduce the density matrix $\rho_{m,n}(t)=u^*_m(t) u_n(t)$.
The evolution equation for the density matrix is obtained from
the equation (\ref{eq:timeevo}) for the amplitudes $u_n(t)$,
\begin{eqnarray}
i\frac{\partial}{\partial t}\rho_{m,n} =
  -\varepsilon(t)\frac{q_e}{\hbar}(m-n)\rho_{m,n}\hspace{12mm}\nonumber\\
   +\frac{\hbar}{2q^2_eL}(\rho_{m+1,n}+\rho_{m-1,n}
        -\rho_{m,n+1}-\rho_{m,n-1}).
\end{eqnarray}
Putting a term describing some kind of dissipatives in the circuit, we should
write the evolution equation for the density matrix as
\begin{eqnarray}
i\frac{\partial}{\partial t}\rho_{m,n}=
  -\varepsilon(t)\frac{q_e}{\hbar}(m-n)\rho_{m,n}
                    \nonumber\\
   +\frac{\hbar}{2q^2_eL}(\rho_{m+1,n}+\rho_{m-1,n}
        -\rho_{m,n+1}-\rho_{m,n-1})    \nonumber\\
   +i\sum_{m',n'}\gamma^{m' n'}_{m n}\rho_{m' n'},\hspace{1cm}
\label{eq:SLiouville}
\end{eqnarray}
which is the stochastic Liouville equation. Eq.(\ref{eq:SLiouville})
is a good start point to study the dissipative effects in the mesoscopic
L-design. The one-parameter off-diagonal dissipative is particularly
interesting, {\it i.e.},
$
\gamma_{m' n'}^{m n}=\gamma (1-\delta_{m,n})\delta_{m,m'}\delta_{n,n'}
$
Further discussions are in progress.

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

In the above, we studied the quantization of  mesoscopic electric
circuit. Differing from the literature in which it is simply treated
as the quantization of a harmonic oscillator, we addressed the importance
of the discreteness of electric charge. Taking the discreteness into
account, we proposed a quantum theory for mesoscopic electric circuit
and give a finite-difference Schr\"{o}dinger equation for mesoscopic
electric circuit.
We used the charge representation and
p-representation. Due to the discreteness of electric
charge, $\hat{p}$ is no longer a current operator.

As the Schr\"{o}dinger equation for LC-design in
p-representation becomes the well known Mathieu
equation, it is solvable. We obtain the wave functions
in terms of Mathieu functions and the energy
spectrum in terms of the eigenvalues of Mathieu equation.
The discussion on uncertainty relation for the charge and current shed some
new light on the knowledge of transitional Heisenberg uncertainty relation.
As further applications  of our theory,
Introducing a gauge field and gauge transformation,
we obtained a formula for the persistent current on
the mesoscopic pure L-design in the presence of the magnetic flux.
As the mesoscopic metal ring is a natural pure L-design,
the formula is certainly valid for the persistent current on
mesoscopic rings.
In our formula, the mass of electrons, the carriers for electric current,
is not involved. It is worthwhile to check that property  by experiments.
The  theory is applied to explain the Coulomb blockade.
By considering a quantum C-design, the Coulomb blockade phenomenon was
formulated as an immediate
consequence. The time evolution of quantum L-design is solved formally,
and the dissipative is introduced with the help of density matrix formulation.

The work is supported by NSFC No. 19675030 and NSFZ No.198024.
\baselineskip=10pt
\parskip=0pt


\section{Introduction}

The cross section of inclusive $J/\psi$ hadroproduction measured in
$p\bar p$ collisions at the Fermilab Tevatron \cite{abe,cdf} turned out to be
more than one order of magnitude in excess of what used to be the best
theoretical prediction, based on the colour-singlet model (CSM).
As a solution to this puzzle, Bodwin, Braaten, and Lepage \cite{bod} proposed
the existence of so-called colour-octet processes to fill the gap.
The central idea is that $c\bar c$ pairs are produced at short distances in
colour-octet states and subsequently evolve into physical (colour-singlet)
charmonia by the nonperturbative emission of soft gluons.
The underlying theoretical framework is provided by nonrelativistic QCD
(NRQCD) endowed with a particular factorization hypothesis, which implies a
separation of short-distance coefficients, which are amenable to perturbative
QCD, from long-distance matrix elements, which must be extracted from
experiment.
This formalism involves a double expansion in the strong coupling constant
$\alpha_s$ and the relative velocity $v$ of the bound charm quarks, and takes
the complete structure of the charmonium Fock space into account.

In the case of inelastic $J/\psi$ photoproduction, NRQCD with
colour-octet matrix elements tuned \cite{cho} to fit the Tevatron data
\cite{abe} predicts \cite{cac} at leading order (LO) a distinct rise in
cross section as $z\to1$, where $z$ is the fraction of the photon energy
transferred to the $J/\psi$ meson in the proton rest frame, which is not
observed by the H1 \cite{aid} and ZEUS \cite{bre} collaborations at HERA.
This colour-octet charmonium anomaly has cast doubts on the validity of the
NRQCD factorization hypothesis \cite{bod}, which seems so indispensible to
interpret the Tevatron data in a meaningful way.

Here, we report on an attempt \cite{kni} to rescue the NRQCD approach by
approximately taking into account dominant higher-order (HO) QCD effects.
The basic idea is as follows.
The predicted excess over the HERA data at $z$ close to unity is chiefly
generated by colour-octet $c\bar c$ pairs in the states ${}^1\!S_0$,
${}^3\!P_0$, and ${}^3\!P_2$ \cite{cac}, where we use the spectroscopic
notation ${}^{2S+1}\!L_J$ to indicate the spin $S$, the orbital angular
momentum $L$, and the total angular momentum $J$.
On the other hand, in hadroproduction at the Tevatron, the contributions from
the colour-octet ${}^1\!S_0$ and ${}^3\!P_J$ states fall off much more
strongly with increasing transverse momentum ($p_T$) than the one due to the
colour-octet ${}^3\!S_1$ state \cite{cho}, which is greatly suppressed in the
quasi-elastic limit of photoproduction \cite{cac}.
Consequently, the nonperturbative matrix elements which are responsible for
the colour-octet charmonium crisis are essentially fixed by the Tevatron data
in the low-$p_T$ regime.
This is precisely where the LO approximation used in Ref.~\cite{cho} is
expected to become unreliable due to multiple-gluon radiation from the
initial and final states.
In Ref.~\cite{can}, this phenomenon was carefully analyzed in a Monte Carlo
framework and found to significantly increase the LO cross section.
In Ref.~\cite{kni}, fits to the latest prompt $J/\psi$ data taken by the
CDF collaboration \cite{cdf} at the Tevatron were performed incorporating this
information \cite{can} on the dominant HO QCD effects.
The resulting HO-improved NRQCD predictions for inelastic $J/\psi$
photoproduction at HERA do not overshoot the H1 \cite{aid} and ZEUS \cite{bre}
data any more.

\section{Theoretical input}

The underlying theoretical framework is explained in Ref.~\cite{kni}.
If $p_T$ is of order $M_{J/\psi}$ or below, we adopt the fusion picture, where
the $c\bar c$ bound state is formed within the primary hard-scattering process.
In the high-$p_T$ regime, we work in the fragmentation picture, where the
$c\bar c$ bound state is created from a single high-energy gluon, charm quark
or antiquark which is close to its mass shell.
We take the renormalization scale $\mu$ and the common factorization scale
$M_f$ to be $\mu=M_f=m_T$, where $m_T=\sqrt{4m_c^2+p_T^2}$ is the $J/\psi$
transverse mass.
We define the starting scale $\mu_0$ of the fragmentation functions (FF's) as
$\mu_0=2m_c=M_{J/\psi}$.
For our LO analysis, we choose CTEQ4L \cite{lai} and GRV-LO \cite{glu} as the
proton and photon PDF's, respectively, and evaluate $\alpha_s$ from the
one-loop formula with $\Lambda^{(4)}=236$~MeV \cite{lai}.
Whenever we include higher orders, we adopt the $\overline{\rm MS}$
renormalization and factorization scheme and employ CTEQ4M \cite{lai}, GRV-HO
\cite{glu}, and the two-loop formula for $\alpha_s$ with
$\Lambda_{\overline{\rm MS}}^{(4)}=296$~MeV \cite{lai}.

Unfortunately, not all ingredients which would be necessary for a fully
consistent NLO analysis are yet available.
In the case of fusion, the NLO corrections to the partonic cross sections are
only known for direct photoproduction in the CSM \cite{kra}.
Furthermore, in the case of fragmentation, the NLO corrections to the FF's
at the initial scale $\mu_0$ are still unknown.
In the case of direct $J/\psi$ photoproduction under typical HERA conditions,
the QCD correction factor $K$ to the inclusive cross section in the CSM was
found \cite{kra} to be as low as 1.2 in the inelastic regime $z\lsim0.9$.
It is plausible that the $K$ factors for the color-octet and resolved-photon
processes should be modest, too.
However, the situation should be very different for inclusive $J/\psi$
hadroproduction at the Tevatron, especially in the low-$p_T$ range, where one
expects substantial HO QCD effects due to multiple-gluon radiation.
Such effects were estimated for the fusion mechanism in Ref.~\cite{can} by
means of the Monte Carlo event generator PYTHIA \cite{sjo} after implementing
therein the relevant colour-octet processes, and they were indeed found to be
very sizeable.
The impact of these effects on the fit to the latest CDF data \cite{cdf} is
demonstrated in Table~\ref{t1}.

\begin{table}
\begin{center}
\caption{Values of the $J/\psi$ matrix elements resulting from the
LO and HO-improved fits to the CDF data \protect\cite{cdf}.
Here, $M_r^{J/\psi}=
\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^1\!S_0]\right\rangle
+(r/m_c^2)
\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!P_0]\right\rangle$.
The first 10 (last 2) of the 11 data points are described in the fusion
(fragmentation) picture.}
\label{t1}
\medskip
\begin{tabular}{|c|c|c|}
\hline\hline
 & LO & HO \\
\hline
$\left\langle{\cal O}^{J/\psi}[\,\underline{1},{}^3\!S_1]\right\rangle$ &
$(7.63\pm0.54)\cdot10^{-1}$~GeV$^3$ & $(1.30\pm0.09)$~GeV$^3$ \\
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!S_1]\right\rangle$ &
$(3.94\pm0.63)\cdot10^{-3}$~GeV$^3$ & $(2.73\pm0.45)\cdot10^{-3}$~GeV$^3$ \\
$M_r^{J/\psi}$ &
$(6.52\pm0.67)\cdot10^{-2}$~GeV$^3$ & $(5.72\pm1.84)\cdot10^{-3}$~GeV$^3$ \\
$r$ & 3.47 & 3.54 \\
$\chi_{\rm DF}^2$ & 7.49/12 & 3.96/12 \\
\hline\hline
\end{tabular}
\end{center}
\end{table}

\section{Predictions for charmonium photoproduction}

We now explore the phenomenological consequences of this HO improvement for
inclusive $J/\psi$ photoproduction in $ep$ collisions at HERA, with beam
energies $E_e=27.5$~GeV and $E_p=820$~GeV in the laboratory frame, assuming
the maximum photon virtuality to be $Q_{\rm max}^2=4$~GeV$^2$.
As in the H1 \cite{aid} and ZEUS \cite{bre} publications, we convert the $ep$
cross sections to averaged $\gamma p$ cross sections by dividing out the
photon-flux factor.
The contribution due to $\psi^\prime$ mesons with subsequent decay into
$J/\psi$ mesons is approximately taken into account by multiplying the
theoretical predictions by an overall factor of 1.15.
The data are mostly concentrated in the low-$p_T$ range, where the fusion
picture should be valid.

In Fig.~1, we compare our LO and HO-improved predictions for the
$p_T^2$, $z$, and $y$ distributions with the H1 \cite{aid} and ZEUS \cite{bre}
data.
Here, $y$ is the $J/\psi$ rapidity in the laboratory frame, which is taken to
be positive in the proton flight direction.
The circumstance that
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^1\!S_0]\right\rangle$ and
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!P_0]\right\rangle$
are not separately fixed by the fit to the CDF data \cite{cdf} induces some
uncertainty
in the colour-octet contributions to the cross sections of direct and resolved
photoproduction and thus also in the total cross section.
This uncertainty is encompassed by the results for
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^1\!S_0]\right\rangle=
M_r^{J/\psi}$ and
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!P_0]\right\rangle=0$ and
those for
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^1\!S_0]\right\rangle=0$ and
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!P_0]\right\rangle
=(m_c^2/r)M_r^{J/\psi}$, which are actually shown in Fig.~1.
We observe that, at LO, the colour-octet contribution of direct
photoproduction is dominant for $z\gsim0.5$.
Thus, it also makes up the bulk of the $p_T^2$ and $y$ distributions,
which are integrated over $0.4<z<0.9$.
This contribution is responsible for the significant excess of the LO
predictions over the experimental results for $d\sigma/dp_T^2$ and
$d\sigma/dz$ at low $p_T$ and high $z$, respectively.
On the other hand, the HO-improved predictions tend to undershoot the data
leaving room for a substantial $K$ factor due to the missing NLO corrections
to the partonic cross sections.
Now, the colour-singlet contribution of direct photoproduction, which is well
under theoretical control \cite{kra}, is by far dominant, except in the
corners of phase space, at $z\lsim0.15$ and $z\gsim0.85$, where the colour-octet
contributions of resolved and direct photoproduction, respectively, take over.
Of course, we should also bear in mind that the predictions shown in
Fig.~1 still suffer from considerable theoretical
uncertainties related to the choice of the scales $\mu$ and $M_f$, the PDF's,
and other input parameters such as $m_c$ and $\Lambda^{(4)}$ \cite{kra}.
From these observations, we conclude that it is premature at this point to
speak about a discrepancy between the Tevatron \cite{abe,cdf} and HERA
\cite{aid,bre} data of inclusive $J/\psi$ production within the framework of
NRQCD \cite{bod}.

\begin{figure}[ht]
\label{f1}
\begin{tabular}{ll}
\parbox{7.5cm}{
\epsfig{file=ptjlo.ps,width=7.5cm}
}
&
\parbox{7.5cm}{
\epsfig{file=ptjnlo.ps,width=7.5cm}
}
\\
\vspace*{1mm}
\parbox{7.5cm}{
\epsfig{file=zjlo.ps,width=7.5cm}
}
&
\parbox{7.5cm}{
\epsfig{file=zjnlo.ps,width=7.5cm}
}
\\
\vspace*{1mm}
\parbox{7.5cm}{
\epsfig{file=yjlo.ps,width=7.5cm}
}
&
\parbox{7.5cm}{
\epsfig{file=yjnlo.ps,width=7.5cm}
}
\end{tabular}
\vspace*{2mm}
\caption[junk]{{\it
The H1 \protect\cite{aid} and ZEUS \protect\cite{bre} data on inelastic
$J/\psi$ photoproduction are compared with the LO and HO-improved predictions.
The total results (solid lines) are built up by the
direct-photon colour-singlet (dotted lines),
direct-photon colour-octet (dashed lines),
resolved-photon colour-singlet (dot-dashed lines), and
resolved-photon colour-octet (dot-dot-dashed lines) contributions.
}}
\end{figure}

\section{Conclusions}

We determined the $J/\psi$ colour-octet matrix elements which appear in the
NRQCD expansion \cite{bod} at leading order in $v$ by fitting the latest
Tevatron data of prompt $J/\psi$ hadroproduction \cite{cdf}.
We found that the result for the linear combination $M_r^{J/\psi}$ of
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^1\!S_0]\right\rangle$ and
$\left\langle{\cal O}^{J/\psi}[\,\underline{8},{}^3\!P_0]\right\rangle$ is
substantially reduced if the HO QCD effects due to the multiple emission of
gluons, which had been estimated by Monte Carlo techniques \cite{can}, are
taken into account.
As an important consequence, the intriguing excess of the LO NRQCD prediction
for inelastic $J/\psi$ photoproduction at $z$ close to unity \cite{cac}
over the HERA measurements \cite{aid,bre} disappears.
We assess this finding as an indication that it is premature to proclaim an
experimental falsification of the NRQCD framework on the basis of the HERA
data.
Although we believe that our analysis captures the main trend of the HO
improvement, we stress that it is still at an exploratory level, since a
number of ingredients which would be necessary for a fully consistent NLO
treatment of inclusive $J/\psi$ hadroproduction and photoproduction are still
missing.

\medskip

\noindent
{\it Achnowledgements.} The author thanks Gustav Kramer for his collaboration
on this work.


\section{Introduction}

Young stellar objects (YSOs) are divided into 4 classes by the spectral-energy
 distribution at infrared to millimeter wave lengths (Andr\'{e} \& Montmerle 1994):
Class 0 and  Class I protostars,
Class II sources (classical  T Tauri stars, CTTS) and
Class III sources (weak-line T Tauri stars, WTTS).
The Class 0 protostars have massive envelopes which collapse towards the central region.
Their age  is  estimated to be $\sim$10$^4$ yr.
The class I protostars have circumsteller disks and dilute dust envelopes.
Their age is estimated to be $\sim$10$^5$ yr.
The Einstein and ROSAT observatories  detected soft X-ray emission from both CTTS and WTTS,
 whose luminosities range from 10$^{28}$ erg s$^{-1}$ to 10$^{30}$ erg s$^{-1}$
(Montmerle et al. 1993, Neuh\"{a}user 1997).
Recently, ROSAT and ASCA discovered X-ray
 emission from Class I protostars (Casanova et al. 1995, Koyama et al. 1996, Grosso et al. 1997, Kamata et al. 1997).
Since X-ray emission has been detected from only a dozen protostars (Carkner et al. 1998),
 it is important to increase the sample size
of X-ray emitting protostars to investigate the origin of their X-ray emission.

The HH24-26 region is located in the L1630 Orion molecular cloud complex at
 a distance of about 460 pc (Chini et al. 1993).
This region is  an active star-forming region containing several YSOs,
  Class 0 sources (HH24MMS, HH25MMS),  Class I sources (SSV63, SSV59) and
 a Class II source (SSV61) (Bontemps et al. 1995).
SSV63 and SSV61 were detected in the 2.2 $\mu$ mapping survey (Strom et al. 1976)
 and later SSV63 was resolved into four near-infrared sources, SSV63E, SSV63W, SSV63NE-1, and SSV63NE-2
 (Moneti \& Reipurth 1995).
Both SSV63E and SSV63W show flat or rising spectra in the near- and mid-infrared
band (Zealy et al. 1992),  and hence belong to Class I.
SSV63NE-1 and SSV63NE-2 are unclassified diffuse sources,
and may be reflection nebulae (Moneti \& Reipurth 1995).
A complex of Herbig-Haro jets (HH24)  was detected around SSV63
 by optical observations  and some of the jets are probably
 associated with SSV63 itself (Jones et al. 1987, Mundt et al. 1991).
In this paper, we report the first detection of X-ray emission
 from the infrared protostar SSV63.

 \section{Observations}

We observed the HH24-26 region  with the ASCA observatory
(Tanaka et al. 1994) for a net exposure time of about 30 ks on 1998 October 2.
The center  coordinates of the pointing position were  R.A.(J2000)$=5^{\rm{h}} 46^{\rm{m}} 27.6^{\rm{s}}$ and Decl.(J2000)$=-0^{\circ} 7' 58.7''$.
ASCA has four identical X-ray telescopes (XRT), which achieve
 large effective area at energies up to 10 keV (Serlemitsos et al. 1995).
Two Solid state Imaging Spectrometers (SIS0, SIS1) and
 two Gas Imaging Spectrometers (GIS2, GIS3) are located at each focus
 of an XRT.
Each SIS has four CCD chips, which have an 11' $\times$ 11' field of view each.
These instruments cover the energy range from 0.4 keV to 10 keV (Burke et al. 1991)
 with an energy resolution of  $\sim$160 eV (FWHM) at 6 keV
 and $\sim$100 eV (FWHM) at 1.5 keV at the time of  the present observation.
The GIS  has a  field of view of 40' diameter and covers
 the energy range from 0.7 keV to 10 keV
with an  energy resolution of 0.5 keV (FWHM) at 5.9 keV
(Ohashi et al. 1996).
In this observation, the SIS was operated in the 2-CCD mode, while the GIS
 was in the normal pulse-height  mode with nominal bit assignments.

\section{Results}
\subsection{X-ray image of HH24-26 region}

Figure \ref{image} shows the SIS  soft band (0.7$-$2.0 keV) and hard band (2.0$-$8.0 keV) images
 of the HH24-26 region averaged over the whole observation.
In the figure, the SIS0 and SIS1 images are superimposed
and smoothed by a gaussian function with $\sigma$ of 0.2'.
The symbols with numbers from 1 to 8 in Figure \ref{image} correspond to the sources cited in the caption.
Clearly, we detected two X-ray sources.
The southern source is conspicuous only in the soft band (0.7$-$2.0 keV) image, whereas
the northern source is predominant in the hard band (2.0$-$8.0 keV) image.
Both detections are consistent with point sources within the angular resolution
 of the XRT (half power diameter about 3').
In order to derive accurate positions and count rates of the detected sources,
we performed 2-dimensional fitting to the SIS image using the point spread function of the XRT (Ueda 1996).
The positions determined with ASCA are ($5^{\rm{h}}46^{\rm{m}}7.0^{\rm{s}}$, $-0^{\circ}10'7''$) and ($5^{\rm{h}}46^{\rm{m}}7.3^{\rm{s}}$, $-0^{\circ}12'4''$), with 1 $\sigma$ error of
14''.
Based on the positional coincidence, the southern source is identified with
 the T Tauri star  SSV61.
The northern  source is identified  with the protostars  SSV63E/SSV63W,
 but  the ASCA image resolution does not allow us to resolve X-ray emission
 of the two (hereafter we designate the emission as from SSV63E+W).
The count rates in the 0.7$-$7.0 keV band of SSV61 and SSV63E+W
are 1.4 (1.3$-$1.6) $\times 10^{-2}$ counts s$^{-1}$
(the errors in parentheses indicate a 90\% confidence region)
and 8.2 (7.0$-$9.4) $\times 10^{-3}$ counts s$^{-1}$
 , respectively.
We were not able to detect any evidence of X-ray emission from the Class I protostar SSV59 or
the Class 0 protostars, HH24MMS and HH25MMS.
We derived upper limits (2 $\sigma$) of count rates of these protostars,  2 $\times 10^{-4}$ counts s$^{-1}$, 8 $\times 10^{-4}$ counts s$^{-1}$, and 5 $\times 10^{-4}$ counts s$^{-1}$, respectively for SSV59, HH24MMS, and HH25MMS.

\subsection{Light curve}

Figure 2 shows background-subtracted light curves in the hard (3.0$-$8.0 keV)
 and the soft (0.7$-$2 keV) X-ray bands.
In the figure, all the data from the  four sensors (SIS0, SIS1, GIS2, GIS3)
are combined to obtain the best statistics.
The data were taken from regions of 6' radius for the GIS and 4' radius
 for the SIS, centered at the position of SSV63E+W.
The circle shown in Figure \ref{image} indicates the region where the SIS light curve was derived.
We extracted the background which is constructed
from the whole FOV of each sensor where point sources were excluded, correcting for the positional dependence of the background intensity (Ueda 1996).
Although  the region contains both SSV61 and SSV63E+W,
the hard band light curve should correspond to that of  SSV63E+W
 and the soft band light curve to that of SSV61,
as expected from the image analysis in section 3.1.
We discarded the data in the 2.0$-$3.0 keV band from the hard band light curve
 because X-rays both from SSV63E+W and  SSV61 contribute comparably in the 2.0$-$3.0 keV band,
 as shown in Figure \ref{spec} of section 3.3.

As can be seen from the top panel of Figure \ref{lc}, SSV63E+W showed
 flare-like time variability.
On the other hand, the soft band light curve (bottom panel of Figure \ref{lc})
 showed little intensity  variability,
 suggesting that the T Tauri star SSV61 was relatively stable during the ASCA observation.
The flare from SSV63E+W is characterized by a slow rise followed by an exponential decay.
We fitted the light curve with a model consisting of a linear rise and an exponential decay
to a  constant level of the form,
\begin{equation}
F(t)=\cases{\frac{N}{\tau_{r}}(t-t_{p}+\tau_r)+C,     &  $(t < t_{p})$, \cr
N exp\{\frac{-(t-t_{p})}{\tau_{d}}\}+C,               &  $(t \ge t_{p}).$       }
\end{equation}
This model has five free parameters, flare peak time $t_{p}$,
rise time $\tau_{r}$,
a constant level $C$,
the normalization of the exponential function $N$,
and e-folding decay time $\tau_{d}$.
The solid curve in Figure 2 shows the best-fit model.
From the fit, we derived a peak flux ($N+C$) of 0.18 (0.16$-$0.19) cts s$^{-1}$,
$C$ = 0.02 (0.016$-$0.026)  cts s$^{-1}$,
$\tau_{r}$ = 1.7 (1.2$-$2.2) $\times 10^{4}$ s,
and $\tau_d$ = 1.2 (0.9$-$1.5) $\times 10^{4}$ s.

\subsection{Energy spectrum}

The upper and lower panels in Figure \ref{spec} show background-subtracted spectra
 taken from the  interval A (hereafter termed the ``flare state'')
 and  the interval B (the ``quiescent state'') indicated in Figure \ref{lc}, respectively.
Since the separation between SSV63E+W and SSV61 ($\sim$2') is not
large enough compared with the PSF of  the XRT of ASCA,
it is difficult to obtain spectra from SSV63E+W and SSV61 without contamination from the other
 even if we extract the spectra in a small region around each source.
Hence, we accumulated photons in the same region as used in section 3.2 to obtain better statistics.
The background spectrum is also extracted  from the same region as used in section 3.2.
Since the region includes both the positions of SSV63E+W and SSV61,
the spectra will contain contributions from both the sources.
Both of the spectra in Figure \ref{spec} show
a double peak feature, suggesting the existence of two (hard and soft)
components. Combined with the results from the image analysis in
section 3.1, we assume that the soft and the hard components
correspond to SSV61 and SSV63E+W, respectively. In fact, the soft
component does not show a significant change between the two states,
while the hard component is enhanced during the flare state. This confirms
that the hard component in the spectrum really corresponds to the
emission from the Class I protostar SSV63E+W.

In the hard component during the flare state, an emission line feature is
clearly seen at 6.6 (6.5$-$6.7) keV, which is consistent with the
K$\alpha$ line from He-like iron ions. This suggests that the hard
component during the flare state originates from an optically thin thermal
plasma.  On the other hand, the origin of the soft component from
SSV61 is also considered to be thin thermal coronal emission, which is
typical for T Tauri stars (Montmerle et al. 1993, Neuh\"auther 1997).
Accordingly, we fit the combined GIS and SIS spectra with a model
consisting of two thermal components with different absorption column
densities, in which the high and low temperature components represent
the contribution from SSV63E+W and from SSV61, respectively.
This model of two thermal components was also adopted for the quiescent spectrum,
utilizing the thin thermal spectral
model calculated by Raymond and Smith (1977; hereafter the ``RS model'').
Since the statistics of the spectrum are not good enough to determine
 the abundance of the soft component, we fixed this  at 0.5 solar  for SSV61 because
 the metal abundances of late-type  stars are roughly 0.5 solar (Carkner et al. 1996).
For the hard component, we used an Auxiliary Response Function (ARF) constructed for a point source
located in the center of the region selected (i.e., SSV63E+W), while
for the soft component, we used an ARF constructed for a point source located at the position of SSV61.

Using this ``two-component'' RS model, we obtained
acceptable fits for both of the flare-state and quiescent state spectra.
The temperature and the absorption column density did
not change significantly between the flare state and the quiescent state for
either the hard or soft components. Hence, to make the tightest
constraints on these parameters, we extracted spectra averaged over
the whole observation, and fit them with the same model. We
again obtained acceptable fits with parameters listed in Table \ref{tbl-2}.
We derived plasma temperatures $kT=$ 0.9 (0.7$-$1.2) and 5.0 (3.3$-$7.9) keV,
time-averaged emission measures of 5.3 (3.2$-$9.0) $\times$ 10$^{54}$ cm$^{-3}$
and 14 (9.5$-$25) $\times$ 10$^{54}$ cm$^{-3}$,
and absorption column densities
$N_{\rm H} =$1.3 (0.9$-$1.7) $\times$ 10$^{22}$ cm$^{-2}$ and
 1.5 (1.2$-$1.8) $\times$ 10$^{23}$ cm$^{-2}$, respectively for the low and  high temperature
 components.
We obtained an abundance of 0.45 (0.28$-$0.69) solar for the high temperature component, SSV63E+W.
Finally, to determine the
spectral parameters of SSV63E+W separately during the flare state and the
quiescent state, we repeated the fit by fixing the
parameters of the soft component at the best-fit values obtained from
the whole observation. The results are also listed in Table
\ref{tbl-2}.

\section{Discussion}

Since X-ray emission from Class 0 protostars has not been detected yet,
 it is interesting  whether Class 0 protostars in the HH24-26 region,
  HH24MMS and HH25MMS,  emit X-rays or not.
The present ASCA observation, however, failed to detect X-ray emission from
either HH24MMS or HH25MMS.
We estimate upper limits of X-ray luminosity as
9 $\times$ 10$^{30}$ erg s$^{-1}$ for HH24MMS,
 7 $\times$ 10$^{30}$ erg s$^{-1}$ for HH25MMS, and
 4 $\times$ 10$^{30}$ erg s$^{-1}$ for the Class I protostar SSV59,
assuming a  RS model with  $kT=$ 5 keV, $N_{\rm{H}}=$ 1 $\times$ 10$^{23}$ cm$^{-2}$,
and a metal abundance of 0.5 solar.
We note that these estimates depend on the spectral parameters assumed.

The most remarkable result here is the first detection of X-ray emission from a Class I protostar SSV63E+W
in Orion.
The significant features of  the X-ray emission are :
(1) a large X-ray luminosity of about 1 $\times$ 10$^{32}$ erg s$^{-1}$ in the 0.5$-$10 keV band
and a high temperature of  $kT$ = 5.2 ( $>$ 2.3)  keV during the quiescent state,
(2) an X-ray flare with total energy release of $5 \times 10^{36}$ erg with an elevated temperature
of $kT \sim 6$ keV, and (3) a large absorption column density of  $N_{\rm{H}}=$ 1.1 (1.0$-$1.3) $\times$ 10$^{23}$ cm$^{-2}$.

SSV63E+W has a high temperature plasma even during the quiescent state,
which is comparable to that during the flare state.
Such a high temperature is not observed in TTSs during the quiescent state.
A cluster of Class I protostars in the R CrA molecular cloud also showed
 a high temperature of $kT \sim 7$ keV  during the quiescent state (Koyama et al. 1996).
Thus, some of the Class I sources have high temperature plasmas of $\sim$ 6 keV.
In order to know the mechanism of the quiescent X-ray emission,
we compare the L$_{X}$/L$_{bol}$ ratio of SSV63E+W with those of T Tauri stars.
The bolometric luminosity of SSV63E is estimated to be 22.4 L$_{\odot}$ (Berrilli et al. 1989, Zealey et al. 1992).
Although that of SSV63W is highly uncertain,
assuming  the X-ray emission of SSV63E+W comes from only SSV63E, we obtained L$_{X}$/L$_{bol}$ ratio of
 1.2 $\times$ 10$^{-3}$ in the quiescent state.
Since the typical L$_{X}$/L$_{bol}$ ratio of T Tauri stars is 5 $\times$ 10$^{-4}$  and the maximum is 3 $\times$ 10$^{-3}$ (Preibisch et al. 1998), our value of the L$_{X}$/L$_{bol}$ ratio of SSV63E+W is comparable to that of T Tauri stars.
Thus the X-ray emission mechanism from protostars during the quiescent state
has to produce
 higher temperature plasmas of $\sim$ 6 keV, but not larger  L$_{X}$/L$_{bol}$ ratios than those of T Tauri stars.
The enhanced disk-magnetosphere interaction (Shu et al. 1997) may provide such high temperature plasmas in protostars.

The hard ($kT \sim 6$ keV) X-ray flare of SSV63E+W can probably be explained by strong magnetic activity,
by analogy with the X-ray flares of TTSs.
Such a hard X-ray flare was also observed  from a protostar
in the R Cr A molecular clouds (Koyama et al. 1996).
Hayashi et al. (1996) suggested a model for hard X-ray flares in
protostars, based on a closed magnetic loop connecting the central
star and disk. This model can explain the high temperature of X-ray emitting protostars during the flare state.
From the decay time analysis (van den Oord \& Mewe 1989, Montmerle 1990) of the observed flare,
we estimated the electron density $n_e$ $\sim 1 \times$ 10 $^{11}$ cm$^{-3}$,
the volume  of the emission region $V \approx E.M./n_e^2$ $\sim 2 \times 10^{33}$ cm$^{3}$,
the typical length of the  region $d = V^{\frac{1}{3}}$ $\sim 1 \times 10^{11}$ cm,
and  the magnetic field strength $B$ $\sim$ 200 gauss for the flare region in SSV63E+W.
These parameters constrain the model of X-ray flare of the protostar.

\acknowledgments
The authors are grateful to Dr. Hugh Hudson and the referee,
Dr. Thierry Montmerle for their stimulating discussions.
H.O. achnowledges the award of Research Fellowships of the Japan Society for the
 Promotion of Science for Young Scientists.

\newpage