Publications

Order by: publication date number of citations

• Critical energy landscape of linear soft spheres

  Silvio Franz, Antonio Sclocchi, Pierfrancesco Urbani

  SciPost Phys. 9, 012 (2020) · published 23 July 2020 | Toggle abstract · pdf

  We show that soft spheres interacting with a linear ramp potential when overcompressed beyond the jamming point fall in an amorphous solid phase which is critical, mechanically marginally stable and share many features with the jamming point itself. In the whole phase, the relevant local minima of the potential energy landscape display an isostatic contact network of perfectly touching spheres whose statistics is controlled by an infinite lengthscale. Excitations around such energy minima are non-linear, system spanning, and characterized by a set of non-trivial critical exponents. We perform numerical simulations to measure their values and show that, while they coincide, within numerical precision, with the critical exponents appearing at jamming, the nature of the corresponding excitations is richer. Therefore, linear soft spheres appear as a novel class of finite dimensional systems that self-organize into new, critical, marginally stable, states.

  9 citations

• Emergent XY* transition driven by symmetry fractionalization and anyon condensation

  Michael Schuler, Louis-Paul Henry, Yuan-Ming Lu, Andreas M. Läuchli

  SciPost Phys. 14, 001 (2023) · published 11 January 2023 | Toggle abstract · pdf

  Anyons in a topologically ordered phase can carry fractional quantum numbers with respect to the symmetry group of the considered system, one example being the fractional charge of the quasiparticles and quasiholes in the fractional quantum Hall effect. When such symmetry-fractionalized anyons condense, the resulting phase must spontaneously break the symmetry and display a local order parameter. In this paper, we study the phase diagram and anyon condensation transitions of a $\mathbb{Z}_2$ topological order perturbed by Ising interactions in the Toric Code. The interplay between the global (``onsite'') Ising ($\mathbb{Z}_2$) symmetry and the lattice space group symmetries results in a non-trivial symmetry fractionalization class for the anyons, and is shown to lead to two characteristically different confined, symmetry-broken phases. To understand the anyon condensation transitions, we use the recently introduced critical torus energy spectrum technique to identify a line of emergent 2+1D XY* transitions ending at a fine-tuned (Ising$^2$)* critical point. We provide numerical evidence for the occurrence of two symmetry breaking patterns predicted by the specific symmetry fractionalization class of the condensed anyons in the explored phase diagram. In combination with large-scale quantum Monte Carlo simulations we measure unusually large critical exponents $\eta$ for the scaling of the correlation function at the continuous anyon condensation transitions, and we further identify lines of (weakly) first order transitions in the phase diagram. As an important additional result, we discuss the phase diagram of a resulting 2+1D Ashkin-Teller model, where we demonstrate that torus spectroscopy is capable of identifying emergent XY/O(2) critical behaviour, thereby solving some longstanding open questions in the domain of the 3D Ashkin-Teller models. To establish the generality of our results, we propose a field theoretical description capturing the transition from a $\mathbb{Z}_2$ topological order to either $\mathbb{Z}_2$ symmetry broken phase, which is in excellent agreement with the numerical results.

  9 citations

• Ab initio electron-lattice downfolding: Potential energy landscapes, anharmonicity, and molecular dynamics in charge density wave materials

  Arne Schobert, Jan Berges, Erik G. C. P. van Loon, Michael A. Sentef, Sergey Brener, Mariana Rossi, Tim O. Wehling

  SciPost Phys. 16, 046 (2024) · published 9 February 2024 | Toggle abstract · pdf

  The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the simulation of electronic properties and nuclear motion-in particular molecular dynamics (MD) simulations. Three different downfolding strategies based on constraining, unscreening, and combinations thereof are benchmarked against full density functional calculations for selected charge density wave (CDW) systems, namely 1H-TaS$_2$, 1T-TiSe$_2$, 1H-NbS$_2$, and a one-dimensional carbon chain. We find that the downfolded models can reproduce potential energy surfaces on supercells accurately and facilitate computational speedup in MD simulations by about five orders of magnitude in comparison to purely ab initio calculations. For monolayer 1H-TaS$_2$ we report classical and path integral replica exchange MD simulations, revealing the impact of thermal and quantum fluctuations on the CDW transition.

  9 citations

