Axion Misalignment Relic Density Calculator

Introduction: Axion Dark Matter from the Misalignment Mechanism

Axion misalignment turns an initially displaced field value into a cold relic population once the Hubble rate drops below the axion mass. In the standard story, the angle θ_i is nearly frozen at early times, then the field begins to oscillate as the universe cools, and those oscillations behave like pressureless matter. This calculator uses that picture to estimate the present-day relic density from f_a and θ_i, along with the corresponding axion mass and a simple underproduced, matches, or overproduced label relative to the familiar Ω_DM h² ≈ 0.12 benchmark.

To keep the estimate aligned with the standard axion-mass relation, the calculator uses the widely cited chiral-perturbation-theory scaling between the axion mass and decay constant:

This inverse relation is why larger decay constants correspond to lighter axions. The misalignment abundance is then approximated by

The calculator is therefore most useful as a quick small-angle estimate in a standard radiation-dominated cosmology. In this implementation, θ_i enters quadratically and f_a enters with the 7/6 power, so a modest change in either input can shift the answer by a large factor. The observed dark matter density corresponds to Ω_DM h² ≈ 0.12, so the result also tells you whether the chosen axion scenario under-produces, roughly matches, or over-produces dark matter relative to that target. By entering f_a and θ_i, the calculator returns the axion mass, the relic density, and a classification indicating whether misalignment alone could account for all of the dark matter.

The misalignment estimate depends on assumptions about the early universe. The simplest reading of the formula is the one used here: Peccei–Quinn symmetry is broken before inflation, the initial angle is homogeneous across the observable patch, and the universe follows the usual radiation-dominated expansion until oscillations begin. If symmetry breaking happens after inflation, then strings and domain walls add more axion production and the relic abundance can move by order-unity factors or more. Likewise, near the top of the potential, anharmonic corrections become important and the small-angle approximation starts to understate the final density. For that reason, this calculator should be treated as a baseline guide rather than the last word on any specific axion model.

Axion cosmology connects the invisible early universe to observable particle-physics searches. The same decay constant that feeds the misalignment estimate often appears in the couplings probed by haloscopes, helioscopes, and astrophysical studies, so a cosmological number is often the first step in narrowing a viable mass range. If a model is tuned to produce all of the dark matter, the result gives you a sense of how much room exists before the abundance becomes too large. If a model is intended to make only part of the dark matter, the calculator shows how far below the benchmark the scenario lands. That kind of back-of-the-envelope filtering is useful long before a full numerical scan.

The energy density stored in the oscillating axion field is often described as a coherent zero mode. Once oscillations begin, the energy redshifts like cold matter, which is why the misalignment mechanism is so compelling as a dark matter source. The compact scaling used here hides a great deal of physics: the onset of oscillations occurs when 3H(t) ≈ m_a, the field amplitude is set by the initial angle, and the late-time abundance follows from how the comoving energy evolves afterward. More elaborate treatments can include temperature-dependent masses, non-standard expansion histories, entropy injection, or detailed lattice-inspired corrections. Those additions matter in precision work, but they are not needed for a quick planning estimate.

The table below illustrates what the calculator returns for two representative axion misalignment choices:

These examples show the two main levers in the standard misalignment picture. Increasing f_a makes the axion lighter, but it also raises the relic density through the displayed scaling because the oscillations begin later and the field stores more comoving energy. Decreasing θ_i is usually the more efficient way to suppress the abundance because the angle enters as θ_i². That is why a small adjustment to the initial angle can have an effect comparable to a substantial move in the decay constant. If the benchmark target matters to your model, the table makes it easy to see which side of the line you are on.

For practical model building, the calculator is best used to compare a few hand-picked scenarios. Try a value of f_a near the scale your model prefers, then vary θ_i to see how quickly the abundance changes. If the answer lands far above Ω_DM h² ≈ 0.12, the model may need a smaller angle, additional dilution, or a different production history. If the result is far below the benchmark, the axion would be only a subcomponent of dark matter unless some other production mechanism fills the gap. This is the sort of quick diagnostic that can save time before moving to a more detailed cosmology code.

