2026-09-30 | | Total: 175
Reliable multistep quantum lattice-Boltzmann evolution requires controlling computational faults. When ideal collision preserves encoded mass and momentum exactly, changes in those charges can provide syndromes for selected faults at the post-collision checkpoint. A fault that shifts one population then leaves a unique mass-momentum residual labeled by the lattice velocity. A coherent reference register stores the expected charges, so the same test remains valid while the physical fluid charges vary across the lattice. Single bit flips on population number registers occupy disjoint syndrome sectors and satisfy the Knill-Laflamme condition. A state-vector demonstration on eighteen data qubits of the two-dimensional nine-velocity lattice (D2Q9) recovers that single-bit family to double-precision roundoff. Conservation still leaves a fifteen-dimensional kinetic nullspace on the three-dimensional nineteen-velocity lattice (D3Q19). Exact multiple-relaxation-time (MRT) streaming analysis ranks those charge-preserving modes by the order at which they return to density or momentum. A classical D3Q19 decaying-flow simulation at moderate Reynolds number then injects about $1.75\times10^4$ selected single-population shifts per realization. Recovery of that alphabet returns the trajectory to floating-point roundoff. The quantum statements assume an ideal charge-preserving collision and a reliable reference. The result is an inner recovery map for a stated charge-changing alphabet together with a kinetic classification of the unresolved sector.
Multimode optomechanical scattering is ubiquitous in free-space imaging protocols, but its radiation pressure counterpart is challenging to detect with extended objects. Here we demonstrate radiation pressure actuation of a nanomechanical membrane using dynamically structured light from a spatial mode demultiplexer (SPADE). Using single- and two-mode illumination, we directly image the position-dependent intermodal coupling between Hermite-Gauss (HG) modes as mediated by the membrane's motion. As an application, we use multimode radiation pressure feedback to simultaneously cool two nearly degenerate membrane modes. Our results are applicable to diverse mechanical systems and serve as a stepping stone to SPADE-based quantum optomechanics.
We introduce and numerically demonstrate receiver-informed transmitter (RITX) waveform design for 200 Gb/s PAM-4 C-band direct-detection PON. Using 50 GHz optoelectronics, RITX achieves N1/E1 loss classes with FFE-Free ONU, enabling 40 km reach and 48.1%/66.7% network-wide complexity savings at 1:64/1:128 splits
A simple model gives good agreement between observed emissions of CO2 into the atmosphere from the combustion of fossil fuels and observed increases of the atmospheric concentrations C of CO2 and decreases of the isotope delta value delta13C. If fossil fuel emissions were to cease the model predicts that: (1) the difference between the CO2 concentration C and the preindustrial value, Co = 280 ppm, would decay exponentially with time constant of 60 y; (2) the difference between the concentration-delta-value product C x delta13C and the preindustrial value Co x delta13Co = (280 ppm) x (-6.5%%), would decay exponentially with a time constant of 8 y. The decay times are different because of isotopic exchange, largely driven by photosynthesis and respiration on the land.
We report the experimental observation of the $^{\text{133}}$Cs 6S$_{\text{1/2}}$ to 6D$_{\text{5/2}}$ single-frequency two-photon transition at low excitation power using an optical nanofiber embedded in a cold atom cloud. We investigate the 917 nm fluorescence spectra modified by the surface-induced interaction of the nanofiber, which exhibits a pronounced asymmetric broadening with a distinct red-shifted tail. We present a systematic theoretical description of the two-photon transition process of atoms near a nanofiber surface, considering atomic angular momentum coupling, guided-mode decay rate, and surface-induced van der Waals interactions. Owing to the tight spatial confinement provided by the guided evanescent field, efficient nonlinear excitation is realized at excitation powers down to tens of microwatts, corresponding to a reduction of approximately four orders of magnitude relative to typical free-space implementations. Our work reveals the interplay among waveguide-modified radiative dynamics, surface-induced interactions, and the spatial distribution of atoms around the nanofiber, and provides insights into further studies on nonlinear optical transitions and surface-mediated atomic dynamics in waveguide quantum electrodynamics systems.
Emergence, features of a whole that no part shows alone, is central to complex systems but rarely measured so as to locate its source. Viewing a system structurally, as a network of interacting parts, we evaluate emergence as what observing the coupled whole reveals beyond, or erases from, the combined observations of the parts. For observations that build what they see route by route, we show that every emergent feature is produced by its own route across the interface between parts, a route that traces it to the components and interactions responsible. We show the number of all such boundary-crossing paths bounds the total emergence discrepancy, and equals the emergence of measures that count routes. The bound is computed locally around the interface, validated through network examples, and in a nervous system the routes show that interneurons relay most of the emergence. Our work associates emergence with path structure in the network system that can be anticipated, attributed, and engineered.
Turbulent flows are highly chaotic, which makes their instantaneous states difficult to predict. Data assimilation aims to reduce this predictive uncertainty by synthesizing the predictions of a physical solver with partial observations of the system. Traditional data assimilation methods such as ensemble and variational approaches can incur substantial computational cost through evaluations of high-fidelity solvers and, for variational methods, adjoint calculations. Recently, generative diffusion models have been explored in the solution of inverse problems in fluid mechanics, demonstrating their potential for data assimilation. However, existing diffusion-based methods often use free parameters instead of the prescribed observation-error covariance to control the balance between the prior and likelihood, departing from the Bayesian formulation of data assimilation. We propose a framework for principled generative data assimilation in which observations are assimilated in the clean-state space and the prescribed observation-error covariance enters explicitly through an ensemble Kalman update. In the Lorenz-63 system, the proposed method yields a lower trajectory error and equation residual than the other diffusion-based posterior sampling methods considered, although its ensemble underestimates the posterior uncertainty. In two-dimensional Rayleigh-Bénard convection, the generated fields recover the statistics of the direct numerical simulation reference at large scales, with deviations confined to the smaller scales. These results suggest that the framework can produce physically plausible turbulent flow fields, indicating its potential for more challenging turbulent flow applications.
Network reconstruction from partial information is usually performed at the same resolution level at which constraints are observable. This becomes problematic when only coarse-grained information is available, while the relevant process occurs at a finer scale. Here we employ the multiscale model of weighted networks introduced in a companion paper and turn it into a probabilistic framework for reconstructing weighted networks across arbitrary aggregation levels. The model is built to preserve its functional form under coarse-graining, so that global parameters calibrated on an observable aggregate layer can be transferred to finer layers without refitting. We test the method on two empirical systems with different aggregation mechanisms. In the International Trade Network, countries are aggregated into geographic macro-regions and the observed coarse-grained layer is used to infer the underlying country-level network. In the Dutch production network, sectoral flows are reconstructed across the hierarchical industrial classification, using coarser sectoral layers to infer finer ones. In both geographical and sectoral settings, we benchmark our genuinely multiscale reconstruction method against a state-of-the-art weighted reconstruction model calibrated directly at the target resolution, and therefore using additional information available at the same (finer) scale at which performance is evaluated. Remarkably, despite this informational disadvantage, our method recovers the fine-scale binary structure with high accuracy, improving precision, specificity, accuracy and maximum degree-error diagnostics, while remaining nearly equivalent in sensitivity.
A droplet breaks up in tens of microseconds, and a recording captures perhaps a dozen frames. The states in between cannot be recovered without repeating the experiment, and simulating them is too costly to sweep an operating envelope. Yet they are present in the corpus as a whole: a campaign spanning a device's actuation range produces morphologies resembling those any single recording missed. Exploiting that requires a representation that is low-dimensional, invertible, and faithful to shape rather than to sampling. Proper orthogonal decomposition supplies the first two but measures distance in sampled coordinates; manifold learning supplies the third but no map back to a shape; elastic shape analysis supplies a shape metric but no reduced coordinates. SHROM composes all three. Interfaces are represented by their square-root velocity functions, a neighbour graph is built over that shape space, and an autoencoder is trained to reconstruct while penalising latents in which graph neighbours are not latent neighbours. The demonstration uses 301,539 inkjet droplet contours and four filmed break-up sequences. The graph term does not improve reconstruction. It determines whether position in the latent carries meaning: clustering the latent of an otherwise identical model recovers the shape-space partition at chance level (ARI = 0.063 +/- 0.059), and at 0.781 +/- 0.049 with the term active. Waveform parameters predict the full contour at R^2 = 0.878 +/- 0.026. Negative results are reported in the same terms. The graph metric proved immaterial across five choices, and interpolation error saturates at the reconstruction limit, so a plain autoencoder leads that task. The main limitation is interpolation across a topology change: no component of the regularised latent holds both a single-component and a post-break-up shape, whereas an unregularised autoencoder mixes them freely.
We introduce chiral rotational wavepackets: coherent superpositions of a few low-lying rotational states of an achiral molecule for which the expectation value of a molecular axis traces a three-dimensional chiral curve in the laboratory frame. The handedness is a dynamical property of the quantum state rather than of the equilibrium structure, and we show that it can be imprinted on linear rotors as well as on spherical and symmetric tops. We derive the action of space inversion, motion reversal, and their product on such wavepackets in the Wigner $D^J_{MK}$ basis, and define the rotational enantiomers of a given wavepacket. To quantify the handedness of the orientation trajectory we introduce two derivative-based pseudoscalar measures: a geometrical measure built from the curvature and torsion of the trajectory, which is even under motion reversal, and an origin-referenced measure built from the binormal vector, which is odd under motion reversal. The latter is shown to be proportional to the geodesic curvature of the trajectory projected onto the unit sphere, which links it, via the Gauss--Bonnet theorem, to the solid angle enclosed by the projected trajectory and to the parallel-transport holonomy of a molecule-fixed vector. Finally we show that a nonzero geometrical measure is a sufficient condition for the wavepacket to be truly chiral in the sense of Barron, while wavepackets with real expansion coefficients are only falsely chiral.
Optical anapoles have attracted significant interest as a promising platform for enhancing light-matter interactions at the nanoscale. An ideal anapole is a localized excitation of electromagnetic fields that does not produce any far-field radiation. In practice, anapoles are not ideal due to the parasitic scattered fields associated with higher-order multipoles. Here, we use the recently established exact current multipole expansion and the exact anapole condition to design an ideal anapole in which radiation due to the multipole orders up to the electric and magnetic octupoles is either eliminated or strongly suppressed, reducing the total radiated power by orders of magnitude. We also propose and study realistic photonic structures in which such anapoles can be excited by a uniform linearly polarized plane wave through second-harmonic generation that involves the off-diagonal elements of the second-order susceptibility tensor. In this case, the anapole condition is reached off-resonantly by adjusting the polarization angle of the incident field, which allows for the anapoles to be excited over a wide spectral range in scatterers of different sizes. Using this approach, we achieve scattering suppression on the order of $10^{-4}$ in photonic structures that are half a wavelength in size. The method offers the possibility to realize localized electromagnetic excitations that are truly isolated from the surroundings.
Multicomponent methods treat both electrons and nuclei quantum mechanically using the standard machinery of electronic structure theory. Previous work has demonstrated that quantitatively describing electron-nuclear correlation in multicomponent many-body perturbation theory methods requires large electronic basis sets with high-angular momentum atomic orbitals, and calculations on large systems can therefore become impractical. Similarly, the slow convergence of the electron-electron correlation energy with respect to highest atomic-orbital angular momentum is a well-known issue in standard electronic structure theory in which only the electrons are treated quantum mechanically, and numerous explicitly correlated methods have been introduced to mitigate this poor convergence. Motivated by the success of single-component explicitly correlated methods, we generalize the single-component MP2-F12/3C(FIX) approach for electron-electron correlation to electron-proton correlation to reduce the need for large electronic basis sets in multicomponent calculations. We also implement the complementary auxiliary basis set (CABS) singles correction for the electrons and protons. Our calculations show that the multicomponent MP2-F12 electron-proton correlation energy is within a few percent of the complete-basis-set (CBS) estimate obtained by QZ-5Z inverse-cubic extrapolation with PB4-D held fixed at every electronic basis-set cardinality tested, with the aug-cc-pVDZ basis set recovering 96-98% of the extrapolated CBS electron-proton correlation energy and exceeding the conventional multicomponent MP2 value at the aug-cc-pV5Z level. These results indicate that explicitly correlated methods can substantially reduce the electronic basis-set requirements of multicomponent calculations.
Adsorption and desorption redistribute soluble surfactant between a liquid film and its free surface, but cannot change the total inventory in a closed periodic domain. We audit this requirement in nonlinear weighted-residual models for soluble-surfactant-laden falling films. Although a previous correction restores conservation of the nonlinear surface Marangoni flux, a separate inconsistency remains in the depth-integrated bulk diffusion. We derive the resulting mass imbalance and construct explicit positive periodic concentration fields that exhibit it. The missing contribution arises from the streamwise concentration gradient appearing in the normal diffusive flux at a sloping interface. Including this contribution restores local conservation of the coupled bulk-surface transport system while leaving the hydrodynamic equations, the leading-order adsorption-desorption closure, and the linear stability properties of the uniform equilibrium unchanged. In an independent no-slip calculation using the parameters of an earlier travelling-wave study, the reference model loses $0.0228832\%$ of its initial surfactant mass by $t=4500$, whereas the corrected run has an absolute relative mass change of $1.2\times10^{-15}$. The main hump and capillary ripples remain closely comparable at $t\approx4500$. This benchmark therefore shows that agreement with travelling-wave calculations does not, by itself, guarantee that a reduced model respects the conservation laws inherited from the underlying transport problem. Although the drift is small in this case, the corrected equations conserve the reduced total inventory for every smooth periodic solution. We also distinguish this exact conservation property from higher-order compatibility with the full curved-interface boundary conditions.
The compressible MHD energy cascade in the solar wind is shaped by Alfvénicity, plasma fluctuations, and proton temperature, yet how these factors regulate the relative weight of compressible and incompressible energy transfer across heliocentric distances remains poorly understood. We apply compressible and incompressible exact relations to a large sample of 2-hour intervals from WIND (1~au), Parker Solar Probe ($<0.25$~au), and Solar Orbiter (0.3--1~au), classifying intervals as Alfvénic ($|σ_c|>0.75$) or non-Alfvénic ($|σ_c|<0.25$) and by wind speed, and examine the Yaglom-like and purely compressible components of the compressible cascade. Velocity fluctuation amplitude is the strongest predictor of the compressible cascade rate $\langle|\varepsilon_c|\rangle$, followed by density and magnetic field fluctuations. $\langle|\varepsilon_c|\rangle$ and $\langle|\varepsilon_i|\rangle$ are tightly correlated in all regimes, with their ratio organized by Alfvénicity and wind speed: non-Alfvénic and Alfvénic fast winds cluster along $\langle|\varepsilon_c|\rangle=\langle|\varepsilon_i|\rangle$, while Alfvénic slow wind shows a systematic excess consistent with internal energy contributions. The largest $\langle|\varepsilon_c|\rangle$ systematically occur in hotter intervals, with a two-decade increase from cool to hot events in Alfvénic slow and non-Alfvénic wind, a trend carried mainly by the purely compressible term, while the dominant Yaglom-like term controls the overall cascade magnitude. Alfvénicity and wind speed thus govern the efficiency of turbulent energy transfer across heliocentric distances, and the tight coupling between compressible and incompressible rates points to the need for compressible, multi-fluid descriptions of the solar wind dissipation budget.
We present an equivalent-circuit analysis that simultaneously considers the photoluminescence and electrical responses of valley photovoltaic devices. Valley photovoltaics (VPV) are a novel concept for hot-carrier solar cells, where intervalley scattering of electrons in the absorbing material helps maintain the conduction-band electron populations at metastable satellite valleys, which are at energies higher than the conduction-band minimum. We measure intensity- and voltage-dependent photoluminescence and electrical signals of a VPV device, which has previously shown S-shaped current-voltage characteristics and thus low efficiency. We present an equivalent circuit that captures the S-shaped current-voltage curves and simultaneously describes our photoluminescence and electric measurements. The equivalent-circuit model indicates that a reverse diode causes inefficient extraction of electrons and thus the S shape. We use a Poisson/drift-diffusion model to demonstrate that the reverse diode can be either explained by the valley-scattering process itself or a heterojunction barrier. The reverse diode must be eliminated before VPV devices can achieve high efficiency.
Correlations between optical modes, ranging from classical fluctuations to quantum entanglement, are a cornerstone of modern optics and emerging quantum technologies. However, probing these correlations on the natural timescale of electronic motion - the attosecond regime - remains a major challenge. Here, we generate and characterize correlated attosecond extreme-ultraviolet (XUV) emission by perturbing gas-phase high-harmonic generation with two-mode bright squeezed vacuum. The signal and idler fields imprint their correlated amplitude and phase fluctuations onto two families of harmonics, thereby transferring these correlations from the infrared to broadband XUV modes. Even harmonics serve as an intrinsic attosecond phase-sensitive probe of these correlations through interference between indistinguishable XUV pathways. Shot-resolved covariance measurements reveal the amplitude correlations, while even-harmonics resolve the phase correlations with sub-cycle precision. Time-domain reconstruction of the correlated emission reveals a transition from maximal fluctuations to near-silence within only 30 as. Whereas the signal and idler driven attosecond pulse trains exhibit excess fluctuations individually, their joint reconstruction shows strongly suppressed fluctuations, reflecting their shared correlations. These results establish a framework for generating and probing quantum correlations on attosecond timescales, opening the door to quantum attosecond science.
Arrays of neutral atoms excited to Rydberg levels have emerged as one of the most promising platforms for quantum computation and simulation. With the hope to outperform classical devices, the number of atoms has been increased by orders of magnitude. However, the interaction time, i.e., the maximum accumulated time during which an atom interacts with its neighbours, has been limited to a few microseconds only. This restrains the number of gates per atom or prevents the simulation of long-time dynamics of quantum many-body systems. Here, we observe the spin-exchange interaction between two spin 1/2s encoded in laser-trapped circular Rydberg atoms over more than $60μs$ and 40 spin-oscillation periods, improving the state of the art by an order of magnitude. These unprecedented timescales allow us to record for weak atomic trapping a collapse and revival of the spin-exchange oscillation contrast induced by spin-motion coupling. Our results constitute the first observation of this coupling with laser-trapped Rydberg atoms. To counteract its detrimental effect for quantum simulation, we demonstrate a novel dynamical decoupling method that prevents the collapse of the spin oscillations. This method exploits the advantages of our recently developed hybrid platform that enables the measurement and the optical manipulation of long-lived circular Rydberg atoms with auxiliary Rydberg atoms. This work opens a direct route to long-duration quantum simulation of strongly-interacting many-body systems.
Molecular mixing models for the conditional scalar diffusion (CSD) term in the filtered mass density function (FMDF) transport equation used in Large Eddy Simulation (LES-FMDF) approach are evaluated against direct numerical simulation (DNS) of spatially evolving jets. A database of DNS of turbulent mixing slot jets of $C_{7}H_{16} / O_2$ and $C_{7}H_{16} / Air$ is developed. The formulation includes the compressible form of the governing equations, a generalized multicomponent diffusion model, a cubic real gas equation of state and realistic property models. Simulations are conducted over a wide range of initial pressures ($1 atm < P_0 < 100 atm$) and jet width based Reynolds numbers of 850 and 1300. FMDF of mixture fraction at various filter widths are obtained from the simulation at several spatial locations within the flow. The CSD term in the exact transport equation of FMDF is calculated from the DNS. An \textit{a priori} analysis of Interaction by Exchange of Mean (IEM), Modified Curl (MC) and Mapping Closure (MAPPING) mixing models for CSD is conducted. At high pressures, significant deviations are observed between the CSD measured from the DNS and those predicted by the models. The source of the deviations is determined to be the mixing frequency used in the models. The significance of the mixing frequency and its variation with change in pressure and diffusion models is studied. The primary objective of this work is to assess the applicability of aforementioned mixing models to flows at large pressures and determine the causes for deviation in prediction.
We introduce an exact algebraic architecture that evaluates an arbitrary finite Fourier series using a two-mode ($2 \times 2$) linear optical circuit, with the only tunable components being single-mode phase shifters encoding the function argument. We prove that such a circuit must exist for every Fourier series and derive an analytical method for its construction based on spectral factorisation. The resulting optical system exhibits an $\mathcal{O}(N)$ depth for an $N$-harmonic expansion, executing function evaluations in the passive optical time-of-flight. Finally, we validate our claims numerically, demonstrating that even for sequences with thousands of Fourier terms, our proposed circuit construction correctly synthesises continuous and discontinuous nonlinear functions. Our architecture thus provides a universal, deterministic foundation for single-variable nonlinear optical computing on integrated photonic platforms.
Second-quantized neural-network quantum states have achieved accurate molecular energies, but extending them across molecular geometries requires a shared representation of the geometry-dependent wavefunction coefficients. We introduce geometry-conditioned foundation neural-network quantum states for molecular electronic structure in second quantization. A single autoregressive model learns a family of ground states from sparse anchor geometries and provides wavefunctions at untrained geometries without further optimization. Orbital alignment matches orbital identities and transports their phases, establishing an aligned orbital basis across geometries. Frozen energies reach chemical accuracy at every untrained query geometry for N$_2$, CO, and H$_4$. On additional molecular paths, the energy-trained wavefunctions yield dipoles, quadrupoles, and natural occupations without property labels. Across three paired N$_2$ training seeds, orbital alignment lowers the mean absolute energy error over all untrained query geometries from 34-37 mHa to 0.049-0.085 mHa. At approximately 1 mHa mean absolute error, frozen evaluation reduces the per-geometry cost by $986\times$ relative to independent optimization, yielding an estimated $25.8\times$ end-to-end GPU-cost reduction on a 161-point N$_2$ grid.
We investigate the transient dynamics of cavity-EIT-based all-optical switching in a system of trapped ions coupled to an optical cavity through numerical simulations. In contrast to steady-state analysis, the time-dependent response provides direct insight into the switching speed, transient dynamics, and achievable switching contrast. We consider three distinct switching schemes and systematically compare their performance as a function of the relevant system parameters. The switching contrast is evaluated from the time-dependent cavity output and used to characterize the performance of each scheme. For the four-level N-type system, we find distinct trade-offs between switching through suppression of the cavity-EIT signal and switching through resonance shifting, with the former exhibiting a faster temporal response. The corresponding three-level scheme can achieve near-unity switching contrast with substantially shorter switching pulses, although it does not provide an independent switching control field. These results establish the performance limits and trade-offs of different cavity-EIT-based switching schemes and provide guidelines for optimizing their operation for applications such as high-speed optical gating, frequency-selective photon routing, frequency-multiplexed quantum communication, etc.
A PIN-diode-reconfigurable quasi-Yagi antenna fed by a coplanar waveguide (CPW) is presented for X-band and Ku-band frequency reconfiguration. PIN diodes are strategically introduced between the driven and parasitic elements to control the effective surface-current paths and thereby modify the operating frequency characteristics. By optimizing the diode locations, number, and switching states, two distinct operating modes are achieved. In Mode I, dual-band operation covering the X-band and Ku-band is obtained, with simulated impedance bandwidths of 9-11 GHz and 12.2-18.8 GHz. In Mode II, a simulated continuous impedance bandwidth from 10 to 14.2 GHz is achieved. The corresponding simulated peak gains are approximately 8 and 7 dBi for Mode I and Mode II, respectively. A prototype is fabricated and experimentally characterized in terms of its S-parameters, and the measured results confirm the frequency-reconfiguration behavior of the proposed antenna. The effects of the number and placement of the PIN diodes on the impedance characteristics are further investigated to clarify their roles in frequency reconfiguration and design complexity. The proposed antenna provides a compact CPW-fed quasi-Yagi configuration with PIN-controlled X-band and Ku-band frequency reconfiguration.
The stimulated Raman shortcut-to-adiabatic passage (STIRSAP) technique has been demonstrated to accelerate adiabatic population transfer processes. Traditionally applied to three-level atomic systems, this method enables rapid and efficient population transfer in a short interaction time. In this work, we extend this approach to a four-level N -type system driven by three laser fields. Under the three-photon phase matching condition, the first-order Doppler effect can be completely eliminated, making this configuration more suitable for applications in frequency metrology and quantum information processing. By adiabatically eliminating the intermediate state, the four-level system is effectively reduced to a two-level model, enabling the application of counterdiabatic driving without requiring additional couplings. This approach allows for fast and efficient population transfer between metastable states. We analyze the influence of laser intensity peaks and detuning, and we show that the transfer time is significantly shorter than that achieved with the STIRAP technique. Moreover, we also show 1 that an optimal coupling strength between the ground and metastable states further minimizes the operation time.
We review recent theoretical and empirical advances that recast fish schooling within a statistical-physics framework, emphasizing how heterogeneous and time-dependent interactions govern collective motion. Building on extensions of Vicsek-type models to complex and weighted social networks, we show that topology and link strength qualitatively alter flocking stability and critical thresholds, with empirical weights typically reducing global alignment. High-resolution trajectory inference reveals selective, nonreciprocal responses: individuals preferentially attend to faster neighbors, producing transient, speed-induced leadership rather than fixed hierarchies. At the mesoscopic scale, schools exhibit avalanche-like turning cascades with scale-free size and duration statistics, aftershock clustering, and an Omori-type temporal decay with a short memory. Together, these findings support a unified picture in which collective order and critical-like fluctuations coexist: alignment provides stability while near-critical variability preserves responsiveness. This synthesis highlights fish schools as a tractable experimental model for how living collectives organize, transmit and process information across scales.
The Laboratory for High Energy Physics and the Institute of Anatomy at the University of Bern are establishing a pre-clinical proton therapy research facility at the Bern University Hospital (Inselspital). The facility centers on the Bern Medical Cyclotron, an 18 MeV compact cyclotron designed for microampere-range beam currents for radioisotope production. The Bern Medical Cyclotron can be adapted as a highly accessible and customizable platform for pre-clinical proton therapy research, enabling advanced modalities such as ultra-high dose rate (FLASH) radiotherapy and spatially fractionated radiation therapy. Human keratinocytes (HaCaT) and mouse melanoma (B16-F10) cell lines were irradiated with 8.14(29) MeV protons at conventional dose rates and with 225 kV X-rays at doses from 0 to 8 Gy. Clonogenic survival assays were analyzed using the linear-quadratic model, and relative biological effectiveness values were derived by comparing the proton and X-ray doses required for 0.1 survival. Proton irradiation yielded steeper survival curves than X-rays, indicating higher cytotoxicity. B16-F10 cells were most sensitive, with a relative biological effectiveness of 1.34 at 0.1 survival, while HaCaT cells showed a relative biological effectiveness of 1.21. These results confirm enhanced proton effectiveness in cell killing and demonstrate that the Bern Medical Cyclotron provides a reliable platform for radiobiological studies. Establishing in vitro proton irradiation capability marks a key milestone toward developing a comprehensive pre-clinical radiobiology research facility in Bern. Ongoing upgrades will further refine dosimetry and proton delivery, paving the way for future FLASH and proton spatially fractionated radiation therapy studies.