• Inferring effective couplings with restricted Boltzmann machines

  Aurélien Decelle, Cyril Furtlehner, Alfonso De Jesús Navas Gómez, Beatriz Seoane

  SciPost Phys. 16, 095 (2024) · published 5 April 2024 | Toggle abstract · pdf

  Generative models offer a direct way of modeling complex data. Energy-based models attempt to encode the statistical correlations observed in the data at the level of the Boltzmann weight associated with an energy function in the form of a neural network. We address here the challenge of understanding the physical interpretation of such models. In this study, we propose a simple solution by implementing a direct mapping between the Restricted Boltzmann Machine and an effective Ising spin Hamiltonian. This mapping includes interactions of all possible orders, going beyond the conventional pairwise interactions typically considered in the inverse Ising (or Boltzmann Machine) approach, and allowing the description of complex datasets. Earlier works attempted to achieve this goal, but the proposed mappings were inaccurate for inference applications, did not properly treat the complexity of the problem, or did not provide precise prescriptions for practical application. To validate our method, we performed several controlled inverse numerical experiments in which we trained the RBMs using equilibrium samples of predefined models with local external fields, 2-body and 3-body interactions in different sparse topologies. The results demonstrate the effectiveness of our proposed approach in learning the correct interaction network and pave the way for its application in modeling interesting binary variable datasets. We also evaluate the quality of the inferred model based on different training methods.

  9 citations

• Automatic contraction of unstructured tensor networks

  Adam S. Jermyn

  SciPost Phys. 8, 005 (2020) · published 15 January 2020 | Toggle abstract · pdf

  The evaluation of partition functions is a central problem in statistical physics. For lattice systems and other discrete models the partition function may be expressed as the contraction of a tensor network. Unfortunately computing such contractions is difficult, and many methods to make this tractable require periodic or otherwise structured networks. Here I present a new algorithm for contracting unstructured tensor networks. This method makes no assumptions about the structure of the network and performs well in both structured and unstructured cases so long as the correlation structure is local.

  9 citations

• Majorana bound states in encapsulated bilayer graphene

  Fernando Peñaranda, Ramón Aguado, Elsa Prada, Pablo San-Jose

  SciPost Phys. 14, 075 (2023) · published 18 April 2023 | Toggle abstract · pdf

  The search for robust topological superconductivity and Majorana bound states continues, exploring both one-dimensional (1D) systems such as semiconducting nanowires and two-dimensional (2D) platforms. In this work we study a 2D approach based on graphene bilayers encapsulated in transition metal dichalcogenides that, unlike previous proposals involving the Quantum Hall regime in graphene, requires weaker magnetic fields and does not rely on interactions. The encapsulation induces strong spin-orbit coupling on the graphene bilayer, which opens a sizeable gap and stabilizes fragile pairs of helical edge states. We show that, when subject to an in-plane Zeeman field, armchair edges can be transformed into p-wave one-dimensional topological superconductors by contacting them laterally with conventional superconductors. We demonstrate the emergence of Majorana bound states (MBSs) at the sample corners of crystallographically perfect flakes, belonging either to the D or the BDI symmetry classes depending on parameters. We compute the phase diagram, the resilience of MBSs against imperfections, and their manifestation as a 4$\pi$-periodic effect in Josephson junction geometries, all suggesting the existence of a topological phase within experimental reach.

  9 citations

• Collective Monte Carlo updates through tensor network renormalization

  Miguel Frías-Pérez, Michael Mariën, David Pérez García, Mari Carmen Bañuls, Sofyan Iblisdir

  SciPost Phys. 14, 123 (2023) · published 22 May 2023 | Toggle abstract · pdf

  We introduce a Metropolis-Hastings Markov chain for Boltzmann distributions of classical spin systems. It relies on approximate tensor network contractions to propose correlated collective updates at each step of the evolution. We present benchmark computations for a wide variety of instances of the two-dimensional Ising model, including ferromagnetic, antiferromagnetic, (fully) frustrated and Edwards-Anderson spin glass instances, and we show that, with modest computational effort, our Markov chain achieves sizeable acceptance rates, even in the vicinity of critical points. In each of the situations we have considered, the Markov chain compares well with other Monte Carlo schemes such as the Metropolis or Wolff's algorithm: equilibration times appear to be reduced by a factor that varies between $40$ and $2000$, depending on the model and the observable being monitored. We also present an extension to three spatial dimensions, and demonstrate that it exhibits fast equilibration for finite ferro- and antiferromagnetic instances. Additionally, and although it is originally designed for a square lattice of finite degrees of freedom with open boundary conditions, the proposed scheme can be used as such, or with slight modifications, to study triangular lattices, systems with continuous degrees of freedom, matrix models, a confined gas of hard spheres, or to deal with arbitrary boundary conditions.

  8 citations

