Papers

Context. Ultraluminous X-ray sources (ULXs) are thought to be powered, in many cases, by super-Eddington accretion onto compact objects. While soft X-ray lags have been detected in several ULXs, hard lags remain rare and poorly understood. Aims. We investigate the temporal and energy dependence of X-ray lags in NGC 7456 ULX-1 to constrain their physical origin and probe the super-Eddington accretion flow. Methods. We analyzed the two deepest XMM-Newton observations, taken in 2018 and 2023. Hard (1-10 keV) and soft (0.3-1 keV) light curves were cross-correlated using adaptive-binning techniques optimized for Poissonian low-count data. We measured lags over consecutive 10 ks intervals and investigated their energy dependence. Spectra were modeled using thermal and Comptonization models. Results. We detect significant hard X-ray lags in both observations, with the hard emission delayed by $\sim 10^3$ s during phases of rapid flux variability. The delays are primarily driven by the lowest-energy photons. Spectral modeling indicates a Comptonization-dominated flow comprising a cooler, extended outer region and a hotter, compact inner flow embedded in an optically thick wind. We interpret the delays as the combined effect of inward propagation of accretion-rate fluctuations and photon diffusion within the dense outflow. Fluctuations first enhance the soft-emitting outer regions and then propagate toward the hotter inner flow, where photons undergo stronger Comptonization before escaping with a kilosecond delay. The small inferred inner emitting radius disfavors an intermediate-mass black hole accretor. Conclusions. The sign, amplitude, and energy dependence of the delays disfavor standard reverberation. Propagation-driven variability coupled with radiative transfer in optically thick winds appears to play a major role in shaping the timing properties of super-Eddington accretion flows.

Simulating microstructure evolution requires quantum-mechanical accuracy and mesoscopic reach in length and time scales, a combination that no current method achieves. Classical phase-field models provide this reach, but their accuracy is limited by phenomenological free energies and mobilities. Here we develop a framework for learning ab initio phase-field models, where the mesoscopic equation is not postulated but derived from a Mori-Zwanzig projection of molecular dynamics onto species-density fields under explicit assumptions. The nonlocal free energy and mobility left unspecified by this equation are parametrized by neural networks and learned from short molecular dynamics trajectories generated with machine-learning interatomic potentials of ab initio accuracy. We demonstrate the framework on an iron-boron melt and on hydrogen-helium mixtures under planetary conditions. For iron-boron, the model shows that the melt at the FeB$_4$ composition is spinodally unstable at ambient pressure but stabilized at 10 GPa, offering a thermodynamic rationale for why FeB$_4$ has been synthesized only under high pressure. For hydrogen-helium, the model predicts the immiscibility boundary and captures droplet nucleation and growth in helium-rain simulations of a column corresponding to 2.2 million atoms, far beyond the scale of atomistic modeling at comparable accuracy. Trained across compositions and conditions, such models could provide a mesoscopic counterpart to ab initio molecular dynamics.

Multiple-choice questions are a staple of educational assessment due to their efficiency, but their diagnostic ability is inherently limited. Traditional psychometric techniques infer misconceptions by fitting models to predetermined data patterns, making them ineffective for analysing small, novel, or non-standardised samples where such prior trends are absent. To address this, we introduce a new distractor analysis method based on Information Theory, which treats the randomness in incorrect answers not as noise but as a measurable signal in student answer patterns. By examining the entropy of distractor choices alongside overall accuracy, our approach offers a sample-sized-independent framework for classifying performance. Cross-referencing low-scoring item-level response data against independently documented misconceptions on a standard mechanics inventory shows that, for items identified as most concentrated (lowest degree of randomness), a majority of incorrect responses select the exact response coded to that item's named, interview-validated misconception. Furthermore, items independently flagged elsewhere as diagnostically unreliable are exactly the items the measure identifies as having the highest degree of randomness. The pedagogically actionable cases are those with low degree of randomness in a class's incorrect answers, signalling a shared, re-teachable misconception, whereas a high degree of randomness chiefly serves as the null case against which that concentration is judged. Our analysis is intended as a fast, single-administration triage step: a practical, theory-driven complement to established distractor-analysis methods that helps instructors flag which items warrant closer, more resource-intensive investigation.

