arXiv:2607.13006v2 Announce Type: replace
Abstract: A growing family of indices scores how predictable a series is from its spectrum. Practitioners increasingly read these scores as answering a different question: whether \emph{adding context}, a longer lookback, a retrieval plug-in, or a pretrained model, will help. These are not the same question. The value of context is a property of the operating point, not of the series. Any index built from the power spectrum is invariant under phase randomization, whereas the beyond-second-order value that retrieval and foundation models supply is not, because a phase-randomized series is asymptotically Gaussian. We state this as an impossibility result and isolate it with surrogate pairs that fix the spectrum and the marginal by construction. We then give a label-free, configuration-level diagnostic, the coverage deficit, whose principal term measures beyond-spectrum structure as the gain of analog over linear prediction. On seven benchmarks the prediction holds: window-keyed retrieval's value collapses across surrogate pairs (ECL median $+33\%\!\to\!-35\%$, $p{<}10^{-40}$) while every spectral index stays frozen; a foundation model's value splits into a surviving second-order part and a small beyond-linear margin that collapses; a longer linear window's value survives. Leave-one-dataset-out, the structure term predicts the sign of beyond-spectrum value where the spectral indices trail it, and the reverse holds for the second-order mechanism. We introduce no new forecaster; the contribution is the distinction, a controlled comparison, and a diagnostic for the deployment decision. Code: https://github.com/KurbanIntelligenceLab/SINE
Science Journals
arXiv:2607.13010v2 Announce Type: replace
Abstract: Dermatological practice routinely involves measuring and tracking lesion size, morphology and texture, as critical components of wound or skin cancer screening, monitoring and diagnosis. To accomplish this task, practitioners often image the skin surface with commonly available off-the-shelf camera sensors. This has led to an overwhelming research focus on 2D methods while these objectives naturally benefit from 3D information. In this paper, we demonstrate that dense monocular 3D reconstructions, metric scale measurements and rich surface normal texture estimates are achievable for both dermoscopic and macroscopic cases without the need for additional hardware or multiple captures. We present DermDepth, the first single-view metric scale 3D model for the dermatological domain and D-Synth, the first synthetic dermoscopic dataset with pixel-perfect 3D information. Our experiments show training DermDepth on D-Synth corrects metric scale error from over 16x to under 1.1x for real dermoscopic data, while preserving geometric quality and increasing texture richness. Fine-tuning on a small amount of real clinical samples generalizes our method across three real-world benchmarks spanning the few mm to hundred cm range, diverse skin-tones, chronic wound cases and produces measurements broadly consistent with disease size reported in medical literature. All code, data and models are available at https://github.com/hectorcarrion/dermdepth.
arXiv:2502.18963v3 Announce Type: replace-cross
Abstract: Exceptional points (EPs) are remarkable spectral degeneracies in a non-Hermitian system's parameter space, where both eigenvalues and eigenstates coalesce. Here, we show that in non-Hermitian molecular chiral systems the position of EPs in the parameter space is enantiomer-specific. First, we show that encircling the EP of one enantiomer drives robust topological population transfer in the chiral molecule while its mirror twin remains unaffected, offering a new route for selective chiral control. Second, we reveal how resonant excitation of EPs in chiral molecules can amplify weak chiral effects, offering an alternative approach to the enhancement of chiral interactions. Third, we demonstrate that a twisted chiral fiber immersed in a liquid solution of chiral molecules exhibits topologically different behavior depending on the solution's enantiomeric excess, offering a new approach to the detection of molecular chirality. Our results combine high enantiosensitivity with topological robustness in chiral discrimination and control, paving the way for new approaches in the exploration of non-Hermitian and chiral phenomena.
arXiv:2507.08138v2 Announce Type: replace-cross
Abstract: We present the Conservative Matrix Field (CMF) as a tool for the analysis and computation of D-finite functions. We use conservative matrix fields to establish asymptotic properties of families of linear forms in periods, including (but not limited to) multivariate Mellin integrals, via a discrete Levinson-type framework due to Benzaid and Lutz. Finally, we present an experimental analysis of the families of linear forms generated by these objects and formalize the resulting observations as conjectures on their continuous asymptotic and arithmetic properties.
arXiv:2509.15849v2 Announce Type: replace-cross
Abstract: We consider polaritons formed by hybridizing particle-hole excitations of an insulating phase with a cavity photon at sub-gap frequencies, where absorption is suppressed. The strength of the hybridization is driven by the Van Hove singularity in the joint density of states (JDOS) at the band gap: the stronger the singularity, the more a photon is hybridized with the interband transitions. In order to increase the singularity and thus the polariton hybridization without absorption, we propose to engineer a non-parabolic momentum dispersion of the bands around the gap in order to implement a high-order Van Hove singularity (HOVHS) in the JDOS. Ultracold atoms in tunable optical lattices are an ideal platform to engineer two-dimensional gapped phases with non-trivial band dispersions at the gap. Moreover, the intrinsic non-interacting nature of polarized fermionic atoms prevents the emergence of sub-gap excitations, which are common in solid-state systems and could otherwise spoil the absence of absorption below the gap. Our findings identify band-engineering at the gap edge as a promising route for polariton control with applications in quantum-nonlinear optics.
arXiv:2510.01943v3 Announce Type: replace-cross
Abstract: Quasar-convex functions form a broad nonconvex class with applications to linear dynamical systems, generalized linear models, and Riemannian optimization, among others. Current nearly optimal algorithms work only in affine spaces due to the loss of one degree of freedom when working with general convex constraints. Obtaining an accelerated algorithm that makes nearly optimal $\widetilde{O}(1/(\gamma\sqrt{\varepsilon}))$ first-order queries to a $\gamma$-quasar convex smooth function \emph{with constraints} was independently asked as an open problem in Mart\'inez-Rubio (2022); Lezane, Langer, and Koolen (2024). In this work, we solve this question by designing an inexact accelerated proximal point algorithm that we implement using a first-order method achieving the aforementioned rate and, as a consequence, we improve the complexity of the accelerated geodesically Riemannian optimization solution in Mart\'inez-Rubio (2022). We also analyze projected gradient descent and Frank-Wolfe algorithms in this constrained quasar-convex setting. To the best of our knowledge, our work provides the first analyses of first-order methods for quasar-convex smooth functions with general convex constraints.
arXiv:2510.14317v2 Announce Type: replace-cross
Abstract: Column generation and branch-and-price (B&P) are leading mathematical optimization methods for large-scale exact optimization, iterating between solving a master problem and a pricing problem. Due to the difficulty of discrete optimization, high-performance column generation often relies on a custom pricing algorithm built specifically to exploit the problem's structure. This bespoke nature of the pricing solver makes column generation a problem-specific method and hinders the use of generic implementations across a wide range of problems. We show that domain-independent dynamic programming (DIDP), a model-based paradigm for dynamic programming, can be used as a generic pricing solver. We develop new modeling features and a solving algorithm for DIDP to achieve better performance in typical pricing problems. We demonstrate that in four problem classes, our implementations of B&P, with pricing by DIDP, empirically outperform an existing automated B&P solver and B&P with pricing by mixed-integer programming or constraint programming.
arXiv:2510.15824v2 Announce Type: replace-cross
Abstract: This article considers an online version of conformal inference, called adaptive conformal inference [ACI] and introduced by Gibbs and Cand\`es (2021): prediction sets are issued sequentially, after observing features and before the outcomes are revealed. These sets are evaluated both in terms of validity (the fraction of rounds where the outcome was lying in the prediction set) and efficiency (the average lengths of the prediction sets). The two criteria point to different directions (validity favors larger sets). We also target a wide range of scenarios, with exchangeable data and arbitrary data (lack of any stochastic guarantees) as two extremes. A series of existing strategies for ACI typically guarantee that empirical coverage converges to the desired level for arbitrary sequences, but they generally lack simultaneous efficiency guarantees. To provide a unified study, we first formulate ACI as a repeated two-player game with finite action sets and vector-valued payoffs encoding validity and efficiency. Building on this reformulation, we introduce a strategy based on Blackwell approachability and on its opportunistic extension by Bernstein et al. (2014) that ensures validity while adapting the efficiency of the prediction intervals to the underlying degree of stochasticity of the opponent player. The resulting guarantee is "best of many worlds": it recovers the relevant efficiency guarantees in exchangeable and adversarial settings, and provides guarantees in intermediate settings that arise in typical applications such as the forecasting of time series.
arXiv:2511.01680v4 Announce Type: replace-cross
Abstract: Social scientists are increasingly turning to unstructured datasets to unlock new empirical insights, e.g., estimating descriptive statistics of or causal effects on quantitative measures derived from text, audio, or video data. In many settings, unsupervised analysis is of primary interest, in that the researcher does not want to (or cannot) manually pre-specify all important aspects of the unstructured data to measure; they are interested in "discovery." This paper proposes a general and flexible framework for pursuing such discovery from unstructured data in a statistically principled way. The framework leverages recent methods from the literature on AI interpretability to map unstructured data points to high-dimensional, sparse, and interpretable "concept embeddings"; computes statistics from these concept embeddings for testing interpretable, concept-by-concept hypotheses; performs selective inference on these hypotheses using algorithms validated by new results in high-dimensional central limit theory, producing a selected set ("discoveries"); and both generates and evaluates human-interpretable natural language descriptions of these discoveries. The proposed framework has few researcher degrees of freedom, is robust to data snooping and other post-selection inference concerns, and facilitates fast and inexpensive sensitivity analysis and replication. Applications to recent descriptive and causal analyses of unstructured data in empirical economics are explored.
arXiv:2511.13592v2 Announce Type: replace-cross
Abstract: The existing method of GS-PowerOpt solves the non-convex optimization problem of the form $\max_{\boldsymbol{x} \in \mathbb{R}^d} f(\boldsymbol{x})$ through maximizing a Gaussian-smoothed surrogate $F_{N,\sigma}(\boldsymbol{\mu}) = \mathbb{E}_{\boldsymbol{x}\sim\mathcal{N}(\boldsymbol{\mu},\sigma^2 I_d)}[e^{N f(\boldsymbol{x})}]$. We analyze the role of the smoothing radius $\sigma>0$ and identify a limitation of the fixed-$\sigma$ design used in GS-PowerOpt. Specifically, $\sigma$ induces an inherent exploration--refinement tradeoff: a larger $\sigma$ improves global exploration and finite-time surrogate optimization, but may distort the location of the surrogate maximizer; in contrast, a smaller $\sigma$ better preserves local structure but can weaken gradient signals away from high-value regions.
To address this limitation, we propose GS-PowerHP, a power-smoothed homotopy method with an incrementally decaying $\sigma$ schedule. The proposed mechanism uses larger smoothing radii in early iterations to maintain informative gradient signals when the iterate is far from high-value regions, and gradually decreases $\sigma$ to improve local refinement near the maximizer. We provide theoretical results showing that this decaying schedule improves the exploration--refinement tradeoff of fixed-$\sigma$ power smoothing. Empirically, GS-PowerHP consistently outperforms the fixed-$\sigma$ baseline and exhibits robust performance across different optimization tasks, including adversarial attacks on ImageNet ($d=150{,}528$), where it substantially improves over other smoothing-based zeroth-order methods.
arXiv:2601.02291v3 Announce Type: replace-cross
Abstract: We investigate the angle of repose ${\theta}_r$ of dense suspensions of colloidal silica particles ($d = 2$ $\mu m$ to $7$ $\mu m$) in water-filled microfluidic rotating drum experiments, to probe the crossover between the thermal (colloidal) and athermal (granular) regimes. For the smallest particles, thermal agitation promotes slow creep flows, and piles always flatten completely regardless of their initial inclination angle, resulting in ${\theta}_r = 0$. Above a critical particle size, piles of colloids stop flowing at a finite angle of repose, which increases with particle size but remains below the minimal value expected for athermal frictionless granular materials: $0 < {\theta}_r < {\theta}_{ath} \approx 5.8{\deg}$. We quantify the arrest dynamics as a function of the gravitational P\'eclet number $Pe_g$, which characterizes the competition between particle weight and thermal agitation. Our measurements are consistent with a recent rheological model [Billon et al., Phys. Rev. Fluids 8, 034302, 2023], in which the arrested state stems from a crossover between glass and jamming transitions as the granular pressure in the pile increases relative to the thermal pressure.
arXiv:2601.07926v2 Announce Type: replace-cross
Abstract: We propose a general method for studying systems that display excitations with arbitrarily low energy in their low-temperature phase. We argue that in a rectangular right prism geometry, with longitudinal size much larger than the transverse size, correlations decay exponentially (at all temperatures) along the longitudinal dimension, but the scaling of the correlation length with the transverse size carries crucial information from which the lower critical dimension can be inferred. The method is applied in the particularly demanding context of Ising spin glasses at zero magnetic field. The lower critical dimension and the multifractal spectrum for the correlation function are computed from large-scale numerical simulations. Several technical novelties (such as the unexpectedly crucial performance of Houdayer's cluster method or the convenience of using open - rather than periodic - boundary conditions) allow us to study three-dimensional prisms with transverse dimensions up to $L=24$ and effectively infinite longitudinal dimensions down to low temperatures. The value that we find for the lower critical dimension turns out to be in agreement with expectations from both the Replica Symmetry Breaking theory and the Droplet model for spin glasses. We argue that our novel setting holds promise in clarifying which of the two competing theories more accurately describes three-dimensional spin glasses.
arXiv:2601.20496v2 Announce Type: replace-cross
Abstract: Generating dense physical fields from sparse measurements is a fundamental question in sampling, signal processing, and many other applications. State-of-the-art approaches to this problem either rely on spatial statistics that ignore the governing physics, integrate the physics into a multiple-objective optimization process, or require examples of the complete, fully-resolved simulation state during training, which are frequently unavailable outside of synthetic benchmarks. Here, we present a novel alternative that leverages recent advances in the integration of numerical simulators with data-driven models. Namely, we propose a hybrid modeling pipeline that couples Radial Basis Function (RBF) reconstruction with a Neural Network (NN) correction and a Partial Differential Equation (PDE) solver, so that the numerical simulator itself is embedded directly in the training loop of the learned component. Notably, the NN is trained without assuming availability of examples of the fully-resolved simulation state. This is made possible by implementing the PDE solver so that it is end-to-end differentiable, allowing gradients to be backpropagated through the simulation step during training. This grey-box methodology is evaluated on three standard benchmarks from fluid mechanics, where it achieves superior results over statistical and machine-learning-based reconstruction methods.
arXiv:2602.12880v2 Announce Type: replace-cross
Abstract: Atomic interactions play an important role in the properties of ultracold atomic gases. In single component bosonic systems, its effect is already present at the critical point for the Bose-Einstein condensate phase transition by shifting it to lower temperatures as a consequence of effective repulsion between the atoms. When considering atomic bosonic mixtures, interesting effects arise from the competition between intra- and interspecies interactions such as the miscible-immiscible phase transition and the particular case of self-bounded quantum droplets. In such a scenario, it is natural to expect that these interactions will also affect the critical point of each species composing the mixture. In this paper, we obtain analytical expressions for the critical temperature shift of the phase transition to a Bose-Einstein condensate in the presence of a second species. We treat differently the cases in with the second species is above or below its own critical temperature and apply the obtained relations to the case of a $^{23}$Na-$^{39}$K bosonic mixture which can be realized in current running experimental setups. Our findings can be easily extended to other atomic mixtures trapped by arbitrary conservative traps.
arXiv:2603.17466v5 Announce Type: replace-cross
Abstract: The goal of this study is to introduce a unified framework for simulating random iteration equations (RIE), understood as iteration equations containing random variables. The main idea is to propagate approximations of the full state density from one iteration to the next, rather than estimating it from many repeated pathwise Monte Carlo simulations. The presentation of the RIE modeling framework is conceptually simple based on recent work on static random equations and designed to be accessible. The modeling requirements for RIEs allow for potential nonsmooth nonlinearities and stochasticities in the transfer function. Additionally, the RIE computational strategy for full-density propagation is presented based on iterative likelihood / posterior calculations. As results, illustrative applications of nonlinear random and stochastic differential equation simulations, a new full-density gradient descent method (FDGD) for global optimization under uncertainty and examples of chaotic mappings are presented in order to demonstrate the breadth of the utility of this framework. In total, the character of the presentation is explorative and encourages new applications and theoretical studies.
arXiv:2603.28408v2 Announce Type: replace-cross
Abstract: Energetic particle populations are ubiquitous throughout the Universe. In our solar system, the most prominent sources of energetic particles are solar flares or collisionless shocks often driven by huge eruptions of magnetised plasma called coronal mass ejections (CMEs). Remotely, low energy electrons from the Sun can be observed as solar radio bursts that are produced by accelerated electron beams undergoing beam-plasma interactions. There are still many open questions on the generation of solar energetic particles (SEP): how and where are SEPs accelerated during solar flares and CMEs and how they escape the solar atmosphere? Another important question is: what is the link between the solar radio bursts and the observed SEPs at spacecraft? SKA can provide high-resolution radio images combined with spectroscopic observations to determine the acceleration time, trajectory and escape of low energy electrons from the solar corona. The synergy between SKA and current space missions will help investigate solar activity and energetic particles across a wide range of wavelengths and particle energies. Particle data from spacecraft can be used to make a connection between radio bursts and SEPs by comparing SEP inferred injection times and energies to those of electrons generating radio bursts at the Sun. Radio observations in turn can be used to distinguish between flare and shock acceleration since different radio bursts pinpoint towards different energetic processes. Since the acceleration region and origin of SEPs of various properties is still largely debated, radio observations have the potential to be an invaluable tool in unraveling these processes.
arXiv:2604.21206v3 Announce Type: replace-cross
Abstract: We study the magnetic-field dependence of the interactions between two alkaline-earth(-like) Rydberg atoms, ${}^{88}$Sr and ${}^{174}$Yb. Considering the pair of Rydberg states $|ns,{}^3S_1,m_J\rangle$ and $|(n+1)s,{}^3S_1,m_J\rangle$, we show that the effective Hamiltonian takes the form of an {\it XXZ}-type quantum spin model, as in the alkali-atom case [M. Kunimi and T. Tomita, Phys. Rev. A {\bf 112}, L051301 (2025)]. We find that the behavior of the anisotropy parameter for ${}^{174}$Yb at zero magnetic field is significantly different from that for other atomic species. This behavior arises from the interplay of strong spin-orbit coupling and the resulting multichannel redistribution of F\"orster defects in ${}^{174}$Yb. We systematically calculate the interaction parameters of the {\it XXZ} model in the presence of a magnetic field and show that they can be tuned by the field. As applications to quantum many-body problems, we investigate one-dimensional systems in the large-anisotropy regime and show that the folded {\it XXZ} model can be realized in ${}^{174}$Yb systems without fine-tuning of the field. We also investigate two-dimensional square-lattice systems and show that a supersolid phase can emerge in the ground state at the mean-field level.
arXiv:2604.21448v3 Announce Type: replace-cross
Abstract: This work presents a unified viscoelastic-viscoplastic continuum framework for modeling rate-dependent granular flows across regimes. The formulation incorporates two distinct rate-dependent mechanisms, namely micro-inertia and viscoelastic dissipation, within a single continuum description. A central contribution is an explicit link between the coefficient of restitution and a continuum viscosity, derived from an analysis of wave attenuation in granular assemblies, thereby establishing a direct connection between particle-scale collision physics and macroscopic damping. This relation is introduced while retaining inertia-dependent plastic flow governed by the classical $\mu(I)$ rheology. The constitutive model is constructed by meticulously partitioning elastic and plastic responses within the model and corresponding stress-update routine, such that viscous dissipation governs wave propagation and collisional processes without altering the plastic flow rule. The framework is implemented within the material point method to simulate transient processes involving large deformations, material separation, and subsequent reconsolidation. A range of numerical examples, including steady, transient, vibrational, and impact-driven flows, demonstrates that the model captures wave propagation, diffusion, and rate-dependent granular behavior within a unified continuum setting.
arXiv:2605.05378v4 Announce Type: replace-cross
Abstract: Well-known chaotic maps, such as the logistic and tent maps, have been used to generate cryptographically secure pseudorandomness, yet we know of no efforts which attempt to utilize the Gauss continued-fraction map, a known chaotic map, as a starting point for producing quality pseudorandom output. In this paper, we consider the family of $r$-continued-fraction maps, which generalize the Gauss map, and use them to generate pseudorandom output which outperforms many standard generators, such as the Mersenne Twister, in statistical quality, as ascertained by the use of the Dieharder, PractRand, and TestU01 suites. In this way, we demonstrate the potential viability of these maps as a starting point for novel generators, and provide practical motivation for further study of the properties of both the exact and finite-precision $r$-continued-fraction maps.
arXiv:2605.14777v2 Announce Type: replace-cross
Abstract: Spectrally multiplexed telecom quantum networks require quantum memories combining efficient storage with programmable frequency addressing. An integrated implementation should therefore unite a native telecom transition, efficient storage, and fast on-chip spectral control. Here we demonstrate a cavity-enhanced memory in an isotopically purified $^{167}\mathrm{Er}^{3+}$-doped thin-film lithium niobate microring. Long-lived hyperfine shelving states enable persistent, high-contrast atomic frequency comb preparation with a single-component lifetime of $277.6(52.6)$~s, while cavity impedance matching yields $23.3(5)\%$ on-chip efficiency for 100-ns storage. The intrinsic electro-optic response enables frequency-selective storage and routing at rates up to 20~MHz. We further store and retrieve time-energy-entangled telecom photons, violating an entanglement-witness bound by more than 11 standard deviations. Our results establish erbium-doped thin-film lithium niobate as a programmable light--matter interface for spectrally multiplexed quantum networks.
arXiv:2605.27682v3 Announce Type: replace-cross
Abstract: We study the problem of determining a matrix whose $k$th multiplicative compound, with $k > 1$, is a prescribed matrix $M$. The cardinality of the set of matrices whose $k$th multiplicative compound equals $M$ is characterized in terms of $\rank(M)$. On the one hand, if $\rank(M)\le 1$, it is shown that there exist infinitely many such matrices for which a complete characterization is determined. On the other hand, if $\rank(M)>1$, then there exists a unique matrix -- up to an overall sign -- whose compound is $M$. An algorithm for finding a matrix whose compound equals $M$ is detailed, and its time complexity is analyzed.
arXiv:2606.02528v2 Announce Type: replace-cross
Abstract: Large language models now power robo-advisors and trading agents, yet whether they carry built-in biases toward specific assets is largely untested. We ask three questions: do LLMs systematically prefer certain financial instruments; can an internal representation with causal leverage over those preferences be identified; and does that representation affect downstream financial decisions?
We develop a three-level audit protocol and apply it to Bitcoin. First, a behavioral audit of nine frontier LLMs shows that Bitcoin's ranking among money-like instruments is frame-dependent: models place it around rank 5 of 8 as "reliable money" but near the top under crisis and autonomous-agent frames, and an attribute-swap experiment shows that rankings track functional properties, not names. Second, we open a model's internals: a search across thousands of sparse-autoencoder features in Gemma 3 identifies a dominant Bitcoin-selective feature. Amplifying it shifts the model toward the asset and suppressing it shifts the model away, even when "Bitcoin" never appears in the prompt. Third, we test financial consequences: amplification raises Bitcoin's portfolio share by 5.2 percentage points while suppression lowers it by 4.6 pp, with amplification reallocating within crypto and suppression cutting total crypto exposure.
We characterize this as bounded behavioral leverage (leverage meaning causal influence over outputs, not financial leverage): an identifiable internal feature can be perturbed to move financial choices, but only within measurable limits. The framework links internal representations to external recommendations, validated with random controls and mechanism boundaries. As LLMs become autonomous financial agents, this is a first step toward a behavioral layer for emerging know-your-agent (KYA) standards: knowing what an agent prefers, and how far that preference can be moved.
arXiv:2606.20735v2 Announce Type: replace-cross
Abstract: In a recent paper [arXiv:2604.24783 (2026)], we have proposed a general approach to treat systems with inherent non-Boltzmann-Gibbsian behaviour. Given the extremely high accuracy of our approach, we have adopted the term hyperstatistics. We have applied such a statistical mechanics approach, i.e., hyperstatistics, to the discharge of a capacitor in a RC series circuit, pumping of $^4$He of a closed cycle cryostat, midrapidity data of $p$-Pb collisions at the LHC, as well as for the distribution of accelerations in turbulent systems. Here, we discuss into more details the ground of hyperstatistics. We demonstrate the versatility of hyperstatistics upon applying it to the velocity autocorrelation function in Brownian motion and also regarding its potential to describe brain dynamics.
arXiv:2607.00507v2 Announce Type: replace-cross
Abstract: Quantum contextuality is a vital non-classical resource, yet illuminating the precise mechanisms through which it enables unconditional computational advantages remains a challenge. We translate graph-theoretic formulations of contextuality into an unconditional quantum memory advantage for formal language recognition. We define a promise problem on an exclusivity graph $G$ where any classical finite automaton respecting exclusivity requires $N \ge \chi(G)$ memory states, whereas a QFA requires a memory of dimension $d = \xi(G)$. The gap between these bounds isolates a structural, information-theoretic incompatibility between classical and quantum descriptions that we term \textit{representational contextuality}. For Boolean orthogonality graphs, this exacts an exponential classical memory penalty ($d=\mathcal{O}(n)$ vs $N=2^{\Omega(n)}$). Finally, we demonstrate a sharp algorithmic phase transition: allowing the classical machine a finite confusability of mutually exclusive events reduces this exponential classical memory cost to $\mathcal{O}(n)$.
arXiv:2607.14017v1 Announce Type: new
Abstract: Bio-heat transfer models play a fundamental role in predicting temperature fields during laser-induced thermal therapy (LITT). Among continuum bio-heat transfer models, the Pennes equation remains the most widely adopted formulation. However, by representing blood perfusion as an isotropic volumetric source, it neglects the directional transport of thermal energy associated with blood flow. The Wulff model overcomes this limitation by incorporating blood averaged enthalpy transport directly into the heat flux. Despite its physical significance, the derivation of the Wulff formulation and the assumptions required to obtain its governing equation remain only briefly discussed in the original work. The aim of the present work is twofold. First, the physical formulation originally proposed by Wulff is revisited in order to clarify its derivation and the assumptions required to obtain the governing equation. In particular, the physical origin of the modified heat flux is investigated and an equivalent formulation with an independent metabolic source term is discussed. Second, the resulting model is assessed against the classical Pennes formulation through representative laser-induced thermal therapy benchmark problem. The numerical results show that accounting for blood flow directionality may substantially alter the predicted temperature field and the extent of thermal damage. A dimensionless analysis further leads to the definition of a biological Peclet number, which quantifies the magnitude of directional blood enthalpy transport relative to thermal diffusion. Finally, a simple extension accounting for venous stasis is proposed, mitigating the excessive temperature rise predicted by the original Wulff formulation