@inproceedings{22002,
  abstract     = {Topological simplification is the process of reducing complexity of a function while maintaining its essential features. Its goal is to find a new filter function, which reorders cells of the input complex in a way which eliminates some persistent homological features, without affecting the rest. We present a new approach to simplification based on the concept of forbidden regions and combinatorial dynamics. It allows us to reorder and cancel critical values, whose cancellation is not possible using existing methods because they are not consecutive in the total order. Each such cancellation takes O(c⋅n) time in the worst case, where c is the number of birth-death pairs and n is the size of the input complex.},
  author       = {Leśkiewicz, Jakub and Furmanek, Bartosz and Lipiński, Michał and Morozov, Dmitriy},
  booktitle    = {42nd International Symposium on Computational Geometry},
  isbn         = {9783959774185},
  issn         = {1868-8969},
  keywords     = {persistent homology, topological simplification, depth posets},
  location     = {New Brunswick, NJ, United States},
  publisher    = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik},
  title        = {{Topological simplification guided by forbidden regions}},
  doi          = {10.4230/LIPIcs.SoCG.2026.72},
  volume       = {367},
  year         = {2026},
}

@article{20971,
  abstract     = {Mountain glaciers are among the natural systems most vulnerable to climate change. However, their interactions with the atmosphere are complex and not fully understood. These interactions can trigger rapid adjustments and climate feedbacks that either amplify or attenuate atmospheric signals, influencing both glacier response and large-scale atmospheric circulation. Observing this functional coupling in nature is challenging because the key processes occur over a wide range of spatial and temporal scales. However, recent advances in observational techniques and modeling have provided new insights into these interactions. In this review, we summarize the current state of knowledge on glacier-atmosphere interactions in high-mountain regions at different scales, and highlight recent advances in observational and numerical modeling. We also highlight important knowledge gaps and outline future research directions to improve the prediction of glacier change in a warming world.},
  author       = {Sauter, T. and Brock, B. W. and Collier, E. and Goger, B. and Groos, A. R. and Haualand, K. F. and Mott, R. and Nicholson, L. and Prinz, R. and Shaw, Thomas and Stiperski, I. and Georgi, A. and Haugeneder, M. and Mandal, A. and Reynolds, D. and Saigger, M. and Sicart, J. E. and Voordendag, A.},
  issn         = {1944-9208},
  journal      = {Reviews of Geophysics},
  number       = {1},
  title        = {{Glacier-atmosphere interactions and feedbacks in high-mountain regions - A review}},
  doi          = {10.1029/2024RG000869},
  volume       = {64},
  year         = {2026},
}

@article{20980,
  abstract     = {Morse decompositions partition the flows in a vector field into equivalent structures. Given such a decomposition, one can define a further summary of its flow structure by what is called a connection matrix. These matrices, a generalization of Morse boundary operators from classical Morse theory, capture the connections made by the flows among the critical structures—such as attractors, repellers, and orbits—in a vector field. Recently, in the context of combinatorial dynamics, an efficient persistence-like algorithm to compute connection matrices has been proposed in Dey, Lipiński, Mrozek, and Slechta [SIAM J. Appl. Dyn. Syst., 23 (2024), pp. 81–97]. We show that, actually, the classical persistence algorithm with exhaustive reduction retrieves connection matrices, both simplifying the algorithm of Dey et al. and bringing the theory of persistence closer to combinatorial dynamical systems. We supplement this main result with an observation: the concept of persistence as defined for scalar fields naturally adapts to Morse decompositions whose Morse sets are filtered with a Lyapunov function. We conclude by presenting preliminary experimental results.},
  author       = {Dey, Tamal K. and Haas, Andrew and Lipiński, Michał},
  issn         = {1536-0040},
  journal      = {SIAM Journal on Applied Dynamical Systems},
  number       = {1},
  pages        = {108--130},
  publisher    = {Society for Industrial & Applied Mathematics},
  title        = {{Computing a connection matrix and persistence efficiently from a morse decomposition}},
  doi          = {10.1137/25m1739406},
  volume       = {25},
  year         = {2026},
}

