Papers

Hypothesis: The presence of a polymer in a droplet plays a significant role in contact line dynamics following droplet impact onto a solid substrate. We hypothesize that the polymer concentration modifies the post-impact wetting behavior. Experiments: In our experiments, we visualize the contact line dynamics of polymer-laden droplets impacting onto a flat sapphire substrate. Using high-speed side and bottom-view imaging, we analyze droplet spreading and retraction over a broad range of Weber numbers (10-700) and polyethylene oxide (PEO) concentrations (0-400 ppm). To systematically characterize the contact line dynamics, we used image processing to quantify the maximum spreading factor, receding velocity, and receding contact angle. Findings: The effect of the polymer on the spreading phase is negligible, as the maximum spreading factor follows classical inertial-capillary scaling, showing that polymer addition produces no measurable change in the maximum spreading under the present conditions. In contrast, the receding phase exhibits a strong dependence on the polymer concentration: increasing the PEO content leads to a significant reduction in the receding velocity and dynamic contact angle, while the spreading dynamics remain water-like. These results are consistent with an additional resistance associated with the receding contact-line region, while the microscopic origin of this resistance remains unresolved. To capture this behavior, we propose a scaling model that relates the normalized retraction velocity to the receding contact angle and impact conditions. The model effectively collapses the experimental data across all investigated concentrations and Weber numbers.

Altermagnets combine compensated magnetic order with momentum-dependent spin splitting, providing spin-polarized electronic states without a net magnetization. Metallic $d$-wave altermagnets are particularly promising because they can support time-reversal-odd spin currents in linear response, yet experimentally validated bulk realizations remain scarce. Here, we use high-throughput density-functional-theory data to relate the magnitude of altermagnetic band splitting to compositional, structural, and DFT-derived magnetic descriptors. An interpretable gradient-boosted model identifies two candidate-prioritization criteria: compact unit cells and magnetic sublattices that sustain sizable local moments. Guided by these trends, we screen tetragonal $A_2XY$ Heusler compounds in space group $P4/mmm$. Among 307 structures, symmetry identifies 169 altermagnetic arrangements, of which 157 remain metallic altermagnets in DFT. Sixteen realize an altermagnetic collinear ground state, of which 15 are dynamically and mechanically stable and 10 also lie on or below the calculated thermodynamic hull. Six candidates exceed the CrSb splitting obtained under the same computational protocol, led by Co$_2$AlSc ($\Delta_\mathrm{max}=2.24$~eV) and Fe$_2$AlGe ($2.07$~eV). A Julliere-model estimate gives a tunneling magnetoresistance of up to $203\%$ at the Fermi level for Co$_2$AlSc. These results identify a chemically tunable family of metallic $d$-wave altermagnets and demonstrate how interpretable machine learning can guide targeted first-principles searches.

The conventional grid-based approach used for storing and evaluating ordinary parton distribution functions (PDFs), as implemented for example in LHAPDF, becomes much more demanding for double parton distributions (DPDs). While a collinear PDF depends on one longitudinal momentum fraction and one factorization scale, an unequal-scale DPD depends on two momentum fractions and two independent factorization scales for every parton-flavor combination. Direct tabulation of this four-dimensional dependence therefore requires considerably more disk space and runtime memory. In this work, we develop a hybrid neural-network method for compact storage and fast evaluation of $y$-independent unequal-scale DPDs. The network compresses the reference grid by reproducing most of its values, while the original values are stored only at points where the neural prediction does not reach the required accuracy. In this way, the neural model captures the bulk of the grid, with a sparse table covering the remaining difficult points. For the GS09-based dense grid studied here, only $0.694\%$ of the active flavor values need to be stored in the sparse table. Compared with the compressed dense grid, the hybrid method reduces the disk footprint by a factor of $10.6$, peak memory use by a factor of $31.6$, and initialization time by about a factor of $10.6$. At the same time, its evaluation throughput reaches about $85\%$ of the dense-grid rate. The method is implemented natively in C++ within PDFxTMDLib and requires neither Python nor PyTorch at runtime. The hybrid scheme therefore offers much lower memory and startup costs, with only a modest reduction in evaluation speed.

