Papers

Recently it has been proposed that comparing the net-baryon ($B$) stopping with net-electric charge ($Q$) stopping can help studies of the baryon stopping mechanism in nuclear collisions. Here we find the $B/Q\times Z/A$ ratio to be very sensitive to the difference between $s$ and $\bar s$ quark rapidity distributions. For mid-rapidity of isobar collisions at 200A GeV, a multi-phase transport (AMPT) model gives slightly more $\bar s$ than $s$, which leads to $B/Q\times Z/A<1$, while the model without the $s-\bar s$ asymmetry gives $B/Q\times Z/A \geq 1$. Comparing Ru+Ru and Zr+Zr isobar collisions, the AMPT (and UrQMD) model gives $B/\Delta Q\times \Delta Z/A<1$ at mid-rapidity at all centralities, which strongly contradicts the recent STAR data. We also find that the $B/\Delta Q\times \Delta Z/A$ ratio is very sensitive to the net-light quark ($u,d$) stoppings, but it is less sensitive to the $s-\bar s$ asymmetry than the $B/Q\times Z/A$ ratio by a factor of 3. These results are expected to help us resolve the stopping puzzle and better understand the baryon stopping mechanism in the future.

We develop a comprehensive Wilson-line formulation of quark Reggeization in QCD. Starting from a next-to-eikonal operator built from a semi-infinite Wilson line and a background-quark insertion, we identify an interpolating operator for the Reggeized quark and derive its nonlinear rapidity evolution using the background-field method. In the dilute regime, its positive-signature component exhibits Regge-pole evolution governed by the quark Regge trajectory, whereas the negative-signature sector displays mixing between quark and gluon degrees of freedom, in accordance with its known Regge-cut structure. In the planar limit, this mixing is suppressed, and Reggeization emerges without an explicit signature projection, recovering signature degeneracy. We illustrate the universality of the construction by extracting the same Reggeized-quark operator from a more general Wilson-line operator describing a gluon-to-quark transition. Furthermore, we extend the formalism to massive quarks, deriving a coordinate-space evolution kernel whose Fourier transform reproduces the massive-quark Regge trajectory. Finally, from the leading operator-mixing structure of the evolution equation and signature arguments, we show that the Reggeized-quark interpolating operator is expected to remain an eigenstate of the rapidity evolution up to next-to-leading logarithmic accuracy.

We prove that, for all $\beta>0$, the Dobrushin-Lanford-Ruelle (DLR) equations for the unit-intensity $\mathsf{Sine}_\beta$ point process have a unique stationary solution under a particular finite electric energy condition. Uniqueness fails if this condition is dropped. This answers a question of Dereudre-Hardy-Leblé-Maïda and gives a canonical statistical physics characterisation of $\mathsf{Sine}_\beta$. The main technical ingredient is a relative entropy estimate which allows us to compare the conditional distribution in a finite box of a DLR solution to a corresponding circular $\beta$ ensemble.

The mean squared displacement (MSD) is a cornerstone in the analysis of diffusion processes in complex media. When the system is heterogeneous and, in particular, when single-particle trajectories are short, it is essential to extract maximal information from each measured trajectory. This is typically done by time-averaging squared increments and examining the scaling of the time-averaged MSD in log-log space. However, classical regression methods perform poorly in this setting because time-averaging introduces correlations aggravated by those inherent to anomalous diffusion. We tackle these limitations by applying a generalized least squares framework, which substantially reduces variance and bias in diffusion parameter estimates, especially for short (around 100 points) and ultra-short (around 10 points) trajectories. The method is fully automated and requires no supervision. Furthermore, it enables prediction of estimation error probability density, which is asymptotically Gaussian, for both classical and enhanced approaches. Leveraging this prediction, we introduce a specialized deconvolution algorithm that reconstructs the underlying particle ensemble structure from experimental data.

