Publications

Order by: publication date number of citations

• Emergent Kitaev materials in synthetic Fermi-Hubbard bilayers

  Daniel González-Cuadra, Alejandro Bermudez

  SciPost Phys. 19, 130 (2025) · published 18 November 2025 | Toggle abstract · pdf

  We investigate the emergence of bond-directional spin-spin interactions in a synthetic Fermi-Hubbard bilayer that can be realized with ultracold fermions in Raman optical lattices. The model exploits synthetic dimensions to couple two honeycomb layers, each corresponding to a different hyperfine atomic state, via Raman-assisted tunneling and, moreover, via an inter-layer Hubbard repulsion due to the cold-atom scattering. In the strong-coupling regime at half filling, we derive effective spin Hamiltonians for the kinetic exchange featuring Kitaev, Heisenberg, off-diagonal exchange ($\Gamma$-couplings), as well as tunable Dzyaloshinskii-Moriya interactions. We identify specific configurations that generate both ferromagnetic and antiferromagnetic Kitaev couplings with various perturbations of relevance to Kitaev materials, providing a tunable platform that can explore how quantum spin liquids emerge from itinerant fermion systems. We analyze the Fermi-liquid and Mott-insulating phases, highlighting a correspondence between Dirac and Majorana quasi-particles, with possible phase transitions thereof. In an extreme anisotropic limit, we show that the model reduces to an inter-layer ribbon in a quasi-1D ladder, allowing for a numerical study of the correlated ground state using matrix product states. We find a transition from a symmetry-protected topological insulator to a Kitaev-like regime characterized by nonlocal string order. Our results establish that cold-atom quantum simulators based on Raman optical lattices can be a playground for extended Kitaev models, bridging itinerant fermionic systems and spin-liquid physics.

• A real-time approach to frequency-mixing spectroscopies: Application to sum and difference frequency generation in two-dimensional crystals

  Mike N. Pionteck, Myrta Grüning, Simone Sanna, Claudio Attaccalite

  SciPost Phys. 19, 129 (2025) · published 18 November 2025 | Toggle abstract · pdf

  We propose a computational framework to extract nonlinear response functions from real-time simulations in the presence of more than one external field. We apply this approach to the calculation of sum frequency generation (SFG) and difference frequency generation (DFG). SFG and DFG are second-order nonlinear processes where two lasers with frequencies $\omega_1$ and $\omega_2$ combine to produce a response at frequency $\omega = \omega_1 ± \omega_2$. Compared with other nonlinear responses such as second-harmonic generation, SFG and DFG allow for tunability over a larger range. Moreover, the optical response can be enhanced by selecting the two laser frequencies in order to match specific electron-hole transitions. To assess the approach, we calculate the SFG and DFG of two-dimensional crystals, h-BN and MoS$_2$ monolayers, from real-time solution of an effective Schrödinger equation. Within the effective Schrödinger equation, one can select from various levels of theory for the effective one-particle Hamiltonian to account for local-field effects and electron-hole interactions. We compare results obtained within the independent particle picture and including many-body effects. Such comparison allows us to identify and characterize excitonic features in the obtained spectra. Additionally, we demonstrate that our approach can also extract higher-order response functions, such as field-induced second-harmonic generation. We provide an example using the h-BN bilayer.

  1 citation

• Three dimensional topological field theories and Nahm sum formulas

  Dongmin Gang, Heeyeon Kim, Byoungyoon Park, Spencer Stubbs

  SciPost Phys. 19, 128 (2025) · published 12 November 2025 | Toggle abstract · pdf

  It is known that a large class of characters of 2d conformal field theories (CFTs) can be written in the form of a Nahm sum. In [D. Zagier, in Frontiers in number theory, physics, and geometry II, Springer, Berlin (2007)], D. Zagier identified a list of Nahm sum expressions that are modular functions under a congruence subgroup of $SL(2,\mathbb{Z})$, which can be thought of as candidates for characters of rational CFTs. Motivated by the observation that the same formulas appear as the half-indices of certain 3d $\mathcal{N}=2$ supersymmetric gauge theories, we perform a general search over low-rank 3d $\mathcal{N}=2$ Abelian Chern-Simons matter theories which either flow to unitary TFTs or $\mathcal{N}=4$ rank-zero SCFTs in the infrared. These are exceptional classes of 3d theories, which are expected to support rational and $C_2$-cofinite chiral algebras on their boundary. We compare and contrast our results with Zagier's and comment on a possible generalization of Nahm's conjecture.

  1 citation

