The observed vortices in each BaTiO3 layer aligned along a specific crystal direction and formed a robust striped pattern. Earlier polar textures of this kind had been observed almost exclusively in lead-containing superlattices comprising ultrathin ferroelectric and dielectric titanate layers on oxide substrates. The reported system therefore extends the material and substrate combinations in which polar topological textures can occur. The authors describe the result as a possible route toward incorporating complex ferroelectric states into silicon-based technologies, but the research does not demonstrate a device-level application.

Mechanical strain from the silicon substrate was identified as an important factor in producing the pattern. The first author, doctoral student Valentin Hevelke, said that this was the first realization of periodically ordered vortices in BaTiO3 superlattices on silicon and that silicon provides the strain needed to induce the reported pattern. The samples also showed reversible temperature-dependent behavior: heating caused a structural phase transition and the polar texture disappeared, while cooling to room temperature caused the polar texture to re-form. The supplied release does not quantify switching speed, electrical control, endurance or operation under applied fields, so it does not establish a practical switching process or device function.

Phase-field modelling by a group at Pennsylvania State University reproduced the coupling between the BaTiO3 and SrTiO3 layers through electrostatic forces and mechanical strain. The simulations also explained how the vortex tubes in the BaTiO3 layers align and provided a model for how the polar-texture pattern develops as strain changes. Overall, the work combined direct experimental characterization of a lead-free, silicon-supported superlattice with modelling of the electrostatic and mechanical effects associated with the observed vortex arrangement.

Subject of Research: Experimental study of polar vortices in lead-free BaTiO3/SrTiO3 superlattices on silicon

Article Title: Patterns of polar vortices in lead-free superlattices on silicon

Article References: Patterns of polar vortices in lead-free superlattices on silicon. Original publication (Bibliographic details unavailable.)

Image Credits: AI Generated

DOI: Not provided

Keywords: polar vortices, BaTiO3, SrTiO3, silicon, superlattices, ferroelectric materials, molecular-beam epitaxy, phase-field modelling, mechanical strain, Patterns, polar, vortices

The classical picture of solute transport descends from work by Sir Geoffrey Taylor in the 1950s, who showed that a dissolved substance carried through a smooth tube spreads in a predictable, Gaussian fashion. This behavior, known as Fickian transport, assumes that mixing is governed by ordinary diffusion superimposed on steady advection, so that a tracer pulse spreads symmetrically around its mean arrival time. In the real subsurface, however, water rarely flows through smooth channels. Rock fractures have rough, uneven walls, variable apertures, and dead-end pockets that generate swirling recirculation zones. Under these conditions the transport becomes non-Fickian: part of the contaminant cloud races ahead of the predicted front, while another part trails behind in a long tail that can persist for years.

This early-arrival and long-tailing behavior is not a mere academic curiosity. It determines how quickly a plume of industrial solvent might reach a drinking-water well, how far seawater intrudes into coastal aquifers, and how reliably repositories for nuclear waste can be sealed against radionuclide migration. The new study tackles the problem head-on by performing a systematic sequence of numerical solute transport simulations in single fractures with flow in the nonlinear regime, deliberately varying four key parameters: the joint roughness coefficient, which quantifies the waviness and asperity of the fracture walls; the inlet velocity of the flowing water; the aperture, which is the width of the opening between the two rock surfaces; and the molecular diffusion coefficient of the dissolved species.

The results reveal a clear hierarchy of influences. Increasing roughness, inlet velocity, or aperture all amplify the early arrival and the long trailing of the contaminant cloud, whereas a larger diffusion coefficient diminishes both effects. In physical terms, rougher and wider fractures at higher flow rates create stronger and more persistent eddies that hijack portions of the flow, while molecular diffusion acts as a smoothing agent that helps solute escape from trapped regions and rejoin the main current. The Peclet number, defined as the product of inlet velocity and aperture divided by the diffusion coefficient, emerges as the master variable: dispersivity increases systematically with this dimensionless ratio, alongside aperture and roughness.

The most striking advance of the work lies at the mesoscopic scale, where the authors turned their attention to the vortices themselves. Within a single rough fracture, the area occupied by an individual vortex follows an exponential distribution, meaning that small eddies are common while large ones are comparatively rare. The kinetic energy of these vortices, by contrast, follows a lognormal distribution, indicating that a few energetic whirls dominate the stored rotational energy of the flow. These statistical fingerprints provide a compact mathematical description of the hidden turbulence that classical transport models simply average away.

