SciPost Submission Page
Polaron spectroscopy of interacting Fermi systems: insights from exact diagonalization
by Ivan Amelio, Nathan Goldman
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Ivan Amelio 
Submission information  

Preprint Link:  https://arxiv.org/abs/2309.07019v1 (pdf) 
Date submitted:  20230914 14:56 
Submitted by:  Amelio, Ivan 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
Immersing a mobile impurity into a manybody quantum system represents a theoretically intriguing and experimentally effective way of probing its properties. In this work, we use exact diagonalization to compute the polaron spectral function in various instances of interacting Fermi settings, within the framework of FermiHubbard models. Inspired by possible realizations in cold atoms and semiconductor heterostructures, we consider different configurations for the background Fermi gas, including charge density waves, multiple Fermi seas and pair superfluids. Our results offer a benchmark for computations based on meanfield approaches and reveal surprising features of polaron spectra, inspiring new theoretical investigations.
Current status:
Reports on this Submission
Anonymous Report 3 on 20231017 (Invited Report)
 Cite as: Anonymous, Report on arXiv:2309.07019v1, delivered 20231017, doi: 10.21468/SciPost.Report.7959
Strengths
1) The authors present results for a bunch of different interacting systems.
2) Despite finitesize effects of exact diagonalization, the results show characteristic features of polaron physics. Thus the paper gives an important benchmark for further investigations.
Weaknesses
1) The method is limited to rather small system sizes. More "fewbody physics" is examined than "manybody physics".
2) Not all the features in the polaron spectra are suffiently discussed and understood.
Report
The article “Polaron spectroscopy of interacting Fermi systems: insights from exact diagonalization” by I. Amelio and N. Goldman presents spectral functions of polarons immersed in differently interacting 2D FermiHubbard lattices obtained from exact diagonalization.
By changing the parameters of the FermiHubbard model, the authors treat four different scenarios. At a spinpolarized background, they first benchmark the method on a noninteracting Fermi gas and then add repulsive Coulomb interactions giving rise to the formation of charge density waves. Next, they use a background of spinful fermions, where they study the effect of local repulsive and attractive Hubbard interactions to simulate the physics of transition metal dichalcogenides (TMD) and a polaron immersed in a superfluid, respectively. Here, they allow for spin imbalance and spindependent coupling between the impurity and the Fermi gas.
Although the method is limited to rather small system sizes, i.e., 38 particles on square and triangular lattices of sizes 4x4 up to 6x6, their spectra show characteristic features of polaron physics, i.e., the attractive and repulsive polaron branches or the signatures of umklapp scattering in presence of Coulomb repulsion. This is a valuable insight for the polaron community. However, their spectra obtain other features like the double line of the attractive polaron branch in a superfluid, which are not adequately explained. It remains an open question whether these are finitesize effects or really give hints to new physics.
Considering this and the following comments and questions, I am not able to yet give a final recommendation for whether to accept this article. To this end, I would kindly ask the authors to address the comments and questions listed under "Requested changes".
Requested changes
1) Exact diagonalization typically comes with an extreme restriction in terms of the system size. In my opinion, in the abstract it should be made clearer that the considered system sizes are extremely small and the manuscript actually deals with fewbody physics rather than manybody physics.
2) To get a better feeling of the numerical complexity, could the authors please specify how large is the considered Krylov space, i.e., what are typical values of M used in the computations?
3) The plots, especially those whose colors are scaled logarithmically, i.e., Fig. 1 and Fig. 2 (b), show a bunch of other lines beside the usual polaron branches. Have the authors investigated whether their number and shape is in relation to number of particles or lattice sites? Also, regarding Fig. 2 (b), how large is the strength of the peak related to umklapp scattering in comparison to the other peaks?
4) For the noninteracting gas and the charge density waves, the authors give results for both the square and the triangular lattice. However, the additional features in Fig. 3 (a) for a triangular lattice are not discussed extensively. Why is the repulsive polaron branch split in two at around $K = 10$ and why does the linearly rising second band have a bifurcation at around $K = 3$? TMD results, i.e., Fig. 4, are exclusively shown for triangular lattices whereas BCS results, i.e., Figs. 5 and 6, are only shown for square lattices. I do not follow the author’s strategy of when they show results for the square lattice and when for the triangular lattice. This is all the more confusing, as in typical TMD experiments the electrons predominantly explore the bands at the band minimum, i.e. for momenta where the dispersion relation is quadratic, and thus lattice effects (band warping etc.) does not matter. For a transparent and complete description, the authors should consider to present the other lattice’s results in an appendix or they just show results for the square lattice after showing once both in Fig. 1. Moreover, a discussion for why band effects should matter at all to simulate TMD physics of Fermi polarons should be included.
5) At the beginning of Sec. 4: I assume that K inside the quantum numbers (K,uparrow) etc. represents the momentum values k. I first confused it with the previously used interaction strength parameter $K = V_0/(a_\mathrm{CDW} E_F)$.
6) At the beginning of page 10, the authors remark that to identify the different lines for the attractive polarons in Fig. 4 they calculated the various binding energies. It would be more convincing if the authors show the binding energies in Fig. 4 as they already do that with cyan lines in Fig. 5.
7) In the Sec. 5 the authors compare their plots several times with those they got from a Chevy ansatz and meanfield calculations of polarons in a BCS superfluid. There, the spectra are given in terms of the density n instead of $E_F/E_B$. On Page 9, it is described that varying $U = (U_\uparrow + U_\downarrow) / 2$ corresponds to varying $E_F/E_B$ and thus is similar to varying the density. Maybe it is just me, but could you make this point clearer. What is the actual expression for the binding energies in terms of $U$? At the end of Sec. 2.1, the authors claim that it is “[their] philosophy th[e] method naturally needs to be complemented with other approaches, such as Chevy ansatz”, but for my taste there has not been made enough effort for a detailed comparison.
8) In the conclusion, it is mentioned that polaron spectroscopy is an effective way of probing manybody systems. Could the authors elaborate more on how the systems presented here can be realized in experiments? There is a great progress in simulating Hubbard models in optical lattices, also fewbody physics is being examined in coldatomic gases. Are there some specific challenges to simulate experimentally the scenarios used in this work?
9) There are two phrasings which are not really precise regarding physics:
a. Page 3: “noninteracting Fermi polaron” refers to polaron immersed in a noninteracting Fermi sea, however the impurityfermionic interaction $U$ needs to be nonzero for a polaron to be formed.
b. Page 6: I stumbled over the phrase “the impurity binds to an object of effective mass $N_\uparrow m_f$, where $m^*_f = \hbar^2/(2t_f)$ is the bare effective mass of a fermion”. I guess you mean “$N_\uparrow m^*_f$” and “bare effective mass of a fermion” is an somewhat confusing expression for “the effective mass of a single fermion” opposed to the “bare mass of a fermion”.
10) The authors are usually very precise when it comes to units. The absorption spectra have a dimension of 1/energy. Are the plotted absorption spectra given in units of $1/E_F$?
11) Finally, there are some defects in orthography and presentation, which leaves the reader with the impression that not enough effort has been made in the final polishing of the draft (with decreasing relevance).
a. Fig. 4(c) is squeezed compared to their counterparts (a) and (b).
b. Figs. 5 and 6 do not provide a colorbar legend.
c. There is an inconsistent way of referring to figures: “panels (ac)” vs. “panels a, b, c” vs. “panels (a) and (b)”
d. Obvious orthographic mistakes: “Feschbach”, “blushift”
e. Inconsistent use of hyphens: “FermiHubbard” vs. “Fermi Hubbard”, “ground state” vs. “groundstate”, “finitesize effect” vs. “finite size effect”, “Chevy ansatz” vs. “Chevyansatz”
f. Specific orthography: “p.e.” is not an official abbreviation, please use “e.g.” instead, in American English “labelled” is written as “labeled”.
g. Usually if a section is not specified with a numbering the word “section” is written in lower case.
Anonymous Report 1 on 20231012 (Invited Report)
 Cite as: Anonymous, Report on arXiv:2309.07019v1, delivered 20231012, doi: 10.21468/SciPost.Report.7933
