Peer Revenue Play credits. Real papers.

Papers

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

Physics Physics

physics Watch Physics (sign in to use)

Thursday 6 August

Fluid Dynamics

Nonisothermal Shock Structure and Universal Flow-Index Thresholds in a Hyperbolic Power-Law Fluid

We investigate whether the two flow-index thresholds previously found for isothermal shock profiles persist when the full nonisothermal dynamics is taken into account in a hyperbolic power-law relaxation model of Rational Extended Thermodynamics. The nonisothermal profile problem is structurally different from its isothermal counterpart: restoring the energy balance determines the temperature along the traveling wave and feeds it back into the pressure, the relaxation production, the temperature-dependent consistency coefficient, and the characteristic structure. Under thermodynamic stability, $p_\theta\ge0$, and strict convexity of the reduced Hugoniot pressure, no nontrivial constant-temperature compressive profile can satisfy the full equations. We derive an exact global characteristic-ordering identity and prove that the positive nonequilibrium characteristic speed has its strict global minimum at the unperturbed upstream state. Consequently, a monotone continuous profile exists for $1<M_0<\Mst$, whereas for $M_0>\Mst$ the Boillat–Ruggeri theorem excludes a $C^1$ profile and any admissible piecewise-smooth connection must contain a subshock. Despite the thermomechanical coupling, the shock-thickness classification remains unchanged: $m=2$ is the weak-shock threshold and $m=1$ the near-critical threshold as $M_0\nearrow\Mst$. The corresponding exponents are constitutive-independent within the present class, while finite limiting values and prefactors depend on the equation of state, internal energy, and temperature-dependent consistency coefficient. For the Tait–Murnaghan example, increasing the reference temperature lowers the critical Mach number when the dimensional viscous–relaxation scale is fixed and, for the thermally thinning law considered, reduces the resolved shock thickness.
Read the full abstract

Geophysics · Cross-listed from hep-ex

Estimating the sensitivity of the IceCube Upgrade to probe the interior of the Earth using atmospheric neutrino oscillations

The IceCube Upgrade is a densely instrumented central region of the IceCube Neutrino Observatory, deployed during the 2025-26 polar season. It will reduce the detector's energy threshold and improve overall reconstruction capabilities for multi-GeV atmospheric neutrinos, which in turn enhance their sensitivity to Earth matter effects as they traverse through the deep Earth. In this study, we describe the potential of the IceCube Upgrade to observe Earth matter effects on atmospheric neutrinos and estimate the detector's sensitivity to probe key features of the Preliminary Reference Earth Model by utilizing these observations. We highlight the IceCube Upgrade's capability to estimate the mass of the Earth and verify the non-homogeneous distribution of matter density within the Earth. We also estimate the IceCube Upgrade sensitivity to measure the correlated densities of the Earth layers while incorporating constraints from the mass and moment of inertia of the Earth. Neutrino-based results would be independent and complementary to the seismic and gravitational measurements.
Read the full abstract

