Papers

Energetic particle induced geodesic acoustic modes (EGAMs) are axisymmetric perturbations of the radial electric field in tokamak plasmas. They are driven unstable by the phase space nonuniformity of a population of energetic particles (EP). In this paper, the nonlinear oscillation in the amplitude of the energetic-particle induced geodesic acoustic modes is studied by means of the gyrokinetic particle-in-cell code ORB5. Similarities are discussed with the beam-plasma instability (BPI), where a Langmuir wave is driven unstable by phase space nonuniformity of a population of energetic electrons. A similar scaling of the nonlinear oscillation frequency as a function of the mode amplitude is found for the EGAMs and for the BPI, confirming that their nonlinear dynamics is strongly determined by the same physical mechanisms. As a product of this study, a novel diagnostics is proposed for the evaluation of the EGAM intensity in tokamak plasmas.

A computational scheme, which utilizes neural networks, was developed to predict properties of molecular-based materials from chemical structures. The method uses a set of simple algorithms to encode the structure and composition of organic molecules directly into numerical vectors, which is used as input for neural networks. Backpropagation type neural networks are then used to correlate these numeric inputs with a set of desired properties. Calculated results for a series of hydrocarbons, hydrofluorocarbons, and crown ethers demonstrate average accuracies of 0.2-8.1% with maximum deviations of 16-20% for a broad range of thermodynamic, physical, and physical-chemical characteristics (heat capacity, enthalpy, heat of evaporation, boiling point, density, refractive index, stability constants, etc.). In addition, a number of physical and mechanical properties were estimated for polymeric materials and compared with regression analysis. Based on the neural network capabilities of formulating accurate quantitative structure property relationships, a technique called computational synthesis is suggested for performing materials design.

Optical trapping and micromanipulation rely on carefully shaped electromagnetic fields to exert precise forces and torques on microscopic particles. Despite their widespread application in biology and nanotechnology, the absolute physical limits of trapping performance, specifically the maximum achievable optical force and trap stiffness, have not yet been rigorously quantified. This work establishes a general theoretical framework to determine these fundamental bounds in the dipole approximation. By relating the optical force and stiffness to a local Taylor expansion of the electromagnetic field at the particle location, we formulate the performance limit as a solution to a quadratically constrained quadratic program. To evaluate these bounds, we employ two complementary approaches. First, we utilize a complete basis of vector spherical wave functions to determine the absolute theoretical limits of optical force and stiffness permitted by Maxwell's equations in free space, revealing Pareto-optimal trade-offs between stable confinement and directional force. Second, we introduce an aperture-based formulation that restricts the incident fields to those realizable by finite planar apertures. This yields device-consistent bounds directly applicable to experimental setups which rely mostly on electromagnetic beams. The finding that optimized aperture fields can outperform standard Gaussian beams by removing the severe axial bottleneck is particularly important. By comparing these two regimes, we identify the specific spatial modes that contribute to stable trapping and quantify the performance trade-offs inherent to physical beam shaping. This dual framework provides provably optimal bounds for power-normalized optical tweezers and serves as a rigorous benchmark for evaluating realistic beam designs.

We study the cumulants of work in a slow isothermal process for gapped systems. Using the Martin-Siggia-Rose-De Dominicis-Janssen (MSRDJ) formalism and the properties of connected correlation functions, we show that in this process, the $n$-th cumulant of work scales as $1/T^{n-1}$ , where $T$ is the time duration. This result holds generally for arbitrary smooth protocols. Furthermore, we derive the coefficients of the cumulants from equilibrium properties. These coefficients are found to be relevant to thermodynamic geometric tensors.

We study a Mach-Zehnder interferometer (MZI)-based optoelectronic magnonic parametric oscillator (OEMPO) incorporating a YIG-loaded magnonic branch and a tunable phase-shifter branch, enabling systematic investigation of adaptive interferometric oscillator dynamics under distributed phase perturbations. Through analysis of nondegenerate OEPO mode pairs and frequency-pulling behavior, the loop free spectral range (FSR) and effective delay time were quantitatively extracted. Despite the nominally frequency-pinned parametric operation, weak frequency pulling and OEPO mode softening were observed, revealing an additional adaptive interferometric degree of freedom introduced by the MZI architecture. By comparing local and global sampling configurations, we demonstrate that the YIG branch behaves predominantly as a local dispersive resonant subsystem governed by the complex magnonic susceptibility, whereas the phase-shifter branch primarily controls the global interferometric redistribution geometry. Nevertheless, coherent recombination and adaptive regeneration within the loop produce finite cross-coupling between the two branches, resulting in partially synchronized interferometric dynamics and branch-dependent adaptive redistribution. Quantitative complex-Lorentzian analysis further reveals substantial phase-to-amplitude conversion and distinct differences between the OEO and OEPO regimes: the phase-pinned OEPO favors strongly dispersive local YIG response, while the frequency-adaptive OEO exhibits more mixed absorptive–dispersive behavior due to spectral relaxation through frequency pulling. Broadly, the present platform establishes a versatile framework for exploring adaptive nonlinear interferometric physics, coherent phase redistribution, and branch-dependent synchronization phenomena in hybrid magnonic-photonic oscillator systems.