Structure formation has been used to test the viability of the IR cutoff in the THDE model, and the future event horizon was found to be a viable IR cutoff (Phys. Rev. D 113, 123531 (2026)). Motivated by this, in this paper, using the combined Pantheon+ SN Ia, OHD, DESI DR2, and CMB distance prior datasets, we constrain the model parameters for the HDE, THDE, and BHDE models, and use parameter sets from the MCMC samples to analyze the evolution of $f\sigma_{8}(z)$ and $S_{8}(z)$ for these models. We find that the evolutionary curves of both $f\sigma_{8}(z)$ and $S_{8}(z)$ for HDE and BHDE are nearly indistinguishable, while the curves for the THDE model lie lowest. The $\chi^{2}$ statistic shows that the THDE model yields the lowest $\chi^{2}$ value for both $f\sigma_{8}$ and $S_{8}$ measurements. Finally, we construct the ($f\sigma_{8}$,$S_{8}$) plane as a diagnostic to analyze the trajectories of these models, and find that the trajectory of the THDE model can achieve the current ($f\sigma_{8}$,$S_{8}$) data region. These results demonstrate that structure formation can distinguish THDE from HDE and BHDE, and suggest that structure formation can serve as a useful probe for distinguishing HDE models and can be extended to other HDE scenarios.

The paper proposes a generalization of the fully implicit energy- and charge-conserving electromagnetic particle-in-cell method to the case where the lightest type of plasma particle (electrons) is described in the drift-kinetic approximation. This allows us to remove the very strict time step limitation of this method requiring to resolve the gyrorotation of electrons. Since the drift-kinetic model is only applicable to light particles whose contribution to plasma polarization is small, this model can be further simplified by neglecting the electron polarization drift and including in Maxwell's equations, in addition to the current of gyrocenters, only the magnetization current. In order for such a hybrid model to retain conservative properties in finite-difference form, a method of self-consistent interpolation of $\nabla B$ in the mirror force from grid to particle and the magnetization vector from particle to grid is proposed. Unlike existing drift-kinetic models, this model is not limited to considering small perturbations near a given equilibrium, and therefore allows one to study the formation of plasma equilibria in regimes with a finite ratio of plasma to magnetic field pressure. Testing of the parallel code implemented in C++ using the PETSc library confirmed the fulfillment of the finite-difference laws of energy and charge conservation, as well as the ability of the model (in the drift-kinetic version for all types of particles) to correctly reproduce the diamagnetic effect and longitudinal ion-acoustic wave.

Structured light can generate electronic dc currents with azimuthal winding, yet the rules governing their winding order $m$ remain unclear. Here we identify sector-resolved winding selection rules using graphene as a clean two-dimensional platform. By decomposing the current response into local and gradient sectors at the current-operator level, we show that the winding order is determined not by the optical orbital angular momentum $\ell$ alone, but by the angular structure of the corresponding current operators together with projection onto the azimuthal direction. For scalar Laguerre-Gaussian beams, linear polarization yields $m=|\ell\pm1|$ in the local sector and $m=|\ell|,|\ell\pm2|$ in the gradient sector, whereas circular polarization with helicity $\sigma=\pm1$ selects $m=|\ell-\sigma|$ in the local sector and $m=|\ell+2\sigma|$ in the gradient sector. Numerical time-evolution calculations verify these rules and further show that helicity can select an $m=0$ branch, producing an azimuthally uniform circulating current whose radial profile determines the axial magnetic field $B_z(z)$. Our results provide a sector-resolved organizing principle for classifying and controlling structured-light-driven dc currents with tailored winding structures.

Assuming post-quantum one-way functions, we construct device-independent (DI) oblivious transfer (OT) and bit commitment: honest parties use only trusted classical computation to operate untrusted quantum devices, which may share arbitrary entanglement and behave non-IID. Security is simulation-based against quantum polynomial-time adversaries and composes sequentially with efficient simulators. One protocol skeleton serves both, in two regimes. With isolated laboratories and coordinate-local measurements in the honest receiver's device, it tolerates a constant rate of honest-device faults. With polylogarithmically many qubits of adaptive leakage between the laboratories and arbitrary joint measurements, it tolerates an inverse-polylogarithmic rate. Each elementary DI call uses a fresh, isolated batch of polylogarithmically many device coordinates, and total device use in the compiled OT protocol is polynomial. The commitment has efficient simulators against both parties and yields DI coin tossing with abort. Because OT is complete for secure computation, the construction yields a DI protocol, with abort, for every efficiently computable classical functionality on a fixed number of parties, secure against static corruption of any proper subset of them. The commitment's extractor changes a public parity relation through classical equivocation and leaves the device execution, hence its leakage, unchanged. A commit-and-prove functionality, disjoint audits, and an affine consistency check link the certified correlations to ideal OT. Sender security rests on a selector-aware parallel-repetition bound for the Magic Square game, which we derive from the two-round threshold theorem of Kundu and Tan.

