{"id":"37782222-bd1d-4042-8e7e-54be3850714d","arxiv_id":"1908.08789","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Uniaxial strain in ZrSiS preserves its flat infrared optical conductivity under compression, triggers a Lifshitz transition under ~2 GPa tension that weakens screening, and tunes the dispersion of ~20 eV low-loss hyperbolic plasmons.","lead":"This paper uses first-principles calculations to predict how squeezing or stretching the nodal-line semimetal ZrSiS changes its optical properties. It finds that the material's infrared response stays flat under compression but loses spectral weight under tension, and that its high-energy plasmons, including a predicted deep-UV hyperbolic regime, can be tuned by strain.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The GGA-PBE placement of the quadratic band is unvalidated at the meV scale, so the predicted ~2 GPa Lifshitz transition and the abstract's stated critical stress are quantitatively insecure.","rationale":"The reader's weakest assumption explicitly identifies the position of the quadratic band near the Fermi level as fragile, and I agree that this is the most load-bearing premise. The abstract's headline tensile prediction (Lifshitz transition at ~2 GPa weakening interband screening) is a quantitative Fermi-surface claim, and the paper's own detailed numbers (1.3 and 3.4 GPa) already disagree with the '2 GPa' used in the abstract and conclusions. That internal mismatch signals that the strain axis is not carefully pinned. The underlying GGA-PBE band structure is calibrated only against low-energy aggregate optical properties, which do not constrain a single quadratic band position to the needed tens-of-meV accuracy, and the neglected SOC is not obviously irrelevant at a Fermi-surface topology change. A single computational check with HSE06 or G0W0 plus SOC at a few tensile strains would settle whether the transition strain survives. I am not rejecting the paper: the low-energy zero-strain results are well reproduced, and the high-energy hyperbolic prediction is also uncertain but less central to the strain-tuning claim. The reader's CONDITIONAL verdict is appropriate, so no verdict change is recommended.","tokens_in":15155,"tokens_out":11057,"duration_ms":116122,"concrete_test":"Repeat the DFT band-structure calculation at tensile strains uzz = 0, -1%, -2%, -3%, and -4% using a more accurate electronic-structure method (HSE06 hybrid functional or G0W0) with spin-orbit coupling included, and determine the strain at which the quadratic band along Z–R crosses epsilon_F and at which the hole pockets merge or separate. If the critical strain moves by more than about 1% (approximately 1 GPa) relative to the paper's P1 and P2 values, or if no Fermi-surface topology change occurs below 5% tension, then the abstract's '2 GPa Lifshitz transition' and the associated infrared spectral-weight reduction are not robust. If such a calculation is too expensive, an alternative check is to compare the unstrained GGA-PBE band position of the quadratic band against high-resolution ARPES data, requiring agreement to within ~20 meV.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central tensile-strain prediction in Sec. IV and the abstract is a Lifshitz transition at ~2 GPa that weakens interband screening and reduces infrared spectral weight. That prediction rests on the position of a quadratic electron band along Z–R relative to the Fermi level. The paper uses GGA-PBE without spin-orbit coupling (Sec. II A) and calibrates only low-energy aggregate quantities: the unscreened plasma frequency, the flat optical conductivity, and epsilon_inf (Sec. III A). None of these is sensitive to a tens-of-meV placement of this quadratic band. The text is internally inconsistent about the critical stress: the abstract and conclusions say 2 GPa, while Sec. IV and Fig. 8 report two transitions at P1 = 1.3 GPa and P2 = 3.4 GPa. Because a Lifshitz transition is a Fermi-surface topology change at exactly epsilon_F, an error of a few tens of meV—comparable to the SOC gaps the authors neglect—can shift the critical strain by more than 1% or remove the transition entirely. The stated justification for neglecting SOC (low temperatures, <20 meV) is not convincing for a zero-energy Fermi-surface property. Thus the quantitative strain-tuning claim is currently unsupported without a band-structure check at the relevant energy scale.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles calculations of the optical conductivity, dielectric function, and plasmonic response of the nodal-line semimetal ZrSiS, focusing on the effect of uniaxial strain. Using GGA-PBE band structures and a scalar RPA dielectric function with a Drude correction (Sec. II), the authors reproduce several low-energy experimental quantities for the pristine compound: the in-plane unscreened plasma frequency (3.15 vs 2.88 eV), the flat infrared conductivity (~7000 vs 6600 Ohm^-1 cm^-1), and epsilon_inf ~ 9 vs the experimental ~7.8. The main claims are that the frequency-independent infrared conductivity is robust under uniaxial compression up to 10 GPa but its flat region narrows with increasing compression; that tensile stress around 2 GPa induces a Lifshitz transition that weakens interband screening and reduces infrared spectral weight; and that the high-energy response hosts low-loss ~20 eV plasmons with strongly anisotropic dispersion, including a possible hyperbolic regime in the deep ultraviolet. The strain dependence of the high-energy plasmon dispersion is also discussed.","tokens_in":15307,"tokens_out":5365,"duration_ms":54513,"significance":"If the predictions hold, the paper identifies ZrSiS as a mechanically tunable natural hyperbolic material for ultraviolet plasmonics and connects the strain response to the chemical-pressure trend across the ZrSiX family (X=S, Se, Te). The manuscript has clear strengths: the zero-strain calculations are calibrated against measured optical quantities, the Drude contribution is handled explicitly, and the central predictions are concrete and falsifiable by EELS and infrared spectroscopy under uniaxial stress. The high-energy and tensile-strain predictions, however, rest on approximations that are not validated at the relevant energy scale, which is the main source of uncertainty in the paper.","major_comments":[{"comment":"The abstract and Conclusions state that 'Upon uniaxial tensile stress of 2 GPa, the Fermi surface undergoes a Lifshitz transition,' but Sec. IV identifies two Lifshitz transitions at P1 ≈ 1.3 GPa and P2 ≈ 3.4 GPa (with tension as negative stress in Fig. 8). The manuscript should reconcile this discrepancy: if the headline claim refers to the P1 transition, the abstract should say so and mention the second transition; if it refers to a single transition near 2 GPa, that conflicts with the two computed critical stresses. As written, the central quantitative claim is internally inconsistent.","section":"Abstract; Sec. IV, Fig. 8"},{"comment":"The prediction that tensile strain of order 2 GPa triggers a Lifshitz transition and weakens interband screening depends on the energy position of the quadratic electron band along Z–R relative to the Fermi level. This band is obtained from GGA-PBE without spin-orbit coupling, and the calibration quantities in Sec. III A—unscreened plasma frequency, flat optical conductivity, and epsilon_inf—are aggregate low-energy properties that are insensitive to the band position at the meV scale. Because the stated SOC-induced gaps are up to ~20–30 meV, an error of that size could shift the critical strain substantially or remove the transition altogether. The authors should provide a direct check of this band's position (e.g., ARPES comparison, SOC-included band structure, or a hybrid-functional calculation) or explicitly qualify the predicted critical stress as an estimate with large uncertainty.","section":"Sec. IV, Figs. 7-8"},{"comment":"The claims of low-loss ~20 eV plasmons and a hyperbolic regime are based on the signs of epsilon_xx(omega) and epsilon_zz(omega) computed in the scalar RPA with local field effects neglected and using GGA-PBE band structures. The experimental calibration in Sec. III A covers only the low-energy (<= 2 eV) response; no independent check is presented for the ~20 eV region, where the permittivity signs are the load-bearing input for the hyperbolic dispersion. The authors should discuss the expected influence of local field effects, self-energy corrections, and SOC on the 20 eV permittivity, or compare with existing EELS or vacuum-UV optical data. Alternatively, the hyperbolic-regime claim should be explicitly labeled as a prediction that remains to be confirmed.","section":"Sec. III B, Fig. 5"}],"minor_comments":[{"comment":"The strain sign convention is confusing: the text defines uzz as positive for compression and negative for tension, while the critical stresses are printed as P1 = -1.3 GPa and P2 = -3.4 GPa. Please adopt a single, clearly stated convention and apply it consistently in the text and figure captions.","section":"Sec. IV; Figs. 7 and 8"},{"comment":"The expression for sigma_{2,intra}(omega) appears to be missing a factor of omega; the standard relation is sigma_2(omega) = omega [1 - epsilon_1(omega)]/(4 pi). Please check and correct the formula.","section":"Sec. II B, Eq. (8)"},{"comment":"There are several typographical errors in this section: 'Simliar' should be 'Similar', 'hyprobolic' should be 'hyperbolic' (twice), and 'electronmagnetic' should be 'electromagnetic'.","section":"Sec. III B"},{"comment":"In the discussion of the nonsymmorphic Dirac node, 'Andreas et al.' should be 'Topp et al.' (Ref. [13]) for accuracy.","section":"Sec. IV"},{"comment":"The values of the dispersion coefficient A in Eq. (10) are not reported; given the paper's emphasis on strain tuning of the plasmon dispersion, providing A and its strain dependence would make the claim quantitative.","section":"Sec. IV, Eq. (10) and Fig. 10"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a careful first-principles study that reproduces the measured zero-strain optics of ZrSiS and then pushes two new claims — a strain-driven Lifshitz transition that weakens interband screening, and a ~20 eV hyperbolic plasmon regime. The low-energy calibration is genuinely good: unscreened plasma frequency 3.15 vs 2.88 eV, flat conductivity ~7000 vs 6600, ε∞≈9 vs 7.8. That anchors the method where it matters for the infrared claim. And the results are computed, not fitted; the only tuned parameter is broadening η, and they test its effect. No circularity.\n\nWhat's actually new: the strain dependence of the flat conductivity region, the tensile Lifshitz transition, and the deep-UV hyperbolic plasmons. Prior work measured zero-strain optics and connected the ZrSiX family via chemical pressure; this adds the strain axis.\n\nThe soft spots are where the paper makes its loudest claims. The hyperbolic regime at ~20 eV rests on scalar RPA without local field effects, computed from GGA-PBE bands without SOC, and there is no independent check in that energy range. That doesn't kill it, but it means the UV prediction is a plausibility argument, not a settled result.\n\nThe bigger problem is the Lifshitz transition. The abstract and conclusions say '2 GPa', but Sec. IV and Fig. 8 report P1=1.3 GPa and P2=3.4 GPa. Those are two distinct topology changes; the text does not say which one causes the screening reduction. This is not a typo-level issue, because the whole tensile-strain story hangs on the quadratic band crossing the Fermi level at a specific stress, and that crossing is exactly the kind of meV-scale feature that GGA-PBE with no SOC can misplace. The authors justify neglecting SOC as a low-temperature/low-frequency effect, but a Lifshitz transition is a zero-energy Fermi-surface property, so that justification doesn't hold. A tens-of-meV error could shift the critical stress by more than 1% or remove the transition. They calibrate aggregate low-energy quantities, but none of those is sensitive to the quadratic band's position within tens of meV.\n\nWho should read this: anyone working on ZrSiX optics, topological semimetal plasmons, or natural hyperbolic materials. The infrared part is solid enough to cite; the UV and Lifshitz parts should be treated as predictions needing follow-up.\n\nRecommendation: send it to peer review. A serious referee should ask for (a) a consistent statement of the critical stress, (b) a band-structure check of the quadratic band position at the relevant energy scale (e.g., hybrid functional or GW, and an estimate of SOC-induced shifts), and (c) a caveat on the scalar-RPA/no-LFE approximation in the UV. With those addressed, it would be a useful paper.","headline":"A solid, well-calibrated DFT+RPA study with interesting strain predictions, but the headline Lifshitz-transition stress is internally inconsistent and rests on meV-scale band physics the method doesn't anchor.","tokens_in":15934,"tokens_out":2635,"would_cite":true,"duration_ms":25936,"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":"ZrSiS, a layered nodal-line semimetal, keeps its flat infrared conductivity under compression up to 10 GPa, undergoes a Lifshitz transition near 2 GPa tension that lowers infrared spectral weight, and hosts low-loss ~20 eV plasmons whose…","keywords":["nodal-line semimetal","ZrSiS","optical conductivity","plasmon","hyperbolic material","strain engineering","Lifshitz transition","first-principles calculation"],"falsifier":"Measure the electron energy loss spectrum of a ZrSiS crystal under uniaxial strain: the paper predicts a sharp low-loss peak near 20 eV with anisotropic quadratic dispersion, and a separate check is that the in-plane and out-of-plane permittivities obtained from optical ellipsometry should have opposite signs in a ~0.6 eV window near 20 eV; if the 20 eV loss peak is absent or the permittivity signs are the same, the high-energy plasmon and hyperbolic claims fail. For the Lifshitz transition, quantum oscillation or ARPES measurements on a crystal under ~2 GPa tensile stress should show the hole pockets merging near $k_z=\\pi/c$ and the quadratic band crossing the Fermi level along Z-R.","tokens_in":14856,"feed_emoji":"🔬","tokens_out":7030,"duration_ms":65144,"temperature":0.7,"pith_summary":"This paper uses first-principles calculations to ask how external strain changes the optical response of the nodal-line semimetal ZrSiS. It predicts that the material's signature flat infrared conductivity survives uniaxial compression up to 10 GPa, but the flat region narrows as the load increases. Under tensile stress near 2 GPa, a Lifshitz transition reshapes the Fermi surface and weakens interband screening, cutting the infrared spectral weight. At high energy the same calculations predict low-loss plasmons near 20 eV with anisotropic dispersion, and the strongly anisotropic dielectric response implies a type-I hyperbolic regime for deep-ultraviolet plasmons. The practical payoff is a single air-stable layered material whose infrared and ultraviolet optical properties could be tuned mechanically.","feed_headline":"ZrSiS survives 10 GPa and gains tunable UV plasmons","feed_subtitle":"First-principles study predicts a 2 GPa Lifshitz switch and low-loss hyperbolic plasmons near 20 eV.","key_machinery":"The argument is carried by the wave-vector- and frequency-dependent dielectric function computed in the random phase approximation from first-principles band structures, with interband transitions treated from dipole matrix elements in the long-wavelength limit and intraband transitions added as a standard Drude term. Plasma excitations are identified with the zeros of the real part of the dielectric function, and the energy loss function $\\operatorname{Im}[-1/\\epsilon(q,\\omega)]$ locates which of these modes are low-loss. The hyperbolic regime is diagnosed from the product $\\epsilon_{xx}(\\omega)\\,\\epsilon_{zz}(\\omega) < 0$, which makes the constant-frequency surface in Eq. (11) a hyperboloid rather than a sphere. The same dielectric function, evaluated at $q \\to 0$ with a finite damping parameter, produces the infrared conductivity and the screened plasma frequency that calibrate the calculation against experiment.","core_discovery":"The central claim is that ZrSiS combines a strain-tolerant low-energy optical fingerprint with strain-tunable high-energy plasmonics. The frequency-independent in-plane optical conductivity, a known fingerprint of this nodal-line semimetal, remains flat under uniaxial compression up to 10 GPa, although the flat window shrinks and the spectral weight grows by roughly 50% at 5% compression. Tensile strain acts differently: near 1.3-3.4 GPa the Fermi surface changes topology through two Lifshitz transitions, and the resulting drop in interband screening reduces the infrared spectral weight and raises the screened plasma frequency from about 1.0 eV to about 1.3 eV at 4% tension. In the high-energy region, the calculations show a weakly damped plasmon near 20 eV whose quadratic dispersion is strongly anisotropic, and where the product of the in-plane and out-of-plane permittivities is negative over a ~0.6 eV window, the defining condition for type-I hyperbolic plasmons. The strain response of this high-energy mode is small in frequency but large in dispersion, so strain acts as a tuning knob rather than a switch.","pith_inferences":["Beyond the paper's own claims, the predicted tensile Lifshitz transition offers a clean test of the chemical-pressure picture: stretched ZrSiS should optically resemble ZrSiSe or ZrSiTe, and the ~2 GPa transition stress could be measured by quantum oscillations or ARPES under a four-point bending setup.","Beyond the paper's own claims, the 20 eV hyperbolic window, if verified by EELS, would make ZrSiS one of the few natural (unstructured) hyperbolic materials in the deep ultraviolet; applications like subwavelength imaging and thermal emission engineering would follow without nanofabrication.","Beyond the paper's own claims, since the scalar RPA neglects local-field effects and spin-orbit coupling, the most decisive numerical check is a full dielectric-matrix or Bethe-Salpeter calculation at 15-25 eV; if the sign of either permittivity changes there, the hyperbolic regime would move or disappear."],"forward_implications":["Uniaxial compression up to 10 GPa leaves the flat infrared conductivity intact, so ZrSiS-based infrared elements could tolerate large mechanical loads without losing their broadband response.","Tensile stress near 2 GPa acts as a mechanical switch: the Lifshitz transition weakens interband screening and lowers the infrared spectral weight, with the effect saturating after the Fermi-surface reconstruction around 3.4 GPa.","Electron energy loss measurements should find a sharp, low-loss peak near 20 eV whose dispersion is quadratic and markedly different in-plane and out-of-plane.","The negative product of permittivities near 5 and 20 eV implies type-I hyperbolic plasmon propagation in windows of about 0.6 eV; the 20 eV window is the practically relevant one because the 5 eV mode is strongly damped.","Strain of up to 5% changes the 20 eV plasmon frequency by only a few percent, but can change its out-of-plane dispersion coefficient by roughly 30%, giving a quantitative strain-tuning handle."],"supporting_citations":[{"why":"Supplies the experimental infrared conductivity, flat-region value, and screening constant used to calibrate the low-energy calculation.","marker":"[33]"},{"why":"Explains the flat optical conductivity as a combination of intraband and interband transitions, the mechanism the paper analyzes under strain.","marker":"[44]"},{"why":"Establishes the optical properties of the ZrSiX family and the chemical-pressure trend with c/a, the comparison for strained ZrSiS.","marker":"[45]"},{"why":"Supplies the power-law scaling for interband conductivity that motivates the frequency-independent behavior in Dirac-like systems.","marker":"[39]"},{"why":"Identifies ZrSiS as a nodal-line semimetal with two kinds of Dirac nodal lines, the band-structure context for the optical analysis.","marker":"[12]"},{"why":"Reports the lower screening constant of ZrSiTe, which supports the chemical-pressure interpretation of stretched ZrSiS.","marker":"[76]"},{"why":"Motivates the strain study by showing a possible topological phase transition in ZrSiS under high pressure.","marker":"[36]"},{"why":"Describes a four-point bending setup that could apply the tensile stresses the paper predicts to access the Lifshitz transition.","marker":"[79]"}],"fun_headline_variants":["Strain-tolerant ZrSiS keeps its infrared flatline","Tensile strain triggers Lifshitz switch in ZrSiS","ZrSiS: 10 GPa compression, still flat IR response","UV plasmons in ZrSiS tune with strain","Hyperbolic plasmons and Lifshitz transitions in strained ZrSiS"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the calculated band structure and scalar dielectric function are accurate enough that the quadratic band sits at the right energy near the Fermi level and the permittivity signs near 20 eV are correct; a shift of a few tens of meV would move or erase the tensile Lifshitz transition, and wrong signs would make the hyperbolic regime an artifact of the calculation.","fun_headline_variants_meta":{"raw":{"variants":["Strain-tolerant ZrSiS keeps its infrared flatline","Tensile strain triggers Lifshitz switch in ZrSiS","ZrSiS: 10 GPa compression, still flat IR response","UV plasmons in ZrSiS tune with strain","Hyperbolic plasmons and Lifshitz transitions in strained ZrSiS"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000602,"raw_usage":{"total_tokens":2809,"prompt_tokens":944,"completion_tokens":1865,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1777}},"tokens_in":560,"tokens_out":1865,"duration_ms":13787,"temperature":1.0,"reasoning_tokens":1777,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:30:02.124384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the electron energy loss spectrum of a ZrSiS crystal under uniaxial strain: the paper predicts a sharp low-loss peak near 20 eV with anisotropic quadratic dispersion, and a separate check is that the in-plane and out-of-plane permittivities obtained from optical ellipsometry should have opposite signs in a ~0.6 eV window near 20 eV; if the 20 eV loss peak is absent or the permittivity signs are the same, the high-energy plasmon and hyperbolic claims fail. For the Lifshitz transition, quantum oscillation or ARPES measurements on a crystal under ~2 GPa tensile stress should show the hole pockets merging near $k_z=\\pi/c$ and the quadratic band crossing the Fermi level along Z-R.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental infrared conductivity, flat-region value, and screening constant used to calibrate the low-energy calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Explains the flat optical conductivity as a combination of intraband and interband transitions, the mechanism the paper analyzes under strain."},{"cited_title":"P.; Nateprov, A","cited_arxiv_id":null,"evidence_quote":"Establishes the optical properties of the ZrSiX family and the chemical-pressure trend with c/a, the comparison for strained ZrSiS."},{"cited_title":"B.; van Heumen, E.; Carbone, F.; van der Marel, D","cited_arxiv_id":null,"evidence_quote":"Supplies the power-law scaling for interband conductivity that motivates the frequency-independent behavior in Dirac-like systems."},{"cited_title":"M.; Ali, M","cited_arxiv_id":null,"evidence_quote":"Identifies ZrSiS as a nodal-line semimetal with two kinds of Dirac nodal lines, the band-structure context for the optical analysis."},{"cited_title":"L.; Mao, Z","cited_arxiv_id":null,"evidence_quote":"Reports the lower screening constant of ZrSiTe, which supports the chemical-pressure interpretation of stretched ZrSiS."},{"cited_title":"A.; Yerger, C","cited_arxiv_id":null,"evidence_quote":"Motivates the strain study by showing a possible topological phase transition in ZrSiS under high pressure."},{"cited_title":"D.; Lifshitz, E","cited_arxiv_id":null,"evidence_quote":"Describes a four-point bending setup that could apply the tensile stresses the paper predicts to access the Lifshitz transition."}],"review_version":1}