Papers

We determine the optimal information–disturbance tradeoff for estimating an unknown spatial direction encoded in two antiparallel spins. Rotational covariance reduces the optimization over all instruments to a finite-dimensional Choi problem: a positive seed operator obeys one trace constraint for each irreducible sector of the input representation, while both the directional score and the operation fidelity are linear functionals of this seed. For two antiparallel spin-$1/2$ particles, whose physical representation decomposes as $0\oplus1$, we derive the two-multiplier dual problem and characterize the optimal instrument from the kernel vectors of the dual slack operator. The optimal operation is a covariant filter with scalar–vector coherence and is generally not a convex interpolation between the identity channel and a measure-and-reprepare strategy. At maximum information we recover the Gisin–Popescu score, but the least disturbing output state is optimized independently, giving a smaller disturbance than both the parallel-spin benchmark and antiparallel measure-and-reprepare. We also formulate the parallel benchmark and, as a central extension of the method, treat antiparallel spin-coherent states of arbitrary spin $j$. In this case the signal coherently occupies all sectors $\ell=0,\ldots,2j$ of $j\otimes j$, the endpoint information is governed by nearest-neighbor sector coherences, and the endpoint disturbance is obtained from an explicit finite block-diagonal eigenvalue problem.

Quasi-optical coupling serves as the critical interface in terahertz (THz) heterodyne receiver systems, enabling efficient transfer of incident radiation to photomixers through a focusing element and a planar microwave antenna. With recent advances in nanofabrication, planar dielectric metalenses have emerged as promising alternatives to conventional refractive optics due to their compactness and scalability. However, unlike conventional elliptical silicon lenses that are often treated as nearly ideal optical components, the focusing characteristics of metalenses, including both phase and amplitude, strongly depend on the local deflection angle across the aperture, creating an urgent need to quantitatively understand the coupling between a dielectric metalens and a planar antenna. In this work, we present a quasi-optical coupling analysis between a planar Si metalens and a logarithmic spiral antenna integrated with a THz superconducting NbN hot-electron bolometer (HEB) mixer operating at 1.63 THz using spherical-coordinate vectorial integration. By combining the angular radiation profile of the spiral antenna with the propagated complex electric-field profile from metalens numerical simulations, the calculated coupling efficiency accounts for angular power distribution, phase-front matching, and polarization-dependent vectorial overlap. The calculated coupling efficiency is then directly correlated with experimentally measured double-sideband receiver noise temperatures through comparison with a conventional elliptical Si lens measured under the same receiver configuration. The analysis establishes a quantitative relationship between metalens focusing efficiency, vectorial antenna coupling, and receiver noise temperature, providing guidance for optimizing metalens design and improving the overall performance of metalens-integrated THz heterodyne receivers.

Recently we showed that minimum weight decoding in the (6.6.6 planar) colour code is NP-hard. However, it remained an open question as to whether it was possible to approximate the minimum weight decoding arbitrarily closely in polynomial time. In this paper we prove that it is possible: for any $\varepsilon>0$ there is an polynomial time algorithm that, given a syndrome, can find an error-set generating that syndrome whose weight is at most $1+\varepsilon$ times the weight of the minimum weight decoding. As a consequence we see that, for any $\varepsilon>0$, there is a polynomial time algorithm that can correct all errors of weight up to $(1-\varepsilon)d/2$ in the distance $d$ colour code (so almost up to the theoretical $d/2$ limit). The polynomial we give is impractically large, but it does open the door for sensible polynomial time algorithms that approximate minimum weight decoding and, in particular, shows that approximate decoding is not NP-hard.

We propose a new weak-form Physics-Informed Neural Network approach (named INI-VPINN). INI-VPINN naturally incorporates Neumann boundary and interface conditions into the variational formulation. It removes the need for additional loss terms or multiple subdomain networks. This framework employs compact support weighting functions and integration by parts to implicitly impose flux and continuity constraints. In this way, it implicitly ensures physical consistency across material boundaries. The proposed method is tested on Poisson and Laplace problems with sharp interfaces and complex geometries. Results show that, compared with several other Physics Informed Neural Networks-based formulations, the INI-VPINN consistently achieves higher accuracy, smoother and faster convergence. The proposed framework provides a general approach for solving multimaterial problems with complex geometries and mixed Neumann-Dirichlet boundary conditions using neural networks. The implementation is publicly available in a GitHub repository.