We discuss mechanisms of $D \bar{D}$ pair production in $e^+ e^-$ collisions. We consider the $D^0 {\bar D}^0$-continuum mechanism with the $t/u$-channel $D^*$ vector-meson exchanges, where the $D^* D \gamma$ coupling constants are derived from the $D^* \to D \gamma$ decays. For the $D^+ D^-$ channel, we consider mainly the $D^{\pm}$ exchanges. We find a relatively large continuum contribution for the $D^0 {\bar D}^0$ channel and a significantly smaller one for $D^+ D^-$, concluding that the broad bump at $M_{D^0 {\bar D}^0} = 3.8$ GeV observed by the Belle and BaBar Collaborations has a continuum origin rather than corresponding to the broad $\chi_{c0}(3860)$ resonance. We also investigate the resonant production of $D \bar{D}$ pairs, including $\chi_{c0}(3860)$ and $\chi_{c2}(3930)$. The $\chi_{c2}(3930)$ meson is a good candidate for the first excited $P$-wave charmonium $\chi_{c2}(2P)$ state. This resonance can decay into both $D \bar{D}$ channels; however, its branching fractions are not well known at present. Comparing our model result with $\Gamma(\chi_{c2}(2P) \to \gamma \gamma) = 0.544$ keV (corresponding to the Buchmüller-Tye potential) to the BaBar data, we find $B(\chi_{c2}(3930) \to D \bar{D}) = 0.58 \pm 0.13$. We present results for differential distributions in several variables for Belle II kinematics. Finally, we shortly mention predictions for the ultraperipheral heavy-ion reaction $Pb Pb \to Pb Pb D \bar{D}$ for low invariant mass ($M_{D \bar{D}} < 4.3$ GeV) at the LHC.

The photostability of DNA nucleobases relies on ultrafast nonradiative relaxation following ultraviolet excitation. In thymine, the competition between internal conversion (IC) and intersystem crossing (ISC) determines the balance between ultrafast energy dissipation and potentially long-lived photochemistry. However, the precise molecular mechanisms linking early singlet-state electronic and structural evolution to the formation of long-lived triplet dark states remain unresolved. Here, we investigate the gas-phase photodynamics of thymine using time-resolved multi-center X-ray photoelectron spectroscopy (tr-mc-XPS) supported by state-of-the-art ab initio quantum calculations. By simultaneously probing multiple carbon 1s core levels, we track site-specific electronic and nuclear dynamics. We discover that the long-lived dark state is not initially a pure triplet; rather, it exists as a mixture of singlet and triplet states for the first 240 ps before transitioning fully to the $^3\pi\pi^*$ state. On the femtosecond timescale, we resolve atom-specific vibrational coherence at $\sim$730 cm$^{-1}$ (a ring-breathing mode with a period of 46 fs) that survives the initial $^1\pi\pi^*$ to $^1n\pi^*$ IC. Furthermore, we identify a critical dual role for the methyl group (C9H$_3$). During singlet relaxation, the methyl group acts as an inertia spectator, and its mass hinders direct internal conversion to the ground state. Upon $^3\pi\pi^*$ state formation, however, state-dependent hyperconjugation couples the methyl group to the pyrimidine ring, turning it into a highly sensitive electronic reporter of ISC. These findings establish tr-mc-XPS as a powerful approach for disentangling the complex, multi-timescale photodynamics governing nucleobase stability.