Crucially, the researchers show how these vortices distort the concentration field in a mechanistically intuitive way. Recirculating eddies squeeze and narrow the main flow channel, forcing the bulk of the advective transport through a restricted corridor. Solvent that drifts into the eddy periphery is detained in low-concentration zones, and this partitioning between a fast mainstream and sluggish recirculation pockets produces exactly the two-faced signature of non-Fickian transport: a leading edge that arrives early and a trailing tail that refuses to die away. The vortices are thus not a secondary detail of fracture hydraulics but the primary engine of anomalous transport.

From these observations the team constructed two complementary predictive models. At the macroscopic level, they expressed the hydrodynamic dispersivity as a function of the Peclet number, the aperture, and the joint roughness coefficient, yielding a dispersivity model that engineers can evaluate from measurable fracture properties. At the mesoscopic level, they showed that dispersivity increases with the mean vortex area and with both the mean and the standard deviation of the vortex kinetic energy, and they reformulated the hydrodynamic dispersivity in terms of the vortex distribution parameters, the aperture, and the diffusion coefficient. The resulting macroscopic-mesoscopic model links the statistical anatomy of the eddy field directly to the spreading coefficient used in large-scale transport calculations.

The authors validated their formulations by running additional simulation cases and comparing the predicted dispersivity against existing models, reporting that the new expressions reproduce the simulated transport behavior with superior accuracy. This kind of verification matters because dispersivity is the single most consequential parameter in practical contaminant-transport assessment, and it has historically been obtained by curve-fitting tracer breakthrough curves rather than derived from first principles. A model that can estimate dispersivity from fracture roughness, aperture, and flow conditions, or from vortex statistics obtainable through imaging and computation, promises substantially more reliable forecasts of plume evolution in fractured rock.

The implications extend naturally beyond a single fracture. Fractured aquifers and petroleum reservoirs are networks of intersecting fractures, and the behavior of a network inherits the physics of its individual segments. By establishing a defensible basis for predicting transport in one rough fracture under nonlinear flow, the study lays the groundwork for scaling up to full fracture networks, a step the authors explicitly identify as the foundation of their approach. Prior work by other groups has documented eddy-driven tailing and inertial mixing at fracture intersections, and the new distribution-based framework offers a vocabulary for integrating those observations into quantitative prediction.

For the broader community of hydrogeologists, environmental engineers, and regulatory scientists, the message is that the invisible architecture of swirling water inside a crack controls the visible fate of a pollutant miles downstream. Roughness that once seemed like an inconvenient detail now has a documented, parameterized role, and the vortices that roughness breeds can be counted, measured, and folded into transport models. As groundwater contamination, coastal salinization, and nuclear waste isolation grow more pressing worldwide, tools that convert mesoscopic eddy statistics into macroscopic dispersivity estimates may become indispensable for protecting the water that hundreds of millions of people drink.

Subject of Research: Non-Fickian solute transport of contaminants in rough rock fractures and its link to vortex dynamics and hydrodynamic dispersivity

Article Title: Macroscopic and mesoscopic characteristics of non-Fickian solute transport in rough fractures

Article References: Li, C., Wang, Z., Liu, J., Qiao, L., & Yang, J. (2026). Macroscopic and mesoscopic characteristics of non-Fickian solute transport in rough fractures. Hydrogeology Journal. https://doi.org/10.1007/s10040-026-03175-3

Image Credits: AI Generated

DOI: 10.1007/s10040-026-03175-3

Keywords: groundwater contamination, solute transport, rough fractures, non-Fickian transport, vortices, hydrodynamic dispersivity, Peclet number, joint roughness coefficient, numerical simulation, hydrogeology, seawater intrusion, nuclear waste disposal

The distinction between what is measured and what is real comes down to the physics of vortices. Type-II superconductors are characterized by a magnetic penetration depth that exceeds the superconducting coherence length, which means magnetic fields penetrate the material as quantized whirlpools of current called vortices. Below a lower critical field, no vortices enter; above an upper critical field, they crowd the sample and superconductivity dies entirely. When a current flows, each vortex feels a Lorentz force. If that force exceeds the pinning that defects in the material provide, the vortices begin to move, dissipating energy and heating the sample until it flips into its normal, resistive state. The critical current density measured this way, denoted Jc, therefore reflects the defect landscape of a particular sample rather than any intrinsic property of the superconducting condensate.

