Search NASA⌕ Search

SEARCH · Search NASA

Results for “Self-Organized Criticality”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An atomistic theory of nucleation: Self-organization via non-equilibrium work and fluctuations

For this work, insights from non-equilibrium statistical mechanics, highlighting the role of work and fluctuations at the microscale, are applied toward the development of a fundamental, rigorous and purely atomistic theory of nucleation. Nanoscale fluctuations in order, density and heat influence the local nucleation rate by orders of magnitude, necessitating their inclusion through a modern approach. Coarse-graining over the underlying Hamiltonian dynamics allows derivation of a microscale expression for the nucleation rate in terms of a classical path integral over far from equilibrium trajectories and their associated work. Second law violating states at the microscale, as found from the dynamics of small critical nucleation clusters, contribute exponentially to the observable macroscale nucleation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Early self-organization of fission gas bubble superlattice formation in neutron-irradiated monolithic U-10Mo fuels

Self-organization of defect superlattices in far-from-equilibrium systems presents a promising way to mitigate swelling concerns in nuclear materials. The gas bubble superlattice (GBS) is a highly ordered, three-dimensional complex defect structure that can retain fission gasses in Uranium-Molybdenum (U-Mo) fuels. Transmission electron microscopy (TEM) investigation of monolithic U-10Mo fuel irradiated to 1.15 × 10 21 fissions/cm 3 and 1.30 × 10 21 fissions/cm 3 revealed that early-stage ordering preferentially occurs at the grain boundaries (GB) and that the critical bubble size for complete ordering is ~3 nm. Once formed at the GB, the GBS extends towards the grain interior; however, the spread in distance from the GB varies likely depending on the type and strength of the GB sink. TEM results also showed a possible correlation between the growth and evolution of the intragranular disordered bubbles and large dislocation networks. The fission product distribution in and outside of the GBS was also investigated confirming the presence of xenon in the GBS, as well as other fission products including cesium, barium, lanthanum, and cerium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Polarization and dynamic phases of aligning active matter in periodic obstacle arrays

Here, we numerically examine a system of monodisperse self-propelled particles interacting with each other via simple steric forces and aligning torques moving through a periodic array of obstacles. Without obstacles, this system shows a transition to a polarized or aligned state for critical alignment parameters. In the presence of obstacles, there is still a polarization transition, but for dense enough arrays, the polarization is locked to the symmetry directions of the substrate. When the obstacle array is made anisotropic, at low densities the particles can form a quasi-isotropic state where the system can be polarized in any of the dominant symmetry directions. For intermediate anisotropy, the particles self-organize into a coherent lane state with one-dimensional polarization. In this phase, a small number of highly packed lanes are adjacent to less dense lanes that have the same polarization, but lanes further away can have the opposite polarization, so that global polarization is lost. For the highest anisotropy, hopping between lanes is suppressed, and the system forms uniformly dense uncoupled but polarized lanes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Efficient Clustering of Software Vulnerabilities using Self Organizing Map (SOM)

The common vulnerabilities and exposures (CVE) database was created with a mission to ``identify, define, and catalog publicly disclosed cybersecurity vulnerabilities''. This rich body of information can be used to enable rapid and efficient response to secure and defend cyber operations and protect critical cyber infrastructure. The main goal of this paper is to develop a visual analytics tool to enable deep analysis of CVEs using unsupervised clustering techniques. We enhance our analysis by first mapping CVEs to hierarchical-classes in Common Weakness Enumeration (CWE) using information in the National Vulnerability Database (NVD). Both the mapping and the numerical representation of CVEs are enabled by V2W-BERT, which uses natural language processing of the extensive information in NVD to generate a large tabular database of 137,226 CVE entries from 1999 to 2020, where each CVE is represented by a vector of 768 numerical features. The vectorized data is processed by Self-Organizing Maps (SOM), which is an unsupervised machine learning technique for dimensionality reduction, visual representation and clustering. Using a Torus map of 6417 units, we achieve ~10-fold data compression of ~140k CVEs using SOM. The trained map is further clustered using standard K-means clustering into 138 clusters of CVEs. We conducted a brief investigation of the rich mapping of CVEs to best-matching-units to K-means clusters, as well as CVEs to CWEs. For example, this novel mapping provided insight into the role of CWE-59 and CWE-264 in several CVEs that is otherwise hard to explore in the original data. We conclude that our this novel approach will not only enable deep analysis of the complex relationships between CVEs and CWEs, but also a mechanism to quickly respond to and design mitigation actions for rapidly evolving vulnerabilities that have not been mapped to existing CWEs.

