Schwinger Pair Production Rate Calculator

What this calculator estimates

The Schwinger effect is one of the striking predictions of quantum electrodynamics: a strong enough electric field can destabilize the vacuum and create charged particle pairs. In the most familiar case, those particles are electron-positron pairs. This calculator evaluates the leading Schwinger pair-production rate per unit volume per unit time for an idealized constant, spatially uniform electric field. In practical terms, it helps you answer a focused question: if a particle has mass m and charge q, and the field strength is E, how large is the tunnelling rate predicted by the leading term of the Schwinger formula?

That focus matters because the result is not a general laboratory yield estimate. The output is a local rate density, written in cubic metres and seconds. If you later want an approximate total number of pairs in a finite setup, you would still need an effective interaction volume, an effective duration, and a judgment about whether the constant-field approximation is meaningful over that region. The calculator intentionally does the narrow physics step first: it turns your chosen mass, charge, and field into a rate, a critical field scale, and the dimensionless ratio E/E_c.

The default mass and charge values correspond to an electron. The electric-field box is left for you to choose on purpose, because the answer is extraordinarily sensitive to field strength. At low fields the rate is so suppressed that it is effectively zero for numerical work; near the critical field it rises explosively. That sharp transition is exactly why people use a calculator here instead of mental arithmetic.

Inputs and how to choose them

The first input is Particle Mass m (kg). Enter the particle rest mass in SI kilograms. For the standard electron case, the prefilled value is 9.10938356 × 10⁻³¹ kg. If you are exploring hypothetical charged particles or comparing how the rate changes with mass, remember that the exponent contains m². Doubling the mass therefore pushes the critical field up strongly and makes pair creation much harder.

The second input is Charge q (C). The script uses the magnitude of the charge, so you may enter either a signed charge or a positive magnitude; the production rate depends on |q| rather than the sign. For an electron, the magnitude is 1.602176634 × 10⁻¹⁹ C. Larger charge magnitude lowers the critical field and tends to increase the rate, because the electric field couples more strongly to the particle.

The third input is Electric Field E (V/m). This must be positive in the model used here. The field strength does two jobs at once. It appears in the prefactor as E², which would already make the rate grow with field, but the much more dramatic effect is hidden in the exponential term. When E is well below the critical field, the exponent is so negative that the result is essentially zero. Once E approaches E_c, the suppression weakens quickly and the rate can become enormous.

A good habit is to think in orders of magnitude before you press the button. For electrons, the critical field is roughly 1.3 × 10¹⁸ V/m. If you enter something like 10¹⁴ V/m, you should expect the calculator to report an almost vanishing rate. If you enter something on the order of 10¹⁸ V/m, you should expect a dramatic jump. That quick expectation check makes it easier to catch unit slips such as confusing kV/m with V/m or omitting a power of ten.

The formula used by the calculator

The page computes the leading term of the Schwinger production rate in SI units. Written explicitly, the calculator uses

It also reports the critical field

Those two expressions explain nearly all of the behaviour you will see on the page. The prefactor grows only polynomially, while the exponential term depends on the ratio E_c/E. That is why the output can look almost binary when you scan across field strengths: far below critical, production is crushed; near critical, production turns on rapidly.

In more complete treatments, the Schwinger rate appears as a sum over terms with integer index n. The calculator deliberately keeps only the leading n = 1 contribution, which is the dominant piece in the simple constant-field picture and matches the JavaScript on this page. That makes the tool useful for first-pass estimates, while still leaving room for more detailed analysis if you care about higher-order terms, pulse shape, magnetic fields, or backreaction.

Worked example with electron values

Suppose you keep the default electron mass and charge and enter an electric field of 1.0 × 10¹⁸ V/m. For electrons, the critical field is about 1.3 × 10¹⁸ V/m, so the ratio E/E_c is a bit below one. That already tells you the exponent will no longer be catastrophically large in magnitude. Instead of an astronomically tiny tunnelling probability, you are in a regime where the suppression has weakened enough for the prefactor to matter.