We find the optimal initial states of a linearized cavity-optomechanical system for estimating the single-photon coupling by quantum Fisher information. In the red-detuned (beam-splitter) regime the interaction exchanges excitations between the two modes, and for pure inputs the Fisher information is four times the variance of the excitation-exchange operator, reducing the optimization to variance maximization with a closed-form solution. With the mechanical mode in a Fock state the probe is degenerate: any state saturating the photon-number constraint is optimal, and a Fock reference with many phonons enhances the per-photon sensitivity linearly. For an arbitrary fixed mechanical reference it becomes a phase-matched squeezed vacuum, from an exact bound on that coherence obtained by a Lagrange-dual argument. Jointly optimizing both modes at fixed total excitation number, the quantum Fisher information obeys a Heisenberg bound quadratic in the energy, saturated by equal superpositions of the two extremal eigenstates of the generator, equivalently a two-mode NOON state after a balanced beam splitter; balanced product Fock states reach only a linear enhancement. Exact diagonalization in truncated Fock spaces confirms all the results, the familiar single-photon input being the simplest saturating case. Three extensions are worked out: the blue sideband has an unbounded two-mode-squeezing spectrum, so the entangled optimum becomes impossible and the balanced product Fock state is exactly optimal at definite energy; under cavity loss the information accumulation turns linear on the cavity-lifetime timescale, at a renewal-optimal duration of about two and a half lifetimes; and for a state-of-the-art quantum-coherent coupling experiment the per-hour relative precision lies between one part in a million and ten parts in a million, number and thermal references being the most robust.

Ringdown spectroscopy probes black hole (BH) dynamics, and exceptional points (EPs) can reveal underlying degrees of freedom. We develop a procedure to identify EPs in ringdowns using a time-weighted Jordan term and a splitting bound. For relative Jordan excitation proportional to separation, we derive a sixth-order squared distance between split and EP waveforms after fitting the resonance center. Independent constraints on these quantities can restore fourth order. Complete waveforms show that omitting ordinary responses can strengthen EP evidence and spuriously tighten splitting bounds.

We analyze thermalization and eigenstate thermalization in the random free fermion model, a quadratic Hamiltonian with a Gaussian Orthogonal Ensemble hopping matrix, by an exact random-matrix calculation. For a single-particle observable $A$ we compute the ensemble average and fluctuations of $\operatorname{Tr}(A\rho(t))$: we derive the relaxation envelope $g^2(t/\tau_\lambda)$ with $\tau_\lambda = \hbar/(2\eta\sqrt{N})$ and its oscillatory $|t|^{-3}$ tail, both properties of the semicircular limiting density, and we show that the variance of the fluctuations is $O(1/N)$, so that the relaxation profile self-averages at fixed scaled time. Incorporating particle-number conservation, we prove a factorization theorem for the energy-shell statistics of the many-body eigenstates: the energy-independent ETH structure follows from the statistical independence of the eigenvector matrix and the spectrum, and at fixed realization the energy constraint enters only through a covariance subtraction of relative order $1/N$. A controlled eigenvector–eigenvalue correlation deformation $H(\varepsilon)=H_0+\varepsilon H_1$ then generates, at first order in $\varepsilon$, the energy dependence of the diagonal energy-resolved slope and of the off-diagonal spectral function, a single-particle diagnostic of the ETH structure in the allowed particle-hole sector. Correlation-function fluctuations and the single-particle spectral form factor finally identify the subleading $O(1/N)$ late-time fingerprints of two-level spectral correlations. The model thus separates eigenbasis chaos from spectral chaos and identifies the statistical origin of energy-independent versus energy-dependent ETH.

We show that there exists a theoretically well-motivated Higgs sector in which the conventional heavy-Higgs production mechanisms are suppressed. In this limit, the $b$-associated production channel assumes the role of the primary production mechanism. We demonstrate that the democratic three Higgs-doublet model (3HDM) provides a viable beyond the SM (BSM) framework in which this limit can be realized. We further show that existing $b$-associated searches provide the most sensitive LHC probe of this 3HDM scenario in the up-phobic limit and derive the corresponding constraints on its parameter space. We also present projections for future collider searches, demonstrating their potential to probe further regions of the model's up-phobic parameter space.

