# Liénard–Wiechert potential

Liénard–Wiechert potentials describe the classical electromagnetic effect of a moving electric point charge in terms of a vector potential and a scalar potential in the Lorenz gauge. Built directly from Maxwell's equations, these potentials describe the complete, relativistically correct, time-varying electromagnetic field for a point charge in arbitrary motion, but are not corrected for quantum-mechanical effects. Electromagnetic radiation in the form of waves can be obtained from these potentials.

These expressions were developed in part by Alfred-Marie Liénard in 1898 and independently by Emil Wiechert in 1900[1] and continued into the early 1900s.

## Implications

The study of classical electrodynamics was instrumental in Einstein's development of the theory of relativity. Analysis of the motion and propagation of electromagnetic waves led to the special relativity description of space and time. The Liénard–Wiechert formulation is an important launchpad into more complex analysis of relativistic moving particles.

The Liénard–Wiechert description is accurate for a large, independent moving particle, but breaks down at the quantum level.

Quantum mechanics sets important constraints on the ability of a particle to emit radiation. The classical formulation, as laboriously described by these equations, expressly violates experimentally observed phenomena. For example, an electron around an atom does not emit radiation in the pattern predicted by these classical equations. Instead, it is governed by quantized principles regarding its energy state. In the later decades of the twentieth century, quantum electrodynamics helped bring together the radiative behavior with the quantum constraints.

## Universal Speed Limit

The force on a particle at a given location r and time t depends in a complicated way on the position of the source particles at an earlier time tr due to the finite speed, c, at which electromagnetic information travels. A particle on Earth 'sees' a charged particle accelerate on the Moon as this acceleration happened 1.5 seconds ago, and a charged particle's acceleration on the Sun as happened 500 seconds ago. This earlier time in which an event happens such that a particle at location r 'sees' this event at a later time t is called the retarded time, tr. The retarded time varies with position; for example the retarded time at the Moon is 1.5 seconds before the current time and the retarded time on the Sun is 500 s before the current time. The retarded time can be calculated as:

$t_r=t-\frac{R(t_r)}{c}$

where $R(t_r)$ is the distance of the particle from the source at the retarded time. Only electromagnetic wave effects depend fully on the retarded time.

A novel feature in the Liénard–Wiechert potential is seen in the breakup of its terms into two types of field terms (see below), only one of which depends fully on the retarded time. The first of these is the static electric (or magnetic) field term that depends only on the distance to the moving charge, and does not depend on the retarded time at all, if the velocity of the source is constant. The other term is dynamic, in that it requires that the moving charge be accelerating with a component perpendicular to the line connecting the charge and the observer and does not appear unless the source changes velocity. This second term is connected with electromagnetic radiation.

The first term describes near field effects from the charge, and its direction in space is updated with a term that corrects for any constant-velocity motion of the charge on its distant static field, so that the distant static field appears at distance from the charge, with no aberration of light or light-time correction. This term, which corrects for time-retardation delays in the direction of the static field, is required by Lorentz invariance. A charge moving with a constant velocity must appear to a distant observer in exactly the same way as a static charge appears to a moving observer, and in the latter case, the direction of the static field must change instantaneously, with no time-delay. Thus, static fields (the first term) point exactly at the true instantaneous (non-retarded) position of the charged object if its velocity has not changed over the retarded time delay. This is true over any distance separating objects.

The second term, however, which contains information about the acceleration and other unique behavior of the charge that cannot be removed by changing the Lorentz frame (inertial reference frame of the observer), is fully dependent for direction on the time-retarded position of the source. Thus, electromagnetic radiation (described by the second term) always appears to come from the direction to the position of the emitting charge at the retarded time. Only this second term describes information transfer about the behavior of the charge, which transfer occurs (radiates from the charge) at the speed of light. At "far" distances (longer than several wavelengths of radiation), the 1/R dependence of this term makes electromagnetic field effects (the value of this field term) more powerful than "static" field effects, which are described by the 1/R2 potential of the first (static) term and thus decay more rapidly with distance from the charge.