All authors The IceCube Collaboration: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J.M. Alameddine, S. Ali, N. M. Amin, K. Andeen, C. Argüelles, S. Athanasiadou, S. N. Axani, R. Babu, X. Bai, A. Balagopal V., S. W. Barwick, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, P. Behrens, J. Beise, C. Bellenghi, S. Benkel, S. BenZvi, D. Berley, E. Bernardini, D. Z. Besson, E. Blaufuss, L. Bloom, S. Blot, F. Bontempo, J. Y. Book Motzkin, C. Boscolo Meneguolo, S. Böser, O. Botner, J. Böttcher, J. Braun, B. Brinson, Z. Brisson-Tsavoussis, L. Brusa, R. T. Burley, D. Butterfield, K. Carloni, J. Carpio, S. Chattopadhyay, N. Chau, Y. C. Chen, Z. Chen, D. Chirkin, S. Choi, A. Chubarov, B. A. Clark, D. A. Coloma Borja, A. Connolly, J. M. Conrad, D. F. Cowen, C. De Clercq, J. J. DeLaunay, D. Delgado, T. Delmeulle, S. Deng, P. Desiati, K. D. de Vries, G. de Wasseige, T. DeYoung, J. C. Díaz-Vélez, S. DiKerby, T. Ding, M. Dittmer, A. Domi, L. Draper, L. Dueser, D. Durnford, K. Dutta, M. A. DuVernois, T. Ehrhardt, L. Eidenschink, A. Eimer, C. Eldridge, P. Eller, E. Ellinger, D. Elsässer, R. Engel, H. Erpenbeck, W. Esmail, S. Eulig, J. Evans, P. A. Evenson, K. L. Fan, K. Fang, K. Farrag, A. Fattorini, A. R. Fazely, A. Fedynitch, N. Feigl, C. Finley, D. Fox, A. Franckowiak, S. Fukami, P. Fürst, J. Gallagher, E. Ganster, A. Garcia, M. Garcia, E. Genton, L. Gerhardt, A. Ghadimi, C. Glaser, T. Glüsenkamp, J. G. Gonzalez, S. Goswami, A. Granados, D. Grant, S. J. Gray, S. Griffin, S. Griswold, K. M. Groth, D. Guevel, C. Günther, P. Gutjahr, C. Ha, A. Hallgren, L. Halve, F. Halzen, L. Hamacher, M. Handt, K. Hanson, J. Hardin, A. A. Harnisch, P. Hatch, A. Haungs, J. Häußler, K. Helbing, J. Hellrung, B. Henke, L. Hennig, F. Henningsen, L. Heuermann, R. Hewett, N. Heyer, S. Hickford, A. Hidvegi, C. Hill, G. C. Hill, R. Hmaid, K. D. Hoffman, A. Hollnagel, D. Hooper, S. Hori, K. Hoshina, M. Hostert, W. Hou, M. Hrywniak, T. Huber, K. Hultqvist, K. Hymon, A. Ishihara, W. Iwakiri, M. Jacquart, S. Jain, O. Janik, M. Jansson, M. Jin, N. Kamp, D. Kang, W. Kang, A. Kappes, L. Kardum, T. Karg, A. Karle, A. Katil, M. Kauer, J. L. Kelley, M. Khanal, A. Khatee Zathul, A. Kheirandish, T. Kim, H. Kimku, F. Kirchner, J. Kiryluk, C. Klein, Y. Kobayashi, S. Koch, A. Kochocki, R. Koirala, H. Kolanoski, T. Kontrimas, L. Köpke, C. Kopper, D. J. Koskinen, P. Koundal, M. Kowalski, T. Kozynets, A. Kravka, N. Krieger, J. Krishnamoorthi, T. Krishnan, K. Kruiswijk, E. Krupczak, A. Kumar, E. Kun, N. Kurahashi, C. Lagunas Gualda, L. Lallement Arnaud, M. J. Larson, F. Lauber, J. P. Lazar, K. Leonard DeHolton, A. Leszczyńska, C. Li, J. Liao, C. Lin, Q. R. Liu, Y. T. Liu, M. Liubarska, C. Love, L. Lu, F. Lucarelli, W. Luszczak, Y. Lyu, M. Macdonald, E. Magnus, Y. Makino, E. Manao, S. Mancina, A. Mand, I. C. Mariş, S. Marka, Z. Marka, L. Marten, I. Martinez-Soler, R. Maruyama, J. Mauro, F. Mayhew, F. McNally, K. Meagher, A. Medina, M. Meier, Y. Merckx, L. Merten, J. Mitchell, L. Molchany, S. Mondal, T. Montaruli, R. W. Moore, Y. Morii, A. Mosbrugger, D. Mousadi, E. Moyaux, T. Mukherjee, M. Nakos, U. Naumann, R. Neshat, L. Neste, M. Neumann, H. Niederhausen, M. U. Nisa, K. Noda, A. Noell, A. Novikov, A. Obertacke, V. O'Dell, A. Olivas, R. Orsoe, J. Osborn, E. O'Sullivan, B. Owens, V. Palusova, H. Pandya, A. Parenti, C. Parisel, N. Park, V. Parrish, E. N. Paudel, L. Paul, T. Pernice, T. C. Petersen, J. Peterson, S. Pick, M. Plum, A. Pontén, V. Poojyam, B. Pries, R. Procter-Murphy, G. T. Przybylski, L. Pyras, C. Raab, J. Rack-Helleis, N. Rad, M. Ravn, K. Rawlins, Z. Rechav, A. Rehman, I. Reistroffer, E. Resconi, C. D. Rho, W. Rhode, L. Ricca, B. Riedel, A. Rifaie, E. J. Roberts, S. Rodan, M. Rongen, A. Rosted, C. Rott, T. Ruhe, L. Ruohan, D. Ryckbosch, J. Saffer, D. Salazar-Gallegos, P. Sampathkumar, A. Sandrock, G. Sanger-Johnson, M. Santander, S. Sarkar, M. Scarnera, M. Schaufel, H. Schieler, S. Schindler, L. Schlickmann, B. Schlüter, F. Schlüter, N. Schmeisser, T. Schmidt, F. Schmitt, A. Scholz, F. G. Schröder, S. Schwirn, S. Sclafani, D. Seckel, L. Seen, M. Seikh, S. Seunarine, P. A. Sevle Myhr, R. Shah, S. Shah, S. Shefali, N. Shimizu, M. Shin, B. Skrzypek, R. Snihur, J. Soedingrekso, D. Soldin, P. Soldin, G. Sommani, D. Song, C. Spannfellner, G. M. Spiczak, C. Spiering, J. Stachurska, M. Stamatikos, T. Stanev, T. Stezelberger, T. Stürwald, T. Stuttard, G. W. Sullivan, I. Taboada, S. Ter-Antonyan, A. Terliuk, A. Thakuri, M. Thiesmeyer, W. G. Thompson, J. Thwaites, W. Tian, S. Tilav, K. Tollefson, J. A. Torres, S. Toscano, D. Tosi, A. K. Upadhyay, K. Upshaw, A. Vaidyanathan, N. Valtonen-Mattila, J. Valverde, J. Vandenbroucke, T. Van Eeden, N. van Eijndhoven, L. Van Rootselaar, J. van Santen, J. Vara, F. Varsi, M. Velazquez, M. Venugopal, M. Vereecken, S. Vergara Carrasco, S. Verpoest, D. Veske, A. Vijai, J. Villarreal, C. Walck, A. Wang, E. H. S. Warrick, C. Weaver, A. Weindl, J. Weldert, A. Y. Wen, C. Wendt, J. Werthebach, M. Weyrauch, N. Whitehorn, C. H. Wiebusch, D. R. Williams, L. Witthaus, J. Woodward, G. Wrede, X. W. Xu, J. P. Yanez, Y. Yao, E. Yildizci, S. Yoshida, F. Yu, S. Yu, T. Yuan, S. Yun-Cárcamo, A. Zander Jurowitzki, A. Zegarelli, S. Zhang, Z. Zhang, P. Zhelnin, P. Zilberman, C. Zilleruelo Cañas

Applied Physics

A 28-GHz Varactor-Based RIS With Continuous Phase Control: From Unit-Cell Modeling to Programmable Wavefront Control and Synthesis

This paper presents a 28 GHz varactor-based reconfigurable intelligent surface (RIS) platform with continuous phase control and establishes a unified device-to-system validation framework for programmable wavefront control and synthesis. The proposed RIS comprises 96 independently controlled elements, each employing a single varactor diode, with a board-integrated analog-bias control architecture. It builds on an experimentally validated unit-cell model providing approximately 300{\deg} of continuous reflection-phase tuning at normal incidence. An analytical framework incorporating measured horn illumination, finite phase availability, and unit-cell reflection losses consistently relates device-level characteristics to beamforming performance. It is evaluated via near-field-to-near-field characterization, near-field-to-far-field beam-steering, and far-field-to-far-field wireless-link experiments. The near-field-to-far-field results show close agreement among analytical predictions, full-wave simulations, and measurements, while the far-field-to-far-field response agrees with simulation and a first-order link-budget estimate. Accurate steering is demonstrated for all investigated angles within $\pm$45{\deg} across three azimuthal planes, with maximum deviation of approximately 2{\deg}. The complete prototype, including driver and bias network, draws only 0.85 W with an estimated full-aperture reconfiguration time of approximately 50 ms. Beyond beam steering, the same platform enables experimental investigation of 3-bit, 2-bit, and 1-bit phase quantization and programmable multi-beam wavefront synthesis using a common RF aperture and control architecture. Collectively, these results bridge realistic varactor behavior, analytical modeling, and programmable wavefront synthesis, providing a rigorous basis for developing and experimentally validating continuously tunable millimeter-wave RISs.
Read the full abstract