The Starobinsky model provides a successful realization of inflation and serves as a benchmark for alternative inflationary models. Motivated by its observational success, we investigate the dynamical evolution of the early Universe within this framework. Incorporating radiation production from scalaron decay, we formulate a coupled system of dynamical equations describing scalaron depletion and radiation buildup while ensuring exact energy conservation in the expanding background. In the Friedmann equation, both scalaron and radiation components gravitate, thereby accounting for the backreaction of the generated radiation. Starting with a radiation-free Universe and initial conditions from the standard slow-roll approximations, exact numerical evolution shows that inflationary dynamics is predominantly governed by the scalaron. Despite continuous radiation production, the radiation component remains subdominant, and the slow-roll conditions persist throughout inflation, yielding approximately 60 e-folds. Following inflation, the scalaron enters a damped oscillatory phase accompanied by rapid radiation growth. This transition provides a graceful exit from inflation and a smooth entry into reheating. During reheating, radiation grows continuously as the scalaron energy is depleted, with the radiation energy density ultimately becoming substantially larger than the residual scalaron energy density. This indicates the onset of a hot, radiation-dominated phase, satisfying the Kofman-Yi criterion for successful reheating. Our results demonstrate that the coupled scalaron-radiation dynamics provides a self-consistent realization of graceful exit and the transition to radiation domination in the Starobinsky model.

A novel analytical framework, designated as the JO&Minna method, is presented for interpreting the depth profiles of residual activation induced by heavy-ion beams in thick structural materials. By combining the limiting fragmentation and factorization principles, the JO-factor formalism rigorously links the spatial distribution of projectile-like fragments to the range derivatives and electronic stopping powers. Furthermore, the Minna procedure establishes a mathematical link via profile differentiation and convolution analysis to reconstruct the original energy spectrum and energy spread of the fragments from the measured activity profiles. The methodology is benchmarked using experimental data for $^{46}\text{Sc}$ production from $^{238}\text{U}$ irradiation in aluminum, demonstrating a reliable and straightforward approach for extracting fragment production kinematics and energy distributions in accelerator materials research.

Co-propagating chiral Majorana edge modes are predicted to exist in proximity-coupled quantum anomalous Hall systems. They combine the exciting prospects of edge vortices as flying non-Abelian Ising anyons, with their fully electrical detection through anyon fusion processes. In this setup, electrical transport arises from the interference of a pair of Majorana edge states. By bosonizing the model with an additional replica sector, we generalize previous treatments to fully asymmetric configurations. We analyze a Mach-Zehnder interferometer and show that while the DC conductance does not exhibit signatures of edge-vortex interference, at low frequencies, this interference manifests as an effective capacitance. This quantity is sensitive to the non-Abelian anyon statistics of edge vortices and, in particular, carries signatures of their topological spin and nontrivial conformal dimension. We demonstrate that the thermal length $l_T = \hbar v/(\pi k_B T)$ is the relevant scale at low temperatures $T$, analogous to quasiparticle transport in conventional quantum Hall edge states with velocity $v$. As long as the interferometer arm lengths match within this length scale, transport signatures of the non-Abelian statistics of edge vortices remain visible.

The small scale density distribution of the cold circumgalactic medium (CGM) plays a central role in galaxy evolution, yet it is still poorly understood. We provide new constraints on the cold CGM density distribution by comparing the CGM emission from MUSE observations of quasar (QSO) nebulae at $z \sim 3.8$ with mock observations generated from the new DaLya simulation suite under the assumption of maximal fluorescence. These are further complemented by mock observations generated from the COLIBRE (m5 and m6 resolutions) and HELLO simulations. We first investigate the distribution of the observed and simulated Ly$\alpha$ surface brightness (SB) values, which directly depend on the cold gas density distribution under our assumption, finding stark differences between observations and upper limits provided by simulations with the exception of the COLIBRE higher-resolution run (m5). The observations generally show a much larger covering fraction of regions with relatively large Ly$\alpha$ SB values and suggest that densities reaching at least 10 cm$^{-3}$ must be present in the CGM of QSOs within a radial distance of 50 to 200 ckpc. To understand the origin of this discrepancy, we study the cold CGM density distribution finding in all cases simulated PDFs consistent with skewed lognormal ones with medians $\sim 10^{-2}$ cm$^{-3}$ and clumping factors between 16 and 102. We compare the broadness of these distributions to the one implied by the observed HeII/Ly$\alpha$ CGM ratio, which is independent on the assumption of maximum fluorescence. Assuming lognormal PDFs with median values similar to the simulations, the observed ratio requires a broader PDF with a clumping factor at least 8 times larger than predicted by simulations. These results suggest that higher resolution and/or additional physical mechanisms are needed to reproduce the observed CGM of QSOs at these redshifts.