Using those values, the exponent is approximately −πE_c/E ≈ −4.1, so the exponential factor is only about 0.016. The prefactor is still extremely large in SI units, which pushes the final leading-term estimate into an enormous rate density, on the order of 10⁵⁴ m⁻³s⁻¹. That number should not be read as a direct count in a real apparatus; it is a reminder that once you enter the near-critical regime, the idealized vacuum becomes violently unstable in the model.

Now compare that with a field of 1.0 × 10¹⁶ V/m. The ratio E/E_c drops to roughly 0.0076, so the exponent becomes about −414. The exponential factor is then so tiny that ordinary floating-point arithmetic underflows to zero. When the calculator says the rate is numerically suppressed for the chosen inputs, it is telling you that the theoretical expression is far below the precision of the computation and effectively irrelevant on any practical scale described by this idealized model.

How to read the results panel

After you submit the form, the page returns three values. The first is Γ, the pair-production rate density in m⁻³s⁻¹. The second is the critical field E_c in V/m for the particle you specified. The third is the field ratio E/E_c, a dimensionless number that often gives the quickest intuition. If the ratio is tiny, the rate should be tiny. If the ratio is approaching one, the exponential suppression is weakening rapidly.

For many inputs the rate shown by the script will be either extraordinarily small or extraordinarily large. That is not a bug; it reflects the stiffness of the underlying exponential. A result displayed as approximately zero usually means the chosen field is deep in the suppressed regime. A huge result means the idealized constant-field model is predicting strong vacuum breakdown, at which point omitted effects such as field depletion, pulse structure, and interactions with the environment may become important. In other words, the calculator is excellent for scale awareness, but it is not a full experimental simulation.

Quick comparison table for electron-positron production

The table below uses the electron mass and charge to show how sensitive the leading estimate is to the field strength. The intent is not to replace the form, but to give you intuition before you test your own numbers.

If your own scenario differs from the electron case, the same reading strategy still works. First inspect the critical field for the chosen mass and charge, then compare your applied field to it. The ratio is the bridge between raw input numbers and physical intuition.

Assumptions and limitations that matter

This calculator is intentionally narrow. It assumes a constant, spatially uniform electric field and uses the leading Schwinger term. Real strong-field environments can be pulsed, inhomogeneous, mixed with magnetic fields, or modified by plasma dynamics and backreaction. The calculator also treats the vacuum response locally, which is useful for order-of-magnitude reasoning but does not by itself predict a detector count, beam profile, or spectrum.

The sign of the charge does not change the result because the magnitude |q| is what enters the tunnelling rate. The mass and charge are assumed to be fixed properties of a single particle species, not environment-dependent effective values. The script also uses standard floating-point arithmetic. When the exponent is extremely negative, the computed rate may underflow to zero. That numerical zero should be interpreted as “far too small to matter within this calculation,” not as a statement that the exact mathematical expression is literally equal to zero.

Another limitation is that the page reports only the leading contribution rather than the full infinite series. For many first estimates that is the right compromise: it keeps the interface understandable and the model transparent. Still, once the field becomes large enough for the rate to explode, the simple picture can cease to be the whole story. If you are doing precision theory work or connecting directly to an experiment, you should treat this page as the first stop in a longer analysis rather than the final word.

The safest workflow is straightforward. Choose physically sensible SI inputs, compute the rate, compare the field to the reported critical field, and then ask whether the constant-field idealization is reasonable for your use case. If the answer depends on pulse duration, geometry, or depletion of the external field, you have learned something important: the order-of-magnitude estimate is useful, but the next layer of modelling now matters.

A general calculator view of the same physics

Even with specialized physics underneath, this page is still a calculator in the basic mathematical sense: it maps a set of inputs to a result. That broad view can be helpful when you compare scenarios, because you can ask which input moved and how strongly the output responded.

Many scientific estimates also build from weighted contributions or summed terms, which is why series ideas appear so often in physics and numerical tools.

For the Schwinger problem, the most important lesson from that generic viewpoint is simple: not every input matters equally. The electric field enters in a way that is far more dramatic than a linear scale-up. When you explore the form below, change one parameter at a time and watch how the ratio E/E_c organizes the behaviour.