Fluid Dynamics
Papers
Fresh from arXiv, typeset properly. Highlighted words trade on the market.
Saturday 17 January
Modelling fluid dynamics with machine learning (ML) has advanced rapidly, yet most data driven approaches remain opaque because they rely on complex architectures to capture nonlinear flow behaviour. This lack of interpretability limits the reliability and hinders the understanding of when and why some models succeed or fail. To address this, we present a transparent approach that provides insights into how data-driven fluids dynamics and ML work. This is achieved by training a convolutional neural network (CNN), on data from a simple laminar fluid flow, to behave as an operator that exactly matches the finite-difference numerics, providing a direct link between well established theory and this new world of ML models. Importantly, the model demonstrates strong generalisation capability by reproducing the dynamics for a wide range of distinct and unseen flow conditions within the same flow category. The CNN learns the forward Euler three-point stencil weights, capturing physical principles such as consistency and symmetry despite having only three tuneable weights. This interpretable ML model goes beyond pure numerical training (numCNN),the approach is shown to work when trained on analytical (anCNN) and even molecular dynam ics (mdCNN) data. In some cases, the physics is not captured, and thanks to the simple and interpretable form, these CNNs provide insight into the limits, pitfalls and best practice of data-driven fluid models. Because the approach is based on finite-difference operators, it naturally extends to many structured-grid computational fluid dynamics problems, including turbulent, multiphase and multiscale flows as well as systems beyond the continuum such as molecular dynamics.
Chemical Physics
Full Reaction Pathway Dynamics for Atmospheric Decomposition Reactions: The Photodissociation of H$_2$COO
Branching ratios for fragmentation channels of important meta- and unstable species are essential for a molecular-level characterization of atmospheric chemistry. Here, the molecular product channels for the decomposition dynamics of the smallest Criegee intermediate, H$_2$COO, are quantitatively investigated. Using a high-quality, full-dimensional machine learned potential energy surface (CASPT2/aug-cc-pVTZ), the translational, rotational, and vibrational energy distributions of the CO$_2$+H$_2$, H$_2$O+CO, and HCO+OH fragmentation channels were analyzed to elucidate partitioning of the available energy. The CO$_2$ + H$_2$ product forms through two different pathways that bifurcate after formation of the OCH$_2$O intermediate. Along the direct pathway, CO$_2$ is preferentially vibrationally excited with H$_2$in its vibrational ground state, whereas for the indirect pathway going through formic acid, H$_2$ can populate levels with $v > 0$. For all product channels passing through energized formic acid, the lifetime distributions are described by stretched exponentials with $\beta$ ranging from 1.1 to 1.7. This is a clear signature of non-RRKM effects and suggests that the explicit molecular dynamics needs to be followed for a quantitative and realistic description of the photodissociation dynamics.
Read the full abstract
Physics and Society
Structure of Pitch-Pattern Motifs in Major League Baseball
Baseball consists of two teams alternating between batting and fielding while competing to score runs through sequential pitching events. Recent advances in tracking technology have enabled all Major League Baseball (MLB) clubs to record every pitch with high resolution, yet most quantitative studies have primarily emphasized single-pitch metrics, leaving the role of sequential structure less explored. Here, we examine pitch-pattern motifs of multiple lengths using approximately 12.4 million Statcast pitch recordings from the 2008-2025 MLB regular seasons at two complementary scales. At the macroscale, we quantify pitch-sequence diversity using the Shannon entropy and inverse Simpson index and examine their relationships with earned run average and wins. At the microscale, we compare hit and out frequencies across pitch-pattern motifs. Rather than identifying outcome-determining sequences, we find that motif usage exhibits stable, non-random organization, as reflected in Zipf's and Heaps' laws, while showing limited association with conventional performance measures. While language-like scaling (Zipf's and Heaps' laws) clearly reveals an underlying 'grammar' of MLB pitch sequences, that grammar alone is insufficient to account for performance indicators such as earned run average or wins. These results suggest that sequence-based analyses clarify the structural organization of pitch usage, while also delineating the limits of motif-based approaches for explaining performance without richer contextual information.
Read the full abstract
Chemical Physics
Towards accurate predictions of bond-selective fluorescence spectra
Vibrational-encoded fluorescence spectro-microscopies are emerging as powerful tools for studying molecular vibrations with the unparalleled sensitivity of fluorescence spectroscopy. We recently described one such technique, termed bond-selective fluorescence-detected infrared-excited (BonFIRE) spectro-microscopy. Currently, prospects of BonFIRE towards rational molecular design are limited, but they have the potential to be assisted by computational tools. In this Perspective, we provide a brief overview of the theory of BonFIRE spectroscopy. We then describe a fully automated computational pipeline for calculating BonFIRE spectra, reproducing key features of experimental results. Finally, we highlight a few potential applications of computational methods for vibrational-encoded fluorescence spectro-microscopies and their broader implications for chemistry and biology.
Read the full abstract
Chemical Physics · Cross-listed from cond-mat.mtrl-sci
Transition Metal Dichalcogenide MoS${}_2$: oxygen and fluorine functionalization for selective plasma processing
Low-temperature plasma processing is a promising technique for tailoring transition metal dichalcogenides (TMDs). For chalcogen substitution processing, a key challenge is to identify the ion energy window that enables selective chalcogen removal while preserving the metal lattice. Using ab-initio molecular dynamics (AIMD), we demonstrate that oxygen and fluorine functionalization widen the processing window by significantly lowering the sulfur sputtering energy threshold ($E_{\text{sputt,S}}$) of MoS${}_2$ from $\sim 30$ eV to $\sim 10$ eV via formation of sputtering products such as SO${}_2$ and SF${}_n$. Additionally, we show that experimentally relevant cryogenic temperatures strongly affect $E_{\text{sputt,S}}$. The dependence is confirmed via AIMD and also predicted by a mechanistic parameter-free theory, suggesting that $E_{\text{sputt}}(T)$ generalizes to other TMDs, functionalization, and surface impacts in general. Our results highlight oxygen/fluorine functionalization, ionic impact angle, and material temperature to be key control parameters for selective, damage-controlled chalcogen removal in TMD processing.
Read the full abstract
Optics
Revealing the long-range coupling for multi-dimensional metasurface multiplexer
Metasurface coupling constitutes a fundamental yet intricate electromagnetic interaction that occurs within a lattice of artificial subwavelength unit cells. Despite its prevalence, such coupling is typically ignored in conventional metasurface design frameworks due to the high characterization complexity, leading to suboptimal device performance. Here, we reveal a distinctive long-range coupling that exceeds an order of magnitude compared with the interaction range of evanescent waves, substantially enriching the metasurface design landscapes. This coupling exhibits pronounced graph topological features, and we design a graph neural network (GNN) to accurately abstract its inherent physics. Through strategic enhancement of the coupling effects, the discrete metasurface responses are transformed into continuous states, thereby unlocking diverse multiplexing channels. By further integrating the GNN into an inverse design agent, we tailor the multi-channel global response of metasurface to support simultaneous multiplexing across angle, frequency, and polarization domains. Experimentally, we demonstrate a compact metasurface multiplexer with eight independent channels, showcasing its potential for next-generation vehicular networks. This work establishes a new paradigm for highly integrated multifunctional metasurfaces, with promising prospects for high-density optical storage, information encryption, and high-capacity wireless communication.
Read the full abstract
Optics
Plasmonic Bi-Cavity Nanostructure for Efficient Light Collection and Localization
Tip-enhanced Raman spectroscopy (TERS) typically relies on high-NA excitation to generate a strong axial field at the tip apex, which shortens the working distance and constrains sample geometries. We show that a plasmonic bi-cavity tip, the plasmon-tunable tip pyramid (PTTP), co-tuned in nanopyramid length L and plateau length W, supports a hybrid antenna-cavity mode that funnels energy to the apex under radially polarized, on-axis excitation, even with a dry objective of NA = 0.75. Finite-element simulations identify W as a design-critical parameter that sets an in-plane surface-plasmon-polariton (SPP) Fabry-Pérot-like resonance; co-tuning (L,W) yields a periodic series of maximal apex |E|^2. Experiments on monolayer graphene confirm near-field enhancement and reproduce the characteristic annular TERS point spread function (PSF) with NA = 0.75. Relaxing the NA requirement increases working distance and compatibility with constrained environments, pointing to practical, deployment-ready nano-Raman instrumentation.
Read the full abstract
Instrumentation and Detectors
Study of the energy calibration of the DEAP-3600 detector using $^{22}\mathrm{Na}$ source data and simulations
DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3.3 tons of LAr contained in a spherical acrylic vessel. At a WIMP mass of 100 GeV, DEAP-3600 has a projected sensitivity of $10^{-46}\,\mathrm{cm}^{2}$ for the spin-independent elastic scattering cross section of WIMPs. Radioactive sources have been used for the energy calibration and to test the detector performance. One of the most effective calibration runs has been taken with a $^{22}\mathrm{Na}$ source deployed in a tube located around the DEAP-3600 steel shell. The simultaneous emission of three $\gamma$'s by the source provides an excellent tagging for the $^{22}\mathrm{Na}$ decay. The results concerning the energy response of the detector and the agreement between data and Monte Carlo simulations in DEAP-3600 are investigated in this study.
Read the full abstract
Optics
Integrated nano electro-optomechanical spiking neuron
Neuromorphic computing offers a pathway toward energy-efficient processing of data, yet hardware platforms combining nanoscale integration and multimodal functionality remain scarce. Here we demonstrate a gallium-phosphide electro-optomechanical spiking neuron that integrates optical and electromechanical interfaces within a single nanostructure on a silicon photonic chip operating at telecommunication wavelengths (1550 nm) and exploiting a 3 gigahertz-frequency mechanical mode. Our device displays excitable dynamics, generating optical spikes at its output, as in the spiking activity of neurons and cardiac cells and defined by the calibrated all-or-none response to external perturbations. This dynamic is consistent with the saddle-node on invariant circle scenario and associated features are demonstrated including control of excitable threshold, temporal summation and refractory period. Our device compact footprint and its CMOS-compatible platform make it well suited for edge-computing applications requiring low latency and establish a foundation for versatile brain-inspired optomechanical computing and advanced on-chip optical pulse sources.
Read the full abstract
Chemical Physics
Extension of the CIPSI-Driven CC($P$;$Q$) Approach to Excited Electronic States
We extend the CIPSI-driven CC($P$;$Q$) methodology [K. Gururangan et al., J. Chem. Phys. 155 (2021) 174114], in which the leading higher-than-doubly excited determinants are identified using the selected configuration interaction (CI) approach abbreviated as CIPSI, to excited electronic states via the equation-of-motion (EOM) coupled-cluster (CC) formalism. By examining vertical excitations in CH+ at equilibrium and stretched geometries, adiabatic excitations in CH, and ground- and excited-state potential cuts of water, we demonstrate that the CIPSI-driven CC($P$;$Q$) method converges parent CC/EOMCC singles, doubles, and triples energetics from relatively inexpensive Hamiltonian diagonalizations in CI spaces smaller than the corresponding triples manifolds.
Read the full abstract
Instrumentation and Detectors
High-Voltage Performance Testing in LAr of the PMMA Cathode Connection for the DarkSide-20k Experiment
DarkSide-20k (DS-20k) is a next-generation dual-phase liquid argon (LAr) time projection chamber (TPC) devoted to the direct-detection of dark matter. The detector is currently under construction in Hall-C at the Laboratori Nazionali del Gran Sasso, Italy, at a depth of approximately 3500 m water equivalent. The detector will instrument 49.7 t of low-radioactivity underground LAr contained within an acrylic TPC and is designed to reach a WIMP-nucleon spin-independent cross-section sensitivity down to $10^{-48}\,\mathrm{cm}^{2}$ for a WIMP mass of $0.1\,\mathrm{TeV}/c^{2}$ in a 200 tonne-year run. In DS-20k a uniform electric drift field is established in the active volume to transport ionization electrons toward the electroluminescence region, with the required high voltage delivered to the TPC cathode through a custom cable and stress-cone assembly. At the University of California, Davis, a dedicated test setup was developed to reproduce the DS-20k cathode high-voltage connection in LAr, matching the local electric-field conditions. This work summarizes the results of a comprehensive test campaign validating the operation of the DS-20k cathode HV system in LAr up to $-100$ kV.
Read the full abstract
Friday 16 January
Physics and Society
Critical Transit Infrastructure in Smart Cities and Urban Air Quality: A Multi-City Seasonal Comparison of Ridership and PM2.5
Public transit is a critical component of urban mobility and equity, yet mobility and air-quality linkages are rarely operationalized in reproducible smart-city analytics workflows. This study develops a transparent, multi-source monitoring dataset that integrates agency-reported transit ridership with ambient fine particulate matter PM2.5 from the U.S. EPA Air Quality System (AQS) for four U.S. metropolitan areas - New York City, Chicago, Las Vegas, and Phoenix, using two seasonal snapshots (March and October 2024). We harmonize heterogeneous ridership feeds (daily and stop-level) to monthly system totals and pair them with monthly mean PM2.5 , reporting both absolute and per-capita metrics to enable cross-city comparability. Results show pronounced structural differences in transit scale and intensity, with consistent seasonal shifts in both ridership and PM2.5 that vary by urban context. A set of lightweight regression specifications is used as a descriptive sensitivity analysis, indicating that apparent mobility-PM2.5 relationships are not uniform across cities or seasons and are strongly shaped by baseline city effects. Overall, the paper positions integrated mobility and environment monitoring as a practical smart-city capability, offering a scalable framework for tracking infrastructure utilization alongside exposure-relevant air-quality indicators to support sustainable communities and public-health-aware urban resilience.
Read the full abstract
Physics and Society
D-MODD: A Diffusion Model of Opinion Dynamics Derived from Online Data
We present the first empirical derivation of a continuous-time stochastic model for real-world opinion dynamics. Using longitudinal social-media data to infer users opinion on a binary climate-change topic, we reconstruct the underlying drift and diffusion functions governing individual opinion updates. We show that the observed dynamics are well described by a Langevin-type stochastic differential equation, with persistent attractor basins and spatially sensitive drift and diffusion terms. The empirically inferred one-step transition probabilities closely reproduce the transition kernel generated from the D-MODD model we introduce. Our results provide the first direct evidence that online opinion dynamics on a polarized topic admit a Markovian description at the operator level, with empirically reconstructed transition kernels accurately reproduced by a data-driven Langevin model, bridging sociophysics, behavioral data, and complex-systems modeling.
Read the full abstract
Physics and Society
The Impossibility of Cohesion Without Fragmentation
Accounts of social division and cohesion presuppose relations whose maintenance already depends on positional compatibility. We develop this feasibility layer as a static theory: a positioning event fixes players' positions, and a compatibility function determines which pairs can sustain a relation. Two results follow from one formalism. First, any non-uniform compatibility pattern on the dyads sharing a positional requirement forces fragmentation and cohesion together; the only events that guarantee no collapse are those rendering every pair compatible. The asymmetry between unilateral severance and bilateral confirmation recovers the convention underlying pairwise stability without a behavioral premise, while a separate specialization, in which compatibility partitions the position space, yields assortative survival without choice or opportunity. Second, at the measure level, the survival rate factorizes exactly into integration and cohesion terms. Along any complete relaxation path that relaxes one pair at a time (the generic case under distance thresholds with absolutely continuous positions), the cohesion term is non-monotone once at least three players are present, and in clustered configurations it is bounded above by a sharp asymptotic value that finite systems approach but never attain, with strict crossing of the two terms once the number of clusters exceeds two.
Read the full abstract
Physics and Society
Meeting the Needs of the Global Quantum Science Community: A Call to Action
2025 marks one hundred years since the discovery of quantum mechanics. In the century since then, quantum science has blossomed into a global community composed of academics, engineers, developers, and entrepreneurs. The world is currently in the middle of the so-called second quantum revolution, with increased public awareness of quantum science and technology, and growing investment in both quantum hardware and software applications. However, representation remains low among historically marginalized groups: women, LGBTQ+, BIPOC, and people from the global south make up disproportionately few physicists. There are numerous efforts to improve diversity within quantum science, including through workforce development. But many of the changes enacted at the highest levels have failed to result in real change, as highlighted and discussed in the recent Women For Quantum Manifesto of Values. Here, we seek to echo and amplify the need for real change in the quantum ecosystem, emphasizing intersectionality and a feminist approach that centers the most vulnerable members of the quantum community: young students and researchers, especially those communities historically marginalized from quantum science. This report is our attempt to help quantum communities meet this need; we have conducted a survey of quantum scientists all over the world, and here we include both a preliminary report of our findings and policy suggestions we have built to address them. The primary results of our survey are that, 1) marginalized quantum scientists are experiencing hardships and challenges more than their more privileged peers across all metrics, 2) that this fact is hurting retention of diverse, talented quantum scientists in our field, and 3) quantum EDI is an investment in talent retention and resilience building, which are essential for a thriving, globally competitive quantum ecosystem.
Read the full abstract
Optics
Sub-Doppler rubidium atom cooling using a programmable agile integrated PZT-on-SiN resonator
Programmability and precise control of laser frequency are essential for quantum experiments and applications such as atomic clocks, quantum computers, and cold-atom sensors. Current systems use bulky, power-hungry modulators and frequency shifters which are difficult to integrate and limit portability and scalability. We report an electrically controllable, agile optical frequency source based on a semiconductor laser stabilized to a photonic-integrated, lead zirconate titanate (PZT)-actuated resonator cavity. We demonstrate this approach with precision programmable frequency control of a 780-nm laser that can periodically reference to rubidium spectroscopy followed by fast, programmable, arbitrary frequency tuning sequences for quantum control. We use this approach to demonstrate sub-Doppler cooling of rubidium-87 without any external modulators, achieving atom-cloud temperatures as low as 16 $\mu$K. The device achieves a tuning strength up to 1 GHz/V with 11 MHz modulation bandwidth while consuming only 10 nW of electrical power. This work establishes a route toward compact, low-power, and chip-scale laser systems for next-generation quantum and atomic sensing technologies.
Read the full abstract
Fluid Dynamics · Cross-listed from math.NA
Global Recovery from Local Data: Interior Nudging for 2D Navier-Stokes equations in a Physical Domain
In many real-world applications of data assimilation (DA), the strategic placement of observers is crucial for effective and efficient forecasting. Motivated by practical constraints in sensor deployment, we show that global recovery of the flow field can be achieved using observations available only in a subregion of the domain, possibly far from the boundary. We focus on the two-dimensional incompressible Navier-Stokes equations posed in a bounded physical domain with Dirichlet boundary conditions. Building on the continuous data assimilation framework of Azouani, Olson, and Titi (2014), we rigorously prove that the assimilated solution converges globally to the true solution under suitable conditions on the nudging parameter, spatial resolution, and the geometry of the observation region, specifically, when the maximum distance from any point in the domain to the observational subregion is bounded by a constant multiple of \( \nu^{1/2} \) (in terms of scaling). Our computational results, conducted via finite element methods over complex geometries, support the theoretical findings and reveal even greater robustness in practice. Specifically, synchronization with the true solution is achieved even when the observational subregion lies farther from the rest of the domain than the theoretical threshold permits. Across all three tested scenarios, the local nudging algorithm performs comparably to full-domain assimilation, reaching global accuracy up to machine precision. Interestingly, observational data near the boundary are found to be largely uninformative. This demonstrates that full observability is not necessary: carefully chosen interior observations, even far from the boundary, can suffice.
Read the full abstract
Chemical Physics
Time-dependent density functional theory study of strong-field laser-induced coulomb explosion of the HCl dimer
We present a channel-resolved interpretation of laser-driven Coulomb explosion of the HCl dimer from an ensemble of trajectories. Three dominant outcomes are identified: a minor three-body channel and two four-body channels (sequential and near-simultaneous dissociation of both molecules). The key result is that pathway selection is strongly correlated with the degree of ionization during the laser interaction, which is in turn strongly modulated by laser-molecule orientation. Higher early-time ionization predisposes the system toward near-simultaneous four-body breakup, whereas lower ionization favors sequential and three-body fragmentation; for low-ionization cases, a fragment-resolved charge metric further differentiates three-body and sequential behavior. These charge-dependent trends consistently map onto experimentally accessible observables: the simultaneous mechanism dominates the high-energy tail of the kinetic energy release (KER) spectrum and populate distinct regions of the emission-angle distributions, while sequential events concentrate at lower KER. Overall, early-time charge evolution provides a unifying explanation for channel branching and for the channel-resolved fragmentation signatures.
Read the full abstract
Optics
A wafer-scale ultrasensitive programmable chiroptical sensor
Chiroptical enantioselective sensing is gaining traction across various applications. However, intrinsic molecular chiroptical responses are weak, and existing amplification approaches add synthesis, manufacturing, or operational complexity that limits sensitivity, scalability, and dynamic control. Here, we present a fundamentally new sensing paradigm merging adsorption-driven chirality induction with wafer-scale optical transduction in a programmable heterostructure containing twisted aligned carbon nanotubes (CNTs) and phase change materials (PCMs). Chiral molecules adsorb onto CNTs to form chiroptically active composites that are macroscopically assembled by alignment and rotational stacking, yielding large ultraviolet circular dichroism (CD). We resolve molecule concentration and handedness in a single device without lithography, hotspot delivery, or differential protocols, achieving sub-$\mu$M sensitivity for CD-silent glucose and chiral amino acids enabled by $>10^5\,\mathrm{M^{-1}}$ adsorption constants. We validate adsorption using molecular dynamics simulations, reproduce experimental results using chiral transfer matrix simulations, and realize sensor programmability by tuning the PCM layer. This platform enables cost-effective in-situ enantiomer monitoring in aqueous environments.
Read the full abstract
Optics
Broadband THz spectroscopy system beyond 25 THz using BNA crystals and a tunable single-ring-fiber compressor
We present a terahertz (THz) time-domain spectroscopy (THz-TDS) system which accesses a broadband spectrum, efficiently covering the so-called "new THz gap" between 5 and 15 THz and extending beyond 25 THz. The system exploits nonlinear interactions within the organic crystal BNA (N-benzyl-2-methyl-4-nitroaniline) to generate and detect THz radiation upon excitation by a near-infrared (NIR) pulse centered at 1.03 $\mu$m. To enable broadband THz spectral monitoring, the NIR pulse from a Yb-based solid-state laser undergoes spectral broadening in a gas-filled single-ring hollow-core photonic-crystal fiber, followed by a pulse compression to achieve durations as short as 31 fs. This approach paves the way for broadband spectroscopy in hard-to-access THz regions using widely available near-infrared ultrafast sources.
Read the full abstract
Optics
Ultra-broadband Mid to Long-wave Infrared Spintronic Poisson Bolometer
Infrared detectors have traditionally been divided into two fundamental classes, mid-wave (MWIR, 3-5 um) and long-wave (LWIR, 8-14 um). Integrating MWIR and LWIR within a single device is challenging due to distinct materials, cooling needs, and detection mechanisms, while such integration is critical for improved object recognition, temperature estimation, and environmental sensing. In this work, we demonstrate a Spintronic Poisson (SP) bolometer enabling room-temperature ultra-broadband sensing across 3-14 um. Unlike conventional bolometers that rely on continuous analog signals, the SP bolometer implements a Poisson-counting detection paradigm, encoding temperature in discrete stochastic events, which turns thermal noise from a limitation into the basis of the estimator itself. We fabricate the SP bolometer using a spintronic transduction layer integrated with a plasmonic nanoantenna array to enhance broadband infrared absorption. Using spintronic transduction, the device achieves the noise-equivalent temperature difference (NETD, thermal sensitivity metric) of 80-100 mK at 300 K, surpassing uncooled detectors and approaching cooled technologies. This work establishes a statistical detection paradigm for room-temperature infrared sensing with broad application potential.
Read the full abstract
Medical Physics
Dose-LET Interactions Predict Capsular Contracture After Proton Postmastectomy Radiation Therapy
Pencil beam scanning (PBS) proton therapy provides highly conformal dose distributions that are increasingly leveraged for postmastectomy radiation therapy (PMRT) to reduce cardiopulmonary exposure. However, implant-based reconstruction in the setting of PMRT remains vulnerable to capsular contracture, and biological mechanisms of possible high linear energy transfer (LET) in PBS have not been well characterized. A retrospective case-control study was conducted on consecutive breast cancer patients who underwent mastectomy followed by implant-based reconstruction and proton PMRT (50 Gy in 25 fractions) between 2015 and 2021. Dose-LET volume histograms (DLVHs) were calculated for peri-implant tissue (5-mm shell around the implant). Generalized linear mixed-effects regression (GLMER) was employed to identify DLVH indices significantly associated with capsular contracture. Spearman correlation analysis was used to eliminate redundance. DLVCs were derived from receiver operating characteristic (ROC) analysis and validated using support vector machine (SVM)-based normal tissue complication probability (NTCP) model. Eight capsular contracture and 16 matched controls patients were analyzed. Three independent and significant DLVH indices were identified(p<0.01). The corresponding DLVCs were: V(55.8 Gy[RBE=1.1], 2.2 keV/{\mu}m) < 0.0033%, V(50.3 Gy[RBE=1.1], 5.4 keV/{\mu}m) < 0.0017%, and V(32.8 Gy[RBE=1.1], 0.9 keV/{\mu}m) > 96.98%. The SVM-based NTCP model achieved an area under the ROC curve (AUROC) of 0.867, with 91.7% accuracy, 87.5% sensitivity, and 93.8% specificity. Capsular contracture following proton PMRT is significantly associated with the synergistic interplay between dose and LETd in peri-implant tissue. The derived DLVCs provide actionable dosimetric constraints that can be integrated into treatment planning to minimize capsular contracture risk in proton PMRT.
Read the full abstract
Fluid Dynamics
Coupled two-phase flow and surfactant/PFAS transport in porous media with angular pores: From pore-scale physics to Darcy-scale modeling
Two-phase surfactant-laden flow and transport in porous media are central to many natural and engineering applications. Surfactants alter two-phase flow by modifying interfacial tension and wettability, while two-phase flow controls surfactant transport pathways and interfacial adsorption. These coupled processes are commonly modeled using Darcy-type two-phase flow equations combined with advection–dispersion–adsorption transport equations, with capillary pressure–saturation relationships scaled by the Leverett $J$-function. However, the Leverett $J$-function idealizes porous media as bundles of cylindrical tubes and decouples interfacial tension and wettability, limiting its ability to represent angular pore geometries and interfacial tension–wettability coupling effects. We present a modeling framework that explicitly incorporates pore angularity and interfacial tension–wettability coupling into Darcy-scale surfactant-laden flow and transport models. Two-phase flow properties are derived for angular pores, upscaled across pore size distributions, and formulated as explicit and closed-form expressions. These upscaled relationships are integrated into a coupled flow–transport model to simulate transient two-phase flow and surfactant transport. Results reveal a nonlinear and nonmonotonic dependence of two-phase flow properties on pore angularity, pore size distribution, and interfacial tension. Example simulations of water flow and PFAS migration in unsaturated soils indicate that surfactant-induced flow effects on PFAS leaching are generally minor under typical conditions, whereas pore angularity strongly controls water flow, interfacial area, and PFAS retention. Overall, the proposed framework provides a more physically grounded approach for modeling two-phase surfactant-laden flow and transport in porous media.
Read the full abstract
Instrumentation and Detectors
AllShowers: One model for all calorimeter showers
Accurate and efficient detector simulation is essential for modern collider experiments. To reduce the high computational cost, various fast machine learning surrogate models have been proposed. Traditional surrogate models for calorimeter shower modeling train separate networks for each particle species, limiting scalability and reuse. We introduce AllShowers, a unified generative model that simulates calorimeter showers across multiple particle types using a single generative model. AllShowers is a continuous normalizing flow model with a Transformer architecture, enabling it to generate complex spatial and energy correlations in variable-length point cloud representations of showers. Trained on a diverse dataset of simulated showers in the highly granular ILD detector, the model demonstrates the ability to generate realistic showers for electrons, photons, and charged and neutral hadrons across a wide range of incident energies and angles without retraining. In addition to unifying shower generation for multiple particle types, AllShowers surpasses the fidelity of previous single-particle-type models for hadronic showers. Key innovations include the use of a layer embedding, allowing the model to learn all relevant calorimeter layer properties; a custom attention masking scheme to reduce computational demands and introduce a helpful inductive bias; and a shower- and layer-wise optimal transport mapping to improve training convergence and sample quality. AllShowers marks a significant step towards a universal model for calorimeter shower simulations in collider experiments.
Read the full abstract
Fluid Dynamics · Cross-listed from cond-mat.soft
Universal wrinkling dynamics of a sheet on viscous liquid
We investigate the wrinkling dynamics of a thin elastic sheet that is indented or compressed while floating on a viscous liquid. We show that the deformation speed controls the dynamics, leading to a wrinkle wavelength significantly smaller than that selected under quasistatic compression. Once active compression ceases, the wrinkles coarsen until their wavelength relaxes toward the equilibrium value. We develop a theoretical model coupling Stokes flow in the liquid to elastic bending of the sheet, which quantitatively predicts both the initial wavelength selection and its subsequent coarsening. We demonstrate that the same mechanism governs two dimensional and axisymmetric geometries, thereby extending classical static wavelength selection laws to dynamic situations. Although developed from controlled laboratory experiments, the model captures a generic viscous-elastic coupling and applies broadly to thin elastic films interacting with viscous environments, including the formation of surface wrinkles in pahoehoe lava flows.
Read the full abstract
History and Philosophy of Physics · Cross-listed from cond-mat.stat-mech
100 Glorious Years of the Ising Model
This is an editorial article based on the reseaches on the Ising model over the last 100 years.
Classical Physics
Conformal Symmetry and the Thermal Effects of Acceleration in Classical Physic
An accelerating Rindler frame in Minkowski spacetime acting for a finite time interval is used to carry a box of particles or waves between two relativistic inertial frames. The finite spatial extent of the box allows treatment of the equations of motion for particles or for waves, while the Rindler acceleration provides a substitute for scattering to test for thermal equilibrium. In the case of equilibrium for relativist particles, the Juttner distribution is derived. For relativistic waves, a full derivation of the Planck spectrum including zero-point radiation is obtained within classical theory. For relativistic waves, relativistic behavior and conformal symmetry are crucial. It is emphasized that the classical two-point correlation function for classical zero-point radiation depends upon the geodesic separation between the spacetime points and is independent of the coordinate system choice. The classical point of view here does not give any support for the idea that a system in uniform acceleration through classical zero-point radiation finds a thermal system.
Read the full abstract
Applied Physics · Cross-listed from nlin.AO
Temporal Complexity and Self-Organization in an Exponential Dense Associative Memory Model
Dense Associative Memory (DAM) models generalize the classical Hopfield model by incorporating n-body or exponential interactions that greatly enhance storage capacity. While the criticality of DAM models has been largely investigated, mainly within a statistical equilibrium picture, little attention has been devoted to the temporal self-organizing behavior induced by learning. In this work, we investigate the behavior of a stochastic exponential DAM (SEDAM) model through the lens of Temporal Complexity (TC), a framework that characterizes complex systems by intermittent transition events between order and disorder and by scale-free temporal statistics. Transition events associated with birth-death of neural avalanche structures are exploited for the TC analyses and compared with analogous transition events based on coincidence structures. We systematically explore how TC indicators depend on control parameters, i.e., noise intensity and memory load. Our results reveal that the SEDAM model exhibits regimes of complex intermittency characterized by nontrivial temporal correlations and scale-free behavior, indicating the spontaneous emergence of self-organizing dynamics. Notably, such regimes arise over finite intervals of noise intensity rather than at a single critical point, consistent with the concept of extended criticality. Further, the noise intensity range needed to reach the critical region, where self-organizing behavior emerges, slightly decreases as the memory load increases. This study highlights the relevance of TC as a complementary framework for understanding learning and information processing in artificial and biological neural systems, revealing the link between the memory load and the self-organizing capacity of the network.
Read the full abstract
Chemical Physics
The rise and fall of stretched bond errors: Extending the analysis of Perdew-Zunger self-interaction corrections of reaction barrier heights beyond the LSDA
Incorporating self-interaction corrections (SIC) significantly improves chemical reaction barrier height predictions made using density functional theory methods. We present a detailed, orbital-by-orbital analysis of these corrections for three semi-local density functional approximations (DFAs) situated on the three lowest rungs of the Jacob's Ladder of approximations. The analysis is based on Fermi-Löwdin Orbital Self-Interaction Correction calculations performed at several steps along the reaction pathway from the reactants (R) to the transition state (TS) to the products (P) for four representative reactions selected from the BH76 benchmark set. For all three functionals, the major contribution to self-interaction corrections of the barrier heights can be traced to stretched bond orbitals that develop near the TS configuration. The magnitude of the ratio of the self-exchange-correlation energy to the self-Hartree energy (XC/H) for a given orbital is introduced as an indicator of one-electron self-interaction error. For the exact, but unknown density functional, XC/H = 1.0 for all orbitals, while for the practical DFAs studied here, XC/H spans a range of values. The largest values are obtained for stretched or strongly lobed orbitals. We show that significant differences in XC/H for corresponding orbitals in the R, TS, and P configurations can be used to identify the major contributors to the SIC of barrier heights and reaction energies. Based on such comparisons, we suggest that barrier height predictions made using the SCAN meta-generalized gradient approximation may have attained the best accuracy possible for a semi-local functional using the Perdew-Zunger SIC approach.
Read the full abstract
Biological Physics · Cross-listed from q-bio.BM
Principles of Client Enrichment in Multicomponent Biomolecular Condensates
Biomolecular condensates are commonly organized by a small number of scaffold molecules that drive phase separation together with client molecules that do not condense on their own but become selectively recruited into the dense phase. A central open question is how client recruitment feeds back on scaffold interactions to determine condensate composition. Here we address this problem in a reconstituted focal adhesion system composed of focal adhesion kinase (FAK) and phosphorylated p130Cas (Cas) as scaffolds and the adaptor protein paxillin (PXN) as a client. We show that both FAK phosphorylation and PXN recruitment produce a common compositional response in which FAK becomes enriched while Cas is depleted within the condensate. To interpret these observations, we develop two complementary theoretical descriptions. First, within a two-component Flory-Huggins framework, we show that phosphorylation can be captured by either strengthening heterotypic FAK-Cas interactions or increasing the effective number of interaction-relevant segments on FAK, both of which bias partitioning toward FAK-rich condensates. Second, we introduce a minimal three-component Flory-Huggins theory without an explicit solvent and map it onto an effective two-component description, demonstrating that client recruitment renormalizes homotypic and heterotypic scaffold interactions. Analytical predictions for the location of the critical point are tested in reconstituted multicomponent systems through PXN addition, showing that client recruitment alone tunes proximity to criticality and reshapes condensate composition. Together, our results reveal distinct yet convergent physical routes by which post-translational modification and client recruitment control scaffold composition in multicomponent condensates.