# Ponderomotive force

Classical motion of a trapped ion in a radiofrequency (rf) quadrupole (Paul) trap. A quadrupole electric field is displayed for reference, which oscillates at a given frequency ${\displaystyle \omega }$. The blue line represents the ion path in the transversal (or radial) direction of a linear trap, while the orange line is the secular (slow) motion resulting from the ponderomotive force due to the electric field onto the ion. Micromotion is the fast oscillatory motion around the secular motion [1]

In physics, a ponderomotive force is a nonlinear force that a charged particle experiences in an inhomogeneous oscillating electromagnetic field. It causes the particle to move towards the area of the weaker field strength, rather than oscillating around an initial point as happens in a homogeneous field. This occurs because the particle sees a greater magnitude of force during the half of the oscillation period while it is in the area with the stronger field. The net force during its period in the weaker area in the second half of the oscillation does not offset the net force of the first half, and so over a complete cycle this makes the particle move towards the area of lesser force.

The ponderomotive force Fp is expressed by

${\displaystyle \mathbf {F} _{\text{p}}=}$${\displaystyle -{\frac {e^{2}}{4m\omega ^{2}}}}$${\displaystyle \nabla }$${\displaystyle (E^{2})}$

which has units of newtons (in SI units) and where e is the electrical charge of the particle, m is its mass, ω is the angular frequency of oscillation of the field, and E is the amplitude of the electric field. At low enough amplitudes the magnetic field exerts very little force.

This equation means that a charged particle in an inhomogeneous oscillating field not only oscillates at the frequency of ω of the field, but is also accelerated by Fp toward the weak field direction. This is a rare case where the sign of the charge on the particle does not change the direction of the force ((-e)2=(+e)2).

## Derivation

The derivation of the ponderomotive force expression proceeds as follows.

Consider a particle under the action of a non-uniform electric field oscillating at frequency ${\displaystyle \omega }$ in the x-direction. The equation of motion is given by:

${\displaystyle {\ddot {x}}=g(x)\cos(\omega t),}$

neglecting the effect of the associated oscillating magnetic field.

If the length scale of variation of ${\displaystyle g(x)}$ is large enough, then the particle trajectory can be divided into a slow time motion and a fast time motion:[2]

${\displaystyle x=x_{0}+x_{1}}$

where ${\displaystyle x_{0}}$ is the slow drift motion and ${\displaystyle x_{1}}$ represents fast oscillations. Now, let us also assume that ${\displaystyle x_{1}\ll x_{0}}$. Under this assumption, we can use Taylor expansion on the force equation about ${\displaystyle x_{0}}$, to get:

