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)

Tuesday 13 October · still arriving

Chemical Physics · Cross-listed from cond-mat.mes-hall

WannierSTM: Efficient superconducting scanning-probe spectroscopy from maximally localized Wannier functions

We present a Python-based program that combines maximally localized Wannier functions-based first-principles electronic structure with the nonequilibrium Green's function formalism to enable efficient simulations of scanning tunneling microscopy (STM). The workflow builds upon standard density functional theory packages to construct compact transport Hamiltonians that preserve the properties of realistic materials. Superconductivity is incorporated via an effective pairing potential within the Nambu-spinor formalism, while coupling to the STM tip and substrate electrodes is treated using the wide-band limit approximation. Within this formulation, transport quantities such as current-voltage characteristics and differential conductance are evaluated in a closed analytical form. This avoids explicit numerical energy integrations and enables rapid spectroscopic calculations over broad parameter spaces. The framework provides a unified and computationally efficient approach for modeling superconducting hybrid nanostructures with spin-orbit coupling, magnetic textures, and spin polarization. Representative applications in spin chains, two-dimensional and molecular magnets on a superconductor demonstrate its robust capabilities for both simple model systems and realistic ab initio molecular structures relevant to modern STM experiments.
Read the full abstract

All authors Oscar Moreno Segura, Dorye L. Esteras, Andrei Shumilin, Gonzalo Rivero-Carracedo, Riku Tuovinen, José J. Baldoví

Fluid Dynamics · Cross-listed from math.NA

Space-Time Parallel Direct Numerical Simulation of Rayleigh-Bénard Convection

Rayleigh-Bénard convection (RBC) is a model for turbulent flow and challenging to simulate by direct numerical simulation (DNS). Adequate resolution of all relevant length scales requires very fine grids in space and stability restrictions can dictate small step sizes. In this paper, we investigate the use of parallel spectral deferred correction (SDC) as a high-order time-stepping method to reduce wall-clock times of direct numerical simulation of RBC. Performance results on JUSUF (CPUs) and JUPITER (GPUs) with two different codes show parallel SDC can offer some concurrency along the temporal axis plus larger maximum stable time steps. Taken together, this allows for a reduction of wall-clock time required for simulations over one free-fall time by having to compute fewer steps and being able to mitigate the higher cost per time step through temporal parallelization.
Read the full abstract

All authors Thomas Saupe, Thibaut Lunet, Mathis Bode, Robert Speck, Sebastian Götschel, Daniel Ruprecht

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

Radiation limits on amino acid preservation and chiral measurement at Oxia Planum

Amino acid enantiomeric excess (EE) is a potential biosignature, but its measurement requires sufficient surviving material as well as preserved chirality. We derive a criterion that separates abundance-limited from racemization-limited preservation and apply it to the Mars Organic Molecule Analyzer (MOMA) on the ExoMars Rosalind Franklin rover. The model combines cosmic-ray and radiogenic doses with laboratory destruction kinetics and assumed ratios of interconversion to destruction rates. For an initially homochiral population of 1000~ppmw, the surviving EE is 58% or 25% when the total abundance reaches 1~ppmw for rate ratios of 0.04 or 0.10, respectively. This abundance criterion does not ensure that both enantiomers can be quantified. Applying provisional estimates of 1~ppmw per enantiomer and 10% EE, serine at 2~m requires initial concentrations of 340 and 960~ppmw after cosmic-ray exposures of $10^{7}$ and $10^{8}$~yr, respectively, under transferred fused-silica–perchlorate kinetics and a 1.57-MGy radiogenic dose. These concentrations exceed measured meteoritic serine abundances, but do not predict the initial Oxia inventory. For organics retained in a closed host since deposition 3.9~Ga ago, the modeled radiogenic dose is 1.6–5.6~MGy across the adopted compositions. In the lower-dose case, drilling to 2~m reduces the total dose by only a factor of 1.5 after $10^{7}$~yr. The boundaries depend on mineral kinetics and organic residence time; they distinguish preserved chirality from a measurable excess.
Read the full abstract

Optics

A Photonic Integrated Light Generation, Control, and Delivery Engine for Creating Cold Atoms

The creation of cold atoms underpins quantum sensing and precision metrology, computing, timekeeping, and fundamental science experiments. Miniaturizing the lasers and optical infrastructure, which occupy over 80% of the volume of these experiments, through photonic integration could improve robustness, reliability, and manufacturability while reducing cost and enabling portable space-based precision measurements and quantum tools for non-experts. Yet, there has been limited progress towards integration of the laser and optical infrastructure required to generate cold atoms onto a single chip. We demonstrate a photonic engine that generates, controls, and delivers precision light to create cold atoms in a rubidium 3-dimensional magneto-optical trap (3D-MOT). We integrate on a single chip two frequency-agile self-injection-locked lasers, laser frequency and intensity control, passive beam combining and splitting, signal monitoring, and high-efficiency beam delivery to a rubidium vapor cell. The integrated engine is shown to support multiple functions required for generating cold atoms including stable generation of low-noise cooling and repump light, frequency control for atom spectroscopy and stabilization, real-time laser monitoring, routing for laser offset locking, optical power modulation, and large-area three-beam delivery of 3 mW to the physics package. Fully electrically controlled, the engine enables sustained 3D-MOT operation, trapping more than 1 million atoms at an average temperature of 350 uK and supporting atom-number and temperature measurements. These results demonstrate the first single-chip integration of light generation, control, and delivery for cold-atom generation, offering a path toward smaller, more robust, and manufacturable cold-atom systems for applications including cold-atom clocks and sensors.
Read the full abstract