Physicists have long defined a more fundamental quantity: the thermodynamic depairing current density, Jc*. This is the current at which the energy shift imparted to quasiparticles by the supercurrent equals the superconducting energy gap itself. Beyond this point, it is no longer energetically favorable for electrons to remain bound in Cooper pairs, and the condensate dissociates into normal charge carriers. Unlike Jc, Jc* depends only on microscopic parameters such as the superfluid density and the coherence length, making it a genuine material property and the theoretical ceiling for supercurrents. Yet despite decades of theoretical predictions and indirect hints, directly measuring Jc* has remained out of reach, because any direct current strong enough to approach it sets vortices in motion within nanoseconds, destroying the superconducting state before depairing can be observed.

The new study, published in Nature Physics by E. Wang of the Max Planck Institute for the Structure and Dynamics of Matter in Hamburg and colleagues, exploits a window in time where vortex physics simply cannot keep up. Vortices move at speeds limited to a few tens of kilometres per second, which means that within a two-picosecond pulse they can travel only a few tens of nanometres, nowhere near enough to cross a sample tens of micrometres wide. On these timescales, vortex motion is inertially frozen near the sample edges, dissipation is negligible, and the bulk of the superconductor experiences a pure, dissipationless supercurrent pushed all the way to its intrinsic limit. The experiment effectively outruns the very mechanism that has always masked the depairing threshold.

The experimental platform is an elegant piece of ultrafast engineering. The team fabricated thin films of two archetypal superconductors: niobium nitride (NbN), a disordered s-wave superconductor with a critical temperature of 15 kelvin, and yttrium barium copper oxide (YBCO), a d-wave cuprate superconducting at 85 kelvin. Each film was integrated into a coplanar waveguide with three pairs of photoconductive switches. Femtosecond laser pulses at 515 nanometres strike the middle pair of voltage-biased switches, launching current pulses with a full width at half maximum of roughly two picoseconds. Unbiased switches on either side sample the incoming, reflected and transmitted pulses at controlled time delays, allowing the researchers to reconstruct exactly how much current the superconductor carried and how the electromagnetic field responded.

The baseline measurements confirmed the familiar picture. Conventional direct-current tests yielded critical current densities of approximately 100 gigaamperes per square metre for NbN at 7 kelvin and slightly below 50 gigaamperes per square metre for YBCO at 50 kelvin. Below the critical temperature, the picosecond response was unmistakably superconducting: the reflected pulse showed a small inductive feature proportional to the kinetic inductance of the condensate, and about 90 percent of the incoming pulse transmitted through the sample. Above the critical temperature, the response became resistive, with the pulse partially reflected and partially transmitted, exactly as expected for a normal metal.

Then came the dramatic part. When the researchers drove NbN at 7 kelvin with pulses reaching about 6.5 times the conventional critical current density, the sample’s response suddenly resembled its own normal-state behaviour at 20 kelvin, indicating that superconductivity had been transiently suppressed. Tracking the transmitted peak electric field as a function of the peak current density revealed a sharp, well-defined drop at roughly 2.2 times Jc, corresponding to about 220 gigaamperes per square metre. Below that threshold, the transmitted current was carried by Cooper pairs; above it, the current was carried by normal carriers, and the slope of the response matched the normal state. Between the two regimes, the transmitted current actually decreased with increasing drive, a hallmark of rapid Cooper pair breakup. The local maximum in transmitted current marks the depairing current density Jc*, observed directly for the first time in this setting.

YBCO told a strikingly different story. Instead of a sharp threshold, the d-wave cuprate showed a gradual, continuous suppression of superconductivity beginning well before the conventional critical current was reached, with no local maximum in the transmitted current at all. The difference traces back to the symmetry of the superconducting gap. In an s-wave superconductor like NbN, the gap is isotropic, so every Cooper pair breaks at the same critical current, producing a clean threshold. In a d-wave superconductor like YBCO, the gap varies from zero along the nodal directions to its maximum along the antinodal directions. Pairs in different parts of the Fermi surface depair at different currents, so partial depairing begins immediately and accumulates continuously, and a well-defined Jc* simply does not exist. The picosecond technique thus doubles as a probe of gap symmetry in materials that cannot easily be accessed by conventional spectroscopy.