Turbulence plays a critical role in the evolution of H II regions, yet recovering its statistical properties from observations remains challenging, particularly when centroid velocities are affected by instrumental noise. We present a methodology for recovering the second-order velocity structure function from intermediate spectral-resolution integral-field spectroscopy and apply it to VLT MUSE observations of the Orion Nebula. The plane-of-sky velocity field is characterized by fitting a simple parametric model to the observed structure function, allowing the turbulent velocity variance, correlation length, power-law slope, and noise level to be estimated simultaneously. We investigate the trade-off between instrumental noise and spatial resolution through systematic spatial binning and validate the recovered turbulent parameters against previously analyzed high spectral-resolution KPNO echelle observations. We find that spatial binning improves the empirical structure function provided that the bin size remains below approximately $0.05 r_0$, while the fitted model recovers consistent turbulent parameters across all binning levels. The MUSE results agree closely with those from the KPNO data, demonstrating that intermediate-resolution integral-field spectroscopy can recover reliable turbulence statistics despite its relatively poor velocity resolution. We also identify modest but systematic differences in the turbulent properties traced by emission lines of different ionization potential, reflecting the geometric and ionization structure of the nebula. Our methodology provides a robust framework for extracting turbulent velocity statistics from noisy integral-field spectroscopic observations.

The present work proposes a novel Lagrangian dynamical theory for the velocity gradient tensor (VGT) $\bm{A}\equiv\bm{\nabla}\bm{u}$, formulated on the basis of its canonical real Schur form.... The proposed theory is expected to be useful for understanding the evolution dynamics of the VGT rotational invariants and intrinsic vorticity modes in a variety of vortical flows.

Squashed entanglement is an additive bipartite entanglement measure. For a state $\rho_{AB}$, it is defined as the infimum of all conditional mutual information $I(A:B\mid E)$ evaluated on extensions $\rho_{ABE}$. It was unresolved whether there is a bound on the dimension of the conditioning system, $E$, needed for this optimization. We show that no such cardinality bound is possible. Explicitly, we find that the squashed entanglement for a partially dephased Bell pair on two qubits cannot be achieved for any finite dimensional extension $E$. A crucial ingredient of this proof is a direct sum step, which in a sense, symmetrizes two purifying systems while keeping the conditional mutual information unchanged and the coherence no worse.

Background and purpose: Dose Volume Histogram (DVH) indices remain the main dose-effect metrics for toxicity prediction, but they depend on how contours were made. Inter-observer variability (IOV) is unavoidable and clinically accepted, so the question for automatic segmentation addresses equivalence: do automatic contours produce DVH errors comparable to human IOV, and can these indices predict patient-reported toxicity? Materials and methods: In 429 patients treated with iodine-125 LDR-BT monotherapy, indices from expert manual delineation were compared with a deterministic and a Bayesian nnUNet on a fixed dose distribution, by two-one-sided tests against margins set from CT contouring IOV. Logistic regression gave patient-level thresholds at 90\,% probability of equivalence. In 380 patients, eleven DHV indices were added to a clinical baseline predicting change in International Prostate Symptom Score (IPSS) at six horizons over 5 years, with nested cross-validation, bootstrap intervals and corrected $t$-tests. Results: All cohort-level comparisons of DVH indices were declared equivalent, with no interval consuming more than 44\,% of its equivalence margin. Individual agreement was weaker with equivalence rates from 55.9\,% to 89.7\,%, depending on the DVH index and automatic segmentation model. Thresholds ranged from 0.864 to 0.958 Dice. No DVH block improved IPSS prediction at any horizon. The largest improvement declared by bootstrap intervals was 0.12 IPSS points, well below the minimal clinically important difference, across all learners. Conclusion: Automatic contours matched expert dosimetry within human IOV on average, but individual equivalence requires a volume dependent quality metric threshold definition. Our DVH indices panel provides no significant predictive power for IPSS, irrespective of segmentation source and horizon.

