Showing 1–50 of 74 results for author: Ferrari, F

Insulating and metallic phases in the one-dimensional Hubbard-Su-Schrieffer-Heeger model: Insights from a backflow-inspired variational wave function

Authors: Davide Piccioni, Francesco Ferrari, Michele Fabrizio, Federico Becca

Abstract: The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model, by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspir… ▽ More The interplay between electron-electron and electron-phonon interactions is studied in a one-dimensional lattice model, by means of a variational Monte Carlo method based on generalized Jastrow-Slater wave functions. Here, the fermionic part is constructed by a pair-product state, which explicitly depends on the phonon configuration, thus including the electron-phonon coupling in a backflow-inspired way. We report the results for the Hubbard model in presence of the Su-Schrieffer-Heeger coupling to optical phonons, both at half-filling and upon hole doping. At half-filling, the ground state is either a translationally invariant Mott insulator, with gapless spin excitations, or a Peierls insulator, which breaks translations and has fully gapped excitations. Away from half-filling, the charge gap closes in both Mott and Peierls insulators, turning the former into a conventional Luttinger liquid (gapless in all excitation channels). The latter, instead, retains a finite spin gap that closes only above a threshold value of the doping. Even though consistent with the general theory of interacting electrons in one dimension, the existence of such a phase (with gapless charge but gapped spin excitations) has never been demonstrated in a model with repulsive interaction and with only two Fermi points. Since the spin-gapped metal represents the one-dimensional counterpart of a superconductor, our results furnish evidence that a true off-diagonal long-range order may exist in the two-dimensional case. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Landau-Zener without a Qubit: Unveiling Multiphoton Interference, Synthetic Floquet Dimensions, and Dissipative Quantum Chaos

Authors: Leo Peyruchat, Fabrizio Minganti, Marco Scigliuzzo, Filippo Ferrari, Vincent Jouanny, Franco Nori, Vincenzo Savona, Pasquale Scarlino

Abstract: Landau-Zener-Stückelberg-Majorana (LZSM) interference emerges when the parameters of a $\textit{qubit}$ are periodically modulated across an avoided level crossing. Here, we investigate the occurrence of the LZSM phenomenon in nonlinear multilevel bosonic systems, where the interference pattern is determined by multiple energy levels and cannot be described by a level crossing between only two sta… ▽ More Landau-Zener-Stückelberg-Majorana (LZSM) interference emerges when the parameters of a $\textit{qubit}$ are periodically modulated across an avoided level crossing. Here, we investigate the occurrence of the LZSM phenomenon in nonlinear multilevel bosonic systems, where the interference pattern is determined by multiple energy levels and cannot be described by a level crossing between only two states. We fabricate two superconducting resonators made of flux-tunable Josephson junction arrays. The first device is very weakly nonlinear (the nonlinearity is smaller than the photon-loss rate) and, when a weak driving field is applied, it behaves as a linear resonator, yet shows the same LZSM interference as in a two-level system. Notably, here the interference originates from multiple avoided level crossings of the harmonic ladder. When subjected to a stronger drive, nonlinear effects start playing a role, and the interference pattern departs from the one observed in two-level systems. We demonstrate that, when two or more LZSM interference peaks merge, dissipative quantum chaos emerges. In the second device, where the nonlinearity surpasses the photon-loss rate, we observe additional LZSM interference peaks due to Kerr multiphoton resonances. When described under the light of the Floquet theory, these resonances can be interpreted as synthetic modes of an array of coupled cavities. We derive a simple effective model highlighting the essential features of the entirety of these phenomena. As the control of LZSM in qubit systems led to the implementation of fast protocols for characterization and state preparation, our findings pave the way to better control of nonlinear resonators, with implications for diverse quantum technological platforms. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Candidate quantum spin liquids on the maple-leaf lattice

Authors: Jonas Sonnenschein, Atanu Maity, Chunxiao Liu, Ronny Thomale, Francesco Ferrari, Yasir Iqbal

Abstract: Motivated by recent numerical studies reporting putative quantum paramagnetic behavior in spin-$1/2$ Heisenberg models on the maple-leaf lattice, we classify Abrikosov fermion mean-field Ansätze of fully symmetric $U(1)$ and $\mathbb{Z}_{2}$ quantum spin liquids within the framework of projective symmetry groups. We obtain a total of $17$ $U(1)$ and $12$ $\mathbb{Z}_{2}$ algebraic PSGs, and, upon… ▽ More Motivated by recent numerical studies reporting putative quantum paramagnetic behavior in spin-$1/2$ Heisenberg models on the maple-leaf lattice, we classify Abrikosov fermion mean-field Ansätze of fully symmetric $U(1)$ and $\mathbb{Z}_{2}$ quantum spin liquids within the framework of projective symmetry groups. We obtain a total of $17$ $U(1)$ and $12$ $\mathbb{Z}_{2}$ algebraic PSGs, and, upon restricting their realization via mean-field Ansätze with nearest-neighbor amplitudes (relevant to the studied models), only 12 $U(1)$ and 8 $\mathbb{Z}_{2}$ distinct phases are obtained. We present both singlet and triplet fields for all Ansätze up to third nearest-neighbor bonds and discuss their spinon dispersions as well as their dynamical spin structure factors. We further assess the effects of Gutzwiller projection on the equal-time spin structure factors, and identify a $U(1)$ Fermi surface spin liquid whose structure factor most closely reproduces the one obtained from pseudo-fermion functional renormalization group calculations. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 28 pages, 21 figures, 7 tables

Journal ref: Phys. Rev. B 110, 014414 (2024)

A contact map method to capture the features of knot conformations

Authors: Neda Abbasi Taklimi, Franco Ferrari, Marcin Radosław Piątek, Luca Tubiana

Abstract: Inspired by recent advances in the chromosome capture techniques, a method is proposed to study the structural organization of systems of polymers rings with topological constraints.To this purpose, the system is divided into compartments and a simple condition is provided in order to determine if two compartments are in contact or not. Next, a set of contact matrices $\bar T_{ab}$ is defined that… ▽ More Inspired by recent advances in the chromosome capture techniques, a method is proposed to study the structural organization of systems of polymers rings with topological constraints.To this purpose, the system is divided into compartments and a simple condition is provided in order to determine if two compartments are in contact or not. Next, a set of contact matrices $\bar T_{ab}$ is defined that count how many times during a simulation a compartment $a$ was found in contact with a non-contiguous compartment $b$ in conformations with a given energy or temperature. Similar strategies based on correlation maps have been applied to the study of knotted polymers in the recent past. The advantage of the present approach is that is coupled with the Wang-Landau algorithm. Once the density of states is computed, it is possible to generate the contact matrices at any temperature. This gives an immediate overview over the changes of phases that polymer systems undergo. The information on the structure of knotted polymers and links stored in the contact matrices is the result of averaging hundred of billions of conformations and visualized by means of colormaps. The obtained color patterns allow to identify the main properties of the structure of the system under investigation at any temperature. The method is applied to detect the structural rearrangements following the phase transitions of a knotted polymer ring and a circular polycatenane composed by four rings in a solution. It is shown that the colormaps have a finite number of patterns that can be clearly associated with the different phases of these systems. Colormaps also bring new knowledge, for instance predicting the average number of tails appearing in the conformations of the considered polymers at a given temperature. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 27 pages, 11 figures, plain LaTex + RevTex4-1

Field sources in a planar anisotropic CPT-odd gauge model

Authors: L. H. C. Borges, A. F. Ferrari, P. H. O. da Silva, F. A. Barone

Abstract: In the present paper we study some new classical Lorentz violating effects in planar electrodynamics due to the presence of stationary point-like field sources. Starting from the Carroll-Field-Jackiw model defined in (3+1) dimensions, which belongs to the electromagnetic CPT-odd sector of the Standard Model Extension (SME), we apply the dimensional reduction procedure obtaining a (2+1)-dimensional… ▽ More In the present paper we study some new classical Lorentz violating effects in planar electrodynamics due to the presence of stationary point-like field sources. Starting from the Carroll-Field-Jackiw model defined in (3+1) dimensions, which belongs to the electromagnetic CPT-odd sector of the Standard Model Extension (SME), we apply the dimensional reduction procedure obtaining a (2+1)-dimensional model that encompasses an electromagnetic sector composed of the Maxwell-Chern-Simons electrodynamics, a pure scalar sector described by a massless Klein-Gordon field, and a mixed sector where the background vector mediates contributions involving both the scalar and the gauge fields. For all the sectors of this planar theory, we explore some physical phenomena that arise from the interactions between external sources. Specifically, we obtain perturbative effects up to second order in the background vector related to the presence of both electric and scalar planar charges and Dirac points. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 26 pages, 13 figures

Spin-phonon interactions on the kagome lattice: Dirac spin liquid versus valence-bond solids

Authors: Francesco Ferrari, Federico Becca, Roser Valenti

Abstract: We investigate the impact of the spin-phonon coupling on the S=1/2 Heisenberg model on the kagome lattice. For the pure spin model, there is increasing evidence that the low-energy properties can be correctly described by a Dirac spin liquid, in which spinons with a conical dispersion are coupled to emergent gauge fields. Within this scenario, the ground-state wave function is well approximated by… ▽ More We investigate the impact of the spin-phonon coupling on the S=1/2 Heisenberg model on the kagome lattice. For the pure spin model, there is increasing evidence that the low-energy properties can be correctly described by a Dirac spin liquid, in which spinons with a conical dispersion are coupled to emergent gauge fields. Within this scenario, the ground-state wave function is well approximated by a Gutzwiller-projected fermionic state [Y. Ran, M. Hermele, P.A. Lee, and X.-G. Wen, Phys. Rev. Lett. 98, 117205 (2007)]. However, the existence of U(1) gauge fields may naturally lead to instabilities when small perturbations are included. Since phonons are ubiquitous in real materials, they may play a relevant role in the determination of the actual physical properties of the kagome antiferromagnet. We perform a step forward in this direction, including phonon degrees of freedom (at the quantum level) and applying a variational approach based upon Gutzwiller-projected fermionic Ansätze. Our results suggest that the Dirac spin liquid is stable for small spin-phonon couplings, while valence-bond solids are obtained at large couplings. Even though different distortions can be induced by the spin-phonon interaction, the general aspect is that the energy is lowered by maximizing the density of perfect hexagons in the dimerization pattern. △ Less