Single-Molecule Force Spectroscopy (SMFS) provides unprecedented insights into biomolecular mechanics, yet the high-throughput generation of force-extension trajectories creates a severe data curation bottleneck. Identifying rare molecular unbinding events within thousands of noise-dominated curves traditionally relies on tedious, non-scalable manual auditing. Here, we present a system-agnostic, interpretable deep learning framework tailored to overcome extreme class imbalance in automated SMFS triage. Utilizing 1D-to-2D rasterized geometric matrices, we deployed a modified ResNet18 architecture governed by an asymmetric Focal Loss objective function. We evaluated this framework on the complex mechanical unfolding pathways of the R. champanellensis cellulosome. Under hyper-imbalanced test conditions where the target interaction constituted only 1.34% of the dataset (13 true events out of 970 traces), the model achieved an overall accuracy of 0.9196 and a remarkable True Positive Rate (Recall) of 0.9231. By implementing an empirically calibrated dual-threshold triage system, the pipeline automatically discarded 880 unambiguous background noise traces , reducing the manual curation workload by over 90% while safely preserving high-value rare data. Finally, Gradient-weighted Class Activation Mapping (Grad-CAM) visually validated that the network's decisions are firmly anchored in the relevant geometric features of the force curves, specifically localizing on the structural unbinding regions, effectively mitigating 'black-box' skepticism. Built for free cloud-based execution, this open-source tool democratizes scalable, highly precise molecular discovery across the biophysics community.

We construct KMS states for the Nelson model, the spinless Pauli–Fierz model, and the Pauli–Fierz model with spin by functional integral representations, and study point-source models after removal of the infrared and ultraviolet cutoffs. If the physical test-function space can distinguish the zero mode, the absence of off-diagonal long-range order, the vanishing of the zero-mode form, the vanishing of the condensate density, the order-parameter criterion, and the triviality of the BEC directions and the BEC ideal are equivalent. We also describe the infrared quotient and the BEC ideal in the resolvent algebra uniformly for the three models, and formulate an operator-algebraic sufficient condition for separately given spatially translation-invariant KMS states.

Multicomponent vapor-deposited organic glasses are essential in organic electronic applications, but achieving controlled component segregation at the nano- and mesoscale remains a challenge, hindering the rational development of high-performance devices. In this study, we investigate binary organic semiconductor mixtures of TPD (N,N'-Bis(3-methylphenyl)-N,N'-diphenylbenzidine) and TCTA (Tris(4-carbazoyl-9-ylphenyl)amine). Despite being miscible in the bulk liquid state, the co-deposited glassy films of these two organic semiconductors exhibit a range of segregation behaviors, from homogenous to clearly phase-separated structures. We employed differential scanning calorimetry and resonant soft X-ray scattering (RSoXS) to study the component segregation behavior and used the National Institute of Standards and Technology RSoXS Simulation Suite, paired with Atomic Force Microscopy, to interpret the energy-dependent RSoXS spectra. Our results indicate that component segregation in co-deposited TPD-TCTA films is due to a kinetically-arrested nucleation-and-growth mechanism, in contrast to the segregation mechanism of a previously reported TPD-DO37 (disperse orange 37) mixture which is strongly immiscible in bulk. This work provides a demonstration of tunable molecular aggregation in organic semiconductor glasses, enabling access to a continuum of morphologies from homogeneously mixed to segregated phases.

