2026-10-20 | | Total: 179
The gut microbiota shapes host metabolism, development and immunity, yet how the host regulates its composition in the presence of complex ecological interactions remains poorly understood. Recent experiments that measured T cell reactivity to individual strains in germ-free mice inoculated with a defined microbial community found that T cell recognition across microbial strains is sparse but broad, with a few polyspecific clones recognizing antigens shared across taxonomically diverse strains. Reanalyzing these experiments, we show that these T cell responses are associated with abundance changes of the strains they target. To interpret these findings, we develop a model linking the T cell repertoire to microbial abundances, in which a repertoire is allocated to minimize microbe-driven harm at the cost of maintaining an immune response. The model rationalizes the observed sparsity of gut T cell reactivity and shows that this allocation can emerge from competition between T cells for microbe-derived antigens. Finally, we find that while immune recognition can regulate microbial abundances, ecological interactions between strains may interfere with this mechanism. Together, these results provide a quantitative framework for immune control of gut microbial composition and identify principles governing the organization of T cell repertoires in the gut.
Dual-phase xenon time projection chambers (TPCs) are the leading technology for WIMP dark matter searches and have also demonstrated competitive sensitivities to neutrinoless double beta decays (0$νββ$) of \ce{^{136}Xe}. However, the large energy disparity between these two signals requires different detector optimizations. In particular, 0$νββ$\ searches in liquid xenon desire efficient rejection of closely located multiscatters of background gamma rays near the 0$νββ$\ energy, a capability achieved only in the vertical dimension in existing dual-phase xenon TPCs. We propose a highly customized, wide-angle optical system that can be conveniently integrated into large xenon TPCs and provide mm-scale transverse position resolution for MeV signals by imaging gas electroluminescence light. This approach preserves a TPC's low-energy sensitivity and optimizes a single experiment for both dark matter and 0$νββ$\ searches. This paper discusses the design, simulated performance, and implications of the imaging system for future liquid xenon-based 0$νββ$\ searches.
The ground state of a molecule coupled to an optical cavity contains virtual photons, and squeezing of the photonic field has been found in exactly solvable models, perturbation theory, and numerical calculations. Here, two exact sum rules for the field quadratures show that, within the dipole approximation, the vector-potential fluctuations of a linearly polarized mode never exceed their vacuum value in the ground state of any molecule at any coupling strength. Consequently, the field is squeezed whenever the mode couples to a dipole-allowed transition from a nondegenerate ground state, and the cavity-induced ground-state energy shift cannot increase with the cavity frequency at a fixed coupling vector, which is an exact constraint on approximate methods. Moreover, the incoherent photon statistics, which exclude the coherent part of the field, are never sub-Poissonian. At weak coupling, an exact resolvent expression reduces the squeezing and the ratio of the anomalous to the normal photon moment to linear-response properties of the isolated molecule, with the ratio set by a mean excitation energy. Time-dependent density functional theory with CAM-B3LYP and PBE0 reproduces exact and QED coupled-cluster references for both quantities to within 5% at 45 benchmark points. A survey of 29 molecules gives mean excitation energies of 9.5-21 eV at a cavity frequency of 2 eV and shows that the most polarizable molecules squeeze the field the most.
Pair-loaded relativistic turbulence is expected in compact high-energy astrophysical environments and may sustain plasma heating. Using self-consistent three-dimensional radiative particle-in-cell simulations, we identify the mechanism of ion heating mediated by pair-enhanced radiative feedback in relativistic turbulence. We derive the theoretical criterion for the transition to a regime dominated by radiative feedback and show that pair loading promotes the growth and persistence of large, connected optically thick structures. These long-lived structures sustain intense localized ion heating, repeatedly generating prominent hot-ion clumps and substantially enhancing the mean ion temperature. Our results establish a physical link between pair creation and sustained ion heating, providing a theoretical basis for investigating intermittent hot-ion production in black-hole coronae and gamma-ray burst outflows.
Foundation machine-learned interatomic potentials (MLIPs) promise broad transferability, yet their out-of-domain reliability is hard to establish. Here, we use pure and O/H-containing disordered carbon spanning multiple phases, competing hybridisations and porous frameworks to stress-test six popular foundation MLIPs. We first screen all six on density functional theory (DFT) energy and force errors and binding curves, then subject four to melt-quench-anneal molecular dynamics and evaluate the resulting structural observables. Binding curves and O/H-containing trajectories reveal unphysical extrapolation and dynamical failures that low DFT energy and force errors fail to reveal. Of the models tested, MACE-MH-1 is the most consistently robust, with low DFT errors and stable MD. Despite differing architectures, OMAT-based models share structural biases, highlighting the importance of training coverage. All four models generate structures of nanoporous carbon that reproduce the principal experimental and specialised carbon potential signatures. Incorporating O/H leaves short-range order largely unchanged but shifts pore-size distributions towards larger pores. These results emphasise the need for application-specific validation and extend molecular dynamics modelling of nanoporous carbon towards experimentally realistic structures.
With energy-related carbon dioxide emissions remaining at high levels, decarbonizing commercial building operations across key economies is critical for achieving global climate targets. This study assesses future operational carbon trajectories, mitigation potential, and the driving factors of this potential in commercial building operations across 16 countries from 2024 to 2060 by combining a Kaya-based emission model and Monte Carlo uncertainty analysis. The key findings indicate that (1) under the business-as-usual scenario, China's commercial building operational carbon emissions are projected to peak at 1150 (+-120) million tonnes of carbon dioxide (MtCO2) in 2033 (+-5 years), while the remaining countries exhibit differentiated carbon lock-in levels, with the highest level in the United States at 470 (+-89) MtCO2. (2) The total emission reduction potential increases from 352 MtCO2 in 2030 to 1119 MtCO2 in 2050 and stabilizes at 1123 MtCO2 in 2060, with China, the United States, and Japan contributing 84.7%-88.9% of the total potential. (3) By 2060, the mitigation potential from energy intensity reduction is concentrated mainly in China and the United States, reaching 52 and 32 MtCO2, respectively; end-use decomposition further identifies space heating and appliances as key targets for decarbonizing commercial building operations. In summary, this study characterizes future operational carbon trajectories and clarifies the drivers and end-use priorities of mitigation potential in commercial building operations, providing cross-country evidence for differentiated strategies for long-term commercial building decarbonization.
Discrete adjoint methods give the sensitivities that aerodynamic design and data assimilation require, but they remain hard to obtain for flow solvers that run on graphics processors. This paper presents CuAD, a source-to-source automatic differentiation tool that reads CUDA kernels and writes their tangent and adjoint versions as ordinary CUDA source. Its reverse mode is organised around the GPU thread and exploits the structure of flow-solver kernels: pure accumulations need no tape, the remaining intermediate values are recomputed rather than stored, and adjoint updates are atomic only where several threads reach the same memory. We evaluate CuAD on four kernels from simple flow models: a viscous Burgers residual, a two-dimensional advection-diffusion residual, a quasi-one-dimensional Euler residual and a Biot-Savart summation over point vortices. All generated derivatives agree with finite differences and pass adjoint dot-product tests to round-off, while the adjoint kernels cost two to six times the primal kernels. Compared with the compiler-based tool Clad under the same compiler, CuAD gives identical results at identical cost on the stencil kernels and is 16 to 470 times faster on the summation kernel because it stores no tape. A steady Burgers study shows that in single precision the adjoint gradient of a finite-difference Jacobian reaches errors of 45 to 75 per cent once the condition number grows to $2.5\times10^{5}$, whereas the CuAD Jacobian keeps the error below $5\times10^{-5}$. Four applications complete the evaluation: adjoint fields of the Burgers shock, the sensitivity of a scalar sensor in a cellular vortex, the inverse design of a transonic nozzle whose area is recovered to $2\times10^{-8}$, and an unsteady adjoint of the roll-up of a vortex sheet over 400 time steps.
Strong light-matter coupling has been proposed to control singlet-exciton fission, yet reported effects remain controversial due to the prevalence of cavity-induced optical artifacts in transient absorption measurements. Here, we investigate singlet fission in dielectric microcavities containing rubrene by probing in the near-infrared region, outside the cavity stopband, to eliminate mode-filtering artifacts. Across contrasting cavity quality factors ($Q = 90$ and $1150$) and a range of exciton-photon detunings ($+20$ to $-110$ meV), we find that the underlying kinetics are virtually identical to uncoupled films. Although minor spectral changes occur, transfer matrix simulations reveal these stem purely from mirror interference rather than modified photophysics. Transient electron paramagnetic resonance further confirms that high-spin triplet dynamics remain unperturbed. Together, these findings establish a widely applicable experimental framework for evaluating molecular photophysics in microcavities, revealing that singlet fission remains unaltered in strongly coupled rubrene.
Artificial-intelligence systems face a growing data-movement bottleneck and require higher interconnect bandwidth without a proportional increase in electrical power consumption or the number of laser sources and fibre connections. Moving optical interfaces closer to computing chips reduces electrical transmission loss, while demanding optical engines that combine high channel capacity, broadband modulation and simple packaging. Although wavelength multiplexing relieves pressure on the fibre budget, integrating shared multi-wavelength sources with dense wavelength routing and broadband modulation in a compact format remains challenging. Here we demonstrate a monolithic transmitter that brings these functions together on a single material platform, enabled by the low birefringence and native Pockels response of thin-film lithium tantalate. The device separates, independently modulates and recombines eight wavelength channels between one optical input and output. This architecture retains a two-fibre interface as the channel count grows, directly addressing the packaging burden of parallel optical links. Measurements across the eight channels establish an aggregate line rate of 1.6 Tbit s$^{-1}$. Operation with a shared optical frequency comb achieves an estimated net data rate of 0.58 Tbit s$^{-1}$ across four channels, demonstrating compatibility with a common multi-wavelength source. Scaling this architecture into multiple wavelength-multiplexed engines offers a path towards 10 Tbit s$^{-1}$ optical interfaces with high bandwidth per unit package edge, positioning thin-film lithium tantalate for the next generation of co-packaged optical interconnects.
The heat pathways in 3D-ICs have long been a labyrinth lacking experimental investigations into spatial distribution of thermal properties, severely limiting study of internal thermal management inside chips. For the first time, we established a methodology for 3D thermal characterizations of advanced 3D-ICs, implemented through a systematic workflow that enables layer-resolved thermal property extraction and the construction of a location-specific thermal database. We demonstrated its application to an advanced sample featuring 2.2-μm-diameter TSVs and ~2-μm-pitch face-to-face hybrid bonding (HB). Time-domain and steady-state thermoreflectance were used to extract the thermal conductivities (TCs) of Al traces, TSV, BEOL interconnects and HB pads, indicating HB structure enables efficient thermal transport with 9 W/m-K TC. Further, we mapped and visualized the TC distribution of TSV structures and RDL Cu wires with small step sizes down to 0.2 μm. Utilizing these results, we proposed an experiment-based thermal model to reveal that co-evolution of interconnect and substrate TCs is essential for thermal design and fabrication towards effective cooling optimization.
We derive discrete moment equations exactly equivalent to the lattice Boltzmann method (LBM) and identify the general moment closure shared by kinetic methods based on the discrete-velocity Boltzmann equation. This closure provides a route to correct lower-order dynamics without requiring accurate reconstruction of the higher-order moments, thereby substantially broadening the scope of existing closure strategies that seek increasingly accurate representations of these moments. These advances establish a unified moment framework bridging LBM and moment methods. The framework enables efficient simulations across all Knudsen regimes, overcoming the prohibitive complexity of existing moment methods at high Knudsen numbers. It further provides a common basis for integrating and extending key features of diverse kinetic methods.
Microwave imaging could monitor breast tumor response to neoadjuvant chemotherapy, but between-visit changes in breast positioning and normal tissue may confound tumor-related changes. We use eight patient-specific dielectric models from longitudinal MRI of two patients at four time points. For each model, 24-port S-parameters were simulated from 1 to 3 GHz, and Jacobians computed for tumor permittivity and loss, fibroglandular tissue (FGT) permittivity, and breast rotation and translation. 1 mm positioning changes produced changes 20 to 50 times larger than a 10% tumor permittivity change above 2 GHz. Nevertheless, tumor sensitivity patterns were weakly aligned with positioning patterns, tumor-FGT similarity exceeded 0.8 at only some visits, and the Jacobians remained well conditioned. Tumor responses are thus small but separable; monitoring will require positioning correction, frequency selection, and noise-aware detectability analysis.
Traveling wave (TW) structures are widely used in linear or circular accelerators for beam acceleration. Unlike the standing wave cavities, the overall accelerating voltage cannot be measured by the probes installed in the TW structures. This creates difficulties measuring and regulating the accelerating voltage and phase felt by the beam when it passes through the structure. In this paper, a model of TW structures is developed, allowing us to estimate the accelerating voltage and phase from the measurement of the RF drive power input to the TW structure. The beam loading effects, i.e., the beam-induced accelerating voltage and phase in the TW structure, are also included in the model. With the TW structure model, control methods are proposed to stabilize the accelerating voltage and phase felt by the beam based on feedback, feedforward, and optimization. The validity of the TW structure model is confirmed with the test results at the SwissFEL with improved beam stability under the same RF feedback controller. This work provides a firm basis for developing the multiple-bunch operation strategy for SwissFEL required by the future upgrade with more beamlines.
A limited understanding of how evaporation couples to crystallization pathways hampers control over the formation of solution-processed perovskite photoactive layers. Identifying these mechanisms-and learning how to select and tune them-is essential for achieving the desired crystal morphology and, ultimately, improving solar-cell performance and stability. To investigate crystallization pathways in evaporating thin films, we develop a simple and general zero-dimensional population-balance model that couples evaporation, chemical reactions, and crystallization. We use the model to explore a range of parameters for a solvent-solute system that follows a two-step crystallization pathway: solvent evaporation first induces formation of a crystalline precursor, which subsequently recrystallizes into the final phase. The model predicts the evolution of the film composition during drying and provides quantitative LaMer curves. Based on these results, we introduce a vocabulary of "supply limitations" and "rate limitations" that enables the interpretation of crystallization kinetics in terms of the governing process conditions.
In conventional coupled oscillator lattices and laser arrays, expanding the active volume inevitably causes spectral crowding and degrades spatial uniformity, posing a long-standing barrier to high-power single-mode operation. We uncover a dissipative selection mechanism in sparse star networks that decouples modal selectivity from lattice size constraints. By engineering dissipative auxiliary couplings between a central hub and N peripheral resonators, the system establishes an exact star-Laplacian symmetry. This enforces a target eigenstate that is immune to all lossy channels while maintaining equal amplitude across all main nodes. All competing eigenstates break this symmetry, driving the auxiliary loss and opening a persistent, size-independent modal loss gap that prevents multi-mode competition. Realistic microring simulations validate the robustness of this dark-state isolation under electromagnetic dispersion and mechanism extension to arbitrary arrays is allowed when meets the system's exceptional point. By exploiting graph topology to bypass multi-mode competition, this work establishes a general route toward scalable high-power single-supermode laser arrays and robust state selection in coupled non-Hermitian systems.
This work presents the development and assessment of a hybrid framework coupling numerical simulation, experimental samples, data assimilation (DA) and machine learning (ML) to investigate the flow within an axial compressor stage. The methodology first implements a Multi-Fidelity Ensemble Kalman Filter (MFEnKF) which integrates experimental PIV data with a two-fidelity numerical simulation setup, combining steady-state RANS on a fixed grid with unsteady hybrid RANS/LES simulations on a moving grid. This data assimilation step successfully reduces discrepancies between experiments and CFD realizations. Besides, it provides additional high fidelity samples when compared with the available PIV data. This enriched database is used to train machine learning (ML) tools, with the objective to obtain black-box forcing terms to be included in the discretized form of the Navier-Stokes equations. These models are expected to preserve the accuracy of the DA results while not longer requiring the input of experimental data. The neural network architecture used here maps local physical features such as the velocity, the velocity gradients and the turbulent viscosity, to obtain the volume forcing terms. The analysis of the performance of such ML models shows that global accuracy is preserved but localized increased errors in specific regions of the compressor, such as close to the blade tips, can exceed those observed from the baseline RANS simulations. To resolve this issue, a Mixture of Experts (MoE) method is employed, which classifies and combines the trained ML tools for different flow regions. The discrepancy with experimental data of this method is systematically lower than the one observed for the baseline RANS solution and very similar to the one obtained with the MFEnKF.
Modularity maximization is widely used to identify community structure in complex networks. In many applications, however, community detection must also incorporate prior knowledge and domain-specific requirements beyond network topology. We study constrained modularity maximization with restrictions on the number of communities, must-link and cannot-link relations, and general resource constraints. While fast heuristics can provide high-quality solutions, they generally lack optimality guarantees, whereas exact mathematical-programming approaches can become computationally demanding as the network size and constraint complexity increase. We develop a unified optimization framework offering complementary heuristic and exact solution methods. These include a constrained Louvain heuristic for feasible solution generation, a column-generation heuristic with graph-partitioning-based pricing, and exact branch-and-price algorithms using mixed-integer linear programming and sparse mixed-integer quadratic programming formulations for the pricing subproblem. We develop a cutting-plane procedure providing independent upper bounds to assess solution quality when exact optimization cannot be completed. Experiments on 11 unconstrained networks with up to 1,005 nodes and 30 constrained instances with up to 379 nodes demonstrate the complementary roles of these methods. The column-generation heuristic rapidly improves the constrained Louvain solutions across the tested constraint classes, while branch-and-price using the sparse quadratic formulation solves 29 of the 30 constrained instances to optimality within a one-hour time limit. This formulation becomes particularly effective on larger, sparse, undirected networks by substantially reducing the size of the exact pricing model. Overall, the framework guides method selection to balance computational time, solution quality, and optimality guarantees.
Purpose: Current eye models often assume a fixed spherical retinal curvature of 12 mm, which does not reflect the anatomical variability observed in vivo. This study proposes an extension incorporating the stochastic relationship between axial length and three-dimensional (3D) retinal curvature, enabling more realistic 3D eye models. Methods: Axial length and 3D retinal curvature data from 94 subjects, split into five subsets for cross-validation, were used to construct two probabilistic models: a Single Gaussian Model (SGM) and two-component Multi Gaussian Mixture Model (MGMM). Both were converted into conditional forms to estimate retinal curvature for a specific axial length, and evaluated against measured data distributions. Results: All models produced radii of curvature comparable to measured data (median difference <0.2 mm). Non-conditionally, the MGMM described the distribution slightly better than SGM, with higher density overlap (>94% vs. >86%), but similar variance (worst f-test (p-value) 1.05 (0.71) vs. 1.05 (0.73)). Conditioned on axial length, both models performed similarly with overlap exceeding 93% and the worst f-test 1.02 (0.84). Given its simplicity and robustness, the conditional SGM was selected, achieving similar validation results: median difference (<0.10 mm), overlap (>92%), and a non-significant worst f-test (1.02 (0.85)). Conclusions: The conditional SGM enables accurate generation of anatomically realistic 3D retinal curvatures based on axial length, supporting improved stochastic eye modeling for visual optics and population-based simulations of healthy eyes. Translational relevance: Replacing the fixed 12 mm retina assumption with a population-based stochastic model improves anatomical accuracy of eye models used in optical design and treatment planning.
Predictive modeling of inertial-confinement-fusion (ICF) implosions increasingly requires resolving ion-kinetic effects beyond the reach of radiation-hydrodynamics, without sacrificing radiation transport. We present, to our knowledge, the first coupling of a kinetic-ion plasma model to thermal radiation transport. The iFP code - a conservative, spherically symmetric, 1D2V hybrid kinetic-ion/fluid-electron Vlasov-Fokker-Planck solver - is extended with a grey P1 radiation model, coupled through a high-order/low-order (HOLO) fixed-point solver and advanced on a shared adaptive moving mesh. Opacities are mixed on the fly from the LANL TOPS library, exploiting the per-species information intrinsic to a kinetic code. We verify the implementation against the optically thick Marshak-wave and Mach-45 radiative-shock benchmarks, and apply it to a NIF polar-direct-drive exploding-pusher capsule. Increasing physics fidelity systematically improves agreement with the measured neutron yield and burn width; the improvement is driven by kinetically enhanced fuel leakage and ballistic preheating, and radiation transport proves essential for quantitative agreement. This work establishes a proof-of-principle capability for higher-fidelity kinetic-plasma/radiation modeling of ICF.
High-accuracy computation of nonsmooth dynamics is a challenging problem in systems with multiple frictional contacts. This study proposes a general physics-informed neural network framework, termed the mixed leading variable physics-informed neural network (MLV-PINN). By using the governing equations and contact-friction laws as physical constraints, this framework provides common solution approaches for computing nonsmooth dynamics not only in discrete systems but also in continuum systems. Specifically, the network represents accelerations, normal contact forces and friction margins as independent outputs and determines them simultaneously, without first reducing the dynamics to discrete velocity-impulse complementary subproblems. A composite physics-informed loss is constructed from the equations of motion in their original form, contact-friction complementary conditions and a friction-margin consistency relation. To improve stick-slip state identification during PINN optimisation, a mixed leading variable combining relative tangential velocity and acceleration is introduced to reformulate the tangential complementary conditions. Four benchmark simulations cover one-dimensional (1-D) nonsmooth dynamics involving single-contact stick-slip, 1-D nonsmooth dynamics involving multi-contact stick-slip, two-dimensional (2-D) nonsmooth dynamics with coupled normal contact-separation and tangential stick-slip, and 2-D nonsmooth dynamics in a continuum system involving moving contact on an elastic beam. Benchmark comparisons show substantial reductions in displacement and velocity errors relative to the corresponding conventional methods at large time steps. These results, together with an ablation study and tests involving friction-velocity slopes, long-time integration and loss weights, illustrate the computational accuracy and support the broad applicability of the proposed framework.
The establishment of metrics and control of fundamental noise sources in quantum sensing is a long-standing goal in quantum science and technology. We report on the observation of polariton-mediated optical noise in a thermal Rydberg atom vapor and demonstrate its active control and suppression using a radio-frequency field in Rydberg electrometry. The noise exhibits structured amplification and reduction across the Rydberg electromagnetically induced transparency resonance. This reveals a nonlinear slow-light feedback mechanism in which the polariton medium stores and releases probe-field energy. We show that the reference radio-frequency field used in Rydberg heterodyne receivers can also tune the medium into an operating regime where high electromagnetic-field responsivity coincides with reduced optical noise. This operating-point control enables optimization of the electric-field sensitivity within the constraints imposed by Rydberg-polariton physics. These results establish polariton noise as a previously unrecognized intrinsic noise source in thermal Rydberg vapors and the corresponding control method as a crucial aspect of continued efforts to achieve fundamental sensitivity limits in Rydberg electrometry.
We describe a class of electromagnetic Gaussian modes with spherical spatiotemporal structure that remains invariant during propagation through anomalously dispersive media. The angular profile of the modes is given by spherical harmonics, whose quantization axis can be freely chosen. As a result, the modes support spatiotemporal orbital angular momentum (OAM) configurations with longitudinal and transverse components set by the orientation of the quantization axis. Superpositions of these modes generalize the familiar solution of a chirped pulse that simultaneously focuses and compresses to arbitrary angular structure. The modes are derived through a series of unitary transformations that map the dispersive paraxial wave equation onto a Schrödinger equation with a 3D isotropic harmonic potential. These transformations facilitate the identification of conserved operators for mixed longitudinal and transverse OAM.
In this study, we investigate the feasibility of accelerating small probes up to relativistic energies via magnetocentrifugal acceleration by pulsars. Specifically, we consider a broad range of pulsars and examine probes with radii ranging from $10$ to $1000$ nm. We show that the probes may attain Lorentz factors of up to approximately $25$. Furthermore, we investigate the detectability of such probe swarms and demonstrate that, if a highly advanced civilization deploys probes with a total mass of approximately one billion tonnes for acceleration, their collective synchrotron emission could be detectable by the FAST radio telescope at a distance of 1 kpc. The study also discusses the potential resemblance between the emission from such probe swarms and fast radio bursts (FRBs), and outlines directions for future research.
Memory capacity and bandwidth are becoming major constraints in neural-network computing, motivating parameter-efficient representations that reduce model storage and data movement. Tensor-network methods already enable substantial compression of neural networks, but their structured operations are often poorly matched to conventional processors optimized for dense linear algebra. Such representations can be implemented in photonic hardware, although this typically requires optical architectures tailored to a particular factorization. Here we propose a different approach based on programmable photonic meshes. Compact electronic parameters determine the interferometer phases, while optical interference produces a dense, generally high-rank transformation. Since the parameterization remains electronic, its dimensionality and structure can be changed without modifying the optical circuit. Numerical experiments in end-to-end training demonstrate that the proposed photonic layers match or outperform strongly compressed tensor-train and brickwork layers, often with fewer trainable parameters. Experiments in parameter-efficient fine-tuning of a language model also show advantages over established digital adapters. The reduced trainable state lowers persistent memory requirements and offers the potential to reduce data movement and bandwidth demands.
The interaction between human infrastructures and external events results in loss of functionality or even complete system disruption through a snowball effect of cascading-failures. Their complex interplay challenges their prevention, affecting the human population that rely on their functionality. Understanding how these networks respond to external perturbations is crucial for predicting and mitigating the impacts of disruptive events. In particular, meteorological events, while not preventable, can have significant impacts on human infrastructure. Here, we propose a model coupling the dynamics of transportation systems to the dynamics of external events considered as fields, to capture at the same time the critical behavior of the former and the level of influence of the latter. The impact of rainfalls on real-world railway networks provides a compelling case study, showcasing local endemic and system-wide disruptions around a critical value of the rainfall strength. Remarkably, localized perturbations on the wider network reveals the local weaknesses of the global infrastructure, and we assess them in terms of system perturbability and corresponding societal risk. Our framework opens the door to the quantitative analysis of criticalities in essential human infrastructure imputable to vulnerabilities that are set to intensify due to climate change.