The optimised tokamak-stellarator hybrid concept (Henneberg and Plunk 2024) has the potential to combine tokamak and stellarator advantages to achieve magnetically confined fusion. These compact quasi-axisymmetric designs can have a low aspect ratio and large plasma volume, good particle confinement, and relatively simple coils. Previous work showed that such magnetic configurations can in principle be reproduced by a single type of non-planar "banana coil" alongside the conventional tokamak coilset (Henneberg and Plunk 2025). In this work, we optimise banana coils while also considering engineering constraints beyond simple geometric measures. We quantify the characteristic geometries of force-optimised banana coils and the magnetic fields they generate, and analyse the mechanisms by which forces may be reduced through optimisation.

Ion exchange membranes are useful for a wide range of applications, including water desalination, fuel cells, and aqueous batteries. Accordingly, a variety of models for ion exchange membranes has been proposed, each emphasizing different aspects that govern their static and dynamic properties. By reviewing these models, we identify key physical contributions and beneficial modeling strategies. Based on these insights, we derive a thermodynamically consistent model by combining mass-action site occupation with mean-field electrostatic interactions along the polymer backbone. In this derivation, we explicitly account for multicomponent electrolytes at elevated concentrations. The parameters of the resulting model relate closely to those of other models, but gain consistency and interpretability through the underlying derivation. A discussion of the model parameters highlights redundancies and linkages between quantities that are commonly treated independently. Comparison to experimental data shows that both static and dynamic membrane properties are reproduced with good accuracy by the presented model. This makes it a promising basis for theory-driven membrane optimization and supports the tailored design of ion exchange membranes for various technologies.

A fourth-order conservative adaptive multiresolution average-interpolating wavelet upwind scheme is proposed for compressible flows governed by hyperbolic conservation laws. A family of asymmetric average-interpolating wavelets with upwind properties is constructed for conservative finite volume discretization, while symmetric average-interpolating wavelets are employed for multiresolution decomposition and reconstruction of physical variables in the adaptive procedure. Since both the conservative discretization and the adaptive multiresolution representation are constructed from cell-average quantities, the proposed scheme preserves strict conservation during both numerical evolution and adaptive cell redistribution. Unlike hybrid adaptive wavelet methods that use wavelets mainly for data compression and mesh adaptation, the present adaptive wavelet upwind scheme utilizes average-interpolating wavelet multiresolution approximation to reconstruct the interface values directly for numerical flux evaluation, thereby avoiding additional ghost-cell marking and reconstruction near coarse–fine mesh interfaces. The boundary variation diminishing reconstruction is incorporated at the finest resolution level to achieve non-oscillatory shock-capturing capability. Numerical tests demonstrate that the proposed scheme achieves the expected fourth-order accuracy, maintains conservation errors close to machine precision, and controls numerical errors around the prescribed threshold. The proposed method also sharply captures shock waves and contact discontinuities without spurious oscillations and resolves multiscale smooth structures through a sparse adaptive representation. These results indicate that the proposed scheme provides an efficient, conservative, and reliable approach for high-resolution simulations of compressible flows.

Improving the resolution of telescope systems will provide the opportunity to study new physical phenomena in previously unobserved environments. Spatial mode de-multiplexing (SPADE) based imaging is a promising and rapidly evolving technique for pushing the resolution of optical telescopes beyond the diffraction limit. A key application of this technique is for near-optimal hypothesis testing for the presence of secondary and extended sources in the sub-diffraction regime. We present the first demonstration of a binary-SPADE based hypothesis testing instrument deployed on-sky. In our proof-of-principle experiment, based on mode demultiplexing with a double clad fiber coupler, we demonstrate detection of a binary star system separated below the diffraction limit. We perform measurements in the photon-starved regime where no image can be formed by traditional direct imaging. We find the scaling of the system's type II error rate (the “binary source miss" chance) was heavily limited by unbalanced loss in our double-clad fiber coupler when compared to the idealized quantum limits. Despite this the evaluated type II error is always lower than a perfect direct imaging measurement. We expect that if this instrument is scaled to larger aperture telescope systems the effects of atmospheric turbulence will further degrade this system's performance.