We introduce a unified theoretical framework for single-plasmon transport in one-dimensional nanowires, bridging the quantized electromagnetic Green's tensor formalism with effective non-Hermitian Hamiltonian models. This approach naturally incorporates propagating surface plasmon polaritons, high-order modes dissipative channels, and intrinsic losses. We investigate both the stationary regime and the spatio-temporal dynamics of a single-plasmon pulse travelling through an atomic chain coupled to a dispersive nanowire. We analyze modal contributions to reflection and transmission spectra for quantum emitter coupled to a silver nanowire, a configuration proposed as a single-plasmon transistor, and we demonstrate that optimized multi-emitter systems offer significant advantages. In case of one quantum emitter coupled to a silver nanowire at telecom wavelengths, we predict a single-plasmon transmittivity down to 7% under realistic conditions, and an atomic qubit population of 12%. Extension to multi-emitter systems using a Löwdin orthogonalization procedure enables a consistent treatment of collective interactions. We show that optimized positioning with just five emitters enhances plasmon modulation, achieving a transmittivity of 2% but also reduces coupling losses to one-third compared to the single-emitter case. Our results establish a robust foundation for analyzing and designing plasmonic waveguide quantum electrodynamics systems.

We consider the $SO(d)$ and $O(d)$ invariant sectors of the bosonic $d$-matrix harmonic oscillator with $U(N)$ gauge symmetry. The micro-canonical degeneracy $\mathcal{Z}( N , d , k )$ for fixed energy $k$ is expressed as a pairing between an $N$-dependent vector and a $d$-dependent vector in the space of partitions of the integer $k$. This pairing formula is derived by counting invariant words in multi-matrix variables $X^i_{j,a}$, using properties of Clebsch-Gordan multiplicities (Kronecker coefficients) for the symmetric group $S_k$, Schur-Weyl duality and harmonic analysis on the homogeneous space $U(d)/SO(d)$. Analytic formulae for large $N$ and $k$ with $ k \le N $ are obtained using group integrals over $U(N)$ and $SO(d)$ (or $ O(d)$). The micro-canonical heat capacity in this regime is negative and turns positive, at a critical value $k_{\rm crit}$, due to finite $N$ modifications to the counting, thus forming what we denote as a characteristic caloric fold in the $ E $ versus $T$ curve. Data from the pairing formula is well fitted by $k_{\rm crit} \sim { N^2 \over 4 }$ for small values of $d$. A derivation of this large $N$ formula is given using a matrix model approximation and semi-classical analysis of the eigenvalue density. The large $N,d$ limit of the degeneracies reveals a key role for ribbon graph combinatorics. The caloric fold is also notably a property of black hole thermodynamics in anti-de-Sitter spaces. We propose the spherically symmetric \(SO(d)\) and \(O(d)\) invariant sectors of \(d\)-matrix quantum mechanics as tractable matrix systems for capturing key features of dual descriptions of black-hole thermodynamics.

Can a general-purpose large language model design molecules with the precision of a seasoned chemist? Current LLM-based frameworks answer this question with scalar feedback loops - generate, score, reject - that amount to informed trial-and-error. Here we show that replacing a single number with the full physicochemical rationale from first-principles calculations transforms the LLM from a stochastic sampler into a causal reasoner. Our system couples retrieval-augmented generation with a self-reflection module that feeds orbital energies, atomic charges, and electron densities - rather than compressed scores - back into the design loop. On HOMO-LUMO gap targets from 2.0 to 5.0 eV, this structure-property-relationship (SPR) reflection achieves a deviation as low as 0.0014 eV with a 100% success rate under the SPR+RAG configuration, consistently outperforming scalar-feedback and non-reflective baselines in median and mean deviation. The framework generalizes seamlessly to dipole-moment design, synthetic accessibility optimization, and molecular docking, and proves robust across 7 distinct LLM backbones. These results establish a new paradigm: when the model understands not only that a molecule fails, but why, iterative molecular design becomes genuinely mechanistic.

Millimeter-wave Adaptive Optics (MAO) is essential for high-precision large-aperture submillimeter telescopes, requiring real-time compensation of wavefront errors by capturing them as spatially-discrete excess path length (EPL) fluctuations. This paper presents a unified control-theoretic framework for the EPL compensation problem. We first model the optical drive system as a plant where input commands relate to measured EPL through a first-order system representing mechanical response delay and a measurement matrix characterizing the actuator-to-sensor coupling. We mathematically formulate the control task as an asymptotic disturbance suppression problem, specifically targeting low-frequency disturbances such as thermal and wind-induced deformations. Second, we propose an Anti-Windup Proportional-Integral (AWPI) control law. By employing a decoupling strategy, the design is reduced to a loop-shaping problem for decoupled scalar sensitivity functions, ensuring both stability margins and asymptotic disturbance suppression of constant-valued disturbances. The anti-windup mechanism is integrated to maintain control continuity during the recovery from saturation, preventing undesirable discontinuities in the drive command. Third, we introduce practical operational tools: a manual focus adjustment scheme that allows observer intervention without interfering with the feedback loop, and the cosine similarity index to quantify the suppressibility of specific Zernike modes. Numerical simulations, incorporating a three-axis secondary reflector drive and five-point EPL measurements, demonstrate direction-dependent disturbance rejection and the suppression of von Karman-modeled wind turbulence, validating the effectiveness of the proposed framework for real-world telescope applications.