Panchal, Khyati↗
DOE OSTI2026

Emerging hierarchical dislocation structures: Insights from scanning electron microscopy-electron backscatter diffraction in situ tensile testing and multifractal analysis

Understanding the evolution of dislocation structures during plastic deformation is critical for predicting the mechanical performance of metallic materials. In this work, we applied in situ scanning electron microscopy/electron backscatter diffraction tensile testing combined with multifractal (MF) analysis to assess deformation-induced dislocation structure evolution in solution-annealed 304 L stainless steel, both in its as-received and neutron-irradiated states (5.4 displacements per atom). The analysis of kernel average misorientation patterns revealed the formation of hierarchical dislocation arrangements that exhibit clear MF scaling behavior. Despite pronounced visual differences between nonirradiated and irradiated specimens—most notably, the appearance of dislocation channels after irradiation—the singularity spectra suggest that both conditions give rise to similar underlying hierarchical structures. MF analysis provides a quantitative measure of the spatial complexity and self-organization of dislocation patterns, highlighting the accelerated emergence and evolution of the dislocation structures in irradiated polycrystalline materials, as well as the limitation of their spatial extent. The findings indicate that irradiation not only modifies microstructure but also alters correlation-driven dislocation organization. More generally, they demonstrate that MF analysis is a powerful tool for probing mesoscale deformation mechanisms.

Dislocation structures↗
DOE OSTI2025

The Power of DESI for Photometric Redshift Calibration: A Case Study with KiDS-1000

Accurate redshift estimates are a critical requirement for weak lensing surveys and one of the main uncertainties in constraints on dark energy and large-scale cosmic structure. In this paper, we study the potential to calibrate photometric redshift (photo-z) distributions for gravitational lensing using the Dark Energy Spectroscopic Instrument (DESI). Since beginning its science operations in 2021, DESI has collected more than 50 million redshifts, adding about one million monthly. In addition to its large-scale structure samples, DESI has also acquired over 256k high-quality spectroscopic redshifts (spec-zs) in the COSMOS and XMM and VVDS fields. This is already a factor of 3 larger than previous spec-z calibration compilations in these two regions. Here, we explore calibrating photo-zs for the subset of KiDS-1000 galaxies that fall into joint self-organizing map (SOM) cells overlapping the DESI COSMOS footprint using the DESI COSMOS observations. Estimating the redshift distribution in KiDS-1000 with the new DESI data, we find broad consistency with previously published results while also detecting differences in the mean redshift in some tomographic bins with an average shifts of Delta Mean(z) = -0.028 in the mean and Delta Median(z) = +0.011 in the median across tomographic bins. However, we also find that incompleteness per SOM cell, i.e., groups of galaxies with similar colors and magnitudes, can modify n(z) distributions. Finally, we comment on the fact that larger photometric catalogs, aligned with the DESI COSMOS and DESI XMM and VVDS footprints, would be needed to fully exploit the DESI dataset and would extend the coverage to nearly eight times the area of existing 9-band photometry.