All authors Andrei Isichenko, Calisto Ruiz, Chuan Zhong, Nick Montifiore, Sam J. Harding, Meiting Song, Chenghang Liu, Nitin Indukuri, Karl D. Nelson, Daniel J. Blumenthal

Atmospheric and Oceanic Physics

Cumulus parameterization controls simulated rainfall more than microphysics, but not at heavy thresholds: a 3 * 2 factorial WRF study over the Western Himalaya

Rainfall forecasts over the Western Himalaya remain poor, and the relative contributions of the cumulus and microphysics parameterizations to that error are not well constrained.

Computational Physics · Cross-listed from cs.LG

Deep Mesoscopic Dynamics

Microscopic simulations resolve spatial patterns but are feasible only over small domains and short timescales, whereas bulk descriptions erase those patterns through spatial averaging. Mesoscopic fields preserve this spatial organisation, but their effective equations are known only for selected microscopic systems and scaling limits. We introduce Deep Mesoscopic Dynamics (DeepMesoD), which uses neural networks to learn an effective mesoscopic field equation from microscopic trajectories at a prescribed coarse-graining scale. For the nonconserved dissipative dynamics studied here, DeepMesoD learns the free-energy functional with a spatial interaction kernel and a state-dependent dissipation operator constrained by the fluctuation-dissipation relation. We first validate the construction in the Ising-Kac model, where the learned mesoscopic equation recovers the known hydrodynamic limit and rollouts remain stable far beyond the training horizon. We then treat two systems without closed-form mesoscopic descriptions. For a two-dimensional Heisenberg ferromagnet governed by atomistic spin dynamics, DeepMesoD learns a mesoscopic model whose predictions match atomistic references up to 1 ns, and a rollout with an unseen initial condition remains bounded up to 25 ns. For layered CrBr$_3$, the learned kernel captures the suppression and sign reversal of interlayer coupling at a stacking fault whose location and strength are not supplied as inputs. Controlled changes to the learned interaction show how fault number and placement shape multistep magnetic responses, complementing our magneto-photoluminescence measurements. Across the two atomistic systems, DeepMesoD achieves a three-to-four-order-of-magnitude reduction in wall-clock time compared with atomistic spin dynamics.
Read the full abstract

All authors Yue Zhao, Shiqi Wu, Yihang Ouyang, Sergey Grebenchuk, Magdalena Grzeszczyk, Xiaoxu Zhao, Kostya S. Novoselov, Maciej Koperski, Qianxiao Li

Optics

On-chip optical isolators

Optical isolation suppresses back-reflections that destabilize lasers, produce parasitic interference, and increase crosstalk in photonic systems. Although free-space and fibre-optic isolators are mature, no broadly adopted on-chip technology yet combines strong isolation, low insertion loss, broad bandwidth, compact footprint, and scalable fabrication. The difficulty is twofold: Lorentz reciprocity must be broken, and every relevant backward spatial, spectral, and polarization channel must be suppressed. This Review develops a common framework for assessing proposed solutions. We first distinguish pairwise nonreciprocity from intensity-asymmetric transmission and from true isolation using the constitutive relations and scattering matrix. We then compare three physical routes to an integrated nonreciprocal response: magneto-optical bias, externally driven space-time modulation, and nonlinear optical processes. Topological and non-Hermitian systems are treated as cross-cutting design frameworks that can route, protect, or enhance an underlying nonreciprocal response, but do not by themselves break reciprocity. Finally, we examine fundamental constraints and practical trade-offs among isolation, insertion loss, bandwidth, footprint, modal coverage, drive requirements, and fabrication compatibility. The resulting mechanism-based perspective connects reciprocity-breaking physics to the requirements of scalable photonic integrated circuits.
Read the full abstract

All authors Abanoub Mikhail, Nikolay Solodovchenko, Sen Xu, Fan Zhang, Maksim Tumashov, Mingzhao Song, Sergei Makarov, Andrey Bogdanov

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

Colossal Spin–Lattice–Charge Coupling in Noncollinear Ferrimagnetic NbMnX$_2$ ($X$ = S, Se, Te) Monolayers

