In this sequence of posts I will present a probabilistic solution of the classical Dirichlet problem using Brownian motion. I will follow the structure of the book Green, Brown and Probability and Kai-Lai Chung with some little changes and somewhat more explanation. We assume knowledge of “advanced calculus”, real variables and some probability theory. No knowledge of Dirichlet problem and Brownian motion is assumed. We hope that the posts are accessible to undergraduates.

**Convention**: For convenience, we will work solely on the Euclidean space . It turns out that the theory is much the same for all spaces with . The theory for is quite different, as we will see.

We begin by recalling the classical Dirichlet problem. Let be a *domain* in . By definition, this means that is a non-empty connected open set in . A real-valued function defined on is said to be *harmonic* on if if the second derivatives of are continuous on (this we denote by writing , and

on . The operator is called the Laplacian.

**Classical Dirichlet Problem:** Let be a continuous function on , the boundary of . Find a function with the following properties:

i) is continuous on .

ii) is harmonic on .

iii) on .

The function is said to be a *solution* of the Dirichlet problem.

**An example in electrostatics.** To put things into perspective, let us describe a situation where this equation arises. Recall that in electrostatics, the force between two point charges and equals

,

where is the distance between the charges and is some constants depending on the unit. (This is called Coulomb’s Law.) The charge (say at ) creates a potential in space, given by

.

Suppose now that we have a distribution of charges on the surface of some conductor, two parallel metal plates. The distribution can be described in terms of a density function , where ranges over the conductor. The charge creates a potential

.

Here the integration is over the conductor (two surfaces in this case). Here is the key fact:

*The potential function is harmonic outside the conductor.*

To understand this at least intuitively, you should check that is harmonic outside the singularity . Then integrating still gives a harmonic function (**Exercise:** Make this precise under appropriate conditions). Suppose that the potential on the plates are fixed. Then the problem of determining the potential between the plates can be modeled as a Dirichlet problem.

Now we have a well posed mathematical problem, namely the Dirichlet problem. (I will omit the word “classical” from now on.) In dealing with a boundary value problem like this, we would like to ask the following questions:

1) Given the domain and a boundary continuous function , can we always solve the Dirichlet problem?

2) If a solution exists, is it unique?

3) Are there ways, or formulas, for computing solutions?

You may already know the answers from a previous course on partial differential equations or private reading. The purpose in this notes is to obtain and explain these results *from the perspective of Brownian motion*. There are very deep connection behind these two seemingly unrelated theories. We prefer not to spell out the relationship at the moment…

This ends the introduction. In the next post we will review Green’s identity and define Green’s function. After that we will give a hint of why Brownian motion is relevant, if at all.