Skip to main content
Physics LibreTexts

5.15: Poisson’s and Laplace’s Equations

  • Page ID
    97136
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    The electric scalar potential field \(V({\bf r})\), defined in Section 5.12, is useful for a number of reasons including the ability to conveniently compute potential differences (i.e., \(V_{21} = V({\bf r}_2) - V({\bf r}_1)\)) and the ability to conveniently determine the electric field by taking the gradient (i.e., \({\bf E}=-\nabla V\)). One way to obtain \(V({\bf r})\) is by integration over the source charge distribution, as described in Section 5.13. This method is awkward in the presence of material interfaces, which impose boundary conditions on the solutions that must be satisfied simultaneously. For example, the electric potential on a perfectly conducting surface is constant1 – a constraint which is not taken into account in any of the expressions in Section 5.13 (this fact is probably already known to the reader from past study of elementary circuit theory; however, this is established in the context of electromagnetics in Section 5.19.)

    In this section, we develop an alternative approach to calculating \(V({\bf r})\) that accommodates these boundary conditions, and thereby facilitates the analysis of the scalar potential field in the vicinity of structures and spatially-varying material properties. This alternative approach is based on Poisson’s Equation, which we now derive.

    We begin with the differential form of Gauss’ Law (Section 5.7):

    \[\nabla \cdot {\bf D} = \rho_v \nonumber \]

    Using the relationship \({\bf D}=\epsilon{\bf E}\) (and keeping in mind our standard assumptions about material properties, summarized in Section 2.8) we obtain

    \[\nabla \cdot {\bf E} = \frac{\rho_v}{\epsilon} \nonumber \]

    Next, we apply the relationship (Section 5.14):

    \[{\bf E} = -\nabla V \nonumber \] yielding \[\nabla \cdot \nabla V = - \frac{\rho_v}{\epsilon} \nonumber \]

    This is Poisson’s Equation, but it is not in the form in which it is commonly employed. To obtain the alternative form, consider the operator \(\nabla \cdot \nabla\) in Cartesian coordinates:

    \[\begin{aligned} \nabla \cdot \nabla & = \left[ \frac { \partial } { \partial x } \hat { \mathbf { x } } + \frac { \partial } { \partial y } \hat { \mathbf { y } } + \frac { \partial } { \partial z } \hat { \mathbf { z } } \right] \cdot \left[ \frac { \partial } { \partial x } \hat { \mathbf { x } } + \frac { \partial } { \partial y } \hat { \mathbf { y } } + \frac { \partial } { \partial z } \hat { \mathbf { z } } \right] \\[4pt] & = \frac { \partial ^ { 2 } } { \partial x ^ { 2 } } + \frac { \partial ^ { 2 } } { \partial y ^ { 2 } } + \frac { \partial ^ { 2 } } { \partial z ^ { 2 } } \\[4pt] &= \nabla ^ { 2 } \end{aligned} \nonumber \]

    i.e., the operator \(\nabla \cdot \nabla\) is identically the Laplacian operator \(\nabla^2\) (Section 4.10). Furthermore, this is true regardless of the coordinate system employed. Thus, we obtain the following form of Poisson’s Equation:

    \[\boxed{ \nabla^2 V = - \frac{\rho_v}{\epsilon} } \label{m0067_ePoisson} \]

    Poisson’s Equation (Equation \ref{m0067_ePoisson}) states that the Laplacian of the electric potential field is equal to the volume charge density divided by the permittivity, with a change of sign.

    Note that Poisson’s Equation is a partial differential equation, and therefore can be solved using well-known techniques already established for such equations. In fact, Poisson’s Equation is an inhomogeneous differential equation, with the inhomogeneous part \(-\rho_v/\epsilon\) representing the source of the field. In the presence of material structure, we identify the relevant boundary conditions at the interfaces between materials, and the task of finding \(V({\bf r})\) is reduced to the purely mathematical task of solving the associated boundary value problem (see “Additional Reading” at the end of this section). This approach is particularly effective when one of the materials is a perfect conductor or can be modeled as such a material. This is because – as noted at the beginning of this section – the electric potential at all points on the surface of a perfect conductor must be equal, resulting in a particularly simple boundary condition.

    In many other applications, the charge responsible for the electric field lies outside the domain of the problem; i.e., we have non-zero electric field (hence, potentially non-zero electric potential) in a region that is free of charge. In this case, Poisson’s Equation simplifies to Laplace’s Equation:

    \[\boxed{ \nabla^2 V = 0 ~~\mbox{(source-free region)} } \label{m0067_eLaplace} \]

    Laplace’s Equation (Equation \ref{m0067_eLaplace}) states that the Laplacian of the electric potential field is zero in a source-free region.

    Like Poisson’s Equation, Laplace’s Equation, combined with the relevant boundary conditions, can be used to solve for \(V({\bf r})\), but only in regions that contain no charge.


    This page titled 5.15: Poisson’s and Laplace’s Equations is shared under a CC BY-SA 4.0 license and was authored, remixed, and/or curated by Steven W. Ellingson (Virginia Tech Libraries' Open Education Initiative) .