by Ruth Williams
In the late 1970s, Kai Lai Chung began investigating connections between probability and the (reduced) Schrödinger equation: \begin{equation}\label{RW-reds} \tfrac{1}{2} \Delta u(x) + q(x)u(x) = 0 \quad \hbox{for } x\in \mathbb{R}^d, \end{equation} where \( q \) is a real-valued Borel-measurable function on \( \mathbb{R}^d \), and \( \Delta \) is the \( d \)-dimensional Laplacian. His work in this area extended over the next 15 years or so. It included collaborations with several colleagues and students, and inspired the work of others. His book, From Brownian Motion to Schrödinger’s Equation [6], written with Zhongxin Zhao, is a compilation and refinement of much of the research conducted in this area up through 1994.
In the following, I will describe some of the background and early advances in this research involving connections with Brownian motion. A complementary article written by Michael Cranston, which also appears in this volume, focuses on related developments involving connections with conditioned Brownian motion. My account is not meant to be exhaustive, but rather to provide a sample of some of the intriguing aspects of the topic and to illustrate the pivotal role that Kai Lai Chung played in some of the developments. My description is necessarily influenced by my own personal recollections.
Background
Stimulated by Feynman’s [e1] proposed “path integral” solution of the complex time-dependent Schrödinger equation, Kac [e2], [e3] considered, for a Borel-measurable function \( q:\mathbb{R}\to \mathbb{R} \) satisfying \( q\leq 0 \), the following multiplicative functional of one-dimensional Brownian motion \( B \): \[ e_q(t) = \exp\Bigl( \int_0^t q(B_s)\, ds \Bigr) \quad\text{for } t\geq 0, \] This functional can also be defined for suitable Borel-measurable functions \( q: \mathbb{R}^d \to \mathbb{R} \), and for \( B \) a \( d \)-dimensional Brownian motion or even a \( d \)-dimensional diffusion process. Such functionals are now called Feynman–Kac functionals.
Consider a continuous bounded function \( q:\mathbb{R}^d\to \mathbb{R} \), and a continuous bounded function \( g:\mathbb{R}^d\to\mathbb{R} \). If \( \psi:[0,\infty)\times \mathbb{R}^d\to \mathbb{R} \) is a continuous bounded function, with continuous partial derivatives \( \partial \psi/\partial t \), \( \partial \psi/\partial x_i \), and \( \partial^2 \psi/\partial x_i\partial x_j \) for \( i,j=1,\ldots, d \), defined on \( (0,\infty)\times \mathbb{R}^d \), and such that the following time-dependent Schrödinger equation holds: \begin{equation} \label{RW-teds} \frac{\partial \psi (t,x)}{\partial t} = \frac {1}{2} \Delta \psi(t,x) + q(x) \psi(t, x) \quad\text{for } t > 0, x\in \mathbb{R}^d, \end{equation} with initial condition \( \psi(0,x) = g(x) \) for \( x\in\mathbb{R}^d \), and where \begin{equation} \Delta \psi (t,x) =\sum_{i=1}^d \frac{\partial^2 \psi}{\partial x_i^2} (t,x), \end{equation} then it can be shown (for example, by using Itô’s formula), that \begin{equation} \label{RW-expk} \psi(t, x) = E^x \bigl[ e_q(t) g(B_t)\bigr] \quad\text{for } t\geq 0, x\in \mathbb{R}^d. \end{equation} Here, \( E^x \) denotes the expectation operator under which \( B \) is a \( d \)-dimensional Brownian motion starting from \( x \).
Kac [e2], [e3] was interested in \eqref{RW-teds} when \( d=1 \) and \( q\leq0 \). However, rather than considering this equation directly, he worked with a reduced Schrödinger equation similar to \eqref{RW-reds} (with \( q-s \) in place of \( q \), and \( d=1 \)) obtained by formally taking the Laplace transform (with parameter \( s \)) in equation \eqref{RW-teds} to eliminate the time variable \( t \). Under mild conditions, for example when \( q \) is bounded and continuous in addition to being non-positive, Kac [e3] showed that the Laplace transform of the right member of \eqref{RW-expk} with \( g=1 \) satisfies this reduced Schrödinger equation.
In the late 1950s and early 1960s, in developing a potential theory for Markov processes, Dynkin (see [e5], Chapter XIII, §4, Theorem 13.16) and others, considered expressions of the form \begin{equation} \label{RW-fkX} E^x\Bigl[\exp\Bigl(\int_0^{\tau} q(X_s)\,ds\Bigr) f(X_{\tau})\Bigr] \quad\text{for } x\in \overline D, \end{equation} where \( X \) is a diffusion process in \( \mathbb{R}^d \), \( \tau=\inf\{ s > 0: X_s\notin D\} \) is the first exit time of \( X \) from a bounded domain \( D \) in \( \mathbb{R}^d \), \( \overline D \) is the closure of \( D \), \( f \) is a continuous function defined on the boundary of \( D \), and \( q\leq0 \) is Hölder-continuous and bounded on \( D \). Here, \( P^x \) and \( E^x \) denote respectively probability and expectation operators for \( X \) starting from \( x\in \overline D \). The domain \( D \) is regular if \[ P^x (\tau =0 ) =1 \quad \hbox{for each } x\in \partial D. \] Under suitable assumptions on \( X \) and assuming that \( D \) is regular, Dynkin showed that expressions of the form \eqref{RW-fkX} yield continuous functions on \( \overline D \) that satisfy the equation \begin{equation}\label{RW-diffreds} {\mathcal L} u(x) + q(x) u(x) = 0 \quad \hbox{for } x\in D, \end{equation} with continuous boundary values given by \( f \), where \( \mathcal L \) is the infinitesimal generator of the diffusion process \( X \). The assumption that \( q\leq 0 \) implies that the action of the stopped Feynman–Kac functional, \begin{equation} \label{RW-FKS} e_q(\tau) = \exp\Bigl( \int_0^{\tau} q(X_s)\,ds\Bigr), \end{equation} is to “kill” the diffusion process at a state dependent exponential rate given by \( -q \) up until the stopping time \( \tau \). The assumption that \( q \) is non-positive ensures that the mean value of the stopped Feynman–Kac functional is always finite; in fact, it is bounded by one. In contrast, Khasminskii [e4] considered the case where \( q \) is non-negative. In this case, the action of the stopped Feynman–Kac functional can be interpreted as “creating mass” at a state-dependent exponential rate given by \( q \) up until the stopping time \( \tau \). Accordingly, the expression in \eqref{RW-fkX} can fail to be finite if the domain is sufficiently large. Indeed, the results of Khasminskii [e4] imply that (under similar conditions to those imposed by Dynkin except that \( q\geq 0 \)) the expression in \eqref{RW-fkX} is finite for all \( x\in \overline D \) if and only if there is a continuous function \( u \) that is strictly positive on \( \overline D \) and satisfies \eqref{RW-diffreds}.
It was not until the work of Chung [2] that probabilistic solutions of \eqref{RW-reds} for general (signed) \( q \) became an object of considerable interest. The question of how the oscillations of such a \( q \) affect the behavior of the stopped Feynman–Kac functional \eqref{RW-FKS} is an intriguing one; in particular, killing of mass in some locations may cancel creation of mass in others. The next section describes some of K. L. Chung’s investigations on stopped Feynman–Kac functionals with general \( q \).
Feynman–Kac gauge and positive solutions of the Schrödinger equation
One-dimensional diffusions
Kai Lai Chung initiated his research on stopped Feynman–Kac functionals in [2] by considering a one-dimensional diffusion process (i.e., continuous strong Markov process) \( X \), and the functional \eqref{RW-FKS} with bounded Borel-measurable \( q:\mathbb{R}\to \mathbb{R} \), and \[ \tau=\tau_b\equiv \inf \{ t > 0: X_t =b\} \quad\text{for }b\in \mathbb{R} .\] Assume that \( P^x(\tau_b < \infty ) =1 \) for each \( x\in \mathbb{R} \) and \( b\in \mathbb{R} \), and define \[ v(x,b) = E^x\Bigl[\exp\Bigl(\int_0^{\tau_b} q(X_s) \, ds \Bigr) \Bigr] \quad\text{for } x\in \mathbb{R}, b\in \mathbb{R}. \] Two fundamental properties of \( v \) are that \( 0 < v(x,b) \leq \infty \) for all \( x, b \in \mathbb{R} \), and that \[ v(a,b ) v(b,c) = v(a, c) \] for any \( a < b < c \) or \( c < b < a \) in \( \mathbb{R} \). The latter follows from the strong Markov property. These properties lead to the following:
For fixed \( b\in \mathbb{R} \),
if \( v(x, b) < \infty \) for some \( x < b \), then \( v(x, b) < \infty \) for all \( x < b \),
if \( v(x, b) < \infty \) for some \( x > b \), then \( v(x,b) < \infty \) for all \( x > b \).
Chung [2] introduced the following measures of finiteness of \( v \): \begin{align*} \alpha &= \inf\{ b\in \mathbb{R}: v(x, b) < \infty \ \hbox{ for all } x > b \}, \\ \beta&= \sup \{ b\in \mathbb{R}: v(x,b) < \infty\ \hbox{ for all } x < b\}, \end{align*} and showed the next two results.
\begin{align*} \alpha &= \sup\{ b\in \mathbb{R}: v(x, b) =\infty \ \hbox{ for all } x > b \},\\ \beta &= \inf \{ b\in \mathbb{R}: v(x, b) =\infty \ \hbox{ for all } x < b \}. \end{align*} Furthermore, if \( \beta \in \mathbb{R} \), then \( v(x, \beta ) =\infty \) for all \( x < \beta \), and, if \( \alpha \in \mathbb{R} \), then \( v(x, \alpha ) = \infty \) for all \( x > \alpha \).
The following conditions (a)–(c) are equivalent.
\( \beta =+\infty \).
\( \alpha =-\infty \).
For all \( a, b\in \mathbb{R} \), \( v(a,b)\,v(b,a) \leq 1 \).
Following on from this, Chung and Varadhan [1] proved the next two theorems for \( q:\mathbb{R}\to \mathbb{R} \) bounded and continuous, and \( X \) a one-dimensional Brownian motion.
Fix \( b\in \mathbb{R} \). Suppose that \( v(x,b) < \infty \) for some (and hence all) \( x < b \). Then, \( u(x)=v(x , b) \) is twice continuously differentiable for \( x\in (-\infty , b) \), continuous for \( x\in(-\infty, b] \), and \( u \) satisfies the reduced Schrödinger equation: \begin{equation}\label{RW-schrod} \tfrac{1}{2} u^{\prime\prime}(x) + q u(x) =0 \quad\text{for } x \in (-\infty , b), \end{equation} with the boundary condition \( u(b) = 1 \).
The following conditions are equivalent.
There is a twice continuously differentiable, strictly positive function \( u \) satisfying \eqref{RW-schrod} with \( b=+\infty \).
\( \beta =+\infty \).
\( \alpha =-\infty \).
For all \( a, b\in \mathbb{R} \), \begin{equation} \label{RW-uab} v(a,b)\,v(b,a) \leq 1. \end{equation}
There is some pair of distinct real numbers \( a \), \( b \) such that \eqref{RW-uab} holds.
The equivalence of (a), (b) and (c) is an analogue of Khasminskii’s [e4] results, but with \( d=1 \), \( D=(-\infty, \infty) \) and \( q \) allowed to change sign. Further discussion of this one-dimensional case for Brownian motion can be found in Chapter 9 of the book by Chung and Zhao [6].
Multidimensional Brownian motion
Chung soon moved on to consider multidimensional Brownian motions and domains of finite Lebesgue measure in the work [3] with K. Murali Rao. This paper appeared in the proceedings of the first “Seminar on Stochastic Processes”, held at Northwestern University in 1981. This series of annual conferences was initiated by K. L. Chung, E. Çinlar and R. K. Getoor. The “Seminars” have grown in size over the years, but the novel format of a few invited talks, with ample time reserved for less formal presentations and discussions, has persisted and is one of the attractions of these annual meetings held over two and a half days.
The paper [3] was a significant advance. In particular, it contained the first “gauge theorem.” It is stated in its original form below, and then some generalizations are mentioned.
For this, assume that \( B \) is a \( d \)-dimensional Brownian motion (\( d\geq 1 \)), \( P^x \) and \( E^x \) are respectively probability and expectation operators for \( B \) starting from \( x\in \mathbb{R}^d \), \( q:\mathbb{R}^d \to \mathbb{R} \) is a bounded Borel-measurable function, \( D \) is a domain in \( \mathbb{R}^d \) with closure \( \overline D \) and boundary \( \partial D \), \( m \) is \( d \)-dimensional Lebesgue measure, and \( f:\partial D\to \mathbb{R} \) is a bounded Borel-measurable function with \( f\geq 0 \). Let \[ \tau^{}_D =\inf\{ t > 0: B_t \notin D\}, \] be the first exit time of \( B \) from \( D \). Define \begin{equation} \label{RJW-defU} u(D, q, f; x) = E^x \Bigl[ \exp\Bigl( \int_0^{\tau^{}_D } q(B_s)\,ds\Bigr) f(B_{\tau^{}_D}) ; \tau^{}_D < \infty\Bigr] \quad\text{for } x\in \overline D. \end{equation} When \( m(D) < \infty \), \( P^x(\tau^{}_D < \infty) =1 \) for all \( x\in\overline D \) (see [6], Theorem 1.17) and the qualifier \( \tau^{}_D < \infty \) in \eqref{RJW-defU} may be omitted. The following is Theorem 1.2 in Chung and Rao [3].
Suppose the domain \( D \) satisfies \( m(D) < \infty \). If \( u(D, q, f;\,\cdot\,) \not\equiv \infty \) in \( D \), then it is bounded in \( \overline D \).
While visiting Chung at Stanford University in the early 1980s, Zhongxin Zhao [e7] (see also [6], Theorems 5.19 and 5.20) extended this result by relaxing the assumptions on \( q \) and \( D \). In particular, he showed that the conclusion of Theorem 6 continues to hold if the boundedness condition on \( q \) is relaxed to simply require that \( q:\mathbb{R}^d \to \overline{\mathbb{R}} \) is a Borel-measurable function satisfying \begin{equation} \label{RJW-Kato} \lim_{\alpha \downarrow 0} \Bigl[ \sup_{x\in \mathbb{R}^d} \int_{|y-x|\leq \alpha} \bigl|G(x-y) q(y) \bigr| \,dy \Bigr ] =0, \end{equation} where \( \overline{\mathbb{R}}=[-\infty, \infty] \) is the extended real line and, for \( x\in\mathbb{R}^d \), \[ G (x) = \begin{cases} |x|^{2-d} & \hbox{if } d\geq 3,\\ \ln \dfrac{1}{|x|} & \hbox{if } d=2, \\ |x| & \hbox{if } d=1. \\ \end{cases} \] The set of Borel-measurable functions \( q:\mathbb{R}^d\to\overline{\mathbb{R}} \) satisfying \eqref{RJW-Kato} is called the Kato class (on \( \mathbb{R}^d \)), and is usually denoted by \( K_d \) or \( J \). Various properties of these functions, as well as analytic properties of associated weak solutions of the reduced Schrödinger equation \eqref{RJW-Sch}, are described in an extensive paper of Aizenman and Simon [e6] which appeared shortly after the work [3] of Chung and Rao. The paper [e6] also describes connections between the stopped Feynman–Kac functional and weak solutions of \eqref{RJW-Sch} under a spectral condition (see Theorem 9 below).
Neil Falkner also visited Chung at Stanford in the early 1980s. During this time, Falkner [e8] proved a gauge theorem, when conditioned Brownian motion is used in place of Brownian motion, for bounded Borel-measurable \( q \) and sufficiently smooth bounded domains \( D \). (Zhao subsequently used conditioned Brownian motion in his work [e7].) For more details on the work in [e8] and a discussion of subsequent generalizations, see the article by Michael Cranston in this volume, and Chapter 7 of the book by Chung and Zhao [6].
The function \( u(D, q, 1;\,\cdot\,) \) obtained by setting \( f\equiv 1 \) is called the gauge (function) for \( (D, q) \), and we say that \( (D,q) \) is gaugeable if this function is bounded on \( D \).
Under the assumptions of Theorem 6 and assuming \( (D,q) \) is gaugeable, a second key result in the paper of Chung and Rao [3] provides sufficient conditions for \( u(D, q, f;\,\cdot\,) \) to be a twice continuously differentiable solution of the reduced Schrödinger equation in \( D \), \begin{equation} \label{RJW-Sch} \tfrac{1}{2} \Delta u(x) + q(x) u(x) = 0 \quad \hbox{for } x\in D, \end{equation} with continuous boundary values given by \( f \). As is usual in the theory of elliptic partial differential equations, to ensure two continuous derivatives for \( u \), in dimensions two and higher, one imposes a stronger condition on \( q \) than simple continuity. For example, locally Hölder-continuous functions are often used. For \( d=1 \), let \( {\mathcal C}_1(D) \) denote the set of bounded continuous functions \( h: D \to \mathbb{R} \) and, for \( d\geq 2 \), let \( {\mathcal C}_d(D) \) denote the set of bounded continuous functions \( h: D\to \mathbb{R} \) such that, for each compact set \( K\subset D \), there are strictly positive, finite constants \( \alpha, M \) such that \[ |h(x) -h(y) | \leq M |x-y|^\alpha \quad\text{for all } x, y \in K. \] The following theorem is proved in Chung and Rao [3] for \( d\geq2 \). For \( d=1 \), they impose local Hölder continuity on \( q \) to obtain the result, but this condition can be relaxed by invoking a suitable analytic lemma for a Green potential, as was shown in Chung’s book Lectures from Markov Processes to Brownian Motion [4] (see Proposition 4 of Section 4.7). Note for this that, for \( d=1 \), \( m(D) < \infty \) implies that \( D \) is a bounded interval.
Let \( D \) be a regular domain in \( \mathbb{R}^d \) satisfying \( m(D) < \infty \). Suppose that \( q\in \mathcal C_d(D) \) and \( f:\partial D\to\mathbb{R} \) is bounded and continuous. Assume that \( (D, q) \) is gaugeable, that is, \( u(D, q, 1; \,\cdot\,) \not=\infty \) in \( D \). Then, \( u=u(D, q, f;\,\cdot\,) \) defined by \eqref{RJW-defU} on \( \overline D \) is twice continuously differentiable in \( D \), continuous and bounded on \( \overline D \), satisfies \eqref{RJW-Sch} in \( D \), and has \( u=f \) on \( \partial D \). Furthermore, \( u \) is the unique twice continuously differentiable solution of \eqref{RJW-Sch} that is continuous and bounded on \( \overline D \) and agrees with \( f \) on the boundary \( \partial D \).
This theorem has been generalized to situations where \( q \) is a Kato-class function, and \eqref{RJW-Sch} is interpreted in the weak sense of partial differential equation theory (see [6], Section 4.4)
Note that under the assumptions of the theorem above, if \( f\geq 0 \), then \( u(D, q, f;\,\cdot\,) \) is a non-negative solution of \eqref{RJW-Sch} and, if \( f > 0 \) on \( \partial D \), then \( u(D, q, f ;\,\cdot\,) > 0 \) on \( \overline D \). One naturally expects there to be some relation between the existence of such positive solutions of \eqref{RJW-Sch} and the sign of \[ \lambda (D,q) = \sup_\varphi \Bigl[ \int_D \bigl\{ -\tfrac{1}{2}|\nabla \varphi(x) |^2 + q(x) \varphi(x)^2\bigr\}\,dx \Bigr] , \] where the supremum is over all \( \varphi :D \to \mathbb{R} \) such that \( \varphi \) is infinitely continuously differentiable in \( D \), has compact support in \( D \), and satisfies \( \int_D\varphi(x)^2 dx =1 \). The quantity \( \lambda (D, q) \) is the supremum of the spectrum of the operator \( \frac{1}{2} \Delta + q \) on \( L^2(D) \) (see [6], Proposition 3.29). Indeed, there is a sharp relationship provided by the following theorem (see Theorem 4.19 of [6] for a proof).
Let \( D \) be a domain in \( \mathbb{R}^d \) satisfying \( m(D) < \infty \) and \( q \) be a Kato-class function. Then, \( (D, q) \) is gaugeable if and only if \( \lambda(D, q) < 0 \).
For bounded domains \( D \), Aizenman and Simon proved in Theorem A.4.1 of [e6] the “if” part of this theorem when \( \lambda (D, q) < 0 \), \( u(D, q, f;\,\cdot\,) \) is a weak solution of \eqref{RJW-Sch}, its boundary values are given by \( f \), and they are assumed continuous if \( f \) is continuous and \( D \) is regular.
The work of Chung and Rao [3] was the seed for much subsequent work on connections between the probabilistically defined quantity \eqref{RJW-defU}, gauges, and solutions of the reduced Schrödinger equation \eqref{RJW-Sch}. Besides continuing his own work in this area, in the 1980s Chung had two students, Elton P. Hsu [5], [e10], [e11] and Vassilus Papanicolaou [e12], who worked on probabilistic representations for other boundary-value problems associated with the reduced Schrödinger equation. A conjecture of Chung on equivalent conditions for finiteness of the gauge in terms of finiteness of \( u(D,q, f; \,\cdot\,) \), when \( f \) is a non-negative function that is positive only on a suitable subset of the boundary, stimulated my work [e9] (as a student of Chung) and then Neil Falkner’s [e8] (as a visitor at Stanford). Falkner’s work used conditioned Brownian motion, which became an object of intense interest, in its own right and for its connections with gauge theorems. For more on this fascinating subject, see the accompanying article by Michael Cranston. Other generalizations have also occurred, especially ones involving more general Markov processes than Brownian motion. The works related to this are too numerous to mention here.
Finally, on a personal note, I would like to thank Kai Lai Chung for the pleasure of our collaborations and for the many lively discussions I have enjoyed with him over the years.