@article{21015,
  abstract     = {Early embryo geometry is one of the most invariant species-specific traits, yet its role in ensuring developmental reproducibility and robustness remains underexplored. Here we show that in zebrafish, the geometry of the fertilized egg—specifically its curvature and volume—serves as a critical initial condition triggering a cascade of events that influence development. The embryo geometry guides patterned asymmetric cell divisions in the blastoderm, generating radial gradients of cell volume and nucleocytoplasmic ratio. These gradients generate mitotic phase waves, with the nucleocytoplasmic ratio determining individual cell cycle periods independently of other cells. We demonstrate that reducing cell autonomy reshapes these waves, emphasizing the instructive role of geometry-derived volume patterns in setting the intrinsic period of the cell cycle oscillator. In addition to organizing cell cycles, early embryo geometry spatially patterns zygotic genome activation at the midblastula transition, a key step in establishing embryonic autonomy. Disrupting the embryo shape alters the zygotic genome activation pattern and causes ectopic germ layer specification, underscoring the developmental significance of geometry. Together, our findings reveal a symmetry-breaking function of early embryo geometry in coordinating cell cycle and transcriptional patterning.},
  author       = {Mishra, Nikhil and Li, Yuting I and Hannezo, Edouard B and Heisenberg, Carl-Philipp J},
  issn         = {1745-2481},
  journal      = {Nature Physics},
  pages        = {139--150},
  publisher    = {Springer Nature},
  title        = {{Geometry-driven asymmetric cell divisions pattern cell cycles and zygotic genome activation in the zebrafish embryo}},
  doi          = {10.1038/s41567-025-03122-1},
  volume       = {22},
  year         = {2026},
}

@article{21018,
  abstract     = {In this paper, we review recent results on stability and instability in logarithmic Sobolev inequalities, with a particular emphasis on strong norms. We consider several versions of these inequalities on the Euclidean space, for the Lebesgue and the Gaussian measures, and discuss their differences in terms of moments and stability. We give new and direct proofs, as well as examples and discuss the stability of a logarithmic uncertainty principle. Although we do not cover all aspects of the topic, we hope to contribute to establishing the state of the art.},
  author       = {Brigati, Giovanni and Dolbeault, Jean and Simonov, Nikita},
  issn         = {2730-9657},
  journal      = {La Matematica},
  publisher    = {Springer Nature},
  title        = {{Logarithmic Sobolev Inequalities: A review on stability and instability results}},
  doi          = {10.1007/s44007-025-00180-y},
  volume       = {5},
  year         = {2026},
}

@article{21132,
  abstract     = {We unify the variational hypocoercivity framework established by D. Albritton, S. Armstrong, J.-C. Mourrat, and M. Novack [2], with the notion of second-order lifts of reversible diffusion processes, recently introduced by A. Eberle and the second author [30]. We give an abstract, yet fully constructive, presentation of the theory, so that it can be applied to a large class of linear kinetic equations. As this hypocoercivity technique does not twist the reference norm, we can recover accurate and sharp convergence rates in various models. Among those, adaptive Langevin dynamics (ALD) is discussed in full detail and we show that for near-quadratic potentials, with suitable choices of parameters, it is a near-optimal second-order lift of the overdamped Langevin dynamics. As a further consequence, we observe that the Generalised Langevin Equation (GLE) is also a second-order lift, as the standard (kinetic) Langevin dynamics are, of the overdamped Langevin dynamics. Then, convergence of (GLE) cannot exceed ballistic speed, i.e. the square root of the rate of the overdamped regime. We illustrate this phenomenon with explicit computations in a benchmark Gaussian case.},
  author       = {Brigati, Giovanni and Lörler, Francis and Wang, Lihan},
  issn         = {1937-5077},
  journal      = {Kinetic and Related Models},
  pages        = {34--55},
  publisher    = {American Institute of Mathematical Sciences},
  title        = {{Hypocoercivity meets lifts}},
  doi          = {10.3934/krm.2025020},
  volume       = {20},
  year         = {2026},
}