Blanco, Diana [UC, Santa Cruz; UC, Santa Cruz, Ins↗
DOE OSTI2026

Proximity Magnetism in Mn(Bi,Sb) 2 Te 4 –(Bi,Sb) 2 Te 3 /MnTe Natural Heterostructures

Magnetic topological insulators and their heterostructures provide significant opportunities to couple band topology with a nontrivial spin configuration for enhanced spintronic device performance, as well as designing magnetoelectric systems and functionalities. We find that Mn interdiffusion from MnTe when interfaced with (Bi,Sb) 2 Te 3 stabilizes as self-organized Mn(Bi,Sb) 2 Te 4 septuple lamellae among alternating (Bi,Sb) 2 Te 3 quintuple layers, as observed using scanning transmission electron microscopy and depth-sensitive polarized neutron reflectometry. We further demonstrate a valuable combination of magnetic and topological orders in these naturally formed Mn(Bi,Sb) 2 Te 4 –(Bi,Sb) 2 Te 3 heterostructures, which are exchange-coupled with MnTe. Magnetotransport experiments and quantum magnetism simulations reveal that, above its own Néel temperature T N ∼ 20 K, Mn(Bi,Sb) 2 Te 4 mediates the exchange field leading to an anomalous Hall effect at the (Bi,Sb) 2 Te 3 /MnTe interface, with an enhanced interfacial T N exceeding 200 K, approaching that of the bulk MnTe. This magnetic interface, in turn, allows a robust and deterministic spin–orbit torque switching without an external magnetic field at a low critical current density of 3 × 10 5 A cm –2 . The antiferromagnetically coupled architecture of Mn(Bi,Sb) 2 Te 4 –(Bi,Sb) 2 Te 3 /MnTe, featuring magnetic and topological proximity effects across a chalcogenide backbone, is rich in fundamental interface physics and holds the potential for practical applications in spintronics.

magnetic topological insulator↗
DOE OSTI2022

Staircase formation by resonant and non-resonant transport of potential vorticity

The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Link statistics of dislocation network during strain hardening

Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.

Crystal plasticity↗
DOE OSTI2025

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2023

Uncertain role of clouds in shaping summertime atmosphere-sea ice connections in reanalyses and CMIP6 models

Downwelling longwave radiation (DLR) driven by the atmospheric and cloud conditions in the troposphere is suggested to be a dominant factor to determine the summertime net surface energy budget over the Arctic Ocean and thus plays a key role to shape the September sea ice. We use reanalyses and the self-organizing map (SOM) method to distinguish CMIP6 model performance in replicating the observed strong atmosphere-DLR connection. We find all models can reasonably simulate the linkage between key atmosphere variables and the clear sky DLR but behave differently in replicating the atmosphere-DLR connection due to cloud forcing. In ERA5 and strongly coupled models, tropospheric high pressure is associated with decreased clouds in the mid- and high-levels and increased clouds near the surface. This out-of-phase structure indicates that DLR cloud forcing is nearly neutral, making the clear sky DLR more important to bridge JJA circulation to late-summer sea ice. In MERRA-2 and weakly coupled models, tropospheric clouds display a vertically homogeneous reduction; the cloud DLR is thus strongly reduced due to the cooling effect, which partially cancels out the clear sky DLR and makes the total DLR less efficient to translate circulation forcing to sea ice. The differences of cloud vertical distribution in CMIP6 appear to be differentiated by circulation related relative humidity. Therefore, a better understanding of the discrepancy of different reanalyses and remote sensing products is critical to comprehensively evaluate simulated interactions among circulation, clouds, sea ice and energy budget at the surface in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Unified Photometric Redshift Calibration for Weak Lensing Surveys Using the Dark Energy Spectroscopic Instrument

The effective redshift distribution n(z) of galaxies is a critical component in the study of weak gravitational lensing. Here, we introduce a new method for determining n(z) for weak lensing surveys based on high-quality redshifts and neural-network-based importance weights. Additionally, we present the first unified photometric redshift calibration of the three leading stage-III weak lensing surveys, the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS), with state-of-the-art spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We verify our method using a new, data-driven approach and obtain n(z) constraints with statistical uncertainties of the order of $σ_z$ ~ 0.01 and smaller. Our analysis is largely independent of previous photometric redshift calibrations and, thus, provides an important cross-check in light of recent cosmological tensions. Overall, we find excellent agreement with previously published results on the DES Y3 and HSC Y1 data sets, while there are some differences on the mean redshift with respect to the previously published KiDS-1000 results. We attribute the latter to mismatches in photometric noise properties in the COSMOS field compared to the wider KiDS self-organizing map-gold catalog. At the same time, the new n(z) estimates for KiDS do not significantly change estimates of cosmic structure growth from cosmic shear. Finally, we discuss how our method can be applied to future weak lensing calibrations with DESI data.