We study equilibrium clustering in a finite one-dimensional lattice gas of $L$ sites with periodic boundary conditions, as a minimal model for adsorption and binding on small ring-like substrates. Using a grand-canonical formulation with nearest-neighbor coupling, we derive exact finite-size expressions for the mean occupancy, the mean number of domain walls, and the mean number of clusters. Building on exact $k$-site correlation functions, we further derive expressions for the mean number of clusters of size $k$ and for two complementary size statistics: the cluster-size distribution, and the site-weighted cluster-size distribution. These observables characterize how spatial organization changes across attractive (cooperative) and repulsive (anticooperative) interactions, and highlight finite-size and parity-dependent effects of the underlying lattice, the latter being particularly pronounced near half filling in small systems. To access larger lattices without enumerating all $2^L$ microstates, we also develop a cluster-based combinatorial formulation in which configurations are classified by cluster counts and sizes, reducing the effective state space to a set whose size scales with integer partitions, $\approx e^{\sqrt{L}}$, rather than with $\approx e^{L}$. Taken together, our results provide exact benchmarks for finite periodic systems and suggest experimentally relevant cluster observables that complement occupancy-based measures of cooperativity, with particular relevance for binding on ring-like substrates for biological assemblies.

We derive an effective orbital-lattice model with quantum $S=1$ degrees of freedom for transition-metal compounds, providing a microscopic understanding of cluster formation driven by the cooperative interplay of spin, orbital, and lattice degrees of freedom. Motivated by the trimerized phases observed in LiVS$_2$ and LiVO$_2$, we consider a triangular-lattice three-orbital system with two electrons per site occupying the threefold-degenerate $t_{2g}$ manifold. Starting from a multiorbital Kanamori-Hubbard Hamiltonian, we project the low-energy sector onto the local $S=1$ triplet manifold, in which two electrons occupy different orbitals according to Hund's coupling. The resulting effective model exhibits exchange networks whose geometry is determined by the orbital configuration. However, the orbital-driven exchange interactions alone do not stabilize the experimentally observed trimer phase. We find that by incorporating ionic lattice displacements that modulate transfer integrals and induce bond-dependent exchange couplings on shortened and elongated bonds, the phase competition is qualitatively altered, leading to the robust stabilization of a trimerized ground state within a fully quantum-mechanical framework. We further show that a simplified orbital-lattice model, in which the spin-exchange energy is replaced by effective bond energies, faithfully reproduces the essential ground-state properties of the microscopic model. This reduced description enables large-scale finite-temperature simulations and reveals a rich sequence of thermal phase transitions, including first-order, second-order, and Kosterlitz-Thouless transitions into distinct spin-, orbital-, and lattice-ordered phases.

We extend the potential-based formulation of self-dual and duality-symmetric gauge theories developed in \cite{Chakrabarti:2025gyt} to curved spacetime. The gravitational coupling is encoded by the same metric-dependent linear map of Sen's formalism. The parent new action retains the higher-form $h$-gauge symmetry, and different gauge choices reproduce either Sen's flux-based formulation or a potential-based description in which the shadow sector decouples directly at the level of the action. We also derive explicit form of the gravitational map for the toroidally reduced theory in $D=4n$ and for a polyform formulation in arbitrary dimension, both of which apply equally to ours as well as Sen's formalism. As non-trivial checks, we reproduce the Álvarez-Gaumé-Witten gravitational anomalies of the two-dimensional chiral boson without fermionisation, and of the ten-dimensional self-dual five-form. Finally, on $\mathrm{AdS}_5\times S^5$ the physical on-shell action yields the non-vanishing boundary contribution required by holography without adding a boundary term by hand. We also establish a local on-shell identification with the textbook supergravity RR four-form potential.