All authors Spandan Manna, Florian Reher, Karim El Isa, Amar Al-Bassam, Dirk Heberling

Fluid Dynamics

Emergence of chaos and fractality in the basin boundary of subcritical shear flow

From a dynamical systems perspective of subcritical transition in shear flows, the basin boundary separating the laminar and turbulent attractors, along with the edge state that governs the long-term dynamics on that boundary, are of fundamental interest. As the Reynolds number is increased from small values, a multiplicity of simple exact coherent structures (ECS) appear in phase space, of which one often undertakes the role of the edge state. At higher values of the Reynolds number, however, the dynamics on the basin boundary and the edge state are sensitive to initial conditions, as shown by a wealth of numerical and experimental studies. The mechanism behind this transition remains unclear. To address this, we examine the subcritical regime of Taylor-Couette flow using a minimal computational box and reveal a generic mechanism whereby the edge state becomes chaotic. The first step in this process involves the formation of a heteroclinic tangle between a travelling-wave-type ECS and an independently engendered chaotic saddle. The interaction causes the basin boundary to incorporate the saddle, thus inheriting its fractal structure, while the edge state itself remains the simple ECS. The transition of the edge state to chaos requires a second step: the formation of a heteroclinic cycle involving the ECS and the chaotic saddle. In consequence, the observation of a simple non-chaotic edge state at given values of the parameters is no guarantee that the same situation will hold at nearby values.
Read the full abstract

All authors Baoying Wang, Roger Ayats, Kengo Deguchi, Alvaro Meseguer, Fernando Mellibovsky

Chemical Physics

Density-Functional Excited-State Gradients and Nonadiabatic Couplings on a Consumer GPU from a Contraction-DAG

Nonadiabatic dynamics needs an excited-state gradient and an interstate nonadiabatic coupling matrix element (NACME) at every nuclear geometry, and a double-hybrid functional's accuracy has been unavailable for the coupling. We report the first analytic derivative NACME for a double-hybrid excited state—deferred in the original hh-TDA method and supplied for hybrids only by Yu et al.—derived, with the hole-hole and particle-particle Tamm–Dancoff (\hhTDA/\ppTDA) gradients and NACMEs, as a single reverse-mode transpose of one contraction graph closed under a non-symmetric atomic-orbital-direct $J/K$ kernel. Its double-hybrid excitation energy lowers the vertical-excitation mean absolute deviation from bare-\hhTDA\ $0.86$ to $0.47$~eV and removes the $+0.53\!\rightarrow\!+0.05$~eV over-excitation bias, improving seven of ten states while over-correcting the ionic $\pi\pi^*$ states—the expected perturbative-doubles failure, reported not trimmed. Every coupling is validated to $\sim\!10^{-4}$ against an independent literal many-electron wavefunction-overlap oracle that shares no code path with the method, and is physically meaningful at the ammonia $n\!\rightarrow\!\sigma^*$ covalent conical intersection, where the \hhTDA/\ppTDA manifolds recover the $F\!-\!2$ seam and adiabatic linear-response TDDFT gives $\tau\!\equiv\!0$ by construction. Gradients, NACMEs, and the double-hybrid coupling all run device-resident and AO-direct through one shared Cholesky-decomposed $J/K$ engine within the 8\,GB of a consumer RTX~4060 (a profile-guided $\sim\!10^2\times$ launch collapse preserving double-precision bit-identity)—placing on a commodity desktop card a correlated excited-state derivative capability that has until now required datacenter hardware.
Read the full abstract

Instrumentation and Detectors

Gamma Neutron Radioactive Source Identification in Water Cherenkov Detectors

Water Cherenkov Detectors (WCDs) are a robust technology widely used in astrophysics, high energy physics, and recently nuclear security applications. They detect high energy interactions through the Cherenkov light emitted by charged particles traveling faster than the speed of light in water. In this work, we demonstrate the feasibility of gamma-neutron discrimination in WCDs using a combined methodology that integrates statistical analysis with machine learning techniques. The experimental setup employs different shielding configurations to isolate gamma and neutron contributions from a \textsuperscript{241}AmBe source, while \textsuperscript{60}Co and \textsuperscript{137}Cs sources are used to establish a signal to energy calibration. A statistical analysis based on a $3\sigma$ significance criterion is used to define energy thresholds, enabling a linear relationship between the measured charge spectrum and the deposited energy. Building on this calibration, pulse shape information is further exploited through machine learning methods to improve event classification. An ensemble model based on a soft-voting strategy combining a Bagging classifier, CatBoost, and a Multilayer Perceptron was trained on detector signals acquired under different shielding conditions, achieving an accuracy of 0.816 and an area under the Receiver Operating Characteristic (ROC) curve. The combined approach demonstrates that statistical thresholding provides a physically grounded discrimination baseline across the full energy range, while machine learning enhances classification performance at higher energies by leveraging pulse level information. This integrated strategy improves radiation identification capabilities in water Cherenkov detectors, with potential applications in nuclear security and radiation detection.
Read the full abstract

All authors A. Núñez Selin, C. Sarmiento Cano, H. Asorey, I. Sidelnik

Applied Physics · Cross-listed from astro-ph.EP

Asteroid Disruption and Deflection Simulations for Multi-Modal Planetary Defense

Planetary defense from asteroids via deflective means alone does not offer viable solutions in terminal scenarios where there is little warning time before impact. The PI method of planetary defense enables operation in terminal interdiction modes where there is little warning time prior to impact, but can also operate in the same extended time scale interdiction modes as made possible by traditional deflection techniques, which results in a versatile, multi-modal planetary defense capability. The method is also practical and cost-effective since it relies solely on launch vehicles and penetrator materials already available today, and thus presents itself as a logical and competitive option for planetary defense. As per the PI method, we investigate the effectiveness of rubble pile asteroid disruption and deflection via hypervelocity impacts with 10:1 aspect ratio cylindrical tungsten penetrators. We present the results of an ongoing simulation campaign dedicated to investigating the PI method, using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian-Eulerian (ALE) hydrodynamics code ALE3D run with the High-End Computing Capability (HECC) at NASA Ames Research Center. We model heterogeneous rubble pile asteroids with a distribution of spherical boulders of varying initial yield strengths set within a weak binder material. We find that rubble pile asteroids of this type in the 20 - 100 meter-class can be effectively mitigated via 20 km/s impacts with 100 - 1000 kg penetrators via the coupling of the penetrator kinetic energy into the bulk material of the asteroid.
Read the full abstract