@article{21164,
  abstract     = {Global emission inventories often fail to capture the complexities of vehicular pollution in regions with unique fuel mixes, such as Brazil’s extensive biofuel use, leading to significant uncertainties in atmospheric modeling. This study presents a century-long (1960–2100) bottom-up vehicular emission inventory for Brazil, leveraging locally derived emission factors. Our estimates reveal substantial discrepancies in magnitude, timing, and speciation of non-CO2 pollutants (CO, NMHC, PM2.5) compared to leading global inventories (EDGAR, CEDS, CAMS), highlighting critical inaccuracies in widely used data sets. More critically, future projections under Shared Socioeconomic Pathways (SSPs) uncover a novel positive feedback mechanism: rising temperatures significantly enhance vehicular evaporative nonmethane hydrocarbon (NMHC) emissions. This temperature-dependent increase and subsequent NMHC oxidation to CO2 suggest an overlooked pathway that could amplify climate warming and air pollution globally, particularly after a breakpoint around 2050 (p < 0.05). While historical emissions peaked in the 1990s–2000s, nonexhaust PM becomes increasingly important. Air quality simulations using our inventory in the MUSICA model show good regional PM2.5 agreement but highlight challenges in resolving local primary pollutant peaks. This comprehensive inventory provides crucial data for Brazil and uncovers globally relevant climate–chemistry interactions, urging a re-evaluation of regional specificities in global emission assessments.},
  author       = {Ibarra-Espinosa, Sergio and Dias de Freitas, Edmilson and Gaubert, Benjamin and Lichtig, Pablo and Ropkins, Karl and da Silva, Iara and Martins Pereira, Guilherme and Schuch, Daniel and Nascimento, Janaina and Hoinaski, Leonardo and Martins, Leila Droprinchinski and Gavidia-Calderón, Mario and Vara-Vela, Angel and Toledo de Almeida Albuquerque, Taciana and Ynoue, Rita Yuri and Diez, Sebastian and Mera, Zamir and Casallas Garcia, Alejandro and Vallejo, Fidel and Diaz, Valeria and Pedruzzi, Rizzieri and Abrutzky, Rosana and Franco, Marco A. and Huneeus, Nicolas and Jorquera, Hector and Belalcázar-Cerón, Luis Carlos and Rojas, Néstor Y. and de Fatima Andrade, Maria and Emmons, Louisa and Brasseur, Guy},
  issn         = {1520-5851},
  journal      = {Environmental Science &amp; Technology},
  publisher    = {American Chemical Society},
  title        = {{A century of vehicular emissions in Brazil: Unveiling the impacts of unique fuel mix on air quality}},
  doi          = {10.1021/acs.est.5c08400},
  year         = {2026},
}

@article{21217,
  abstract     = {This study investigates the mechanisms driving clustered convection and the breakdown of the Intertropical Convergence Zone (ITCZ) over the Western Pacific Warm Pool using high‐resolution cloud‐resolving simulations and machine‐learning sensitivity experiments. Results show that ITCZ breakdown episodes, marked by spatially homogeneous convection and weakened meridional moisture gradients, are triggered primarily by anomalous moisture advection linked to the equatorial Rossby‐wave activity. While large‐scale moisture advection regulates the background convective state strongly, it is the surface and low‐level meridional winds that dominate transitions between clustered and random convection. Simulations demonstrate that moisture alone can sustain convective clustering, but breakdown episodes are more persistent and widespread when coupled with southerly meridional advection. These findings confirm that wave‐driven advection acts as a regulatory mechanism, periodically disrupting convective clustering and reshaping the meridional moisture gradient. This modulation of organization by wave‐induced breakdown events is critical for understanding tropical convection variability and its implications for the climate system.},
  author       = {Casallas Garcia, Alejandro and Mark Tompkins, Adrian and Muller, Caroline J},
  issn         = {1477-870X},
  journal      = {Quarterly Journal of the Royal Meteorological Society},
  publisher    = {Wiley},
  title        = {{Moisture and wind effects of Rossby waves on Western Pacific Intertropical Convergence Zone breakdown events}},
  doi          = {10.1002/qj.70131},
  year         = {2026},
}