We propose a novel framework that simultaneously addresses the origin of Dirac neutrino masses and the nature of self-interacting dark matter (SIDM). The model is based on an $A_{4}$ modular symmetry to ensure the Diracness of neutrinos as well as the stability of the DM. The neutrino sector realizes a Dirac Inverse Seesaw mechanism where the smallness of the neutrino mass is governed by the vacuum expectation value (VEV) of a singlet scalar $\phi$. This same scalar couples to a vector-like fermion DM candidate, inducing a tiny Majorana mass splitting that renders the DM pseudo-Dirac and inelastic. Crucially, the scalar also acts as a light mediator for DM self-interactions, potentially solving the small-scale structure problems of Cold DM. In light of the recent 248 keV nuclear-recoil event, LZ230616, observed by LUX-ZEPLIN (LZ) DM direct detection experiment, we demonstrate that our inelastic SIDM parameter space naturally accommodates this signal via endothermic scattering kinematics. Furthermore, the spontaneous breaking of the dark parity required for this inelasticity produces a network of cosmological domain walls. We show that the explicit symmetry breaking needed to safely annihilate these walls generates a stochastic gravitational wave background. The non-holomorphic modular symmetry reduces the free parameters, correlating neutrino observables, addressing DM phenomenology, $\Delta N_{\rm eff}$, and gravitational-wave signatures.

The recent convergence of 6G wireless systems and Tiny Machine Learning (TinyML) has driven the need for on-device intelligence in edge networks, where ultra-low latency, stringent energy budgets, and tight compute constraints demand novel architectures. Lightweight deep models efficiently extract local patterns but fail to capture global dependencies, while attention mechanisms do so at the expense of energy and computational cost. To bridge this gap, we introduce Quantumer, a hybrid TinyML–quantum framework that integrates multi-scale dilated convolutions and scaled dot-product attention within a lightweight transformer architecture, employing a two-stage transfer learning pipeline from Quantum Pre-Training (Quantumer-Q) to Classical Fine-Tuning (Quantumer-C). We also present QuantiblentLayer, a four-qubit variational circuit that maps compact traffic representations into measurement-based Hilbert-space features using trainable rotations and cyclic entangling operations. The circuit is used only during offline pre-training as a nonlinear embedding teacher and is removed before Quantumer-C deployment, leaving a fully classical inference model without runtime quantum execution. By transferring these quantum-assisted embeddings into an energy-efficient, lightweight transformer, Quantumer achieves strong detection performance with minimal compute and memory overhead on resource-constrained edge devices. The intrusion detection system (IDS) is used as a case study and evaluated on the Edge-IIoTset, TON IoT, and WUSTL-IIoT-2021 datasets. Quantumer-Q achieves competitive compact-model performance with 105.86K parameters, 0.4038 MB memory usage, 0.5525 MB model size, and 5.5646 MFLOPs; the INT8 Raspberry Pi 4 deployment obtains 16.8413 ms latency with a 0.6493 MB footprint. These results support training-time quantum-assisted representation learning for compact edge-deployable IDS.

We investigate secure communication over a classical multiple-access wiretap channel (MAC-WTC), specifically exploring the benefit of shared entanglement between the transmitters. Under the strict semantic security criterion, we derive an achievable rate region and a regularized expression for the secrecy capacity. We further establish a single-letter upper bound on the secure capacity of the MAC-WTC with entangled transmitters. Our results demonstrate that entanglement strictly enlarges the secrecy capacity of MAC-WTC compared to sharing only classical correlated randomness, a finding we illustrate using a pseudo-telepathy game example. Finally, we establish new strong soft-covering lemmas for the output statistics of multipleaccess channels (MACs) with entangled transmitters. Our results generalize existing results for non-entangled systems.

A multimessenger method is developed for scheduling optical spectroscopy of ultracompact double white dwarf binaries from Laser Interferometer Space Antenna (LISA) data. The orbital parameters inferred from the gravitational wave (GW) signal are used to predict the argument of latitude $u(t)$, which sets the phase dependence of the radial velocity. With $u(t)$ known in advance, spectra can be taken at selected orbital phases. The GW and spectroscopic descriptions are linked by $u(t)=\Phi_{\rm GW}(t)/2+\delta_{e,\omega}(t)+\mathcal{O}(e^2)$, where $\Phi_{\rm GW}$ is the GW phase, $e$ the eccentricity, and $\delta_{e,\omega}$ the equation of the centre to $\mathcal{O}(e)$, which vanishes for a circular orbit. Radial velocities are regressed on $\cos u$; for double-lined binaries, the two slopes determine the semi-amplitudes $K_1$ and $K_2$. With the inclination and orbital period measured by LISA, both component masses follow from the semi-amplitudes through Kepler's third law, and the resulting chirp mass with the GW amplitude gives the source distance. We derive waveforms to first order in the Laplace–Lagrange parameters; the circular quadrupole waveform is recovered as $e\to0$. Using simulated LISA data for one circular and one eccentric system, we forecast how long $u(t)$ remains accurate enough for phase-resolved spectroscopy. Masses derived from the velocity amplitudes are free of the biases that tides, mass transfer, or acceleration along the line of sight introduce when $\dot f_{\rm GW}$ alone is used. For a four-year LISA observation, scheduling at quadrature remains viable for roughly a decade after the observation ends. Four years after the observation midpoint, the ephemeris predicts the orbital phase to $0.3$~s for the eccentric system and $0.9$~s for the circular one, compared with $1$–$30$~s for optical ephemerides propagated over the same baseline.