Submitted 18 April, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Journal ref: Phys. Rev. B 109, 165133 (2024)

$\require{mhchem}$Quantum paramagnetism in the decorated square-kagome antiferromagnet $\ce{Na6Cu7BiO4(PO4)4Cl3}$

Authors: Nils Niggemann, Nikita Astrakhantsev, Arnaud Ralko, Francesco Ferrari, Atanu Maity, Tobias Müller, Johannes Richter, Ronny Thomale, Titus Neupert, Johannes Reuther, Yasir Iqbal, Harald O. Jeschke

Abstract: $\require{mhchem}$The square-kagome lattice Heisenberg antiferromagnet is a highly frustrated Hamiltonian whose material realizations have been scarce. We theoretically investigate the recently synthesized $\ce{Na6Cu7BiO4(PO4)4Cl3}$ where a Cu$^{2+}$ spin-$1/2… ▽ More $\require{mhchem}$The square-kagome lattice Heisenberg antiferromagnet is a highly frustrated Hamiltonian whose material realizations have been scarce. We theoretically investigate the recently synthesized $\ce{Na6Cu7BiO4(PO4)4Cl3}$ where a Cu$^{2+}$ spin-$1/2$ square-kagome lattice (with six site unit cell) is decorated by a seventh magnetic site alternatingly above and below the layers. The material does not show any sign of long-range magnetic order down to 50 mK despite a Curie-Weiss temperature of $-212$ K indicating a quantum paramagnetic phase. Our DFT energy mapping elicits a purely antiferromagnetic Hamiltonian that features longer range exchange interactions beyond the pure square-kagome model and, importantly, we find the seventh site to be strongly coupled to the plane. We combine two variational Monte Carlo approaches, pseudo-fermion/Majorana functional renormalization group and Schwinger-Boson mean field calculations to show that the complex Hamiltonian of $\ce{Na6Cu7BiO4(PO4)4Cl3}$ still features a nonmagnetic ground state. We explain how the seventh Cu$^{2+}$ site actually aids the stabilization of the disordered state. We predict static and dynamic spin structure factors to guide future neutron scattering experiments. △ Less

Submitted 20 December, 2023; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: 5 pages, 4 figures. Supplement: 7 pages, 11 figures v2: minor revisions in text and figure labels

Journal ref: Phys. Rev. B 108, L241117 (2023)

Quantum spin liquids on the diamond lattice

Authors: Aishwarya Chauhan, Atanu Maity, Chunxiao Liu, Jonas Sonnenschein, Francesco Ferrari, Yasir Iqbal

Abstract: We perform a projective symmetry group classification of spin $S=1/2$ symmetric quantum spin liquids with different gauge groups on the diamond lattice. Employing the Abrikosov fermion representation, we obtain $8$ $SU(2)$, $62$ $U(1)$ and $80$ $\mathbb{Z}_{2}$ algebraic PSGs. Constraining these solutions to mean-field parton Ansätze with short-range amplitudes, the classification reduces to only… ▽ More We perform a projective symmetry group classification of spin $S=1/2$ symmetric quantum spin liquids with different gauge groups on the diamond lattice. Employing the Abrikosov fermion representation, we obtain $8$ $SU(2)$, $62$ $U(1)$ and $80$ $\mathbb{Z}_{2}$ algebraic PSGs. Constraining these solutions to mean-field parton Ansätze with short-range amplitudes, the classification reduces to only $2$ $SU(2)$, $7$ $U(1)$ and $8$ $\mathbb{Z}_{2}$ distinctly realizable phases. We obtain both the singlet and triplet fields for all Ansätze, discuss the spinon dispersions, and present the dynamical spin structure factors within a self-consistent treatment of the Heisenberg Hamiltonian with up to third-nearest neighbor couplings. Interestingly, we find that a zero-flux $SU(2)$ state and some descendent $U(1)$ and $\mathbb{Z}_{2}$ states host robust gapless nodal loops in their dispersion spectrum, owing their stability at the mean-field level to the projective implementation of rotoinversion and screw symmetries. A nontrivial connection is drawn between one of our $U(1)$ spinon Hamiltonians (belonging to the nonprojective class) and the Fu-Kane-Mele model for a three-dimensional topological insulator on the diamond lattice. We show that Gutzwiller projection of the 0- and $π$-flux $SU(2)$ spin liquids generates long-range Néel order. △ Less

Submitted 19 October, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: Editors' Suggestion. 36 pages, 9 figures, 9 tables

Journal ref: Phys. Rev. B 108, 134424 (2023)

Steady-state quantum chaos in open quantum systems

Authors: Filippo Ferrari, Luca Gravina, Debbie Eeltink, Pasquale Scarlino, Vincenzo Savona, Fabrizio Minganti

Abstract: We introduce the notion of steady-state quantum chaos as a general phenomenon in open quantum many-body systems. Classifying an isolated or open quantum system as integrable or chaotic relies in general on the properties of the equations governing its time evolution. This however may fail in predicting the actual nature of the quantum dynamics, that can be either regular or chaotic depending on th… ▽ More We introduce the notion of steady-state quantum chaos as a general phenomenon in open quantum many-body systems. Classifying an isolated or open quantum system as integrable or chaotic relies in general on the properties of the equations governing its time evolution. This however may fail in predicting the actual nature of the quantum dynamics, that can be either regular or chaotic depending on the initial state. Chaos and integrability in the steady state of an open quantum system are instead uniquely determined by the spectral structure of the time evolution generator. To characterize steady-state quantum chaos we introduce a spectral analysis based on the spectral statistics of quantum trajectories (SSQT). We test the generality and reliability of the SSQT criterion on several dissipative systems, further showing that an open system with chaotic structure can evolve towards either a chaotic or integrable steady state. We study steady-state chaos in the driven-dissipative Bose-Hubbard model, a paradigmatic example of out-of-equilibrium bosonic system without particle number conservation. This system is widely employed as a building block in state-of-the-art noisy intermediate-scale quantum devices, with applications in quantum computation and sensing. Finally, our analysis shows the existence of an emergent dissipative quantum chaos, where the classical and semi-classical limits display an integrable behaviour. This emergent dissipative quantum chaos arises from the quantum and classical fluctuations associated with the dissipation mechanism. Our work establishes a fundamental understanding of the integrable and chaotic dynamics of open quantum systems and paves the way for the investigation of dissipative quantum chaos and its consequences on quantum technologies. △ Less

Submitted 29 November, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: 23 pages, 12 figures

Variational Benchmarks for Quantum Many-Body Problems

Authors: Dian Wu, Riccardo Rossi, Filippo Vicentini, Nikita Astrakhantsev, Federico Becca, Xiaodong Cao, Juan Carrasquilla, Francesco Ferrari, Antoine Georges, Mohamed Hibat-Allah, Masatoshi Imada, Andreas M. Läuchli, Guglielmo Mazzola, Antonio Mezzacapo, Andrew Millis, Javier Robledo Moreno, Titus Neupert, Yusuke Nomura, Jannes Nys, Olivier Parcollet, Rico Pohle, Imelda Romero, Michael Schmid, J. Maxwell Silvester, Sandro Sorella , et al. (8 additional authors not shown)

Abstract: The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, ide… ▽ More The continued development of novel many-body approaches to ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. Here we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide the most extensive curated dataset of variational calculations of many-body quantum systems to date, identifying cases where state-of-the-art numerical approaches show limited accuracy, and novel algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods towards quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 25 pages, 5 figures

A $j_{\rm eff}=1/2$ Kitaev material on the triangular lattice: The case of NaRuO$_2$

Authors: Aleksandar Razpopov, David A. S. Kaib, Steffen Backes, Leon Balents, Stephen D. Wilson, Francesco Ferrari, Kira Riedl, Roser Valenti

Abstract: Motivated by recent reports of a quantum disordered ground state in the triangular lattice compound NaRuO$_2$, we derive a $j_{\rm eff}=1/2$ magnetic model for this system by means of first-principles calculations. The pseudospin Hamiltonian is dominated by bond-dependent off-diagonal $Γ$ interactions, complemented by a ferromagnetic Heisenberg exchange and a notably antiferromagnetic Kitaev term.… ▽ More Motivated by recent reports of a quantum disordered ground state in the triangular lattice compound NaRuO$_2$, we derive a $j_{\rm eff}=1/2$ magnetic model for this system by means of first-principles calculations. The pseudospin Hamiltonian is dominated by bond-dependent off-diagonal $Γ$ interactions, complemented by a ferromagnetic Heisenberg exchange and a notably antiferromagnetic Kitaev term. In addition to bilinear interactions, we find a sizable four-spin ring exchange contribution with a strongly anisotropic character, which has been so far overlooked when modeling Kitaev materials. The analysis of the magnetic model, based on the minimization of the classical energy and exact diagonalization of the quantum Hamiltonian, points toward the existence of a rather robust easy-plane ferromagnetic order, which cannot be easily destabilized by physically relevant perturbations. △ Less

Submitted 17 July, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: v2 includes supplementary information, a more detailed discussion of selected aspects and correction of small errors

Journal ref: npj Quantum Mater. 8, 36 (2023)

Projective symmetry group classification of Abrikosov fermion mean-field ansätze on the square-octagon lattice

Authors: Atanu Maity, Francesco Ferrari, Ronny Thomale, Saptarshi Mandal, Yasir Iqbal

Abstract: We perform a projective symmetry group (PSG) classification of symmetric quantum spin liquids with different gauge groups on the square-octagon lattice. Employing the Abrikosov fermion representation for spin-$1/2$, we obtain $32$ $SU(2)$, $1808$ $U(1)$ and $384$ $\mathbb{Z}_{2}$ algebraic PSGs. Constraining ourselves to mean-field parton ansätze with short-range amplitudes, the classification red… ▽ More We perform a projective symmetry group (PSG) classification of symmetric quantum spin liquids with different gauge groups on the square-octagon lattice. Employing the Abrikosov fermion representation for spin-$1/2$, we obtain $32$ $SU(2)$, $1808$ $U(1)$ and $384$ $\mathbb{Z}_{2}$ algebraic PSGs. Constraining ourselves to mean-field parton ansätze with short-range amplitudes, the classification reduces to a limited number, with 4 $SU(2)$, 24 $U(1)$ and 36 $\mathbb{Z}_{2}$, distinct phases. We discuss their ground state properties and spinon dispersions within a self-consistent treatment of the Heisenberg Hamiltonian with frustrating couplings. △ Less