@article{21275,
  abstract     = {DNA methylation is a primary layer of epigenetic modification that plays a pivotal role in the regulation of development, aging, and cancer. The concurrent activity of opposing enzymes that mediate DNA methylation and demethylation gives rise to a biochemical cycle and active turnover of DNA methylation. While the ensuing biochemical oscillations have been implicated in the regulation of cell differentiation, their functional role and spatiotemporal dynamics are unknown. In this work, we demonstrate that chromatin-mediated coupling between these local biochemical cycles can lead to the emergence of phase-locked domains, regions of locally synchronized turnover activity, whose coarsening is arrested by genomic heterogeneity. We introduce a minimal model based on stochastic oscillators with constrained long-range and nonreciprocal interactions, shaped by the local chromatin organization. Through a combination of analytical theory and stochastic simulations, we predict both the degree of synchronization and the typical size of emergent phase-locked domains. We qualitatively test these predictions using single-cell sequencing data. Our results show that DNA methylation turnover exhibits surprisingly rich spatiotemporal patterns that may be used by cells to control cell differentiation.},
  author       = {Olmeda, Fabrizio and Gupta, Misha and Bektas, Onurcan and Rulands, Steffen},
  issn         = {2835-8279},
  journal      = {PRX Life},
  publisher    = {American Physical Society},
  title        = {{Spatiotemporal patterns of active epigenetic turnover}},
  doi          = {10.1103/89bj-79g5},
  volume       = {4},
  year         = {2026},
}

@article{21344,
  abstract     = {Tropospheric ozone has the potential to become an increasingly pressing public health issue in Bogotá, Colombia, due to rising concentrations across the city driven by complex interactions among emissions, meteorology, and urban structure. This study presents a comprehensive spatiotemporal analysis of ozone levels from 2013 to 2023 and assesses the associated health burden using mortality data from the same period. Results reveal a consistent upward trend in ozone concentrations, particularly in northern, western, and southern localities, with seasonal peaks linked to biomass burning and photochemical conditions. Mortality analysis, based on the Global Exposure Mortality Model, estimates that 18.3% of all deaths among individuals aged 25 and older are attributable to long-term ozone exposure. The highest burdens are found in densely populated and socioeconomically vulnerable areas such as Kennedy, Suba, and Ciudad Bolívar, with the elderly being the most affected. Building on these findings, we developed a machine learning prediction model for ozone using a convolutional merge with a long-short term memory network architecture trained on air quality and meteorological variables. The model demonstrated strong predictive performance (mean Rho=0.86, RMSE=3.5 μg/m3) across monitoring stations (17 with at least 35000 data points), supporting its potential application in real-time early warning systems across Bogotá. This integrated approach highlights the importance of localized air quality management, combining epidemiological assessment with predictive modeling. The findings underscore the urgency of implementing region-specific mitigation strategies and improving monitoring infrastructure to reduce health risks from ozone exposure in Bogotá’s rapidly growing urban environment.},
  author       = {Bustos, Daniela and Garcia, Diana and Rojas, Nestor Y. and Lopez-Barrera, Ellie A. and Peña-Rincon, Carlos and Casallas Garcia, Alejandro},
  issn         = {2509-9434},
  journal      = {Earth Systems and Environment},
  publisher    = {Springer Nature},
  title        = {{Ozone trends and mortality risk: The growing need for machine learning predictions in Bogotá, Colombia}},
  doi          = {10.1007/s41748-026-01052-3},
  year         = {2026},
}

@article{21385,
  abstract     = {We prove that the average size of a mixed character sum (math. formular) (for a suitable smooth function w) is on the order of √x for all irrational real θ satisfying a weak Diophantine condition, where χ is drawn from the family of Dirichlet characters modulo a large prime r and where x 6 r. In contrast, it was proved by Harper that the average size is o(√x) for rational θ. Certain quadratic Diophantine equations play a key role in the present paper. },
  author       = {Wang, Victor and Xu, Max},
  issn         = {1473-7124},
  journal      = {Proceedings of the Royal Society of Edinburgh: Section A Mathematics},
  pages        = {1--15},
  publisher    = {Cambridge University Press},
  title        = {{Average sizes of mixed character sums}},
  doi          = {10.1017/prm.2026.10123},
  year         = {2026},
}

