Publications

Order by: publication date number of citations

• Chaos in the vicinity of a singularity in the three-body problem: The equilateral triangle experiment in the zero angular momentum limit

  Hugo D. Parischewsky, Alessandro A. Trani, Nathan W. C. Leigh

  SciPost Phys. Core 6, 016 (2023) · published 21 March 2023 | Toggle abstract · pdf

  We present numerical simulations of the gravitational three-body problem, in which three particles lie at rest close to the vertices of an equilateral triangle. In the unperturbed problem, the three particles fall towards the center of mass of the system to form a three-body collision, or singularity, where the particles overlap in space and time. By perturbing the initial positions of the particles, we are able to study chaos in the vicinity of the singularity. Here we cover both the singular region close to the unperturbed configuration and the binary-single scattering regime where one side of the triangle is very short compared to the other two. We make phase space plots to study the regular and ergodic subsets of our simulations and compare them with the outcomes expected from the statistical escape theory of the three-body problem. We further provide fits to the ergodic subset to characterize the properties of the left-over binaries. We identify the discrepancy between the statistical theory and the simulations in the regular subset of interactions, which only exhibits weak chaos. As we decrease the scale of the perturbations in the initial positions, the phase space becomes entirely dominated by regular interactions, according to our metric for chaos. Finally, we show the effect of general relativity corrections by simulating the same scenario with the inclusion of post-Newtonian corrections to the equations of motion.

  4 citations

• Testing the mechanism of lepton compositeness

  Vincenzo Afferrante, Axel Maas, René Sondenheimer, Pascal Törek

  SciPost Phys. 10, 062 (2021) · published 10 March 2021 | Toggle abstract · pdf

  Strict gauge invariance requires that physical left-handed leptons are actually bound states of the elementary left-handed lepton doublet and the Higgs field within the standard model. That they nonetheless behave almost like pure elementary particles is explained by the Fr\"ohlich-Morchio-Strocchi mechanism. Using lattice gauge theory, we test and confirm this mechanism for fermions. Though, due to the current inaccessibility of non-Abelian gauged Weyl fermions on the lattice, a model which contains vectorial leptons but which obeys all other relevant symmetries has been simulated.

  4 citations

• Quantum engineering of a synthetic thermal bath for bosonic atoms in a one-dimensional optical lattice via Markovian feedback control

  Ling-Na Wu, André Eckardt

  SciPost Phys. 13, 059 (2022) · published 21 September 2022 | Toggle abstract · pdf

  We propose and investigate a scheme for engineering a synthetic thermal bath for a bosonic quantum gas in a one-dimensional optical lattice based on Markovian feedback control. The performance of our scheme is quantified by the fidelity between the steady state of the system and the effective thermal state. For double-well and triple-well systems with non-interacting particles, the steady state is found to be an exact thermal state, which is attributed to the fact that the transfer rates between all pairs of coupled eigenstates satisfy detailed balance condition. The scenario changes when there are more lattice sites, where the detailed balance condition does not hold any more, but remains an accurate approximation. Remarkably, our scheme performs very well at low and high temperature regimes, with the fidelity close to one. The performance at the intermediate temperature regime (where a crossover into a Bose condensed regime occurs) is slightly worse, and the fidelity shows a gentle decrease with increasing system size. We also discuss the interacting cases. In contrast to the non-interacting cases, the scheme is found to perform better at a higher temperature. Another difference is that the minimal temperature that can be engineered is nonzero and increases with the interaction strength.

  4 citations

• A convenient Keldysh contour for thermodynamically consistent perturbative and semiclassical expansions

  Vasco Cavina, Sadeq S. Kadijani, Massimiliano Esposito, Thomas L. Schmidt

  SciPost Phys. 15, 209 (2023) · published 27 November 2023 | Toggle abstract · pdf

  The work fluctuation theorem (FT) is a symmetry connecting the moment generating functions (MGFs) of the work extracted in a given process and in its time-reversed counterpart. We show that, equivalently, the FT for work in isolated quantum systems can be expressed as an invariance property of a modified Keldysh contour. Modified contours can be used as starting points of perturbative and path integral approaches to quantum thermodynamics, as recently pointed out in the literature. After reviewing the derivation of the contour-based perturbation theory, we use the symmetry of the modified contour to show that the theory satisfies the FT at every order. Furthermore, we extend textbook diagrammatic techniques to the computation of work MGFs, showing that the contributions of the different Feynman diagrams can be added to obtain a general expression of the work statistics in terms of a sum of independent rescaled Poisson processes. In this context, the FT takes the form of a detailed balance condition linking every Feynman diagram with its time-reversed variant. In the second part, we study path integral approaches to the calculation of the MGF, and discuss how the arbitrariness in the choice of the contour impacts the final form of the path integral action. In particular, we show how using a symmetrized contour makes it possible to easily generalize the Keldysh rotation in the context of work statistics, a procedure paving the way to a semiclassical expansion of the work MGF. Furthermore, we use our results to discuss a generalization of the detailed balance conditions at the level of the quantum trajectories.

  4 citations

