math problems

[Mohamed Bana <mbana.lists@googlemail.com>]

[-- Attachment #1: Type: text/plain, Size: 765 bytes --]

hi guys,

i generally use this file to test if math is "working", with the recent
update i get a lot of errors such as;

! Missing number, treated as zero.
<to be read again>
                      $
\@@dobig ...o #1\bodyfontsize {}\right .\n@space $
                                                     }}
\@mt@defaultBigl ...\puremathcomm {open}{\Big {#1}
                                                     }
l.99 \pauli = \dirac^2 = \Bigl(
                                  \sum_{j=1}^2
\sigma_j\big(-i\PD{}{x^j}-a_j\bi...

?


with texlive 2008 it works just fine, i've attached the output of the
TL2008.  i'm certain that the same file was compiling just fine with (i
think) luatex 0.31 (or pre 0.31).

the content is from Mikael Persson's thesis.

Mohamed



[-- Warning: decoded text below may be mangled, UTF-8 assumed --]
[-- Attachment #2: mickep-math.tex --]
[-- Type: text/x-tex; name="mickep-math.tex", Size: 12903 bytes --]

\usemodule[bib]

\def\mathbb#1{{\blackboard #1}}
\def\pauli{\mathfrak{P}}
\def\mathfrak#1{{\fraktur #1}}


% Matriser av olika typ.
% Paranthesis
\definemathmatrix
[pmatrix]
[left={\left(\,},right={\,\right)}]
% Brackets
\definemathmatrix
[bmatrix]
[left={\left[\,},right={\,\right]}]
% Curly braces
\definemathmatrix
[Bmatrix]
[left={\left\{\,},right={\,\right\}}]
% vertical bars
\definemathmatrix
[vmatrix]
[left={\left\vert\,},right={\,\right\vert}]
% double vertical bars
\definemathmatrix
[Vmatrix]
[left={\left\Vert\,},right={\,\right\Vert}]

\def\PD#1#2{\frac{\partial #1}{\partial #2}}
\def\ann{\mathscr{Q}^{\vphantom{*}}}
\def\mathscr#1{{\gothic #1}}
\def\cre{\mathscr{Q}^*}
\def\dirac{{\mathfrak{D}}}
\def\pform{\mathfrak{p}}
\def\ed{\mfunction{\,d}}
\definemathcommand [dom] [nolop] {\mfunction{Dom}}
\def\bz{\bar{z}}
\def\psip{{\psi_+}}
\def\psim{{\psi_-}}

% \eqref
\definereferenceformat[eqref][left=(,right=)]

\defineenumeration[problem]
        [text=Problem,
        location=serried,
        width=fit,
        indenting=first,
        distance=0.5em,
	way=bysection,
        ]

\subsubsection[ksec:paulitvad]{The Pauli operator in $\mathbb{R}^2$}

A charged spin $1/2$ particle is described by the Pauli Hamiltonian, which acts in $L_2(\mathbb{R}^2)\otimes \mathbb{C}^2$, and is formally defined as
\placeformula[keq:pauli]
\startformula
\pauli = 
\startpmatrix
\NC H-\frac{g}{2}B \NC 0 \NR
\NC 0 \NC H+\frac{g}{2}B\NR
\stoppmatrix.
\stopformula
Here $H$ is the two-dimensional Schrödinger Hamiltonian $H=(-i\nabla -\vec{a})^2$, $B$ is the magnetic field (In two dimensions we identify the two-form and the coefficient function), and $g$ is the {\em gyromagnetic ratio}. We identify the real point $(x^1,x^2)$ with the complex number $z=x^1+ix^2$, and denote a scalar potential of $B$ by $W$,
\startformula
-\Delta W = B.
\stopformula
We set $\Pi_j = -i\PD{}{x^j}-a_j$ and
\startformula
 \ann = \Pi_1 - i \Pi_2,\quad \cre = \Pi_1 + i\Pi_2,
\stopformula
and note that
\placeformula[keq:komm]
\startformula
\ann\cre=\cre\ann+2B,\quad H=\cre\ann+B=\ann\cre-B.
\stopformula
From~\eqref[keq:pauli] and~\eqref[keq:komm] we get
\placeformula[keq:paulig]
\startformula
\pauli = 
\startpmatrix
\NC \cre\ann-\frac{g-2}{2}B \NC 0 \NR
\NC 0 \NC \ann\cre+\frac{g-2}{2}B\NR
\stoppmatrix.
\stopformula
The number $\frac{g-2}{2}$ is called the {\em anomaly factor} of the magnetic moment. Experiments give an anomaly factor of $0.00159$ for the electron~\cite[bovo]. We assume that $g=2$, which is the simplest case. Thus, the Pauli Hamiltonian we study in this thesis is formally defined by
\placeformula[keq:pauliw]
\startformula
\pauli = 
\startpmatrix
\NC \cre\ann \NC 0 \NR
\NC 0 \NC \ann\cre\NR
\stoppmatrix.
\stopformula
The Pauli operator can be written as the square of the Dirac operator
\placeformula[keq:paulia]
\startformula
\pauli = \dirac^2 = \Bigl(\sum_{j=1}^2 \sigma_j\big(-i\PD{}{x^j}-a_j\big)\Bigr)^2 = 
\startpmatrix
\NC 0 \NC \cre \NR
\NC \ann \NC 0\NR
\stoppmatrix^2
\stopformula
from which it follows that it is a non-negative operator. Now let us be more precise about the domains. As in the case of the Schrödinger Hamiltonian there is a problem in defining the Pauli Hamiltonian if the magnetic field is too singular. The quadratic form corresponding to $\pauli$ is given by
\placeformula[keq:pform]
\startformula
\pform(\psi,\psi)= \int_{\mathbb{R}^2} \Big|\sum_{j=1}^2 \sigma_j\big(-i\PD{}{x^j}-a_j\big)\psi\Big|^2\ed m(x).
\stopformula

If $\vec{a}\in L_{2,\text{loc}}(\mathbb{R}^2)\otimes \mathbb{R}^2$ then $\pform(\psi,\psi)$ makes sense for $\psi\in C_0^\infty(\mathbb{R}^2)\otimes \mathbb{C}^2$. We define the {\em minimal} Pauli form $\pform_{\text{min}}$ as
\startformula
\startalign
\NC \dom(\pform_{\text{min}}) \NC = C_0^\infty(\mathbb{R}^2)\otimes \mathbb{C}^2;\NR
\NC \pform_{\text{min}}(\psi,\psi) \NC = \pform(\psi,\psi),\quad \psi\in\dom(\pform_{\text{min}}).\NR
\stopalign
\stopformula
It is closable and thus a self-adjoint operator $\pauli_{\text{min}}$ can be defined. We also define the {\em maximal} Pauli form $\pform_{\text{max}}$ as
\placeformula
\startformula
\startalign
\NC \dom(\pform_{\text{max}}) \NC = \bigl\{\,\psi\in L_2(\mathbb{R}^2)\otimes \mathbb{C}^2\bigm| \pform(\psi,\psi)<\infty\,\bigr\};\NR[keq:pformmax]
\NC \pform_{\text{max}}(\psi,\psi) \NC = \pform(\psi,\psi),\quad\psi\in\dom(\pform_{\text{max}}).\NR
\stopalign
\stopformula
In the presence of AB solenoids, $\vec{a}$ does not belong to $L_{2,\text{loc}}(\mathbb{R}^2)\otimes \mathbb{R}^2$. It was proved in~\cite[so] that the Pauli form can not be defined on smooth compactly supported $\psi$ via~\eqref[keq:pform] in this case. The way out of this is to redefine the Pauli form $\pform$ by an expression that makes sense even in this more singular case. This is done in~\cite[ervo] by writing the operators $\ann$ and $\cre$ as
\placeformula[keq:anncre]
\startformula
\ann =-2i e^{W}\PD{}{\bar{z}} e^{-W}\quad\text{and}\quad \cre = -2i e^{-W}\PD{}{z} e^{W},
\stopformula
and noting that the quadratic form
\placeformula[keq:pformw]
\startformula
\pform(\psi,\psi) = 4\int_{\mathbb{R}^2} \Big|\PD{}{\bz}\left(e^{-W}\psip\right)\Big|^2e^{2W}+\Big|\PD{}{z}\left(e^{W}\psim\right)\Big|^2 e^{-2W}\ed m(x),
\stopformula
$\psi=(\psip,\psim)^t$ makes sense even with this more singular field. If $\pform$ is defined on a maximal domain in the same way as in~\eqref[keq:pformmax], it yields a self-adjoint operator even with this singular field, usually called the {\em maximal} Pauli operator. The forms in~\eqref[keq:pform] and~\eqref[keq:pformw] coincides for more regular fields.


\page

\placeformula[keq:pauli]
\startformula
\pauli = 
\startpmatrix
\NC H-\frac{g}{2}B \NC 0 \NR
\NC 0 \NC H+\frac{g}{2}B\NR
\stoppmatrix.
\stopformula

\placeformula
\startformula
-\Delta W = B.
\stopformula
We set $\Pi_j = -i\PD{}{x^j}-a_j$ and
\startformula
 \ann = \Pi_1 - i \Pi_2,\quad \cre = \Pi_1 + i\Pi_2,
\stopformula

\placeformula[keq:komm]
\startformula
\ann\cre=\cre\ann+2B,\quad H=\cre\ann+B=\ann\cre-B.
\stopformula

\placeformula[keq:paulig]
\startformula
\pauli = 
\startpmatrix
\NC \cre\ann-\frac{g-2}{2}B \NC 0 \NR
\NC 0 \NC \ann\cre+\frac{g-2}{2}B\NR
\stoppmatrix.
\stopformula

\placeformula[keq:pauliw]
\startformula
\pauli = 
\startpmatrix
\NC \cre\ann \NC 0 \NR
\NC 0 \NC \ann\cre\NR
\stoppmatrix.
\stopformula

\placeformula[keq:paulia]
\startformula
\pauli = \dirac^2 = \Bigl(\sum_{j=1}^2 \sigma_j\big(-i\PD{}{x^j}-a_j\big)\Bigr)^2 = 
\startpmatrix
\NC 0 \NC \cre \NR
\NC \ann \NC 0\NR
\stoppmatrix^2
\stopformula

\placeformula[keq:pform]
\startformula
\pform(\psi,\psi)= \int_{\mathbb{R}^2} \Big|\sum_{j=1}^2 \sigma_j\big(-i\PD{}{x^j}-a_j\big)\psi\Big|^2\ed m(x).
\stopformula

\startformula
\startalign
\NC \dom(\pform_{\text{min}}) \NC = C_0^\infty(\mathbb{R}^2)\otimes \mathbb{C}^2;\NR
\NC \pform_{\text{min}}(\psi,\psi) \NC = \pform(\psi,\psi),\quad \psi\in\dom(\pform_{\text{min}}).\NR
\stopalign
\stopformula

\placeformula
\startformula
\startalign
\NC \dom(\pform_{\text{max}}) \NC = \bigl\{\,\psi\in L_2(\mathbb{R}^2)\otimes \mathbb{C}^2\bigm| \pform(\psi,\psi)<\infty\,\bigr\};\NR[keq:pformmax]
\NC \pform_{\text{max}}(\psi,\psi) \NC = \pform(\psi,\psi),\quad\psi\in\dom(\pform_{\text{max}}).\NR
\stopalign
\stopformula

\placeformula[keq:anncre]
\startformula
\ann =-2i e^{W}\PD{}{\bar{z}} e^{-W}\quad\text{and}\quad \cre = -2i e^{-W}\PD{}{z} e^{W},
\stopformula

\placeformula[keq:pformw]
\startformula
\pform(\psi,\psi) = 4\int_{\mathbb{R}^2} \Big|\PD{}{\bz}\left(e^{-W}\psip\right)\Big|^2e^{2W}+\Big|\PD{}{z}\left(e^{W}\psim\right)\Big|^2 e^{-2W}\ed m(x),
\stopformula





\placeformula
\startformula
H = (-i\nabla-\vec{a})^2.
\stopformula
It is an unbounded operator, so one should be more specific about the domain $\dom(H)$. To be able to define $H$ on $C^\infty_0(\mathbb{R}^{n})$ it is sufficient that
\placeformula[keq:areq]
\startformula
\vec{a}\in L_{4,\text{loc}}(\mathbb{R}^{n})\otimes \mathbb{R}^{n},\quad \div\vec{a}\in L_{2,\text{loc}}(\mathbb{R}^{n}).
\stopformula
This follows by expanding $H$ as
\startformula
H = -\Delta + i\div \vec{a}+2i\vec{a}\cdot\nabla+\vec{a}\cdot\vec{a}.
\stopformula

\placeformula[keq:sform]
\startformula
h(u,u)=\int_{\mathbb{R}^n} \left|(-i\nabla-\vec{a})u\right|^2\ed m(x).
\stopformula
Assuming that $\vec{a}\in L_{2,\text{loc}}(\mathbb{R}^n)\otimes \mathbb{R}^n$, we can define two forms $h_{\min}$ and $h_{\max}$ as
\startformula
\startalign
\NC \dom(h_{\text{min}}) \NC = C_0^\infty(\mathbb{R}^n);\NR
\NC h_{\text{min}}(u,u)\NC =h(u,u),\quad u\in\dom(h_{\text{min}});\NR
\stopalign
\stopformula

\placeformula
\startformula
\startalign
\NC \dom(h_{\text{max}}) \NC = \bigl\{\,u\in L_2(\mathbb{R}^n)\bigm| h(u,u)<\infty\,\bigr\}\NR
\NC h_{\text{max}}(u,u)\NC =h(u,u),\quad u\in\dom(h_{\text{max}}).\NR
\stopalign
\stopformula

\placeformula
\startformula
e^{if}(-i\nabla-\vec{a}_1)^2 e^{-if} = (-i\nabla -\vec{a}_2)^2.
\stopformula

\placeformula
\startformula
B(x)=2\pi\alpha \delta(x) \ed x^1\wedge \ed x^2.
\stopformula

\placeformula
\startformula
\vec{a}(x)=\frac{\alpha}{|x|^2}\bigl(-x^2,x^1\bigr).
\stopformula

\placeformula
\startformula
\startalign
\NC h(u,u) \NC =\int_{\mathbb{R}^2} |(-i\nabla-\vec{a})u|^2\ed m(x) \geq \int_{|x|<1}|(-i\nabla-\vec{a})u|^2\ed m(x)\NR
\NC \NC = \int_{|x|<1}|\vec{a}|^2\ed m(x) = \int_{|x|<1}\frac{\alpha^2}{|x|^2}\ed m(x) = +\infty\NR
\stopalign
\stopformula

\placeformula
\startformula
u(re^{i\theta})\sim c_{-\alpha}r^{-\alpha}+c_{\alpha-1}r^{\alpha-1}e^{-i\theta}+c_{\alpha}r^{\alpha}+c_{1-\alpha}r^{1-\alpha}e^{-i\theta},\quad r\searrow 0, 
\stopformula

\placeformula
\startformula
N\bigl(\Lambda_q\pm\lambda,\mu_{\pm},H(\vec{a}_0,\pm V)\bigr) \sim \frac{|\log\lambda|}{\log|\log\lambda|},\quad\lambda \searrow 0,
\stopformula

Moreover, if $W_\infty=0$, then they are able to prove a similar formula for the higher Landau levels,
\startformula
\lim_{j\to\infty} \bigl(\pm j!(\lambda_{j,q}^{\pm}-\Lambda_q\bigr)^{1/j} = \frac{B_0}{2}\kap(K)^2,\quad q=0,1,\ldots.
\stopformula

\placeformula
\startformula
\lambda_{1,q}^+\geq\lambda_{2,q}^+\geq 
\cdots,\qquad \lambda_{1,q}^-\leq\lambda_{2,q}^-\leq \cdots
\stopformula
be the eigenvalues of $H(\vec{a},\pm V)$ in $(\Lambda_{q},\Lambda_{q+1})$ (for~$+$) and $(\Lambda_{q-1},\Lambda_q)$ (for~$-$), and let $\kap(K)$ be the logarithmic capacity of the set $K$, see~\cite[lan]. Then, if the magnetic scalar potential can be written as $W=-\frac{B_0}4|z|^2+W_\infty$, where $W_\infty$ is bounded, and $V$ has support in a compact $K$, $V\geq c>0$ on $K$, it holds that
\startformula
\lim_{j\to\infty} \bigl(\pm j!(\lambda_{j,0}^{\pm}-\Lambda_0)\bigr)^{1/j} = \frac{B_0}{2}\kap(K)^2.
\stopformula
Moreover, if $W_\infty=0$, then they are able to prove a similar formula for the higher Landau levels,
\startformula
\lim_{j\to\infty} \bigl(\pm j!(\lambda_{j,q}^{\pm}-\Lambda_q\bigr)^{1/j} = \frac{B_0}{2}\kap(K)^2,\quad q=0,1,\ldots.
\stopformula

\placeformula
\startformula
\startalign
\NC N(-\infty,\Lambda_0,H)\NC \leq 2n,\NR
\NC N(\Lambda_q,\Lambda_{q+1},H)\NC \leq qn,\quad q=0,1,\ldots.\NR
\stopalign
\stopformula

\placeformula
\startformula
\lim_{j\to\infty} \bigl(j!(\lambda_{j,q}^{+}-\Lambda_q)\bigr)^{1/j} = \frac{B_0}{2}\kap(K)^2,\quad q=0,1,\ldots.
\stopformula

\placeformula
\startformula
\lim_{j\to\infty} \bigl(j!(\Lambda_q-\lambda_{j,q}^{-})\bigr)^{1/j} = \frac{B_0}{2}\kap(K)^2,\quad q=0,1,\ldots.
\stopformula

\placeformula[keq:gpauli]
\startformula
\pauli = 
\startpmatrix
\NC \vphantom{\frac{g-2}{2}B}\cre\ann \NC 0 \NR
\NC 0 \NC \vphantom{\frac{g-2}{2}B}\ann\cre\NR
\stoppmatrix
+
\startpmatrix
\NC -\frac{g-2}{2}B \NC 0 \NR
\NC 0 \NC \frac{g-2}{2}B\NR
\stoppmatrix,
\stopformula

\placeformula
\startproblem (This problem was proposed to me by Prof. Ari Laptev in private communication) Study Schrödinger operators with singular potentials in higher dimensions. For example, one could try to describe the self-adjoint extensions of the Schrödinger operator in $\mathbb{R}^{2d}$ with magnetic one-form
\startformula
a(x) = \sum_{j\neq k} \Phi_{j,k}\frac{x^{2k}-x^{2j}}{|z^j-z^k|^2}\ed x^{2j-1} + \Phi_{j,k} \frac{x^{2j-1}-x^{2k-1}}{|z^j-z^k|^2}\ed x^{2j}
\stopformula
initially defined on $C_0^\infty\bigl(\mathbb{R}^{2d}\setminus \{z^j=z^k\}_{j\neq k}\bigr)$.

The main difficulty is that the vector potential is singular, not just in one point as in two dimensions, but in all hyperplanes $z^j=z^k$. In two dimensions the extensions were described by certain singular boundary terms at the singular points. 
\stopproblem




[-- Attachment #3: mickep-math.log --]
[-- Type: text/x-log, Size: 4322 bytes --]

This is pdfTeX, Version 3.1415926-1.40.9 (Web2C 7.5.7) (format=cont-en 2009.3.11)  13 MAR 2009 21:37
entering extended mode
\write18 enabled.
 (/home/mbana/context-minimals/tex/texmf-context/web2c/natural.tcx)
**mickep-math.tex \emergencyend
(./mickep-math.tex

ConTeXt  ver: 2009.03.11 09:45 MKII  fmt: 2009.3.11  int: english/english

system          : cont-new loaded
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/cont-new.tex
systems         : beware: some patches loaded from cont-new.tex
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/cont-new.mkii)
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/cont-mtx.tex))
system          : cont-fil loaded
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/cont-fil.tex
loading         : Context File Synonyms
)
system          : cont-sys.rme loaded
(/home/mbana/context-minimals/tex/texmf-context/tex/context/user/cont-sys.rme
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/type-tmf.tex)
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/type-siz.tex)
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/type-one.tex)
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/type-akb.tex))
system          : mickep-math.top loaded
(./mickep-math.top)
bodyfont        : 12pt rm is loaded
language        : patterns nl->texnansi:texnansi->1->2:2 nl->ec:ec->2->2:2 fr->
texnansi:texnansi->3->2:2 fr->ec:ec->4->2:2 de->texnansi:texnansi->5->2:2 de->e
c:ec->6->2:2 it->texnansi:texnansi->7->2:2 it->ec:ec->8->2:2 pt->texnansi:texna
nsi->9->2:2 pt->ec:ec->10->2:2 hr->ec:ec->11->2:2 pl->pl0:pl0->12->2:2 pl->ec:e
c->13->2:2 pl->qx:qx->14->2:2 cs->il2:il2->15->2:2 cs->ec:ec->16->2:2 sk->il2:i
l2->17->2:2 sk->ec:ec->18->2:2 sl->ec:ec->19->2:2 ru->t2a:t2a->21->2:2 gb->ec:e
c->22->2:2 us->ec:ec->23->2:2 agr->agr:agr->24->2:2 da->ec:ec->25->2:2 sv->ec:e
c->26->2:2 af->ec:ec->27->2:2 nb->ec:ec->28->2:2 nn->ec:ec->29->2:2 deo->ec:ec-
>30->2:2 es->ec:ec->35->2:2 ca->ec:ec->36->2:2 la->ec:ec->37->2:2 ro->ec:ec->38
->2:2 tr->ec:ec->39->2:2 fi->ec:ec->41->2:2 hu->ec:ec->42->2:2 loaded
specials        : dvips loaded
\openout3 = `mickep-math.tui'.

\openout0 = `mickep-math-mpgraph.mp'.