Near-threshold $J/\psi$ production is sensitive to scalar and traceless spin-two chromoelectric structure in the nucleon. We study whether present integrated and differential measurements can separate these contributions and constrain their relative phase. The GlueX integrated cross section does not determine a unique scalar contribution: production and elastic scalar terms are strongly correlated, while flatter non-forward profiles drive the fitted scalar strength to zero. In a 38-point large-$\xi$ differential analysis, the scalar term is collinear at fixed kinematics with the retained C-GFF direction, leading to nearly degenerate solutions and strong profile dependence. A common single-channel final-state factor preserves the scalar-spin-two relative phase, whereas coupled-channel rescattering can modify it. Its determination therefore requires polarization-sensitive observables together with process-specific complex helicity amplitudes.

This paper focuses on the quasinormal modes of a charged loop quantum black hole that preserves diffeomorphism covariance. We have examined the influence of the black hole charge Q and the quantum parameter {\zeta} on the quasinormal modes of neutral massless scalar field, neutral massive scalar field and massless charged scalar field. For the neutral massless scalar field, our work shows that as in the case of the corresponding uncharged black hole that has been examined by other authors, the turn-on of {\zeta} renders the real parts of the overtones a tendency to approach zero. Increasing Q can inhibit this tendency. For the neutral massive scalar field, when Q and {\zeta} are turned on, arbitrarily long-lived modes still exist. We also observe a non-monotonic behavior in the fundamental mode as the field mass varies for the first time. For the massless charged scalar field, our study shows that the asymptotic behaviors of the quasinormal modes of charged field and uncharged field are different as the extremal charge Qis approached. We also find that the real part of the fundamental mode tends to zero at some critical value eQfor a charged scalar field with field charge e.

We report the full sub-Kelvin characterization of the BAUSCIA high-frequency gravitational-wave antenna, based on cryogenic Bulk Acoustic Wave resonators coupled to SQUID amplifiers for signal readout. The end-to-end performance analysis of the experimental chain, from the sensor to the DAQ, demonstrates that resonator thermal noise dominates the noise budget at temperatures well below 100 mK for all investigated modes. The measurements also verify proper thermalization of the resonator modes and rule out significant SQUID backaction on the strain sensor, confirming that the system operates close to its fundamental limits with comparable resonator thermal noise and additive readout noise. From this characterization, we derive projected strain sensitivities in the range $1$–$10 \times 10^{-21}\,\mathrm{Hz}^{-1/2}$ for several resonant frequencies between 5 and 20~MHz, with corresponding fundamental thermal limits in the range $0.7$–$1.4 \times 10^{-21}\,\mathrm{Hz}^{-1/2}$. Equipped with three resonators and synchronized to an absolute global time reference for multi-site coincidences, the antenna is ready to enter a pilot science run.

Generative augmentation is a standard remedy for label scarcity in physiological signal classification, but existing quantum generative models start from uninformative noise, ignoring class structure that is already available. We propose Conditional Quantum Flow Matching (CQFM): a single 306-parameter circuit, conditioned on both flow time and class label, transports a compact class-conditional prior toward the target distribution. Quantum flow matching as published is unconditional, so this is to our knowledge the first conditional one, and the first EEG augmentation on a parameterized quantum circuit. A nonnegative spectral embedding removes the need for tomography at readout. On BCI Competition IV-2a, starting from a prior rather than noise is worth $+5.1$ accuracy points over QuDDPM (9/9 subjects), though at that operating point a class-conditional Gaussian matches CQFM. Where the prior fails the transport earns its keep: given one transferred from other subjects it regains $+7.2$ TSTR points (9/9).