@article{21504,
  abstract     = {Selecting an appropriate divergence measure is a critical aspect of machine learning, as it directly impacts model performance. Among the most widely used, we find the Kullback–Leibler (KL) divergence, originally introduced in kinetic theory as a measure of relative entropy between probability distributions. Just as in machine learning, the ability to quantify the proximity of probability distributions plays a central role in kinetic theory. In this paper, we present a comparative review of divergence measures rooted in kinetic theory, highlighting their theoretical foundations and exploring their potential applications in machine learning and artificial intelligence.},
  author       = {Auricchio, Gennaro and Brigati, Giovanni and Giudici, Paolo and Toscani, Giuseppe},
  issn         = {1793-6314},
  journal      = {Mathematical Models and Methods in Applied Sciences},
  publisher    = {World Scientific Publishing},
  title        = {{From kinetic theory to AI: A rediscovery of high-dimensional divergences and their properties}},
  doi          = {10.1142/S0218202526410010},
  year         = {2026},
}

@article{21657,
  abstract     = {We compare three global kilometer-scale models (ICON, IFS and NICAM) to clarify the advantages and challenges of high-resolution global weather and climate modeling, using different approaches to represent convection, from fully parameterized to fully explicit. Our analysis focuses on tropical precipitation characteristics spanning a wide range of spatio-temporal scales—including the diurnal cycle, extreme precipitation, convective organization, and the Madden-Julian Oscillation (MJO)—along with interactions between convection and the thermodynamic environment. All three models commonly show weaker convective organization with smaller precipitation cells than observed, though the strength of the bias varies by model. This diversity is introduced by differences in the representation of (a) convective initiation affected by the convective sensitivity to moisture and (b) tropospheric moistening associated with deep convection. Models with stronger thermodynamic-convection coupling increase environmental moisture near convection, thereby enhancing convective organization. This has important upscale effects on the MJO; while IFS and NICAM capture its eastward propagation well, ICON has difficulty reproducing it. The amplitudes and phases of precipitation diurnal cycles over land show much greater disagreement among the models than over ocean, influenced by how convection is initiated. Biases in rain evaporation and cold pool formation hinder the propagation of mesoscale convection, leading to errors such as the misrepresentation of nocturnal convection moving off the coast of Sumatra in IFS and ICON. These results highlight the importance of thermodynamic-convection coupling in realistically simulating tropical convection across scales. To improve this coupling, kilometer-scale models require better representation of the interaction between resolved convection and three-dimensional turbulent mixing.},
  author       = {Takasuka, Daisuke and Becker, Tobias and Bao, Jiawei},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {3},
  publisher    = {Wiley},
  title        = {{Precipitation characteristics and thermodynamic-convection coupling in global kilometer-scale simulations}},
  doi          = {10.1029/2025MS005343},
  volume       = {18},
  year         = {2026},
}

@article{21661,
  abstract     = {Model checking undiscounted reachability and expected-reward properties on Markov decision processes (MDPs) are key for the verification of systems that act under uncertainty. Popular algorithms are policy iteration and variants of value iteration; in tool competitions, most participants rely on the latter. These algorithms generally need worst-case exponential time. However, the problem can equally be formulated as a linear programme, solvable in polynomial time. In this paper, we give a detailed overview of today’s state-of-the-art algorithms for MDP model checking with a focus on performance and correctness. We highlight their fundamental differences, and describe various optimizations and implementation variants. We experimentally compare floating-point and exact-arithmetic implementations of all algorithms on three benchmark sets using two probabilistic model checkers. Our results show that (optimistic) value iteration is a sensible default, but other algorithms are preferable in specific settings. This paper thereby provides a guide for MDP verification practitioners—tool builders and users alike.},
  author       = {Hartmanns, Arnd and Junges, Sebastian and Quatmann, Tim and Weininger, Maximilian},
  issn         = {1433-2787},
  journal      = {International Journal on Software Tools for Technology Transfer},
  keywords     = {Quantitative model checking, Markov decision process, Linear programming, Value iteration, Policy iteration},
  publisher    = {Springer Nature},
  title        = {{The revised practitioner’s guide to MDP model checking algorithms}},
  doi          = {10.1007/s10009-026-00848-y},
  year         = {2026},
}