The interaction of two qubits with a common bosonic reservoir is encoded by a matrix-valued spectral density $\mathbf{J}(\omega)$, whose diagonal entries describe the local spectra, while the off-diagonal entries encode cross-correlations. Even when \(\mathbf{J}(\omega)\) has rank one, a frequency-independent dark channel need not exist because the family \(\{\mathbf{J}(\omega)\}_{\omega}\) may have a trivial common kernel. We term the frequency-dependent rotation of the bright and dark directions spectral twisting and quantify it through the Fubini–Study speed $\tau(\omega)$ of the bright spectral projector. We analyze how twisting modifies two-qubit dynamics and quantify the loss of dark-state protection through the leakage \(P_{\mathrm{leak}}(t)\). Comparisons with untwisted asymmetric reservoirs and rotating-wave dynamics, together with detuned and finite-temperature calculations, distinguish spectral twisting from coupling asymmetry, counter-rotating processes, and thermal absorption. For resonant qubits tuned to the crossing of the two local spectra $\omega_\times$, the singlet is locally dark at the transition frequency but couples to off-resonant components whose bright directions are rotated. In the weak-twisting regime, the fixed-time leakage scales as $P_{\mathrm{leak}}(t)\propto[\omega_\times\tau(\omega_\times)]^2$. We test this prediction for mismatched Drude-Lorentz spectra using nonperturbative hierarchical equations of motion generalized to cross-correlated bath forces. These results provide a geometric framework for dark-state engineering in structured reservoirs.

A change in the geometry of the Hilbert space of a quantum system interacting with a nonadiabatic electromagnetic field is found. To identify its physical origin, the evolution of orthogonality to non-orthogonality of the states of the quantum system has been investigated based on the matrix elements of the quantum states in three main physical cases: a quantum system in zero external field (closed quantum system), a quantum system in slowly changing (adiabatic) field, and a quantum system in rapidly changing (nonadiabatic) field. The quantum states are orthogonal in the first two cases and change to non-orthogonal in a nonadiabatic field, following continuously the parameters of the field. A time-dependent Fubini-Study metric of the Hilbert space can be determined based on the matrix elements of the quantum states in the nonadiabatic case. Such a dynamical metric of the Hilbert space in quantum mechanics resembles the dynamical metric of spacetime in general relativity.

We investigate the critical non-Hermitian skin effect (CNHSE) and the coupling-driven morphing of scale-free localization (SFL) in two Hatano-Nelson chains coupled locally at two bulk sites separated by an inclusive distance d. At weak interchain coupling, competition between the opposite skin accumulations of the two chains generates a geometrically selected critical branch with a complex loop-like spectrum and scale-free skin states localized near the physical boundaries. Remarkably, this SFL branch persists as the coupling becomes strong, although its underlying mechanism and spatial profile change qualitatively. Strong local hybridization generates high-energy impurity states and effectively divides the low-energy Hilbert space into geometry-selected inner and outer segments, with the former supporting SFL states near the interchain-coupled sites and the latter retaining conventional skin states. The number and localization scale of the SFL states are controlled by the separation between the coupling sites. Consequently, increasing the interchain coupling continuously relocates the SFL states from the physical boundaries to the coupling-defined internal interfaces. Our results demonstrate that bilocal coupling can induce a CNHSE and further drive its scale-free states into a segmentation-induced localization regime.

