Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Physics

physics Watch Physics (sign in to use)

Wednesday 29 April

Chemical Physics · Cross-listed from cond-mat.mtrl-sci

A Theoretical Investigation of the Thermal and Photochemical Mechanisms of Ethylbenzene Dehydrogenation on Rutile TiO$_{2}$(110)

This master's thesis investigates the thermal and photochemical dehydrogenation of ethylbenzene (EB) to styrene on the rutile TiO$_{2}$(110) surface. A dual-methodological quantum chemical approach is used for this investigation. While industrial styrene production is energy-intensive, photocatalysis on semiconductor materials offers a promising alternative under significantly milder conditions. To elucidate the underlying mechanisms, this study employs density functional theory (DFT-PBE-D3) for geometry optimization and high-level multi-reference methods (SA-CASSCF) to accurately describe the electronic complexity of excited states and radical intermediates. The investigation reveals that, on the stoichiometric surface, both thermal and photochemical pathways are dominated by proton-coupled electron transfer (PCET). The wavelength dependence observed in the literature is explained by how the system navigates electronic manifolds. 343 nm irradiation leads to rapid relaxation into the ground state, where high kinetic barriers persist. In contrast, 257 nm excitation enables the system to persist in higher excited states (S1/T2). This allows the reaction to bypass the rate-determining ground-state barrier. Furthermore, the study demonstrates that surface oxidation causes a fundamental mechanistic shift. On oxidized surfaces, pre-adsorbed oxygen radicals (O$_{Ti}$) enable direct hydrogen atom transfer (HAT), which is more efficient than PCET on reduced surfaces. This "hydrogen scavenger" effect explains the significant increase in styrene yield. This work underscores the necessity of multi-reference treatments for complex surface reactions and provides a fundamental understanding of how surface stoichiometry and photon energy govern photocatalytic efficiency.
Read the full abstract

Fluid Dynamics

Impulse-driven capillary detachment

Capillary interfaces subjected to impulsive forcing arise in many natural and technological systems, yet the pathway by which rapid substrate motion is converted into droplet detachment remains unclear. Here we study this process in a controlled setting: a liquid droplet resting on a taut wire that is plucked and suddenly released. The resulting transverse wave imparts a brief inertial forcing at the droplet base, initiating rapid stretching that precedes sheet formation and jet breakup. We show that the maximum extension prior to detachment is set by the mechanical work transmitted from the wire through capillary traction at the three-phase contact line, balanced by viscous dissipation during filament extension. This energetic balance identifies the contact line as the pathway by which mechanical impulse is converted into capillary deformation and governs impulsive droplet detachment.
Read the full abstract

All authors Dilip Kr. Maity, Sandip Dighe, Nilamani Sahoo, Tadd Truscott

Plasma Physics

Nitrogen-induced ELM suppression and confinement improvement in the EAST tokamak with a full metal wall

This paper reports the achievement of an ELM-free H-mode regime with confinement improvement enabled by nitrogen (N2) seeding on the Experimental Advanced Superconducting Tokamak (EAST) with a full metal wall. Following N2 injection, large Edge-Localized Mode (ELM) bursts are completely suppressed, while global energy confinement is significantly enhanced, with the H98 factor increasing from approximately 0.9 to 1.2. A distinct edge coherent mode (ECM), localized at the pedestal foot (psi_N ~ 0.99), is identified using O-mode Poloidal Correlation Reflectometry and AXUV diagnostics. This mode operates within a frequency range of 20-50 kHz with a poloidal wavenumber of k_theta ~ 0.54 cm^-1. Linear gyrokinetic simulations performed with the CGYRO code reveal a dominant instability that quantitatively matches the experimental measurements. Detailed scans of parameters identify this mode as a Dissipative Trapped Electron Mode (DTEM). The energy and particle transport driven by this pedestal-foot DTEM effectively regulates the edge gradients, preventing the pedestal from crossing the Peeling-Ballooning stability boundary and sustaining a stationary ELM-free state. These findings provide a physical basis for an integrated scenario to maintain high confinement and protect plasma-facing components in future steady-state fusion reactors.
Read the full abstract

All authors Jingyan Hu, Peng Shi, Chu Zhou, Jinyue Liu, Gongshun Li, Kangning Geng, Yiren Zhu, Kedong Li, Hailin Zhao, Xiang Jian, Ge Zhuang

Atomic Physics

Precision Spectroscopy of 2S-nS Transitions in Atomic Hydrogen: A Determination of the Proton Charge Radius

We present absolute frequency measurements of 2S_{1/2}-nS_{1/2} two-photon transitions with n = 8, 9, and 10 in a cryogenic beam of atomic hydrogen. Each transition has been measured with a fractional uncertainty of approximately 2.6*10^(-12). Combining the results from this work and the 1S_{1/2}-2S_{1/2} transition frequency, we extract a root-mean-square proton radius of r_p = 0.8433(31) fm and a Rydberg frequency of cR_{\infty} = 3,289,841,960,252.9(9.7) kHz. These are in good agreement with the CODATA 2022 recommended values.
Read the full abstract

All authors R. G. Bullis, W. L. Tavis, M. R. Weiss, J. Orellana Cisneros, A. J. Cheeseman, U. D. Jentschura, D. C. Yost

Instrumentation and Detectors

Comparison of Silvaco and Synopsys TCAD Predictions Including the Perugia Radiation Damage Model in Silicon Pixel Detectors for the HL-LHC

At the High Luminosity Large Hadron Collider (HL-LHC), silicon pixel detectors will be exposed to radiation fluences about 5 to 10 times larger than those experienced by the current innermost pixel layers up to today. Due to radiation damage to bulk of pixel detectors, leakage current and depletion voltage will increase significantly over time, posing severe constraints on operating conditions, with important modifications to the electric field profile due to radiation induced deep defects in silicon bulk. It is important to have reliable predictions for all observables - such as leakage current level and breakdown voltage - after irradiation, in order to estimate operational voltage values. In this paper, the predictions of Silvaco and Synopsys TCAD device simulations are compared when the surface and bulk defects and traps of the “New University of Perugia radiation damage model” are included by studying observables like leakage current, depletion and breakdown voltage, along with electric field and trap occupancy profiles. The results are quite promising regarding leakage current, depletion voltage, electric field and trap statistics, at two distinct reference temperatures and fluences.
Read the full abstract

All authors M. Bomben, T. Croci, K. Aouadj, A. Fondacci, F. Moscatelli, A. Morozzi, D. Passeri

Computational Physics · Cross-listed from math.NA

Drift-Free Conservative Dynamics from Quantized Interaction Rules

Conservation laws are conventionally discretized through floating-point flux evaluation, with invariants obtained by cancellation of approximate interface contributions and admissible weak solutions selected by reconstruction and Riemann solvers. Here we introduce an operator-level formulation in which conservative dynamics is realized as an exact discrete interaction rule on a quantized state space. The update is defined by an antisymmetric integer-transfer operator, which enforces conservation exactly at the arithmetic level and eliminates round-off drift from the primitive evolution. For scalar laws, monotone order-preserving transfers select admissible shock structures within the primitive update, rather than through flux reconstruction. Numerical experiments show that the interaction rule preserves high-frequency transport near the Nyquist limit and maintains sharply localized discontinuities in Burgers dynamics. The same construction extends to multidimensional problems and systems of conservation laws through oriented, vector-valued integer transfers. The results show that exact integer-transfer dynamics can suppress cumulative transport drift while preserving entropy-shock localization in nonlinear conservative evolution.
Read the full abstract

Biological Physics · Cross-listed from cond-mat.soft

Coexistence of patterned phases in chemically active multicomponent mixtures

Chemically active mixtures exhibit complex patterns that emerge from the interplay of physical interactions and reactions among components. Individually, these two processes are well-understood: Physical interactions can give rise to phase separation, whereas reactions can form reaction-diffusion patterns. To understand the combination of both processes, we identify a Lyapunov functional for a class of chemical reactions. By minimizing this functional, we identify a generalized Gibbs phase rule that governs the number of coexisting patterns, and we demonstrate that complex patterns can be created by the modular combination of independent phases. Our theory unveils complex stationary patterns in chemically active mixtures and provides a framework for analyzing more complex systems.
Read the full abstract

All authors Chengjie Luo, Yicheng Qiang, Guido L. A. Kusters, David Zwicker

Physics and Society

People, Places & Things: Network topology & motifs of R&D missions

Challenge-led R and D programs increasingly assemble heterogeneous people, organizations, funders, projects, and technical outputs around defined missions. Yet program evaluation often describes these systems through project lists, output counts, or retrospective case narratives. This article develops a typed network framework for representing R and D program architecture directly. We model programs as networks of people, places, and things: researchers, program directors, institutions, funders, publications, patents, projects, and citations. Applied to ARPA-E project impact sheets from the agency's first decade, the framework reconstructs 23 program-induced networks and an agency-level composed network. We show that R and D programs have an analysable topology: a typed arrangement of people, institutions, funders, projects, publications, patents, and citations that can be reconstructed, compared, and monitored. The analysis shows that programs can be compared by their local structural patterns, that cross-program overlap is concentrated more in recurring institutions than in individual researchers, and that program fingerprints differ across thematic areas. The article contributes to network science by extending topological analysis to R and D program systems, a class of governed, typed, and output-generating networks that has not been systematically represented in existing innovation-network work.
Read the full abstract

All authors Henry C. W. Price, Martin Ho, Tim S. Evans, Eoin O'Sullivan

Chemical Physics · Cross-listed from quant-ph

Amplitude Encoding of Slater-Type Orbitals via Matrix Product States: Efficient State Preparation and Integral Evaluation on Quantum Hardware

Slater-type orbitals (STOs) provide the physically correct description of atomic wavefunctions but have been largely replaced by Gaussian-type orbitals in computational chemistry due to the lack of closed-form multi-center integrals. We present a systematic study of amplitude encoding of STOs on quantum computers using matrix product states (MPS). For one-dimensional orbital functions of the form $p_d(x) e^{-\zeta x}$, we derive analytical MPS constructions with constant bond dimension $\chi = d + 1$, requiring $O(n)$ classical and quantum resources for $n$-qubit registers with no grid sampling. We demonstrate a complete one-electron integral pipeline – overlap, kinetic energy, and nuclear attraction – in one dimension, validating the overlap and kinetic energy on IBM Heron processors at 5~qubits with 0.67% hardware-induced error using Zero-Noise Extrapolation. In three dimensions, we compute multi-center overlap integrals between 1s and 2s orbitals in Cartesian coordinates with 0.02% discretization error at 18~qubits. A systematic entanglement analysis reveals that the MPS bond dimension of three-dimensional STOs in Cartesian coordinates saturates with increasing grid resolution – reaching $\sim$138 for the hydrogen 1s orbital at 12~qubits per coordinate – establishing bounded encoding complexity rather than the exponential scaling initially expected. The SVD truncation threshold provides a practical resource parameter, reducing the bond dimension to 39 at threshold $10^{-6}$ with negligible accuracy loss. These results map the entanglement landscape for amplitude encoding of atomic orbitals and establish MPS-based state preparation as a viable path toward exact STO basis sets on quantum computers.
Read the full abstract

Atmospheric and Oceanic Physics · Cross-listed from math.AP

A mathematical study of an elastic-viscous-plastic sea-ice model with the Kelvin-Voigt rheology

Motivated by the elastic-viscous-plastic (EVP) sea-ice model [E. C. Hunke and J. K. Dukowicz, J. Phys. Oceanogr., 27, 9 (1997), 1849–1867], which is used in large-scale numerical climate simulations, we proposed in [D. W. Boutros, X. Liu, M. Thomas and E. S. Titi, arXiv:2505.03080 (2025)] the use of the inviscid Voigt regularisation for the constitutive (stress-tensor) relation and proved the global well-posedness of the resulting model. The EVP model treats sea ice as a non-Newtonian fluid. In turn, elastic-viscous-plastic solids often involve a Kelvin-Voigt viscosity in terms of the strain rate. Therefore, in the present work we formulate an elastic-viscous-plastic sea-ice model with a Kelvin-Voigt regularisation in terms of the strain rate. In other words, we introduce the Voigt regularisation in the momentum balance rather than in the constitutive relation (for the stress tensor). We then prove the local well-posedness for the Kelvin-Voigt EVP model with the advection term, in the momentum balance, and the global well-posedness in the absence of the advection term (following a very standard approximation in the latter case). A crucial component of the proof of these results, is a new $L^\infty$-estimate for the stress tensor which relies on the damping structure. Note that, both with and without the advection term, we are able to handle the case of viscosity coefficients without a cutoff from above, which remains a major open problem for the closely related Hibler sea-ice model. We are also able to prove the existence of solutions for much less regular initial data compared to our previous paper on the Voigt-EVP model.
Read the full abstract

All authors Daniel W. Boutros, Xin Liu, Marita Thomas, Edriss S. Titi

Fluid Dynamics

Coherent structures in Newtonian and viscoelastic turbulent planar jets

The addition of a small amount of long-chain polymers confers viscoelastic properties to Newtonian flows. The resulting non-Newtonian solution now exhibits different dynamics, such as enhanced mixing at low Reynolds, where elastic instabilities can trigger elastic turbulence even though inertial turbulence is absent. Here, we study this phenomenon in viscoelastic planar jets and, in particular, we do it from the perspective of coherent structures to understand how elastic turbulence is triggered and sustained, which remain barely explored in this setup. We introduce the spatio-temporal Koopman decomposition for extracting the dominant flow patterns, and we compare them with those from Newtonian planar jets at high Reynolds number. Global flow structures are similar between jets, with low-frequency streaks and high-frequency wave packets dominating the turbulent dynamics. However, structures are strikingly different in the near field, where elasticity-driven streaks affect the dynamics in the potential core of the viscoelastic planar jet, modifying the bulk flow and interacting with the flow instability. The analysis of the polymer field reveals stretched polymer filaments and centre-mode structures, which support the implication of the near-field streaks on sustaining elastic turbulence in three-dimensional viscoelastic planar jets.
Read the full abstract

All authors Christian Amor, Adrián Corrochano, Giovanni Soligo, Soledad Le Clainche, Marco Edoardo Rosti

Fluid Dynamics

Scaling in Supersonic Turbulence: Energy Spectra and Fluxes using High-Fidelity Direct Numerical Simulations

Supersonic turbulence is vital to astrophysical and high-speed engineering flows, yet its energy transfer mechanisms remain poorly understood. We present high-resolution ($1024^3$) direct numerical simulations (DNS) of forced compressible turbulence across a range of turbulent Mach numbers ($M_t = 0.2$ to $3.0$). Using the GPU-accelerated solver DHARA with a seventh-order, low-dissipation Targeted Essentially Non-Oscillatory (TENO) scheme, we resolve both fine-scale eddies and sharp shock fronts. Our results reveal a fundamental shift in the energy cascade in the supersonic regime. As $M_t$ increases, the rotational kinetic energy spectrum steepens from a Kolmogorov-like $k^{-5/3}$ scaling toward a Burgers-like $k^{-2}$ scaling. Conversely, the compressive energy spectrum becomes shallower, deviating from Burgers scaling. We show that these spectral modifications are driven by a dominant cross-scale transfer of energy from solenoidal to compressive modes within the inertial range, alongside significant contributions from pressure dilatation. Scaling laws for the root-mean-square compressive velocity ($U_C$) and compressive energy flux ($\Pi_C$) are found to mirror classical Burgers turbulence. Finally, we show that while energy injection rates depend on forcing type rather than Mach number, increased $M_t$ leads to decreased rotational dissipation and increased compressive dissipation and pressure dilatation. These findings elucidate intermodal energy cascade mechanisms, advancing our understanding of energy transfers in supersonic turbulence.
Read the full abstract

All authors Harshit Tiwari, Dhananjay Singh, Mahendra K. Verma, Rajesh Ranjan

Physics and Society

Digital Epidemiology with Awareness-Based Event-Triggered Migration in Networked Cyber-Physical Systems

Understanding how human mobility and information propagation influence the course of an epidemic remains a key challenge in digital epidemiology. In this work, we develop a new awareness-based, event-triggered epidemic model embedded within a networked Cyber-Physical System (CPS). In our framework, disease transmission and the dissemination of epidemic-related information evolve together on two interconnected layers. In detail, the physical layer models disease spread through human movement between two types of locations - residences and transfer stations - forming a bipartite metapopulation network. This structure captures the rendezvous effect, which reflects how gatherings in shared locations contribute to infection spread. The cyber layer represents the flow of information through digital communication networks. We introduce an event-triggered migration regulation mechanism, whereby individuals adapt their movement patterns based on local awareness thresholds, leading to a decentralized control process embedded within the network. Using a microscopic Markov chain approach (MMCA), we derive the epidemic threshold analytically and validate our results through extensive Monte Carlo simulations. Our findings show that event-triggered migration effectively suppresses the overall spread of the disease and lowers infection peaks - especially in heterogeneous populations and densely connected gathering points. These results demonstrate the potential of CPS-based epidemic models to enable real-time, awareness-driven interventions and to inform the design of decentralized control strategies that leverage digital communication dynamics.
Read the full abstract

All authors Yusheng Li, Minyu Feng, Liang-jian Deng, Matjaž Perc, Jürgen Kurths

Fluid Dynamics

Reduced-order modeling of a viscoelastic turbulent jet with hybrid machine learning models

Adding flexible polymers to a Newtonian solvent confers complex properties to the resulting solution. The additional complexity substantially increases the computational cost of numerical simulations, which often makes them prohibitively expensive. Here, we propose hybrid reduced-order models to accelerate simulations of viscoelastic turbulent jets. The model combines modal decompositions with deep networks: we use proper orthogonal decomposition to obtain a compact representation of the data, and a neural network is trained to predict the mode coefficients in the low-dimensional space. Results show that the hybrid model effectively captures the long-term behavior of the viscoelastic jet, that we demonstrate by computing relevant statistics of the jet. While small models are capable of predicting large-scale dynamics more than one-step at a time, thus facilitating greater accelerations, larger models are mandatory for forecasting smaller-scale dynamics, with skip connections the most effective strategy for deeper and generalizable models. The proposed methodology underpins the potential of hybrid approaches for compact and robust reduced-order models of viscoelastic turbulent jets.
Read the full abstract

All authors Christian Amor, Adrián Corrochano, Marco Edoardo Rosti, Soledad Le Clainche

Fluid Dynamics

Revisiting the mixing length scaling in pressure-gradient turbulent boundary layers via a symmetry approach

A century after Prandtl's mixing length hypothesis, full-profile scaling of the mixing length in pressure-gradient turbulent boundary layers (PG TBLs) remains debated, especially for adverse pressure gradients (APGs). This work presents a symmetry-based analytical model for the mixing length in equilibrium APG TBLs by extending the structural ensemble dynamics theory and coupling a two-layer total shear stress model. The framework unifies the inner layer, logarithmic region, half-power-law transition zone, and wake region with an invariant Karman constant. A critical Clauser parameter is identified, above which the logarithmic layer shrinks and transitions to the half-power-law scaling. The wake-region mixing-length parameter is analytically formulated, and the viscous sublayer and buffer layer thicknesses are determined self-consistently without ad hoc fitting. With only one finite-Reynolds-number correction parameter determined by the maximum shear stress, the model accurately predicts full profiles of mixing length, mean velocity, and Reynolds shear stress, validated against extensive numerical and experimental data. This work provides a unified, physically consistent framework for mixing-length scaling in PG TBLs and clarifies the transition mechanism from the log law to the half-power law under strong APG. It also enables assessment of the invariance of the logarithmic law and Karman constant using the full-profile scaling law of the mixing length.
Read the full abstract

Chemical Physics

Seniority-zero Quadratic Canonical Transformation Theory

We propose a method to solve the Schrödinger equation for systems with static/strong electron correlation using Hamiltonian transformations. Building on our previous work on seniority-zero canonical transformation theory, which seeks a unitary transformation that maps the Hamiltonian into the seniority-zero space, this method presents an alternative way of evaluating the Baker–Campbell–Hausdorff (BCH) expansion based on quadratic canonical transformation theory. The extension aims to relax the small-generator constraint by allowing approximate four-body contributions in the expansion, thus expanding the class of excitations previously allowed in SZ-LCT, where only approximate three-body operators were retained. Numerical tests reveal that the seniority-zero quadratic canonical transformation method (SZ-QCT) delivers good accuracy, with most errors within chemical accuracy. In particular, SZ-QCT shows sub-millihartree errors in cases where larger generators are needed to recover the residual dynamic correlation. The computational scaling of SZ-QCT is the same as that of SZ-LCT, $\mathcal{O}(N^8/n_c)$, where $n_c$ is the number of cores available for the computation
Read the full abstract

Fluid Dynamics · Cross-listed from cond-mat.soft

Viscous Settling of Bravais Unit-Cells

We study experimentally and theoretically the Stokesian settling of a well-known class of porous shapes: Bravais lattice unit-cells, whose porosity we vary controllably by changing their lattice spacing. In our experiments, conducted in a square cuboidal container with its long-axis aligned along gravity, we find that the settling speed U and the solid fraction {\phi} of these lattice units obey a power-law relationship U $\propto$ {\phi}^{\gamma} , with an exponent {\gamma} = 0.43 independent of their shape. To understand the observed scaling exponent, we analytically and numerically investigate the settling of the simple cubic structure under different approximations. We find that the walls of the container, though far from the sinking object, have a defining effect. Our Stokesian boundary integral simulations show that the Faxen's boundary correction captures the wall-effects accurately and enables us to discount the wall-effect from the experimental data, yielding a power-law exponent {\gamma} = 0.30 for settling in an unbounded domain. The power-law relating sinking speed and porosity is a step towards predictively understanding the sedimentation fluxes of complex objects in the clouds and the oceans. However, the applicability of this universal scaling to irregular and biologically richer aggregates found in nature remains an open direction.
Read the full abstract

All authors Sebastian Bürger, Harshit Joshi, S Ganga Prasath, Rahul Chajwa, Rama Govindarajan

Tuesday 28 April

Physics and Society

Achieving Generational Peace in Mali through Intergenerational Mean-Field-Type Game-based Incentives

This article develops an intergenerational mean-field-type game (MFTG) to model Mali's and neighbouring countries multi-actor conflict ecosystem, which includes formal state forces, traditional hunters, nonstate militias, jihadists, criminal networks, civil societies, and international proxies. Each decision-maker (agent, a group of agents or representative agent) is defined by a type, state, information structure, and action, with payoffs dependent not only on individual decisions but also on the evolving distribution of all agents' profiles. The model reveals that cycles of violence can persist across multiple generations due to the embedded presence of retaliatory types such as revenger child-soldiers whose trauma-conditioned best-responses favor conflict, and whose behavior reinforces intergenerational transmission of violence. The model also captures the strategic exploitation of institutional fragility by war entrepreneurs who profit from sustained instability through arms sales, militia contracting, and unregistered market mediation. These actors inject minimal resources to trigger profitable escalations, turning latent tensions into self-reinforcing violence economies. We show that in the absence of structural counterincentives, peaceful strategies are non-absorbing, and violence remains dynamically rewarding for war entrepreneurs. However, by embedding incentive-compatible, information-adaptive transfers directly into instantaneous payoffs, rewarding verifiable peacebuilding and penalizing aggression, it is possible to shift the mean-field-type equilibrium distribution intergenerationally toward more peaceful types and drive systemic de-escalation. We also discuss about the funding and the real implementation of such mechanisms in the field.
Read the full abstract

Applied Physics

A micromechanical frequency reference with parts-per-trillion holdover stability

Microelectromechanical (MEMS) resonators are widely used in timekeeping applications, and recent advances in fabrication, materials, and encapsulation technology have advanced their potential as high stability frequency references. However, for holdover applications that require the highest levels of long-term frequency stability, compact vapor atomic clocks remain dominant. In this work, we demonstrate a 268 MHz MEMS clock that achieves record fractional frequency stability of ~8 parts-per-trillion at an averaging time of 8 hours, competitive with chip-scale atomic clocks. We achieved this using a single-crystal silicon electrostatic resonator that has no currently known intrinsic drift mechanism and is protected from the environment with a wafer-level encapsulation. We specifically identify gain variations in the sustaining electronics as the dominant limitation in conventional phase-locked oscillator architectures – originating from temperature sensitivity and drifts in the electronic components – and overcome this by implementing a frequency-locked loop architecture based on dual-frequency resonance tracking (DFRT). This novel approach removes the specific gain of the supporting electronics as a frequency determining variable in the oscillator. When combined with dual-mode tracking and ratiometric temperature stabilization of the resonator, this approach enables a dramatic enhancement to long-term frequency stability and establishes gain-insensitive DFRT locking as a general paradigm for high-stability MEMS clocks.
Read the full abstract

All authors Jie Yan, Jintark Kim, Rakibul Islam, Jiawei Yang, Karim Elmeligy, Alkim Bozkurt, Thomas W. Kenny, Pavan K. Hanumolu, Gaurav Bahl

Accelerator Physics

Proposal for applying the novel gas-dynamic ion-beam extraction and bunching technique to the cryogenic stopping cells at FAIR

To enhance the quality of pulsed cold ion beams extracted from the two cryogenic stopping cells at FAIR (i.e., the one currently used in the FRS and the one under development for the Super-FRS at FAIR), we propose using a novel gas-dynamic ion beam extraction and bunching technique as an alternative to the radiofrequency quadrupole (RFQ) method. This technique allows for 100% ion transmission by placing a short stack of thin cylindrical electrodes behind the extraction RF carpet. Detailed gas-dynamic and ion trajectory computer simulations demonstrate that implementing this proposal will enable the achievement of world-record emittance values for ion beams in a wide mass range. The results of these simulations are presented and discussed.
Read the full abstract

General Physics

Pairwise-comparison-valued cosurfaces: a projective framework for multi-scale relational structures

We introduce cosurfaces with values in the group \(\PC_n(H)\) of \(H\)-valued reciprocal pairwise comparison matrices. The composition law is covariant on upper triangular coefficients and contravariant on lower triangular coefficients, which makes \(\PC_n(H)\) a natural target for oriented gluing constructions. Starting from a directed family of finite oriented discretizations, we define finite configuration spaces, coarse-graining maps induced by ordered refinements, and the associated universal projective limit. This yields a multi-scale organization of local comparative data in which global objects are reconstructed only through compatibility across scales. In the stochastic setting, projectively compatible probability laws define a cylindrical semantics on the limit space. We also introduce inconsistency observables, interpreted as discrete curvature-type defects measuring obstructions to global coherence. The resulting framework is simultaneously geometric, algebraic, and probabilistic, and suggests a foundational perspective on relational structures built from local comparisons rather than absolute observables.
Read the full abstract

Chemical Physics

Thermodynamic Properties of Diatomic Molecules from the Frost-Musulin Potential

In this study, we present a quantum-statistical analysis of H$_2$ and LiH diatomic molecules within the Frost–Musulin potential framework. By combining the analytical bound-state approach to the radial Schrödinger problem with the near-equilibrium Pekeris representation, we obtain a validated rotation-vibration spectrum that reproduces a physically consistent ordering of energy levels. These bound states are subsequently combined with standard translational and rotational ideal gas contributions to construct the total partition function and the corresponding thermodynamic observables of the ground state. The resulting formulation captures the Gibbs free energy deviation function for both molecules with high quantitative accuracy and provides chemically plausible trends for heat capacity and enthalpy increase over a wide temperature range. At the same time, residual errors become increasingly pronounced in derivative-sensitive quantities, particularly at high temperatures; this indicates that the dominant limitations now stem not from the local bound-state spectrum itself, but from the neglect of inelastic rotational, continuity contributions and dynamics close to dissociation. Consequently, the present results define the potential model as a compact and analytically tractable representation of the bound region, recovering a significant portion of the observed thermochemistry whilst also delineating the regime where more comprehensive molecular statistical mechanics is required.
Read the full abstract

All authors Mohammadjavad Parsanasab, Reza Khordad, Meysam Asadipour, Ahmad Ghanbari, Vatan Badalov

Chemical Physics

Novel Chemical Pathways for the Formation of Nucleobase Precursors via Benzene {\pi}-Bond Addition to HCN

We propose a simple and efficient pathway for the formation of precursors to core nucleobases in DNA and RNA using a suite of computational chemistry methods. Benzene, which is thermochemically stable in N2- or CO2-dominated atmospheres, could have formed via upper-atmospheric photochemistry or surface lightning and accumulated on the early Earth or Mars. However, nitrogen insertion into the benzene ring to form pyrimidine and purine is widely considered to be challenging. We propose that nitrogen incorporation occurred through HCN 1,4-cycloaddition to benzene's {\pi}-system, followed by a C2H2 fragmentation mechanism, as confirmed by quantum chemistry calculations. This pathway, potentially facilitated by photochemistry at the ocean surface or episodic impact events on local reservoirs, can lead to pyrimidine formation, which can further react with NH3 and HCN to produce purine. Extending this pathway to early Mars, our photochemical model simulates heterocyclic compound formation under cold, dry surface conditions that favor high benzene and HCN concentrations but lack liquid water. We thus propose that organics formed during dry phases may have later dissolved into surface waters during wet phases and become concentrated as ocean sediments. This result supports Mars Sample Return efforts focused on ancient aqueous environments likely to retain prebiotic signatures.
Read the full abstract

All authors Jeehyun Yang, Danica J. Adams, Renyu Hu, Yuk L. Yung

Applied Physics · Cross-listed from quant-ph

Naturally Resonant Emitters: Approaching Fundamental Antenna Limits

Antenna miniaturization remains a critical technological challenge across frequency scales - from microwave RF links in phones and wearables to VLF for underwater-to-air communications and ionospheric probing. At deeply subwavelength scales conventional antennas require complex and lossy matching circuits due to absent intrinsic material resonances, motivating resonant electrically small emitters (ESEs) like mechanical resonators and quantum emitters. Here, we extend the theory of electrically small antennas (ESAs) to this broader ESE class, deriving the fundamental efficiency limit for a unit volume emitter at given frequency and bandwidth. Our figure of merit (FOM) - quantifying proximity to this limit - enables direct comparison across ESE types, frequencies, bandwidths and scales. We demonstrate its utility using public data from ELF and VLF Navy facilities alongside two mechanical ESEs reported in literature. The measurements reveal that mechanical antennas operate near theoretical FOM limit, questioning claims of possible further orders-of-magnitude gains. A naturally resonant emitter is still subject to the Chu-Harrington limit (CHL) under its standard assumptions. Indeed, we derive novel CHL-dictated constraints on atomic ESE properties: lower bound on excited-state lifetime and an upper bound on transition dipole moment.
Read the full abstract

Physics and Society

Power, Depletion and Energy Quality Model of Thermo-industrial Civilization

The current thermo-industrial civilization is critically dependent on fossil fuel energy sources. An intuitive model capturing the interplay between economic activity, physical power consumption, depletion and energy quality is presented.

Optics

Taming Randomness in Random Lasers: Programmable Disorder for Active Control of Random Lasing via Electric-Field-Directed Assembly of Nanowires

Random lasing exploits multiple scattering to provide optical feedback without conventional resonant cavities, enabling simplified architectures that are readily integrated into compact photonic platforms such as wearable sensors and lab-on-chip devices. However, the same disorder that enables cavity-free lasing also makes it challenging to control and tune the emission properties. Here, an electrically reconfigurable random-lasing platform based on dielectrophoretic assembly of chaining silver nanowires suspended in a dye gain medium is reported. An applied electric field across patterned quadrupole electrodes induces nanowire chaining and programmable alignment, enabling real-time reconfiguration of the disorder landscape. Based on the electrically driven disorder state transitions, tunable random-lasing characteristics, including reduced lasing thresholds, modulation of emission intensity, and control of the polarization state have been demonstrated. Simulations further indicate that chaining enhances scattering relative to absorption, providing more efficient radiative feedback, and the orientation of the nanowire network governs the polarization dependence of the system. These results establish a route to actively modulate random lasing through controllable disorder and point toward adaptive, reconfigurable photonic light sources and sensing systems.
Read the full abstract

All authors Jinkai Yang, Kumudu N. Ranasinghe, Lei Kang, Jennifer R. Decker, Cheng-Yu Wang, Douglas H. Werner, Christine D. Keating, Zhiwen Liu

Computational Physics

Mixture of Experts Framework in Machine Learning Interatomic Potentials for Atomistic Simulations

First-principles atomistic simulations are essential for understanding complex material phenomena but are fundamentally limited by their computational cost. While Machine Learning Interatomic Potentials (MLIPs) have drastically improved cost for a given accuracy, their inference cost remains a bottleneck for massive systems or long timescales. To address this, we introduce a multifidelity "Mixture-of-Experts" framework based on the E(3)-equivariant Allegro architecture. Our method spatially partitions the simulation domain into a chemically complex region (e.g., reactive interfaces) and a simple region (e.g., bulk lattice), assigning models of varying capacity to each. Among the challenges in such static domain decomposition, the mechanical mismatch between models at the interface is particularly critical, as it can generate artificial stress fields and instability. We address this challenge with a co-training strategy in which the loss function includes agreement constraints – penalties on per-atom energy and force discrepancies between models evaluated on shared bulk environments – forcing the independent models to learn a consistent physical description of the bulk material. We validate this approach on a realistic Pt+CO catalytic system, demonstrating that the co-trained models maintain exact energy conservation, align their bulk mechanical response (e.g., equation of state and bulk modulus), and achieve predictive accuracy comparable to a full high-fidelity simulation at more than twice the computational speed.
Read the full abstract

All authors Gabriel de Miranda Nascimento, Marc L. Descoteaux, Laura Zichi, Chuin Wei Tan, William C. Witt, Nicola Molinari, Sriteja Mantha, Daniil Kitchaev, Mordechai Kornbluth, Karim Gadelrab, Charles Tuffile, Boris Kozinsky

Chemical Physics

Solv-eze: Automated Placement of Explicit Water Molecules Using 3D-RISM

Molecular dynamics (MD) simulations are widely used to study biological systems, where water molecules often play a critical role in protein-ligand interactions. In conventional MD preparation protocols, water molecules are typically added from a pre-equilibrated solvent box and removed using conservative steric cutoffs, an approach that can eliminate important interfacial waters that are often not recovered during equilibration due to kinetic barriers limiting exchange with bulk solvent. In this work, we present an automated and computationally efficient method for placing water molecules around biomolecular solutes using three-dimensional reference interaction site model (3D-RISM) solvent density distributions. By identifying regions of high solvent probability, the method generates physically meaningful initial hydration structures without requiring extended sampling or specialized techniques such as grand canonical Monte Carlo (MC) or hybrid MC/MD approaches, and will be released as an update to AmberTools 26, enabling seamless integration into standard MD preparation pipelines. We validated the approach on a diverse set of protein-ligand complexes with crystallographically resolved bridging waters, showing that the method reproduced over 80% of experimentally observed bridging waters and 85% of buried waters not accessible to the bulk. Subsequent energy minimization of both crystallographic and predicted waters further improved agreement. Overall, this method enables more accurate and practical initialization of interfacial hydration, improving the reliability of MD simulations with modest computational cost relative to routine system preparation.
Read the full abstract

All authors Felipe Silva Carvalho, Steven Ramsey, George M. Giambaşu, Tom Kurtzman, Tyler Luchko

Space Physics · Cross-listed from astro-ph.SR

Parker Solar Probe Observations of Compound Reconnection Exhaust Boundaries and Mirror-Mode Structures in the Near-Sun Heliospheric Current Sheet

Magnetic reconnection is a fundamental physical process that can drive rapid conversion of magnetic energy into plasma bulk flows, thermal heating, and particle acceleration in space and astrophysical plasmas. Classical reconnection theory predicts that the Alfvenic reconnection exhausts are bounded by pairs of slow-mode shocks. However, identifying and characterizing these shocks through in situ spacecraft observations remains a challenge. Here we report Parker Solar Probe (PSP) observations of a reconnection exhaust embedded in the heliospheric current sheet (HCS) at a heliocentric distance of 12.2 R_O. The reconnection exhaust is bounded on both boundaries by compound magnetic structures rather than a pair of pure slow shocks. Each boundary consists of a rapidly evolving, steep inner slow shock, whose Mach numbers and shock-normal angles change significantly within several minutes, and an outer, gradual compound structure which comprises a slow shock and a rotational discontinuity. These slow shocks are quasi-perpendicular and are accompanied by enhanced proton perpendicular heating. Deep within the reconnection exhaust, high perpendicular temperature together with large plasma beta trigger mirror instability and generate mirror-mode structures. These observations provide new insights into the structure of reconnection exhaust boundaries and their role in energy conversion in the near-Sun plasma.
Read the full abstract

All authors Weijie Sun, Tai Phan, Jia Huang, Yi-Hsin Liu, James A. Slavin, Orlando Romeo, Mingzhe Liu, Vassilis Angelopoulos, Ali Rahmati, Davin Larson, Nehpreet Walia, Stuart Bale, Marc Pulupa, Jiutong Zhao, Roberto Livi

Fluid Dynamics · Cross-listed from math.NA

An Asymptotic-Preserving Dual Formulation Finite-Volume Method for the Thermal Rotating Shallow Water Equations

We propose a new second-order asymptotic-preserving (AP) dual formulation finite-volume (DF-FV) method for the thermal rotating shallow water (TRSW) equations. The TRSW system models geophysical flows characterized by horizontal temperature/density variations, exhibiting multi-scale dynamics due to the coexistence of fast rotational waves and slower advective processes. To efficiently address challenges associated with the multiscale nature of the TRSW system, we follow the DF-FV framework and develop a DF-FV method, in which both the conservative and nonconservative (primitive) forms of the equations are simultaneously solved, allowing the method to exploit the complementary strengths of each representation across different flow regimes. The primitive formulation is better suited for preserving the correct asymptotic behavior in nearly thermal quasi-geostrophic (TQG) regimes characterized by a low Rossby number, while the conservative formulation is essential for robust shock capturing in high-Rossby-number regimes, in which nonconservative discretizations may fail to converge to physically relevant weak solutions.
Read the full abstract

All authors Alina Chertock, Alexander Kurganov, Lorenzo Micalizzi, Nan Zhang