Visualize charged particle trapping and escape in an oscillating electric field (Paul trap) in real time. Adjust AC voltage, frequency, and DC offset to explore Mathieu equation stability.
Presets
Field Parameters
Particle Parameters
Controls
Show electric field
Show particle trails
Run simulation
Statistics
Results
0
Trapped
0
Escaped
0.00
Avg KE
0.0
Time [s]
Mathieu Stability Diagram (bottom-right)
Check whether the current a and q parameters fall inside the stable region (blue area).
Trap
Paul Trap Physics:
Electric potential φ(x,y,t) = [U − V₀cos(Ωt)](x²−y²)/(2r₀²).
Field components Eₓ = −(U−V₀cosΩt)x/r₀², Eᵧ = +(U−V₀cosΩt)y/r₀².
Substituting into the equation of motion yields the Mathieu equation: d²u/dτ² + (a − 2q·cos2τ)u = 0.
Here a = 4qU/(mΩ²r₀²) and q = 2qV₀/(mΩ²r₀²).
Particles are trapped when (a, q) lies inside the first stability region (Nobel Prize 1989).
What exactly is a Paul trap? It sounds like something from sci-fi, but you're saying it can hold a single particle in mid-air?
🎓
Basically, it's a device that uses oscillating electric fields to trap charged particles, like ions, in a small region of space. In practice, it doesn't use physical walls. Instead, the rapidly changing electric field creates a dynamic "potential well" that the particle can't escape. Try moving the RF Voltage slider in the simulator above—you'll see how a stronger field creates a tighter, more stable trap.
🙋
Wait, really? How can an oscillating field trap something? Wouldn't it just push the particle back and forth until it flies out?
🎓
Great intuition! That's the clever part. For a static field, you're right—it would just push the particle away. But by switching the field direction millions of times per second, the particle experiences a net restoring force toward the center, like a ball in a spinning bowl. A common case is trapping ions for quantum computing. In the simulator, adjust the RF Frequency to see how changing the speed of oscillation affects the particle's stability.
🙋
So what's the DC Offset slider for? If the AC field is doing the trapping, why do we need a DC component?
🎓
Excellent question. The DC field creates a static bias that shapes the trap. It can make the trap more stable in one direction but less stable in another. For instance, in mass spectrometers, this is used to selectively eject ions based on their mass. In the simulator, set the DC Offset to a positive value and watch—you'll see the particle's motion become unstable in one plane, causing it to escape. It's a perfect demo of tuning trap stability.
Physical Model & Key Equations
The motion of a charged particle in a Paul trap is governed by the Mathieu equation, which comes from applying Newton's second law to the forces from the oscillating quadrupole electric field.
Here, $u$ represents the particle's position in either the x or y direction. The dimensionless parameters $a_u$ (related to the DC voltage) and $q_u$ (related to the AC voltage and frequency) determine whether the solution is stable (particle trapped) or unstable (particle escapes). Time is scaled to the RF drive frequency.
The stability parameters $a$ and $q$ are what you directly control with the simulator sliders. They are defined in terms of the trap's physical parameters.
Where $e$ and $m$ are the particle's charge and mass, $U_{dc}$ is the DC offset, $V_{rf}$ is the RF voltage amplitude, $\Omega$ is the RF angular frequency, and $r_0$ is a characteristic trap dimension. The particle is trapped only for specific ($a, q$) regions in "stability diagrams"—which is what you're exploring when you adjust the sliders.
Frequently Asked Questions
First, try lowering the AC voltage (q value). If it is too high, it becomes unstable. Next, set the DC voltage (a value) close to 0, and adjust the frequency to enter the stability region of the Mathieu diagram (the triangular region in the a-q plane).
The stability region is the state where the particle trajectory does not diverge and remains near the trap center. Observe whether the particle is trapped for a long time (e.g., more than 100 cycles) while changing parameters. It is efficient to refer to the theoretical curves of the Mathieu diagram.
It is useful for exploring the optimal combination of AC voltage, frequency, and DC component before experiments. In particular, since the stability region of the Mathieu equation can be visually confirmed, it improves the efficiency of parameter adjustment for trapping particles in actual equipment.
a represents the strength of static confinement due to DC voltage, and q represents the strength of dynamic confinement due to AC voltage. The orbital stability of the particle is determined by the combination of the two, and trapping is possible only in the region where the solution of the Mathieu equation is bounded (the stability region).
Real-World Applications
Quantum Computing: Paul traps are the workhorse for ion-trap quantum computers. Individual ions (like Ytterbium) are trapped and laser-cooled to near absolute zero, where their internal energy states become qubits. The precise control offered by the trap's electric fields allows for quantum logic operations.
Mass Spectrometry: In a Quadrupole Mass Filter (a linear Paul trap), a specific combination of DC and RF voltages is applied. Only ions with a specific mass-to-charge ratio ($m/e$) have stable trajectories and pass through to the detector; all others collide with the rods and are filtered out.
Precision Frequency Standards (Atomic Clocks): The most accurate clocks in the world use an optical transition in a single trapped ion (like Aluminium). The Paul trap isolates the ion from disturbing environmental effects, allowing its "ticking" frequency to be measured with extraordinary precision.
Fundamental Physics Research: Paul traps are used to study charged particle dynamics, test quantum mechanics, and even simulate complex many-body physics. They provide a pristine, controllable environment to observe phenomena that are difficult to see elsewhere.
Common Misconceptions and Points to Note
First, you might think "the higher the AC voltage, the stronger the trapping", but this is a misconception. While increasing the AC voltage (the q parameter) does strengthen the force pulling the particle back toward the center, it also increases the amplitude of its oscillations. Beyond a certain threshold, the particle's oscillation diverges and it escapes almost instantly. For example, with the DC component (a) at 0, setting q above 0.9 typically makes the particle unstable. Parameter tuning is all about "applying the right amount of force"; the golden rule is to aim for a point inside the stability region.
Next, remember that simulations represent an ideal environment. In actual experiments, poor vacuum conditions can cause ions to lose energy through collisions with residual gas, preventing trapping. Furthermore, minor imperfections in electrode geometry or power supply noise cannot be ignored. Using this tool while imagining the gap between "theoretically it should trap" and "it doesn't work on the actual device" will help you develop practical intuition.
Finally, do not treat the "DC component as an optional extra". The DC voltage (a) plays a crucial role in balancing the trapping strength in the x and y directions. For instance, you'll find that changing a from 0 to a slightly positive value (e.g., 0.1) alters the shape of the stability region and narrows the range of q values that allow trapping. This is precisely the principle of "mass selectivity" used to filter ions of a specific mass. Understand that DC and AC are two wheels of the same cart.
Set RF voltage (V0) between 100–500 V and RF frequency (Omega) from 1–10 MHz using the sliders.
Adjust the DC offset voltage (U) from −50 to +50 V to create axial confinement asymmetry.
Input the charge-to-mass ratio (Q/M) in C/kg; typical values: singly ionized argon = 1.53×10⁷ C/kg, doubly ionized neon = 3.84×10⁷ C/kg.
Click "Start Simulation" to launch particle trajectory visualization in the trap's quadrupole field.
Monitor stability diagram overlay—green zones indicate stable trapping; red indicates escape.
Worked Example
Paul trap confining a single Ar⁺ ion: V0 = 250 V, Omega = 3.2 MHz, U = +10 V, Q/M = 1.53×10⁷ C/kg. The ion oscillates with secular frequency ωsec ≈ 45 kHz in the radial direction and 52 kHz axially. Electrode diameter = 4 mm. Increasing U to +25 V shifts the axial stability boundary; exceeding this destabilizes axial confinement within 15 microseconds. Decreasing V0 to 150 V narrows the stable parameter region, causing radial escape.
Practical Notes
Stability region (Mathieu q-parameter) depends on V0/(Omega²×r₀²): mass spectrometry applications use q ≈ 0.3–0.7 for broad mass range, while ion traps for quantum simulation operate at q ≈ 0.1 for deep confinement.
DC offset U breaks radial-axial symmetry; use small U values (|U| < 30 V) for symmetric trapping of multiple ions.
For heavy particles (lower Q/M), increase V0 proportionally to maintain secular frequency above space-charge collision rates (>10 kHz typical).
Cooling gas damping (He buffer) reduces oscillation amplitude by ~40% per RF cycle in practical systems; pure vacuum simulations neglect this.