\openout0 = `mpgraph.mp'.

systems         : system commands are enabled
language        : language en is active
specials        : loading definition file tpd
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/spec-tpd.tex
specials        : loading definition file fdf
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/spec-fdf.tex
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/spec-fdf.mkii)
)
(/home/mbana/context-minimals/tex/texmf-context/tex/context/base/spec-tpd.mkii)
)
system          : module bib loaded
(/home/mbana/context-minimals/tex/texmf-context/tex/context/bib/t-bib.tex
publications    : loading formatting style from bibl-apa
(/home/mbana/context-minimals/tex/texmf-context/tex/context/bib/bibl-apa.tex))
subsubsection   : 1 The Pauli operator in $\mathbb {R}^2$
(./mickep-math.tuo
system          : new version of utility file, second pass needed
)
references      : unknown reference [][keq:pauli]
references      : unknown reference [][keq:komm]
publications    : warning: cite argument bovo is unknown on 170
! Missing number, treated as zero.
<to be read again> 
                   $
\@@dobig ...o #1\bodyfontsize {}\right .\n@space $
                                                  }}
\@mt@defaultBigl ...\puremathcomm {open}{\Big {#1}
                                                  }
l.99 \pauli = \dirac^2 = \Bigl(
                               \sum_{j=1}^2 \sigma_j\big(-i\PD{}{x^j}-a_j\bi...

? X
 
Here is how much of TeX's memory you used:
 2089 strings out of 257043
 34980 string characters out of 1271064
 2692006 words of memory out of 4125384
 41891 multiletter control sequences out of 10000+100000
 166255 words of font info for 80 fonts, out of 2000000 for 5000
 486 hyphenation exceptions out of 8191
 35i,10n,33p,409b,619s stack positions out of 10000i,500n,10000p,4000000b,50000s
No pages of output.
PDF statistics:
 0 PDF objects out of 1000 (max. 8388607)
 0 named destinations out of 1000 (max. 131072)
 1 words of extra memory for PDF output out of 10000 (max. 10000000)




[-- Attachment #4: Type: text/plain, Size: 487 bytes --]

___________________________________________________________________________________
If your question is of interest to others as well, please add an entry to the Wiki!

maillist : ntg-context@ntg.nl / http://www.ntg.nl/mailman/listinfo/ntg-context
webpage  : http://www.pragma-ade.nl / http://tex.aanhet.net
archive  : https://foundry.supelec.fr/projects/contextrev/
wiki     : http://contextgarden.net
___________________________________________________________________________________