We revisit the effective coupling of the scalaron to gauge fields in $f(R)$ gravity minimally coupled to the Standard Model, focusing on the scalaron decay into two photons. Treating the scalaron as an intrinsic component of the Jordan-frame metric, we analyse quantum effects arising from its coupling to the energy–momentum tensor. In this framework, the trace anomaly contributes to the scalaron–gauge boson interaction. Using Fujikawa's approach, we obtain the trace-anomaly contribution, associated with the transformation of matter fields between the Jordan and Einstein frames, first in QED and then in the full Standard Model. The resulting effective scalaron–photon interaction agrees with direct perturbative calculations that include the trace anomaly and differs from the result obtained when only the classical expression for the trace of the energy–momentum tensor is taken into account in the scalaron interaction. In the limit where the scalaron mass is much smaller than the masses of particles circulating in the loop, the diagrammatic contribution from the classical expression of the trace of the energy–momentum tensor cancels the anomaly-induced term, leading to a vanishing effective coupling and a strong suppression of the scalaron decay rate into photons. These results clarify the origin of discrepancies in the literature concerning the effective scalaron coupling to massless gauge fields. They stem from inequivalent prescriptions for incorporating quantum loop effects in $f (R)$ gravity, leading to different scalaron–gauge field interactions with direct implications for scalaron dark-matter phenomenology.

The planar Hall effect (PHE) is an anisotropic magnetotransport response generated by coplanar electric and magnetic fields. We investigate the PHE in single- and bilayer two-dimensional electron gases (2DEGs) with Rashba spin-orbit coupling and identify two distinct mechanisms: Zeeman coupling and a band geometric channel. In the Zeeman channel, an in-plane magnetic field distorts the Rashba spin-orbit-coupled band dispersion and generates anisotropic carrier velocities, producing a finite PHE. In an asymmetric Rashba bilayer, interlayer electronic delocalization generates finite planar Berry curvature and orbital magnetic moment components, giving rise to a band geometric PHE channel. Using semiclassical Boltzmann transport theory, we calculate the chemical potential and angular dependence of the planar Hall conductivity for both mechanisms. Symmetry analysis shows that the leading response is quadratic in the magnetic field and exhibits the characteristic $\pi$-periodic angular dependence. For the parameter regime considered here, the Zeeman-induced contribution dominates, while the band geometric channel provides a distinct symmetry-allowed contribution unique to asymmetric Rashba bilayers. Our results reveal microscopic origins of anisotropic magnetotransport in spin-orbit-coupled two-dimensional materials.

The production of a top-antitop quark pair in association with a $W$ boson constitutes one of the heaviest final states currently studied at the Large Hadron Collider (LHC) at CERN. Measurements of its production rate have consistently exceeded Standard Model predictions. Owing to the complexity of the two-loop amplitudes entering the double-virtual correction, next-to-next-to-leading-order (NNLO) QCD calculations for this process have so far employed dynamical approximations for the two-loop contribution. We present NNLO QCD predictions based, for the first time, on a direct computation of the required two-loop amplitudes in the generalised leading-colour limit.

Understanding quantum dynamics in curved spacetime is a central challenge at the intersection of quantum mechanics and gravity. Anti-de-Sitter (AdS) spacetime plays a pivotal role in the context of the AdS/CFT correspondence, which relates gravitational dynamics in the AdS bulk to a conformal field theory (CFT) living on its boundary. Despite its foundational importance, direct experimental access to dynamical quantum phenomena in Lorentzian AdS spacetime has so far remained out of reach. Here, we report the first experimental emulation of fermionic wave packet dynamics in Lorentzian AdS spacetime using a photonic platform. By mapping the Dirac equation in curved spacetime onto the propagation of light in engineered wave\-guide arrays, we directly observe gravitational confinement of relativistic wave packets and resolve their center-of-mass motion in real time. We identify a characteristic superposition of slow geodesic oscillations governed solely by spacetime curvature and fast Zitterbewegung arising from relativistic particle–antiparticle interference. While the geodesic frequency is independent of fermion mass, the Zitterbewegung frequency exhibits a distinct joint dependence on mass and curvature, revealing a curvature-induced modification of relativistic quantum dynamics. Our results provide the first quantitative experimental access to fermionic bulk dynamics in emulated AdS$_2$ spacetime with Lorentzian signature. This establishes a scalable analog platform that may potentially be used for exploring dynamical aspects of holography.