Pushing the brightness limit of electron sources requires simultaneously minimizing the intrinsic emittance and maximizing the accelerating field. Alkali-antimonide photocathodes exhibit excellent properties at low fields, yet their high-field photoemission physics remains poorly understood owing to acute vacuum sensitivity. Here, we demonstrate robust photoemission from alkali-antimonide photocathodes in a radio-frequency gun at peak fields exceeding 100 MV/m, with quantum efficiency maintained above 1% for more than two weeks, enabling the first systematic measurements of their high-field brightness limits. The measured field dependence of the intrinsic emittance provides direct insight into photoemission physics at high fields. Near-threshold photoemission further reduces the intrinsic emittance, increasing the attainable brightness, and reveals constraints imposed by surface roughness. This work opens the high-field regime for advanced photocathodes, extending the brightness limits of electron sources.

We investigate the gravitational $D$-form factors of the nucleon and pion within an extended linear sigma model (eLSM) that incorporates both two-quark and four-quark states, as well as nucleon fields. Within the scalar-meson-dominance picture, we find that meson mixing between two-quark and four-quark states suppresses the nucleon $D$-form factor and leads to agreement with the lattice-QCD results. Interestingly, a quantitatively similar suppression was previously found in the parity doublet model (PDM) as a consequence of a nonzero chiral invariant component of the nucleon mass. This suggests that the effects of meson mixing and the chiral invariant mass are difficult to distinguish from the nucleon $D$-form factor alone. We therefore examine the pion $D$-form factor in the chiral limit. While its value in the forward limit remains fixed at $-1$ in both the eLSM and the PDM, meson mixing modifies its momentum-transfer dependence. In contrast, the chiral invariant mass does not directly affect the pion $D$-form factor. These results suggest that a combined analysis of the nucleon and pion $D$-form factors can provide further insight into the mass decomposition of the nucleon and the quark composition of mesons.

We analyze the correlations of the vacuum state of a charged scalar field by constructing localized observables in two disjoint regions. We consider a time-dependent homogeneous electric field and examine the time evolution of correlations by relating it with the particle-antiparticle creation through the Schwinger effect. Our results show that the electric field increases the mixedness of the states and enhances total correlations, while the distillable entanglement among the regions decreases, vanishing at a finite time for sufficiently intense pulses within the non-perturbative pair-production regime. We also discuss the effect of anisotropy induced by the external background, showing that the distribution of correlations acquires a nontrivial directional dependence.

We have theoretically investigated the stochastic switching dynamics of a charged microparticle in a quadrupole trap near the principal parametric resonance at atmospheric pressure. We show that this behavior arises from a balance of thermal noise-activated escape and dissipative return between linear and nonlinear dynamical regimes of the particle motion. The switching follows an effective Gibbs distribution with an effective potential barrier between the regimes and effective system temperature, which are both experimentally accessible. Because the effective barrier depends on particle mass and charge, and since thermal noise reveals rather than obscures the transition, this exponential dependence could provide the foundation to precision noise-aware detection schemes. The system further offers a highly adjustable platform for studying activation dynamics in non-equilibrium parametric systems.

Peptide sequence governs both the structures formed through supramolecular assembly and the dynamics by which they emerge, but predicting either requires resolving slow collective processes among many interacting molecules. Molecular dynamics (MD) provides microscopic insight into these processes, yet the long timescales of assembly and the vast peptide sequence space make systematic exploration computationally demanding. We introduce SupraTITO, a transferable generative molecular dynamics (GenMD) framework for supramolecular systems, demonstrated through peptide self-assembly. SupraTITO learns transferable implicit transfer operators (TITO) conditioned on peptide sequence, molecular topology, and periodic geometry, allowing configurations to be propagated over physical intervals much longer than an MD integration step. On a comprehensive dipeptide benchmark, SupraTITO generalizes to held-out sequences and reproduces sequence-dependent structures and dynamics while maintaining molecular integrity over long rollouts. Compared with direct ensemble prediction trained on the same trajectory data, SupraTITO more accurately reproduces assembly structures while also resolving their temporal evolution. The learned dynamics generalize across peptide concentrations, including dilute conditions not represented during training. These results extend transferable GenMD to collective dynamics in periodic supramolecular systems and provide a foundation for modeling related processes beyond peptide assembly.