Many of the most powerful and elegant models in physics are grounded in symmetries. In electrodynamics, for example, geometric symmetries govern the observable effects of light-matter interactions. However, for man-made objects, exact symmetries are rarely met and tiny deviations are common. Nonetheless, even approximate symmetries keep many symmetry-derived rules effectively intact. However, as we will show here, this is not universally true. We demonstrate that an incremental violation of the symmetry of a carefully designed system can produce an optical response maximally different from the unbroken symmetry case. To do so, we exploit symmetry-protected quasi-bound states in the continuum (qBICs). Specifically, we design a four-fold rotationally symmetric metasurface composed of nearly dual-symmetric meta-atoms that supports a pair of spectrally aligned electric and magnetic qBICs. At normal incidence, symmetry forbids helicity-preserving reflection. However, for arbitrarily small deviations from normal incidence, the strong resonant enhancement associated with the qBICs overcomes the near-symmetry suppression and enables perfect helicity-preserving reflection. This rapidly emerging violation of symmetry-rules reveals a fundamental intricacy when it comes to treating near-symmetric systems. At the same time, our work opens the door to novel applications in metrology and sensing.

The widely known Yin-Yang symbol (Taijitu) is based on nested circles of different radii whose areas are colored black and white such that the interface traces an $\mathcal{S}$-shaped curve. We address the question of how this symbol can be related to physical phenomena such as daytime and nighttime duration and the annual seasons. Using a simple dynamic model of daytime duration, we introduce the excess daytime fraction and reconstruct the symbol using the latitude of Zurich. In particular, we explain how the black and white areas are linked to the stability of Yin or Yang predominance. We further demonstrate that the Golden and Silver Ratios found in the geometry of the symbol carry meaning with respect to the Gregorian calendar. Finally, we construct an alternative Yin-Yang symbol using logarithmic spirals with the Golden Ratio as the growth parameter. The didactical quantitative derivation of the Yin-Yang symbol and its grounding in real-world observations can be regarded as a novel perspective on this iconic pattern.

The crack closure effect of a low-alloyed steel subjected to low-cycle fatigue loading has been characterized at two different imposed strain amplitudes. Two techniques (non-continuous DIC H-DIC and DCPD) have been employed in this aim, leading to similar conclusions. Thus, it is shown that the crack does not remain completely closed during a part of the compressive portion of the fatigue cycle for both applied loadings. The crack opening strains, combined with the cyclic stress-strain curve of the material allowed to determine an equivalent cyclic opening stress. This confirmed that the crack opening stresses decrease with the applied maximum stress when subjected to tension-compression loading, as commonly found in the literature.

We present a detailed timing analysis of the transition from the hard-intermediate state (HIMS) to the soft-intermediate state (SIMS) in MAXI J1820+070 using NICER observations. This transition is marked by a sharp drop of the broadband noise across both the soft and hard X-ray bands, the disappearance of the Type-C quasi-periodic oscillation (QPO), the quenching of the steady, optically thick, compact jet, the appearance of a Type-B QPO, and the detection of discrete, optically thin, radio ejections. For the first time, we detect a QPO at 3.5-5.9 Hz in the 2-12 keV power density spectrum of MAXI J1820+070 roughly half a day before the transition, which appears to evolve smoothly into the Type-B QPO observed immediately after the transition. The location of this additional QPO component in the broadband rms vs. QPO frequency plot is consistent with that of the Type-B QPOs in GX 339-4 and GRO J1655-40, suggesting a possible connection between this additional QPO in the HIMS and the Type-B QPO in the SIMS. This result, together with recent findings in Swift J1727.8-1613, suggests that QPOs with these characteristics can emerge prior to the HIMS-to-SIMS transition and are not confined exclusively to the SIMS. If this additional QPO feature is the precursor of the Type-B QPO in the SIMS, its presence before the transition, whereas the bright discrete, optically thin, radio ejections appear at the transition, would imply that there may be no direct physical connection between the Type-B QPO and the discrete radio ejections. Our results further suggest a link between the disappearance of the Type-C QPO, the drop of the broadband noise, and the emergence of discrete radio ejections at the HIMS-to-SIMS transition. We speculate that the simultaneous presence of such a QPO, non-harmonically related to the Type-C QPO in the HIMS, could be compatible with a spine-sheath outflow structure.