All authors Alexander N. Cohen, Philip Lubin, Darrel Robertson, Mark Boslough, Sasha Egan, Brin Bailey, Angela M. Stickle, Elizabeth A. Silber, Peter Meinhold, Dharv Patel

Optics

The Quantitative and dynamical analysis of anisotropic dispersive optical solitons

In this article, we explore the integrable aspects of generalized non-linear fluid equation possessing mixed spatial and time evolutions. This non-linear model has been acknowledged as integrable tool for the non-linear dynamical character- ization of dense material such as in optical fibre, metallic ionic fluid. This model widely has been applied in physical interpretations of electromagnetic propaga- tion through perturbed charged plasma. Here, we use few integrable approaches and attain a broad range of solitary solutions involve deep analytical analysis. These results are furnished with computational simulations, which enrich under- standing about the accurate pictorial flow of energy. We also present a bifurcation analysis for canonical coordinates to examine the impact of parameter variations on non-linear behavior of the systems under consideration. Subsequently, we per- form stability analysis to evaluate the robustness of optical solitons features in the vicinity of equilibrium points. These results establish the reliable and accurate framework for analyzing non linear response of materials for energy propagation by varying physical parameters associated with this non-linear equation model.
Read the full abstract

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

Multi-wavelength synthesis of a flux rope-trapped mini-prominence eruption and post-flare coronal rain

Small-scale eruptive phenomena in the solar corona including miniature flux ropes and associated cool plasma condensations are not fully understood despite increasing high-resolution observations. We perform forward modeling based on a 2.5D MHD simulation capturing homologous flux rope eruptions, in-situ condensation leading to a mini-prominence, and subsequent post-flare coronal rain. Synthetic diagnostics are obtained using optically-thin EUV and UV emissions, and non-LTE radiative transfer treatment for the H$\alpha$ line. The synthetic EUV emission reveals the flux ropes as bright rim-like structures. The corresponding UV diagnostic shows bright region, which is co-spatial with the dark core due to embedded cool plasma ($\sim$ tens of kK) inside the flux rope, identifying an erupting mini-prominence. Spectral synthesis of Si IV 1402.77 A indicates an upward motion of the mini-filament, and reveals the presence of two predominant velocity components during eruption. At a later stage, thermal instability in post-flare arcades produces coronal rain with temperatures of $\approx 10^4$ K. The EUV diagnostics reveal brightening at the downstream of the rain blob, indicating localized heating associated with compressional effects. The H$\alpha$ spectral synthesis shows enhanced absorption signatures and red-shifted profiles corresponding to downflows of the coronal rain blobs up to $\approx 23$ km s$^{-1}$, whereas the Si IV 1402.77 A spectral profile shows the maximum downflow velocity of $\approx50$ km s$^{-1}$, highlighting the evidence of thermodynamic and kinematic structuring within the falling rain blobs. The synthetic diagnostics provide clear, multi-wavelength signatures that can guide future high-resolution observations, and highlight the importance of small-scale reconnection-driven processes in shaping the multi-thermal structure (between MK to kK) of the solar corona.
Read the full abstract

All authors Samrat Sen, Alexander G. M. Pietrow, Veronika Jercic, Patrick Antolin

Chemical Physics

Interference between multiple photoionization pathways in chiral molecules: Converging continuum results in Gaussian bases

An accurate description of photoionization observables is a central challenge for theoretical models of molecular photoionization. Within standard electronic-structure approaches, the continuum states are represented by unoccupied Hartree-Fock orbitals expanded in Gaussian basis sets. Since these basis sets are optimized for bound states, computed observables may exhibit a noticeable basis-set dependence. Here, we augment these basis sets with diffuse functions and investigate the convergence of photoelectron circular dichroism (PECD), anisotropy parameters, and the forward backward ionization time delay in the multiphoton ionization of randomly oriented chiral molecules. All observables converge systematically with the number of added diffuse functions, resolving previously observed basis-set discrepancies and indicating that the augmented basis sets provide a more accurate representation of the intermediate continuum states. In particular, our fully ab initio calculations yield forward-backward time delays in qualitative agreement with recent measurements, for which previous theoretical descriptions relied on empirical modeling. Diffuse augmentation thus provides a computationally efficient and transferable route to converged Gaussian-basis calculations for molecular multiphoton photoionization.
Read the full abstract

All authors Muhammad Sakhi, Alexander Blech, Corbin Allison, Christiane P. Koch, Loren Greenman

Physics and Society

Game-Theoretic Inverse Reinforcement Learning for Modeling Competitive Human Driving: A Cut-in Prediction Study

Capturing the strategic decision-making inherent in competitive human driving is critical for autonomous vehicle safety and traffic simulation. This study demonstrates that game-theoretic Inverse Reinforcement Learning (IRL) provides a robust framework for this challenge. We present a comprehensive analysis comparing data-driven IRL models against an established physics-based game-theoretic approach for predicting aggressive, safety-critical cut-in lane changes. Using the high-fidelity highD dataset, we systematically develop and evaluate a series of IRL models with increasing feature complexity. Our results reveal significant advantages: the best-performing IRL models achieve an overall prediction accuracy exceeding 75 percent while maintaining a Cut-In precision up to 51 percent and recall up to 49 percent. This represents a significant improvement over the established physics-based benchmark, which achieved only 4.4 percent precision in these high-stakes scenarios. The analysis reveals a clear trade-off: incorporating granular, instantaneous features yields higher precision, while adding temporal consistency features maximizes recall. These findings suggest that IRL-based models can effectively bridge the gap between microscopic driver intent and macroscopic safety outcomes, providing a more reliable foundation for modeling interactions in mixed-autonomy environments.
Read the full abstract

Atmospheric and Oceanic Physics · Cross-listed from astro-ph.IM

Atmospheric Light Pollution by Proposed Reflect Orbital Space Mirrors

ReflectOrbital is planning to launch a $18\times18$\,m mirror in space at a $\sim$600\,km altitude to illuminate a $2.5$\,km radius circular patch on Earth. The pilot satellite, called EARENDIL-1, is a pathfinder for a constellation of $\sim$50,000 larger mirrors, each with $54\times54$\,m in diameter, with the explicit goal to provide space-borne illumination on Earth. We calculate the sky brightness caused by the Rayleigh scattering and aerosol scattering of the incoming beam from one such satellite and by the light reflected back by the illuminated patch, for various atmospheric properties and ground albedos under cloud-free conditions. We demonstrate that the light pollution caused by one of these satellites is very significant, altering the nighttime environment up to distances of $\sim$30\,kms. The results show that for an observer within the beam of a single 54\,m satellite, the mirror will appear as a $-16.7$\,mag point-like source, i.e.~about 4 magnitudes brighter than the full moon. The diffuse sky background will be similar to the dusk sky shortly after sunset; too bright to see even the brightest stars. From a distance of 14\,km, the glow from a single mirror will exceed the luminance of the full moon sky for majority of the sky. From a distance of 34\,km the sky will still appear brighter than the moon-lit sky in the direction of the beam. If 400 such mirrors illuminate the same patch simultaneously, the glow will be obvious from 80\,kms.
Read the full abstract

Biological Physics · Cross-listed from cond-mat.stat-mech

A Universal Control Budget for First-Passage Kinetics

The first-passage time is the natural observable of reaction completion, yet how its mean responds to a rate change has lacked a general constraint. We show that the logarithmic sensitivity of the mean first-passage time of any finite Markov chain to any rate is bounded by one in magnitude, and that these sensitivities sum to -1. Together the two laws form a conserved control budget: speeding completion through some transitions must be paid for by others, and a coordinated change shifts the completion time only as far as the budget allows. Raising an activation barrier or shifting the depth of a well moves many rates at once, yet neither can shift the completion time further than a single rate could. The budget caps kinetic-proofreading discrimination at the checkpoint count, and prices it in sensitivity to substrate concentration.
Read the full abstract

Optics

Structured coherence: A modern perspective on optical coherence as a resource

Optical coherence is a well-established branch of physical optics in which the statistical properties of fluctuating optical fields are described in terms of correlation functions over continuous spatial and temporal degrees of freedom (DoFs). Nevertheless, in any practical setting, only discrete DoFs are ever accessible. In these scenarios, the modes are fixed, stable, and deterministic, and partial coherence arises from random relative complex amplitudes, a configuration we refer to as 'structured coherence'. Advances in structured coherence have recently unveiled new conceptual ground in optical communications and information processing in which partial coherence may be preferable to full coherence, which we call a 'coherence advantage'. We formulate structured coherence in terms of coherence matrices to investigate these recent theoretical and experimental breakthroughs. We first review optical fields characterized by a binary DoF via 2x2 Hermitian coherence matrices and introduce key concepts that take on new significance for larger-dimensional modal sets. Next, we examine the structured coherence of two binary DoFs, which can be described by 4x4 coherence matrices. We introduce the concept of coherence rank, entropy swapping, and optical cross-purity. In the perspective outlined here, coherence is viewed as a 'resource', which can be exchanged between DoFs, concentrated into a DoF or into particular modes, or spread over the DoFs. We then examine larger-dimensional modal sets, which allow for more versatile applications in optical information processing. The formulation presented here lends itself particularly to the manipulation of partial optical coherence in integrated photonic platforms, thereby opening myriad avenues for novel fundamental investigations of structured coherence and exploiting the coherence advantage in optical communications and information processing.
Read the full abstract

Chemical Physics

Nuclear-Spin Statistical Weights from Young Diagrams

Nuclear-spin statistical weights and selection rules in molecular spectroscopy are governed by the permutation symmetry of identical nuclei, although rigid and semi-rigid molecules typically realize only a proper subgroup of the full symmetric group. Building on the Schur-Weyl framework of Schmiedt, Jensen, and Schlemmer, we extend this approach to molecular geometries whose rovibrational motion realizes cyclic and dihedral permutation subgroups. By combining the Schur-Weyl decomposition of the $n$-spin Hilbert space with the Kraśkiewicz-Weyman-Adin-Roichman major-index branching rule for the restriction $S_N \downarrow C_m$, we obtain a fully combinatorial method for determining nuclear-spin symmetry species and statistical weights. Each standard Young tableau corresponds to a definite cyclic character determined by its major index, with reflections yielding the associated dihedral symmetry, so that nuclear-spin species follow directly from tableau data without projection operators or case-specific constructions. Applications to $XeOF_4$, $SF_4$, benzene, deuterated benzene, and the tropylium cation demonstrate that the method applies uniformly to realistic molecular point groups and reproduces established statistical weights while providing a transparent spin-resolved structure. The approach supplies the symmetry information required for rovibrational line-intensity modeling and provides a practical, automatable framework for nuclear-spin symmetry analysis in polyatomic molecules.
Read the full abstract

Medical Physics

An open-source framework for predicting ultrasound neuromodulation: bridging tissue elastomechanics and neuron firing dynamics

Transcranial focused ultrasound is a non-invasive neuromodulation modality with millimetre-scale resolution, but its biophysical mechanism of action remains unresolved. Exposure is conventionally specified by transducer surface or derated focal pressure, quantities only indirectly related to what matters for therapy: which neurons fire and through which pathway. We address this gap with an end-to-end computational framework that maps a transcranial acoustic field to per-voxel neural firing maps registered to anatomy. The pipeline couples heterogeneous nonlinear full-wave acoustic propagation, viscoelastic shear-wave propagation, Pennes bioheat diffusion, a bilayer-mechanics conversion from tissue strain to membrane tension, and a multi-compartment Hodgkin-Huxley neuron carrying mechanosensitive, cavitation-coupled, calcium-coupled, thermosensitive, astrocytic-gliotransmitter, and mechanosensitive-synaptic pathways. Six candidate mechanisms are implemented as interchangeable components on a shared neuron model, so their firing predictions can be compared directly on the same field, and every numerical parameter is classified by source and bracketed by sensitivity analysis. We demonstrate the framework on a theta-burst sonication delivered through a micro-CT human-skull specimen targeting the left dorsal anterior cingulate cortex, predicting a focal firing zone of approximately 8,500 mm^3 at a focal thermal rise well within ITRUSST consensus safety envelopes. The principal output is a per-voxel firing-volume map resolved jointly with the acoustic, elastic, and thermal field histories that drive it, giving spatially resolved, falsifiable predictions that are testable against high-density extracellular recordings and support parameter estimation, cell-type-resolved mechanism identification, and quantitative safety assessment for ultrasound neuromodulation.
Read the full abstract

Fluid Dynamics

Three-layer water flows: Dirichlet-Neumann operators and approximations

The object of investigation in this paper are the nonlinear equations of motion for two-dimensional inviscid water flows with piecewise constant density stratification in a three-layer fluid with a flat bottom, a free surface and two interfaces. We establish a Hamiltonian formulation for the nonlinear governing equations in this setup. The Hamiltonian of the system and the equations of motion of the surface and of the interfaces are expressed with the help of the Dirichlet-Neumann (DN) operators, which are introduced for each of the layers. Then, the linear equations for small amplitudes of the elevation of the surface and of the interfaces in the leading order are derived from which a bi-cubic equation for the dispersion relation is obtained, whose solutions are analysed. The six real solutions for the possible propagation speeds (three positive, related to right-moving waves and three negative, related to left-moving waves) have magnitudes of different order. Upper and lower bounds for the previously mentioned roots are also given in terms of the coefficients of the equation. Subsequently, approximate formulae for the propagation speeds are derived. The importance of the DN operators is further illustrated in a separate analysis of the three-layer model with flat surface (rigid lid). The full nonlinear evolution equations are expressed again in terms of the DN operators, and the equations in the linear regime and the weakly nonlinear propagation regime (the Boussinesq approximation) are derived by a proper expansion of the DN operators. Limits to the two-layer free surface model are obtained as well. The obtained results are applicable to internal waves in lakes and in the ocean as well as to laboratory experiments with three superimposed fluid layers.
Read the full abstract

Computational Physics · Cross-listed from astro-ph.IM

NestyNet. II. Coherent Function-Space Posteriors from Scientific Neural Surrogates (or How to Avoid Expensive MCMC)

Scientific analyses increasingly use flexible neural networks, but their thousands of correlated parameters make it challenging to interpret the associated uncertainties. Here we develop a low-dimensional posterior for the fitted function itself, for scientific neural surrogates trained with second-order optimization. Linearizing the fitting procedure with respect to the randomized residual rows gives a measurement-to-function transport, the linear map, assembled from the converged Jacobians and Gauss–Newton curvature, that carries measurement perturbations into the function perturbations that refitting would produce. Its leading singular functions define coherent deformation modes. Independent Gaussian coefficients then generate smooth function draws, so that any derived quantity, including those requiring derivatives or integrals of the draw, inherits the posterior. The construction distinguishes repeated-experiment covariance from the local Gauss–Newton/Laplace posterior and propagates both to correlated quantities of scientific interest. The result is conditional on the fit's declared choices (architecture, hyperparameters, active set, and optimization branch), and every fit is certified as converged by checking that a further optimization step would change the fitted predictions by less than a chosen small fraction of the measurement errors. Our primary example is an 800-parameter phase-space distribution function fit for a mock stellar disk. Four uncertainty coordinates, two orders of magnitude fewer than the fitted parameters and stable under refinement of the force basis, capture 99% of the vertical-force posterior variance, and 4000 coherent draws propagate through the force, total-density, surface-density, and frequency calculations in 0.8s. The method provides a highly efficient route to uncertainty propagation for derivative-dependent scientific inference.
Read the full abstract

All authors Rodrigo Ibata, Wassim Tenachi, Foivos Diakogiannis, Neil Ibata, Anirudh Shankar

History and Philosophy of Physics

Newton's Second Law: A Theoretical Identity Derived from the Principle of Excluded Perpetual Motion and the Weak Equivalence Principle

For more than three centuries, Newton's Second Law (F=ma) has governed mechanics with undisputed success, yet its epistemic authority has rested on empirical adequacy alone. That empirical contingency vanishes under two physical principles - the Principle of Excluded Perpetual Motion (PEPM) and the Weak Equivalence Principle (WEP) - from which the law emerges as a structural necessity of admissible mechanics. Under the PEPM constraint, Suppes' operational measurement protocol grounds gravitational mass and force as independent primitives. The fusion of Stevin's static and Galileo's kinematic inclined-plane analyses establishes F/a as a well-defined, object-intrinsic quantity - the operational definition of inertial mass. The WEP compels the equivalence of inertial and gravitational mass without presupposing Newton's Second Law. Substituting this equivalence into the definition of inertial mass yields the law as a theoretical identity between independently defined quantities. The derivation is anchored in gravity because weight uniquely bridges statics and kinematics, yet the result is independent of the underlying force mechanism. The derivation carries structural consequences: modified-inertia formulations of MOND are inadmissible under the PEPM-WEP constraint; the location-invariance of the Kibble-balance realization of the SI kilogram is secured on first principles; and under the PEPM alone, torsion-balance and free-fall experiments constitute equally direct tests of the WEP. More fundamentally, Newton's Second Law is not an irreducible axiom of mechanics, but a structural consequence of deeper physical principles: energy conservation and the Einstein Equivalence Principle.
Read the full abstract

Chemical Physics

A Phase Space Electronic Structure View of The Solid State

We develop a phase space electronic structure view of the solid state, where electronic bands are parameterized by both nuclear position $\mathbf{Q}$ and nuclear momentum $\mathbf{\Pi}$ (rather than the standard practice of only $\mathbf{Q}$). Within this phase space, non-Born-Oppenheimer electronic structure context, we prove a nuclear wavevector ($\mathbf{q}$)-dependent version of Nafie's equality relating the change in electronic momentum as a function of nuclear momentum, $\partial \langle\mathbf{p}\rangle/\partial \mathbf{\Pi}_q$, to the change in electronic position as a function of nuclear position, $\partial \langle\mathbf{r}\rangle/\partial \mathbf{Q}_q$. Furthermore, we show how information about electron-phonon interactions, specifically nuclear-induced electronic momentum, can be extracted from simple band structure calculations – without Berry curvature calculations.
Read the full abstract

Physics Education

Using LLMs to Detect Growth in Computational Thinking in Introductory Physics

As computation becomes more central to physics education, creating scalable methods to assess authentic computational thinking (CT) in students remains a critical challenge. While student-written responses capture nuanced reasoning, they are difficult to evaluate at scale. In this study, we investigated the use of Large Language Models (LLMs) to analyze students' written explanations of computational physics problems on a pre- and post- semester survey. By first establishing a human-coded baseline, grounded in CT literature, we identified significant growth in Data Practices and Computational Problem-Solving Practices. When given the same responses, an LLM successfully mirrored the human evaluations and scaled up the detection of these key trends across a large dataset. Notably, both human raters and the LLM struggled to reliably evaluate more complex constructs such as Systems Thinking. Overall, this study demonstrates that LLMs offer a viable method to scale the evaluation of students' CT in large-enrollment physics courses
Read the full abstract

All authors Sean Savage, Anand Shanker, Grace Michlitsch, N. Sanjay Rebello

Fluid Dynamics

VENUSS: a unified finite-element model of solidifying lava

The development of a solid rind or carapace at the surface of lava flows and domes results in a transition in deformation mechanism from dominantly viscous to elastic or plastic. This transition has a significant impact on the rate and style of emplacement, including on the construction of channelized flows, over-steepened margins, and flow advance due to lava breakouts. These processes are particularly important in subaqueous, subglacial, and extraterrestrial environments in which cooling is accelerated, requiring models specifically calibrated for these environments. We present a new numerical model, Viscous-Elastic Numerically Unified Solver for Solidifying flows (VENUSS), for cooling and solidifying free surface flows. The model couples a viscous fluid interior with an elastic shell whose thickness grows in response to cooling. As a demonstration of the impact of including a solidified crust in the flow model, we show that a dome-like shape fed from below with an elastic shell coupled to the basal topography results in more lateral expansion and less vertical uplift than a comparable highly-viscous rind, demonstrating the need for lateral transfer of stress in solid layers to accurately interpret and predict dome deformation.
Read the full abstract

All authors Janine Birnbaum, Einat Lev, Marc W. Spiegelman, Jackie E. Kendrick, Yan Lavallee

Optics · Cross-listed from eess.SP

Numerical Model of a Multiple-Input-Multiple-Output Distributed Acoustic Sensing System with Joint Phase and Birefringence Estimation

In this work, we introduce and experimentally validate a numerical model for a Multiple-Input-Multiple-Output Distributed Acoustic Sensing (MIMO-DAS) system that accounts for dynamic perturbations of fiber birefringence and of the common optical phase of the backscattered signal (or polarization-averaged phase, shared by both polarization tributaries). The MIMO-DAS system probes the fiber using polarization-multiplexed constant-power coded sequences that are suited for coexistence of DAS with WDM data transmission over the same fiber. We study the effect of both axisymmetric and anisotropic events on the two quantities. We demonstrate, through numerical simulations and lab experiments, the estimation of effective birefringence magnitude in static conditions, and the joint estimation of common phase and effective birefringence magnitude in the case of dynamic longitudinal strain and anisotropic transverse strain. This allows for event discrimination and increased sensitivity to disturbances that act transversely on the fiber, since polarization will be responsive to perturbations that break cylindrical symmetry, while the phase will strongly respond to longitudinal strain.
Read the full abstract

All authors Diane Prato, Mehran Mokhtari Sheramin, Renaud Gabet, Elie Awwad

Atomic Physics · Cross-listed from cond-mat.quant-gas

Quantum transport in a non-Hermitian 1D synthetic lattice: from quantum Zeno reflection to near-perfect absorption

We experimentally explore the absorption of a propagating wavepacket impinging upon a dissipative region in an engineered quantum system. We employ a 1D non-Hermitian synthetic lattice with an abrupt interface between dissipative and non-dissipative subchains, using the states of the electronic ground-state hyperfine manifold in a $^{87}$Rb Bose-Einstein condensate as sites. By tuning the dissipation rate, we observe a progression from ballistic propagation, to near-perfect absorption, to quantum Zeno reflection. Guided by numerical simulations, we identify that optimal absorption occurs when tunneling and dissipation are properly matched, and find qualitative agreement with an idealized semi-infinite model across all dissipation regimes. Our results establish synthetic lattices as a versatile Quantum simulation platform for dissipation-engineered quantum transport and highlight controlled dissipation as a resource for tailoring quantum dynamics.
Read the full abstract

All authors Emmanuel. D. Mercado-Gutiérrez, W. Xu, E. Gvozdiovas, A. M. Piñeiro, Javier Argüello-Luengo, W. D. Phillips, Q. Zhou, I. B. Spielman

Optics

Differentiable eigendecomposition-free RCWA for full-tensor anisotropic photonics

Full-tensor anisotropy transforms rigorous coupled-wave analysis (RCWA) into a large, fully coupled non-Hermitian eigenproblem, making eigendecomposition expensive and difficult to differentiate. We introduce a differentiable, eigendecomposition-free RCWA framework for spatially patterned media with fully coupled permittivity tensors, using boundary fields rather than internal eigenmodes as the layer representation. A boundary-value cascade constructs scattering operators directly from these fields, enabling automatic differentiation and efficient GPU execution. Benchmarks against finite-element and transfer-matrix solutions show close agreement in scattering responses, while automatic-differentiation gradients agree with finite differences and enable topology optimization. At 529 Fourier harmonics, layer construction is 22.8 times faster than conventional eigendecomposition on the same GPU. Our framework offers a general computational route toward scalable forward modeling and inverse design in anisotropic photonic systems.
Read the full abstract

Fluid Dynamics

From expansion to collapse: Bubble and continuum multiscale modeling in open-system magmas

