Papers

We construct an explicit binary linear system game that separates $C_{qa}$, the closure of the set of finite dimensional quantum correlations, from $C_{qc}$, the set of commuting operator correlations. The game admits a perfect commuting operator strategy, while every correlation in $C_{qa}$ has a success probability strictly less than 1. This provides a concrete counterexample to Tsirelson's problem in its approximation form. The defining system has 1417152 equations in 1889684 variables, with exactly three nonzero coefficients per equation and a single nonzero entry on the right hand side. We also compute the classical value exactly. The complete system is specified in Lean 4, and the separation and exact classical value are formalized using Mathlib.

Strong gravitational lensing provides a unique tool of studying intrinsically faint high-redshift active galactic nuclei (AGN). In particular, low-luminosity AGN (LLAGN) can probe black-hole growth in low-mass galaxies that current X-ray surveys cannot otherwise access. We investigate the nature and physical properties of a faint X-ray source triply imaged by the galaxy cluster SDSS J2243$-$0935, and determine whether the source hosts an AGN. We combine archival Hubble Space Telescope imaging and Chandra X-ray observations with the published strong-lensing models of SDSS J2243$-$0935. We constrain the source redshift geometrically, derive magnification maps, and infer the physical properties of the host galaxy. The geometric lensing analysis gives a source redshift of $z=2.14^{+0.12}_{-0.17}$, consistent with the photometric redshift estimate, and a total magnification of $\mu=78\pm7$. The source is triply lensed and detected only in the HST F160W band, with an intrinsic magnitude of $m_{\rm AB}=28.2\pm0.2$. The rest-frame optical emission implies a stellar mass of $8.6\lesssim\log(M_\star/M_\odot)\lesssim9.4$. The unobscured star formation rate is constrained to $\mathrm{SFR}_{\rm UV}<0.009\,M_\odot\,\mathrm{yr}^{-1}$, although substantial dust obscuration remains possible. Chandra detects X-ray point-source counterparts to all three HST-detected lensed images. The combined X-ray spectrum contains $100\pm11$ net counts and is consistent with a moderately absorbed power law, with $\Gamma \simeq 1.48$ and $N_{\rm H}=(5.5\pm4.4)\times10^{22}\,\mathrm{cm}^{-2}$. After correcting for gravitational magnification, we obtain an intrinsic rest-frame $2$–$10$ keV luminosity of $L_{\rm X}=1.2\pm0.3\times10^{42}\,\mathrm{erg\,s^{-1}}$. The X-ray emission is likely dominated by an AGN hosted by a low-mass galaxy at $z\simeq2.14$, making this system a strongly lensed LLAGN at cosmic noon (abridged).

Scientific dynamics forecasting is often framed as an architecture choice, although deployment is also determined by observed history, rollout feedback, compute budget, physical objective, and test distribution. We formulate protocol-dependent model selection and introduce ProtocolMatch, a compute-matched, validation-selected, and failure-preserving evaluation framework. On driven quantum-spin dynamics, we compare recurrent, patched-attention, causal-attention, and low-rank linear predictors across three independently generated datasets. The causal-attention–recurrence ordering reverses as the training set grows within a fixed two-spin task, while a linear predictor has the lowest mean error in the six-spin local-observable comparison. Restricting observed history worsens every refreshed-history view but improves every closed-loop view in the four-spin study. A latest-state MLP has lower error than persistence on every dataset under state refresh across all five cells, yet its closed-loop rank varies by system and includes finite explosive errors. Physical penalties improve targeted consistency without reliably improving prediction error, and in-distribution intervals lose most coverage after a driving-frequency shift. Thus scientific model selection should return a predictor with its protocol and report accuracy, physical validity, and shifted-distribution reliability separately.