The rapid neutron-capture process (r-process) is widely regarded as the dominant mechanism responsible for the synthesis of heavy elements in the universe, yet its astrophysical sites remain an open question. Recent studies suggest that the high-entropy, rapidly expanding baryonic material ejected by magnetar giant flares may provide favorable conditions for r-process nucleosynthesis, while the late-time gamma-ray emission observed from the magnetar SGR 1806-20 offers direct observational support for this scenario. In this work, we perform nuclear reaction network simulations to investigate the nucleosynthesis yields of magnetar giant flares and to characterize the associated nuclear gamma-ray line emission arising from the radioactive decay of heavy nuclei. The nuclei synthesized in magnetar giant flares are found to be mainly distributed near the first and second r-process abundance peaks. Owing to this nuclide composition, the gamma-ray opacity is found to be strongly energy-dependent with the opacity in the keV band exceeding that in the MeV band by approximately three orders of magnitude. The nuclear gamma-ray emission is dominated by MeV photons at early times and gradually extends toward the sub-MeV and keV bands as time progresses, thereby offering a diagnostic of heavy element enrichment in the ejecta. The gamma-ray spectrum exhibits a peak near 1 MeV with major contributions from $^{88}$Kr and $^{92}$Sr, whose radioactive decays produce several bright gamma-ray lines with fluxes exceeding $\sim10^{-8}$ erg cm$^{-2}$ s$^{-1}$, making them the most promising lines for detection by MeV gamma-ray detectors. Because magnetar giant flares occur in the Galaxy at a rate roughly three orders of magnitude higher than neutron star mergers and their gamma-ray lines are accessible to current MeV instruments, they offer new and valuable science opportunities for MeV gamma-ray astronomy.

We present a robust assessment of cosmological constraints on the sum of neutrino masses ($\sum m_\nu$) when relaxing the standard assumption of purely adiabatic primordial initial conditions. Allowing for a neutrino density isocurvature (NDI) component alongside the adiabatic mode, we analyse the latest CMB-SPA combination (Planck 2018, ACT DR6, and SPT-3G), DESI DR2 baryon acoustic oscillation data, and the DES Year 5 supernova sample. Within the $\Lambda$CDM model, the 95% upper limit weakens only marginally from $\sum m_\nu < 0.052$ eV (purely adiabatic) to $< 0.057$ eV (including NDI), with the NDI amplitude consistent with zero. In the CPL dynamical dark energy model, the adiabatic limit is $< 0.111$ eV, shifting to $< 0.115$ eV with NDI, yet the isocurvature mode remains undetected. While these limits are robust against the inclusion of isocurvature perturbations, they are highly sensitive to both the assumed dark energy equation of state and the prior lower bound on $\sum m_\nu$. Notably, the adiabatic $\Lambda$CDM limit of $0.052$ eV lies below the minimum sum required by the normal neutrino mass hierarchy ($0.05878$ eV), indicating that this bound is an artifact of the statistical prior extending to zero. Imposing a physically motivated hierarchy-informed prior raises the limit to $< 0.092$ eV. Our results demonstrate that current data show no evidence for NDI modes and that the inferred neutrino mass upper limit is robust against this extension, but a definitive, model-independent bound requires addressing prior dependencies and dark energy uncertainties. This work provides the first joint constraint on $\sum m_\nu$ and NDI using the full CMB-SPA+DESI DR2+DES dataset.

We calculate the one-loop corrections to the quark TMD in the projectile light-cone gauge using the background field formalism, with the Mandelstam-Leibbrandt (ML) prescription for the extra singularity present in the light-cone gauge propagator. We use the pure rapidity regulator for rapidity divergences. The Collins-Soper-Sterman (CSS) evolution equations are obtained after one-loop renormalization of the quark TMD in this gauge. We discuss how the structure of the rapidity divergences and the double-logarithmic contributions to the CSS resummation compares with the analogous calculation performed in the target light-cone gauge, and discuss the implications for the connection between the TMD factorization and Color Glass Condensate frameworks.