Submitted 27 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 28 pages, 8 figures, 5 tables

Journal ref: Phys. Rev. B 107, 134438 (2023)

Detecting topological phases in the square-octagon lattice with statistical methods

Authors: Paul Wunderlich, Francesco Ferrari, Roser Valentí

Abstract: Electronic systems living on Archimedean lattices such as kagome and square-octagon networks are presently being intensively discussed for the possible realization of topological insulating phases. Coining the most interesting electronic topological states in an unbiased way is however not straightforward due to the large parameter space of possible Hamiltonians. A possible approach to tackle this… ▽ More Electronic systems living on Archimedean lattices such as kagome and square-octagon networks are presently being intensively discussed for the possible realization of topological insulating phases. Coining the most interesting electronic topological states in an unbiased way is however not straightforward due to the large parameter space of possible Hamiltonians. A possible approach to tackle this problem is provided by a recently developed statistical learning method [T. Mertz and R. Valentí, Phys. Rev. Research 3, 013132 (2021)], based on the analysis of a large data sets of randomized tight-binding Hamiltonians labeled with a topological index. In this work, we complement this technique by introducing a feature engineering approach which helps identifying polynomial combinations of Hamiltonian parameters that are associated with non-trivial topological states. As a showcase, we employ this method to investigate the possible topological phases that can manifest on the square-octagon lattice, focusing on the case in which the Fermi level of the system lies at a high-order van Hove singularity, in analogy to recent studies of topological phases on the kagome lattice at the van Hove filling. △ Less

Submitted 4 May, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 9 pages, 6 figures

Journal ref: EPJ Plus 138, 336 (2023)

On the thermal properties of knotted block copolymer rings

Authors: Neda Abbasi Taklimi, Franco Ferrari, Marcin Radosław Piątek, Luca Tubiana

Abstract: The thermal properties of coarse grained knotted polymers containing two kinds of monomers $A$ and $B$ fluctuating in a solution are investigated on a simple cubic lattice using the Wang-Landau MC algorithm. These knots have a more complex phase diagram than knots formed by homopolymers, including the possible presence of metastable states. Two different setups are considered: i) charged block cop… ▽ More The thermal properties of coarse grained knotted polymers containing two kinds of monomers $A$ and $B$ fluctuating in a solution are investigated on a simple cubic lattice using the Wang-Landau MC algorithm. These knots have a more complex phase diagram than knots formed by homopolymers, including the possible presence of metastable states. Two different setups are considered: i) charged block copolymers in a ion solution and ii) neutral copolymers with the $A$ monomers above and the $B$ monomers below the theta point. A precise interpretation of the peaks observed in the plots of the specific heat capacity is provided. In view of possible applications in medicine and the construction of intelligent materials, it is also shown that the behavior of copolymer rings can be tuned by changing both their monomer configuration and topology. We find that the most stable compact states are formed by charged copolymers in which very short segments with $A$ monomers are alternated by short segments with $B$ monomers. In such knots the transition from the compact to the expanded state is very fast, leading to a narrow and high peak in the specific heat capacity which appears at very high temperatures. The effects of topology allow to tune the radius of gyration of the knotted polymer ring and to increase or decrease the temperatures at which the observed phase transitions or rearrangements of the system occur. While we observe a general fading out of the influence of topology in longer polymers, our simulations have captured a few exceptions to this rule. △ Less

Submitted 5 September, 2022; originally announced September 2022.

Comments: 35 pages, 19 figures

Static and dynamical signatures of Dzyaloshinskii-Moriya interactions in the Heisenberg model on the kagome lattice

Authors: Francesco Ferrari, Sen Niu, Juraj Hasik, Yasir Iqbal, Didier Poilblanc, Federico Becca