In the LEPS2 facility constructed at SPring-8, hadron photoproduction experiments have been carried out by using a photon beam with high linear polarization in the tagged energy range of 1.3-2.4 GeV. The BGOegg and Solenoid experiments are alternately running with large-acceptance detector setups, which have been prepared by focusing on the measurement of gamma-rays and charged particles, respectively. These two projects have covered a wide range of hadron physics research, including the search for in-medium meson mass modification, the baryon resonance spectroscopy, the studies on exotic hadrons, and so on. This article reviews the physics programs, the obtained results, and the future prospects in the LEPS2 facility.

We analyze how an arbitrary observation mask on a finite periodic ring of equally spaced sensors affects the posterior covariance, Fourier-mode coupling, and normalized posterior trace in Gaussian-process reconstruction. Under a rotationally stationary prior and homogeneous independent measurement noise, complete observation gives independent scalar posterior formulas for the Fourier modes. For an arbitrary mask, however, the matrix $Q_M=F D_MF^\ast$ is generally non-diagonal; its off-diagonal entries are finite Fourier components of the realized mask and couple modes in the posterior precision. Consequently, masks with the same unavailable-channel fraction can have different normalized posterior traces because their geometries differ. A dimensionless 64-channel synthetic benchmark illustrates this finite-matrix effect. A circumferential array of equally spaced wall-mounted microphones at a fixed axial station of a circular fan or compressor duct provides one concrete mechanical-engineering interpretation: failed, saturated, corrupted, or dropped-out channels form the observation mask. The analysis is a finite-dimensional reference calculation under the stated rotational-stationarity and common-noise assumptions, not a performance claim for a nonuniform or unequally instrumented operating duct.

We present a systematic framework to extract QCD thermodynamic properties at finite baryon density from ultra-central heavy-ion collisions, where the impact parameter is nearly zero and volume fluctuations are strongly suppressed. By mapping the measured multi-particle state of each event to an effective homogeneous fireball through the conservation of total energy, total entropy and net baryon number, a near-thermal state with variations is realized from event to event, which allows thermodynamic response relations to be tested. Using STAR Beam Energy Scan data on charged particle multiplicity and net proton fluctuations, we show that the thermalization condition is satisfied within uncertainties across collision energies, validating the interpretation of ultra-central events as small deviations from a thermally equilibrated state. We then generalize the response relations to finite baryon chemical potential, expressing the response coefficients in terms of the speed of sound and susceptibilities. Our results provide a baseline for using ultra-central collisions to probe the QCD equation of state and phase structure, including the search for the critical endpoint.

Nonreciprocal many-body dynamics can develop a swap phase characterized by spontaneous global oscillations far from equilibrium. To determine how local nonreciprocal interactions predict global oscillations and how broadly this prediction applies, we compare mean-field theory, independent-node theory, and nearest-neighbor spin-coupling theory for a two-spin nonreciprocal Ising model on networks of fixed coordination number (\(z=4\)). Nearest-neighbor spin-coupling theory retains nine independent joint edge probabilities. Monte Carlo benchmarks show that it predicts the oscillation amplitude, angular momentum, and period on edge-swapped networks with a 2.1% aggregate mean absolute relative error. We then extend the square lattice to degree-preserving edge-swapped networks and find sharply different behavior. With increasing system size, the oscillation amplitude and angular momentum decay strongly on the periodic square lattice but remain nearly size independent on edge-swapped networks. These results show that identical local parameters do not determine the macroscopic phase. The contrasting finite-size trends identify network geometry as a key determinant of collective oscillations at fixed local interaction parameters.

We propose a new product entropy, defined as the Rényi (or von Neumann) entropy of a normalized product operator constructed from two density matrices. We establish a duality showing that the SVD entanglement entropy of a subsystem for two pure states is exactly equivalent to the product entropy of the complementary subsystem. This connection provides both a transparent physical interpretation in terms of the spectral diversity of the subsystem state product and a computationally efficient route that bypasses the reduced transition matrix. For low-lying eigenstates, we derive analytical expressions for the subsystem product entropy between the ground state and primary excitations in two-dimensional conformal field theories, explicitly verified against the critical Ising chain. In quantum quench dynamics, the quasiparticle picture yields time evolution in the scaling limit: following a global quench, the subsystem product entropy exhibits distinct sequences of thermalization and revivals, whereas under a local operator quench, it develops characteristic plateaus whose constant values are determined by the inserted operator. Extensive numerical calculations on the critical Ising chain confirm the analytical predictions with excellent accuracy.