### Existence and uniqueness of the retarded time

#### Existence

The retarded time is not guaranteed to exist in general. For example, if, in a given frame of reference, an electron has just been created, then at this very moment another electron does not yet feel its electromagnetic force at all. However, under certain conditions, there always exists a retarded time. For example, if the source charge has existed for an unlimited amount of time, during which it has always travelled at a speed not exceeding $v_M < c$, then there exists a valid retarded time $t_r$. To see this, consider the function $f(t') = |\mathbf{r} - \mathbf{r}_s(t')| - c(t - t')$. At the present time, $t' = t$, we have $f(t') = |\mathbf{r} - \mathbf{r}_s(t')| - c(t - t) = |\mathbf{r} - \mathbf{r}_s(t')| \geq 0$. The derivative $f'(t')$ is given by

$f'(t') = \frac{\mathbf{r} - \mathbf{r}_s(t_r)}{|\mathbf{r} - \mathbf{r}_s(t_r)|} \cdot (-\mathbf{v}_s(t')) + c \geq c - \left|\frac{\mathbf{r} - \mathbf{r}_s(t_r)}{|\mathbf{r} - \mathbf{r}_s(t_r)|}\right| \, |\mathbf{v}_s(t')| = c - |\mathbf{v}_s(t')| \geq c - v_M > 0$

By the mean value theorem, $f(t - \Delta t) \leq f(t) - f'(t)\Delta t \leq f(t) - (c - v_M)\Delta t$. By making $\Delta t$ sufficiently large, we can force this to be negative, i.e., at some point in the past, $f(t') < 0$. By the intermediate value theorem, there exists an intermediate $t_r$ with $f(t_r) = 0$, the defining equation of the retarded time. Intuitively, as the source charge moves back in time, the cross section of its light cone at present time expands faster than it can recede, so eventually it must reach the point $\mathbf{r}$. Note that this is not necessarily true if the source charge's speed is allowed to be arbitrarily close to $c$, i.e., if for any given speed $v < c$ there was some time in the past when the charge was moving at this speed. In this case the cross section of the light cone at present time approaches the point $\mathbf{r}$ as we travel back in time but does not necessarily ever reach it.

#### Uniqueness

For a given point $(\mathbf{r}, t)$ and trajectory of the point source $\mathbf{r}_s(t')$, there is at most one value of the retarded time $t_r$, i.e., one value $t_r$ such that $|\mathbf{r} - \mathbf{r}_s(t_r)| = c(t - t_r)$. To see this, suppose there are two retarded times $t_1$ and $t_2$, with $t_1 \leq t_2$. Then, $|\mathbf{r} - \mathbf{r}_s(t_1)| = c(t - t_1)$ and $|\mathbf{r} - \mathbf{r}_s(t_2)| = c(t - t_2)$. Subtracting gives $c(t_2 - t_1) = |\mathbf{r} - \mathbf{r}_s(t_1)| - |\mathbf{r} - \mathbf{r}_s(t_2)| \leq |\mathbf{r}_s(t_2) - \mathbf{r}_s(t_1)|$ by the triangle inequality. Unless $t_2 = t_1$, this then implies that the average velocity of the charge between $t_1$ and $t_2$ is $|\mathbf{r}_s(t_2) - \mathbf{r}_s(t_1)|/(t_2 - t_1) \geq c$, which is impossible. The intuitive interpretation is that we can only ever "see" the point source at one location/time at once unless it travels at least at the speed of light to another location. As the source moves forward in time, the cross section of its light cone at present time contracts faster than the source can approach, so it can never intersect the point $\mathbf{r}$ again.

We conclude that, under certain conditions, the retarded time exists and is unique.

## Equations

### Definition of Liénard–Wiechert potentials

The Liénard–Wiechert potentials $\varphi$ (scalar potential field) and $\mathbf{A}$ (vector potential field) are for a source point charge $q$ at position $\mathbf{r}_s$ traveling with velocity $\mathbf{v}_s$:

$\varphi(\mathbf{r}, t) = \frac{1}{4 \pi \epsilon_0} \left(\frac{q}{(1 - \mathbf{n} \cdot \boldsymbol{\beta}_s)|\mathbf{r} - \mathbf{r}_s|} \right)_{t_r}$

and

$\mathbf{A}(\mathbf{r},t) = \frac{\mu_0c}{4 \pi} \left(\frac{q \boldsymbol{\beta}_s}{(1 - \mathbf{n} \cdot \boldsymbol{\beta}_s)|\mathbf{r} - \mathbf{r}_s|} \right)_{t_r} = \frac{\boldsymbol{\beta}_s(t_r)}{c} \varphi(\mathbf{r}, t)$

where $\boldsymbol{\beta}_s(t) = \frac{\mathbf{v}_s(t)}{c}$ and $\mathbf{n} = \frac{\mathbf{r} - \mathbf{r}_s}{|\mathbf{r} - \mathbf{r}_s|}$.

### Corresponding values of electric and magnetic fields

We can calculate the electric and magnetic fields directly from the potentials using the definitions:

$\mathbf{E} = - \nabla \varphi - \dfrac{\partial \mathbf{A}}{\partial t}$ and $\mathbf{B} = \nabla \times \mathbf{A}$

The calculation is nontrivial and requires a number of steps. The electric and magnetic fields are (in non-covariant form):

$\mathbf{E}(\mathbf{r}, t) = \frac{1}{4 \pi \epsilon_0} \left(\frac{q(\mathbf{n} - \boldsymbol{\beta})}{\gamma^2 (1 - \mathbf{n} \cdot \boldsymbol{\beta})^3 |\mathbf{r} - \mathbf{r}_s|^2} + \frac{q \mathbf{n} \times \big((\mathbf{n} - \boldsymbol{\beta}) \times \dot{\boldsymbol{\beta}}\big)}{c(1 - \mathbf{n} \cdot \boldsymbol{\beta})^3 |\mathbf{r} - \mathbf{r}_s|} \right)_{t_r}$

and

$\mathbf{B}(\mathbf{r}, t) = \frac{\mu_0}{4 \pi} \left(\frac{q c(\boldsymbol{\beta} \times \mathbf{n})}{\gamma^2 (1-\mathbf{n} \cdot \boldsymbol{\beta})^3 |\mathbf{r} - \mathbf{r}_s|^2} + \frac{q \mathbf{n} \times \Big(\mathbf{n} \times \big((\mathbf{n} - \boldsymbol{\beta}) \times \dot{\boldsymbol{\beta}}\big) \Big)}{(1 - \mathbf{n} \cdot \boldsymbol{\beta})^3 |\mathbf{r} - \mathbf{r}_s|} \right)_{t_r} = \frac{\mathbf{n}(t_r)}{c} \times \mathbf{E}(\mathbf{r}, t)$

where $\boldsymbol{\beta}(t) = \frac{\mathbf{v}_s(t)}{c}$, $\mathbf{n}(t) = \frac{\mathbf{r} - \mathbf{r}_s(t)}{|\mathbf{r} - \mathbf{r}_s(t)|}$ and $\gamma(t) = \frac{1}{\sqrt{1 - |\boldsymbol{\beta}(t)|^2}}$ (the Lorentz factor).

Note that the $\mathbf{n} - \boldsymbol{\beta}$ part of the first term updates the direction of the field toward the instantaneous position of the charge, if it continues to move with constant velocity $c\boldsymbol{\beta}$. This term is connected with the "static" part of the electromagnetic field of the charge.

The second term, which is connected with electromagnetic radiation by the moving charge, requires charge acceleration $\dot{\boldsymbol{\beta}}$ and if this is zero, the value of this term is zero, and the charge does not radiate (emit electromagnetic radiation). This term requires additionally that a component of the charge acceleration be in a direction transverse to the line which connects the charge $q$ and the observer of the field $\mathbf{E}(\mathbf{r}, t)$. The direction of the field associated with this radiative term is toward the fully time-retarded position of the charge (i.e. where the charge was when it was accelerated).

## Derivation

In the case that there are no boundaries surrounding the sources, the retarded solutions for the scalar and vector potentials (SI units) of the nonhomogeneous wave equations with sources given by the charge and current densities $\rho(\mathbf{r}, t)$ and $\mathbf{J}(\mathbf{r}, t)$ are in the Lorenz gauge (see Nonhomogeneous electromagnetic wave equation)

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi \epsilon_0}\int \frac{\rho(\mathbf{r}', t_r')}{|\mathbf{r} - \mathbf{r}'|} d^3\mathbf{r}'$

and

$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \int \frac{\mathbf{J}(\mathbf{r}', t_r')}{|\mathbf{r} - \mathbf{r}'|} d^3\mathbf{r}'$

where $t_r' = t - \frac{1}{c}|\mathbf{r} - \mathbf{r}'|$ is the retarded time.

For a moving point charge whose trajectory is given as a function of time by $\mathbf{r}_s(t')$, the charge and current densities are as follows:

$\rho(\mathbf{r}', t') = q \delta^3(\mathbf{r'} - \mathbf{r}_s(t'))$
$\mathbf{J}(\mathbf{r}', t') = q\mathbf{v}_s(t') \delta^3(\mathbf{r'} - \mathbf{r}_s(t'))$

where $\delta^3$ is the three-dimensional Dirac delta function and $\mathbf{v}_s(t')$ is the velocity of the point charge.

Substituting into the expressions for the potential gives

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi\epsilon_0} \int \frac{q \delta^3(\mathbf{r'} - \mathbf{r}_s(t_r'))}{|\mathbf{r} - \mathbf{r}'|} d^3\mathbf{r}'$
$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \int \frac{q\mathbf{v}_s(t_r') \delta^3(\mathbf{r'} - \mathbf{r}_s(t_r'))}{|\mathbf{r} - \mathbf{r}'|} d^3\mathbf{r}'$

These integrals are difficult to evaluate in their present form, so we will rewrite them by replacing $t_r'$ with $t'$ and integrating over the delta distribution $\delta(t' - t_r')$:

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi \epsilon_0} \iint \frac{q\delta^3(\mathbf{r'} - \mathbf{r}_s(t'))}{|\mathbf{r} - \mathbf{r}'|} \delta(t' - t_r') \, dt' \, d^3\mathbf{r}'$
$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \iint \frac{q\mathbf{v}_s(t') \delta^3(\mathbf{r'} - \mathbf{r}_s(t'))}{|\mathbf{r} - \mathbf{r}'|} \delta(t' - t_r') \, dt' \, d^3\mathbf{r}'$

We exchange the order of integration:

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi \epsilon_0} \iint \frac{\delta(t' - t_r')}{|\mathbf{r} - \mathbf{r}'|} q\delta^3(\mathbf{r'} - \mathbf{r}_s(t')) \, d^3\mathbf{r}' dt'$
$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \iint \frac{\delta(t' - t_r')}{|\mathbf{r} - \mathbf{r}'|} q\mathbf{v}_s(t') \delta^3(\mathbf{r'} - \mathbf{r}_s(t')) \, d^3\mathbf{r}' dt'$

The delta function picks out $\mathbf{r}' = \mathbf{r}_s(t')$ which allows us to perform the inner integration with ease. Note that $t_r'$ is a function of $\mathbf{r}'$, so this integration also fixes $t_r = t_r(\mathbf{r}_s(t'), t')$.

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi \epsilon_0} \int q\frac{\delta(t' - t_r')}{|\mathbf{r} - \mathbf{r}_s(t')|} dt'$
$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \int q\mathbf{v}_s(t') \frac{\delta(t' - t_r')}{|\mathbf{r} - \mathbf{r}_s(t')|} \, dt'$

The retarded time $t_r'$ is a function of the field point $(\mathbf{r}, t)$ and the source trajectory $\mathbf{r}_s(t')$, and hence depends on $t'$. To evaluate this integral, therefore, we need the identity

$\delta(f(t')) = \sum_i \frac{\delta(t' - t_i)}{|f'(t_i)|}$

where each $t_i$ is a zero of $f$. Because there is only one retarded time $t_r$ for any given space-time coordinates $(\mathbf{r}, t)$ and source trajectory $\mathbf{r}_s(t')$, this reduces to:

\begin{align}\delta(t' - t_r') = \frac{\delta(t' - t_r)}{\frac{\partial}{\partial t'}(t' - t_r')|_{t' = t_r}} =& \frac{\delta(t' - t_r)}{\frac{\partial}{\partial t'}(t' - (t - \frac{1}{c}|\mathbf{r} - \mathbf{r}_s(t')|))|_{t' = t_r}} = \frac{\delta(t' - t_r)}{1 +\frac{1}{c}(\mathbf{r} - \mathbf{r}_s(t'))/|\mathbf{r} - \mathbf{r}_s(t')|\cdot (-\mathbf{v}_s(t')) |_{t' = t_r}}=\\ &= \frac{\delta(t' - t_r)}{1 - \boldsymbol{\beta}_s \cdot (\mathbf{r}-\mathbf{r}_s)/|\mathbf{r}-\mathbf{r}_s|}\end{align}

where $\boldsymbol{\beta}_s = \mathbf{v}_s/c$, and $\boldsymbol{\beta}_s$ and $\mathbf{r}_s$ are evaluated at the retarded time, and we have used the identity $|\mathbf{x}|' = \hat{\mathbf{x}} \cdot \mathbf{v}$. Finally, the delta function picks out $t' = t_r$, and

$\varphi(\mathbf{r}, t) = \frac{1}{4\pi \epsilon_0} \left(\frac{q}{|\mathbf{r}-\mathbf{r}_s| (1 - \boldsymbol{\beta}_s \cdot (\mathbf{r}-\mathbf{r}_s)/|\mathbf{r}-\mathbf{r}_s|)}\right)_{t_r} = \frac{1}{4\pi \epsilon_0} \left(\frac{q}{(1-\mathbf{n}\cdot \boldsymbol{\beta}_s)|\mathbf{r}-\mathbf{r}_s|}\right)_{t_r}$
$\mathbf{A}(\mathbf{r}, t) = \frac{\mu_0}{4\pi} \left(\frac{q\mathbf{v}}{|\mathbf{r}-\mathbf{r}_s| (1 - \boldsymbol{\beta}_s \cdot (\mathbf{r}-\mathbf{r}_s)/|\mathbf{r}-\mathbf{r}_s|)}\right)_{t_r} = \frac{\mu_0 c}{4\pi} \left(\frac{q\boldsymbol{\beta}_s}{(1-\mathbf{n}\cdot \boldsymbol{\beta}_s)|\mathbf{r}-\mathbf{r}_s|}\right)_{t_r}$

which are the Liénard–Wiechert potentials.

### Lorenz gauge, electric and magnetic fields

In order to calculate the derivatives of $\varphi$ and $\mathbf{A}$ it is convenient to first compute the derivatives of the retarded time. Taking the derivatives of both sides of its defining equation (remembering that $\mathbf{r_s} = \mathbf{r_s}(t_r)$):

$t_r + \frac{1}{c}|\mathbf{r}-\mathbf{r_s}|= t$

Differentiating with respect to t,

$\frac{d t_r}{d t} + \frac{1}{c}\frac{d t_r}{d t}\frac{d |\mathbf{r}-\mathbf{r_s}|}{d t_r}= 1$
$\frac{d t_r}{d t} \left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right) = 1$
$\frac{d t_r}{d t} = \frac{1}{\left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right)}$

Similarly, Taking the gradient with respect to $\mathbf{r}$ gives

${\boldsymbol \nabla} t_r + \frac{1}{c}{\boldsymbol \nabla} |\mathbf{r}-\mathbf{r_s}| = 0$
${\boldsymbol \nabla} t_r + \frac{1}{c} \left({\boldsymbol \nabla} t_r \frac{d |\mathbf{r}-\mathbf{r_s}|}{d t_r} + \mathbf{n}\right) = 0$
${\boldsymbol \nabla} t_r \left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right) = -\mathbf{n}/c$
${\boldsymbol \nabla} t_r = -\frac{\mathbf{n}/c}{\left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right)}$

It follows that

$\frac{d |\mathbf{r}-\mathbf{r_s}|}{d t} = \frac{d t_r}{d t}\frac{d |\mathbf{r}-\mathbf{r_s}|}{d t_r} = \frac{- \mathbf{n}\cdot{\boldsymbol \beta}_s c}{\left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right)}$
${\boldsymbol \nabla} |\mathbf{r}-\mathbf{r_s}| = {\boldsymbol \nabla} t_r \frac{d |\mathbf{r}-\mathbf{r_s}|}{d t_r} + \mathbf{n} = \frac{\mathbf{n}}{\left(1 - \mathbf{n}\cdot{\boldsymbol \beta}_s\right)}$

These can be used in calculating the derivatives of the vector potential and the resulting expressions are

\begin{align}\frac{d \varphi}{d t} =& -\frac{q}{4\pi\epsilon_0}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^2}\frac{d}{d t}\left[(|\mathbf{r}-\mathbf{r_s}|(1-\mathbf{n}\cdot{\boldsymbol \beta}_s)\right]\\ =& -\frac{q}{4\pi\epsilon_0}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^2}\frac{d}{d t}\left[|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\right]\\ =& -\frac{q c}{4\pi\epsilon_0}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\left[- \mathbf{n}\cdot{\boldsymbol \beta}_s + {\beta_s}^2 - (\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s /c \right]\end{align}
\begin{align}{\boldsymbol \nabla}\cdot\mathbf{A} =& -\frac{q}{4\pi\epsilon_0 c}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^2}\big({\boldsymbol \nabla} \left[\left(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\right)\right]\cdot{\boldsymbol \beta}_s - \left[\left(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\right)\right]{\boldsymbol \nabla}\cdot{\boldsymbol \beta}_s\big)\\ =& - \frac{q}{4\pi\epsilon_0 c}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\cdot\\ &\left[(\mathbf{n}\cdot {\boldsymbol \beta}_s) - {\beta}_s^2(1-\mathbf{n}\cdot {\boldsymbol \beta}_s) - {\beta}_s^2\mathbf{n}\cdot {\boldsymbol \beta}_s + \left((\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s/c\right)(\mathbf{n}\cdot {\boldsymbol \beta}_s) + \big(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\big)(\mathbf{n}\cdot \dot {\boldsymbol \beta}_s/c)\right] \\=&\frac{q}{4\pi\epsilon_0 c}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\left[\beta_s^2 - \mathbf{n}\cdot{\boldsymbol \beta}_s - (\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s/c\right]\end{align}

These show that the Lorenz gauge is satisfied, namely that $\frac{d \varphi}{d t} + c^2 {\boldsymbol \nabla}\cdot\mathbf{A} = 0$.

Similarly one calculates:

${\boldsymbol \nabla}\varphi = -\frac{q}{4\pi\epsilon_0}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\left[\mathbf{n}\left(1-{\beta_s}^2 + (\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s/c\right) - {\boldsymbol \beta}_s(1-\mathbf{n}\cdot{\boldsymbol \beta}_s)\right]$
$\frac{d\mathbf{A}}{dt} = \frac{q}{4\pi\epsilon_0}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\left[{\boldsymbol \beta}_s\left(\mathbf{n}\cdot{\boldsymbol \beta}_s-{\beta_s}^2 + (\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s/c\right) + |\mathbf{r}-\mathbf{r_s}|\dot {\boldsymbol \beta}_s (1-\mathbf{n}\cdot{\boldsymbol \beta}_s)/c\right]$

By noting that for any vectors $\mathbf{u}$, $\mathbf{v}$, $\mathbf{w}$:

$\mathbf{u}\times(\mathbf{v}\times\mathbf{w}) = (\mathbf{u}\cdot\mathbf{w})\mathbf{v}- (\mathbf{u}\cdot \mathbf{v})\mathbf{w}$

The expression for the electric field mentioned above becomes

\begin{align}\mathbf{E}(\mathbf{r}, t) =& \frac{q}{4 \pi \epsilon_0} \frac{1}{|\mathbf{r} - \mathbf{r}_s|^2(1 - \mathbf{n} \cdot {\boldsymbol \beta}_s)^3}\cdot \\ &\left[\left(\mathbf{n} - {\boldsymbol \beta}_s\right)(1-{\beta_s}^2) + |\mathbf{r} - \mathbf{r}_s|(\mathbf{n} \cdot \dot{\boldsymbol \beta}_s/c) (\mathbf{n} - {\boldsymbol \beta}_s) - |\mathbf{r} - \mathbf{r}_s|\big(\mathbf{n} \cdot (\mathbf{n} - {\boldsymbol \beta}_s)\big) \dot{\boldsymbol \beta}_s/c \right]\end{align}

which is easily seen to be equal to $-{\boldsymbol \nabla}\varphi - \frac{d\mathbf{A}}{dt}$

Similarly ${\boldsymbol \nabla}\times\mathbf{A}$ gives the expression of the magnetic field mentioned above:

\begin{align}{\mathbf{B}} =& {\boldsymbol \nabla}\times\mathbf{A} = -\frac{q}{4\pi\epsilon_0 c}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^2}\big({\boldsymbol \nabla} \left[\left(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\right)\right]\times{\boldsymbol \beta}_s - \left[\left(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\right)\right]{\boldsymbol \nabla}\times{\boldsymbol \beta}_s\big)\\ =& - \frac{q}{4\pi\epsilon_0 c}\frac{1}{|\mathbf{r}-\mathbf{r_s}|^2\left(1-\mathbf{n}\cdot{\boldsymbol \beta}_s\right)^3}\cdot\\ &\left[(\mathbf{n}\times {\boldsymbol \beta}_s) - ({\boldsymbol \beta}_s\times {\boldsymbol \beta}_s)(1-\mathbf{n}\cdot {\boldsymbol \beta}_s) - {\beta}_s^2\mathbf{n}\times {\boldsymbol \beta}_s + \left((\mathbf{r}-\mathbf{r_s})\cdot \dot {\boldsymbol \beta}_s/c\right)(\mathbf{n}\times {\boldsymbol \beta}_s) + \big(|\mathbf{r}-\mathbf{r_s}|-(\mathbf{r}-\mathbf{r_s})\cdot{\boldsymbol \beta}_s\big)(\mathbf{n}\times \dot {\boldsymbol \beta}_s/c)\right] \\=& -\frac{q}{4 \pi \epsilon_0 c} \frac{1}{|\mathbf{r} - \mathbf{r}_s|^2(1 - \mathbf{n} \cdot {\boldsymbol \beta}_s)^3}\cdot \\ &\left[\left(\mathbf{n}\times{\boldsymbol \beta}_s\right)(1-{\beta_s}^2) + |\mathbf{r} - \mathbf{r}_s|(\mathbf{n} \cdot \dot{\boldsymbol \beta}_s/c) (\mathbf{n}\times {\boldsymbol \beta}_s) + |\mathbf{r} - \mathbf{r}_s|\big(\mathbf{n} \cdot (\mathbf{n} - {\boldsymbol \beta}_s)\big) \mathbf{n}\times\dot{\boldsymbol \beta}_s/c \right] = \frac{\mathbf{n}}{c}\times\mathbf{E} \end{align}