• Testing Kubo formula on a nonlinear quantum conductor driven far from equilibrium, via power exchanges

  Zubair Iftikhar, Jonas Müller, Yuri Mukharsky, Philippe Joyez, Patrice Roche, Carles Altimiras

  SciPost Phys. 19, 127 (2025) · published 11 November 2025 | Toggle abstract · pdf

  We present an experimental test of Kubo formula performed on a nonlinear quantum conductor, a Superconductor-Insulator-Superconductor tunnel junction, driven far from equilibrium by a DC voltage bias. We implement the proposal of Lesovik & Loosen [G. B. Lesovik and R. Loosen, J. Exp. Theor. Phys. Lett. 65, 295 (1997)] and demonstrate experimentally that it is possible to extract both the emission and absorption noise of the conductor by measuring the power it exchanges with a linear detection circuit whose occupation is tuned close to vacuum levels. We then compare their difference to the real part of the admittance which is independently measured by coherent reflectometry, finding that Kubo formula holds within experimental accuracy. Last, we show theoretically that the spectral density of power exchanged between a quantum conductor and its linear detection circuit follows a Lesovik & Loosen like formula, even in the presence of strong detection back-action. This result applies as long as the conductor acts as a current source for the detection circuit and the detection circuit is not singular.

• Non-invertible symmetries of two-dimensional non-linear sigma models

  Guillermo Arias-Tamargo, Chris Hull, Maxwell L. Velásquez Cotini Hutt

  SciPost Phys. 19, 126 (2025) · published 11 November 2025 | Toggle abstract · pdf

  Global symmetries can be generalised to transformations generated by topological operators, including cases in which the topological operator does not have an inverse. A family of such topological operators are intimately related to dualities via the procedure of half-space gauging. In this work we discuss the construction of non-invertible defects based on T-duality in two dimensions, generalising the well-known case of the free compact boson to any Non-Linear Sigma Model with Wess-Zumino term which is T-dualisable. This requires that the target space has an isometry with compact orbits that acts without fixed points. Our approach allows us to include target spaces without non-trivial 1-cycles, does not require the NLSM to be conformal, and when it is conformal it does not need to be rational; moreover, it highlights the microscopic origin of the topological terms that are responsible for the non-invertibility of the defect. An interesting class of examples are Wess-Zumino-Witten models, which are self-dual under a discrete gauging of a subgroup of the isometry symmetry and so host a topological defect line with Tambara-Yamagami fusion. Along the way, we discuss how the usual 0-form symmetries match across T-dual models in target spaces without 1-cycles, and how global obstructions can prevent locally conserved currents from giving rise to topological operators.

  1 citation

• Characterizing the hadronization of parton showers using the HOMER method

  Benoît Assi, Christian Bierlich, Philip Ilten, Tony Menzo, Stephen Mrenna, Manuel Szewc, Michael K. Wilkinson, Ahmed Youssef, Jure Zupan

  SciPost Phys. 19, 125 (2025) · published 11 November 2025 | Toggle abstract · pdf

  We update the HOMER method, a technique to solve a restricted version of the inverse problem of hadronization – extracting the Lund string fragmentation function $f(z)$ from data using only observable information. Here, we demonstrate its utility by extracting $f(z)$ from synthetic Pythia simulations using high-level observables constructed on an event-by-event basis, such as multiplicities and shape variables. Four cases of increasing complexity are considered, corresponding to $e^+e^-$ collisions at a center-of-mass energy of $90≥v$ producing either a string stretched between a $q$ and $\bar{q}$ containing no gluons; the same string containing one gluon $g$ with fixed kinematics; the same but the gluon has varying kinematics; and the most realistic case, strings with an unrestricted number of gluons that is the end-result of a parton shower. We demonstrate the extraction of $f(z)$ in each case, with the result of only a relatively modest degradation in performance of the HOMER method with the increased complexity of the string system.