• Quantum-critical properties of the one- and two-dimensional random transverse-field Ising model from large-scale quantum Monte Carlo simulations

  Calvin Krämer, Jan Alexander Koziol, Anja Langheld, Max Hörmann, Kai Phillip Schmidt

  SciPost Phys. 17, 061 (2024) · published 23 August 2024 | Toggle abstract · pdf

  We study the ferromagnetic transverse-field Ising model with quenched disorder at $T = 0$ in one and two dimensions by means of stochastic series expansion quantum Monte Carlo simulations using a rigorous zero-temperature scheme. Using a sample-replication method and averaged Binder ratios, we determine the critical shift and width exponents $\nu_\mathrm{s}$ and $\nu_\mathrm{w}$ as well as unbiased critical points by finite-size scaling. Further, scaling of the disorder-averaged magnetisation at the critical point is used to determine the order-parameter critical exponent $\beta$ and the critical exponent $\nu_{\mathrm{av}}$ of the average correlation length. The dynamic scaling in the Griffiths phase is investigated by measuring the local susceptibility in the disordered phase and the dynamic exponent $z'$ is extracted. By applying various finite-size scaling protocols, we provide an extensive and comprehensive comparison between the different approaches on equal footing. The emphasis on effective zero-temperature simulations resolves several inconsistencies in existing literature.

  4 citations

• Carrollian $\mathscr Lw_{1+\infty}$ representation from twistor space

  Laura Donnay, Laurent Freidel, Yannick Herfray

  SciPost Phys. 17, 118 (2024) · published 22 October 2024 | Toggle abstract · pdf

  We construct an explicit realization of the action of the $\mathscr Lw_{1+∞}$ loop algebra on fields at null infinity. This action is directly derived by Penrose transform of the geometrical action of $\mathscr Lw_{1+∞}$ symmetries in twistor space, ensuring that it forms a representation of the algebra. Finally, we show that this action coincides with the canonical action of $\mathscr Lw_{1+∞}$ Noether charges on the asymptotic phase space.

  4 citations

• Kernel polynomial method for linear spin wave theory

  Harry Lane, Hao Zhang, David Dahlbom, Sam Quinn, Rolando D. Somma, Martin Mourigal, Cristian D. Batista, Kipton Barros

  SciPost Phys. 17, 145 (2024) · published 26 November 2024 | Toggle abstract · pdf

  Calculating dynamical spin correlations is essential for matching model magnetic exchange Hamiltonians to momentum-resolved spectroscopic measurements. A major numerical bottleneck is the diagonalization of the dynamical matrix, especially in systems with large magnetic unit cells, such as those with incommensurate magnetic structures or quenched disorder. In this paper, we demonstrate an efficient scheme based on the kernel polynomial method for calculating dynamical correlations of relevance to inelastic neutron scattering experiments. This method reduces the scaling of numerical cost from cubic to linear in the magnetic unit cell size.

  4 citations

• Codebase release 1.2 for YASTN

  Marek M. Rams, Gabriela Wójtowicz, Aritra Sinha, Juraj Hasik

  SciPost Phys. Codebases 52-r1.2 (2025) · published 26 February 2025 | Toggle abstract · src

  We present an open-source tensor network Python library for quantum many-body simulations. At its core is an Abelian-symmetric tensor, implemented as a sparse block structure managed by a logical layer on top of a dense multidimensional array backend. This serves as the basis for higher-level tensor network algorithms operating on matrix product states and projected entangled pair states. An appropriate backend, such as PyTorch, gives direct access to automatic differentiation (AD) for cost-function gradient calculations and execution on GPU and other supported accelerators. We show the library performance in simulations with infinite projected entangled-pair states, such as finding the ground states with AD and simulating thermal states of the Hubbard model via imaginary time evolution. For these challenging examples, we identify and quantify sources of the numerical advantage exploited by the symmetric-tensor implementation.

  4 citations

• Generative learning for the problem of critical slowing down in lattice Gross-Neveu model

  Ankur Singha, Dipankar Chakrabarti, Vipul Arora

  SciPost Phys. Core 5, 052 (2022) · published 2 December 2022 | Toggle abstract · pdf

  In lattice field theory, Monte Carlo simulation algorithms get highly affected by critical slowing down in the critical region, where autocorrelation time increases rapidly. Hence the cost of generation of lattice configurations near the critical region increases sharply. In this paper, we use a Conditional Generative Adversarial Network (C-GAN) for sampling lattice configurations. We train the C-GAN on the dataset consisting of Hybrid Monte Carlo (HMC) samples in regions away from the critical region, i.e., in the regions where the HMC simulation cost is not so high. Then we use the trained C-GAN model to generate independent samples in the critical region. We perform both interpolation and extrapolation to the critical region. Thus, the overall computational cost is reduced. We test our approach for Gross-Neveu model in 1+1 dimension. We find that the observable distributions obtained from the proposed C-GAN model match with those obtained from HMC simulations, while circumventing the problem of critical slowing down.

  4 citations

• A numerical approach for calculating exact non-adiabatic terms in quantum dynamics

  Ewen D. C. Lawrence, Sebastian F. J. Schmid, Ieva Čepaitė, Peter Kirton, Callum W. Duncan

  SciPost Phys. 18, 014 (2025) · published 14 January 2025 | Toggle abstract · pdf

  Understanding how non-adiabatic terms affect quantum dynamics is fundamental to improving various protocols for quantum technologies. We present a novel approach to computing the Adiabatic Gauge Potential (AGP), which gives information on the non-adiabatic terms that arise from time dependence in the Hamiltonian. Our approach uses commutators of the Hamiltonian to build up an appropriate basis of the AGP, which can be easily truncated to give an approximate form when the exact result is intractable. We use this approach to study the AGP obtained for the transverse field Ising model on a variety of graphs, showing how the different underlying graph structures can give rise to very different scaling for the number of terms required in the AGP.

  4 citations