\documentclass[reqno]{amsart}
\usepackage{hyperref}
\AtBeginDocument{{\noindent\small
{\em Electronic Journal of Differential Equations},
Vol. 2007(2007), No. 142, pp. 1--9.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu (login: ftp)}
\thanks{\copyright 2007 Texas State University - San Marcos.}
\vspace{9mm}}
\begin{document}
\title[\hfilneg EJDE-2007/142\hfil Ambarzumian's theorem for trees]
{Ambarzumian's theorem for trees}
\author[R. Carlson, V. Pivovarchik\hfil EJDE-2007/142\hfilneg]
{Robert Carlson, Vyacheslav Pivovarchik} % in alphabetical order
\address{Robert Carlson \newline
University of Colorado at Colorado Springs,
Colorado Springs, CO 80921, USA}
\email{rcarlson@uccs.edu}
\address{Vyacheslav Pivovarchik \newline
South-Ukrainian State Pedagogical University,
Staroportofrankovskaya Ul. 26, 65020 Odessa, Ukraine}
\email{v.pivovarchik@paco.net \quad vnp.@dtp.odessa.ua}
\thanks{Submitted April 12, 2007. Published October 24, 2007.}
\subjclass[2000]{34B45}
\keywords{Inverse eigenvalue problem; quantum graph}
\begin{abstract}
The classical Ambarzumian's Theorem for Schr\"odinger operators
$-D^2 + q$ on an interval, with Neumann conditions at the endpoints,
says that if the spectrum of $(-D^2+q)$ is the same as the spectrum of
$(-D^2)$ then $q=0$. This theorem is generalized to Schr\"odinger operators
on metric trees with Neumann conditions at the boundary vertices.
\end{abstract}
\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\section{Introduction}
In 1929, Ambarzumian \cite{Ambar} published
the following foundational result in inverse spectral
theory.
\begin{theorem}[Ambarzumian] \label{thmo}
If $q(x)$ is real-valued and continuous function,
and the spectrum of the boundary problem
\begin{gather*}
-y''+ q(x)y=\lambda y, \\
y'(0) = y'(1) = 0,
\end{gather*}
is $k^2\pi ^2$, for $k = 0,1,2,\dots $, then $q(x) = 0$.
\end{theorem}
As Borg \cite{Borg1,Borg2} established, the case $q=0$ turns out
to be exceptional. In general the spectra of two Sturm-Liouville
problems are needed to determine the potential. However, some
generalizations of Ambarzumian's theorem were obtained in
\cite{Chak,Chern,Chern2,Harrell}. In particular, it is known
\cite{Harrell} that one may relax the original hypotheses,
assuming that $q \in L^1[0,1]$, and
\[
\lim_{k \to \infty} (\lambda _k - k^2\pi ^2) = 0.
\]
This work establishes an extension of Ambarzumian's theorem to trees, where
each (directed) graph edge $e_1,\dots ,e_n$ is identified with the
interval $[0,1]$. Functions on the graph may then be identified with
an ordered $n$-tuple of functions on $[0,1]$,
With these identifications, the eigenvalue equation
$-Y'' + Q(x) Y = \lambda Y$ for the tree
with the real integrable diagonal matrix potential
$Q = {\rm diag}[q_1,\dots , q_n]$
is a system of $n$ scalar equations
\begin{equation} \label{eeqn}
-y_j'' + q_j(x)y_j = \lambda y_j, \quad q_j \in L^1[0,1].
\end{equation}
Eigenfunctions for $-Y'' + Q(x) Y = \lambda Y$ on a graph are
(nonzero) solutions of \eqref{eeqn} which satisfy vertex
(boundary) conditions. For trees, the Neumann condition $y'(v) =
0$ is assumed to hold at the boundary vertices, which have one
incident edge. At interior vertices, with more than one incident
edge, solutions of \eqref{eeqn} are required to be continuous and
satisfy the Kirchhoff condition, that the derivatives sum to zero
in outward pointing local coordinates. For a star shaped graph a
version of Ambarzumian's theorem was obtained in \cite{Pivo}. Our
generalization is as follows.
\begin{theorem} \label{main}
Suppose $T$ is a finite tree with all edges of length $1$.
For $r = 1,2,\dots $, let $\{ m_r \}$ be a sequence of integers with
$\lim_{r \to \infty} m_r = \infty $.
If the set of eigenvalues for \eqref{eeqn} satisfying the tree vertex
conditions
is nonnegative, and contains a subsequence $\{ \lambda _r \} $ with
\begin{equation} \label{noshift}
\lim_{r \to \infty} \bigl ( \lambda _r - (\pi m_r)^2 \bigr ) = 0,
\end{equation}
then $q_j(x)=0$ a.e., for each $j=1,2,\dots , n$.
\end{theorem}
The following observations may provide motivation for the
subsequent arguments. Consider the equation $-Y'' = \lambda Y$ on
a tree $T$ whose edges have length $1$. If $v$ is a boundary
vertex, and $x$ is the distance in $T$ to $v$, then for $m =
0,1,2,\dots $ the function $\cos(m\pi x)$ will be continuous with
derivative $0$ at all the vertices. Thus for the case $Q(x) = 0$
the spectrum contains a subsequence $(m\pi)^2$, and these
eigenvalues are easily seen to be simple.
If $Q \not= 0$, a number of arguments \cite{Harrell}, \cite[p. 33-35]{Pos}
suggest that one should expect a spectral shift given by $\int _T Q$.
The hypothesized behavior \eqref{noshift} suggests that
$\int _T Q = 0$.
A previously known argument \cite{Harrell,LG}
would then show that positivity of all eigenvalues forces $Q$ to be $0$.
The proof of Theorem \ref{main} uses an analysis of the asymptotic
behavior of eigenfunctions for \eqref{eeqn} when the eigenvalues
are close to $(\pi m_r)^2$. This material, along with a review of
differential equations on graphs, is provided in Section 2.
Section 3 contains three results, which together comprise the
proof of Theorem \ref{main}. The expected result $\int _T Q = 0$
is first established. Second, a version of an argument of
\cite{LG} is used to complete the proof, except for a technical
issue about the form domain for $-D^2 + Q$ when some $q_j$ is
integrable, but not square integrable. The third result is
included to resolve this technical issue.
\section{Background Material}
Let $T$ be a metric tree with $n$ edges, all of length 1, and
an interior vertex chosen as the root.
Local coordinates for the edges identify each edge with $[0,1]$
so that the local coordinate increases as the distance to the root decreases.
This means that each boundary vertex has local coordinate $0$,
the root has local coordinate $1$ on each of its incident edges,
and all other interior vertices $v$ have one outgoing edge, with
local coordinate $0$ for $v$, while the local coordinate for $v$
is $1$ on the remaining incoming incident edges.
The given choice of local coordinates for $T$ provides a convenient form
for the vertex conditions. For an edge $e_j$ incident on a boundary vertex,
the Neumann condition is
\[
y_j'(0)=0.
\]
For the root we impose the continuity conditions
\[
y_j(1)=y_k(1)
\]
for all incident edges $e_j$ and $e_k$ and the Kirchhoff condition
\[
\sum _j y_j'(1)=0,
\]
the sum taken over all edges $e_j$ incident on the root.
For all other interior vertices $v$ with incoming edges $e_j$ and outgoing edge $e_k$
the continuity conditions are
\[
y_j(1)=y_k(0),
\]
and the Kirchhoff condition is
\[
y'_k(0)=\sum _{j} y_j'(1).
\]
Denote by $c_j(x,\lambda )$ the solution of \eqref{eeqn} which
satisfies the conditions $c_j(0,\lambda )-1=c_j'(0, \lambda )=0$
and by $s_j(x, \lambda )$ the solution of \eqref{eeqn} which
satisfies the conditions $s_j(0,\lambda ) = s_j'(0,\lambda ) - 1 =
0$. These solutions and their $x$ derivatives satisfy the integral
equations
\begin{gather*}
c(x,\lambda ) = \cos (\sqrt{\lambda }x) + \int_0^x
\frac{\sin(\sqrt{\lambda }(x-t))}{\sqrt{\lambda }} q(t)
c(t,\lambda ) \,dt , \\
c'(x,\lambda ) = -\sqrt{\lambda }\sin (\sqrt{\lambda }x)
+ \int_0^x \cos(\sqrt{\lambda }(x-t)) q(t) c(t,\lambda ) \,dt , \\
s(x,\lambda ) = \frac{\sin (\sqrt{\lambda }x)}{\sqrt{\lambda }}
+ \int_0^x \frac{\sin(\sqrt{\lambda }(x-t))}{\sqrt{\lambda }} q(t)
s(t,\lambda ) \,dt , \\
s'(x,\lambda ) = \cos (\sqrt{\lambda }x) + \int_0^x \cos(\sqrt{\lambda }(x-t))
q(t) s(t,\lambda ) \,dt .
\end{gather*}
As a consequence, for $\lambda > 0$ one obtains the well know estimates
\begin{equation} \label{estimates1}
c(x,\lambda ) = \cos(\sqrt{\lambda }x) + O(1/\sqrt{\lambda}), \quad
s(x,\lambda ) = \frac{\sin(\sqrt{\lambda }x)}{\sqrt{\lambda}}
+ O(1/\lambda ).
\end{equation}
These results are established in \cite[p. 13]{Pos} for $q \in L^2[0,1]$, but the same ideas
work for $q \in L^1[0,1]$.
Suppose that $Y(x,\lambda )$ is a vector function whose components $y_j$
satisfy \eqref{eeqn},
and which is given the graph $L^2$ norm,
\[
\| Y(x,\lambda ) \| ^2 = \sum_{j=1}^n \int_0^1 |y_j(x,\lambda )|^2 \,dx .
\]
Each $y_j(x,\lambda )$ may be written as a linear combination
\[
y_j(x,\lambda ) = A_j(\lambda )c_j(x,\lambda ) + B_j(\lambda )s_j(x,\lambda ) .
\]
\begin{lemma} \label{bound}
Suppose $q_j \in L^1[0,1]$ and $\| Y(x,\lambda ) \| = 1$.
Then there is a $\lambda _0 > 0$ and a constant $C$ such that
\[
|A_j(\lambda )| \le C, \quad |B_j(\lambda )/\sqrt{\lambda }|
\le C, \quad \lambda \ge \lambda _0.
\]
\end{lemma}
\begin{proof}
The bound on $\| Y(x,\lambda ) \|$ gives a bound for each
$y_j(x,\lambda )$, so
\begin{align*}
& \int_0^1 | A_j(\lambda )c_j(x,\lambda ) +
B_j(\lambda ) s_j(x,\lambda )|^2 \,dx \\
&= \int_0^1 | A_j|^2c_j^2 + \frac{|B_j|^2}{\lambda } \lambda s_j^2
+ (A_j \frac{\overline{B_j}}{\sqrt{\lambda }}
+ \overline{A_j}\frac{B_j}{\sqrt{\lambda }}) c_j \sqrt{\lambda} s_j \le 1 .
\end{align*}
Simple trigonometric integrals and \eqref{estimates1} give
\[
\int_0^1 c_j^2 = \frac{1}{2} + O(\lambda ^{-1/2}) , \quad
\int_0^1 \lambda s_j^2 = \frac{1}{2} + O(\lambda ^{-1/2}) , \quad
\int_0^1 \sqrt{\lambda }s_jc_j = O(\lambda ^{-1/2}) .
\]
Since
\[
2|\overline{A_j}\frac{B_j}{\sqrt{\lambda }}| \le |A_j|^2
+ \frac{|B_j|^2}{\lambda } ,
\]
we have
\[
|A_j|^2[\frac{1}{2} - O(\lambda ^{-1/2})]
+ \frac{|B_j|^2 }{\lambda }[\frac{1}{2} - O(\lambda ^{-1/2})]
\le 1,
\]
which provides the desired estimate.
\end{proof}
Some refinements of the estimates \eqref{estimates1} are available
if $x = 1$ and $\lambda $ is close to $m^2\pi ^2$.
A recent reference is \cite{Carlson04},
but the ideas have a long history, \cite[p. 11]{Borg1}.
Define
\[
[q_j] = \int_0^1 q_j(t) \,dt ,\quad\text{and}\quad
\omega _j = \sqrt{\lambda - [q_j}] .
\]
\begin{lemma}
Suppose that for positive integers $m$ we have $ \lambda = m^2\pi ^2 + O(1)$.
For edge indices $j = 1,\dots ,n$, the following estimates hold uniformly
if $\int_0^1 |q_j(t)| \,dt $ is bounded:
\begin{equation} \label{refest}
\begin{gathered}
c_j(1,\lambda ) = \cos(\omega _j) - (-1)^m2^{-1}m^{-1} \int_0^1 \sin(2m\pi t)
q_j(t) \,dt + O(m^{-2}),
\\
c_j'(1,\lambda ) = -\omega _j\sin(\omega _j) + (-1)^m2^{-1} \int_0^1
\cos(2m\pi t) q_j(t) \,dt + O(m^{-1}), \\
s_j(1,\lambda ) = \omega _j^{-1}\sin(\omega _j) + (-1)^m2^{-1}m^{-2}
\int_0^1 \cos(2m\pi t) q_j(t) \,dt + O(m^{-3}), \\
s_j'(1,\lambda ) = \cos(\omega _j) + (-1)^m2^{-1}m^{-1} \int_0^1
\sin(2m\pi t) q_j(t) \,dt + O(m^{-2}).
\end{gathered}
\end{equation}
\end{lemma}
Suppose $\lambda $ has the form $ \lambda = m^2\pi ^2 + O(1)$, and
$m \to \infty $. Using \eqref{refest},
the Riemann-Lebesgue lemma and elementary trigonometric identities give
\begin{equation} \label{refest2}
\begin{gathered}
c_j((1, \lambda )= (-1)^m + o(1), \\
c_j'(1, \lambda ) = (-1)^m \frac{1}{2}\int_0^1 q_j(x) \,dx + o(1), \\
s_j(1, \lambda ) = o(m^{-1}), \\
s_j'(1, \lambda )= (-1)^m + o(1).
\end{gathered}
\end{equation}
At the end of the proof we will use some Hilbert space operator
theory, so it will be helpful to recall \cite{Kuchment} that on
finite graphs the formal Schr\"odinger operators $-D^2 + Q$ with
real-valued edge potentials $q_j(x)$ in $L^1[0,1]$, together with
the interior and boundary vertex conditions described above, are
associated with a self adjoint operator $\mathcal{L}$, with compact
resolvent, acting componentwise on the Hilbert space $\oplus _j
L^2(e_j)$ of square integrable functions on the edges. These facts
may be established by using the variation of parameters formula to
construct the resolvent. The assumption that all eigenvalues are
nonnegative in Theorem \ref{main} then implies that $\mathcal{L}$ is
a nonnegative operator.
\section{Proof of Theorem \ref{main}}
This proof will be split into several parts.
\begin{lemma}
With the hypotheses of Theorem \ref{main},
\begin{equation} \label{qbarsum}
\sum_j [q_j] = \sum_j \int_0^1 q_j(t) \,dt = 0 .
\end{equation}
\end{lemma}
\begin{proof}
The sequence $\{ m_r \}$ has either an infinite subsequence of even integers,
or an infinite subsequence of odd integers.
We will assume the first case holds,
and use this subsequence instead of the original sequence with no change of
notation.
If an odd subsequence were used instead, the leading terms in \eqref{refest2}
would change sign, but this has no significant impact on the proof.
Suppose that $\{ Y(x,\lambda _r ) \} $ is a sequence of eigenfunctions
for \eqref{eeqn}
with norm $1$. Write the components $y_j(x,\lambda _r)$ as a linear combination
\[
y_j(x,\lambda _r ) = A_j(\lambda _r)c_j(x,\lambda _r)
+ B_j(\lambda _r)s_j(x,\lambda _r) .
\]
Recall that the coefficients $A_j(\lambda _r)$ and
$B_j(\lambda _r)/\sqrt{\lambda _r}$ are
bounded sequences by Lemma \ref{bound}. The proof now proceeds in three steps.
First, consider the value of $y_j$ for edges incident on a vertex $v$.
If $j$ is the index of an outgoing edge, then
\[
y_j(0,\lambda _r) = A_j(\lambda _r).
\]
If $j$ is the index of an incoming edge, then by \eqref{refest2}
\[
y_j(1,\lambda _r) = A_j(\lambda _r)c_j(1,\lambda _r)
+ B_j(\lambda _r)s_j(1,\lambda _r) = A_j(\lambda _r ) + o(1).
\]
The continuity of $Y$ at the vertex $v$ thus implies
\begin{equation} \label{Aeq}
A_j(\lambda _r) = A_k(\lambda _r) + o(1), \quad r \to \infty,
\end{equation}
for all edges $j,k$ incident on $v$. Since the tree $T$ is connected,
\eqref{Aeq} extends to all edges $j,k$.
Second, considering the root vertex to be at the top of the tree $T$, say that an edge
$e_j$ is below an edge $e_k \not= e_j$ if a path from $e_j$ to the root passes through $e_k$.
Label each vertex $v$ with the combinatorial distance from the root,
and label the edges with the larger of the vertex labels on the edge.
Let $M$ be the maximum label.
If the vertex $v$ has label $M-1$, then all its incoming edges $e_k$
join $v$ to a boundary vertex.
If $e_j$ is the outgoing edge for $v$, then the derivative condition at $v$ gives
\[
B_j(\lambda _r) = \sum_k A_k(\lambda _r) c_k'(1,\lambda _r) =
\frac{1}{2} \sum_k A_k(\lambda _r) [q_k] + o(1) .
\]
For vertices $v$ with label $M-2$, and outgoing edge $e_j$,
the derivative condition at $v$ and \eqref{refest2} give
\begin{equation} \label{Beqn}
B_j(\lambda _r) = \sum_k [A_k(\lambda _r) c_k'(1,\lambda _r)
+ B_ks_k'(1,\lambda _r)]
= \frac{1}{2} \sum_l A_l(\lambda _r) [q_l] + o(1) ,
\end{equation}
where the last sum is over edges $e_l$ which are below $e_j$.
This pattern continues as we reduce the vertex label until we reach the root.
The derivative condition at the root now gives
\[
0 = \sum_j A_l(\lambda _r) [q_j] + o(1) =
A_1(\lambda _r) \sum_j [q_j] + o(1) ,
\]
the sum taken over all edges $e_j$.
Third, from \eqref{Beqn} and the extension of \eqref{Aeq} to all edges
$j,k$, we conclude that
$\liminf_{r \to \infty} |A_j(\lambda _r)| > 0$ for all $j$,
since otherwise we would have
\[
\liminf_{r \to \infty} |y_j(x,\lambda _r)| = 0
\]
uniformly in $x$, and
the condition $\| Y (\lambda _r) \| = 1$ would be violated.
In particular, $A_1(\lambda _r)$ is bounded away from zero as $r \to \infty $,
so we obtain
$\sum_j [q_j] = 0$.
\end{proof}
The next step in the proof of Theorem \ref{main} uses an argument
from \cite{LG}, which is extended in \cite{Harrell}. Recall
\cite[pp. 318]{Kato} that the nonnegative operator $\mathcal{L}$ has
an associated nonnegative form
\[
\mathbf{t}[f,g] = \langle \mathcal{L}f,g \rangle
= \sum_j \int_0^1 (\mathcal{L}f_j)\overline{g_j}
= \sum_j \int_0^1 (f_j' \overline{g_j'} + q_jf_j \overline{g_j}) ,
\]
defined for $f,g$ in the domain of $\mathcal{L}$.
The closure of this form is
\[
\mathbf{t}[f,g] = \langle \mathcal{L}^{1/2}f,\mathcal{L}^{1/2}g \rangle ,
\]
with the form domain equal to the domain of $\mathcal{L}^{1/2}$,
the nonnegative square root of $\mathcal{L}$.
\begin{lemma} \label{lem3.2}
With the hypotheses of Theorem \ref{main}, if the function $1$ is
in the domain of $\mathcal{L}^{1/2}$, then $q_j(x) = 0$ almost everywhere
for all edges $e_j$.
\end{lemma}
\begin{proof}
Recalling that the lowest eigenvalue $\lambda _0$ of $\mathcal{L}$ is nonnegative,
we apply the minimax principle for our differential operator $\mathcal{L}$.
For $Y$ in the domain of $\mathcal{L}^{1/2}$ we have
\begin{equation} \label{vareqn}
\begin{aligned}
0 \le \lambda _0
&= \min _{||Y||=1}\langle \mathcal{L}^{1/2}Y,\mathcal{L}^{1/2}Y \rangle \\
&= \min _{||Y||=1}\Big(\sum
_{j=1}^{n}\int_0^1 y_j'\overline{y_j'}dx + \sum _{j=1}^{n}\int_0^1
q_j(x)|y_j|^2dx\Big).
\end{aligned}
\end{equation}
The constant function $Y_0 =(1,1,\dots ,1)^t$ is in the domain of
$\mathcal{L}^{1/2}$, and the combination
of \eqref{vareqn} and \eqref{qbarsum} shows that $0$ is an
eigenvalue of $\mathcal{L}$ with eigenfunction $Y_0$. The proof is
completed by noticing that the equation $\mathcal{L}Y_0 = 0$ gives
$q_j(x)=0$ a.e., for all edges $e_j$.
\end{proof}
Let $\mathcal{D}_0$ denote the domain of the self adjoint operator
$\mathcal{L}_0 = -D^2$ associated with the continuity and Kirchhoff
boundary conditions described above. That is \cite{Carlson98},
$\mathcal{D}_0$ consists of those vector functions
$G = (g_1,\dots ,g_n)$ with $g_j$, $g_j'$ absolutely continuous,
with $g_j'' \in
L^2[0,1]$, and which satisfy the continuity and Kirchhoff boundary
conditions. The function $1$ is in the domain of $\mathcal{L}_0$.
Since the domain of $\mathcal{L}^{1/2}$ includes the domain of
$\mathcal{L}$, the function $1$ will be in the domain of $\mathcal{L}^{1/2}$,
thus yielding the conclusion of Theorem \ref{main}, if the domain
of $\mathcal{L}$ is $\mathcal{D}_0$, which will hold under the stronger
hypothesis $q_j \in L^2[0,1]$.
A more elaborate argument is required for $q_j \in L^1[0,1]$.
Define the symmetric form $\mathbf{t}_1$ with domain $\mathcal{D}_0$,
\[
\mathbf{t}_1[F,G] = \sum_{j=1}^n \int_0^1 \Bigl ( f_j'\overline{g_j'}
+ q_jf_j\overline{g_j} \Bigr ) .
\]
The proof of Theorem \ref{main}
will be completed with the next result, whose proof closely
follows the argument in \cite[pp. 345-346]{Kato}, where additional
details are provided.
\begin{theorem} \label{thm2}
The symmetric form $\mathbf{t}_1$ is closable. The self adjoint operator
associated with the
closure of $\mathbf{t}_1$ is $\mathcal{L} = -D^2 + Q$,
with a domain $\mathcal{D}$ consisting of vector functions
$F:[0,1] \to \mathbb{C} ^n$
whose components $f_j$ satisfy the following conditions:
\begin{itemize}
\item[(i)] $f_j$ and $f_j'$ are absolutely continuous, and
$-f_j'' + q_j(x)f_j \in L^2[0,1]$,
\item[(ii)] $f_j$ satisfies the Neumann boundary conditions $f_j'(v) = 0$
at boundary vertices,
and the continuity conditions and Kirchhoff conditions at interior vertices.
\end{itemize}
\end{theorem}
\begin{proof}
In addition to $\mathbf{t}_1$, define a second symmetric form with
domain $\mathcal{D}_0$,
\[
\mathbf{t}_0[F,G] = \sum_{j=1}^n \int_0^1 f_j'(x)\overline{g_j'(x)} \,dx.
\]
The form $\mathbf{t}_0$ is nonnegative, and so closable.
For any $\epsilon > 0$ there is a $\delta > 0$ such that
\[
|f_j(x)|^2 \le \epsilon \int_0^1 |f_j'(x)|^2 \,dx
+ \delta \int_0^1 |f_j(x)|^2 \,dx .
\]
This implies that $\mathbf{t}_1$ is the sum of $\mathbf{t}_0$ and a form
with $\mathbf{t}_0$
bound $0$, so that \cite[p. 320]{Kato} $\mathbf{t}_1$ is bounded below
and closable, and the closures of $\mathbf{t}_0$ and $\mathbf{t}_1$
have the same domain $\mathcal{D}_1$.
Let $x$ and $y$ lie on edges incident on the vertex $v$.
Integrating from $x$ to $v$ and then from $v$ to $y$ on the respective edges,
the continuity of $F \in \mathcal{D}_0$ at $v$ gives
\[
F(y) - F(x) = \int_x^y F'(t) \,dt.
\]
This formula, and the implied continuity at $v$, extends to functions
$F \in \mathcal{D}_1$.
The closure of $\mathbf{t}_1$ has \cite[p. 331]{Kato} an
associated self adjoint operator $\mathcal{L}$, with domain $\mathcal{D}$.
If $F \in \mathcal{D}$ and $\mathcal{L}F = H$, then
\[
\mathbf{t}_1[F,G] = \langle \mathcal{L}F,G\rangle
= \langle H,G \rangle , \quad G \in \mathcal{D}_1.
\]
For functions $z_j$ such that
\[
z_j'(x) = h_j(x) - q_j(x)f_j(x),
\]
a rearrangement of terms from $\mathbf{t}_1[F,G] = \langle H,G \rangle $
gives
\begin{equation} \label{ibp}
\begin{aligned}
\sum_{j=1}^n \int_0^1 f_j'(x)\overline{g_j'(x)} \,dx
&= \sum_{j=1}^n \int_0^1 z_j'(x) \overline{g_j(x)} \,dx \\
& = \sum_{j=1}^n \Bigl ( z_j(1)\overline{g_j}(1) -
z_j(0)\overline{g_j}(0) \Bigr )
- \sum_{j=1}^n \int_0^1 z_j(x) \overline{g_j'(x)} \,dx.
\end{aligned}
\end{equation}
If the components of $G$ satisfy $g_j(0) = g_j(1) = 0$, then
\[
g_j(x) = \int_0^x g_j'(t) \,dt , \quad \int_0^1 g_j'(t) \,dt = 0.
\]
The functions $g_j'$ have dense span in the orthogonal complement
of the constants in $L^2[0,1]$,
so \eqref{ibp} gives
\begin{equation} \label{preac}
f_j' + z_j = c_j.
\end{equation}
This equation shows that $f_j'$ is absolutely continuous, and
\[
- f_j'' + q_jf_j = h_j.
\]
Since $\mathcal{D} \subset \mathcal{D}_1$,
functions $F \in \mathcal{D}$ are continuous at all vertices $v$,
as noted above.
Relaxing the constraints on functions $G \in \mathcal{D}_1$,
equation \eqref{ibp} gives
\[
\sum_{j=1}^n [c_j - z_j(1)]\overline{g_j}(1)
- \sum_{j=1}^n [c_j - z_j(0)]\overline{g_j}(0) = 0,
\]
or from \eqref{preac}
\[
\sum_{j=1}^n f_j'(1)\overline{g_j}(1)
- \sum_{j=1}^n f_j'(0)\overline{g_j}(0) = 0,
\]
If $G$ vanishes at each vertex except $v$, the continuity of $G$ at $v$
gives the desired derivative conditions.
For instance, if $v$ has one outgoing edge $e_k$ and incoming edges $e_j$,
we get
\[
f'_k(0) = \sum_{j} f'_j(1) .
\]
This shows that every function in $\mathcal{D}$ has the desired properties,
and one may show as in \cite[p. 328]{Kato} that these properties
characterize $\mathcal{D}$.
\end{proof}
\subsection*{Acknowledgements}
This work was partly supported by Grant UKM2-2811-OD-06 from the
U.S. Civilian Research and Development Foundation.
V. Pivovarchik is grateful to the University of Colorado at Colorado Springs
for its hospitality.
\begin{thebibliography}{10}
\bibitem{Ambar} V. A. Ambarzumian; \"Uber eine Frage der Eigenwerttheorie,
Zeitschrift f\"ur Physik, \textbf{53} (1929),690-695.
\bibitem{Borg1} G. Borg;
Eine Umkehrung der Sturm-Liouvilleschen Eigenwertaufgabe.
\emph{Acta Math.}, 78, 1946, 1-96.
\bibitem{Borg2} G. Borg; Uniqueness theorems in the spectral theory of
$y''+(\lambda -q(x))y=0$, Proc. 11th Scandinavian Congress of
Mathematicians, Johan Grundt Tanums Forlag Oslo, (1952), 276-287.
\bibitem{Carlson98} R. Carlson;
Adjoint and self adjoint differential operators on graphs.
\emph{Electronic J. Differential Equations}, 1998 (1998), No. 6,
1-10,
\bibitem{Carlson04} R. Carlson;
A Spectral Transform for the Matrix Hill's Equation,
\emph{Rocky Mountain Journal of Mathematics}, 34(3), 2004.
\bibitem{Chak} N. K. Chakravarty and S. K. Acharaya;
\emph{On an extension of the theorem of V. A. Ambarzumyan},
Proc. Royal Soc. Edinburg \textbf{A 110} (1988), 79--84 .
\bibitem{Chern} H. H. Chern and C. L. Shen, \emph{On the $n$-dimensional
Ambarzumyan's theorem}, Inverse Problems \textbf{13} (1997), 15--18.
\bibitem{Chern2} H. H. Chern, C. K. Law, and H. J. Wang, \emph{Corregendum
to: Extensions of Ambarzumyan's Theorem to general boundary conditions},
J. Math. Analysis and Applications \textbf{309} (2005), 764--768 .
\bibitem{Harrell} E.M. Harrel II;
\emph{On the extension of Ambarzumian's
inverse spectral theorem to compact symmetric spaces}, Amer. J.
Math. \textbf{109}(1987), No. 5, 787--795.
\bibitem{Kato} T. Kato;
Perturbation Theory for Linear Operators,
Springer-Verlag, 1995
\bibitem{Kuchment} P. Kuchment, \emph{Quantum graphs: I. Some basic structures},
Waves in Random Media \textbf{14} (2004), S107--S128.
\bibitem{LG} B. Levitan and M. Gasymov, \emph{Determination of a differential equation
by two of its spectra}, Russ. Math. Surveys \textbf{19} (2), 1--64 (1964).
%
\bibitem{Marchenko} V. A. Marchenko;
Sturm-Liouville operators and applications,
Birkh\"auser Verlag, Basel, 1986.
\bibitem{Pivo} V. Pivovarchik;
Ambarzumian's theorem for the Sturm-Liouville boundary value problem on
star-shaped graph. Funktsional. Anal. Prilozh. (Russian), \textbf{39}
(2005), No. 2, 78-81: English translation in Funct. Anal. Appl.
\textbf{39} (2005), No. 2, 148-151.
\bibitem{Pos} J. P\"oschel and E. Trubowitz;
Inverse Spectral Theory, Academic Press, Orlando, 1987
\end{thebibliography}
\end{document}