We study measurement-based feedback control of a damped harmonic oscillator, motivated by optomechanical and interferometric position measurements, in which feedback is synthesized from a displacement measurement. For a band-limited derivative controller, the loop shapes the mechanical susceptibility and reinjects measurement noise, creating a tradeoff between stabilization and noise amplification. We present a compact classical formulation that isolates this tradeoff in the closed-loop displacement spectrum and identifies the parameters governing resonance suppression versus noise-driven actuation. We then extend the same topology to a quantum-limited position detector, enforcing the imprecision-backaction constraint and introducing finite measurement bandwidth that regularizes high-frequency noise reinjection and adds phase lag in the feedback path. In the high-$Q$ regime we obtain an analytic design rule for the optimal feedback gain at fixed controller and measurement bandwidths and compare it with full-spectrum calculations. The comparison shows that a high mechanical quality factor alone does not guarantee quantitative agreement: the accuracy of the near-resonant approximation depends strongly on the controller and measurement bandwidths. Because finite measurement bandwidth changes both the magnitude and phase of the feedback loop, the analysis separates the near-resonant damping contribution from imprecision-noise reinjection and shows that positive near-resonant feedback damping requires the product of the controller cutoff and measurement bandwidth to exceed the square of the mechanical resonance frequency. Finally, we present a dimensionless map of the gain-optimized high-$Q$ occupation and compare its predictions with full-spectrum optimization. All analytical gains evaluated in this map satisfy the full closed-loop stability criterion for these parameters.

We study the retrieval dynamics of an $n$-body asymmetric dense associative memory model, in which sequential retrieval is implemented by shifting one pattern index in the higher-order Hebbian interaction. The model exhibits sequential retrieval of the stored patterns, and the retrieval state corresponds to a stable limit-cycle attractor. The storage capacity is defined as the maximum loading rate for which this limit cycle remains stable. Using generating functional analysis, we derive the dynamical equations describing the sequential retrieval state. The generating functional analysis reveals that the retarded self-interaction term vanishes in the asymmetric model, and the resulting overlap dynamics is reproduced by an Amari-Maginu-type signal-to-noise analysis. Also, we find that for the interaction order $n\ge 3$, the normalized storage capacity decreases monotonically with increasing interaction order, in contrast to symmetric Modern Hopfield models.

In this review we give an overview about dileptons as probes for the properties of the strongly interacting hot and dense matter created in heavy-ion collisions at various beam energies. As penetrating probes they leave the fireball unaffected from final-state interactions and thus provide a space-time-evolution weighted average of the in-medium properties of hot and dense QCD matter.

Mixed-halide perovskite nanocrystal electrodes fail in supercapacitors through field-driven halide segregation. In undoped $\mathrm{CsPbBr_2I}$, this appears as a capacitance that climbs to 218% of its first-cycle value by cycle 1378 and then collapses to 39% by cycle 2500, with complete loss of the iodine signal from the cycled electrode. This work tests whether sodium incorporation suppresses that failure mode. Na-doped $\mathrm{CsPbBr_{3-x}I_x}$ ($x = 0, 1, 2$) nanocrystals were prepared by ligand-assisted reprecipitation at a Na/Pb precursor ratio of 1.25:1.00 and compared with undoped analogues in 0.1 M tetrabutylammonium tetrafluoroborate in anhydrous dichloromethane. Sodium contracts the pseudocubic lattice parameter of $\mathrm{CsPbBr_3}$ from $5.908 \pm 0.029$ to $5.851 \pm 0.019$ $\mathring{\mathrm{A}}$ after correction for specimen displacement. Specific capacitance at 0.3 A g$^{-1}$ rises for every composition, from 42 to 75, from 63 to 96, and from 56 to 84.5 F g$^{-1}$. Na-$\mathrm{CsPbBr_2I}$ gives the lowest charge-transfer resistance at 175 $\Omega$ and the highest ion diffusion coefficient at $1.6 \times 10^{-16}$ m$^2$ s$^{-1}$, and power-law exponents between 0.33 and 0.45 at all potentials examined show that charge storage is limited by ion transport through the pore network rather than by the interfacial process. In the sodium-containing (Na-$\mathrm{CsPbBr_2I}$) electrode, the capacitance rise reaches only 115% at cycle 600, no collapse follows, and 97% is retained at 2500 cycles, the excess above the first-cycle value being reduced by a factor of 7.9. Iodine is retained at unchanged binding energy, and the Br:I ratio measured by elemental mapping is 2.30 after cycling against 2.25 before.

