Axions, a dark matter candidate, are usually produced when their mass switches on uniformly. But if the mass appears only inside bubbles from a cosmic phase transition, the process changes. Simulations reveal two regimes: fast bubble merging boosts the axion abundance; slow expansion suppresses it via spatial gradients, like a dimmer switch. A new formula unifies these cases and matches simulations. This also reshapes the small-scale matter power spectrum, altering predictions for dense axion miniclusters, which could be observed.
Dark matter remains one of the greatest mysteries of modern physics. Back when Vera Rubin observed galaxy rotation, she proved that visible matter is insufficient to keep stars from flying apart. Among candidates for the invisible mass, the axion stands out for its elegance: solving the strong CP problem of the Standard Model, it naturally fills the universe with cold matter. In the classic scenario, the axion field, like an orchestra awaiting the conductor, freezes after the Big Bang and then smoothly begins to oscillate when the expansion rate matches its mass. But this symphony can sound different if the conductor's baton sweeps not smoothly but abruptly — via a first-order phase transition in the early universe.
Precisely such a scenario was considered by scientists inspired by Stephen Hawking's ideas on quantum fluctuations and Georges Lemaître's on the expanding universe. In their lattice simulation, bubbles of true vacuum nucleated and expanded, like shock waves from thunder, carrying with them a new value of the axion mass. The key parameter turned out to be the ratio of the phase transition duration to the natural oscillation period of the field. If the bubbles meet and collide before the field would have started to swing (fast transition), the mass turns on synchronously throughout the volume — as if the conductor waits for all musicians to take their places and only then gives the downbeat. The delay results in a loud chord: the axion “forgets” to dilute due to expansion, and its relic density skyrockets. In the opposite limit of a slow transition, bubble walls act as dampers — they excite shock waves that quench the oscillation amplitude, reducing the final density. This bubble displacement mechanism radically changes the picture.
The results break the usual predictions for axion cosmology. Instead of a single branch on the “mass–coupling constant” diagram, a bifurcation appears: fast transitions shift the allowed parameter boundary toward smaller decay constants, while slow ones shift it toward larger ones. This immediately affects the interpretation of JWST telescope data on dark matter miniclusters and constraints from black hole superradiance. Moreover, the same phase transition must have produced a stochastic gravitational wave background that future detectors will be able to “hear” — an independent channel for testing the model. The violation of entropy during bubble nucleation adds drama: nonequilibrium effects of the early universe cease to be an annoying obstacle and become a tool.
The prospects are fascinating. Work is already underway to adapt the calculations to the QCD axion, whose mass depends on temperature, and to include the Peccei–Quinn gauge symmetry, where the phase transition may be dynamic. The semi-analytic formula for the relic density derived in the paper opens the way to quick parameter estimation without cumbersome simulations. In essence, we are learning to read the score of the cosmic symphony: the moment when the conductor swept the baton determines the entire melody of the dark sector. And now we know that this instant might not have been a smooth crescendo but a timpani strike, whose echoes still reverberate.
🎯 If the phase transition had occurred a little later, axions could have supercooled, and we would live in a completely different universe — with much denser miniclusters that might have already collapsed into black holes.