A system has a nonreciprocity when the light propagating through it in opposite directions encounters different optical adventures. This effect lies in the origin of various applications, including signal processing, noise reduction, unidirectional propagation, and sensing. Here, we experimentally demonstrate the detection of nonreciprocity in plasmonic metasurfaces and a method to determine its direction. We use the Inverted Leakage Radiation Microscopy setup to show that the non- reciprocity is linked to a false chirality of the medium. The setup behaves as an equivalent of the classical Sagnac loop providing an interference pattern at the image plane. The non-reciprocity is then detected by examining the topological order of the resulting image. Our system represents a compact and versatile method to detect false chirality in plasmonics and suggests a variety of potential areas of implementation in nanophotonics, chemistry and applied physics.

Vacancy ordered double perovskites (VODPs) have been known for more than 100 years, but their exact crystal structures and magnetic ground state remained controversial till date. Anhydrous VODP iridates ($A_2$Ir$X_6$; $A$ = K$^+$, (NH$_4$)$^+$, Cs$^+$, Rb$^+$, and $X$ = F$^-$, Cl$^-$, Br$^-$), having $face-centered~cubic$ (fcc) antiferromagnetic ground state, offer a chemically well-defined platform to investigate the influence of on-site Coulomb repulsion, crystal field, spin-orbit coupling, and competing Heisenberg and Kitaev interactions. Interestingly, among these, the Na-based VODP iridates seem to possess an unusual hydration dependent structural and magnetic uncertainty. As a result, a cubic Na$_2$IrCl$_6$ system having fcc antiferromagnetism could not be stabilized ever. Here, for the first time, we report the synthesis and comprehensive characterization of anhydrous Na$_2$IrCl$_6$ single crystals having fcc structure with antiferromagnetic transition temperature of 1.91 K. Thermogravimetric analysis establishes the absence of lattice water and also demonstrates thermal stability of the crystals up to 453 K. The cubic phase reported here contrasts with previously described 'anhydrous' monoclinic Na$_2$IrCl$_6$ and highlights the strong structural and physical property relationship in Na-based VODP iridates. This new material offers an exciting playing ground for exploring another potential perfect $J_{eff}$ = $\frac{1}{2}$ octahedra and the consequent strong Kitaev interaction within the nearest-neighbor (nn) pairs of Ir ions.

Average atom (A-atom) potentials provide a mean-field description of a chemically disordered alloy and are used to predict the properties of solid solutions without short-range order. Such potentials are usually averaged from an existing interatomic potential and are therefore only as accurate as the parent model. Accurate interatomic potentials are themselves difficult to parameterize and can require large training datasets. Here we benchmark a recently developed formalism that computes an exact A-atom potential directly from a linear atomic cluster expansion (ACE). We first fit a linear ACE to Fe-W data generated with an embedded atom method (EAM) potential. The resulting A-atom potential reproduces the properties of the disordered phase computed from an explicit random supercell and from a conventional A-atom potential averaged from the same EAM potential. We then fit a linear ACE to energies and forces computed from electronic structure calculations for about 1500 small Mo-Nb structures with an average of 7 atoms per structure. The A-atom potential derived from this ACE reproduces the DFT elastic constants, lattice parameter, mixing enthalpy, and Bain path of special quasirandom structures. Of the three chemical site bases, only the occupation basis also reproduces the DFT surface energies of the alloy. For Mo-Nb, the A-atom potential also predicts that the ideal solution entropy outweighs the destabilizing vibrational contribution to the finite-temperature free energy of the disordered phase. These benchmarks show that the properties of disordered alloys can be recovered from small training datasets when the chemical site basis is chosen carefully.