As increasingly-damaging floods can strike within hours of a storm and in ungauged reaches, hourly network-wide simulation has become critical societal infrastructure. Here we demonstrate a multi-timescale physics-embedded learning model which outperforms the United States' operational system and surpasses AI-based systems at flood peaks. Covering more than 800,000 river reaches of the conterminous U.S., the model elevates median hourly Nash-Sutcliffe efficiency at 2,831 gauges to 0.683 from 0.461 for the operational National Water Model v3.0, and narrows flood-peak timing errors from 7-8 hours to 4.5-6 hours. dHBV2.0MTS-MC captures 33% more >=50-year floods than NWM3.0 and 159% more than an operational LSTM baseline. Against recent AI models, overall hourly skill is comparable while rare-flood accuracy is distinctly higher, with relative peak-magnitude error reduced by 34% for >=100-year floods. It combines long-term hydrologic context, short-term shocks, and infiltration excess to resolve extraordinary hourly peaks not visible on a daily plot. Process-model parameters and hourly discharge are produced for every reach, seamlessly covering the continent at 7.2 km2 median resolution. This candidate for the next-generation National Water Model sets a new operational accuracy level for national-scale flood prediction.

A theory for the general circulation of the atmosphere must be based on a theory of its macroturbulence. A central component is the near-surface eddy heat flux, which is closely tied to the mass transport of the surface branch of the circulation. Existing scaling theories are largely formulated within the two-layer quasi-geostrophic framework, which assumes supercritical states and an inverse cascade of kinetic energy—conditions that are frequently not satisfied in planetary atmospheres, including Earth's. Nevertheless, such theories have shown some empirical success, including that the Rhines scale is associated with the mixing length even in the absence of any inverse energy cascade. Here, we analyze hundreds of idealized, dry general circulation model simulations spanning wide ranges of rotation rates, meridional temperature gradients, vertical stratification, and seasonality. Most simulations reside in the marginally critical regime where strongly nonlinear scaling theories fail. We propose a new theory that captures both domain-averaged and local eddy heat flux behavior across all simulations. The theory uses scaling arguments for the zonal momentum balance and accounts for kinematic effects on eddy mixing. The scaling depends on a non-dimensional thermal Rossby number and implicitly accounts for the role of nonlinear eddy momentum fluxes across jets. We demonstrate the skill of the proposed scaling in explaining global near-surface temperature distributions across a wide range of climate states in a simple energy balance model. Implications for general circulation theory and extending the scalings to moist atmospheres are discussed.

AETHER-P\textsuperscript{3} Nowcast v1.0 is a machine-learning-based global thermospheric neutral-density model developed for low-Earth-orbit applications and prepared for onboarding to NASA's Community Coordinated Modeling Center (CCMC). The model provides pointwise neutral-density estimates together with predictive uncertainty using a deep evidential regression framework driven by causal solar, solar-wind, geomagnetic, spatial, temporal, and empirical-model inputs. This report documents the released model configuration, training-data construction, software traceability, output products, validation strategy, benchmark performance, and known limitations. The training archive combines accelerometer- and mission-derived density observations from CHAMP, GRACE-A, GOCE, Swarm-C, and GRACE-FO spanning 2000–2023. More than 40 million eligible 30-s observations are available, but the archive is strongly dominated by consecutive quiet-time measurements. To preserve coverage of physically important regimes, the final 1.67-million-sample training set is constructed using deterministic regime-aware sampling that progressively subsamples quiet conditions while retaining all available extreme-condition observations. Validation uses temporally disjoint chronological blocks with exclusion guards to reduce information leakage. Evaluation across quiet, moderate, and extreme conditions shows competitive performance relative to HASDM, JB2008, NRLMSISE-00, and available WAM-IPE cases, while also identifying limitations associated with sparse training coverage, mission-dependent density products, and condition-dependent uncertainty calibration. The report provides a reproducible technical description of the AETHER-P\textsuperscript{3} Nowcast v1.0 research release and its current CCMC onboarding configuration.

Precision metrology, quantum information, and quantum-controlled chemistry studies all have many implementations in trapped ion systems that require a precise understanding of the motional dynamics of ions. Between the well-studied Coulomb-crystal and gaseous regimes lies the intermediate liquid-like regime whose properties and boundaries have largely eluded experimental investigation. Here, we track the entire liquid-to-gas transition process by leveraging direct velocity measurements using our newly developed ion-trap with integrated velocity map imaging. We demonstrate that the colder the ensemble, the more viscous the ion motion becomes until radial localization emerges, which marks the onset of the liquid regime at a temperature warmer than commonly assigned. At even hotter temperatures, we observe an abrupt change in the heating rate power-law, which indicates the maximum density threshold. The enhanced sensitivity of direct velocity measurements reveals the full transition process from a fluid with short-range order to a weakly-interacting gas.