For students and researchers, the appeal of the misalignment mechanism is that it links simple inputs to rich physics. A single angle and a single decay constant control the abundance, but those inputs encode the state of the field, the symmetry-breaking scale, and the timing of oscillation onset. The calculator compresses that story into a handful of numbers, which is useful when you want intuition rather than a full derivation. It is also a good reminder that cosmology can be highly sensitive to initial conditions: two models with similar masses can give very different relic densities once the starting angle changes. That sensitivity is part of what makes axion dark matter both challenging and interesting.

In summary, the axion misalignment mechanism remains one of the cleanest ways to relate particle theory to a cosmological abundance. By estimating the relic density from the decay constant and initial angle, this calculator gives you a fast way to test whether a model sits near the dark-matter target or drifts away from it. Whether you are screening parameter choices, comparing papers, or just building intuition about how a displaced field becomes a cold relic, the tool provides a focused first pass. It is intentionally simple, but the physics it summarizes reaches all the way from the Peccei–Quinn scale to the structure of the present-day universe.

Documenting Your Axion Estimate

Once the axion misalignment calculator shows a relic density, use the copy button to save the mass, abundance, and classification for later comparison. Keeping f_a, θ_i, m_a, and Ω_ah² together in your notes makes it much easier to compare one cosmology against another.

Sharing the copied output with collaborators is especially helpful when several groups are exploring the same axion parameter space. Others can reproduce the exact angle and decay constant, check whether they are using the same misalignment assumptions, and decide whether a different early-universe history should be tested next.

How to use this axion calculator

1. Enter Axion Decay Constant f a (GeV) using a GeV value that matches the axion model you want to test.
2. Enter Initial Misalignment Angle θ i (rad) in radians, since the abundance scaling depends on the angle exactly as written.
3. Run the estimate and compare it with a second axion scenario before treating the result as your final conclusion.

Formula: how the axion misalignment estimate is built

The calculator follows the relations shown above rather than a generic algebraic placeholder. In this axion-specific setting, f_a controls the mass scale inversely and also appears in the relic-density scaling, while θ_i enters quadratically. That means the initial angle is often the fastest way to move the answer, but changes in the decay constant still matter because they affect both the axion mass and the onset of oscillations. Keep f_a in GeV and θ_i in radians so the estimate uses the same conventions as the equations above.

Worked example: compare two axion misalignment scenarios

As a concrete axion misalignment check, enter f_a = 1×10¹² GeV and θ_i = 1 rad. The calculator returns m_a ≈ 5.7×10⁻⁶ eV and Ω_ah² ≈ 0.18, which is above the usual 0.12 dark-matter benchmark. If you keep the same formulas but lower the decay constant to 5×10¹¹ GeV and the angle to 0.5 rad, the result becomes m_a ≈ 1.14×10⁻⁵ eV and Ω_ah² ≈ 0.020. That second case is clearly underproduced, which makes it easy to see how strongly the abundance responds to the initial angle and the symmetry-breaking scale.

Limitations and assumptions for axion misalignment estimates

This calculator is a planning estimate for the standard axion misalignment picture, not a full numerical treatment of every cosmological wrinkle. It assumes a homogeneous initial angle, the small-angle approximation shown above, and a radiation-dominated universe until the axion field starts oscillating. If your model includes post-inflation symmetry breaking, strings or domain walls, entropy injection, an early matter-dominated era, or large anharmonic corrections near θ_i ≈ π, the simple abundance estimate can shift noticeably. Results therefore depend on accurate inputs and on whether the standard-history assumptions match the scenario you are studying. It does not replace model-specific review, dedicated simulation, or source material that handles the early universe in more detail.