Photometric decomposition is often used to study galaxy structure by modelling surface brightness distribution in its components. However, most of the relevant codes are designed for single-band analysis. I present IMFIT-MW: a lightweight, open source Python tool that extends IMFIT capabilities to simultaneously model multiwavelength images. IMFIT-MW is able to set shared parameters between filters and to impose smooth variation of them with wavelength. IMFIT-MW achieves this by using IMFIT to produce model brightness distributions for given parameters, whereas hyperparameters defining their wavelength dependence are fitted by an external optimizer. The code is available on GitHub (https://github.com/IVChugunov/IMFIT-MW). I provide a few examples of use, showing that IMFIT-MW matches GALFITM in single-thread performance and can have advantage of multithreading. In a broader sense, I demonstrate the possibility to extend IMFIT features even without modifying its source code.

We investigate the dynamics of quantum ergotropy in a multipartite system of two-level atoms interacting with a quantized field. The role of multipartite quantum mutual information and quantum relative entropy is analyzed to understand how correlations and the distinguishability between the passive state and the product of local states affect the extractable work. Numerical results for systems containing two to five atoms show that the maximum ergotropy does not increase linearly with the system size, while the intervals of zero ergotropy are gradually modified and suppressed for larger systems. We also study the effect of the average thermal photon number and different initial atomic states. Increasing the thermal photon number mainly reduces the relative-entropy contribution, whereas the initial state strongly influences both the magnitude and persistence of ergotropy. Among the considered states, the partially entangled state provides the largest and more persistent extractable work, while the mixed GHZ state gives the smallest ergotropy. These results highlight the importance of multipartite correlations, thermal effects, and state preparation in controlling quantum work extraction and may be useful for future studies of quantum batteries and quantum power.

In the Minimal Supersymmetric Standard Model with heavy sfermions, a sub-hundred-GeV Bino-like neutralino can thermally achieve the observed dark matter abundance only through $h$- or $Z$-resonant annihilation. Since the couplings governing these resonances simultaneously determine the spin-independent and spin-dependent scattering off nucleons, reproducing the relic density imposes irreducible lower envelopes on the scattering cross sections. We evaluate these bounds including the complete one-loop corrections to these vertices, which enhance the direct-detection rates by up to $50\%$ and $30\%$, respectively, adopt the FLAG-2024 nucleon matrix elements with their most adverse $3\sigma$ shifts, and confront the predictions with a likelihood-based recast of the LZ data. We find that the LZ-2024 results exclude both resonant scenarios at conservative one-sided significances exceeding $5.9\sigma$ and $4.1\sigma$, respectively, irrespective of whether the neutralino constitutes all or only part of the dark matter, and independently of LHC electroweakino searches.

Quantum sources of light are of paramount importance for the development of photonic quantum technologies. The growing demand for bright, efficient, and spectrally engineered photonic sources for applications in quantum sensing, spectroscopy, and imaging has motivated the development of tailored nonlinear crystals. Here, we report on the development and characterization of chirped, apodized KTP crystals specifically engineered for quantum spectroscopy. Our source exhibits an unprecedented bandwidth in a highly non-degenerate configuration, demonstrating its potential to support quantum spectroscopy protocols over broad spectral ranges while enabling significantly shorter measurement times.

We present the angular BAO analysis using DESI DR1 data as a complementary approach to DESI's primary 3D BAO analysis. We perform the analysis using the two-point angular correlation function $w(\theta)$ in tomographic redshift shells of $\delta z = 0.02, 0.05$, spanning the redshift range $0.4 < z < 2.1$ using the DESI tracers: LRG1, LRG2, LRG3+ELG1, ELG2 and QSO. The main benefit is that it requires no fiducial cosmology to convert redshifts into distances and yields a direct measurement of the transverse dilation parameter $\alpha_{\perp}(z)$. We fit each tracer with a single shared $\alpha_{\perp}(z)$ using templates constructed in a fiducial cosmology and a mock-based covariance from 1000 EZmock realisations. We measure $\alpha_\perp$ to $2.5$ – $4.2\%$ precision in each tracer bin. Compared to the $D_M/r_{\rm d}$ measurements obtained from the DESI DR1 3D analysis for the same tracers, our results are consistent within the statistical precision. We find no evidence for a shift in the transverse BAO scale between the 2D and 3D analyses. Flat $\Lambda$CDM fits to $\alpha_\perp(z)$ alone give $\Omega_{\rm m} = 0.198^{+0.045}_{-0.068}$, $1.6\sigma$ below the DESI 3D value, and the pull is dominated by the $z =1.49$ QSO bin, whose $\alpha_\perp(z=1.49)=1.058\pm0.028$ sits $\sim 2\sigma$ above the fiducial. Restricting to LRGs and ELGs only, our $\Omega_{\rm m}$ and $H_0 r_{\rm d}$ agree with the 3D transverse-only fit over the same bins to better than $0.1\sigma$. With an additional BBN prior, the full 2D sample gives $H_0 = 63.9^{+1.9}_{-3.3}\, {\rm km\,s^{-1}\,Mpc^{-1}}$. This pipeline is implementable on forthcoming DESI data releases, where the increase in volume and sampling density will determine if the QSO offset persists and will also tighten 2D BAO constraints.

Transient supplies of nutrients to the surface ocean, both natural and artificial, stimulate blooms of phytoplankton and the formation of organic matter, driving air-sea gradients and uptake of CO$_2$. However, quantifying the associated carbon budget remains challenging, as it requires tracking the coupled biophysical evolution of water masses while they are transported, stretched and diluted. Here, we present an idealized three-dimensional model that describes biomass production and carbon dynamics within a Lagrangian patch in the ocean. The framework reproduces observed biogeochemical patterns from an artificial fertilization experiment and provides integrated metrics for the local carbon budget. Utilizing large ensembles of simulations, we examine the sensitivity of patch-scale primary production and carbon uptake to biochemical and physical factors. Our results show that patch dilution can enhance the ecosystem response, and how carbon uptake is sensitive to initial injected area, horizontal divergence and vertical diffusivity. In the context of renewed interest in ocean fertilization strategies for climate mitigation, our approach can thus provide quantitative tools to assess their efficacy and potential.

We construct chiral extensions of the cotangent Yangian using vertex algebroids. At zero level, the construction yields a vertex algebra with a compatible coproduct, whose Zhu reduction recovers the cotangent Yangian of Abedin and Niu. We also construct vertex algebra extensions at non-zero levels and study level-additive coproducts. The zero-level vertex algebra admits two natural quasi-conformal structures. We compute the corresponding sphere one-point conformal blocks and equip each block space with an associative product induced by the coproduct. For the first quasi-conformal structure, the blocks space can be identified with the dual of a certain Lie algebroid homology; for the second, we describe the block algebra and its classical limit in terms of mini-twistor geometry.

Recent cosmological constraints on the sum of neutrino masses have reached, and in some analyses fallen below, the minimum value implied by the normal mass ordering, suggesting a possible tension with neutrino-oscillation measurements. We investigate whether this discrepancy can be alleviated in a mass-varying-neutrino (MaVaN) cosmology in which a canonical pseudo-Nambu-Goldstone boson (pNGB) drives late-time cosmic acceleration and controls the masses of massive neutrinos through an exponential coupling. Using current cosmic microwave background, baryon acoustic oscillation, and type Ia supernova data, we reconstruct the neutrino-mass evolution and show that the cosmologically inferred effective neutrino mass can fall below the present-day normal-ordering floor while remaining consistent with oscillation constraints. The same interaction produces a distinctive late-time signature, with 96% of the posterior expansion histories exhibiting a crossing of the apparent phantom divide, $w_{\rm app}=-1$, at $z_{\rm cross}=0.90^{+0.28}_{-0.13}$ at 68% C.L. The underlying scalar field nevertheless remains canonical, with $w_\phi \geq -1$. Our results therefore show that neutrino mass variation can simultaneously alleviate the apparent neutrino-mass tension and generate the phantom-crossing behavior favored by current cosmological data, without introducing a phantom field or violating the null energy condition.

A comprehensive (start-to-end) simulation was performed regarding the formation of a front-tilted laser pulse and the subsequent acceleration of non-relativistic electrons moving along a reflective grating structure. The laser pulse tilt was formed using another diffraction grating – either reflective or transmissive. The pulse had a spectral range (FWHM) of 803–863 nm and an initial duration of approximately 17 fs. Acceleration under normal laser pulse incidence was also calculated for comparison. Research results indicate that the use of a femtosecond laser pulse with a front, tilted by means of a reflective diffraction grating, significantly enhances the efficiency of laser acceleration of non-relativistic electrons compared to the case of a pulse at normal incidence. Furthermore, it was found that using a transmissive diffraction grating to shape the tilted pulse front also increases acceleration efficiency relative to normal incidence on the accelerating grating, although the acceleration rate is lower than when a reflective diffraction grating is used.

The radiation-ionic effect was recently reported in 3 mol% Gd-doped CeO$_2$ (3GDC): under $^{60}$Co gamma irradiation at 35 Gy/min, the grain-boundary (GB) ionic resistance of a bulk pellet decreased 680-fold at room temperature, because electrons trapped at the GB core lower the space-charge potential. The effect has been demonstrated only with photons. Gd has the largest thermal-neutron capture cross section of any stable element and is also the dopant that creates the GB barrier in GDC, so the same mechanism could link GB impedance to thermal-neutron sensing; whether it does, and under what conditions, has not been examined. Here the thermal-neutron response of 3GDC is modelled with the gamma-calibrated Mott-Schottky/Seager framework, with unmeasured parameters carried as intervals and propagated jointly. Each capture deposits about 50 keV by conversion and Auger electrons within 1 $\mu$m, the remaining 8 MeV leaving as prompt photons; Gd segregation to the GBs adds a factor of only 1-2 because the electron range exceeds the grain size. A thermal flux of $8\times10^{9}$ cm$^{-2}$s$^{-1}$ (6-12$\times10^{9}$) matches the GB-effective generation rate of the 35 Gy/min calibration, and the predicted room-temperature GB resistance ratio is 110 (70-170) at $10^{9}$ and 970 (600-1500) at $10^{10}$ cm$^{-2}$s$^{-1}$, falling to 2 at 61 (55-68) and 81 (73-92) $^\circ$C. The neutron case differs from gamma irradiation in three respects: capture events are discrete, so the continuous-source approximation depends on the unmeasured trapped-charge lifetime; self-shielding (mean free path 0.27 mm) makes the response thickness dependent; and capture photons absorbed in the pellet add a gamma term comparable to the local term for pellets of 0.3 mm or more. The intervals are conditional on the model and exclude its structural error. The conditions under which the prediction can be tested are given.

Paying particular attention to the different types of gluons, collinear, soft and Glauber, we argue that, in its standard formulation, the CGC does not contain the collinear sector required for a meaningful comparison with TMD factorization. The CGC Wilson lines usually identified as part of the proton matrix element should instead be interpreted as the equivalent of a soft factor. Restoring the missing collinear exchanges, we obtain a new definition of small-$x$ PDFs from an explicit calculation of the DIS cross section using CGC Feynman rules. The resulting expression can be cast in a form resembling TMD factorization without being equivalent to it. Finally, we discuss the evolution of these small-$x$ PDFs with respect to the rapidity cutoff $\Lambda$ and the Bjorken variable $x_b$.

Mediators coupled predominantly to the electron remain weakly constrained above the kinematic reach of the $B$~factories, a region that is difficult to access at hadron colliders and only partially covered by existing $e^+e^-$ data. We point out that the asymmetric $\mu^+e^-$ mode of the proposed $\mu$TRISTAN collider, operating at $\sqrt{s}=346$~GeV, is well suited to this window, and quantify its reach for two minimal realisations: an axion-like particle and a $Z'$ gauge boson, treated as effective interactions with the electron current. Both are radiated off the electron line and decay to an electron-positron pair, giving a common four-lepton signature whose Standard Model counterpart is under good control. A resonance search in the electron-positron invariant mass, performed for mediator masses between $5.5$ and $300$~GeV, yields projected $95\%$~CL sensitivities to couplings of order $10^{-5}$ at the light end of this range, improving on existing bounds by more than four orders of magnitude for the pseudoscalar and by up to two orders of magnitude for the vector, and providing the only constraint on the heaviest masses considered.

The security of a secret sharing protocol is compromised if one or more dishonest participants can reconstruct the secret by deviating from the prescribed protocol, particularly when such malicious behavior remains undetected. Semi-Quantum Secret Sharing (SQSS) is a variant of secret sharing in which the participants possess only classical capabilities, such as preparing and measuring qubits in the computational ($Z$) basis, while a quantum-capable third party assists the dealer in generating and distributing quantum shares of a classical secret. The limited quantum capabilities of the participants in SQSS protocols may introduce vulnerabilities to collusion attacks, enabling an assisting quantum third party, in collaboration with a dishonest classical participant, to jointly recover the secret even without detection. In this work, we propose a novel SQSS protocol that eliminates this vulnerability and achieves information-theoretic security against collusion attacks. We further prove that the proposed protocol is secure against a broad class of external and internal attacks. Finally, a comparative analysis demonstrates that our construction advances existing SQSS protocols by simultaneously optimizing qubit efficiency, involves a classical dealer, and resilience against DCNA attacks using basic Bell states.

Optically exciting a typical low-dimensional semiconductor produces a nodeless \(s\)-type exciton as its lowest-energy bound state. Here, we explain how Coulomb phase matching stabilises a quantum-vortex exciton as the lowest bound state in flat Chern bands. Using a prototypical Yin–Yang kagome lattice, we demonstrate that increasing the spin–orbit coupling drives a transition of the lowest exciton from a quantum-vortex state to a zero-winding state, while leaving the band Chern numbers unchanged. We trace this behaviour to a gauge-invariant combination of exciton phase differences and Bloch-overlap phases that can reduce the Coulomb energy. By comparing states with identical wavefunction amplitudes, we isolate this phase contribution and show that it is sufficient to drive the reversal in exciton ordering. We classify exciton topology by the momentum-space wavefunction winding relative to the conduction–valence Chern-number difference,connecting the competing states to circular-polarisation selection rules. We then show that the resulting reordering can be tracked directly in the optical absorption spectrum, through the transition from a dark vortex exciton to a bright zero-winding exciton. These findings establish a microscopic link between exciton topology, Coulomb binding, and optical response that extends beyond the information contained in the underlying electronic band Chern numbers.

The stability of the Berezinskii–Kosterlitz–Thouless phase against long-range interactions remains an open problem. We address this question in the two-dimensional XY model with couplings decaying as $r^{-2-\sigma}$, by developing a real-space renormalization-group treatment near the Gaussian manifold of quasi-long-range ordered (QLRO) states. We retain the full many-body vortex interaction and show that, at leading order in the fugacity, screening by small vortex–antivortex pairs preserves its Gaussian-average representation without an expansion in the long-range coupling. Including spin waves perturbatively then yields a local stability criterion: a stable QLRO interval exists only for $\sigma>7/4$. We also investigate a Wilson–Fisher-like fixed point at finite coupling for $\sigma<7/4$, for which the retained flow predicts $\nu^{-1}\sim\sqrt{7/4-\sigma}$ as it merges with the Gaussian line.

We construct unitary, local, 3+1-dimensional conformal field theories with finitely many fields whose thermal state breaks a continuous global symmetry at every temperature. They are finite-$N$ versions of the Chaudhuri–Choi–Rabinovici theories, whose planar circle of fixed points has an ordered arc. The circle is lifted at order $1/N^2$, and an infinite family of fixed points on this arc reaches arbitrarily weak coupling, where the superfluid is under parametric control. We also show that weakly gauging the symmetry gives a superconductor over an exponential range of temperatures.

The lattice dynamics of disordered systems (amorphous solids, polymers and glasses) are often distinct from crystalline analogues. Examples include the "two-level systems" and the Boson peak. Here, we show that even the venerable material amorphous silicon offers new wrinkles: (1) realistic networks studied with DFT-quality interatomic potentials show the existence of harmonic modes exhibiting anomalously large root-mean-square variation in atomic positions, and (2) these dynamics determine the processes of heat transport. We correlate the atomic dynamics with transport using the site-projected thermal conductivity [C. Ugwumadu et. al Phys. Status Solidi B 263, e202500316 (2026)], and conjecture that these observations are relevant to other disordered systems.

The near-infrared (NIR) scintillation signal of xenon has the potential to enhance the performance of xenon-based detectors, but a robust understanding of its properties is a prerequisite for its application. In this work, we measure the NIR time response for $\alpha$-particle interactions in gaseous xenon in 12 wavelength bands spanning approximately $(950-1620)\,$nm. The time response is decomposed into a fast (ns-scale) component and a dominant slow ($\mu$s-scale) component. The fast-component spectrum coincides with strong atomic emission lines, predominantly 6p$\,\to\,$6s, whereas the slow-component spectrum follows the broad continuum centred around $1.3\,\mu$m. The pressure dependence of the continuum can be described by a two-step kinetic model in which the emitting state is formed by a three-body process and depopulated by a two-body collisional process, with density-independent terms in both rates. Interpreting the slower of these as radiative decay yields an intrinsic lifetime of $\sim3\,\mu$s. Theoretical calculations identify a neutral-excimer transition $0_u^- \rightarrow 0_g^-$ and an ionic-dimer transition $1/2_{u} \rightarrow 1/2_g$ as plausible candidates for the continuum emission, with the ionic-dimer transition showing better overall agreement with the measurements.

Non-equilibrium electron and phonon populations govern heat generation and dissipation in nanoscale devices, yet their direct characterization remains challenging. Anti-Stokes Raman scattering provides a sensitive probe of phonon populations and local temperatures, although its interpretation becomes more complex when electronic and phononic populations are driven out of equilibrium. Here, we investigate the gate-dependent anti-Stokes Raman response of monolayer graphene and find that the anti-Stokes response deviates from that expected from an equilibrium phonon population. Using a spatially localized microdroplet spectroelectrochemical platform, dielectric back-gates, and chemical doping, we observe a strong and reversible suppression of the G-mode anti-Stokes intensity as the Fermi level is tuned away from the Dirac point, while the Stokes intensity remains essentially unchanged. By contrast, the G-mode frequency and linewidth of both the Stokes and anti-Stokes branches evolve with doping as expected from the established Kohn-anomaly picture. The anti-Stokes/Stokes intensity ratio decreases by several-fold over the measured range and displays a pronounced, gate-dependent power response with maximum sensitivity at charge neutrality. Comparison with previous studies of non-equilibrium carrier dynamics in doped graphene suggests an important role of Fermi-level-dependent electronic phase space, which can influence both hot-carrier-driven phonon populations and electronically mediated Stokes-anti-Stokes scattering. Our findings establish anti-Stokes Raman scattering as a sensitive probe of electron-phonon coupling and non-equilibrium phonon dynamics in graphene.

In a recent paper, a mathematical analogy was found between two physically distinct classes of two-dimensional free boundary problems in inviscid fluid mechanics. One class involves free surface flows with constant vorticity in the bulk but without boundary surface tension, the other comprises free surface flows without vorticity in the bulk but with non-zero boundary surface tension (capillarity). Several examples were given of pairs of physically-distinct free boundary problems that amount to solving the same mathematical problem, albeit with the possibility of a single constraint arising for the capillary flows. The present work shows that this analogy extends to another pair of physically-distinct problems: a hollow vortex equilibrium in a shear flow with a third order irrotational strain and a steady bubble with surface tension in a simple straining flow. For both problems, we construct an isolated exact solution when a parameter $\gamma=0$ and provide numerical solutions for non-zero $\gamma$. An interesting feature is that non-trivial conformal mapping identities emerge from the class of problems with capillarity that apply equally for the hollow vortex problems when the parameters are appropriately constrained.

Any two orthogonal quantum states can be perfectly distinguished via local operations and classical communication (LOCC), yet whether all orthogonal multipartite pairs incur an identical operational cost remains an open question. Here, we address the resource consumption required for local state discrimination across multipartite pair of states. We demonstrate that distinct pairs of orthogonal multipartite maximally entangled states require unequal amounts of quantum resources to achieve local distinguishability. Extending our framework beyond individual pairs, we analyze larger sets of both locally distinguishable and indistinguishable multipartite orthogonal states, revealing a nontrivial hierarchy in their entanglement consumption. These results show that the resource cost of LOCC discrimination is state-dependent, providing insight into the fine-grained structure of multipartite nonlocality.

Artificial intelligence (AI) is becoming an increasingly useful tool across the experimental sciences, including angle-resolved photoemission spectroscopy (ARPES), which routinely produces large, multidimensional datasets of electronic structure. Recent advances in AI and machine learning (ML) have opened new opportunities across the entire ARPES workflow, from automated sample preparation and real-time data acquisition to post-experiment data analysis and comparison with theoretical calculations. Despite this progress, a comprehensive review of ML applications, their capabilities, and reliability across the different stages of ARPES workflow is still lacking. In this review, we first introduce ML methods that are most relevant to experimentalists working in condensed matter physics and materials science. We then follow the ARPES workflow, reviewing existing ML applications at each step and discussing their advantages, limitations and potential for future development. We also examine the current ARPES data landscape, where several open databases are available but remain relatively small and fragmented compared with large, shared datasets such as ImageNet. Given these limitations, we suggest that the community focus on sharing pretrained models that can be further trained, adapted to specific tasks, and redistributed, while working toward a larger and standardized open ARPES dataset repository. Finally, we discuss our perspectives on the future of AI within the ARPES workflow using a six-level framework of laboratory automation, highlighting the opportunities and challenges in moving toward a fully autonomous, self-driving ARPES laboratory.

Low-dimensional pnictogen chalcohalides (MChX; M = Bi, Sb; Ch = S, Se; X = I, Br) are promising non-toxic semiconductors for solar energy conversion, yet the power conversion efficiencies of MChX solar cells remain below 10%, far from their Shockley-Queisser limit of $\sim$30%. Efficient charge extraction depends critically on the band-edge alignment between the absorber and the charge-selective contacts, which is a surface property. Here, we combine first-principles calculations with device-level modelling to investigate the surface energetics and band alignments across the MChX family. An exhaustive sampling of crystal orientations, including low-symmetry facets, identifies the (011) and (010) surfaces as the most stable terminations in all eight compounds. Their formation energies differ by only $\approx 0.01$-$0.03\mathrm{J\,m^{-2}}$, yet their band edges are shifted rigidly with respect to each other by up to $0.8$eV. Drift-diffusion simulations of BiSBr show that the coexistence of such facets can reduce the open-circuit voltage by up to $0.6$V, even in an otherwise ideal, defect-free absorber. Surface orientation thus emerges as a hidden design parameter for charge extraction. This study may help to explain the gap between the theoretical and experimental efficiencies of MChX solar cells, and it provides new design principles for MChX-based photovoltaic and photocatalytic technologies.

Extreme debris disks are signposts of violent rocky collisions, but it remains unclear whether their mid-infrared variability traces transient dust clouds or longer-lived post-impact reservoirs. We analyze a 19 yr mid-infrared record of the young solar analog NGC 2547 ID8, combining four spectra from 2007-2024 with 3-5 micron monitoring through 2026. The optically thin mineral surface remains compositionally stable: broad fits spanning multiple grain sizes give 550-577 K, and a common-basis decomposition with T_surf fixed at 560 K is dominated by 1 micron olivine/Mg-olivine and 0.1 micron forsterite. The W2-equivalent 4.5 micron excess, traced across six recurrent high states, is instead driven mainly by changes in the hot-base continuum component and defines a 956.8 d infrared recurrence timescale with cycle-to-cycle variations in timing and amplitude. This combination supports the interpretation of ID8 as an early post-impact stage in which material on eccentric orbits repeatedly interacts with a warm rocky reservoir, preserving the impact's mineral fingerprint while the interaction geometry evolves. Material on eccentric orbits can persist for decades, and possibly longer, while interacting with debris disks. Longer-baseline, higher-cadence infrared monitoring of ID8 will further constrain the stability of its recurrence timescale and the detailed evolution of its infrared variability.

Transmission silicon detectors are widely used in $\Delta E$–$E$ telescope configurations for particle identification and timing in nuclear physics experiments. Variations in signal rise time can introduce time-walk effects that degrade time-of-flight measurements, as recently observed in light-ion production experiments at the Neutrons For Science (NFS) facility at GANIL using the Medley setup. This work investigates the influence of particle energy and detector non-uniformity on signal rise time.

Understanding the non-equilibrium dynamics of quantum many-body systems plays a central role in modern many-body physics. Intriguing insights emerge from investigating how quantum resources, which are essential for achieving unambiguous quantumness, are generated and spread under chaotic quantum dynamics. In fermionic systems, non-Gaussianity constitutes an important resource beyond entanglement, characterizing the deviation of generic quantum states from free-fermion states. In this Letter, we propose using the replica twist entropy, defined through the expectation value of a replica rotation operator with a tunable angle $\alpha$, as a natural probe of non-Gaussianity in pure fermionic states. After deriving its behavior for typical many-body states, we establish a general framework for analyzing the replica twist entropy in solvable Sachdev-Ye-Kitaev models in the large-$N$ limit. Applying this framework to the dynamics of the thermofield double state, we uncover a novel dynamical phase diagram in the parameter space of $\alpha$ and evolution time $t$, exhibiting spontaneous $Z_2$ symmetry breaking at $\alpha=\pi/4$ and universal cusps as a function of $\alpha$ in the long-time regime. Our results pave the way for studying non-Gaussianity in strongly correlated fermionic systems and reveal a new class of dynamical transitions probed by non-Gaussianity.

We use the MUSE-ALMA Haloes sample to build a census of the condensed baryons (stars and cold gas) to study the multiphase interstellar medium and circumgalactic medium (CGM). The sample includes 79 gas-rich galaxies at 0.2<z<1.2 located near strong HI Ly-alpha absorbers along quasar sightlines within projected separations of ~5-270 kpc and radial velocity separations of up to 500 km/s. HST spectroscopy characterises the absorbers, while HST imaging, MUSE and ALMA observations trace the stars, ionised gas and molecular gas. We find the following results. (i) HI masses inferred from optical properties show broad agreement with those derived from the CGM absorber properties for impact parameters <100 kpc, although group environments can significantly influence these estimates. (ii) The molecular gas fraction declines more steeply with stellar mass than the HI fraction. This suggests that stars build up more rapidly than the cold gas reservoir, while massive galaxies either consume their molecular gas more efficiently or convert HI into H2 less efficiently. (iii) Galaxies with high stellar surface densities exhibit the highest baryon fractions. This may reflect baryonic feedback, which can produce flatter inner dark matter profiles in high surface density systems, or aperture effects, where the ionised gas kinematic tracer probes only the innermost, baryon-dominated regions in compact sources, but extends further into the dark matter halo in sources with extended baryonic distributions. (iv) A few galaxies are consistent with being the baryon-dominated descendants of the low mass z~6 galaxies detected by JWST, while many appear to be substantially more dark matter dominated, lower mass systems with progenitors beyond JWST detection limits. This study provides a multiscale, multiphase view of the gas in and around galaxies, allowing us to probe the baryon budget across cosmic time.