Near-Earth Objects (NEOs) can appear as faint streaks in long-exposure astronomical images. Detecting these streaks across diverse observatories requires methods that remain reliable despite differences in image quality, orientation, sky background, and noise. However, existing detectors are commonly evaluated using data from only one source, providing limited evidence of cross-source generalization. We introduce NEO-Bench, a multi-source benchmark containing 8,376 images from five astronomical-image datasets. The sources include the Hubble Space Telescope, a Stellina smart telescope, the United Arab Emirates Meteor Monitoring Network, a TETRA1 telescope using a Celestron C14 with Fastar, and the Roboflow Asteroid dataset. We converted the data to a common YOLO format, audited a sample of labels, and defined within-source and leave-one-source-out evaluation protocols. We evaluated four approaches: Hough, Radon, Gaussian PSF, and YOLO26L. Leave-one-source-out F1 decreased in 14 of 20 image-level method-source pairs and 13 of 20 IoU@0.50 localization pairs. Across the datasets categorized as medium or hard, F1 decreased in 11 of 12 image-level pairs and 9 of 12 localization pairs. These results show that cross-source performance remains inconsistent and that reliable generalization across astronomical imaging sources remains an open challenge. The benchmark, code, and data are publicly available.

Population synthesis based on the grid of pre-computed evolutionary sequences for close binaries is applied to the problem of formation of the Galactic population of stripped helium stars, the remnants of $5 \lessapprox M/M_\odot \lessapprox 24$ primary components of close binaries after RLOF. Stellar rotation is taken into account. For the current star formation rate in the Galaxy $2\,M_\odot yr^{-1}$, the estimate of the number of 1 to $7\,M_\odot$ objects with $T_{eff} \gtrapprox $25000K ranges from 14700 to 28500, while for the number of objects with masses from $2\,M_\odot$ to $7\,M_\odot$ it is 3200 to 5500, depending on the assumptions about formation of common envelopes and their outcomes. Other factors influencing the estimate are adopted star formation rate, distributions of close binaries over initial masses of primary components, mass ratios of components, orbital periods. Observational selection may severely limit the number of detectable objects. Angular velocities of rotation of companions to stripped helium stars are 60 to 80 per cent of the critical one. Using the models of initially homogeneous helium stars as approximations for the remnants after mass loss results in the underestimate of pre-supernovae masses.

We investigate a fractional Bose Einstein condensate dark matter scenario in which an ultralight bosonic component forms a self gravitating core embedded in a predominantly cold dark matter halo. The stationary GPP system solved numerically beyond the TF approximation leads to energetically bound equilibrium configurations and good phenomenological fits to the selected rotation curves from 30 SPARC galaxies.

The LUX-ZEPLIN (LZ) Collaboration has reported a $248$ keV nuclear-recoil candidate for which endothermic dark matter (DM) gives some of the largest local significances. We study Solar-capture constraints on three interpretations: a thermal Higgsino, a thermal pseudo-Dirac fermion with off-diagonal vector interactions, and neutron-philic endothermic spin-dependent scattering through $\mathcal{O}_4$. The canonical full-density thermal Higgsino is excluded: its LZ-preferred splitting near $0.38\,\mathrm{MeV}$ lies well below the splitting required to suppress Solar capture sufficiently to satisfy the IceCube upper limits on DM annihilation in the Sun, $\delta\gtrsim0.51$-$0.56\,\mathrm{MeV}$. For the pseudo-Dirac benchmark parameters that fit the LZ event, we find $C_\odot=1.53\times10^{20}\,\mathrm{s}^{-1}$, while two-state kinetics limits the fixed-orbit annihilation rate to $\Gamma_A\lesssim4.4\times10^{14}\,\mathrm{s}^{-1}$, about $1.7\times10^5$ below the IceCube upper limit for the $b\bar b$ channel. A semi-analytic treatment indicates that re-excitation cycles further cool the captured population, although a full phase-space calculation is required for its final distribution. For neutron-philic $\mathcal{O}_4$ scattering at $m_\chi=1$ TeV and $\delta=300$ keV, finite-temperature capture gives $C_\odot=3.98\times10^{17}\,\mathrm{s}^{-1}$, with nuclear-structure uncertainties giving an envelope $(1.77$-$9.72)\times10^{17}\,\mathrm{s}^{-1}$. Under the equilibrium assumption, the upper edge remains about a factor $154$ below the IceCube upper limit for the $b\bar b$ channel. We use $b\bar b$ only as a soft-hadronic proxy because an exact constraint requires the model-specific annihilation spectrum. Solar capture therefore excludes the thermal-Higgsino interpretation but not endothermic explanations generically.