Understanding a nonlinear dynamical system from time series requires not only reproducing its trajectories, but also identifying a simple representation that preserves its essential dynamical structure. Almost-linear recurrent neural networks (AL-RNNs) are piecewise-linear RNNs in which only a subset of units use ReLU nonlinearities, so that nonlinear capacity is explicitly controlled by the number of ReLU units. Their activation patterns define linear regions, represented as symbols, whose observed transitions form a symbolic transition graph. However, directly training AL-RNNs with few ReLU units to realize minimal dynamical representations can be unreliable. We ask whether an AL-RNN with more ReLU units can instead be trained first and systematically reduced to a minimal dynamical representation. We introduce a fixed-point-guided hierarchical reduction procedure that progressively linearizes selected ReLU units, merging neighboring linear regions and graph nodes while preserving distinct symbols containing fixed points (FPs). The resulting reduction tree defines a hierarchy of progressively simpler candidates. Each reduced candidate is initialized from the parent parameters and retrained under guidance from the parent dynamics. We also prove that reproducing $Q$ distinct fixed points requires at least $Q$ FP-containing symbols, providing a certificate of symbol-level minimality when this bound is attained. On the 3-scroll Chua system, direct training with the theoretical minimum of three ReLU units achieves high-fidelity minimal realizations in only 20% of seeds, whereas our learn-reduce-retrain strategy increases the seed-macro success rate to approximately 71% at the same final nonlinear capacity. These results show that redundant nonlinear capacity can serve as a scaffold for discovering and realizing minimal dynamical representations.

Non-abelian infrared radiation can be seen as both a practical ingredient for QCD phenomenology and a window on the long-distance structure of gauge theories. During the last decade, a novel viewpoint was developed on infrared phenomena: soft currents are seen as arising from asymptotic symmetries realised on the celestial sphere, and soft colour correlations similarly emerge from a two-dimensional conformal field theory on the sphere. In this contribution, we briefly review the status of this approach. We begin with the standard bulk picture, presenting the infrared factorisation of massless scattering amplitudes, and then we show how universal infrared information gets mapped to marked points on the celestial sphere. In particular, soft colour-dipole singularities are reproduced to all orders by a correlator of vertex operators in a free-boson theory on the sphere. We then proceed to study soft-gluon currents, and show how bulk results constrain the structure of the celestial theory beyond the free-field limit. For single emission, we note that all logarithms in dipole-like loop contributions can be absorbed into the choice of scale for the strong coupling. For multiple emission, we show how the emerging mixed-helicity celestial operator product expansion must be weighted by energy fractions, breaking holomorphic factorisation and obstructing associativity under successive collinear limits. These results illustrate a remarkable connection between the familiar physics of infrared QCD and the dynamics of a two-dimensional theory on the celestial sphere.

I present a greatly improved version of SOLARPROP, an open source framework for building models of charge-sign dependent solar modulation of cosmic rays in the heliosphere. The transport equation of charged particles is solved by converting it to a set of stochastic differential equations (SDEs). The code now includes the full set of SDEs in two spatial dimensions, allowing for magnetic fields with latitudinal components, and it has been generalized to accommodate models in three spatial dimensions. Due to code parallelization and optimization, a speedup of up to three orders of magnitude with respect to the original version is achieved. Python bindings have been added for ease of use. The code has been validated against the results of several model calculations found in the literature. While excellent agreement is found in many cases, important differences are observed for two model implementations. Possible root causes for the differences have been investigated.

Ultrashort light pulses can drive transient photocurrents in solids, the basis of modern terahertz (THz) emitters. Ferroelectrics have recently emerged in this context as sources that require neither the bias voltage of photoconductive antennas nor the magnetic field of spintronic emitters, yet the ultrafast photocurrent generating the emission has remained debated. Here, we resolve this question in epitaxial BiFeO$_3$ thin films with engineered ferroelectric domain configurations, from single-domain to periodic stripe patterns. The stripe-domain film emits a THz field four times stronger than the single-domain configurations, reaching about half the output of a metallic spintronic reference emitter while remaining unsaturated at our highest excitation fluences. Resolving the emitted THz field in amplitude and phase as a function of the pump light polarization and the crystal orientation, we separate the coexisting photocurrents by their symmetry. This identifies the ultrafast screening of the ferroelectric polarization by photocarriers as the dominant emission mechanism in the in-plane-polarized films, enhanced in multidomain films by the strong built-in fields at the domain walls, while an out-of-plane polarization favors the bulk photovoltaic current. Our results establish the ferroelectric domain configuration as a design parameter setting both the strength and the microscopic origin of light-driven THz emission.

We study the impact of higher-order QCD corrections on Higgs production via weak-boson fusion with subsequent decays of the Higgs boson to bottom-quark pairs. We consider realistic fiducial cuts and observe large corrections in a fixed-order calculation. They arise from the interplay of the particular fiducial cuts considered and the radiative corrections to the Higgs decay subprocess. Switching to a parton-shower-matched description of the decay captures a significant portion of the corrections already at the lowest order and improves the perturbative convergence of the calculation.