Lange, J. U. [American Univ., Washington, DC (Unit↗
DOE OSTI2022

Achievements of actively controlled divertor detachment compatible with sustained high confinement core in DIII-D and EAST

The compatibility of efficient divertor detachment with high-performance core plasma is vital to the development of magnetically controlled fusion energy. The joint research on the EAST and DIII-D tokamaks demonstrates successful integration of divertor detachment with excellent core plasma confinement quality, a milestone towards solving the critical plasma–wall-interaction (PWI) issue and core-edge integration for ITER and future reactors. In EAST, actively controlled partial detachment with T et,div ~ 5 eV around the strike point and H 98 > 1 in different H-mode scenarios including the high β P H-mode scenario have been achieved with ITER-like tungsten divertor, by optimizing the detachment access condition and performing detailed experiments for core-edge integration. For active long-pulse detachment feedback control, a 30 s H-mode operation with detachment-control duration being 25 s has been successfully achieved in EAST. DIII-D has achieved actively controlled fully detached divertor with low plasma electron temperature (T et,div ≤ 5 eV across the entire divertor target) and low particle flux (degree of detachment, DoD > 3), simultaneously with very high core performance (β N ~3, β P > 2 and H 98 ~ 1.5) in the high βP scenario being developed for ITER and future reactors. The high-β P high confinement scenario is characterized by an internal transport barrier (ITB) at large radius and a weak edge transport barrier (ETB, or pedestal), which are synergistically self-organized. Both the high-β P scenario and impurity seeding facilitate divertor detachment. The detachment access leads to the reduction of ETB, which facilitates the development of an even stronger ITB at large radius in the high β P scenario. Thus, this strong large radius ITB enables the core confinement improvement during detachment. In conclusion, these significant joint DIII-D and EAST advances on the compatibility of high confinement core and detached divertor show a great potential for achieving a high-performance core plasma suitable for long-pulse operation of fusion reactors with controllable steady-state PWIs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An electronic microemulsion phase emerging from a quantum crystal-to-liquid transition

Strongly interacting electronic systems often exhibit a complicated phase diagram that results from the competition between different quantum ground states. One feature of these phase diagrams is the emergence of microemulsion phases, where regions of different phases self-organize across multiple length scales. The experimental characterization of these microemulsions can pose considerable challenges, as the long-range Coulomb interaction microscopically mingles with the competing states. Here we observe the signatures of the microemulsion between an electronic Wigner crystal and an electron liquid in a MoSe 2 monolayer using cryogenic reflectance and magneto-optical spectroscopy. We find that the transition into this microemulsion state is marked by anomalies in exciton reflectance, spin susceptibility and umklapp scattering, establishing it as a distinct phase of electronic matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Lifelike behavior of chemically oscillating mobile capsules

Inspired by the self-organization of unicellular species into multicellular organisms, we use theory and simulation to design a system of mobile, active microcapsules that produce chemicals according to a catalytic reaction network (CRN). In solution, the catalytic reactions generate a force that drives the flow of the surrounding fluid. The convective flow, in turn, transports the capsules to new chemical surroundings with each successive translation. Consequently, the chemical signal produced by a cluster of capsules is critically dependent on the proximity and spatial configuration of the neighboring capsules. If all of the capsules and chemical products involved in the CRN lie within sufficient proximity, then the system develops oscillatory behavior that drives the dynamic self-assembly of capsules into larger clusters that are capable of collective action. Finally, these simulations indicate potential chemo-mechanic mechanisms that enable single cellular units, such as ameba, to organize into multicellular life forms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Directed Self-Assembly of the Organic Semiconductor C8-BTBT-C8 in Anodic Aluminum Oxide Nanopores

Controlling the self-assembly of organic semiconductors at the nanoscale is critical for advancing high-performance electronic and photonic devices, yet it remains challenging due to their intrinsic anisotropic crystallization and sensitivity to processing conditions. Here, in this study, we demonstrate that cylindrical nanoconfinement within anodic aluminum oxide membranes provides a versatile platform to precisely tune the molecular orientation and phase behavior of the prototypical organic semiconductor 2,7-dioctyl[1]benzothieno[3,2-b][1]benzothiophene (C8-BTBT-C8). Combining temperature-dependent high-resolution synchrotron X-ray scattering with optical birefringence measurements, we uncover that confinement geometries (pore diameters 25–180 nm) and surface chemistry govern the emergence of distinct smectic A textures, featuring molecular layers either parallel or perpendicular to the pore axis. The competition between axial and radial smectic layering is modulated by pore size, surface hydrophilicity, and thermal history, enabling reversible control over domain orientations and transitions between liquid crystalline and crystalline states. Notably, nanoconfinement stabilizes the smectic phase over an expanded temperature range compared to bulk, while inducing complex multidomain configurations owing to geometric constraints and anchoring conditions. Our results elucidate fundamental mechanisms by which anisotropic nanoscale confinement directs the self-organization of highly conjugated organic molecules, with implications for optimizing directional charge transport and anisotropic optical responses in organic–inorganic hybrid nanoarchitectures. This study establishes nanoconfinement as a powerful strategy to engineer morphology and functional properties in organic semiconducting materials with nanoscale precision.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Field-Driven Out-of-Equilibrium Collective Patterns for Swarm Micro-Robotics

Soft robotics has been rapidly advancing, offering significant improvements over traditional rigid robotic systems through the use of compliant materials that enhance adaptability and interaction with the environment. However, current approaches face critical challenges, including the reliance on complex “top-down” fabrication techniques and the difficulty of wireless powering and control at the microscale. Swarm robotics introduces a paradigm shift, leveraging collective dynamics to achieve cooperative and adaptable behaviors among multiple robotic units. Inspired by nature, this “bottom-up” approach enables swarm robots to execute task-specific reconfigurations, enhancing flexibility and robustness. Field-driven active colloids emerge as a promising platform for swarm microrobotics, capable of self-propulsion and self-organization into dynamic collective patterns under external field excitation and manipulation. These systems mimic biologically inspired swarm behaviors, such as flocking and vortex formation, providing a versatile foundation for designing innovative swarm microrobots. Here, this review discusses the principles of electric and magnetic field-driven collective self-organization, focusing on the particle dynamics, the emergence of collective swarm patterns, and illustrative examples of functional swarm microrobots. It concludes with future perspectives on harnessing these systems for adaptive, scalable, and multifunctional microrobotic applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Self-consistent dispersal puts tight constraints on the spatiotemporal organization of species-rich metacommunities

Dispersal can be critical to the maintenance of ecosystems as it allows local communities to be recolonized after extinction. However, it remains unclear whether the extinction-mitigating effect of dispersal persists when the number of competing species is large. Based on a spatially explicit mathematical description of metacommunities, we show that when many species coexist, each species operates near its extinction threshold, barely surviving due to dispersal. This has general consequences for spatiotemporal abundance patterns. For short-range dispersal, species organize into fractal spatiotemporal extinction patterns characteristic of a directed percolation phase transition. As species approach their extinction threshold, biodiversity is very sensitive to perturbation, suggesting that dispersal within a metacommunity puts tight constraints on the robustness and evolution of species-rich metacommunities. Biodiversity is often attributed to a dynamic equilibrium between the immigration and extinction of species. This equilibrium forms a common basis for studying ecosystem assembly from a static reservoir of migrants—the mainland. Yet, natural ecosystems often consist of many coupled communities (i.e., metacommunities), and migration occurs between these communities. The pool of migrants then depends on what is sustained in the ecosystem, which, in turn, depends on the dynamic migrant pool. This chicken-and-egg problem of survival and dispersal is poorly understood in communities of many competing species, except for the neutral case—the “unified neutral theory of biodiversity.” Employing spatiotemporal simulations and mean-field analyses, we show that self-consistent dispersal puts rather tight constraints on the dynamic migration–extinction equilibrium. When the number of species is large, species are pushed to the edge of their global extinction, even when competition is weak. As a consequence, the overall diversity is highly sensitive to perturbations in demographic parameters, including growth and dispersal rates. When dispersal is short range, the resulting spatiotemporal abundance patterns follow broad scale-free distributions that correspond to a directed percolation phase transition. The qualitative agreement of our results for short-range and long-range dispersal suggests that this self-organization process is a general property of species-rich metacommunities. Our study shows that self-sustaining metacommunities are highly sensitive to environmental change and provides insights into how biodiversity can be rescued and maintained.

54 ENVIRONMENTAL SCIENCES↗