The LZ collaboration has reported a nuclear-recoil excess near 248 keV with a global significance of $2.6\sigma$. Although halo dark matter with a magnetic dipole-dipole interaction and a TeV-scale mass provides the best fit to the excess among the interactions considered by LZ, it predicts a considerable number of events at lower recoil energies, where no excess is observed. We show that a boosted velocity distribution can alleviate this tension and provide a better fit to the LZ recoil spectrum. Moreover, the boost opens up the possibility of explaining the excess with much lighter dark matter, with masses down to the GeV scale. We demonstrate these features first in a model-independent analysis and then realize them in a concrete dark matter model, in which halo dark matter annihilates into on-shell mediators that subsequently decay into boosted dark-sector particles. Our results demonstrate that boosted dark sector particles provide a viable interpretation of the LZ excess.

We investigate superconductivity in monolayer NbSe$_2$, an Ising superconductor, under an in-plane Zeeman field $h_x$, focusing on the emergence of symmetry selected equal-spin triplet pairing. Using a self-consistent Bogoliubov–de Gennes approach with realistic hopping parameters for monolayer NbSe$_2$ on a triangular lattice, we determine the energetically favored singlet pairing states in a range of chemical potentials $\mu$ for onsite, nearest-neighbor, and next-nearest-neighbor pairing channels, and map the resulting phase diagram in $(h_x,\mu)$ plane. A momentum-resolved analysis reveals distinct dominant contributions in order parameters arise from the surroundings of $\Gamma$, $K$, and $K'$ points of the Brillouin zone. We find that an in-plane magnetic field helps to induce a triplet pairing whose symmetry is determined by the parent singlet state, upon opening the triplet interaction channel. For the nonlocal pairing channels, chiral-$d$ pairing is energetically favored and it induces chiral-$p$ equal-spin triplet component with opposite chirality. We further find that Rashba spin-orbit coupling, relevant to substrate coupling and electrostatic gating, suppresses the superconducting orders and reduces the critical field. Our results establish a direct link between the symmetry of the parent singlet condensate and the emergent triplet superconductivity, highlighting monolayer NbSe$_2$ as a promising platform for field-tunable mixed-parity superconducting states.

Recently, a generalized Landau paradigm has been proposed, in which a broad class of "unconventional" phase transitions lying beyond the original symmetry-breaking framework of Landau are understood as the symmetry-breaking transitions of generalized symmetries. In particular, in (1+1)D, such "topological" or "deconfined" transitions can be mapped to the symmetry-breaking transitions of non-invertible symmetries. How to extract the universal dynamical properties of the transitions from this understanding? In this paper, we follow the Ginzburg-Landau philosophy and develop a field-theoretic description of such critical points in terms of the fluctuations of local order parameters, now of non-invertible symmetries. The Symmetry Topological Field Theory formalism plays a central role in our analysis, allowing us to identify the algebra of local order parameters from the Lagrangian algebra of a (2+1)D topological order. We illustrate the key ideas and the full steps of this generalized Ginzburg-Landau procedure using the transitions between gapped phases with a simple non-invertible symmetry — the Rep($S_3$) symmetry.

