{"id":"2691c281-5c7f-4cb6-b65c-9f08027ef640","arxiv_id":"2608.03853","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper derives a density-independent renormalization of Fermi polaron g factors, g1,2 = g_X ± ξg_c, and predicts a tunable crossover from linear to elliptical polarization when magnetic field and strain compete.","lead":"A theory paper derives how magnetic fields change the energy levels and polarization of Fermi polarons, quasiparticles formed when an exciton interacts with a Fermi sea of electrons, in atomically thin semiconductors. It explains why the magnetic response differs from simple trions and predicts a critical field where the electron sea becomes fully valley-polarized, leaving only one circularly polarized polaron state.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Orbital Landau quantization, neglected by the asserted ω_cτ_e≪1 criterion, undermines the Zeeman-only field dependence on which Eq. (41) and B_crit rest.","rationale":"I read the full manuscript. Eq. (41) follows algebraically from Eq. (39), given the constant-density expression for \\tilde E_F^τ, and the mass-ratio coefficient ξ is a parameter-free consequence of the single-pole self-energy. The reader's weakest-assumption identification matches mine: the Zeeman-only premise is the least secure condition for the central claim. I do not regard the orbital issue as an immediate refutation, because in the strict B→0 limit the linear Zeeman coefficient may be robust and because the author explicitly delimits the regime; but for the experimentally relevant high-mobility samples used in the literature, ω_cτ_e ≫ 1 at fields as low as 0.1 T, so spectra near B_crit and the predicted kink would be modified. I also weighed the two self-reported limitations: the retained beyond-accuracy term in Eq. (24) and the one failed Stokes state in Sec. III B. Both are real but target the strained two-level polarization model, not the unstrained g-factor formula Eq. (41). The proposed Landau-level calculation would settle whether the central claim survives beyond the Zeeman-only assumption. Since the reader already issued CONDITIONAL on essentially this concern, my verdict is UNCHANGED.","tokens_in":21893,"tokens_out":15726,"duration_ms":142913,"concrete_test":"Recompute the attractive Fermi polaron energies from Eq. (10) for the Fig. 2(c,d) parameters, replacing the free-electron density of states in S^{στ}(q), Eq. (13), with the Landau-quantized density of states (with a small phenomenological broadening Γ ≪ ℏω_c) while keeping the same Zeeman shifts. Extract (i) the slope of the intravalley splitting E_{Rr} − E_{Ll} at B → 0 and (ii) the field at which the minority valley becomes empty. If (i) differs from (g_X + ξg_c)μ_B by more than 10% or (ii) shifts from 8 T by more than the typical linewidth 0.2 meV, the Zeeman-only premise is load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central magnetic-field dependence rests on the Zeeman-only premise in Sec. II A, justified by the inequality |eB/(M_e c)|τ_e << 1. For hBN-encapsulated TMDC monolayers, mobilities of 10^3–10^4 cm^2/Vs imply τ_e ~ 0.1–2 ps, so at B = 1 T the cyclotron frequency satisfies ω_c τ_e ~ 10^2–10^3; the criterion is violated in the very samples where Fermi-polaron effects are measured. In the paper's own parameter set (N_e = 8×10^10 cm^-2, m* = 0.4m0, Fig. 2(c,d) and Fig. 4), B_crit = 8 T, E_F ≈ 0.48 meV, while ℏω_c ≈ 2.3 meV; the Landau level spacing exceeds the Fermi energy, so the electron-hole-pair continuum in S^{στ}(q), Eq. (13), and the valley-filling formulas (35)–(37) are not described by a constant 2D density of states. The derivation of \\tilde E_F^τ = E_F ∓ g_c μ_B B/2 and of B_crit = 2E_F/(|g_c|μ_B) is tied to that constant density of states. If Landau quantization is included, the self-energy integral and the population threshold change; the linear g-factor correction may survive in the strict B→0 derivative, but the claimed field dependence, the B_crit kink, and the density-independent renormalization at finite small fields are not controlled. The paper does not estimate τ_e or test its criterion against literature mobility values. Secondary internal flags (the retained beyond-accuracy term in Eq. (24) and the admitted failure of one of the four Stokes states in Sec. III B) concern the strained two-level model rather than Eq. (41) itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Green's-function theory of attractive Fermi polarons (Suris tetrons) in doped TMDC monolayers in a perpendicular magnetic field, with and without uniaxial strain. Starting from a band model with valley-dependent Zeeman shifts, it derives the exciton self-energy, obtains the Fermi polaron energy in the single-pole approximation, and predicts a correlation-induced renormalization of the polaron g-factor: Eqs. (41b,c) give g1 = g_X + ξ g_c for intravalley and g2 = g_X − ξ g_c for intervalley attractive polarons, with ξ a mass-ratio coefficient. It also predicts a critical magnetic field B_crit = 2E_F/(|g_c| μ_B) for full valley polarization, a quadratic-to-linear Zeeman crossover in strained samples, and circular/linear dichroism in absorption and reflection. The central derivation is presented in closed form, and the simplified two-level model is compared with the full Green's-function numerics in Fig. 4.","tokens_in":22178,"tokens_out":6547,"duration_ms":62722,"significance":"If the results hold, the paper provides parameter-free, experimentally testable predictions: the polaron g-factor differs from the bare trion g-factor by a calculable correlation correction of either sign, and the difference is independent of density in the small-field regime. This goes beyond the bare-trin model and connects to known experiments on inter- and intravalley trion g-factors. The derivation uses material inputs g_X, g_c and masses from prior measurements, not the target result, and the analytic formulas are a genuine strength. The finite-field predictions, especially B_crit and the field dependence in Figs. 2-5, are nevertheless tied to a Zeeman-only, constant-density-of-states assumption that is not quantitatively justified in the regime where the predictions are made.","major_comments":[{"comment":"The Zeeman-only premise is asserted with the criterion |eB/(M_e c)| τ_e << 1, but no estimate of τ_e is given and the criterion is not tested against the parameters used in Figs. 2-5. For hBN-encapsulated TMDC monolayers with mobilities of 10^3-10^4 cm^2/Vs, τ_e ~ 0.1-2 ps, so at B = 1 T one has ω_c τ_e ~ 0.1-1 and at B_crit = 8 T, ω_c τ_e ~ 1-10; the inequality is not satisfied. In the paper's own example (N_e = 8×10^10 cm^-2, m* = 0.4m0), E_F ≈ 0.48 meV while ℏω_c ≈ 3.7 meV at B = 8 T, so the constant 2D density of states used for \\tilde E_F^τ and B_crit is not controlled. This affects the valley-filling formulas (35)-(37), the self-energy integral in Eq. (13), and the finite-field splittings and dichroism shown in Figs. 2-5. Only the strict B→0 derivative leading to Eq. (41) is protected from this criticism. Please add a quantitative validity estimate for τ_e and either restrict the finite-field claims to the regime where the criterion holds or incorporate Landau quantization.","section":"Sec. II A, Eqs. (4)-(6) and Sec. III A 2, Eqs. (35)-(37)"},{"comment":"The manuscript states that the additional term ∝ E_F^2 in Eq. (24) 'exceeds the accuracy' of the derivation but is retained for better agreement with numerics. Because Eq. (24) enters the simplified two-level splitting Eq. (54), which is compared with the full numerical solution of Eq. (55) in Fig. 4(a), the agreement shown cannot be interpreted as a clean validation of the analytic model; it may be produced partly by an uncontrolled term. The authors should either remove the beyond-accuracy term and repeat the comparison, or provide a controlled derivation that justifies it.","section":"Eq. (24) and the text immediately after it"},{"comment":"The text admits that 'three out of four states are in reasonable agreement with the model, while one state remains mostly linearly polarized.' Since the simplified two-level Hamiltonian of Eq. (52) is used as the paper's main analytical tool for the strain-Zeeman interplay, this is a clear limitation of the model, not just a numerical detail. The manuscript should state explicitly the regime of validity of Eq. (52) and explain why one of the four states is not captured, rather than leaving this as an unexplained exception in the figure caption.","section":"Sec. III B, Fig. 4(b), Eqs. (52)-(56)"}],"minor_comments":[{"comment":"The statement that the g-factor renormalization 'appears even at E_F → 0' is imprecise because B_crit ∝ E_F: in the limit E_F → 0, the regime |B| < B_crit shrinks to B = 0. The intended meaning is presumably the low-density limit at fixed small field; please rephrase.","section":"Sec. III A 2, text below Eq. (41)"},{"comment":"The quantity S_s is used in Eq. (17) before its definition in Eq. (18); defining S_s before Eq. (17) would improve readability.","section":"Sec. II B, Eq. (17)"},{"comment":"There is a typo in the conclusion: 'Fermi plaron states' should be 'Fermi polaron states'.","section":"Conclusion, Sec. V"},{"comment":"The term 'pseudo-Zeeman splitting' is used for the strain-induced splitting; since the same term appears in the literature for other effects, a brief clarifying definition at first use would help.","section":"Abstract and Sec. III B"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for cond-mat.mes-hall. The g-factor renormalization formula is new and analytically clean, and the derivation is largely self-consistent. The main risk is the orbital-field approximation: the finite-field predictions and B_crit rest on a criterion that appears to fail in high-mobility TMDC samples at the fields shown in the paper's own figures. I would ask the authors to supply a quantitative mobility/time-scale estimate and, if necessary, restrict the claims to the regime where the Zeeman-only approximation is controlled."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper's central result is Eq. (41): attractive Fermi polarons get g factors g1 = g_X + ξ g_c and g2 = g_X - ξ g_c for intra- and intervalley complexes, with ξ a mass-ratio coefficient. That is a parameter-free, physically transparent outcome of the Fermi polaron (Suris tetron) model, and it gives a natural explanation for the experimentally observed differences between intra- and intervalley trion g factors. The derivation is coherent: the self-energy, the trion-pole approximation, and the linear-in-density expansion are consistent, and the numerics match the two-level model well. The strain–magnetic-field crossover and the B_crit kink are nice concrete predictions.\n\nThe main soft spot is the neglect of orbital magnetic effects. The paper justifies it with ω_c τ_e << 1, but does not estimate τ_e. In high-mobility hBN-encapsulated samples, τ_e can be a few ps, and at the B_crit of the paper's own parameter set (8 T, E_F ~ 0.5 meV), ℏω_c is ~2 meV, so the continuum density of states used for the valley-filling formulas is not controlled. The stress-test's claimed ω_c τ ~ 10^2–10^3 seems too aggressive by an order of magnitude or more, but the qualitative point stands: the finite-field predictions (B_crit, the quadratic-to-linear crossover) are conditional on the continuum approximation, while the strict B→0 g-factor renormalization should survive. I would like the paper to state the regime more carefully and ideally compare with a Landau-level treatment in the clean limit.\n\nMinor notes: Eq. (24) retains a term beyond the accuracy of the derivation for better numerical agreement; the author says so, but a referee should ask if that fudge leaks into the quoted formulas. The simplified two-level model also fails to describe one of the four Stokes states, as the paper admits; that limits the intuitive model but not the core g-factor result. Citation practice looks fair; the strain-only predecessor is the author's own published benchmark.\n\nWho is this for? People working on magneto-optics of doped TMDC monolayers and trion g factors. It is a solid theory paper with a clear, useful result and a plausible explanation of data. It deserves a serious referee. My recommendation: engage with it, but demand a careful discussion of the orbital-effect validity range and a clean separation between the robust linear-response result and the model-dependent high-field predictions.","headline":"Fermi polaron g-factor renormalization (g1,2 = g_X ± ξ g_c) is a clean, parameter-free result that likely explains measured intra-/intervalley trion g-factor differences; the high-field predictions are shakier because orbital Landau effects are neglected, but the paper is worth a serious referee.","tokens_in":22805,"tokens_out":10170,"would_cite":true,"duration_ms":78712,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The Zeeman splitting of attractive Fermi polarons in doped monolayer semiconductors is set by two renormalized g factors, g_X plus or minus a mass-ratio term, not by the bare trion g factor.","keywords":["Fermi polaron","Suris tetron","transition metal dichalcogenide monolayers","Zeeman effect","g-factor renormalization","uniaxial strain","valley polarization","circular dichroism"],"falsifier":"Measure the energy and circular polarization of the two attractive Fermi polaron absorption lines in a gate-tunable WSe$_2$ monolayer as a function of $B$ at fixed electron density below $B_{\\rm crit}$; the central prediction fails if either line splits with the bare trion slope $g_X\\mu_B$, if both circularly polarized lines remain above $B_{\\rm crit}$, or if the small-field Zeeman slopes change with electron density.","tokens_in":21570,"feed_emoji":"🧲","tokens_out":12357,"duration_ms":103091,"temperature":0.7,"pith_summary":"The paper predicts that in doped monolayer semiconductors, the Zeeman splitting of attractive Fermi polarons is not the bare trion splitting: correlations with the valley-resolved Fermi sea renormalize the $g$ factor by a mass-ratio coefficient $\\xi$, giving $g_1=g_X+\\xi g_c$ for intravalley and $g_2=g_X-\\xi g_c$ for intervalley polarons in W-based materials. This renormalization is density-independent at small fields, so it survives in the $E_F\\to0$ limit, and it changes sign between the two polaron branches. A critical field $B_{\\rm crit}=2E_F/(|g_c|\\mu_B)$ fully valley-polarizes the electron sea, switching the radiative doublet to a single circularly polarized state. With uniaxial strain added, the Zeeman and strain-induced pseudo-Zeeman splittings compete, producing a quadratic-to-linear Zeeman crossover and continuously tunable elliptical polarization. The theory also computes absorption and reflection spectra, including circular and linear dichroism, giving concrete experimental signatures.","feed_headline":"Polarons in doped monolayers get a Fermi-sea g-factor shift","feed_subtitle":"Doped-monolayer polarons should show Zeeman slopes shifted from bare trions by a mass-ratio term.","key_machinery":"The load-bearing object is the attractive Fermi polaron, a four-particle quasiparticle composed of an exciton plus a Fermi-sea electron-hole pair, described by the Hamiltonian in Eq. (6) and the wavefunction in Eq. (8). Its energy is fixed by the self-consistent equation $E^\\sigma=\\varepsilon^\\sigma_0+\\Sigma^\\sigma(E^\\sigma)$, where the exciton self-energy $\\Sigma$ sums the interaction with electrons in each valley; in the trion-pole approximation this reduces to the compact form in Eq. (45) with the coefficients $\\alpha$, $\\beta$, and the mass-ratio factor $\\xi$ of Eq. (40), which carries the entire $g$-factor renormalization. For strain, the central object is the tensor Green's function $\\hat G_s(E)$ of Eq. (50) and the effective two-level Hamiltonian $H_{1,2}=(-E_{T1,2}+\\xi E_F)\\hat 1+\\frac{\\Delta E^{(s)}_{FP}}{2}\\hat\\sigma_x+\\frac{\\Delta E^{W}_{FP1,2}}{2}\\hat\\sigma_z$, which treats the magnetic field as a $\\sigma_z$ term and the strain as a $\\sigma_x$ pseudomagnetic term on the polaron pseudospin. This Hamiltonian yields the analytic splittings, Stokes parameters, and crossover fields used throughout the paper.","core_discovery":"The central result is Eq. (41): in the small-field regime $|B|<B_{\\rm crit}$, the Zeeman splittings of the two attractive Fermi polaron doublets in W-based monolayers are $\\Delta E^{W}_{FP1,2}=g_{1,2}\\mu_B B$ with $g_1=g_X+\\xi g_c$ and $g_2=g_X-\\xi g_c$, where $\\xi$ is a negative coefficient fixed by the electron-to-exciton mass ratio through Eq. (40). Correlations between the exciton and the Fermi sea thus shift the polaron $g$ factor away from the trion value $g_X$, with opposite signs for the intra- and intervalley branches, and the shift is independent of electron density at small fields. At $B_{\\rm crit}=2E_F/(|g_c|\\mu_B)$ the resident electrons become fully valley polarized, so only one circularly polarized polaron state remains and the density dependence of the splitting has a kink. Under uniaxial strain, the same doublets first shift quadratically in $B$ while remaining nearly linearly polarized, then cross over to linear Zeeman splitting with elliptical polarization. The paper derives these behaviors from the Green's-function self-energy approach and a simplified two-level effective Hamiltonian, and works out their signatures in absorption, reflection, and both circular and linear dichroism.","pith_inferences":["A natural extension is to look for the same mass-ratio factor $\\xi$ in the magnetic-field dependence of the Fermi polaron spectral continuum and oscillator-strength redistribution, since the self-energy that produces Eq. (41) also controls those quantities.","The two-level $\\sigma_x$-versus-$\\sigma_z$ picture implies that sweeping $B$ and strain moves the polaron doublet along a full pseudospin trajectory, so the predicted elliptical polarization could be used as a continuously tunable polarization rotor; the paper itself does not propose this application.","The kink in the density dependence at $B_{\\rm crit}$ could be turned into a spectroscopic measurement of the conduction-band $g_c$ that is independent of transport measurements, by tracking the threshold as electron density is tuned.","If orbital (Landau) effects become visible before $B_{\\rm crit}$ in a given sample, the one-state regime would acquire additional structure; the Zeeman-only prediction should then be tested at fields below the first Landau gap or in samples with short momentum relaxation."],"forward_implications":["In W-based monolayers, the intravalley and intervalley attractive polaron lines should show Zeeman slopes $g_X+\\xi g_c$ and $g_X-\\xi g_c$ rather than a single trion slope, a direct magneto-optical test.","Because the renormalized $g$ factor is density-independent at small fields (and survives $E_F\\to0$), lowering the electron density should not push the slopes back toward the bare trion value as long as $B<B_{\\rm crit}$.","Above $B_{\\rm crit}$, only one circularly polarized polaron state remains in each radiative doublet, so the number of observed transitions drops by half and the splitting versus density shows a kink at the threshold.","Under uniaxial strain, each doublet first exhibits a quadratic Zeeman shift while staying nearly linearly polarized, then crosses to linear Zeeman splitting with elliptical polarization; the crossover field is set by the relative size of strain splitting and $g\\mu_B B$.","Mo-based monolayers have only an intervalley polaron, with $g_{\\rm Mo}=g_X+\\xi g_c$, providing a cleaner single-line test of the same renormalization."],"supporting_citations":[{"why":"supplies the Fermi polaron model for charged excitons in two-dimensional semiconductors and the self-energy formalism used to derive the Zeeman splittings.","marker":"[39]"},{"why":"introduces the exciton plus Fermi-sea electron-hole pair (tetron) wavefunction on which the calculation is built.","marker":"[60]"},{"why":"establishes that trion absorption features must include correlations with the Fermi sea, motivating the deviation from the bare trion g factor.","marker":"[34]"},{"why":"extends many-body exciton-polaron theory to strong magnetic fields, providing the context for valley polarization and high-field behavior.","marker":"[35]"},{"why":"reports experimentally distinct g factors for intra- and intervalley trions in WSe2, the observation the renormalized polaron g factors explain.","marker":"[37]"},{"why":"previous theory of the strain-induced fine structure of Fermi polarons that the strain part of the present model extends to magnetic fields.","marker":"[53]"},{"why":"experimental observation of Fermi polarons under strain-induced pseudomagnetic fields, the regime the paper's magneto-strain predictions address.","marker":"[54]"},{"why":"provides the strain-induced exciton fine-structure Hamiltonian used for the pseudo-Zeeman term.","marker":"[51]"},{"why":"fixes the band g-factor conventions for g_c, g_v, and g_X in the Zeeman shifts of Eq. (4).","marker":"[62]"}],"fun_headline_variants":["Polaron g-factor shifts with Fermi-sea correlations in monolayers","Magnetic-field control of polaron doublets in strained monolayers","Fermi sea shifts polaron Zeeman splitting in doped monolayers","Critical field flips polaron polarization in strained monolayers","Strain and magnetic field tune polaron fine structure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All magnetic-field effects are treated as Zeeman shifts of the bands, with orbital (Landau) effects dropped under the condition $|eB/(M_e c)|\\tau_e\\ll1$; in high-mobility monolayers at low temperature this condition can break at fields of a few tesla, and then the renormalized $g$ factors and the valley-polarization threshold derived here would no longer apply.","fun_headline_variants_meta":{"raw":{"variants":["Polaron g-factor shifts with Fermi-sea correlations in monolayers","Magnetic-field control of polaron doublets in strained monolayers","Fermi sea shifts polaron Zeeman splitting in doped monolayers","Critical field flips polaron polarization in strained monolayers","Strain and magnetic field tune polaron fine structure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1393,"prompt_tokens":1054,"completion_tokens":339,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":252}},"tokens_in":670,"tokens_out":339,"duration_ms":3525,"temperature":1.0,"reasoning_tokens":252,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:45:25.808231+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the energy and circular polarization of the two attractive Fermi polaron absorption lines in a gate-tunable WSe$_2$ monolayer as a function of $B$ at fixed electron density below $B_{\\rm crit}$; the central prediction fails if either line splits with the bare trion slope $g_X\\mu_B$, if both circularly polarized lines remain above $B_{\\rm crit}$, or if the small-field Zeeman slopes change with electron density.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the Fermi polaron model for charged excitons in two-dimensional semiconductors and the self-energy formalism used to derive the Zeeman splittings."},{"cited_title":"Yagodkin, K","cited_arxiv_id":null,"evidence_quote":"introduces the exciton plus Fermi-sea electron-hole pair (tetron) wavefunction on which the calculation is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"establishes that trion absorption features must include correlations with the Fermi sea, motivating the deviation from the bare trion g factor."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"extends many-body exciton-polaron theory to strong magnetic fields, providing the context for valley polarization and high-field behavior."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"experimental observation of Fermi polarons under strain-induced pseudomagnetic fields, the regime the paper's magneto-strain predictions address."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"fixes the band g-factor conventions for g_c, g_v, and g_X in the Zeeman shifts of Eq. (4)."}],"review_version":2}