The researchers took care to rule out alternative explanations. All measurements were performed under zero-field cooling with ambient magnetic fields shielded below one microtesla, ensuring no pre-existing vortices in the sample. Even in the most pessimistic scenario, a vortex could move no more than a few tens of nanometres during a pulse, making full penetration implausible. And the stark contrast between the s-wave and d-wave responses points to a mechanism rooted in the microscopic gap structure rather than in vortex dynamics. On the theory side, simulations based on Bardeen-Cooper-Schrieffer theory in the dirty limit, using the Usadel equation with parameters taken directly from experimental characterization, reproduced the sharp transmission drop in NbN, attributing it to the strong nonlinearity of the superfluid density. Time-dependent Ginzburg-Landau theory, valid only near the critical temperature, could not capture this nonlinearity. The calculated temperature dependence of Jc* also agreed closely with the measurements, with experimental values slightly below theory, plausibly because thermal fluctuations initiate depairing marginally early. At 0.8 times the critical temperature, the conventional critical current was nearly an order of magnitude smaller than the intrinsic depairing current.

The implications extend in two directions. Fundamentally, the technique opens a route to probing superconducting properties, including gap symmetry and condensate stiffness, in regimes that conventional transport cannot touch, and it may eventually help interrogate cuprate superconductors for which no complete microscopic theory exists even in equilibrium. Practically, the results show that the maximum supercurrent a type-II superconductor can sustain is accessible on picosecond timescales, suggesting a platform for generating ultrashort, intense magnetic field pulses and for superconducting electronics that operate closer to their intrinsic current limits. For a field in which the gap between what materials can do and what engineers can measure has persisted for more than half a century, outrunning the vortices may prove to be one of the most consequential ideas in modern superconducting research.

Subject of Research: Picosecond measurement of the intrinsic depairing current density in type-II superconductors

Article Title: Probing picosecond depairing currents in type-II superconductors

Article References: Wang, E., Chavez-Cervantes, M., Satapathy, J., Matsuyama, T., Meier, G., Zhang, X., You, L., Marijanovic, F., Curtis, J. B., Demler, E., & Cavalleri, A. (2026). Probing picosecond depairing currents in type-II superconductors. Nature Physics. https://doi.org/10.1038/s41567-026-03469-z

Image Credits: AI Generated

DOI: 10.1038/s41567-026-03469-z

Keywords: superconductivity, depairing current, type-II superconductors, picosecond pulses, Cooper pairs, vortices, NbN, YBCO, s-wave gap, d-wave gap, ultrafast transport, BCS theory

The story began more than a decade ago, when Hyun Youk, a professor of physics at Illinois, started his laboratory with a broad interest in how complex dynamics can arise from simple deterministic rules, particularly in systems of living cells that interact to form spatial patterns. In 2020, his group computationally explored how cells might communicate by secreting molecules and sensing those released by their neighbors, and they discovered communication modes that matched the ways real cells in nature construct the very same kinds of spatial patterns. That work led the team to a cellular automaton, a discrete grid of individual cells, each occupying one of a finite number of possible states and updating according to fixed rules. Cellular automata are famous for generating astonishingly sophisticated behavior—patterns that appear to move, interact, and even self-replicate—from deceptively simple local rules, which is why they have long served as models of pattern formation in physics, computer science, and biology.

The automaton Youk’s team studied was designed to mimic living tissue. Each cell takes on one of four states, and after every timestep it changes state according to a gene-circuit-like rule that imitates how real cells secrete and sense signaling molecules. The grid also carries periodic boundary conditions, meaning each edge wraps around to meet its opposite edge, much like the classic video game character Pac-Man reappearing on the left side of the screen after exiting the right. When simulated, the automaton starts from a randomly chosen initial configuration and always terminates in exactly one of three final fates: a static configuration of same-state cells, a moving rectilinear wave, or a moving spiral wave. Because the number of possible configurations is finite—though astronomically large—the system is, by the strict mathematical definition, not chaotic. Yet a difference of just a single cell’s state between two initial configurations often produced entirely different fates, with nothing visually obvious to explain why.