The recent observation of an event with large recoil energy at LUX-ZEPLIN (LZ) has found a possible explanation in terms of the inelastic scattering of a pseudo-Dirac particle, for example, the Higgsino of the MSSM. For pseudo-Dirac dark matter, the mass splitting $\delta m$ is correlated with the mass of some heavier states. While in the MSSM rough dimensional considerations seem to indicate mass ranges of these heavy particles hopelessly out of reach even at future colliders, we point out that this is not the case in GUT-induced scenarios, with or without supersymmetry. Working in a non-supersymmetric model based on SU(6), we show that the experimentally favored $\delta m$ implies the presence of a $Z'$ gauge boson with a mass of several TeVs. We derive the predicted $Z'$ signal in the dilepton channel and we confront it with the sensitivity of the HL-LHC and the FCC-$hh$ to $Z'$ resonances. We show that the $Z'$ falls squarely in reach of the discovery potential of both experiments.

This paper presents a study of the Weber-Schafheitlin integral and its connections to random walk theory. The Weber-Schafheitlin integral, defined as the improper integral of a product of Bessel functions, arises naturally in the evaluation of probability distributions for multidimensional random walks. This work bridges pure special function theory with applied probability.

The NASA/JPL CNEOS Fireball Data catalogue provides a near-global record of bright atmospheric entries but does not report per-event velocity covariance. We calibrate its kinematic accuracy using 19 events with independent reference trajectories and use optical, satellite, radar, and infrasound catalogues as source-specific controls. For nine modern high-quality comparisons, the CNEOS-minus-reference speed residual has mean -1.00 km s$^{-1}$ and standard deviation 0.86 km s$^{-1}$; variance decomposition assigns 0.77 km s$^{-1}$ to the CNEOS component. The radiant residuals have a 0.74 degree per-axis core plus one 16.3 degree outlier, so the extreme tail remains poorly determined. Seven of 354 complete CNEOS states are nominally hyperbolic without correction. A +0.9935 km s$^{-1}$ atmospheric-stage shift produces five additional crossings, but these are treated only as sensitivity cases. Polar-IM (2026 April 1) is the strongest post-2018 CNEOS anomaly, whereas independent GLM-LI stereo solutions yield 57.3 +/- 2.0 and 56.7 +/- 3.0 km s$^{-1}$, close to their direction-specific parabolic thresholds. We conclude that CNEOS can identify priority candidates, but secure interstellar classification requires an independent, uncertainty-bearing three-dimensional trajectory for the same event.

The inward particle flux, or particle pinch, is routinely observed in magnetically confined fusion experiments, yet its mechanism is not fully understood. We study the particle pinch in the tokamak edge with the flux-driven global turbulence model GDB. Starting from a flat density profile fueled only near the last closed flux surface, the simulation develops a strong inward particle flux that builds a centrally peaked density profile over $O(10)$ ms. The pinch coexists with the usual outward turbulent heat transport; two distinct mechanisms carry it. In the early stage, when density and temperature gradients oppose each other ($\eta_\alpha=L_n/L_{T_\alpha}<0$), drift-wave turbulence drives the inward flux through electron thermal diffusion. In the late stage, once the density profile has flattened, a persistent inward equilibrium $E\times B$ flux, carried by the poloidally asymmetric density and electrostatic potential, drives the central peaking. While the Pfirsch–Schlüter neoclassical transport sets the amplitude of the up-down asymmetric density, classical theory predicts no net radial flux from this asymmetry at leading order. The observed flux flows in a separate non-ideal channel, opened by parallel resistivity, electron inertia, and electromagnetic induction in the equilibrium electron force balance. This channel shifts the equilibrium potential poloidally against the density by $\delta_s\simeq-0.08\pi$ and supplies the late-phase density build-up. These results show that a flux-driven global edge simulation can self-consistently produce a centrally peaked density profile without ad hoc assumptions, and bear on longstanding edge questions such as density pedestal formation.

Conjugated polymer aggregates exhibit optical decoherence in the presence of strong vibronic coupling and substantial diversity in chemical structure, solid-state organization, and excitonic character. Here, we use coherence-detected and population-detected two-dimensional electronic spectroscopies to examine the temperature dependence of the homogeneous optical linewidth across a series of semiconducting polymers. Despite substantial differences in molecular architecture and solid-state organization, all polymers studied exhibit remarkably weak thermal linewidth scaling over the measured temperature range. Comparison between complementary detection modalities further shows that, while this weak temperature dependence is obust, the absolute homogeneous linewidth depends on the measured observable. This behavior reflects the different ways in which coherence- and population-detected measurements project population relaxation and pure dephasing onto the spectroscopic response. These results establish the weak thermal scaling of optical decoherence across a diverse series of conjugated polymers and show that its experimental manifestation must be interpreted in the context of the detection observable.