The strong electric fields from tightly-focused and ultrashort laser beams have always been discussed as a way to accelerate charged particles without any need for a medium or external cavity. Radially-polarized light is one way to do this, motivated by the emergence of longitudinal electrical fields with tight focusing. However, the laser pulse will generally quickly overtake the electrons under its influence, and every-other half-cycle will decelerate the electrons in effect partially reversing the acceleration. In this work we present the effect of optical aberrations, primarily spherical aberration, and how despite their purely spatial nature they can significantly optimize the net acceleration, and advantageously allow for longer pulses to drive this optical field-based process. We discuss the optical physics responsible for this increase in performance and find optimal aberration profiles using a stochastic algorithm.

We present a weak gravitational lensing mass calibration of 19 Planck Sunyaev–Zel'dovich (SZ) selected galaxy clusters using shape measurements from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Year 3 shape catalog. We measure the stacked weak-lensing signal $\Delta\Sigma(R)$ using per-cluster lensing weights that match the measurement pipeline's stacking scheme, and construct an analytical covariance matrix that includes shape noise and projected large-scale structure contributions. Our primary constraint on the SZ mass bias comes from a forward-modeling approach that integrates over the halo mass function while accounting for the Planck SZ selection function, Eddington bias from log-normal scatter in the SZ mass proxy, and cluster miscentering. Fitting four free parameters, the log mass bias $\ln(1-b)$, the miscentered fraction $f_{\rm mis}$, the offset scale $\sigma_{\rm off}$, and the SZ scatter $\sigma_{\ln M}$, over the radial range $0.5$–$5.0\,h^{-1}\,\mathrm{Mpc}$, we obtain $1-b = 0.73^{+0.10}_{-0.11}$ with $\chi^2/\mathrm{dof} = 5.2/5$ at an effective redshift $z_{\rm eff}\simeq 0.24$. This measurement is consistent with recent weak-lensing calibrations of SZ-selected clusters and supports the picture that significant mass bias corrections are required to reconcile cluster abundance measurements with primary cosmic microwave background constraints on cosmological parameters.

Extending the holographic principle to de Sitter (dS) spacetimes remains one of the most vital open frontiers in quantum gravity, where a microscopic, bottom-up tensor-network framework that relates boundary quantum data to emergent de Sitter spacetime is still lacking. In this work, we first show the emergence of de Sitter spacetime from boundary entanglement by formulating a non-unitary continuous multi-scale entanglement renormalization ansatz (cMERA) for a concrete non-Hermitian critical fermion chain. Within this emergent spacetime, we analyze the associated geodesics and show that they act as extremal Ryu-Takayanagi (RT) surfaces undergoing a smooth timelike-to-null transition. Remarkably, we demonstrate that this continuum trajectory suggests a distinct tensor-network architecture in which the bond-counting contribution naturally truncates at the discrete timelike-to-null transition toward the deep infrared. In the resulting architecture, the null ray along the horizon is represented by zero-cost links, since the associated cut severs no tensor legs. This network structure successfully reproduces the logarithmic scaling of non-unitary critical entanglement entropy, offering a bond-counting picture for the de Sitter RT formula. Our results provide the long-sought dS/(c)MERA correspondence at the level of both emergent spacetime and discrete holographic entanglement.

We present the first X-ray polarization measurements of the PSR/PWN complex within SNR Kes 75. Two $\rm {\sim}\,500\,ks$ IXPE observations were conducted in October/November 2024 and April 2025. The second observation yields a significant phase-average 2-8 keV polarization degree $\rm PD = 9.9\% \pm 2.5\%$ at $\rm PA=36.8^\circ \pm 7.3^\circ$, implying a toroidal field aligned with the PWN symmetry axis. The first epoch, however, has only a polarization upper limit. During this epoch, an additional pulsed component is visible at $\Delta \phi \approx 0.5$, detected at ${\sim}\,3.7\sigma$. An unbinned phase-resolved analysis reveals a high-PD rotating vector model PA sweep at the ${\sim}\,99.5\%$ confidence level, with angles fixed at those inferred from the PWN morphology; this can explain the loss of phase-average polarization. Additional observations are needed to pin down the nature of the anomalous pulse.