Abstract: Motivated by recent experiments on Cs$_2$Cu$_3$SnF$_{12}$ and YCu$_{3}$(OH)$_{6}$Cl$_{3}$, we consider the ${S=1/2}$ Heisenberg model on the kagome lattice with nearest-neighbor super-exchange $J$ and (out-of-plane) Dzyaloshinskii-Moriya interaction $J_D$, which favors (in-plane) ${\bf Q}=(0,0)$ magnetic order. By using both variational Monte Carlo (based upon Gutzwiller-projected fermionic wave f… ▽ More Motivated by recent experiments on Cs$_2$Cu$_3$SnF$_{12}$ and YCu$_{3}$(OH)$_{6}$Cl$_{3}$, we consider the ${S=1/2}$ Heisenberg model on the kagome lattice with nearest-neighbor super-exchange $J$ and (out-of-plane) Dzyaloshinskii-Moriya interaction $J_D$, which favors (in-plane) ${\bf Q}=(0,0)$ magnetic order. By using both variational Monte Carlo (based upon Gutzwiller-projected fermionic wave functions) and tensor-network approaches (built from infinite projected-entangled pair/simplex states), we show that the ground state develops a finite magnetization for $J_D/J \gtrsim 0.03 - 0.04$, while the gapless spin liquid remains stable for smaller values of the Dzyaloshinskii-Moriya interaction. The relatively small value of $J_D/J$ for which magnetic order sets in is particularly relevant for the interpretation of low-temperature behaviors of kagome antiferromagnets, including ZnCu$_{3}$(OH)$_{6}$Cl$_{2}$. In addition, we assess the spin dynamical structure factors and the corresponding low-energy spectrum, by using the variational Monte Carlo technique. The existence of a continuum of excitations above the magnon modes is reported within the magnetically ordered phase, similarly to what has been detected by inelastic neutron scattering on Cs$_{2}$Cu$_{3}$SnF$_{12}$. △ Less

Submitted 20 February, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 15 figures

Journal ref: SciPost Phys. 14, 139 (2023)

Pinch-points to half-moons and up in the stars: the kagome skymap

Authors: Dominik Kiese, Francesco Ferrari, Nikita Astrakhantsev, Nils Niggemann, Pratyay Ghosh, Tobias Müller, Ronny Thomale, Titus Neupert, Johannes Reuther, Michel J. P. Gingras, Simon Trebst, Yasir Iqbal

Abstract: Pinch point singularities, associated with flat band magnetic excitations, are tell-tale signatures of Coulomb spin liquids. While their properties in the presence of quantum fluctuations have been widely studied, the fate of the complementary non-analytic features -- shaped as half-moons and stars -- arising from adjacent shallow dispersive bands has remained unexplored. Here, we address this que… ▽ More Pinch point singularities, associated with flat band magnetic excitations, are tell-tale signatures of Coulomb spin liquids. While their properties in the presence of quantum fluctuations have been widely studied, the fate of the complementary non-analytic features -- shaped as half-moons and stars -- arising from adjacent shallow dispersive bands has remained unexplored. Here, we address this question for the spin $S=1/2$ Heisenberg antiferromagnet on the kagome lattice with second and third neighbor couplings, which allows one to tune the classical ground state from flat bands to being governed by shallow dispersive bands for intermediate coupling strengths. Employing the complementary strengths of variational Monte Carlo, pseudo-fermion functional renormalization group, and density-matrix renormalization group, we establish the quantum phase diagram. The U(1) Dirac spin liquid ground state of the nearest-neighbor antiferromagnet remains remarkably robust till intermediate coupling strengths when it transitions into a pinwheel valence bond crystal displaying signatures of half-moons in its structure factor. Our work thus identifies a microscopic setting that realizes one of the proximate orders of the Dirac spin liquid identified in a recent work [Song, Wang, Vishwanath, He, Nat. Commun. 10, 4254 (2019)]. For larger couplings, we obtain a collinear magnetically ordered ground state characterized by star-like patterns. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 7 pages, 5 figures. Supplemental Material appended

Journal ref: Phys. Rev. Research 5, L012025 (2023)

Charge-density waves in kagome-lattice extended Hubbard models at the van Hove filling

Authors: Francesco Ferrari, Federico Becca, Roser Valentí

Abstract: The Hubbard model on the kagome lattice is presently often considered as a minimal model to describe the rich low-temperature behavior of AV$_{3}$Sb$_{5}$ compounds (with A=K, Rb, Cs), including charge-density waves (CDWs), superconductivity, and possibly broken time-reversal symmetry. Here, we investigate, via variational Jastrow-Slater wave functions, the properties of its ground state when both… ▽ More The Hubbard model on the kagome lattice is presently often considered as a minimal model to describe the rich low-temperature behavior of AV$_{3}$Sb$_{5}$ compounds (with A=K, Rb, Cs), including charge-density waves (CDWs), superconductivity, and possibly broken time-reversal symmetry. Here, we investigate, via variational Jastrow-Slater wave functions, the properties of its ground state when both onsite $U$ and nearest-neighbor $V$ Coulomb repulsions are considered at the van Hove filling. Our calculations reveal the presence of different interaction-driven CDWs and, contrary to previous renormalization-group studies, the absence of ferromagnetism and charge- or spin-bond order. No signatures of chiral phases are detected. Remarkably, the CDWs triggered by the nearest-neighbor repulsion possess charge disproportionations that are not compatible with the ones observed in AV$_{3}$Sb$_{5}$. As an alternative mechanism to stabilize charge-bond order, we consider the electron-phonon interaction, modeled by coupling the hopping amplitudes to quantum phonons, as in the Su-Schrieffer-Heeger model. Our results show the instability towards a tri-hexagonal distortion with $2\times 2$ periodicity, in a closer agreement with experimental findings. △ Less

Submitted 15 August, 2022; v1 submitted 21 April, 2022; originally announced April 2022.

Journal ref: Phys. Rev. B 106, L081107 (2022)

Accuracy of Restricted Boltzmann Machines for the one-dimensional $J_1-J_2$ Heisenberg model

Authors: Luciano Loris Viteritti, Francesco Ferrari, Federico Becca

Abstract: Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 $J_1-J_2$ antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has… ▽ More Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 $J_1-J_2$ antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has a non-trivial sign structure, especially for $J_2/J_1>0.5$, forcing us to work with complex-valued RBMs. Two variational Ansätze are discussed: one defined through a fully complex RBM, and one in which two different real-valued networks are used to approximate modulus and phase of the wave function. In both cases, translational invariance is imposed by considering linear combinations of RBMs, giving access also to the lowest-energy excitations at fixed momentum $k$. We perform a systematic study on small clusters to evaluate the accuracy of these wave functions in comparison to exact results, providing evidence for the supremacy of the fully complex RBM. Our calculations show that this kind of Ansätze is very flexible and describes both gapless and gapped ground states, also capturing the incommensurate spin-spin correlations and low-energy spectrum for $J_2/J_1>0.5$. The RBM results are also compared to the ones obtained with Gutzwiller-projected fermionic states, often employed to describe quantum spin models [F. Ferrari, A. Parola, S. Sorella and F. Becca, Phys. Rev. B 97, 235103 (2018)]. Contrary to the latter class of variational states, the fully-connected structure of RBMs hampers the transferability of the wave function from small to large clusters, implying an increase of the computational cost with the system size. △ Less

Submitted 14 April, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 21 pages, 12 figures, submitted to Scipost Physics

Journal ref: SciPost Phys. 12, 166 (2022)

Ring-o-rings: a new category of supramolecular structures with topologically tunable properties

Authors: Luca Tubiana, Franco Ferrari, Enzo Orlandini

Abstract: Macrochains of topologically interlocked rings with unique physical properties have recently gained considerable interest in supramolecular chemistry, biology, and soft matter. Most of the work has been, so far, focused on linear chains and on their variety of conformational properties compared to standard polymers. Here we go beyond the linear case and show that, by circularizing such macrochains… ▽ More Macrochains of topologically interlocked rings with unique physical properties have recently gained considerable interest in supramolecular chemistry, biology, and soft matter. Most of the work has been, so far, focused on linear chains and on their variety of conformational properties compared to standard polymers. Here we go beyond the linear case and show that, by circularizing such macrochains, one can exploit the topology of the local interlockings to store torsional stress in the system, altering significantly its metric and local properties. Moreover, by properly defining the twist (Tw) and writhe (Wr) of these macrorings we show the validity of a relation equivalent to the Cǎlugǎreanu-White-Fuller theorem $Tw + Wr$=const, originally proved for ribbon like structures such as ds-DNA. Our results suggest that circular structures of topologically linked rings with storable and tunable torsion can form a new category of highly designable multiscale structures with potential applications in supramolecular chemistry and material science △ Less

Submitted 16 December, 2021; originally announced December 2021.

Comments: 10 pages, 8 figures

Statistical learning of engineered topological phases in the kagome superlattice of AV$_3$Sb$_5$

Authors: Thomas Mertz, Paul Wunderlich, Shinibali Bhattacharyya, Francesco Ferrari, Roser Valentí

Abstract: Recent experimental findings have reported the presence of unconventional charge orders in the enlarged ($2 \times 2$) unit-cell of kagome metals AV$_3$Sb$_5$ (A=K,Rb,Cs) and hinted towards specific topological signatures. Motivated by these discoveries, we investigate the types of topological phases that can be realized in such kagome superlattices. In this context, we employ a recently introduce… ▽ More Recent experimental findings have reported the presence of unconventional charge orders in the enlarged ($2 \times 2$) unit-cell of kagome metals AV$_3$Sb$_5$ (A=K,Rb,Cs) and hinted towards specific topological signatures. Motivated by these discoveries, we investigate the types of topological phases that can be realized in such kagome superlattices. In this context, we employ a recently introduced statistical method capable of constructing topological models for any generic lattice. By analyzing large data sets generated from symmetry-guided distributions of randomized tight-binding parameters, and labeled with the corresponding topological index, we extract physically meaningful information. We illustrate the possible real-space manifestations of charge and bond modulations and associated flux patterns for different topological classes, and discuss their relation to present theoretical predictions and experimental signatures for the AV$_3$Sb$_5$ family. Simultaneously, we predict new higher-order topological phases that may be realized by appropriately manipulating the currently known systems. △ Less

Submitted 31 January, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Journal ref: npj Comput Mater 8, 66 (2022)

Pinwheel valence-bond-crystal ground state of the spin-$\frac{1}{2}$ Heisenberg antiferromagnet on the $shuriken$ lattice

Authors: Nikita Astrakhantsev, Francesco Ferrari, Nils Niggemann, Tobias Müller, Aishwarya Chauhan, Augustine Kshetrimayum, Pratyay Ghosh, Nicolas Regnault, Ronny Thomale, Johannes Reuther, Titus Neupert, Yasir Iqbal

Abstract: We investigate the nature of the ground-state of the spin-$\frac{1}{2}$ Heisenberg antiferromagnet on the $shuriken$ lattice by complementary state-of-the-art numerical techniques, such as variational Monte Carlo (VMC) with versatile Gutzwiller-projected Jastrow wave functions, unconstrained multi-variable variational Monte Carlo (mVMC), and pseudo-fermion/Majorana functional renormalization group… ▽ More We investigate the nature of the ground-state of the spin-$\frac{1}{2}$ Heisenberg antiferromagnet on the $shuriken$ lattice by complementary state-of-the-art numerical techniques, such as variational Monte Carlo (VMC) with versatile Gutzwiller-projected Jastrow wave functions, unconstrained multi-variable variational Monte Carlo (mVMC), and pseudo-fermion/Majorana functional renormalization group (PF/PM-FRG) methods. We establish the presence of a quantum paramagnetic ground state and investigate its nature, by classifying symmetric and chiral quantum spin liquids, and inspecting their instabilities towards competing valence-bond-crystal (VBC) orders. Our VMC analysis reveals that a VBC with a pinwheel structure emerges as the lowest-energy variational ground state, and it is obtained as an instability of the U(1) Dirac spin liquid. Analogous conclusions are drawn from mVMC calculations employing accurate BCS pairing states supplemented by symmetry projectors, which confirm the presence of pinwheel VBC order by a thorough analysis of dimer-dimer correlation functions. Our work highlights the nontrivial role of quantum fluctuations via the Gutzwiller projector in resolving the subtle interplay between competing orders. △ Less

Submitted 17 November, 2021; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: Main paper (7 pages, 5 figures, 1 table) + Supplemental Material

Journal ref: Phys. Rev. B 104, L220408 (2021)

Dimerization tendencies of the pyrochlore Heisenberg antiferromagnet: A functional renormalization group perspective

Authors: Max Hering, Vincent Noculak, Francesco Ferrari, Yasir Iqbal, Johannes Reuther

Abstract: We investigate the ground state properties of the spin-$1/2$ pyrochlore Heisenberg antiferromagnet using pseudofermion functional renormalization group techniques. The first part of our analysis is based on an enhanced parton mean-field approach which takes into account fluctuation effects from renormalized vertex functions. Our implementation of this technique extends earlier approaches and resol… ▽ More We investigate the ground state properties of the spin-$1/2$ pyrochlore Heisenberg antiferromagnet using pseudofermion functional renormalization group techniques. The first part of our analysis is based on an enhanced parton mean-field approach which takes into account fluctuation effects from renormalized vertex functions. Our implementation of this technique extends earlier approaches and resolves technical difficulties associated with a diagrammatic overcounting. Using various parton ansätze for quantum spin liquids, dimerized and nematic states our results indicate a tendency for lattice symmetry breaking in the ground state. While overall quantum spin liquids seem unfavorable in this system, the recently proposed monopole state still shows the strongest support among all spin liquid ansätze that we have tested, which is further confirmed by our complementary variational Monte Carlo calculations. In the second part of our investigation, we probe lattice symmetry breaking more directly by applying the pseudofermion functional renormalization group to perturbed systems. Our results from this technique confirm that the system's ground state either exhibits broken $C_3$ rotation symmetry, or a combination of inversion and $C_3$ symmetry breaking. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 14 pages, 5 figures

Journal ref: Phys. Rev. B 105, 054426 (2022)

Gutzwiller-projected states for the $J_1$-$J_2$ Heisenberg model on the kagome lattice: achievements and pitfalls

Authors: Yasir Iqbal, Francesco Ferrari, Aishwarya Chauhan, Alberto Parola, Didier Poilblanc, Federico Becca

Abstract: We assess the ground-state phase diagram of the $J_1$-$J_2$ Heisenberg model on the kagome lattice by employing Gutzwiller-projected fermionic wave functions. Within this framework, different states can be represented, defined by distinct unprojected fermionic Hamiltonians that include hopping and pairing terms, as well as a coupling to local Zeeman fields to generate magnetic order. For $J_2=0$,… ▽ More We assess the ground-state phase diagram of the $J_1$-$J_2$ Heisenberg model on the kagome lattice by employing Gutzwiller-projected fermionic wave functions. Within this framework, different states can be represented, defined by distinct unprojected fermionic Hamiltonians that include hopping and pairing terms, as well as a coupling to local Zeeman fields to generate magnetic order. For $J_2=0$, the so-called U(1) Dirac state, in which only hopping is present (such as to generate a $π$-flux in the hexagons), has been shown to accurately describe the exact ground state [Y. Iqbal, F. Becca, S. Sorella, and D. Poilblanc, Phys. Rev. B 87, 060405 (2013); Y.-C. He, M. P. Zaletel, M. Oshikawa, and F. Pollmann, Phys. Rev. X 7, 031020 (2017)]. Here, we show that its accuracy improves in presence of a small antiferromagnetic super-exchange $J_2$, leading to a finite region where the gapless spin liquid is stable; then, for $J_2/J_1=0.11(1)$, a first-order transition to a magnetic phase with pitch vector ${\bf q}=(0,0)$ is detected, by allowing magnetic order within the fermionic Hamiltonian. Instead, for small ferromagnetic values of $|J_2|/J_1$, the situation is more contradictory. While the U(1) Dirac state remains stable against several perturbations in the fermionic part (i.e., dimerization patterns or chiral terms), its accuracy clearly deteriorates on small systems, most notably on $36$ sites where exact diagonalization is possible. Then, upon increasing the ratio $|J_2|/J_1$, a magnetically ordered state with $\sqrt{3} \times \sqrt{3}$ periodicity eventually overcomes the U(1) Dirac spin liquid. Within the ferromagnetic $J_{2}$ regime, evidence is shown in favor of a first-order transition at $J_2/J_1=-0.065(5)$. △ Less

Submitted 19 September, 2021; v1 submitted 5 August, 2021; originally announced August 2021.

Comments: 9 pages, 6 figures, 1 table

Journal ref: Phys. Rev. B 104, 144406 (2021)

Phase diagram of a distorted kagome antiferromagnet and application to Y-kapellasite

Authors: Max Hering, Francesco Ferrari, Aleksandar Razpopov, Igor I. Mazin, Roser Valenti, Harald O. Jeschke, Johannes Reuther

Abstract: We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings Jhexagon, J and J', and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear Q = 0 magnetic phase, two unusual non-colli… ▽ More We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings Jhexagon, J and J', and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear Q = 0 magnetic phase, two unusual non-collinear coplanar Q = (1/3,1/3) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y3Cu9(OH)19Cl8 is a realization of this model and predict its ground state to lie in the region of Q = (1/3,1/3) order, which remains stable even after inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets. △ Less

Submitted 29 November, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: 13 pages, 8 figures, 6 pages of supplement

Journal ref: npj Comput. Mater. 8, 10 (2022)

Effects of spin-phonon coupling in frustrated Heisenberg models

Authors: Francesco Ferrari, Roser Valenti, Federico Becca

Abstract: The existence and stability of spin-liquid phases represent a central topic in the field of frustrated magnetism. While a few examples of spin-liquid ground states are well established in specific models (e.g. the Kitaev model on the honeycomb lattice), recent investigations have suggested the possibility of their appearance in several Heisenberg-like models on frustrated lattices. An important re… ▽ More The existence and stability of spin-liquid phases represent a central topic in the field of frustrated magnetism. While a few examples of spin-liquid ground states are well established in specific models (e.g. the Kitaev model on the honeycomb lattice), recent investigations have suggested the possibility of their appearance in several Heisenberg-like models on frustrated lattices. An important related question concerns the stability of spin liquids in the presence of small perturbations in the Hamiltonian. In this respect, the magnetoelastic interaction between spins and phonons represents a relevant and physically motivated perturbation, which has been scarcely investigated so far. In this work, we study the effect of the spin-phonon coupling on prototypical models of frustrated magnetism. We adopt a variational framework based upon Gutzwiller-projected wave functions implemented with a spin-phonon Jastrow factor, providing a full quantum treatment of both spin and phonon degrees of freedom. The results on the frustrated $J_1-J_2$ Heisenberg model on one- and two-dimensional (square) lattices show that, while a valence-bond crystal is prone to lattice distortions, a gapless spin liquid is stable for small spin-phonon couplings. In view of the ubiquitous presence of lattice vibrations, our results are particularly important to demonstrate the possibility that gapless spin liquids may be realized in real materials. △ Less

Submitted 14 July, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: 10 pages, 10 figures

Journal ref: Phys. Rev. B 104, 035126 (2021)

Gapless spin liquids in disguise

Authors: Francesco Ferrari, Alberto Parola, Federico Becca

Abstract: We show that gapless spin liquids, which are potential candidates to describe the ground state of frustrated Heisenberg models in two dimensions, become trivial insulators on cylindrical geometries with an even number of legs. In particular, we report calculations for Gutzwiller-projected fermionic states on strips of square and kagome lattices. By choosing different boundary conditions for the fe… ▽ More We show that gapless spin liquids, which are potential candidates to describe the ground state of frustrated Heisenberg models in two dimensions, become trivial insulators on cylindrical geometries with an even number of legs. In particular, we report calculations for Gutzwiller-projected fermionic states on strips of square and kagome lattices. By choosing different boundary conditions for the fermionic degrees of freedom, both gapless and gapped states may be realized, the latter ones having a lower variational energy. The direct evaluation of static and dynamical correlation functions, as well as overlaps between different states, allows us to demonstrate the sharp difference between the ground-state properties obtained within cylinders or directly in the two-dimensional lattice. Our results shed light on the difficulty to detect bona fide gapless spin liquids in such cylindrical geometries. △ Less

Submitted 18 May, 2021; v1 submitted 9 March, 2021; originally announced March 2021.

Comments: 9 pages, 10 figures

Journal ref: Phys. Rev. B 103, 195140 (2021)

Towards a Topological Quantum Chemistry description of correlated systems: the case of the Hubbard diamond chain

Authors: Mikel Iraola, Niclas Heinsdorf, Apoorv Tiwari, Dominik Lessnich, Thomas Mertz, Francesco Ferrari, Mark H. Fischer, Stephen M. Winter, Frank Pollmann, Titus Neupert, Roser Valentí, Maia G. Vergniory

Abstract: The recently introduced topological quantum chemistry (TQC) framework has provided a description of universal topological properties of all possible band insulators in all space groups based on crystalline unitary symmetries and time reversal. While this formalism filled the gap between the mathematical classification and the practical diagnosis of topological materials, an obvious limitation is t… ▽ More The recently introduced topological quantum chemistry (TQC) framework has provided a description of universal topological properties of all possible band insulators in all space groups based on crystalline unitary symmetries and time reversal. While this formalism filled the gap between the mathematical classification and the practical diagnosis of topological materials, an obvious limitation is that it only applies to weakly interacting systems-which can be described within band theory. It is an open question to which extent this formalism can be generalized to correlated systems that can exhibit symmetry protected topological phases which are not adiabatically connected to any band insulator. In this work we address the many facettes of this question by considering the specific example of a Hubbard diamond chain. This model features a Mott insulator, a trivial insulating phase and an obstructed atomic limit phase. Here we discuss the nature of the Mott insulator and determine the phase diagram and topology of the interacting model with infinite density matrix renormalization group calculations, variational Monte Carlo simulations and with many-body topological invariants. We then proceed by considering a generalization of the TQC formalism to Green's functions combined with the concept of topological Hamiltonian to identify the topological nature of the phases, using cluster perturbation theory to calculate the Green's functions. The results are benchmarked with the above determined phase diagram and we discuss the applicability and limitations of the approach and its possible extensions. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 16 pages, 13 figures

Variational wave functions for the spin-Peierls transition in the Su-Schrieffer-Heeger model with quantum phonons

Authors: Francesco Ferrari, Roser Valenti, Federico Becca

Abstract: We introduce variational wave functions to evaluate the ground-state properties of spin-phonon coupled systems described by the Su-Schrieffer-Heeger model. Quantum spins and phonons are treated on equal footing within a Monte Carlo sampling, and different regimes are investigated. We show that the proposed variational Ansatz yields good agreement with previous density-matrix renormalization group… ▽ More We introduce variational wave functions to evaluate the ground-state properties of spin-phonon coupled systems described by the Su-Schrieffer-Heeger model. Quantum spins and phonons are treated on equal footing within a Monte Carlo sampling, and different regimes are investigated. We show that the proposed variational Ansatz yields good agreement with previous density-matrix renormalization group results in one dimension and is able to accurately describe the spin-Peierls transition. This variational approach is neither constrained by the magnetoelastic-coupling strength nor by the dimensionality of the systems considered, thus allowing future investigations in more general cases, which are relevant to spin-liquid and topological phases in two spatial dimensions. △ Less

Submitted 3 October, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 7 pages, 5 figures

Journal ref: Phys. Rev. B 102, 125149 (2020)

Gapless spin liquid and valence-bond solid in the $J_1-J_2$ Heisenberg model on the square lattice: insights from singlet and triplet excitations

Authors: Francesco Ferrari, Federico Becca

Abstract: The spin-$1/2$ $J_1-J_2$ Heisenberg model on the square lattice represents one of the simplest examples in which the effets of magnetic interactions may suppress magnetic order, eventually leading to a pure quantum phase with no local order parameters. This model has been extensively studied in the last three decades, with conflicting results. Here, by using Gutzwiller-projected wave functions and… ▽ More The spin-$1/2$ $J_1-J_2$ Heisenberg model on the square lattice represents one of the simplest examples in which the effets of magnetic interactions may suppress magnetic order, eventually leading to a pure quantum phase with no local order parameters. This model has been extensively studied in the last three decades, with conflicting results. Here, by using Gutzwiller-projected wave functions and recently developed methods to assess the low-energy spectrum, we show the existence of a level crossing between the lowest-energy triplet and singlet excitations for $J_2/J_1 \approx 0.54$. This fact supports the existence of a phase transition between a gapless spin liquid (which is stable for $0.48 \lesssim J_2/J_1 \lesssim 0.54$) and a valence-bond solid (for $0.54 \lesssim J_2/J_1 \lesssim 0.6$), even though no clear sign of dimer order is visible in the correlations functions. These results, which confirm recent density-matrix renormalization calculations on cylindrical clusters [L. Wang and A.W. Sandvik, Phys. Rev. Lett. 121, 107202 (2018)] reconcile the contraddicting results obtained within different approaches over the years. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. B 102, 014417 (2020)

Dynamical properties of Néel and valence-bond phases in the $J_1-J_2$ model on the honeycomb lattice

Authors: Francesco Ferrari, Federico Becca

Abstract: By using a variational Monte Carlo technique based upon Gutzwiller-projected fermionic states, we investigate the dynamical structure factor of the antiferromagnetic $S=1/2$ Heisenberg model on the honeycomb lattice, in presence of first-neighbor ($J_1$) and second-neighbor ($J_2$) couplings, for ${J_2 < 0.5 J_1}$. The ground state of the system shows long-range antiferromagnetic order for… ▽ More By using a variational Monte Carlo technique based upon Gutzwiller-projected fermionic states, we investigate the dynamical structure factor of the antiferromagnetic $S=1/2$ Heisenberg model on the honeycomb lattice, in presence of first-neighbor ($J_1$) and second-neighbor ($J_2$) couplings, for ${J_2 < 0.5 J_1}$. The ground state of the system shows long-range antiferromagnetic order for ${J_2/J_1 \lesssim 0.23}$, plaquette valence-bond order for ${0.23 \lesssim J_2/J_1 \lesssim 0.36}$, and columnar dimer order for ${J_2/J_1 \gtrsim 0.36}$. Within the antiferromagnetic state, a well-defined magnon mode is observed, whose dispersion is in relatively good agreement with linear spin-wave approximation for $J_2=0$. When a nonzero second-neighbor super-exchange is included, a roton-like mode develops around the $K$ point (i.e., the corner of the Brillouin zone). This mode softens when $J_2/J_1$ is increased and becomes gapless at the transition point, $J_2/J_1 \approx 0.23$. Here, a broad continuum of states is clearly visible in the dynamical spectrum, suggesting that nearly-deconfined spinon excitations could exist, at least at relatively high energies. For larger values of $J_2/J_1$, valence-bond order is detected and the spectrum of the system becomes clearly gapped, with a triplon mode at low energies. This is particularly evident for the spectrum of the dimer valence-bond phase, in which the triplon mode is rather well separated from the continuum of excitations that appears at higher energies. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Neural Gutzwiller-projected variational wave functions

Authors: Francesco Ferrari, Federico Becca, Juan Carrasquilla

Abstract: Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzw… ▽ More Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid Helium. More recently, Gutzwiller-projected wave functions, typically constructed from fermionic degrees of freedom, have been employed to examine quantum spin models in the presence of competing interactions, where exotic phases with no spontaneous symmetry breaking and fractional excitations may exist. In this work, we investigate the aforementioned fermionic wave functions supplemented with neural networks, specifically with the so-called restricted Boltzmann machine (RBM), to boost their accuracy and obtain reliable approximations to the ground state of generic spin models. In particular, we apply our neural augmented fermionic construction to the description of both magnetically ordered and disordered phases of increasing complexity, including cases where the ground state displays a non-trivial sign structure. Even though the RBM state is by far more effective for Néel states endowed with a particularly simple sign structure, it provides a significant improvement over the original fermionic state in highly frustrated regimes where a complex sign structure is anticipated, thus marking the path to an understanding of strongly-correlated spin models on the lattice via neural Gutzwiller-projected variational wave functions. △ Less

Submitted 13 September, 2019; v1 submitted 2 June, 2019; originally announced June 2019.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. B 100, 125131 (2019)

A new strategy to microscopic modelling of topological entanglement in polymers based on field theory

Authors: Franco Ferrari

Abstract: In this work a new strategy is proposed in order to build analytic and microscopic models of fluctuating polymer rings subjected to topological constraints. The topological invariants used to fix these constraints belong to a wide class of the so-called numerical topological invariants. For each invariant it is possible to derive a field theory that describes the statistical behavior of knotted an… ▽ More In this work a new strategy is proposed in order to build analytic and microscopic models of fluctuating polymer rings subjected to topological constraints. The topological invariants used to fix these constraints belong to a wide class of the so-called numerical topological invariants. For each invariant it is possible to derive a field theory that describes the statistical behavior of knotted and linked polymer rings following a straightforward algorithm. The treatment is not limited to the partition function of the system, but it allows also to express the expectation values of general observables as field theory amplitudes. Our strategy is illustrated taking as examples the Gauss linking number and a topological invariant belonging to a class of invariants due to Massey. The consistency of the new method developed here is checked by reproducing a previous field theoretical model of two linked polymer rings. After the passage to field theory, the original topological constraints imposed on the fluctuating paths of the polymers become constraints over the configurations of the topological fields that mediate the interactions of topological origin between the monomers. These constraints involve quantities like the cross-helicity which are of interest in other disciplines, like for instance in modeling the solar magnetic field. While the calculation of the vacuum expectation values of generic observables remains still challenging due to the complexity of the problem of topological entanglement in polymer systems, we succeed here to reduce the evaluation of the moments of the Gauss linking number for two linked polymer rings to the computation of the amplitudes of a free field theory. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: 34 pages, no figures

Phases Of Melonic Quantum Mechanics

Authors: Frank Ferrari, Fidel I. Schaposnik Massolo

Abstract: We explore in detail the properties of two melonic quantum mechanical theories which can be formulated either as fermionic matrix quantum mechanics in the new large $D$ limit, or as disordered models. Both models have a mass parameter $m$ and the transition from the perturbative large $m$ region to the strongly coupled "black-hole" small $m$ region is associated with several interesting phenomena.… ▽ More We explore in detail the properties of two melonic quantum mechanical theories which can be formulated either as fermionic matrix quantum mechanics in the new large $D$ limit, or as disordered models. Both models have a mass parameter $m$ and the transition from the perturbative large $m$ region to the strongly coupled "black-hole" small $m$ region is associated with several interesting phenomena. One model, with ${\rm U}(n)^2$ symmetry and equivalent to complex SYK, has a line of first-order phase transitions terminating, for a strictly positive temperature, at a critical point having non-trivial, non-mean-field critical exponents for standard thermodynamical quantities. Quasi-normal frequencies, as well as Lyapunov exponents associated with out-of-time-ordered four-point functions, are also singular at the critical point, leading to interesting new critical exponents. The other model, with reduced ${\rm U}(n)$ symmetry, has a quantum critical point at strictly zero temperature and positive critical mass $m_*$. For $0<m<m_*$, it flows to a new gapless IR fixed point, for which the standard scale invariance is spontaneously broken by the appearance of distinct scaling dimensions $Δ_+$ and $Δ_-$ for the Euclidean two-point function when $t\rightarrow +\infty$ and $t\rightarrow -\infty$ respectively. We provide several detailed and pedagogical derivations, including rigorous proofs or simplified arguments for some results that were already known in the literature. △ Less

Submitted 27 June, 2019; v1 submitted 15 March, 2019; originally announced March 2019.

Comments: 135 pages, 50 figures, 2 tables. Minor corrections and updated list of references in v3. Title changed to match published version

Journal ref: Phys. Rev. D 100, 026007 (2019)

Dynamical structure factor of the $J_1-J_2$ Heisenberg model on the triangular lattice: magnons, spinons, and gauge fields

Authors: Francesco Ferrari, Federico Becca

Abstract: Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials. However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and m… ▽ More Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials. However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and may be directly assessed by inelastic neutron scattering. In this work, we investigate the low-energy properties of the $S=1/2$ Heisenberg model on the triangular lattice, including both nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ super-exchanges, by a dynamical variational Monte Carlo approach that allows accurate results on spin models. For $J_2=0$, our calculations are compatible with the existence of a well-defined magnon in the whole Brillouin zone, with gapless excitations at $K$ points (i.e., at the corners of the Brillouin zone). The strong renormalization of the magnon branch (also including roton-like minima around the $M$ points, i.e., midpoints of the border zone) is described by our Gutzwiller-projected state, where Abrikosov fermions are subject to a non-trivial magnetic $π$-flux threading half of the triangular plaquettes. When increasing the frustrating ratio $J_2/J_1$, we detect a progessive softening of the magnon branch at $M$, which eventually becomes gapless within the spin-liquid phase. This feature is captured by the band structure of the unprojected wave function (with $2$ Dirac points for each spin component). In addition, we observe an intense signal at low energies around the $K$ points, which cannot be understood within the unprojected picture and emerges only when the Gutzwiller projection is considered, suggesting the relevance of gauge fields for the low-energy physics of spin liquids. △ Less

Submitted 19 July, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 12 pages, 9 figures

Journal ref: Phys. Rev. X 9, 031026 (2019)

Spectral signatures of fractionalization in the frustrated Heisenberg model on the square lattice

Authors: Francesco Ferrari, Federico Becca

Abstract: We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio $J_2/J_1$, the ground state undergoes a continuous transition from a Néel antiferromagnet to a $\mathbb{Z}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases… ▽ More We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio $J_2/J_1$, the ground state undergoes a continuous transition from a Néel antiferromagnet to a $\mathbb{Z}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases and highlight the existence of a broad continuum of excitations in the proximity of the spin-liquid phase. The magnon branch, which dominates the spectrum of the unfrustrated Heisenberg model, gradually loses its spectral weight, thus releasing nearly-deconfined spinons, whose signatures are visible even in the magnetically ordered state. Our results show how free spinons emerge across a quantum critical point, providing evidence for the fractionalization of magnons into deconfined spinons. △ Less

Submitted 17 September, 2018; v1 submitted 23 May, 2018; originally announced May 2018.

Comments: 5 pages, 3 figures + supplemental materials

Journal ref: Phys. Rev. B 98, 100405 (2018)

Dynamical structure factor of the $J_1-J_2$ Heisenberg model in one dimension: the variational Monte Carlo approach

Authors: Francesco Ferrari, Alberto Parola, Sandro Sorella, Federico Becca

Abstract: The dynamical spin structure factor is computed within a variational framework to study the one-dimensional $J_1-J_2$ Heisenberg model. Starting from Gutzwiller-projected fermionic wave functions, the low-energy spectrum is constructed from two-spinon excitations. The direct comparison with Lanczos calculations on small clusters demonstrates the excellent description of both gapless and gapped (di… ▽ More The dynamical spin structure factor is computed within a variational framework to study the one-dimensional $J_1-J_2$ Heisenberg model. Starting from Gutzwiller-projected fermionic wave functions, the low-energy spectrum is constructed from two-spinon excitations. The direct comparison with Lanczos calculations on small clusters demonstrates the excellent description of both gapless and gapped (dimerized) phases, also including incommensurate structures for $J_2/J_1>0.5$. Calculations on large clusters show how the intensity evolves when increasing the frustrating ratio and give an unprecedented accurate characterization of the dynamical properties of (non-integrable) frustrated spin models. △ Less

Submitted 7 June, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. B 97, 235103 (2018)

Phase Diagram of Planar Matrix Quantum Mechanics, Tensor, and Sachdev-Ye-Kitaev Models

Authors: Tatsuo Azeyanagi, Frank Ferrari, Fidel I. Schaposnik Massolo

Abstract: We compute the phase diagram of a $\text{U}(N)^{2}\times\text{O}(D)$ invariant fermionic planar matrix quantum mechanics [equivalently tensor or complex Sachdev-Ye-Kitaev (SYK) models] in the new large $D$ limit, dominated by melonic graphs. The Schwinger-Dyson equations can have two solutions describing either a high entropy, SYK black-hole-like phase, or a low entropy one with trivial IR behavio… ▽ More We compute the phase diagram of a $\text{U}(N)^{2}\times\text{O}(D)$ invariant fermionic planar matrix quantum mechanics [equivalently tensor or complex Sachdev-Ye-Kitaev (SYK) models] in the new large $D$ limit, dominated by melonic graphs. The Schwinger-Dyson equations can have two solutions describing either a high entropy, SYK black-hole-like phase, or a low entropy one with trivial IR behavior. In the strongly coupled region of the mass-temperature plane, there is a line of first order phase transitions between the high and low entropy phases. This line terminates at a new critical point which we study numerically in detail. The critical exponents are nonmean field and differ on the two sides of the transition. We also study purely bosonic unstable and stable melonic models. The former has a line of Kazakov critical points beyond which the Schwinger-Dyson equations do not have a consistent solution. Moreover, in both models the would-be SYK-like solution of the IR limit of the equations does not exist in the full theory. △ Less

Submitted 8 February, 2018; v1 submitted 11 July, 2017; originally announced July 2017.

Comments: 14 pages, 5 figures, minor changes in v2

Journal ref: Phys. Rev. Lett. 120, 061602 (2018)

Competition between spin liquids and valence-bond order in the frustrated spin-$1/2$ Heisenberg model on the honeycomb lattice

Authors: Francesco Ferrari, Samuel Bieri, Federico Becca

Abstract: Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor $J_1$ and second-neighbor $J_2$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. Néel order is stabiliz… ▽ More Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor $J_1$ and second-neighbor $J_2$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. Néel order is stabilized for small values of the frustrating second-neighbor coupling. Increasing the ratio $J_2/J_1$, we find strong evidence for a continuous transition to a nonmagnetic phase at $J_2/J_1 \approx 0.23$. Close to the transition point, the Gutzwiller-projected uniform resonating valence bond state gives an excellent approximation to the exact ground-state energy. For $0.23 \lesssim J_2/J_1 \lesssim 0.4$, a gapless $Z_2$ spin liquid with Dirac nodes competes with a plaquette valence-bond solid. In contrast, the gapped spin liquid considered in previous works has significantly higher variational energy. Although the plaquette valence-bond order is expected to be present as soon as the Néel order melts, this ordered state becomes clearly favored only for $J_2/J_1 \gtrsim 0.3$. Finally, for $0.36 \lesssim J_2/J_1 \le 0.5$, a valence-bond solid with columnar order takes over as the ground state, being also lower in energy than the magnetic state with collinear order. We perform a detailed finite-size scaling and standard data collapse analysis, and we discuss the possibility of a deconfined quantum critical point separating the Néel antiferromagnet from the plaquette valence-bond solid. △ Less

Submitted 5 October, 2017; v1 submitted 23 June, 2017; originally announced June 2017.

Comments: 10 pages, 12 figures

Journal ref: Phys. Rev. B 96, 104401 (2017)

Block copolymer knots

Authors: Franco Ferrari

Abstract: An extensive study of single block copolymer knots containing two kinds of monomers $A$ and $B$ is presented. The knots are in a solution and their monomers are subjected to short range interactions that can be attractive or repulsive. In view of possible applications in medicine and the construction of intelligent materials, it is shown that several features of copolymer knots can be tuned by cha… ▽ More An extensive study of single block copolymer knots containing two kinds of monomers $A$ and $B$ is presented. The knots are in a solution and their monomers are subjected to short range interactions that can be attractive or repulsive. In view of possible applications in medicine and the construction of intelligent materials, it is shown that several features of copolymer knots can be tuned by changing the monomer configuration. A very fast and abrupt swelling with increasing temperature is obtained in certain multiblock copolymers, while the size and the swelling behavior at high temperatures may be controlled in diblock copolymers. Interesting new effects appear in the thermal diagrams of copolymer knots when their length is increased. △ Less

Submitted 6 April, 2016; v1 submitted 29 March, 2016; originally announced March 2016.

Comments: 5 pages, 6 figures, LaTeX + RevTeX 4, minor improvements in the text, a new reference added

The Analytic Renormalization Group

Authors: Frank Ferrari

Abstract: Finite temperature Euclidean two-point functions in quantum mechanics or quantum field theory are characterized by a discrete set of Fourier coefficients $G_{k}$, $k\in\mathbb Z$, associated with the Matsubara frequencies $ν_{k}=2πk/β$. We show that analyticity implies that the coefficients $G_{k}$ must satisfy an infinite number of model-independent linear equations that we write down explicitly.… ▽ More Finite temperature Euclidean two-point functions in quantum mechanics or quantum field theory are characterized by a discrete set of Fourier coefficients $G_{k}$, $k\in\mathbb Z$, associated with the Matsubara frequencies $ν_{k}=2πk/β$. We show that analyticity implies that the coefficients $G_{k}$ must satisfy an infinite number of model-independent linear equations that we write down explicitly. In particular, we construct "Analytic Renormalization Group" linear maps $\mathsf A_μ$ which, for any choice of cut-off $μ$, allow to express the low energy Fourier coefficients for $|ν_{k}|<μ$ (with the possible exception of the zero mode $G_{0}$), together with the real-time correlators and spectral functions, in terms of the high energy Fourier coefficients for $|ν_{k}|\geqμ$. Operating a simple numerical algorithm, we show that the exact universal linear constraints on $G_{k}$ can be used to systematically improve any random approximate data set obtained, for example, from Monte-Carlo simulations. Our results are illustrated on several explicit examples. △ Less

Submitted 16 August, 2016; v1 submitted 23 February, 2016; originally announced February 2016.

Comments: 52 pages, 25 figures; v2: a few comments and explanations added

Report number: CERN-TH-2016-023

Black Hole Horizons and Bose-Einstein Condensation

Authors: Frank Ferrari

Abstract: Consider a particle sitting at a fixed position outside of a stable black hole. If the system is heated up, the black hole horizon grows and there should exist a critical temperature above which the particle enters the black hole interior. We solve a simple model describing exactly this situation: a large N matrix quantum mechanics modeling a fixed D-particle in a black hole background. We show th… ▽ More Consider a particle sitting at a fixed position outside of a stable black hole. If the system is heated up, the black hole horizon grows and there should exist a critical temperature above which the particle enters the black hole interior. We solve a simple model describing exactly this situation: a large N matrix quantum mechanics modeling a fixed D-particle in a black hole background. We show that indeed a striking phenomenon occurs: above some critical temperature, there is a non-perturbative Bose-Einstein condensation of massless strings. The transition, even though precisely defined by the presence of the condensate, cannot be sharply detected by measurements made in a finite amount of time. The order parameter is fundamentally non-local in time and corresponds to infinite-time correlations. △ Less

Submitted 29 January, 2016; originally announced January 2016.

Comments: 11 pages, 1 figure

Report number: CERN-TH-2016-021

The topological effect on the Mechanical properties of polymer knots

Authors: Yani Zhao, Franco Ferrari

Abstract: The mechanical properties of polymer knots under stretching in a bad or good solvent are investigated by applying a given force $F$ to a point of the knot while keeping another point fixed. The Monte Carlo sampling of the polymer conformations on a simple cubic lattice is performed using a variant of the Wang-Landau algorithm. The results of the calculations of the specific energy, specific heat c… ▽ More The mechanical properties of polymer knots under stretching in a bad or good solvent are investigated by applying a given force $F$ to a point of the knot while keeping another point fixed. The Monte Carlo sampling of the polymer conformations on a simple cubic lattice is performed using a variant of the Wang-Landau algorithm. The results of the calculations of the specific energy, specific heat capacity and gyration radius for several knot topologies show a general trend in the behavior of short polymer knots with lengths up to seventy lattice units. At low tensile force $F$, knots can be found either in a compact or an extended phase, depending if the temperature is low or high. At any temperature, with increasing values of the force $F$, a polymer knot undergoes a phase transition to a stretched state. This transition is characterized by a strong peak in the heat capacity. There is also a minor peak, which corresponds to a transition occurring at low temperatures when the conformations of polymers in the stretched phase become swollen with increasing temperatures. It is also shown that the behavior of short polymer rings is strongly influenced by topological effects. The limitations in the number of accessible energy states due to topological constraints is particularly evident in knots of small size and such that their minimum number of crossing according to the Rolfsen knot table is high. An example is provided by a cinquefoil knot $5_1$ with a length of only fifty lattice units. The thermal and mechanical properties of knots that can be represented with diagrams having the same minimum number of crossings, are very similar. The size effects on the behavior of polymer knots have been analyzed too. Surprisingly, it is found that topological effects fade out very fast with increasing polymer length. △ Less

Submitted 13 December, 2015; originally announced December 2015.

Comments: 13 pages, 12 figures

Cubic-scaling iterative solution of the Bethe-Salpeter equation for finite systems

Authors: Mathias P. Ljungberg, Peter Koval, Francesco Ferrari, Dietrich Foerster, Daniel Sanchez-Portal

Abstract: The Bethe-Salpeter equation (BSE) is currently the state of the art in the description of neutral electron excitations in both solids and large finite systems. It is capable of accurately treating charge-transfer excitations that present difficulties for simpler approaches. We present a local basis set formulation of the BSE for molecules where the optical spectrum is computed with the iterative H… ▽ More The Bethe-Salpeter equation (BSE) is currently the state of the art in the description of neutral electron excitations in both solids and large finite systems. It is capable of accurately treating charge-transfer excitations that present difficulties for simpler approaches. We present a local basis set formulation of the BSE for molecules where the optical spectrum is computed with the iterative Haydock recursion scheme, leading to a low computational complexity and memory footprint. Using a variant of the algorithm we can go beyond the Tamm-Dancoff approximation (TDA). We rederive the recursion relations for general matrix elements of a resolvent, show how they translate into continued fractions, and study the convergence of the method with the number of recursion coefficients and the role of different terminators. Due to the locality of the basis functions the computational cost of each iteration scales asymptotically as $O(N^3)$ with the number of atoms, while the number of iterations is typically much lower than the size of the underlying electron-hole basis. In practice we see that , even for systems with thousands of orbitals, the runtime will be dominated by the $O(N^2)$ operation of applying the Coulomb kernel in the atomic orbital representation △ Less

Submitted 20 May, 2015; originally announced May 2015.

FQHE on curved backgrounds, free fields and large N

Authors: Frank Ferrari, Semyon Klevtsov

Abstract: We study the free energy of the Laughlin state on curved backgrounds, starting from the free field representation. A simple argument, based on the computation of the gravitational effective action from the transformation properties of Green functions under the change of the metric, allows to compute the first three terms of the expansion in large magnetic field. The leading and subleading contribu… ▽ More We study the free energy of the Laughlin state on curved backgrounds, starting from the free field representation. A simple argument, based on the computation of the gravitational effective action from the transformation properties of Green functions under the change of the metric, allows to compute the first three terms of the expansion in large magnetic field. The leading and subleading contributions are given by the Aubin-Yau and Mabuchi functionals respectively, whereas the Liouville action appears at next-to-next-to-leading order. We also derive a path integral representation for the remainder terms. They correspond to a large mass expansion for a related interacting scalar field theory and are thus given by local polynomials in curvature invariants. △ Less

Submitted 13 December, 2014; v1 submitted 24 October, 2014; originally announced October 2014.

Comments: 14 pages; v3: conformal spin rescaled, minor changes

Journal ref: JHEP12(2014)086

On a field theoretical model of polymeric $2s-$plats and some of its consequences

Authors: Franco Ferrari, Jarosław Paturej, Marcin Piątek, Yani Zhao

Abstract: The field theory approach to the statistical mechanics of a system of N polymer rings linked together is generalized to the case of links that have a fixed number $2s$ of maxima and minima. Such kind of links are called plats and appear for instance in the DNA of living organisms. The topological states of the link are distinguished using the Gauss linking number. This is a relatively weak link in… ▽ More The field theory approach to the statistical mechanics of a system of N polymer rings linked together is generalized to the case of links that have a fixed number $2s$ of maxima and minima. Such kind of links are called plats and appear for instance in the DNA of living organisms. The topological states of the link are distinguished using the Gauss linking number. This is a relatively weak link invariant in the case of a general link, but its efficiency improves when $2s-$plats are considered. It is proved that, if we restrict ourselves to $2s-$plat conformations, the field theoretical model established here is able to take into account also the interactions of topological origin involving three chains simultaneously. It is shown that these three-body interactions have nonvanishing contributions when three or more rings are entangled together, enhancing for instance the attractive forces between monomers. The model can be used to study the statistical mechanics of polymers in confined geometries, for instance when $2s$ extrema of a few polymer rings are attached to membranes. Its partition function is mapped here into that of a multi-layer electron gas. Such quasi-particle systems are studied in connection with several interesting applications, including high-$T_c$ superconductivity and topological quantum computing. At the end an useful connection with the cosh-Gordon equation is shown. △ Less

Submitted 9 July, 2018; v1 submitted 30 September, 2014; originally announced September 2014.

Comments: 35 pages, 8 figures. arXiv admin note: substantial text overlap with arXiv:1304.5980

Ring polymers with topological constraints

Authors: Yani Zhao, Franco Ferrari

Abstract: In the first part of this work a summary is provided of some recent experiments and theoretical results which are relevant in the research of systems of polymer rings in nontrivial topological conformations. Next, some advances in modeling the behavior of single polymer knots are presented. The numerical simulations are performed with the help of the Wang-Landau Monte Carlo algorithm. To sample th… ▽ More In the first part of this work a summary is provided of some recent experiments and theoretical results which are relevant in the research of systems of polymer rings in nontrivial topological conformations. Next, some advances in modeling the behavior of single polymer knots are presented. The numerical simulations are performed with the help of the Wang-Landau Monte Carlo algorithm. To sample the polymer conformation a set of random transformations called pivot moves is used. The crucial problem of preserving the topology of the knots after each move is tackled with the help of two new techniques which are briefly explained. As an application, the results of an investigation of the effects of topology on the thermal properties of polymer knots is reported. In the end, original results are discussed concerning the use of parallelized codes to study polymers knots composed by a large number of segments within the Wang-Landau approach. △ Less

Submitted 2 February, 2014; originally announced February 2014.

Comments: 14 pages, 7 figures, proceedings of the conference ISMC 2013, 15-19 September 2013, Rome, Italy

A study of polymer knots using a simple knot invariant written consisting of multiple contour integrals

Authors: Yani Zhao, Franco Ferrari

Abstract: In this work the thermodynamic properties of short polymer knots (up to 120 segments) defined on a simple cubic lattice are studied with the help of the Wang-Landau Monte Carlo algorithm. The sampling process is performed using pivot transformations starting from a given seed conformation. Both cases of short-range attractive and repulsive interactions acting on the monomers are considered. The pr… ▽ More In this work the thermodynamic properties of short polymer knots (up to 120 segments) defined on a simple cubic lattice are studied with the help of the Wang-Landau Monte Carlo algorithm. The sampling process is performed using pivot transformations starting from a given seed conformation. Both cases of short-range attractive and repulsive interactions acting on the monomers are considered. The properties of the specific energy, heat capacity and gyration radius of several knots are discussed. It is found that the heat capacity exhibits a sharp peak. If the interactions are attractive, similar peaks have been observed also in single open chains and have been related to the transition from a frozen crystallite state to an expanded coil state. Some other peculiarities of the behavior of the analyzed observables are presented, like for instance the increasing or decreasing of the knot specific energy at high temperatures with increasing polymer lengths depending if the interactions are attractive or repulsive. Besides the investigation of the thermodynamics of polymer knots, the second goal of this paper is to introduce a method for distinguishing the topology of a knot based on a topological invariant which is in the form of multiple contour integrals and explicitly depends on the physical trajectory of the knot. The chosen invariant is related to the second coefficient of the Conway polynomial. It has been first isolated from the amplitudes of a Chern-Simons field theory with gauge group SU(N). It is shown that this invariant is very reliable in distinguishing the topology of polymer knots. One of the advantages of the proposed approach is that it allows to reduce the number of samples needed by the Wang-Landau algorithm. Some solutions to speed up the calculations exploiting Monte Carlo integration techniques are developed. △ Less

Submitted 3 July, 2013; v1 submitted 24 June, 2013; originally announced June 2013.

Comments: 18 pages, 15 figures, LaTeX + RevTeX 4, a small misprint in fig. 9 corrected, some new references added

Polymer plats and multicomponent anyon gases

Authors: Franco Ferrari, Jarosław Paturej, Marcin Piątek, Yani Zhao

Abstract: Anyon systems are studied in connection with several interesting applications including high $T_C$ superconductivity and topological quantum computing. In this work we show that these systems can be realized starting from directed polymers braided together to form a nontrivial link configuration belonging to the topological class of plats. The statistical sum of a such plat is related here to the… ▽ More Anyon systems are studied in connection with several interesting applications including high $T_C$ superconductivity and topological quantum computing. In this work we show that these systems can be realized starting from directed polymers braided together to form a nontrivial link configuration belonging to the topological class of plats. The statistical sum of a such plat is related here to the partition function of a two-component anyon gas. The constraints that preserve the topological configuration of the plat are imposed on the polymer trajectories using the so-called Gauss linking number, a topological invariant that has already been well studied in polymer physics. Due to these constraints, short-range forces act on the monomers or, equivalently, on the anyon quasiparticles in a way that closely resembles the appearance of reaction forces in the constrained systems of classical mechanics. If the polymers are homogeneous, the anyon system reaches a self-dual point, in which these forces vanish exactly. A class of self-dual solutions that minimize the energy of the anyons is derived. The two anyon gas discussed here obeys an abelian statistics, while for quantum computing it is known that nonabelian anyons are necessary. However, this is a limitation due to the use of the Gauss linking invariant to impose the topological constraints, which is a poor topological invariant and is thus unable to capture the nonabelian characteristics of the braided polymer chains. A more refined treatment of the topological constraints would require more sophisticated topological invariants, but so far their application to the statistical mechanics of linked polymers is an open problem. △ Less

Submitted 28 April, 2013; v1 submitted 22 April, 2013; originally announced April 2013.

Comments: title changed, several references added, typos corrected, text improved. 42 pages, 5 figures

A numerical technique for preserving the topology of polymer knots: The case of short-range attractive interactions

Authors: Yani Zhao, Franco Ferrari

Abstract: The statistical mechanics of single polymer knots is studied using Monte Carlo simulations. The polymers are considered on a cubic lattice and their conformations are randomly changed with the help of pivot transformations. After each transformation, it is checked if the topology of the knot is preserved by means of a method called pivot algorithm and excluded area (in short PAEA) and described in… ▽ More The statistical mechanics of single polymer knots is studied using Monte Carlo simulations. The polymers are considered on a cubic lattice and their conformations are randomly changed with the help of pivot transformations. After each transformation, it is checked if the topology of the knot is preserved by means of a method called pivot algorithm and excluded area (in short PAEA) and described in a previous publication of the authors. As an application of this method the specific energy, the radius of gyration and heat capacity of a few types of knots are computed. The case of attractive short-range forces is investigated. The sampling of the energy states is performed by means of the Wang-Landau algorithm. The obtained results show that the specific energy and heat capacity increase with increasing knot complexity as in the case of repulsive interactions. The data about the gyration radius allow to estimate the size of the polymer knots at different temperatures. △ Less

Submitted 28 December, 2012; originally announced December 2012.

Comments: 8 pages, 2 figures, LaTeX + appolb.cls, Proceedings of the 25th Marian Smoluchowski Symposium on Statistical Physics, 9--13 September 2012, Cracow

Topological effects in the thermal properties of knotted polymer rings

Authors: Yani Zhao, Franco Ferrari

Abstract: The topological effects on the thermal properties of several knot configurations are investigated using Monte Carlo simulations. In order to check if the topology of the knots is preserved during the thermal fluctuations we propose a method that allows very fast calculations and can be easily applied to arbitrarily complex knots. As an application, the specific energy and heat capacity of the tref… ▽ More The topological effects on the thermal properties of several knot configurations are investigated using Monte Carlo simulations. In order to check if the topology of the knots is preserved during the thermal fluctuations we propose a method that allows very fast calculations and can be easily applied to arbitrarily complex knots. As an application, the specific energy and heat capacity of the trefoil, the figure-eight and the $8_1$ knots are calculated at different temperatures and for different lengths. Short-range repulsive interactions between the monomers are assumed. The knots configurations are generated on a three-dimensional cubic lattice and sampled by means of the Wang-Landau algorithm and of the pivot method. The obtained results show that the topological effects play a key role for short-length polymers. Three temperature regimes of the growth rate of the internal energy of the system are distinguished. △ Less

Submitted 24 September, 2012; v1 submitted 8 September, 2012; originally announced September 2012.

Comments: 7 pages, 12 figures, LaTeX + RevTeX. With respect to the first version, in the second version the text has been improved and all figures are now in black and white