We construct a product spectral decomposition from power-normalised parabolic-cylinder functions. At noninteger real orders, these functions are not square integrable with respect to the corresponding full-line Gaussian measure. Bringing the admissible half-line branches into the same Gaussian space and matching them at the origin gives a corrected orthonormal eigenbasis. A unitary transformation to a weighted probability space makes the eigenfunction corresponding to the lowest eigenvalue constant and permits consistent countable products. The resulting law is non-Gaussian under deformation, and the decomposition recovers classical Wiener chaos when the deformation is removed. We identify the associated closed gradient, adjoint divergence and deformed number operator. The Rodrigues formula covers all orders occurring in the spectral branches, including negative orders. Under deformation, however, changing both branch orders by the same nonzero amount breaks the required relation between them. Fractional evolution defined through spectral calculus preserves the basis. After unitary identification, the one-coordinate fractional evolution operators converge in operator norm to Gaussian and three-dimensional radial Ornstein–Uhlenbeck fractional evolutions at the respective parameter endpoints. Convergence is uniform on compact time intervals away from zero. An application example with inverse stable clock is provided in the accompanying github repository.

The non-standard self-interactions of neutrinos could be induced by (light) scalars that couple primarily to neutrinos. In this work, we investigate the effects of neutrinophilic scalars $\phi$ on muon decays at the tree, 1-loop and 2-loop levels, considering both Dirac and Majorana neutrinos. The most phenomenologically interesting processes are the 1-loop decay $\mu \to e \phi \phi$ and $\phi$-induced decay $\mu \to e \gamma$ at the 2-loop order from the same coupling $g_{\rm eff}$, and the 1-loop corrections of $\phi$ to the SM muon decay and the four-body muon decays $\mu^- \to e^- \nu_\mu \nu_\mu \phi,\, e^- \bar\nu_e \bar\nu_e \phi$, both induced by the coupling $h_{e\mu}$ of $\phi$ to Majorana neutrinos. The current experimental bounds on $\mu \to e \gamma$ provide the most stringent limits on the coupling $g_{\rm eff}$, up to roughly 0.064, for the scalar mass $m_\phi \gtrsim 8.7$ MeV. The precise $G_F$ measurements constrain the coupling $|h_{e\mu}|$ down to roughly $0.027$, with the scalar mass $m_\phi$ constrained up to roughly $5.9$ TeV. For $m_\phi \gtrsim {\cal O}({\rm GeV})$, the $G_F$ measurements are more stringent than other existing laboratory, astrophysical and cosmological constraints. Some other $\phi$-induced processes are also possible, which are, however, heavily suppressed by the tiny neutrino masses and/or the loop factor.

In this paper, we introduce a method to constrain the variation of fundamental physical constants predicted by grand unified theories. This approach utilizes OH sum rules, which exhibit different forms depending on the physical conditions of the OH gas. Employing a least-squares technique, we can directly derive the relative velocity offsets among four OH lines in each ground state. This circumvents the need for spectral profile fitting, which is often unfeasible for spectra with complex profiles. Moreover, this method remains robust even when portions of the spectra are contaminated by radio frequency interference (RFI). Applying this method to a Galactic source G035.1+0.5, we obtained a stringent constraint of $\Delta X/X=(-1.8\pm5.1)\times10^{-6} $, where $X=\mu^{-1}\alpha^{2}g_{p}^{0.54}$. This method is expected to play a crucial role in extracting fundamental physics from large datasets of upcoming Square Kilometre Array (SKA) telescope.