Accurate acoustic modelling of the skull is essential for simulation-guided transcranial focused ultrasound (tFUS), but commonly used skull parameterisation strategies differ in complexity and reported accuracy. This study experimentally compared five k-Wave skull models: two voxel-wise linear mapping models, two three-layer models, and one single-layer fixed-parameter model. Nineteen regions of interest from five historical and two Thiel-embalmed human skulls were tested at 220 kHz, 680 kHz, and 1000 kHz. Bowl-surface source fields were reconstructed using acoustic holography, and simulated intracranial pressure fields were benchmarked against needle-hydrophone measurements. Across frequencies, mean peak-pressure errors ranged from 20% to 31%, whereas intensity errors reached 41% to 77%. Errors in -6 dB focal volume ranged from 11% to 67%, and focal-position discrepancies were typically several millimetres. Simulations generally predicted smaller insertion losses than measured, indicating a tendency to underestimate skull-related attenuation and overestimate transmitted intracranial exposure. The linear mapping model with fixed attenuation gave the lowest frequency-averaged pressure error, but no model showed a consistent advantage across all metrics. These results show that current skull models can reproduce gross intracranial beam patterns while retaining substantial quantitative uncertainty in exposure, focal coverage, and target localisation.

Ecological communities with structured interactions exhibit collective phenomena such as percolation and clustering of occupied sites. While these effects have been documented in experiments and simulations, systematic analytical understanding has remained limited. In this paper, we develop a dynamical theory of these phenomena for competitive ecological systems defined on random interaction graphs. We introduce a discrete version of the generalized Lotka-Volterra model that preserves key macroscopic features of continuous ecological dynamics while enabling analytical treatment. Within this framework, we characterize the emergence of percolating clusters and describe the spatial organization of surviving sites. Our analysis uncovers which equilibria can be reached by the dynamics and shows how this dynamical accessibility governs the onset of clustering and percolation. In doing so, our framework complements classical Lotka-Volterra theory by providing a dynamical perspective on the collective organization of structured communities.

We investigate collisionless power absorption in resonant, low$-$pressure capacitively coupled plasmas (CCPs). In these radio-frequency (RF) discharges, the sheath capacitance almost exactly balances the plasma inductance, driving the total RF discharge voltage down to just a few volts. However, plasma persists not only in this ultra$-$low$-$voltage regime; it also generates ions that strike the electrodes with kinetic energies substantially exceeding the amplitude of the applied RF voltage. This counterintuitive behavior arises from the presence of a pronounced electrostatic potential well of approximately 40 V within the plasma bulk, which confines electrons while simultaneously accelerating ions toward the electrodes. We show that, under these resonant conditions, collisionless electron heating exhibits a fundamentally different behavior from the conventional paradigm of stochastic sheath heating mediated by electron$-$sheath interactions. Instead, the predominant energy transfer mechanism is bulk electron heating in RF electric fields via a primarily collisionless process that emerges from the synergistic action of: (i) a strongly amplified RF electric field within the plasma bulk, (ii) electron oscillatory motion (bouncing) within the plasma potential well, and (iii) electron scattering resulting from collisions with neutral atoms. Collectively, these phenomena give rise to a pronounced high$-$energy tail in the electron energy distribution function and thereby lead to a substantial enhancement of the ionization rates. As the gas pressure rises, the resonance is disrupted. At the same time, the region of maximum power absorption moves from the plasma core toward the edges and the sheath, which is accompanied by the disappearance of the high$-$energy electron population and a corresponding decrease in ionization rates

