Can rooted staggered fermions describe nonzero baryon density at low temperatures?

Szabolcs Borsányi (Department of Physics, Wuppertal University, Gaußstraße 20, D-42119, Wuppertal, Germany) ; Zoltán Fodor (Department of Physics, Wuppertal University, Gaußstraße 20, D-42119, Wuppertal, Germany; Institute for Theoretical Physics, ELTE Eötvös Loránd University, Pázmány P. sétány 1/A, H-1117 Budapest, Hungary; Pennsylvania State University, Department of Physics, State College, Pennsylvania 16801, USA; Physics Department, UCSD, San Diego, California 92093, USA; Jülich Supercomputing Centre, Forschungszentrum Jülich, D-52425 Jülich, Germany) ; Matteo Giordano (Institute for Theoretical Physics, ELTE Eötvös Loránd University, Pázmány P. sétány 1/A, H-1117 Budapest, Hungary) ; Jana N. Guenther (Department of Physics, Wuppertal University, Gaußstraße 20, D-42119, Wuppertal, Germany) ; Sándor D. Katz (Institute for Theoretical Physics, ELTE Eötvös Loránd University, Pázmány P. sétány 1/A, H-1117 Budapest, Hungary) ; et al. - Show all 7 authors

Research on the quantum chromodynamics (QCD) phase diagram with lattice field theory methods is dominated by the use of rooted staggered fermions, as they are the computationally cheapest discretization available. We show that rooted staggered fermions at a nonzero baryochemical potential μB predict a sharp rise in the baryon density at low temperatures and μB3mπ/2, where mπ is the Goldstone pion mass. We elucidate the nature of the nonanalyticity behind this sharp rise in the density by a comparison of reweighting results with a Taylor expansion of high order. While at first sight this nonanalytic behavior becomes apparent at the same position where the pion condensation transition takes place in the phase-quenched theory, the nature of the nonanalyticity in the two theories appears to be quite different: While at nonzero isospin density the data are consistent with a genuine thermodynamic (branch-point) singularity, the results at nonzero baryon density point to an essential singularity at μB=0. The effect is absent for four flavors of degenerate quarks, where rooting is not used. For the two-flavor case, we show numerical evidence that the magnitude of the effect diminishes on finer lattices. We discuss the implications of this technical complication on future studies of the QCD phase diagram.