Bubble growth in silicate melts drives significant volume expansion, which has a first order control on magma transport dynamics. When magmas are exposed to external environments, heat and volatile loss at free surfaces can reverse bubble growth, leading to shrinkage and complex feedbacks between diffusion, rheology, and flow. To resolve how magma flow controls, or is controlled by, bubble expansion, we couple a micro-mechanical model for volatile diffusion into individual bubbles, with a macro-scale thermal evolution and fluid flow of the surrounding magmatic suspension. This two-way coupling captures the co-evolution of bubble size, melt viscosity, and pressure gradients, allowing both growth and resorption to emerge naturally from local conditions. We identify distinct dynamical regimes governed by (i) bubble growth limited by (a) viscous resistance or (b) diffusion at the bubble scale, (ii) viscous transport of the suspension, (iii) outgassing through permeable porous networks and exposed magma-fluid interfaces, and (iv) thermal quenching. Across these regimes, thin, high-viscosity boundary layers arising from temperature and volatile concentration gradients play a central role in modulating flow and bubble evolution. The model is implemented in a flexible, modular numerical framework (Multiscale Vesiculation, Fluid flow, Failure, and Interaction Nonlinear model: MVFFIN) enabling extension to a wide range of systems and applications, including conduit flow and pyroclast evolution. By resolving the interplay between internal bubble dynamics and external boundary conditions, this approach provides a unified framework for understanding multiscale degassing and its impact on magmatic transport and fragmentation.
Read the full abstract

All authors Janine Birnbaum, Fabian B. Wadsworth, Anthony Lamur, Jackie E. Kendrick, Yan Lavallee

Physics and Society

The Structure of Spreading on Temporal Networks

The physics of spreading in static networks is well understood through mappings to percolation. We show that spreading dynamics on temporal networks can analogously be mapped to reachability in temporal event graphs. This provides a theoretical and computational framework for a class of processes, such as variants of the susceptible-infected-susceptible model. Without explicit simulations, through the component analysis of event graphs, we obtain epidemic prevalence and derive epidemic thresholds for temporal networks with arbitrary degree and inter-event time distributions, with significant computational advantages as compared to explicit simulations.
Read the full abstract

Atomic Physics

Dual-Faraday-laser-pumped cesium beam clock with $7.7\times 10^{-13}/\sqrt\tau$ frequency stability

Compact cesium beam clocks are major frequency references for deployable timing systems. However, further improvement of their short-term frequency stability is limited by the clock signal-to-noise ratio (SNR). Although two-laser optical pumping can increase the effective atomic utilization, the achievable clock SNR has long been limited by laser-induced frequency-to-amplitude noise conversion. Here, we demonstrate a compact dual-Faraday-laser-pumped (DFP) Cs beam clock enabled by a low-frequency-noise atom-referenced laser architecture. The intracavity Faraday anomalous dispersion optical filter provides inherent alignment to the Cs D$_2$ resonances, while modulation transfer spectroscopy offers suppressed frequency noise and drift. The resulting laser system supports robust turnkey operation with a Lorentzian linewidth of 2.12 kHz. The DFP Cs clock achieves a clock SNR of 46,365 in a 1-Hz bandwidth and a fractional Allan deviation of $7.7\times 10^{-13}/\sqrt\tau$ , with Hadamard deviation reaching $7.7\times 10^{-15}$ at 10,000 s. This work pushes the fractional frequency stability of a compact Cs beam clock into the $10^{-13}/\sqrt\tau$ regime, providing a pathway toward high-performance Cs frequency references for field-deployable precision timing, navigation, and synchronization.
Read the full abstract

All authors Xiaomin Qin, Suyang Wei, Haijun Chen, Yufei Yan, Qiang Wei, Hangbo Shi, Zhiyang Wang, Zheng Xiao, Zijie Liu, Tiantian Shi, Jingbiao Chen

Atomic Physics

Simulation of positronium laser cooling using the Lindblad master equation

We present a formulation and numerical results for positronium (Ps) laser cooling. The formulation is based on the Lindblad master equation and follows the time evolution of the density matrix of Ps atoms. It therefore accounts for atomic coherence, which is necessary to describe the interaction of Ps with the train of short laser pulses generated by the system developed by Shu $\textit{et al.}$ [K. Shu $\textit{et al.}$, Phys. Rev. A $\textbf{109}$, 043520 (2024)]. Using this formulation, we calculate the time evolution of the populations in each internal and momentum state and thereby quantitatively predict the momentum distribution after laser cooling. We present the representative time evolution of the internal-state populations and momentum distribution, together with a comprehensive scan of the laser parameters used to optimize the cooling efficiency. A prominent feature of the simulated distributions is sub-recoil cooling through velocity-selective coherent population trapping, a coherent effect captured by the quantum-mechanical treatment.
Read the full abstract

Medical Physics

Evaluation of Silicon-Based Photon-Counting CT for Coronary Stenosis Quantification with Realistic Coronary Artery Phantoms

Objective: To quantify the impact of high-resolution deep silicon photon-counting CT (dSi-PCCT) on coronary stenosis quantification in anatomically realistic calcified coronary artery phantoms using Micro-CT as ground truth. Methods: Twelve vessel sections representing four calcification geometries (Type I-IV) and three luminal iodine concentrations (10, 15, and 20 mg/mL) were scanned under static conditions using energy-integrating detector CT (EID-CT), dSi-PCCT, and Micro-CT. Images were registered to Micro-CT and segmented using an automated threshold-based pipeline. The primary analysis compared longitudinal profiles of Micro-CT-referenced percent area stenosis and segmented vessel area. A secondary analysis evaluated ellipse-derived percent area stenosis and percent vessel-area deviation at the maximum-calcification cross-section. Results: dSi-PCCT reduced whole-profile mean absolute error in Micro-CT-referenced percent area stenosis from 3.10% with EID-CT to 1.62% (p=0.027) and reduced segmented vessel-area error from 0.55 to 0.31 mm2 (p=0.001). In the secondary analysis, absolute deviations in ellipse-derived percent area stenosis ranged from 0.1% to 6.5% for dSi-PCCT and from 0.8% to 24.2% for EID-CT (p<0.001). Mean absolute differences in percent vessel-area deviation from Micro-CT were also lower with dSi-PCCT than with EID-CT (15.0% vs 26.6%, p<0.001). Conclusion: Under static, resolution-optimized conditions, dSi-PCCT improved task-based coronary stenosis quantification and vessel delineation relative to EID-CT, supporting further evaluation in dynamic phantoms and clinical CCTA.
Read the full abstract

All authors Maria Jose Medrano, Jed Pack, Stephen Araujo, Grant M. Stevens, Koen Nieman, Ge Wang, Bruno De Man, Adam S. Wang