To test whether this apparent complexity amounted to genuine unpredictability, the researchers turned to machine learning. They gave a suite of algorithms a deceptively simple binary-classification task: given an initial configuration, forecast whether the automaton would end in a static fate or a moving-wave fate. The results were startling. Half the time the algorithms guessed correctly, and half the time they were completely wrong—performance indistinguishable from a coin flip. Here was a deterministic system whose future could not be extracted from its past by any human or machine. “Despite the simplicity of our system, the cells self-organized in a way that no human or machine could initially predict,” said Elinor Kay, an Illinois Physics graduate student and co-author of the study. The finding pointed to a source of unpredictability fundamentally distinct from chaos, one that had gone largely unnoticed in discussions of determinism and forecasting.

The breakthrough came when the team changed how they looked at the system. Instead of visualizing the four cell states as colors, they reimagined them as arrows pointing up, down, left, or right. What emerged was striking: cores of oppositely pointing arrows, completely surrounded by closed loops of arrows tracing circular paths, forming structures the researchers called vortices. Each vortex could be labeled positive, negative, or neutral depending on its orientation and degree of circular tracing. Youk and his former student, lead author Lars Koopmans, found the vortices to be highly dynamic entities. “Lars and I were looking at a lot of simulation movies, and watching them eventually trained our eyes to see these vortices,” Youk recalled. “At first, they were oddities, but then we noticed that their abundances always decreased over time.”

The vortices behaved like Brownian particles, drifting across the lattice, merging, and annihilating one another after an initial period of rapid proliferation. Their fate carried decisive information: if every vortex vanished, the system settled into a static configuration or a rectilinear wave, whereas the survival of even a single vortex guaranteed a spiral wave. The team also noticed that every positive vortex was tethered to a negative one, and vice versa, by an unbroken string of same-state cells, and that positive and negative vortices always appeared and disappeared together. The system’s total charge, therefore, remained exactly zero throughout every simulation—an unexpected conservation law that was nowhere encoded in the automaton’s update rules. Some of this behavior traced back to the periodic boundary conditions: whenever a vortex loop crossed a lattice boundary, it reentered the grid and generated a new vortex core of opposite orientation and charge, explaining both the pairing of charged vortices and the conservation of total charge. “We originally imposed periodic boundary conditions because they made the simulations easier to implement,” Youk said. “But it turns out this seemingly trivial decision is actually very important. Without these boundary conditions, we wouldn’t have unending waves or pairing of oppositely charged vortices.”

Something deeper was at work, rooted in topology. A two-dimensional grid with periodic boundary conditions is topologically equivalent to a torus—the surface of a donut—which can be visualized by gluing the lattice’s opposite edges together. Because of this connectivity, certain strings of same-state cells wrap around the donut and cannot be untangled from its hole; these are known as noncontractible loops. The researchers tested whether these loops could serve as predictors of fate. In simulations destined for static or rectilinear-wave outcomes, the average count of noncontractible loops fluctuated near zero before the final vortex annihilation, dipping to zero with increasing frequency as the deciding moment approached. In principle, observing such rapid drops could signal an impending static or rectilinear fate before it actually occurred. But in simulations bound for spiral-wave fates, the loop count simply fluctuated near one, a condition that still permits either non-spiral outcome. The loops, in other words, could predict some fates but not all.

The decisive predictor emerged from a generalization of these loops that the team calls the winding field, a topological quantity capturing the lattice’s full spatial configuration—precisely the kind of information to which the best-performing machine-learning model, a convolutional neural network, had proven sensitive. Rather than examining only wrapped strings, the winding field considers entire connected regions of same-state cells that wrap around the torus, assigning to each cell a pair of numbers, a winding vector, that quantifies how many times its region winds around the donut. At the start of nearly every simulation, the field was zero. But within the first one percent of each run’s runtime, the field switched on as vortices materialized, then self-organized as the dynamics unfolded, expanding as vortices moved and annihilated. Crucially, the field then changed in one of two characteristic ways: it either abruptly engulfed the entire lattice, producing a static or rectilinear fate, or it expanded gradually, leaving persistent regions of zero winding field and producing a spiral wave. That distinction allowed reliable discrimination among all three fates. When the winding field was fed to the convolutional neural network, its accuracy climbed steadily over the course of each simulation, from no better than random guessing at the start to near perfection at the end. “This shows that information is always present but slowly becomes accessible, which is very exciting because it implies that there’s a greater order just below our grasp,” Kay said.