@article{21755,
  abstract     = {Tropical shallow clouds are a major source of uncertainty in Earth's climate sensitivity, especially through their spatial arrangement, which global climate models do not represent. Efforts to understand their organization have partly relied on classifying observed scenes, identifying four patterns as archetypal regimes. Here we analyze geostationary satellite imagery of the western tropical Atlantic using the L‐function, a tool based on point pattern theory that quantifies cloud organization across spatial scales. Classical examples of the four patterns show distinct L‐function fingerprints, revealing their characteristic clustering and regularity scales and aiding physical interpretation. Yet, when evaluating many scenes at fixed spatial scales, the L‐function distribution lacks the distinct modes expected from discrete regimes. This is corroborated by analyses of other organization indices employing diverse approaches, from inter‐cloud nearest‐neighbor distances to fractal analysis. Implications for the parameterization of mesoscale cloud organization in climate models are discussed.},
  author       = {Biagioli, Giovanni and Mandorli, Giulio and Freischem, Lilli Johanna and Casallas Garcia, Alejandro and Tompkins, Adrian Mark},
  issn         = {1944-8007},
  journal      = {Geophysical Research Letters},
  number       = {8},
  publisher    = {Wiley},
  title        = {{Spatial patterns of shallow clouds: Challenging the concept of defined regimes}},
  doi          = {10.1029/2025gl119921},
  volume       = {53},
  year         = {2026},
}

@article{21760,
  abstract     = {3I/ATLAS is the third interstellar object discovered to date, following 1I/‘Oumuamua and 2I/Borisov. Its unusually high excess velocity and active cometary nature make it a key probe of the Galactic population of icy planetesimals. Understanding its origin requires its past trajectory through the Galaxy to be traced and the possible role of stellar encounters to be assessed, both as a potential origin and a perturber to its orbit. We integrated the orbit of 3I/ATLAS backward in time for 10 Myr, together with a sample of Gaia DR3 stars with high-quality astrometry and radial velocities, to identify close passages within 2 pc. We identify 93 nominal encounters, 62 of which are significant at the 2σ level. However, none of these encounters produced any meaningful perturbation. The strongest perturber Gaia DR3 6863591389529611264 at 0.30 pc and with a relative velocity of 35 km s−1, imparted only a velocity change of ∣Δv∣  ≃  5  ×  10−4 km s−1 to the orbit of 3I/ATLAS. Our results indicate that no stellar flybys within the past 10 Myr and 500 pc contained in Gaia DR3 can account for the present trajectory of 3I/ATLAS or be associated with its origin. We further show that 3I/ATLAS is kinematically consistent with a thin-disk population, despite its large peculiar velocity.},
  author       = {Pérez-Couto, X. and Torres Rodriguez, Santiago and Villaver, E. and Mustill, A. J. and Manteiga, M.},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{3I/ATLAS: In search of the witnesses to its voyage}},
  doi          = {10.3847/1538-4357/ae56ff},
  volume       = {1001},
  year         = {2026},
}

@article{21762,
  abstract     = {Bacteria, like eukaryotes, use conserved cytoskeletal systems for intracellular organization. The plasmid-encoded ParMRC system forms actin-like filaments that segregate low–copy number plasmids. In multicellular cyanobacteria such as Anabaena sp., we found that a chromosomally encoded ParMR system has evolved into a cytoskeletal system named CorMR with a function in cell shape control rather than DNA segregation. Live-cell imaging, in vitro reconstitution, and cryo–electron microscopy revealed that CorM formed dynamically unstable, antiparallel double-stranded filaments that were recruited to the membrane by CorR through an amphipathic helix conserved in multicellular cyanobacteria. CorMR filaments were regulated by MinC, which excluded them from the poles and division plane. Comparative genomics indicated that the repurposing of ParMR and Min systems coevolved with cyanobacterial multicellularity, highlighting the evolutionary plasticity of cytoskeletal systems in bacteria.},
  author       = {Springstein, Benjamin L and Javoor, Manjunath and Megrian, Daniela and Hajdu, Roman and Hanke, Dustin M. and Zens, Bettina and Weiss, Gregor L. and Schur, Florian Km and Loose, Martin},
  issn         = {1095-9203},
  journal      = {Science},
  number       = {6795},
  publisher    = {AAAS},
  title        = {{Repurposing of a DNA segregation machinery into a cytoskeletal system controlling cell shape}},
  doi          = {10.1126/science.aea6343},
  volume       = {392},
  year         = {2026},
}

