SciPost Submission Page

Full counting statistics for interacting trapped fermions

by Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Grégory Schehr

Submission summary

As Contributors: Gregory Schehr
Preprint link: scipost_202106_00033v1
Date submitted: 2021-06-18 17:34
Submitted by: Schehr, Gregory
Submitted to: SciPost Physics
Academic field: Physics
Specialties:
  • Mathematical Physics
Approach: Theoretical

Abstract

We study $N$ spinless fermions in their ground state confined by an external potential in one dimension with long range interactions of the general Calogero-Sutherland type. For some choices of the potential this system maps to standard random matrix ensembles for general values of the Dyson index $\beta$. In the fermion model $\beta$ controls the strength of the interaction, $\beta=2$ corresponding to the noninteracting case. We study the quantum fluctuations of the number of fermions ${\cal N}_{\cal D}$ in a domain $\cal{D}$ of macroscopic size in the bulk of the Fermi gas. We predict that for general $\beta$ the variance of ${\cal N}_{\cal D}$ grows as $A_{\beta} \log N + B_{\beta}$ for $N \gg 1$ and we obtain a formula for $A_\beta$ and $B_\beta$. This is based on an explicit calculation for $\beta\in\left\{ 1,2,4\right\} $ and on a conjecture that we formulate for general $\beta$. This conjecture further allows us to obtain a formula for the higher cumulants of ${\cal N}_{\cal D}$. Our results for the variance in the microscopic regime are found to be consistent with the predictions of the Luttinger liquid theory with parameter $K = 2/\beta$, and allow to go beyond. In addition we present families of interacting fermion models in one dimension which, in their ground states, can be mapped onto random matrix models. We obtain the mean fermion density for these models for general interaction parameter $\beta$. In some cases the fermion density exhibits interesting transitions, for example we obtain a noninteracting fermion formulation of the Gross-Witten-Wadia model.

Current status:
Editor-in-charge assigned



Reports on this Submission

Anonymous Report 3 on 2021-8-26 (Invited Report)

Strengths

1. Many new analytical results supported by numerical evidence
2. Puzzling behaviour of higher order cumulants of FCS is found
3. Nice introduction to the link between RMT and Calogero-Sutherland model including free fermions
4. Intuitive but rigorous derivations.

Weaknesses

Not so many, the authors may consider adding a paragraph or two with estimates of various energy scales when discussing universality of higher cumulants.

Report

The manuscript reports new results on Full Counting Statistics of free
fermions and Calogero-Sutherland particles in their ground state. The authors
give results for cumulants of distribution of number of particles in an
interval (finite or infinite) in the thermodynamic limit. Some results are
exact, obtained using mapping of the ground state probability to the joint
probability of eigenvalues of random matrices ensembles. There are also many
conjectured results supported by numerical evidence.

I find the manuscript very well written and useful as a reference for
calculation (at least as a zero order approximation) for strongly interacting
many body systems in one dimension. The authors do great job in explaining the
mapping between random matrices and quantum particles. The previous results
are clearly referenced.

I found the universality of higher cumulants particularly interesting. I think
it deserves more discussion, so the authors may try to add some explanation
with a back of envelope estimates of the energy scales. They could probably
speculate whether such universality will hold for interacting systems other
than Calogero-Sutherland. I also find Eq. (20) a bit puzzling as there is no
dependence on the interval size, unlike in the variance.

I summary I think this is a very good paper with excellent presentation of
many new interesting results and theoretical techniques. While being slightly
on the side of mathematical statistical physics it can benefit
experimentalists in strongly correlated systems.

  • validity: high
  • significance: high
  • originality: high
  • clarity: high
  • formatting: excellent
  • grammar: excellent

Anonymous Report 2 on 2021-8-25 (Invited Report)

Strengths

explicit results on FCS for interacting fermions

Weaknesses

some of the models considered are very fine-tuned

Report

The manuscript deals with the study of counting statistics in the
ground state of certain interacting 1D fermion systems.
In particular, they consider model systems with two-body interactions,
where a mapping onto a certain random matrix ensemble exists.
The main example is the Gaussian $\beta$ ensemble, which corresponds
to the Calogero-Sutherland model in the sense that the joint PDF
describing the random matrix eigenvalues is identical to the square
of the fermionic many-body wave function under an appropriate variable
substitution. Similar mappings can be found for the circular,
Wishart-Laguerre and Jacobi $\beta$ ensembles, which describe fermions
on various domains with generalized interaction and potential terms.

The main goal of the paper is to find, using previous RMT results
obtained via the powerful Coulomb gas technique, a formula that relates
the fermion counting statistics to the noninteracting ($\beta=2$) case.
In the Gaussian case the authors find the result (14) for the variance
in a box, which applies to $\beta=1,2,4$. This shows that, after a
proper rescaling of the lengths, the main difference w.r.t. the
free-fermion case is a factor $2/\beta$ that multiplies the leading
order logarithmic term. Furthermore, the authors also propose a more
general formula (18) that is expected to hold for each of the ensembles
and for arbitrary $\beta$. A nontrivial conjecture of the subleading
term is presented in (17) and tested numerically later on in Sec. 2.
Results for higher order cumulants are also presented in (20)-(21).
Finally, some arguments on more general RM ensembles are discussed.

I believe that the manuscript contains novel and interesting results,
extending the studies on FCS towards interacting fermionic models.
Even though some of the models discussed are, due to the mapping to
a concrete RMT ensemble, very fine tuned, these could be important
starting points for further studies. I thus recommend the publication
of the manuscript and have only a few questions.

Requested changes

1.
The conjecture on higher order cumulants is not tested numerically.
Is it more demanding to check, e.g. requires larger matrix statistics?

2.
The statement about the variance at the end of Sec. 6 was not clear to me.
Does (18) with an appropriate rescaling and the SAME constant (17) work
also for the models in Table 2?

3.
It is known that for non-interacting fermions the FCS becomes universal
under varying the potential term in appropriate bulk and edge regimes.
Could one expect a similar kind of universality when perturbing the
interaction term away from the ones that correspond to random matrix
ensembles?

4. Typos
- sign convention in (31) and (32) seem to differ
- bottom of p. (18): "matches exactly with the formula (58)"
seems to be the wrong reference

  • validity: high
  • significance: good
  • originality: good
  • clarity: high
  • formatting: excellent
  • grammar: excellent

Anonymous Report 1 on 2021-8-1 (Invited Report)

Report

The authors consider N spinless fermions with long-range interactions and confined by an external potential in one dimension. They focus on systems which admit a mapping to random matrix ensembles, which they use to study the full counting statistics of the particle number in the ground state.

The manuscript contains many results which, to the best of my knowledge, are new, and could be summarized as follows. First, by following an approach pioneered by Sutherland and Calogero, they identify a set of models for which a mapping to random matrix theory (RMT) exists. By construction, for these models the ground-state wave function can be written down exactly, but they are not necessarily integrable. Next, they present a thorough analysis of the full counting statistics in these models, and study: (i) the mean density; (ii) the variance of the number of fermions in a given interval; (iii) higher cumulants beyond the variance. Among the main results, they present a conjecture for the variance which they prove analytically in some special cases and show that it passes highly non-trivial tests based on previous mathematical results and numerical computations. Furthermore, they substantiate a conjecture for the higher cumulants, stating that they are determined solely from the microscopic scale, thus being independent of the size of the interval and in fact universal. This generalizes a conjecture put forward in a previous work for non-interacting fermions.

Overall, I think the paper is very strong. While some of the material is a natural generalization of earlier work, I believe the paper contains a lot non-trivial new results, and represents a significant step forward in this branch of the literature.

Furthermore, the paper is very well written. I believe the authors make a very good job in the first section, guiding the reader through what was previously known and what are the new contributions of this work. I think this part is very easy to read even by non-experts.

Finally, I believe the research presented is timely, and of potential interest for a broad audience.

For these reasons, I recommend publication of the manuscript essentially as is.

I have only one comment about a minor aspect of the work. In Sec. 1.2, after Eq. (9), the authors compare the results they obtain against LDA, and conclude that the latter is only accurate for noninteracting fermions. However, one could expect this to be true by construction, since for the LDA prediction they consider noninteracting fermions. On the other hand, when the model under consideration happens to be integrable one could perform a slightly more sophisticated version of LDA, by combining it with Bethe Ansatz. In this case, one still assumes the local chemical potential to be given by \mu-V(x), but the density is then determined by the corresponding exact Bethe Ansatz solution for the interacting model. For locally-interacting systems, this is known to work in the bulk, so one could expect that the same is true here. The authors may optionally comment on this.

  • validity: -
  • significance: -
  • originality: -
  • clarity: -
  • formatting: -
  • grammar: -

Login to report or comment