• Bootstrapping nonequilibrium stochastic processes

  Minjae Cho

  SciPost Phys. 19, 124 (2025) · published 10 November 2025 | Toggle abstract · pdf

  We show that bootstrap methods based on the positivity of probability measures provide a systematic framework for studying both synchronous and asynchronous nonequilibrium stochastic processes on infinite lattices. First, we formulate linear programming problems that use positivity and invariance property of invariant measures to derive rigorous bounds on their expectation values. Second, for time evolution in asynchronous processes, we exploit the master equation along with positivity and initial conditions to construct linear and semidefinite programming problems that yield bounds on expectation values at both short and late times. We illustrate both approaches using two canonical examples: the contact process in 1+1 and 2+1 dimensions, and the Domany-Kinzel model in both synchronous and asynchronous forms in 1+1 dimensions. Our bounds on invariant measures yield rigorous lower bounds on critical rates, while those on time evolutions provide two-sided bounds on the half-life of the infection density and the temporal correlation length in the subcritical phase.

• Interpolating amplitudes

  Víctor Bresó-Pla, Gudrun Heinrich, Vitaly Magerya, Anton Olsson

  SciPost Phys. 19, 123 (2025) · published 10 November 2025 | Toggle abstract · pdf

  The calculation of scattering amplitudes at higher orders in perturbation theory has reached a high degree of maturity. However, their usage to produce physical predictions within Monte Carlo programs is often precluded by the slow evaluation of two- and higher-loop virtual amplitudes, particularly those calculated numerically. As a remedy, interpolation frameworks have been successfully used for amplitudes depending on up to two kinematic invariants. For amplitude interpolation with more variables, such as the five dimensions of a $2\to 3$ phase space, efficient and reliable solutions are sparse. This work aims to pave the way for using amplitude interpolation in higher-dimensional phase spaces by reviewing state-of-the-art interpolation methods, and assessing their performance on a selection of $2\to 3$ scattering amplitudes. Specifically, we investigate interpolation methods based on polynomials, splines, spatially adaptive sparse grids, and neural networks (multilayer perceptron and Lorentz-Equivariant Geometric Algebra Transformer), all under the constraint of limited obtainable data. Our additional aim is to motivate further studies of the interpolation of scattering amplitudes among both physicists and mathematicians.

• Geometric expansion of fluctuations and average shadows

  Clément Berthière, Benoit Estienne, Jean-Marie Stéphan, William Witczak-Krempa

  SciPost Phys. 19, 122 (2025) · published 7 November 2025 | Toggle abstract · pdf

  Fluctuations of observables provide unique insights into the nature of physical systems, and their study stands as a cornerstone of both theoretical and experimental science. Generalized fluctuations, or cumulants, provide information beyond the mean and variance of an observable. In this paper, we develop a systematic method to determine the asymptotic behavior of cumulants of local observables as the region becomes large. Our analysis reveals that the expansion is closely tied to the geometric characteristics of the region and its boundary, with coefficients given by convex moments of the connected correlation function: the latter is integrated against intrinsic volumes of convex polytopes built from the coordinates, which can be interpreted as average shadows. A particular application of our method shows that, in two dimensions, the leading behavior of odd cumulants of conserved quantities is topological, specifically depending on the Euler characteristic of the region. We illustrate these results with the paradigmatic strongly-interacting system of two-dimensional quantum Hall state at filling fraction $1/2$, by performing Monte-Carlo calculations of the skewness (third cumulant) of particle number in the Laughlin state.

• Automated evaluation of imaginary time strong coupling diagrams by sum-of-exponentials hybridization fitting

  Zhen Huang, Denis Golež, Hugo U. R. Strand, Jason Kaye

  SciPost Phys. 19, 121 (2025) · published 7 November 2025 | Toggle abstract · pdf

  We present an efficient separation of variables algorithm for the evaluation of imaginary time Feynman diagrams appearing in the bold pseudo-particle strong coupling expansion of the Anderson impurity model. The algorithm uses a fitting method based on AAA rational approximation and numerical optimization to obtain a sum-of-exponentials expansion of the hybridization function, which is then used to decompose the diagrams. A diagrammatic formulation of the algorithm leads to an automated procedure for diagrams of arbitrary order and topology. We also present methods of stabilizing the self-consistent solution of the pseudo-particle Dyson equation. The result is a low-cost and high-order accurate impurity solver for quantum embedding methods using general multi-orbital hybridization functions at low temperatures, appropriate for low-to-intermediate expansion orders. In addition to other benchmark examples, we use our solver to perform a dynamical mean-field theory study of a minimal model of the strongly correlated compound Ca$_2$RuO$_4$, describing the anti-ferromagnetic transition and the in- and out-of-plane anisotropy induced by spin-orbit coupling.