We present the mid-infrared variability study of 120 long-lived disk candidates (ages 10-100 Myr) and 13 confirmed "Peter Pan" disks using NEOWISE to constrain the longevity of dynamic disk processes beyond the standard 10 Myr dispersal threshold. We identified 32 sources exhibiting significant variability by combining stochastic flux metrics and the Stetson $J$ index. The majority of these variables belong to the 10-25 Myr age range, suggesting a sharp decline in disk brightness variability after 25 Myr. Morphologically, the variability is characterized by irregular fluctuations with mean amplitudes of $\Delta W2 \approx 0.5$ mag, consistent with variable extinction or inner disk warping induced by embedded planetesimals. We also report a strong correspondence between mid-infrared and optical variability for the Peter Pan disks, with all 9 mid-IR variables confirmed as optical variables in the literature. Our findings reveal a population of 23 new long-lived disk candidates that retain active, gas-rich environments, supporting the existence of a rarefied class of systems that survive well into the debris disk era.

Turbulent fluid simulation is computationally demanding because nonlinear interactions couple a wide range of scales. Quantum computing offers another computational model, but fluid evolution is nonlinear and dissipative whereas closed gate-based dynamics is linear and reversible. The lattice Boltzmann method (LBM) is an attractive discrete target because streaming is a local permutation, yet conventional collision reconstructs a nonlinear equilibrium and applies dissipative relaxation. We present a digital quantum lattice Boltzmann (QLBM) algorithm that keeps that nonlinear timestep in computational-basis registers. Each step writes the updated populations into a clean destination bank, uncomputes the workspace, and resets the obsolete source. The resulting channel is completely positive and trace preserving (CPTP), composes without intermediate measurement, and applies to any lattice Boltzmann stencil that can be evaluated reversibly into a clean bank. Forced D3Q19 homogeneous isotropic turbulence, in direct numerical simulation (DNS) and Smagorinsky large-eddy simulation (LES) at Reynolds number 15000 over one million lattice steps, shows statistical agreement in mass, energy budgets, spectra, and intermittency. The three-dimensional calculations use a digital-register emulator of the circuit. A separate D2Q9 implementation compiles the same contract to reversible gates and agrees bit-exactly with its integer twin. Two-bank storage and serial lattice depth limit near-term feasibility.

I consider vortex solitons in two-dimensional star networks created by several one-dimensional line waveguide arrays (or rays) brought in close proximity. First waveguide in each array belongs to the ring of certain radius and its modal field overlaps with modal fields of waveguides from other rays laying on this ring, as well as with modal fields of waveguides in this line array. This creates a system with discrete rotational symmetry that is capable of supporting various localized linear vortex states due to local coupling between waveguides on the central ring, where line arrays approach each other. Localized vortex states emerge from the continuum of the states extended along the rays of the structure upon variation of the radius of the central ring. I show how different vortex modes with progressively increasing topological charges emerge when the number of rays in structure and, therefore, its discrete rotational symmetry increases. Each such linear localized mode gives rise to a family of thresholdless vortex solitons in nonlinear optical medium, whose topological charges are limited by the discrete rotational symmetry of the system. Vortex solitons in star networks exhibit qualitatively different behavior with increasing power depending on whether the propagation constant of linear vortex that gives rise to a family of solitons is located above or below the band of extended states. Stability properties of vortex solitons nontrivially depend on their topological charge and discrete rotational symmetry of the system, showing very different picture from usual ring-like waveguide arrays.