Noncollinear ferrimagnets combine a net magnetic moment with internal spin degrees of freedom, creating opportunities to couple magnetic, mechanical, thermal, and electronic responses. Using first-principles calculations, symmetry analysis, atomistic spin simulations, and transport modeling, we identify stable NbMnX$_2$ ($X$ = S, Se, Te) monolayers with noncollinear ferrimagnetic ground states and Curie temperatures up to 313 K. Competing exchange, Dzyaloshinskii–Moriya interactions, and magnetic anisotropy produce field-responsive canted states that couple strongly to the lattice and charge carriers. Magnetization rotation changes the elastic modulus of NbMnTe$_2$ by up to 36%, while lattice entropy substantially enhances the near-room-temperature magnetocaloric response of NbMnS$_2$ and NbMnSe$_2$. In NbMnSe$_2$, rotation restores a mirror-protected crossing, continuously closes an approximately 0.13 eV gap, and produces a conductivity-based anisotropic magnetoresistance ratio up to 10,000% at 100 K. The combined results establish noncollinear NbMnX$_2$ monolayers as a multifunctional platform for two-dimensional spintronics, nanomechanics, and thermal control.
Read the full abstract

All authors Peng Wang, Nanshu Liu, Hongkuan Yuan, Xue Jiang, Jijun Zhao

Optics

Theory of deeply subwavelength electromagnetic structuring of semiconductors due to thermal-ionization instability

We develop a theory of ionization-assisted instability that governs the reflection of intense electromagnetic pulses from a semiconductor surface. Local overheating of the material increases the free-carrier density and conductivity, which in turn enhances the absorption of the incident radiation and further amplifies local heating and avalanche ionization. At high pump repetition rates, the overheating accumulates from pulse to pulse, leading to a gradual rise of the local material temperature and free carrier density - a positive feedback loop that drives thermal runaway. However, this effect competes with diffraction of the pump field on the emerging perturbation of the dielectric permittivity or on the surface relief, which reduces the local electric field amplitude in the overheated regions, thus providing negative feedback. We show that the competition between heat diffusion and diffraction near-field weakening creates conditions for the exponential growth of temperature perturbations with a certain wavenumber determined by the pulse repetition rate, pump field, thermal conductivity and other material parameters. The derived analytical expressions accurately describe the spatial period of deeply subwavelength structures observed in recent experiments on silicon structuring by intense terahertz radiation, and reproduce qualitatively dependencies of their period on the pump intensity, the incident wavelength and the number of pump pulses.
Read the full abstract

Optics

Practical Mid-infrared Dual-comb Spectroscopy with Absolute Frequency Referencing Based on a Single-cavity Dual-comb Fiber Laser

We demonstrate a practical mid-infrared (MIR) dual-comb spectroscopy (DCS) method based on a bidirectional dual-comb fiber laser. By exploiting the inherently high relative stability and mutual coherence of the dual-comb laser, we stabilize the MIR dual-comb modes with sufficient precision using only piezoelectric control of the repetition frequencies. This simple configuration enables correction of residual timing jitter and phase noise in interferograms while preserving the absolute optical frequency scale defined by the experimentally obtained comb parameters, without requiring an external continuous-wave laser reference or tight phase-locking. Importantly, the proposed method does not merely produce apparently clean interferograms through signal processing, but enables coherent averaging while retaining both the physically meaningful interferometric waveforms and well-defined absolute frequency information. The proposed system is experimentally validated through MIR spectroscopy of N$_2$O around 3.9 $\mu$m, achieving a signal-to-noise ratio of approximately 690 through coherent averaging of 10,000 interferograms over 65 s, along with the acquisition of an absorption spectrum consistent with the reference data. These results demonstrate a simple and practical approach to absolute-frequency-referenced MIR-DCS.
Read the full abstract

Optics · Cross-listed from quant-ph

Boson Sampling Using Fock States from On-Demand Single Photon Source Arrays: Quantum Complexity & Impact of Dephasing

Boson sampling with array of on-demand single photons (discreet variable Fock states) as resource states provides distinct advantages in complexity scaling compared to single photon resource states mimicked by the states (continuous variable) generated in nonlinear photonics. Recent progress in realization of regular ordered array of mesa top single quantum dots (MTSQDs) enabling on-demand generation of many photon resource states with high single photon purity and indistinguishability [1,2] has opened the possibility of implementing chip-scale Boson sampling with single photons. An important question towards scale up of such circuits is how the dephasing and decay timescale of individual on-demand single photons affect the final probability distribution of the N-photon entangled state. In this work we report a framework that provides direct link between the physical sources of decoherence to the outcome of multiphoton correlation measurements. We explicitly incorporate the symmetrization of the N-photon states while accounting for the shape of the individual wave packets and dephasing. Using this framework, we show predictive estimates of the probability of generation of a 4-photon NOON state using an optical circuit initialized with 4 identical single photon wave packets, mimicking the emitted photons from MTSQD single photon sources reported in [2]. Specifically, we show how the probability of generating various 4-photon entangled states depend on the two-photon interference visibility which is determined by the pure dephasing timescale in relation to the radiative decay timescale. The results provide quantitative guidelines to measure 4 photon and higher number of entangled photon states generated in an optical circuit with on demand single photon sources and provide a link between physical processes such as temperature and the output state of the quantum optical circuits.
Read the full abstract

