## When solitons collide

**Adapted from notes I copied during a course I attended**

Consider a 1+1 dimensional scalar field where x corresponds to the space dimension, and t corresponds to the time dimension. The Lagrangian and the Hamiltonian energy density for the field are given by,

Here, the Einstein summation convention is used, where repeated indices are summed over. is used to denote the four-vector

The evolution of this scalar field is governed by the Euler-Lagrange equation,

Plugging in Lagrangian for the scalar field I’ve written above, the E-L equation reduces to,

where, is the D’Alembertian. ie,

A specific example, where the solutions of the E-L equation have interesting properties, is given by the potential,

The lowest energy solutions which satisfy the E-L equations are known as the **classical vaccum**. In this specific example, they are given by

You can see that these solutions correspond to the minimum of the potential. Additionally, the total energy is also minimized by these solutions, because the contribution of the space and time derivatives to the Kinetic Energy density term is zero.

Apart from these vacuum solutions, there exist **finite energy** solutions of the E-L equations for this potential. From the criteria that the energy must be finite when the Hamiltonian energy density is integrated over x, such a solution must tend to the classical vacuua solutions at the extremities.

Such a solution is known as the *kink solition*, and is given by,

where,

You can check that this solution satisfies the Euler-Lagrange equation by plugging it in. Note that the solution is time-independent. Furthermore, when , . Therefore, the solution interpolates between the classical vacuua solutions.

To understand why it’s called a kink soliton, let’s plot its functional form, and the corresponding Energy density. For the plots below, I’ve set

As you can see, the solution tends to the vaccum solution (here, a=1) at either end and the energy density appears as a stationary wave packet which is like a kink in a rope.

Now, there are certain global symmetries that Lagrangian has, which are absent in the kink soliton solution. For instance, if we set , the Lagrangian in Eq. (1), and consequently the corresponding action, is invariant. However, the solution , is a **new solution** and interpolates the other way (ie, as . This is called the *anti-kink* solition, and has the same energy as the kink soliton.

There are other symmetries where the Lagrangian is invariant which give rise to new solutions. For example, a spatial translation will give a solution , which shifts the location of the center of the kink. A lorentz transformation will give a new *time dependent solution*, , where . This corresponds to a kink which starts out at the origin and moves to the right with a constant velocity.

Therefore, the most general form of the kink soliton can be written as,

This corresponds to a soliton which is initially centered at x_0 and moves with a velocity v to the right. The anti-kink soliton corresponding to this is simply

Now, consider a box which has both a kink and an anti-kink solition placed away from each other and with velocities in opposite directions. What happens to them as they approach each other? They can’t pass each other, otherwise they’d violate the boundary condition (remember that the kink and anti-kink solitons tend to *different* vacuum solutions at either end). So you’d expect them to bounce off each other.

In fact, we can even numerically test this, because we know that they are governed by the Euler-Lagrange equation. If represents both the kink and anti-kink soliton, then initially, one has

then, the E-L equation can be rewritten as two first order differential equations

where, is the momentum density

To solve the system of equations above, I discretized the spatial derivative using the central difference method, and used a fourth-order Runge Kutta solver with corresponding intial and boundary condition discussed above. You can see how the kink and anti-kink soliton bounce off each other in the video below.

Kinky, ain’t it?

Wikipedia has a gallery of images on solitons for the Sine-Gordon model, where a Sine potential is used instead of the quartic potential. There are also a lot of nice websites on the fantastic properties that such non-linear PDE’s exhibit.

__See Also__

Sine-Gordon Solitons and Soliton Collisions:

Pendulum Model

Soliton-Lab Art Gallery

## Leave a Reply