@article{21849,
  abstract     = {The development of complex tissues relies on the precise assignment of cell identity. At the molecular scale, this process depends on the deposition of epigenetic modifications—such as methylation—that are regulated by complex biochemical networks and occur at specific regions on the DNA and chromatin. Here we show that despite the complexity of epigenetic regulation, dynamical scaling and self-similarity of DNA methylation marks emerge in embryonic development. Drawing on single-cell multi-omics experiments, super-resolution microscopy and statistical physics, we demonstrate that these phenomena originate in dynamical feedback between DNA methylation and the formation of nanoscale dynamic chromatin aggregates. These nanoscale processes lead to genome-wide increase in DNA methylation marks following a power law and self-similar correlation functions. Using this framework, we identify methylation patterns that precede gene expression changes in embryonic symmetry breaking. Our work identifies linear sequencing measurements as a laboratory to study mesoscopic biophysical processes in vivo.},
  author       = {Olmeda, Fabrizio and Lohoff, Tim and Kafetzopoulos, Ioannis and Clark, Stephen J. and Benson, Laura and Santos, Fatima and Krueger, Felix and Walker, Simon and Reik, Wolf and Rulands, Steffen},
  issn         = {1745-2481},
  journal      = {Nature Physics},
  publisher    = {Springer Nature},
  title        = {{Scaling and self-similarity in the formation of the embryonic epigenome}},
  doi          = {10.1038/s41567-026-03263-x},
  year         = {2026},
}

@article{21884,
  abstract     = {We show that a randomly perturbed digraph, where we start with a dense digraph Dα and add a small number of random edges to it, will typically contain a fixed orientation of a bounded-degree spanning tree. This answers a question posed by Araujo, Balogh, Krueger, Piga and Treglown and generalizes the corresponding result for randomly perturbed graphs by Krivelevich, Kwan and Sudakov. More specifically, we prove that there exists a constant c=c(α,Δ) such that if 
T is an oriented tree with maximum degree Δ and Dα is an n-vertex digraph with minimum semidegree αn, then the graph obtained by adding cn uniformly random edges to Dα will contain T with high probability.},
  author       = {Morawski, Patryk and Petrova, Kalina H},
  issn         = {1077-8926},
  journal      = {Electronic Journal of Combinatorics},
  number       = {2},
  publisher    = {Electronic Journal of Combinatorics},
  title        = {{Randomly perturbed digraphs also have bounded-degree spanning trees}},
  doi          = {10.37236/13316},
  volume       = {33},
  year         = {2026},
}

@misc{21960,
  abstract     = {Solitons - localized wave packets that travel without spreading - play a central role in understanding transport and properties of nonlinear systems. In quantum many-body systems, however, such robust excitations are typically destroyed by thermalization. Here, we theoretically demonstrate the existence of solitonic excitations in high-energy states of Rydberg atom chains in the regime of strong nearest-neighbor Rydberg blockade. 
These localized wave packets propagate directionally atop a special class of reviving initial states related to quantum many-body scars and are capable of carrying energy. Exhibiting long coherence times, these states constitute a form of non-ergodic quantum dynamics and can be efficiently implemented on Rydberg atom simulators. In this work, in addition to a phenomenological description of solitons, we identify their counterpart in a classical nonlinear dynamical system, demonstrate their potential use in quantum information transfer, and conjecture their relevance for anomalous energy transport reported in numerical studies of Rydberg atom arrays.},
  author       = {Kerschbaumer, Aron},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Research Data: "Quasi-solitons in Rydberg atom chains"}},
  doi          = {10.15479/AT-ISTA-21960},
  year         = {2026},
}