Dark photons are a gauge boson of a hypothetical dark sector, representing one of the most-studied minimal extensions of the Standard Model, with wide-ranging theoretical and observational implications. Here, we consider scenarios in which an initially unpopulated dark photon sector is populated via resonant photon to dark photon conversions. This process leads to observable spectral distortions in the cosmic microwave background (CMB), that can be used to constrain these models. We extend previous spectral distortion studies of the monopole spectrum to anisotropic spectral distortions, using the newly developed Frequency Hierarchy (FH) framework of CosmoTherm. We illustrate the physics by presenting detailed computations of the photon transfer functions and distortion cross power spectra throughout the dark photon parameter space. We find that the dark photon mass explicitly controls the shape (i.e., multipole-dependence) of the signal power spectra, while the overall amplitude of the signal is determined by the kinetic mixing parameter of the model. Using these results, we place complementary limits on the minimal dark photon model using data from Planck, finding that the constraints are only marginally weaker than those obtained with COBE/FIRAS data for the average (monopole) distortion. In addition, we compute the corrections to the standard temperature field, arguing that conversions at redshifts larger than $2\times 10^6$ may add iso-curvature type perturbations, which could lead to novel constraints in regimes where distortion anisotropies thermalize.

In the presence of stratification, turbulent flows exhibit intermittency not only at small scales but also at large scales, comparable to the mean flow, as observed in the atmosphere and oceans. We study such flows through a large parametric exploration using direct numerical simulations of the Boussinesq equations with different forcing types. We examine two Prandtl numbers (1 and 6) and vary the Froude number ($Fr$) over a range of geophysical interest values, $0.01\le Fr \le 1$, corresponding to a variation in terms of the buoyancy Reynolds number ($R_{IB}$) of $0.06\le R_{IB} \le 2300$. We analyze the dependence on $R_{IB}$ of the buoyancy flux ($B_f$), the mixing efficiency, the shear parameters, and the vertical momentum flux. Strongly non-Gaussian tails in the spatio-temporal distribution of the $B_f$ are observed, with kurtosis reaching $\approx 10^2$, indicating the potential for stratified geophysical flows to be characterized by highly variable transport properties along the direction of gravity even under stable stratification. This is associated with long-time intermittent behavior of vertical velocity and temperature at large scale, which produces local turbulence and enhances dissipation and transport. We present evidence that the skewness of $B_f$ increases with $R_{IB}$ as a power-law and saturates in the passive-scalar limit. We also show that the domain-averaged $B_f$ exhibits two distinct trends: logarithmic growth with $R_{IB}$ and approach to a small offset as stratification strengthens. A simple model for the temporal evolution of energy and $B_f$ indicates that the defect between vertical and potential energy drives strong $B_f$ events. This trend directly leads to convective instabilities, the formation of two-dimensional and three-dimensional eddies, and rapid dissipation on a turnover timescale, allowing the energetic cycle to restart-also occurring in bursts.

Magnetic equilibrium reconstruction provides the plasma state estimate required for real-time shape control in tokamaks. We present a fast, physics-informed neural network surrogate of the liuqe equilibrium reconstruction code \cite{liuqe1} for the TCV tokamak at EPFL, achieving inference times below 100~$\bm\mu$s and enabling 10~kHz shape control. The model is trained on around 10,000 TCV discharges spanning the full operational range of plasma shapes. Its modular branch/trunk architecture decouples magnetic measurement encoding from spatial coordinate processing, enabling physics-informed regularization via automatic differentiation of the predicted flux map. The surrogate has been compiled and deployed on the TCV real-time control system, and validated both offline and in real time against the models liuqe-rt and lih, showing comparable accuracy. Closed-loop performance assessed with the real-time software in-the-loop fge \cite{fge1} demonstrates control-equivalent behavior across multiple control strategies.

Accurate stellar ages and masses are essential for interpreting the demographics and physical properties of exoplanets, particularly for intermediate-mass, early-type stars where conventional age indicators are ineffective. Isochrone fitting remains the primary tool for characterising such stars, yet its uncertainties are often underestimated, especially in the presence of rapid rotation and unresolved binarity. We present a population-synthesis framework designed to quantify realistic mass and age uncertainties for intermediate-mass stars (1.4-2.5 M$_{\odot}$), incorporating distributions in rotation rate, mass, metallicity, binarity, inclination, and observational error. Rotational and geometric effects are applied a posteriori to stellar evolutionary models, enabling a continuous treatment of rotation and its impact on effective temperature and luminosity. By comparing synthetic populations against commonly used isochrone grids, we demonstrate that rotation and unresolved companions systematically bias inferred masses and ages, particularly for young stars, and introduce random uncertainties at the $\sim$0.1-M$_{\odot}$ and $\sim$180-Myr level, often exceeding formal fitting errors. The effect is strongest near the zero-age main sequence, where ages are underestimated by a factor of $\geq2$, while for older A stars ($>$10% of their main-sequence lifetime), ages are overestimated by 31% with 27% scatter. Our findings carry important consequences for planet detectability, characterisation, and population studies. We provide a publicly available tool, RAPID, for probabilistic inference of stellar parameters from these synthetic populations, and we demonstrate its application to known exoplanet hosts.