Recent studies have revealed a notable timing feature in several black hole X-ray binaries (BHXBs) during the soft-to-hard transition at the outburst decay. Within a narrow frequency range, the phase lags between high- and low-energy X-ray light curves exhibit a sudden increase, accompanied by a drop in the coherence function. These narrow features have been associated with a quasi-periodic oscillation (QPO) appearing only in the imaginary part of the cross spectrum (CS). This QPO remains undetected in the power density spectrum (PDS) and is known as imaginary QPO. Motivated by these results, we analyse five years of NICER observations of the BHXB AT2019wey during its low-hard state (LHS) and hard-intermediate state (HIMS). We find an imaginary QPO in the CS of AT2019wey, with similar characteristics as those found in other BHXBs, making AT2019wey the fifth BHXB in which such QPOs have been found. As the source hardens, the frequency of the imaginary QPO drops from $\sim$ 5 Hz to $\sim$ 1 Hz, while its phase lag rises from $\sim$ 0.3 rad to $\sim$ 0.7 rad during the HIMS and from $\sim$ 0.5 rad to $\sim$ 0.6 rad during the LHS. During the HIMS, the phase-lag energy spectrum of the imaginary QPO shows a typical U-shaped profile, while the shape changes in the LHS. The rms spectrum of the imaginary QPO rises below $\sim$ 2 keV, peaks at around $\sim$ 2 keV and decreases at higher energies, which may be associated with the presence of a relatively cool corona. We compare the properties of the imaginary QPO with those of the type-B and C QPOs in BHXBs and find a tentative connection to type-C QPOs. Combining the imaginary QPOs detected in AT2019wey with those reported in other sources, we find a systematic increase of QPO phase lags with QPO frequency. However, we cannot conclude whether the phase lags of imaginary QPOs exhibit the inclination dependence previously observed in type-C QPOs.

Higher-order topological phases provide robust corner modes, but their use requires controllable creation, isolation, and transfer of individual modes and their superpositions. Here we demonstrate, using the two-dimensional Benalcazar-Bernevig-Hughes model as an example, that subchiral symmetry provides a general control principle for manipulating topological corner modes. The conventional chiral symmetry decomposes into four subchiral symmetries, each associated with one zero-energy corner mode. By selectively breaking these subsymmetries with controlled intercell hoppings, we reduce the fourfold corner-state manifold step by step to single isolated modes. We further design adiabatic protocols that transfer either a single corner state or a superposition of two corner states between selected corners, while preserving the relative phase in the latter case. Both numerical simulations and IBM quantum-processor implementations show that the proposed protocols can be executed with high fidelity, establishing subchiral symmetry as a route to programmable higher-order topological state manipulation.

Repeated thermal optimization must enforce local temperature limits while solving linear systems whose unknowns change as the constraints activate and release. Zero-extension transfer, which carries correction vectors from one inactive set to the next, assigns zero to the degrees of freedom that a constraint has just released. We construct a spectral reference on the full state domain before optimization and restrict it directly to the current inactive set at each update, so the coarse correction receives reference entries at every released degree of freedom. Spatial and space–time constructions cover steady and transient conjugate heat transfer with prescribed fluid velocity, and a restriction-residual identity yields alignment and deflated-convergence bounds under stated rank and spectral-separation assumptions. In the primary comparisons on body-fitted transformer and bore-in-block models, the resulting reference deflation reduces the median complete optimization time relative to the tested iterative alternatives on the same graphics processor. Complete time includes reference construction. Matched transfer experiments connect the improvement to the reference information retained at released constraints. The lower cost persists on refined Cartesian and bore-in-block meshes, on space–time trajectories and over sequences with changing loads and query order; the margin narrows, and eventually reverses, as a trajectory gains time slabs at a fixed horizon. Full-domain reference retention thus supports repeated evaluation of temperature-constrained heat inputs in engineering thermal models.