Strengths
1 A wide variety of polaron problems is studied, including many problems that have not been studied in detail yet using other methods.
2 These problems are studied in detail and compared with previously known results, in general finding good qualitative agreement.
Weaknesses
1 The authors discuss that finite size effects are strong, but no attempt is made to quantify them. For instance, is there any form of convergence of the polaron spectra for 3x3, 4x4, 5x5 to the continuum result?
2Similarly, can some statement be made about how the discrete lines in Fig.2 converge to a moleculehole continuum?
3 Also, how do the results depend on the filling of the lattice?
Report
The paper studies in detail how polaron spectra can be studied using exact diagonalisation, finding (to my surprise), that even small systems suffice to get realiable estimates for the qualitative features of polaron spectra. Due to the relevance of polaron spectra to a wide variety of fields in condensed matter and ultracold atomic physics, and the many directions that can be explored from here, I can recommend this paper for publication after the points in weaknesses are at least attempted to be answered in more detail.
Requested changes
1 Answer the questions in "weaknesses"
(optional): Polaron's have been recently also probed using Rabi oscillation spectroscopy, see for example https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.125.133401 or https://arxiv.org/abs/2308.05746 or https://arxiv.org/abs/2308.06659. Rabi oscillation's are very challenging to calculate with field theory methods (due to conceptual difficulties), but should be easy to calculate in ED, so that could be an interesting application of this method.