${\displaystyle {\ddot {x}}_{0}+{\ddot {x}}_{1}=\left[g(x_{0})+x_{1}g'(x_{0})\right]\cos(\omega t)}$
${\displaystyle {\ddot {x}}_{0}\ll {\ddot {x}}_{1}}$, and because ${\displaystyle x_{1}}$ is small, ${\displaystyle g(x_{0})\gg x_{1}g'(x_{0})}$, so
${\displaystyle {\ddot {x}}_{1}=g(x_{0})\cos(\omega t)}$

On the time scale on which ${\displaystyle x_{1}}$ oscillates, ${\displaystyle x_{0}}$ is essentially a constant. Thus, the above can be integrated to get:

${\displaystyle x_{1}=-{\frac {g(x_{0})}{\omega ^{2}}}\cos(\omega t)}$

Substituting this in the force equation and averaging over the ${\displaystyle 2\pi /\omega }$ timescale, we get,

${\displaystyle {\ddot {x}}_{0}=-{\frac {g(x_{0})g'(x_{0})}{2\omega ^{2}}}}$
${\displaystyle \Rightarrow {\ddot {x}}_{0}=-{\frac {1}{4\omega ^{2}}}\left.{\frac {d}{dx}}\left[g(x)^{2}\right]\right|_{x=x_{0}}}$

Thus, we have obtained an expression for the drift motion of a charged particle under the effect of a non-uniform oscillating field.

## Time averaged density

Instead of a single charged particle, there could be a gas of charged particles confined by the action of such a force. Such a gas of charged particles is called plasma. The distribution function and density of the plasma will fluctuate at the applied oscillating frequency and to obtain an exact solution, we need to solve the Vlasov Equation. But, it is usually assumed that the time averaged density of the plasma can be directly obtained from the expression for the force expression for the drift motion of individual charged particles:[3]

${\displaystyle {\bar {n}}(x)=n_{0}\exp \left[-{\frac {e}{\kappa T}}\Phi _{\text{P}}(x)\right]}$

where ${\displaystyle \Phi _{\text{P}}}$ is the ponderomotive potential and is given by

${\displaystyle \Phi _{\text{P}}(x)={\frac {m}{4\omega ^{2}}}\left[g(x)\right]^{2}}$

## Generalized ponderomotive force

Instead of just an oscillating field, a permanent field could also be present. In such a situation, the force equation of a charged particle becomes:

${\displaystyle {\ddot {x}}=h(x)+g(x)\cos(\omega t)}$

To solve the above equation, we can make a similar assumption as we did for the case when ${\displaystyle h(x)=0}$. This gives a generalized expression for the drift motion of the particle:

${\displaystyle {\ddot {x}}_{0}=h(x_{0})-{\frac {g(x_{0})g'(x_{0})}{2\omega ^{2}}}}$

## Applications

The idea of a ponderomotive description of particles under the action of a time-varying field has applications in areas like:

The ponderomotive force also plays an important role in laser induced plasmas as a major density lowering factor.

Often, however, the assumed slow-time independency of ${\displaystyle \Phi _{P}}$ is too restrictive, an example being the ultra-short, intense laser pulse-plasma(target) interaction. Here a new ponderomotive effect comes into play, the ponderomotive memory effect.[4] The result is a weakening of the ponderomotive force and the generation of wake fields and ponderomotive streamers.[5] In this case the fast-time averaged density becomes for a Maxwellian plasma: ${\displaystyle {\bar {n}}(x,t)=n_{0}e^{-\Psi }[1+{\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{+\infty }dve^{-v^{2}/2}M(x,v,t)]}$, where ${\displaystyle M(x,v,t):=\int _{-\infty }^{t}d\tau \partial _{\tau }\Psi (x-v(t-\tau ),\tau )}$ and ${\displaystyle \Psi (x,t):={\frac {e}{\kappa T}}\Phi _{P}(x,t)}$.

## References

General
• Schmidt, George (1979). Physics of High Temperature Plasmas, second edition. Academic Press. p. 47. ISBN 978-0-12-626660-3.
Citations
1. ^ D. J. Berkeland; J. D. Miller; J. C. Bergquist; W. M. Itano; D. J. Wineland (1998). "Minimization of ion micromotion in a Paul trap". Journal of Applied Physics. American Institute of Physics. 83 (10): 5025. doi:10.1063/1.367318.
2. ^ Introduction to Plasma Theory, second edition, by Nicholson, Dwight R., Wiley Publications (1983), ISBN 0-471-09045-X
3. ^ V. B. Krapchev, Kinetic Theory of the Ponderomotive Effects in a Plasma, Phys. Rev. Lett. 42, 497 (1979), http://prola.aps.org/abstract/PRL/v42/i8/p497_1
4. ^ H. Schamel and Ch. Sack,"Existence of a Time-dependent Heat Flux-related Ponderomotive Effect",Phys. Fluids 23,1532(1980), doi:10.1063/1.863165
5. ^ U. Wolf and H. Schamel,"Wake-field Generation by the Ponderomotive Memory Effect", Phys. Rev.E 56,4656(1997), doi:10.1103/PhysRevE.56.4656