Fast ionic transport is a defining feature of many solid electrolytes, yet its microscopic origin is not fully understood. In the absence of microscopic insights, the development of next-generation solid-state batteries remains largely empirical. Most existing measurements access either the low-frequency transport response or the high-frequency bound polarization, yet the intermediate mesoscopic frequency regime is where ionic transport emerges. By varying the $\mathrm{AgI}$ concentration ($x$) and performing time-domain terahertz spectroscopy (TDTS) in a prototypical glassy electrolyte $\left(\mathrm{AgI}\right)_{x}\left(\mathrm{AgPO_3}\right)_{(1-x)}$, we reveal this intermediate frequency regime and identify a crossover from bound-current-dominated conduction to conductivity arising from short-range dispersive ionic transport. We find that bound polarization associated with the bond-bending motion of the $\mathrm{P{-}O^- -Ag^+}$ motif is present across compositions but is insufficient to produce ionic transport on its own. Transport emerges only when this polarization is embedded in a sufficiently soft $\mathrm{AgPO_3}$ glassy matrix and accompanied by a high carrier density. These ingredients together take the system from a vibrationally bound response to short-range dispersive motion.

We propose a $N=1$ supersymmetric $SU(6)$ theory with an intermediate $(SU(3)_C\times SU(3)_L\times U(1)_X)$ phase, in which common fields and operators relate the neutrinos, gauge mediation, unification, and symmetry-breaking sectors. The inverse-seesaw Dirac matrix is antisymmetric and thus has rank two, leaving one massless neutrino at leading order. The leading symmetry-allowed correction lifts this zero mode through its projection onto the left and right null vectors of the Dirac matrix, resulting in a quartic suppression by the ratio of the intermediate and ultraviolet scales. The complete unified messenger multiplets preserve the relative one-loop gauge-coupling crossing, while the supersymmetry-breaking spurion entering the messenger sector also determines the scale of lepton-number violation. The coupling governing (331) breaking fixes both the inflationary normalization and the radial mass in the broken vacuum. For the conditional heavy spectrum adopted in the threshold calculations, the gauge couplings remain perturbative. The physical scalar directions transverse to the inflaton are stable within the neutral $((\Phi,\bar\Phi,S))$ subspace. These results provide the analytical scale relations among the inverse seesaw, messenger sector, unified gauge evolution, and inflation, while keeping the assumptions associated with the ultraviolet Higgs and heavy-threshold sectors explicitly.

We study Maxwell-Stefan diffusion with additive friction coefficients $f_{ij}=g_i+g_j$. In mass fractions, the system isolates the constrained pair-friction block; in mole fractions, it is the classical ideal isothermal/isobaric Maxwell-Stefan system at constant total molar concentration. Additivity makes the constrained pair-friction dissipation species-diagonal; conversely, species-diagonality on one interior barycentric constraint space forces a pair-sum representation. At operator level, the positive constrained relaxation operator is a scalar shift of a compression of $G=diag[g_1,\ldots,g_N]$. Its scalar resolvent yields both an explicit constrained inverse and interlacing spectral roots, which form global real-analytic coordinates on the open simplex and whose differentials are left eigen-covectors. In root coordinates the principal part is diagonal, no self-square gradient term occurs, and scalar comparison yields invariant rectangles and separation from the simplex boundary. For regularity we introduce entropy-stabilized one-sided multi-EPD truncations: Euler-Poisson-Darboux entropies cancel mixed quadratic production, while for $N\ge4$ a truncation-weighted mixing-entropy correction supplies transverse coercivity. Caccioppoli and logarithmic estimates, shrinking, and critical mass yield Hölder continuity up to the Neumann boundary. The mixing entropy also symmetrizes the moment system; frozen conormal estimates give spatial Lipschitz bounds. Together with time Hölder control and short-interval maximal regularity, this yields global strong solvability on bounded $C^{2+\alpha}$ domains with $0<\alpha<1$, for each $N\ge2$, $d\ge2$, and $p>d+2$, for all uniformly positive, compatible initial concentrations in the natural trace class. Solutions become classical for positive times and converge exponentially to equilibrium in relative entropy, $L^2$, and $C^1(\bar\Omega)$.