We study the possibility of significant PBH production from a slow first-order phase transition with delayed reheating. Since delayed reheating results in an early matter-dominated phase between percolation and reheating, we developed a peak-theoretic approach to PBH formation during this phase based on the non-Gaussian distribution of overdensity arising from the transition. To obtain the collapse probability, we performed large-scale Monte Carlo simulations and employed the hoop-conjecture criterion. We include tidal-torque terms to investigate the initial spin of the PBHs and find that the average spin parameter is $\mathcal{O}(10^{-3})$. Furthermore, we obtain an emergent overdensity threshold for collapse that depends on the phase transition properties and reheating efficiency. We find that the resulting PBH abundance is extremely sensitive to the reheating efficiency, with order-unity changes in efficiency leading to variations of many orders of magnitude in the collapse fraction. We identify regions of parameter space where the resulting PBHs can account for the entirety of the dark matter abundance. Finally, we also constrain the phase transition and reheating properties from current data on (non-)observations of PBHs.

The inner edge of the dead zone (DZIE) in protoplanetary discs is prone to episodic instability caused by the activation of the magneto-rotational instability (MRI) in the weakly turbulent regions. We show how different magnetic field configurations set the inner disc structure and regulate the morphologies of instability cycles. We performed 2D and 3D radiation hydrodynamic simulations of the regions around the DZIE of a Class II disc over a thousand-year timescale. We implemented MRI activation criteria based on ambipolar and Ohmic diffusion coupled to magnetic field strength profiles comprising stellar and disc components. The properties and consequences of the episodic accretion events are highly sensitive to the magnetic field strength. We recover previously reported behaviour by considering relatively strongly magnetised discs. A new burst mode is revealed, in which the midplane MRI activity is restricted to small radii in the presence of weak magnetic fields. In this narrow mode, the pressure bump at the DZIE does not remain static even during quiescence. A distinct dichotomy between the wide and narrow modes is established by the hydrodynamic (in)stability of the ionisation front. Both modes are additionally separated into a reflaring and a non-reflaring version. Our setup does not lead to the classical thermal instability by hydrogen ionisation. In quiescence, the MRI active region shows a layered structure that converges towards the midplane near the star. Our 3D model reveals the breaking of density features produced in the narrow mode, leading to vortices at radii smaller than 0.5 AU. Coupling MRI activity directly to different magnetic field strengths, rather than using simple temperature thresholds, enables a variety of burst modes. Each mode exhibits characteristic accretion burst signatures and has different consequences for planet formation and migration conditions.

Denial-of-service (DoS) and distributed denial-of-service (DDoS) mitigation requires separating malicious traffic from benign traffic while minimizing disruption to legitimate users. Prior work proposed mapping honeypot traffic partitioning to a weighted MaxCut problem and solving the resulting graphs with variational quantum algorithms. We extend this proof of principle direction with a reproducible event-level honeypot-to-QUBO pipeline, labeled temporal bipartite benchmark graphs with 16, 32, 66, and 110 event nodes, QAOA executions on IBM quantum hardware, classical heuristic baselines, a noiseless matrix product state reference, and a routing overhead analysis across quantum processor architectures. The largest benchmark is a 110-node, 181-edge instance executed on three IBM backends. Our results show that a shallow QAOA can execute real traffic partitioning workloads at the utility scale, while backend architecture and routing overhead affect objective quality, security metrics, and observed runtime. Because simple classical heuristics can solve the current labeled benchmark graphs, these experiments are not a quantum advantage claim. Instead, we deliberately use a fixed, shallow QAOA implementation to enable controlled comparisons across problem sizes and hardware architectures. This work establishes a hardware feasibility and architecture benchmark framework, and demonstrates that MaxCut cost, security quality, routing overhead, and runtime must be reported as separate metrics for cybersecurity relevant quantum optimization.