• Noncoplanar magnetic orders and gapless chiral spin liquid on the kagome lattice with staggered scalar spin chirality

  Fabrizio Oliviero, João Augusto Sobral, Eric C. Andrade, Rodrigo G. Pereira

  SciPost Phys. 13, 050 (2022) · published 5 September 2022 | Toggle abstract · pdf

  Chiral three-spin interactions can suppress long-range magnetic order and stabilize quantum spin liquid states in frustrated lattices. We study a spin-1/2 model on the kagome lattice involving a staggered three-spin interaction $J_\chi$ in addition to Heisenberg exchange couplings $J_1$ on nearest-neighbor bonds and $J_d$ across the diagonals of the hexagons. We explore the phase diagram using a combination of a classical approach, parton mean-field theory, and variational Monte Carlo methods. We obtain a variety of noncoplanar magnetic orders, including a phase that interpolates between cuboc-1 and cuboc-2 states. In the regime of dominant $J_{\chi}$, we find a classically disordered region and argue that it may harbor a gapless chiral spin liquid with a spinon Fermi surface. Our results show that the competition between the staggered three-spin interaction and Heisenberg exchange interactions gives rise to unusual ground states of spin systems.

  8 citations

• Out of equilibrium many-body expansion dynamics of strongly interacting bosons

  Rhombik Roy, Barnali Chakrabarti, Arnaldo Gammal

  SciPost Phys. Core 6, 073 (2023) · published 3 November 2023 | Toggle abstract · pdf

  We solve the Schrödinger equation from first principles to investigate the many-body effects in the expansion dynamics of one-dimensional repulsively interacting bosons released from a harmonic trap. We utilize the multiconfigurational time-dependent Hartree method for bosons (MCTDHB) to solve the many-body Schrödinger equation at high level of accuracy. The MCTDHB basis sets are explicitly time-dependent and optimised by variational principle. We probe the expansion dynamics by three key measures; time evolution of one-, two- and three-body densities. We observe when the mean-field theory results to unimodal expansion, the many-body calculation exhibits trimodal expansion dynamics. The many-body features how the initially fragmented bosons independently spreads out with time whereas the mean-field pictures the expansion of the whole cloud. We also present the three different time scale of dynamics of the inner core, outer core and the cloud as a whole. We analyze the key role played by the dynamical fragmentation during expansion. A Strong evidence of the many-body effects is presented in the dynamics of two- and three-body densities which exhibit correlation hole and pronounced delocalization effect.

  8 citations

• Dynamics of the order parameter statistics in the long range Ising model

  Nishan Ranabhat, Mario Collura

  SciPost Phys. 12, 126 (2022) · published 11 April 2022 | Toggle abstract · pdf

  We study the relaxation of the local ferromagnetic order in the transverse field quantum Ising chain with power-law decaying interactions $1/r^{\alpha}$. We prepare the system in the GHZ state and study the time evolution of the probability distribution function (PDF) of the order parameter within a block of $l$ when quenching the transverse field. The model is known to support long range order at finite temperature for $\alpha \leq 2.0$ . In this regime, quasi-localized topological magnetic defects are expected to strongly affect the equilibration of the full probability distribution. We highlight different dynamical regimes where gaussification mechanism may be slowed down by confinement and eventually breaks. We further study the PDF dynamics induced by changing the effective dimensionality of the system; we mimic this by quenching the range of the interactions. As a matter of fact, the behavior of the system crucially depends on the value of $\alpha$ governing the unitary evolution.

  8 citations