The discovery challenges traditional notions of predictability that sit at the intersection of determinism and uncertainty, and Youk’s laboratory is among the first to probe these ideas. “We didn’t expect to find a new source of unpredictability at all,” Youk admitted, describing the finding as an accidental discovery born of the lab’s general interest in self-organizing systems. “What’s most exciting is the way our results exemplify the rich dynamics and layers of order that can form out of only local rules,” Kay added. Much remains to be understood. The team hopes to explore the system’s topological details more fully and to pin down a rigorous definition of predictability, which they currently define operationally as the ability of a human observer or machine-learning model to forecast fate better than chance. “So far, we haven’t come up with a deep answer to why topology matters so much in our simulations,” Youk concluded. “Rigorously defining predictability and examining its properties are our next goals.” The work was funded by the National Institutes of Health through an NIH-NIGMS R35 grant under Grant No. GM147508 and by the National Science Foundation’s Science and Technology Center for Quantitative Cell Biology under Grant No. 2243257.

Subject of Research: Emergence of predictability in a deterministic, nonchaotic cellular automaton modeling living cell communication

Article Title: The emergence of predictability in an unpredictable system

Article References: The emergence of predictability in an unpredictable system. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: cellular automata, predictability, chaos, topology, emergence, self-organization, machine learning, winding field, vortices, spatial patterns, determinism, Nature Communications

The governing equations of fuzzy dark matter are deceptively compact. The model is described by the Schrödinger–Poisson system, in which a complex wave function evolves under the influence of its own self-gravity. Despite their quantum pedigree, these equations admit a fluid interpretation through the Madelung transformation, which recasts the wave function in terms of a density and a velocity field. In this fluid picture, the familiar equations of hydrodynamics acquire an additional term, the quantum pressure, which counteracts gravity on small scales and introduces a characteristic Jeans scale. Structures larger than this scale grow much as they would in the standard cold dark matter scenario, while smaller structures are suppressed, because the uncertainty principle prevents the bosons from being localized into arbitrarily tight clumps. This suppression of small-scale structure is one of the model’s most distinctive and testable predictions, and it manifests as a sharp cutoff in the linear matter power spectrum at a wavenumber that scales with the square root of the boson mass.

Simulating this quantum fluid is dramatically harder than simulating cold dark matter. Wave-based schemes, which evolve the wave function directly, must resolve not only the de Broglie wavelength but also its rapid oscillations in time. The time-step constraints scale with the square of the spatial resolution, a diffusion-like requirement that makes high-resolution simulations extraordinarily expensive. The review details how the computational cost of resolving density granulation within a halo scales as the fourth power of the boson mass parameter and as the 8/3 power of halo mass, meaning that simulations of massive halos with heavier bosons quickly become infeasible. This is the central reason why most genuine wave-based cosmological simulations to date have been restricted to small boxes of a few megaparsecs, or to relatively light boson masses, while heavier candidates and larger volumes remain out of reach.

Among wave-based methods, the pseudo-spectral approach reigns supreme in accuracy. The global Fourier method applies a fast Fourier transform to the wave function, multiplies by the kinetic operator in momentum space, and transforms back, achieving spectral convergence with errors that decay faster than any algebraic rate as resolution increases. The method is unconditionally stable and conserves mass to machine precision, making it the algorithm of choice for codes such as GAMER, AREPO, PyUltraLight and SPoS. Its Achilles heel is that the discrete Fourier transform assumes a uniform grid with periodic boundaries, which is incompatible with adaptive mesh refinement, the technique that makes cosmological simulations affordable by concentrating resolution only where it is needed. This limitation has spurred the development of alternatives, including finite-difference schemes that work on refined grids but converge far more slowly, and a newer local pseudo-spectral technique based on Fourier continuations with Gram polynomials, which achieves eleventh-order accuracy on non-periodic grids and has been implemented in the GAMER code.

Fluid-based methods offer a compelling escape from the resolution bottleneck. Because the density, velocity and phase fields are smooth in the vast low-density regions outside halos, fluid schemes do not need to resolve the de Broglie wavelength there, allowing much larger simulation volumes. Smoothed particle hydrodynamics, implemented in codes such as AX-GADGET, and meshless finite-volume methods in the GIZMO code both exploit this advantage, and the finite-volume formulation conserves mass and momentum to machine precision. But fluid approaches have a fundamental weakness: at density nodes created by destructive wave interference, the velocity and quantum potential diverge, and the phase field becomes discontinuous. These singular points, which correspond to quantized vortices threading the halos, are notoriously difficult for any fluid scheme to capture, and the review documents persistent discrepancies between SPH simulations and wave-based calculations, including disagreements over whether soliton cores and the soliton–halo relation are recovered correctly.