In practical applications, compressible particle-laden flows in moving geometries involve complex, non-linear, and multi-scale inter-actions with turbulent structures. Resolving these dynamics numerically requires careful algorithmic treatment to accurately predict particle trajectories. This work presents a high-fidelity Euler-Lagrange framework that couples a high-order discontinuous Galerkin spectral element method for the continuous phase with a Lagrangian point-particle tracking scheme. To manage moving and deforming domains, the framework integrates two distinct mesh movement strategies: the arbitrary Lagrangian-Eulerian method for general mesh deformations such as time-resolved particle-induced surface deformations and its special interface case, the sliding mesh approach, uniquely suited for rigid rotational or translational movements. A primary focus is placed on tightly coupling the arbitrary Lagrangian-Eulerian formulation into the temporal evolution step by utilizing radial basis function morphing to capture the non-linear feedback loop between evolving surface topologies and the continuous phase. Concurrently, the framework ensures time- and high-order accurate coupling of the mesh movement algorithms with the dispersed phase. In particular, the proposed algorithm resolves the sliding mesh tracking problem by enforcing strict spatial and temporal accuracy as Lagrangian particles cross non-conforming grid interfaces between adjacent moving zones. The algorithms are rigorously validated against multiple benchmarks and subsequently applied to two challenging compressor rotor applications: the first focusing on solid-particle erosion, and the second featuring an upstream cylindrical wake generator.

According to the recent Wealth Thermalization Hypothesis (WTH) the wealth inequality in the world is described by the Rayleigh-Jeans (RJ) thermal distribution of interacting agents in a society with social stratification. In this concept, the wealth layers of society are associated with energy levels from a nonlinear dynamical system conserving two integrals of motion being total energy and probability norm. This leads to RJ condensation and the formation of a huge poverty phase of low wealth and a tiny oligarchic phase that captures a main part of total society wealth. This RJ phenomenon has similarities with self cleaning in multimode optical fibers and constraint driven condensation in various physical systems. We analyze real Lorenz and Pareto curves for wealth of households in countries and the world, Gross Domestic Product of countries, market capitalization of companies at stock exchange of Hong Kong, Shanghai, London, bitcoin transactions, world trade between countries and show that the WTH theory gives a good description of these curves. On the basis of this comparison we argue that the RJ thermal distribution provides a universal description of wealth inequality in the world.

We study the path integral quantisation of a scalar field on a generic noncommutative deformation of Minkowski space, built as a quantum homogeneous space of a deformed Poincaré group. We show that the procedure depends on the choice of noncommutative functional derivative, and we isolate two natural choices, distinguished by the space in which the Leibniz rule remains undeformed. The first, which carries the undeformed statistics, reproduces the standard scheme and yields n-point functions that break the deformed Poincaré covariance and exhibit the UV/IR mixing of [8]. The second, which adapts the functional calculus to the braided statistics of the fields in the spirit of [20], yields covariant n-point functions free of this mixing. We trace both the covariance breaking and the mixing to a single source, the intertwining of external and loop momenta in the non-planar contributions, and conclude that, in the models considered, the UV/IR mixing of [8] is an artifact of a quantisation that breaks the deformed symmetry rather than a feature of noncommutativity itself. We further disentangle covariance, fixed by the quantisation scheme, from finiteness, fixed independently by the propagator, and illustrate the formalism on the T-Minkowski models, the Euclidean three-dimensional quantum gravity model, and the quantum two-sphere.

We review the mechanism by which loops of matter fields contribute to the graviton self-energy during de Sitter inflation. This self-energy is used to quantum-correct the linearized Einstein equations. A Green's function method is employed to obtain exact 1-loop corrections to the plane wave mode functions of gravitational radiation, subject to the usual ambiguity in the initial state. Conformally coupled matter, which does not experience inflationary particle production, makes only a logarithmic enhancement of the rate at which the imaginary part of the mode function goes to zero after horizon crossing. These corrections can be understood, and even summed up, using a variant of the renormalization group. However, massless, minimally coupled scalars, which experience massive inflationary particle production, induce a much stronger enhancement of the rate at which the real part of the mode function approaches a constant. One interpretation of this effect is as a shift of the inflationary Hubble parameter.

Random states play a foundational role in different branches of modern quantum science. In this work, we study a recently proposed random state ensemble induced from von Neumann entropy through the Bogoliubov-Kubo-Mori (BKM) metric. In particular, we derive an exact yet explicit formula of average entanglement entropy over BKM ensemble. In obtaining the formula, we only make use of properties of normalization constant of the ensemble in the absence of its correlation kernel, contrary to average entropy computation of other ensembles. This new framework paves the way for calculating higher-order cumulants of BKM ensemble beyond the average.