Physics and Society

Replacement distance controls growth at long wavelengths in a coevolving spatial network

When an agent leaves a spatial network, its replacement removes old contacts and may transport population to a new location. The frequency of replacement alone therefore does not determine how turnover changes spatial organization. We study diffusing agents joined by persistent attractive contacts, with the total pull on each agent bounded. At fixed population, successors appear at the departure position, within a surrounding disk, or anywhere in the system. We ask whether small differences between dense and dilute regions grow or fade. The theory separates an exactly calculable diffusion and replacement contribution from an approximate attraction contribution obtained by averaging the contact network. Paired simulations find similar attraction contributions under the three rules at all nine sampled wavelengths. Changing only the replacement location reverses the broadest tested pattern: it grows with replacement at the departure position but fades under both displaced rules. An independent radius scan traces this change from growth to decay. Under displaced replacement, intermediate wavelengths grow between sampled long and short scales with decay. Rates measured in the earliest finite window are compatible with the prediction at long and intermediate wavelengths; discrepancies are already present at short wavelengths. Replacement distance thus helps determine which spatial density contrasts are initially amplified.
Read the full abstract

Atmospheric and Oceanic Physics

Extreme events and statistics in an AI emulator of a global storm-resolving model

We present HiRO-ACEv2 (High Resolution Output Ai2 Climate Emulator), a two-stage framework for AI emulation of the 3~km resolution X-SHiELD global storm-resolving model (GSRM), and highlight its ability to generate large ensembles of data for studying the characteristics of tropical cyclones and extreme events in high-resolution simulations at an inference computational cost over 300 times lower than that of traditional numerical GSRMs. HiRO-ACEv2 expands upon the precipitation-only HiRO-ACEv1 model to generate coherent precipitation, 10~m winds, and sea level pressure outputs, enabling a broader range of events to be simulated at high-resolution. Fine-resolution predictions are generated by running an autoregressive emulator of coarsened X-SHiELD at 100 km, whose multivariate outputs are regionally downscaled to 3~km resolution via a separately trained diffusion model. Particular importance is placed on capturing the statistics of extreme events during the diffusion downscaling step. Our ability to match the tails of distributions and generate realistic-looking extreme events in both the `perfect prediction' case as well as in free-running emulator rollouts is improved via the use of multiple denoisers at inference time and several updates to the data sampling, noise sampling, and loss weighting during training. HiRO-ACEv2 closely matches the high-resolution reference distributions of precipitation and near-surface winds out to their 99.9999th percentile values, and broadly tracks the reference distribution of tropical cyclone intensities.
Read the full abstract

All authors Anna Kwa, W. Andre Perkins, Troy Arcomano, Noah D. Brenowitz, Christopher S. Bretherton, Spencer K. Clark, Lucas M. Harris, Jeremy McGibbon, Michael S. Pritchard, Oliver Watt-Meyer

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

A Parameter-Free Reference Hamiltonian for Tight-Binding Vanadium Oxides

Density-functional tight binding has yet to yield a transferable description of vanadium oxides; for the standard second-order construction and the tested cells, we demonstrate that two of its approximations deviate by electron volts and offset each other in the band-index centroid separation between the V 3d and O 2p bands: the neutral two-centre reference omits electron-volt environmental terms, whereas the Mulliken charge model assigns the compact V 3d shell less than half of the self-consistent potential shift observed in PBE. Instead, we build a tight-binding model where the reference Hamiltonian derives from superposed neutral atoms, the full PBE exchange-correlation potential and frozen-core orthogonality, without adjustable terms, and where the orbital-resolved charge response employs three screening factors obtained from PBE linear response and an unscreened Coulomb range. Relative to the primary PAW reference, the plane-wave Harris-Foulkes energy differences of VO2 are matched by the reference energy at fixed geometry and the stated numerical settings, with signed deviations of +2.6 meV/atom for M1 relative to rutile (an energy difference of 92 meV/atom) and +6.2 and -5.2 meV/atom under 3% volume compression and expansion; the charge response reproduces 14 of 22 fitted and 5 of 8 held-out matched PBE responses within one tolerance, and 28 of the 30 within twice the tolerance, and in rutile the fitted first-radial V 3d shift is -4.39 eV, 0.47 eV more negative than PBE's -3.92 eV, whereas the tested DFTB2 charge model gives -1.53 eV. The convexity of the spinless model at fixed geometry prevents it from spontaneously differentiating symmetry-equivalent sites under the symmetries respected by its implemented operators and sampling, so within this model such differentiation in the Magnéli phases must arise from electron-lattice coupling or an explicit spin or correlation term.
Read the full abstract

Optics · Cross-listed from eess.IV

Background2Background: a self-supervised background removal framework for high dimensional microscopy

High-dimensional microscopy captures features across space, time, and spectrum, but complex backgrounds obscure signals. We present Background2Background (B2B), a self-supervised background removal framework leveraging the backgrounds larger correlation relative to the signal. Local pixel shuffling disrupts signal correlations while preserving background structure, transforming signals into effective noise for Noise2Noise-based background removal. Across fluorescence, stimulated Raman scattering and mid-infrared photothermal imaging, B2B achieves over 40-fold background suppression and reveals hidden signatures.
Read the full abstract

All authors Guangrui Ding, Haonan Lin, Rong Tang, Bethany Weinberg, Jianpeng Ao, Qing Xia, Guorong Hu, Dingcheng Sun, Jiaze Yin, Lei Tian, Ji-Xin Cheng

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

A grand potential for molten-salt electrochemical phase diagrams from first principles

The clean-energy transition is driving unprecedented demand for metals, whose extraction is among the largest industrial sources of carbon emissions. Molten-salt electrolysis offers a path toward carbon-free production of aluminum, iron, titanium, and critical metals. However, corrosive operating conditions and complex, multicomponent feedstocks make experimental development costly and challenging. Here, we introduce a first-principles framework for constructing electrochemical phase diagrams of multicomponent systems to accelerate electrode and electrolyte design and guide the selection of compatible materials and operating conditions for electrochemical cells. We derive a grand potential that generates these diagrams at scale, and supply free energies using first-principles data obtained through a fine-tuned foundation machine-learned interatomic potential. For silicon in molten CaCl$_2$, the construction not only recovers the experimental diagram when using the first-principles framework on a limited set of structures, but also predicts further domains when extended to a larger set of crystal structures. Finally, we demonstrate the predictive power of the approach by screening ten alkali and alkaline-earth halide electrolytes for aluminum electrowinning. Three thermodynamic criteria identify NaF, where cryolite emerges as the stable aluminum-bearing phase, recovering the chemistry that Hall and Héroult established empirically in 1886. The framework provides a first-principles basis for selecting compatible electrodes, electrolytes, and operating conditions toward low-carbon metal production from complex source materials.
Read the full abstract

Atomic Physics · Cross-listed from quant-ph

Fast charging of quantum batteries assisted by nonlinear couplings

Charger-mediated quantum batteries offer promising opportunities for efficient energy transfer and storage at the quantum level, with potential advantages beyond their classical counterparts. Here, we theoretically demonstrate that nonlinear charger-battery interactions can substantially accelerate the charging of quantum batteries. In the weak-driving regime, sufficiently strong nonlinear coupling can dominate over the effect of dephasing and significantly reduce the charging time. In the strong- and optimal-driving regimes, where the optimal driving strength maximizes the steady-state ergotropy in the linear-coupling regime, nonlinear coupling gives rise to a new optimum in the charging time beyond that of the linear-coupling regime, accompanied by nearly critically damped dynamics. Remarkably, even in the absence of external driving, nonlinear coupling alone can substantially reduce the charging time, providing an alternative route toward fast battery charging. Our results establish nonlinear charger-battery coupling as an effective resource for accelerating quantum-battery charging and broaden the scope of interaction-enabled quantum advantages in energy transfer and storage.
Read the full abstract

All authors Y.-J. Shen, Y.-C. Shen, C.-H. Chu, H. H. Jen

Fluid Dynamics

Gas radiation and the gravity limit of flicker in a co-flowing hydrogen jet diffusion flame

Although weak at normal gravity, gas radiation sets the gravity at which a co-flowing hydrogen jet diffusion flame stops flickering. Flamelet-based simulations with unity Lewis numbers, which reproduce the flicker frequency measured by femtosecond coherent anti-Stokes Raman scattering within 4%, are combined with local stability analysis. At one tenth of normal gravity, the adiabatic flame flickers, the radiating flame does not, and switching the emission off restores flicker. Radiation acts through two mechanisms that do not weaken at reduced gravity. The hot outer layer loses buoyancy in proportion to the gas residence time, and a disturbance that displaces the flame also displaces its emission, which damps the disturbance on the lean flank (displaced-emission damping). Freezing the radiative source, which removes only its response to disturbances, slows the downstream recession of the oscillation more than tenfold. Both mechanisms slow the upstream spreading of disturbances. This spreading, competing with the co-flow, sets the onset location and the critical co-flow. Without radiation, the spreading speed scales with the square root of gravity; its radiative reduction does not. At reduced gravity, the flicker limit therefore becomes a straight line relating the square root of the co-flow Froude number to a radiation number, the ratio of buoyant to radiative time scales. The calibrated line places the critical gravity for the experimental co-flow at 0.13 times normal gravity, 2.6 times the adiabatic value. It correctly predicted whether two further flames keep flickering. Beyond the limit, transient wave packets can mimic flicker throughout a short drop-tower test.
Read the full abstract

Plasma Physics

Effects of Coulomb Collisions on the Structures of Fast Magnetosonic Shocks

We present fully kinetic particle-in-cell simulations to demonstrate qualitative effects of binary Coulomb collisions on fast magnetosonic shocks. We find that with sufficiently high collisional frequencies, shock rippling and reformation can be inhibited, leading to stationary, laminar shock fronts. We observe that ion-ion collisions disrupt the motion of particles reflected back upstream that is the key for cyclical reformation. Meanwhile, the electron-ion collisions provide a Spitzer resistivity that can dissipate energy from the instability that causes shock rippling. This result helps us understand the role of collisions in laboratory shocks and possibly the outer heliosphere.
Read the full abstract

All authors Michael B. Farnell III, Fan Guo, Yi-Hsin Liu, Xiaocan Li

Optics

Inverse design of poles and zeros: tailoring phase control in metasurfaces

We develop a physically interpretable inverse-design framework for generalized Huygens metasurfaces, using prescribed scattering poles and zeros to engineer resonance interference and transmission-phase control. Compared with direct optimization of transmission or reflection coefficients, this approach offers a physically transparent route to rationally controlling resonance properties, including linewidth and group delay, while reducing computational cost. The resulting designs demonstrate how near-unity transmission and $2\pi$ phase excursion can arise from resonance interference. These signatures of the Kerker effect extend beyond the conventional alignment of electric and magnetic dipole resonances to other mechanisms, including guided-mode resonances and quasi-bound states in the continuum. We further optimize the Kerker condition across a continuum of refractive-index states in a reconfigurable phase change material metasurface, enabling high-transmission $2\pi$ phase tuning. Our results establish a versatile and practical approach to metasurface phase control through the targeted design of scattering poles and zeros.
Read the full abstract

All authors Junhe Zhang, Seoyoung Kim, Hyun Jung Kim, Juejun Hu

Atmospheric and Oceanic Physics

The hypervariational method: a general, automated technique based on innovations to tune the parameters of the Hybrid covariance matrix

The Hybrid forecast covariance model is an operationally successful way to build a background error covariance matrix. It is constructed as a weighted sum, with scalar weights, of the climatological covariance matrix and a localized ensemble covariance matrix, with a localization matrix typically depending on a localization length. For best performance, the Hybrid needs to be optimized by tuning those weights and, possibly, other parameters, such as the localization length. Here, a systematic and automated method for this tuning operation is introduced. This technique is referred to as hypervariational method, and it is formulated in terms of a cost function, built using a time series of innovations as an essential ingredient. This cost function, which depends on the parameters, is minimized when they have optimal values. Equations for its minimization, based on Newton's method, are derived. Besides the aforementioned, commonly used Hybrid, the hypervariational method is applicable to any known Hybrid (and to any set of parameters on which it depends analytically), including recently developed Hybrids, not necessarily constructed as weighted sums. Here, the hypervariational method is applied to tune the ensemble covariance matrix weight in the seminal-formulation Hybrid (which does not feature localization). This is done in a synthetic experimental setting where a theoretically-derived expression for the optimal value for that weight is known. Moreover, in that experimental setting, a further method, based on a linear fit procedure, to determine that optimal value is possible. Both the hypervariational method and that linear fit procedure are shown to produce an optimal weight value matching (up to a tolerance of $0.1$) the optimal value known from the theory. Furthermore, Newton's minimization in the hypervariational method is shown to converge immediately (at the first recursive step).
Read the full abstract

Accelerator Physics

A modern particle-in-cell code for beam-beam simulations of future colliders

A key challenge in the design of next-generation colliders is predicting the effects of beam-beam interactions and using these predictions to optimize machine performance. Beam-beam interactions are characterized by strong focusing forces, radiative effects, and the generation of background particles. While analytical formulas provide insight into the dynamics of beam-beam interactions, these nonlinear, coupled effects are most accurately described with numerical methods. In this paper, we introduce the open-source particle-in-cell code WarpX as a modern tool for modeling the interaction point in future lepton colliders. WarpX provides high-performance execution on distributed GPU architectures within an open-source and actively maintained code base. We present convergence scans and benchmarks against the established beam-beam codes GUINEA-PIG (v1.2.2) and CAIN (v2.42) for luminosity spectra, photon emission, and secondary-particle production under a wide range of beam parameters corresponding to different collider concepts. Our results show that WarpX (v26.09) reproduces the benchmark beam-beam observables while enabling GPU-portable, extensible simulations. These results establish WarpX as a scalable platform for high-fidelity beam-beam modeling in future collider studies.
Read the full abstract

All authors Arianna Formenti, Remi Lehe, Axel Huebl, Stepan S. Bulanov, Jens Osterhoff, Andrew Myers, Weiqun Zhang, Spencer Gessner, Peter Kicsiny, Anna Kinderman, Bao Nguyen, Jinseo Park, Jean-Luc Vay

Physics and Society

Hurricane recovery gaps lie in reconnection, not local restart

Climate-intensified hurricanes disrupt urban mobility, yet whether recovery is equitable across the rural-urban gradient is unclear. Analyzing device-based mobility for Hurricanes Helene and Milton, we find local activity restarted within approximately 4.5 days in both storms, whereas reconnection to the outside world took 5 days in coastal-urban Florida but 11 in inland-rural Appalachia. The gap lies in external reconnection, not local restart, implying planners should target rural isolation.
Read the full abstract

All authors Angelina Liu, Mengqiu Cao, Sen Pei, Qing Yao

Optics

Non-contact lens metrology from optical flow and differentiable ray tracing

Accurate lens surface metrology is essential to determine the manufactured geometry of refractive optical elements and their deviation from the specified design. We demonstrate a transmission metrology method that reconstructs refractive lens surfaces from the displacement of a wavelet-noise-patterned mask. Lens-induced distortion is represented by a dense displacement field. A differentiable ray tracer fits the lens geometry using a calibrated near-telecentric camera model with an effective finite ray origin. We validate the method on commercial lenses. The recovered geometries agree with nominal specifications, and the plano-aspheric reconstruction is compared with chromatic confocal profilometry.
Read the full abstract

Physics and Society

Collective Benefits and Efficiency in Coupled Echo State Networks

Simple aggregation of independent responses ("wisdom of crowds") generally improves collective performance through statistical facilitation, but interaction among members can sometimes hinder group outcomes. However, under certain circumstances, coordination studies show that interacting groups can outperform their individual members and nominal groups, and interaction can also be more efficient than independence. We present a reservoir computing model to investigate how interaction can enhance joint performance in terms of accuracy and efficiency. Echo state networks were trained independently to identify Japanese speakers from vowel recordings and tested jointly with and without coupled feedback. Model pairs reproduced the dyadic advantage reported by Bahrami et al. (2010) with interaction via pooled-feedback, but not when independent responses were aggregated together. Model dynamics revealed that interaction was beneficial whereas aggregated outputs were prone to diverging from runaway feedback over time. In terms of efficiency, larger groups of smaller interacting models outperformed smaller groups of larger models while equating for total amount of group computation. Our results offer a formal framework for explaining when and why interaction enhances collective intelligence.
Read the full abstract

All authors Polyphony Bruna, Rick Dale, Michael J. Spivey, Christopher T. Kello

Medical Physics

Photoacoustic imaging of human brain function through intact adult skulls

Blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has revolutionized noninvasive imaging of human brain function, but it indirectly senses deoxyhemoglobin through water-proton relaxation and relies on bulky, costly, noisy, enclosed, superconducting-magnet scanners. In contrast, photoacoustic computed tomography (PACT) directly probes hemoglobin absorption using a more compact, lower-cost, quieter, open, magnet-free system. However, functional PACT of the adult human brain has so far been demonstrated only after surgical skull removal, because the skull strongly attenuates optical and acoustic signals. Here, we demonstrate non-invasive transcranial functional PACT (T-fPACT) through the intact adult skull. A densely packed, high-channel-count (3,072-element) hemispherical ultrasonic array provides near-$2\pi$-steradian simultaneous acquisition, enabling high sensitivity and single-shot three-dimensional imaging at 20 Hz that eliminates motion blurring. Participant-specific stabilization further enables measurement of small task-evoked hemodynamic responses over the entire task duration. In healthy adults, T-fPACT localized responses to expected sensorimotor regions, distinguished task-specific motor representations, and yielded consistent localizations across repeated task runs and visits. Comparison with high-field (7 T) fMRI showed spatially corresponding functional regions, while dual-wavelength imaging revealed distinct hemoglobin-weighted responses. These findings overcome the long-standing skull barrier in human brain PACT and establish non-invasive, rapid volumetric photoacoustic imaging of brain function.
Read the full abstract

All authors Yilin Luo, Rui Cao, Manxiu Cui, Karteekeya Sastry, Kay Jann, Yousuf Aborahama, Geng Ku, Jinhua Xu, Xiaofei Luo, Shuai Na, Samuel Davis, Konstantin Maslov, Yang Zhang, Danny J. Wang, Lihong V. Wang

Monday 12 October

Fluid Dynamics

Weber-number scaling and the distance dependence of drop size in full-cone pressure-swirl sprays

A pressure-swirl nozzle converts the feed pressure into a swirling conical sheet that breaks into ligaments and then into drops. The size of the drops governs the deposition and drift of agricultural sprays, the capture of airborne dust and the combustion efficiency of fuel injectors. The standard tool for measuring the drop size is laser diffraction, which returns a spatial, concentration-weighted distribution at the beam station. However, the dependence of the measured drop size on the nozzle-to-beam distance has not been quantified alongside a Weber-number scaling, so prefactors reported at different distances cannot be compared on an equal footing. To this end, we measure the drop-size distribution of two full-cone pressure-swirl nozzles by laser diffraction at a range of distances, pressures and orientations, together with high-speed imaging of the sheet breakup. We find that the volume-median diameter $D_{v50}$ decreases with the measurement distance $L$ as $D_{v50} \sim L^{-0.18}$ beyond the drop-formation region. A joint fit over both nozzles gives $D_{v50}/d_0 \sim \mathrm{We}^{-0.31}\,(L/d_0)^{-0.21}$, where $d_0$ is the orifice diameter and $\mathrm{We}$ the Weber number, which separates the pressure and distance effects. At a fixed distance the two nozzles collapse onto the $\mathrm{We}^{-1/3}$ scaling of Kooij et al. (2018), and the measured distance dependence may account for part of the offset between our near-field prefactor and the far-field value. Thus, a drop size obtained by laser diffraction is meaningful only when the measurement distance is reported with it, and a spray should be characterized at the actual distance between the nozzle and its target.
Read the full abstract

All authors Bilal Ben Jemaa, Vincent He, Marko Trifunović, Daniel Duke, Rishabh V. More

Chemical Physics

TorchSpin: Differentiable Electron Paramagnetic Resonance Simulation in PyTorch

Interpreting electron paramagnetic resonance (EPR) spectra commonly requires repeated forward simulation and optimization of spin Hamiltonian parameters. Established software provides mature simulation and fitting methods, but these calculations are not generally available as differentiable components within machine-learning workflows. We present TorchSpin, an open-source PyTorch framework that ports the organization and numerical conventions of EasySpin into a tensor-based computational environment. Representative calculations spanning continuous-wave EPR, ENDOR, pulse EPR, spin dynamics, trajectory simulation, and magnetometry closely reproduced corresponding EasySpin results. In matched benchmarks, TorchSpin CPU execution was faster than EasySpin for six of nine workloads although performance depended strongly on the calculation. We therefore focused on the capabilities enabled by automatic differentiation in PyTorch. In a Cu(II) powder-spectrum demonstration, feedback propagated through TorchSpin improved the quality of a supervised convolutional neural network (CNN) without using parameter labels during refinement. Models trained from randomly initialized CNNs using only noisy spectra and TorchSpin reconstruction loss achieved median similarities of $\geq$ 0.99. CNNs trained on TorchSpin-generated libraries also transferred without adaptation to independently generated EasySpin spectra with performance comparable to EasySpin-trained controls. Once trained, the CNN produced initial parameter estimates in approximately 1 ms per spectrum while a batch of 512 was processed in approximately 1.3 ms with GPU acceleration. TorchSpin therefore provides both a tested Python implementation of established EPR calculations and a differentiable physical layer for rapid learned inference, spectral refinement, and parameter-label-free training.
Read the full abstract

Physics and Society

Popularity Trajectories and Mixing Patterns of Urban Venues in Two European Cities

Urban venues such as restaurants, shops, bars and parks host the daily life of a city, and their collective rise and fall signals how vibrant a neighbourhood is. The usual assumption is that a venue's fortunes follow its neighbourhood's, so that how wealthy an area is and how many people live there will decide a venue's performance. Yet whether there is a common performance trajectory shared by venues over time is underexplored. In this study, we use 5.5 years of Foursquare data for about 10^4 venues in London and Barcelona. We found six distinguishable shapes recur in both cities, that are (1) Fall, (2) Rise, (3) Rise & Fall, (4) Steady, (5) Ephemeral (an early hit that loses popularity fast) and (6) Underdogs (a slow starter that gains popularity late). However, no recovery shape appears, that is neither city yields a U-shaped fall-and-rise group. Additionally, we found that deprivation and population density are not indicative which shape a venue follows. We also discovered that how early a shape can be recognized depends on when its distinctive segment occurs, not a shape's complexity. Proximity between venues, however, does predict how well-performing a venue is. We uncovered a mixing pattern where venues of the same shape tend to locate together more often than chance within roughly 200 to 600 metres. Finally, our qualitative reading of two contrasting districts, one rising in London and one falling in Barcelona, show that comparable socioeconomic conditions but divergent gentrification histories can generate vastly different popularity trajectories. All in all, our results altogether point at that urban vibrancy emerges at the block- and street- level, and social and historical contexts matter in terms of the characteristics of such patterns.
Read the full abstract

All authors Wang Ngai Yeung, Lakmal Meegahapola, Daniele Quercia, Riccardo Di Clemente

Geophysics

A Multi-frame Analysis of Flow Reversal in Tidal Bore

A tidal bore is a nonlinear step like surge propagating on the top of a river counter-current. The scientific literature reports that the passage of the tidal front may induce local flow reversals and, in some cases, a global reversal of the river current. Here, combining Particle Image Velocimetry (PIV) measurements and 2D numerical simulations based on a Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) method, we show that the flow classification of tidal bores must explicitly take into account the frame of reference. When the river flow does not reverse after the bore passing, this is a consequence of the perspective imposed by the background current (itself resulting from the previous tide and the natural river discharge toward the sea). Depending on the relative magnitudes of the river current and the velocity field induced by the bore, one may observe a global reversal, a local reversal (or both), or no reversal at all. When the tide is spatially inhomogeneous, with a tidal range that varies along the estuary rather than a simple step in water depth, transitions between bore types may occur during the bore passing.
Read the full abstract

All authors Assia Belemkadem, Johan Fourdrinoy, Anthony Beaudoin, Germain Rousseaux