The most promising path forward may be hybrid schemes that combine the strengths of both formulations. In the approach developed by Schwabe and Niemeyer and implemented in the AxioNyx code, collisionless N-body particles carry wave packets on coarse levels, and the full wave function is reconstructed on fine levels using a Gaussian beam method that captures statistically correct interference patterns. A newer fully grid-based hybrid scheme in GAMER, described by Kunkel and collaborators, solves the Hamilton–Jacobi–Madelung fluid equations on coarse levels and switches to the wave formulation on refined levels, with careful phase unwrapping at the fluid–wave interfaces. These hybrid methods, coupled with adaptive mesh refinement, have enabled cosmological zoom-in simulations of individual halos down to the present day, something that pure wave-based schemes could not previously achieve, and they open the door to simulating more massive halos and heavier boson masses than ever before.

The distinctive phenomenology of fuzzy dark matter gives simulations their scientific payoff. Each virialized halo hosts a dense, stable soliton core, the ground-state solution of the halo potential, surrounded by an outer profile that closely resembles the standard Navarro–Frenk–White form. The soliton density profile is redshift-independent and obeys an elegant scaling symmetry, and its properties are linked to the host halo through a debated soliton–halo relation that recent work interprets as thermal equilibrium between the core and the surrounding granules. Beyond the soliton, the halo is permeated by stochastically fluctuating density granules, roughly a quarter of the local de Broglie wavelength in size, generated by constructive and destructive interference. These granules scatter stars and gas, dynamically heating dwarf galaxies, stellar streams and galactic disks, and their effects are now being used to place some of the tightest constraints on the boson mass, with recent analyses of ultrafaint dwarfs suggesting the mass must exceed roughly 2.2 x 10^-21 electronvolts.

Simulations also reveal subtler pitfalls that can silently corrupt results. Spurious halos, seeded by numerical noise along filaments, can form below the half-mode mass in both N-body and genuine fuzzy dark matter simulations, mimicking real objects and potentially biasing constraints derived from the halo mass function. The review shows that the global Fourier method appears largely immune to this artificial fragmentation, while hybrid and N-body schemes produce suspiciously regular clumps along filaments that lack counterparts across different numerical setups. Equally treacherous is insufficient resolution of the de Broglie wavelength in low-density regions, which underestimates flow velocities and delays structure formation, and inadequate resolution within halos, which distorts the quantum pressure and can drive unphysical halo contraction. The review demonstrates that even when a simulated halo’s mass and central soliton profile look correct, the outer profile can still be wrong, a sobering warning that agreement with theory is necessary but not sufficient evidence of numerical accuracy.

To help the community navigate these hazards, the review provides a test bench of standard problems, from Gaussian wave packets and vortex pairs to the Jeans instability, isolated solitons, isolated halos and full cosmological boxes, and it makes publicly available initial condition files for isolated-halo and cosmological simulations so that different codes can be compared on equal footing. Eigenmode methods, which construct halos as superpositions of bound states of a static potential, offer yet another efficient tool for studying interference phenomena without solving the full nonlinear system, while collisionless N-body simulations with fuzzy initial conditions remain valuable for probing large-scale structure and checking convergence. Looking ahead, the review notes that machine learning approaches, such as physics-informed neural networks, may eventually accelerate these calculations, and that the algorithms can be extended to self-interacting, vector and multi-component variants of the model. What remains clear is that simulating fuzzy dark matter has matured into a rich computational discipline in its own right, one whose progress will determine whether this quantum alternative to cold dark matter survives the confrontation with the universe it seeks to explain.

Subject of Research: Numerical simulation methods for fuzzy dark matter, an ultralight bosonic dark matter model exhibiting wave phenomena on galactic scales

Article Title: Fuzzy dark matter simulations

Article References: Schive, H.-Y. (2026). Fuzzy dark matter simulations. Living Reviews in Computational Astrophysics, 12(1), Article 1. https://doi.org/10.1007/s41115-026-00027-5

Image Credits: AI Generated

DOI: 10.1007/s41115-026-00027-5

Keywords: fuzzy dark matter, ultralight bosons, Schrödinger–Poisson equations, soliton cores, density granulation, adaptive mesh refinement, pseudo-spectral methods, quantum pressure, cosmological simulations, halo mass function, vortices, computational astrophysics