Successfully performing direct numerical simulations of polymeric flows remains a major challenge in computational fluid mechanics. In addition to the velocity field, such simulations must resolve polymeric degrees of freedom, often expressed via the conformation tensor, $\mathbf{c}$, which captures the local stretch of polymer molecules. A key difficulty here lies in maintaining the physical requirement $\mathrm{Tr}\, \mathbf{c}>3$, which is not explicitly enforced by the governing equations. Consequently, simulations initiated from physical conditions may silently drift into unphysical states with $\mathrm{Tr}\, \mathbf{c}<0$, indicating a loss of positive-definiteness of the conformation tensor. Existing numerical methods to prevent this are costly, making direct numerical simulations of chaotic polymer flows, such as elastic turbulence, heavily reliant on high-performance computing. Here, we ask whether simulations that violate $\mathrm{Tr}\, \mathbf{c}>3$ can still yield meaningful physical insight into the underlying dynamics. We simulate a model dilute polymer solution driven through a plane channel at low Reynolds number and observe the transition to elastic turbulence. Our simulations exhibit two threshold resolutions: below the first, they become numerically unstable and exhibit a finite-time blow-up; above the second, they maintain positive-definiteness. In between, simulations remain stable and chaotic despite local violations of $\mathrm{Tr}\, \mathbf{c}>3$. Surprisingly, these violations do not affect mid-plane statistics of velocity, its gradients, or polymer stretch, which match results from fully positive-definite simulations. This suggests that resolving flow structures or key flow statistics may not require the extreme resolutions needed to preserve positive-definiteness, potentially lowering computational barriers for studying elastic turbulence.

The centre comonad model provided the first concrete cohesive linear $\infty$-topos, settling an open problem of Schreiber. However, the model is degenerate: the quantum modality annihilates all non-commutative algebras, and the associated linear logic collapses to classical cartesian logic. In this paper we give a complete mathematical diagnosis of this degeneracy. We prove that the centre comonad sends the representable sheaf of a simple non-commutative algebra to the empty presheaf, and that the state space of any such algebra is empty. We then prove that the Day convolution on the classical core is cartesian, forcing the Seely isomorphism to hold trivially and collapsing the linear logic. We isolate the structural reason behind this collapse: whenever the opposite of the classical core is monoidally equivalent to a cartesian monoidal category, any coreflective precomposition comonad will exhibit the same degeneracy. We conclude that a non-degenerate quantum modality must be constructed without precomposition, and we briefly discuss possible directions.

Observations of wide post-interaction binaries show an unexpected feature; orbital eccentricity, which is not understood theoretically. A promising resolution to this long-standing puzzle is eccentric mass transfer (MT). Here the first complete framework for MT in orbits with arbitrary eccentricity, the general mass-transfer (GeMT) model, is confronted with the latest observations of hot subdwarfs of spectral type B (sdB) with main-sequence (MS) companions in wide orbits. SdBs are excellent benchmarks for binary evolution models, since their progenitors provide unique constraints on their formation histories. We isolate the effects of eccentric MT from other physical process and we show that it explains the observed orbital-parameter distributions and their correlations seen in wide sdB+MS binaries. To quantitatively constrain their orbital parameters, stellar evolution and tides should be included in future studies, yet it is encouraging that the GeMT model provides the first mechanism that qualitatively reproduces the observed distributions. We further demonstrate that post-MT eccentricities depend directly on key MT parameters, including transferred mass, accretion efficiency, and angular momentum loss. Given the multitude of eccentric post-MT binaries with components ranging from low- to high-mass stars to compact objects, we propose that post-MT eccentricities offer a new window onto binary evolution, presenting a powerful tool to constrain highly uncertain binary-evolution parameters and mass-transfer formation histories across diverse populations. Post-MT eccentricity should therefore be embraced as a key observable, rather than treated as a problem to be corrected.

Recently, a number of spin chain models have been discovered that are solvable via hidden free-fermionic structures, going beyond the Jordan-Wigner paradigm. However, all examples in the literature displayed degeneracies that grow exponentially with the volume and that are homogeneous in the spectrum (identical degeneracies for all energy levels). In this note we present a model that can be solved by “free fermions in disguise” (FFD), such that the spectrum is free from exponential degeneracies for generic coupling constants. The model can be seen as a particular perturbation of two Ising chains. Alternatively, it can be realized as an interpolation between a standard Jordan-Wigner solvable chain and the original FFD model of Fendley. We used ChatGPT Pro 5.4 and 5.5 as a research assistant; in the Supplemental Material we provide details about the collaboration between the AI and the human author.