Electrides are an unconventional class of materials in which electrons are localized in crystallographic void spaces rather than solely around atomic nuclei, giving rise to appealing properties such as low work functions, strong electron-donating character, and even superconductivity. Here, we use ab initio methods to investigate metal-organic framework (MOF) electrides, a new class of materials that combines the interstitial electrons of electrides with the permanent porosity and chemical tunability of MOFs. These materials host pore-confined electrons: occupied electronic states localized in the pore space and with bands slightly below or crossing through the Fermi level. Using density functional theory calculations, we establish several design rules for stabilizing pore-confined electrons in MOFs via an anion-electron exchange process and identify candidate MOF electrides. As a proof-of-concept, we also demonstrate that the pore-confined electrons can directly facilitate chemical reactions, substantially lowering the activation barrier for H2 dissociation without requiring adsorption at a surface site. We envision that pore-confined electrons in nanoporous materials may enable a fundamentally new type of catalysis in which chemical reactions take place in the pore space, driven by electron-centered active sites.

Constrained binary optimization is a representative class of NP-hard problems in scheduling, resource allocation, and finance. Segmented feasibility-preserving variational quantum algorithms (VQAs) are a promising approach that restricts ideal circuit evolution to feasible assignments and executes short measure-and-reseed segments. The existing boundary runtime enforces feasibility through purification, which discards measurements that violate the constraints. As problem size and noise increase, feasible measurements become rare; if no shot survives, the execution chain terminates. In our 72-qubit graph-coloring experiment on a real-device IBM Heron, 11 of 12 purification chains terminate before completing all segments. We propose SUTURE: syndrome-guided repair, a runtime that repairs infeasible measurements instead of discarding them and is deployable on current quantum devices. SUTURE exploits a parity-check-like structure induced by the problem constraints: violated constraints form a syndrome that detects and localizes corruption and, under suitable structural conditions, often identifies a correction. SUTURE combines three stages: (1) a compile-time profiler that uses the compiled constraints and feasible initialization samples to predict single-flip recovery, with a maximum error of 0.011 across four constraint families; (2) a bounded runtime decoder that replaces purification inside the segmented execution loop; and (3) a regime analysis that identifies when repair is preferable to purification. In simulation up to 120 qubits, SUTURE continues the execution chain far beyond the noise level at which purification collapses. In the same 72-qubit experiment on IBM Heron hardware, SUTURE completes all segments in all 12 runs, and hardware timing experiment, decoding adds 1.6% to execution-path latency and less than 1% to total pipeline latency.

We consider an unknown pure stabilizer state of \(n\) qudits with a fixed prime local dimension \(p\), together with prescribed commuting Pauli observables. Rather than learning the whole state, the learner must recover only complementary stabilizers that commute with the prescribed observables, together with their eigenvalues. Once the prescribed measurement is performed, this information and the observed outcome determine the corresponding conditional stabilizer state. We compare the number \(k\) of copies available to the learner with the number \(n\) of qudits and refer to their ratio \(k/n\) as the copy rate. Suppose that the dimension \(m\) of the requested complementary stabilizer subspace satisfies \(m=\beta n+o(n)\) for a fixed \(0<\beta\leq1\). If the copy rate converges to a value below one, the optimal exact recovery probabilities and verification scores converge to zero, both in the worst case and under the uniform prior. Conversely, if the number \(k\) of copies minus the number \(n\) of qudits tends to positive infinity, the entire stabilizer state can be identified with probability tending to one, and the requested complementary information can then be extracted. Thus, restricting the learning target to the stabilizer information needed after the prescribed measurement yields no copy-rate discount.

We prove the all-genus conifold gap conjecture for the Gromov–Witten (GW) theories of the quintic threefold, other smooth Fermat Calabi–Yau threefold hypersurfaces, and local $\mathbb P^2$, as well as for the FJRW theories associated with these hypersurfaces. Our proof uses Mixed Spin P-field (MSP) theory and analogous master-space constructions. The central object is the cone-vertex theory. We express this theory as a translation of the point CohFT via a formal Laplace transform and prove a universal conifold gap theorem for the resulting translated theories. The master-space construction then transfers the gap from the cone-vertex theory to the GW and FJRW theories.