{
  "_oai": {
    "updated": "2024-03-26T00:30:29Z", 
    "id": "oai:repo.scoap3.org:84263", 
    "sets": [
      "PRD"
    ]
  }, 
  "authors": [
    {
      "raw_name": "Szabolcs Bors\u00e1nyi", 
      "affiliations": [
        {
          "country": "Germany", 
          "value": "Department of Physics, Wuppertal University, Gau\u00dfstra\u00dfe 20, D-42119, Wuppertal, Germany"
        }
      ], 
      "surname": "Bors\u00e1nyi", 
      "given_names": "Szabolcs", 
      "full_name": "Bors\u00e1nyi, Szabolcs"
    }, 
    {
      "raw_name": "Zolt\u00e1n Fodor", 
      "affiliations": [
        {
          "country": "Germany", 
          "value": "Department of Physics, Wuppertal University, Gau\u00dfstra\u00dfe 20, D-42119, Wuppertal, Germany"
        }, 
        {
          "country": "Hungary", 
          "value": "Institute for Theoretical Physics, ELTE E\u00f6tv\u00f6s Lor\u00e1nd University, P\u00e1zm\u00e1ny P. s\u00e9t\u00e1ny 1/A, H-1117 Budapest, Hungary"
        }, 
        {
          "country": "USA", 
          "value": "Pennsylvania State University, Department of Physics, State College, Pennsylvania 16801, USA"
        }, 
        {
          "country": "USA", 
          "value": "Physics Department, UCSD, San Diego, California 92093, USA"
        }, 
        {
          "country": "Germany", 
          "value": "J\u00fclich Supercomputing Centre, Forschungszentrum J\u00fclich, D-52425 J\u00fclich, Germany"
        }
      ], 
      "surname": "Fodor", 
      "given_names": "Zolt\u00e1n", 
      "full_name": "Fodor, Zolt\u00e1n"
    }, 
    {
      "raw_name": "Matteo Giordano", 
      "affiliations": [
        {
          "country": "Hungary", 
          "value": "Institute for Theoretical Physics, ELTE E\u00f6tv\u00f6s Lor\u00e1nd University, P\u00e1zm\u00e1ny P. s\u00e9t\u00e1ny 1/A, H-1117 Budapest, Hungary"
        }
      ], 
      "surname": "Giordano", 
      "given_names": "Matteo", 
      "full_name": "Giordano, Matteo"
    }, 
    {
      "raw_name": "Jana N. Guenther", 
      "affiliations": [
        {
          "country": "Germany", 
          "value": "Department of Physics, Wuppertal University, Gau\u00dfstra\u00dfe 20, D-42119, Wuppertal, Germany"
        }
      ], 
      "surname": "Guenther", 
      "given_names": "Jana N.", 
      "full_name": "Guenther, Jana N."
    }, 
    {
      "raw_name": "S\u00e1ndor D. Katz", 
      "affiliations": [
        {
          "country": "Hungary", 
          "value": "Institute for Theoretical Physics, ELTE E\u00f6tv\u00f6s Lor\u00e1nd University, P\u00e1zm\u00e1ny P. s\u00e9t\u00e1ny 1/A, H-1117 Budapest, Hungary"
        }
      ], 
      "surname": "Katz", 
      "given_names": "S\u00e1ndor D.", 
      "full_name": "Katz, S\u00e1ndor D."
    }, 
    {
      "raw_name": "Attila P\u00e1sztor", 
      "affiliations": [
        {
          "country": "Hungary", 
          "value": "Institute for Theoretical Physics, ELTE E\u00f6tv\u00f6s Lor\u00e1nd University, P\u00e1zm\u00e1ny P. s\u00e9t\u00e1ny 1/A, H-1117 Budapest, Hungary"
        }, 
        {
          "country": "Hungary", 
          "value": "HUN-REN-ELTE Theoretical Physics Research Group, P\u00e1zm\u00e1ny P\u00e9ter s\u00e9t\u00e1ny 1/A, H-1117 Budapest, Hungary"
        }
      ], 
      "surname": "P\u00e1sztor", 
      "given_names": "Attila", 
      "full_name": "P\u00e1sztor, Attila"
    }, 
    {
      "raw_name": "Chik Him Wong", 
      "affiliations": [
        {
          "country": "Germany", 
          "value": "Department of Physics, Wuppertal University, Gau\u00dfstra\u00dfe 20, D-42119, Wuppertal, Germany"
        }
      ], 
      "surname": "Wong", 
      "given_names": "Chik Him", 
      "full_name": "Wong, Chik Him"
    }
  ], 
  "titles": [
    {
      "source": "APS", 
      "title": "Can rooted staggered fermions describe nonzero baryon density at low temperatures?"
    }
  ], 
  "dois": [
    {
      "value": "10.1103/PhysRevD.109.054509"
    }
  ], 
  "publication_info": [
    {
      "journal_volume": "109", 
      "journal_title": "Physical Review D", 
      "material": "article", 
      "journal_issue": "5", 
      "year": 2024
    }
  ], 
  "$schema": "http://repo.scoap3.org/schemas/hep.json", 
  "acquisition_source": {
    "date": "2024-03-26T00:30:23.193826", 
    "source": "APS", 
    "method": "APS", 
    "submission_number": "fae14820eb0711eeae4696b6a0e1ccbd"
  }, 
  "page_nr": [
    12
  ], 
  "license": [
    {
      "url": "https://creativecommons.org/licenses/by/4.0/", 
      "license": "CC-BY-4.0"
    }
  ], 
  "copyright": [
    {
      "statement": "Published by the American Physical Society", 
      "year": "2024"
    }
  ], 
  "control_number": "84263", 
  "record_creation_date": "2024-03-25T14:30:03.173156", 
  "_files": [
    {
      "checksum": "md5:637f7c376b93b2c7de1d3ed7d111b800", 
      "filetype": "pdf", 
      "bucket": "d0af9475-4484-4797-8554-1ac8cfa314dc", 
      "version_id": "88f54b2b-2177-4965-adf7-98dcfad2543f", 
      "key": "10.1103/PhysRevD.109.054509.pdf", 
      "size": 1336173
    }, 
    {
      "checksum": "md5:2784fe0231fd28d0b77845d3bf12ee66", 
      "filetype": "xml", 
      "bucket": "d0af9475-4484-4797-8554-1ac8cfa314dc", 
      "version_id": "45cd0397-5d72-4d88-9a72-80b2798d5005", 
      "key": "10.1103/PhysRevD.109.054509.xml", 
      "size": 173299
    }
  ], 
  "collections": [
    {
      "primary": "HEP"
    }, 
    {
      "primary": "Citeable"
    }, 
    {
      "primary": "Published"
    }
  ], 
  "arxiv_eprints": [
    {
      "categories": [
        "hep-lat"
      ], 
      "value": "2308.06105"
    }
  ], 
  "abstracts": [
    {
      "source": "APS", 
      "value": "Research on the quantum chromodynamics (QCD) phase diagram with lattice field theory methods is dominated by the use of rooted staggered fermions, as they are the computationally cheapest discretization available. We show that rooted staggered fermions at a nonzero baryochemical potential <math><msub><mi>\u03bc</mi><mi>B</mi></msub></math> predict a sharp rise in the baryon density at low temperatures and <math><msub><mi>\u03bc</mi><mi>B</mi></msub><mo>\u2273</mo><mn>3</mn><msub><mi>m</mi><mi>\u03c0</mi></msub><mo>/</mo><mn>2</mn></math>, where <math><msub><mi>m</mi><mi>\u03c0</mi></msub></math> is the Goldstone pion mass. We elucidate the nature of the nonanalyticity behind this sharp rise in the density by a comparison of reweighting results with a Taylor expansion of high order. While at first sight this nonanalytic behavior becomes apparent at the same position where the pion condensation transition takes place in the phase-quenched theory, the nature of the nonanalyticity in the two theories appears to be quite different: While at nonzero isospin density the data are consistent with a genuine thermodynamic (branch-point) singularity, the results at nonzero baryon density point to an essential singularity at <math><msub><mi>\u03bc</mi><mi>B</mi></msub><mo>=</mo><mn>0</mn></math>. The effect is absent for four flavors of degenerate quarks, where rooting is not used. For the two-flavor case, we show numerical evidence that the magnitude of the effect diminishes on finer lattices. We discuss the implications of this technical complication on future studies of the QCD phase diagram."
    }
  ], 
  "imprints": [
    {
      "date": "2024-03-25", 
      "publisher": "APS"
    }
  ]
}
Published on:
25 March 2024
Publisher:
APS
Published in:
Physical Review D , Volume 109 (2024)
Issue 5
DOI:
https://doi.org/10.1103/PhysRevD.109.054509
arXiv:
2308.06105
Copyrights:
Published by the American Physical Society
Licence:
CC-BY-4.0

Fulltext files: