{"id":"148d7c1f-6469-4583-8f9b-782872ecc3c4","arxiv_id":"1908.06494","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"First-principles calculations show YH3 hosts a pressure-tunable pseudo nodal surface semimetal that transitions through a shrinking nodal ring to a trivial insulator above 31 GPa.","lead":"This paper uses computer simulations to show that the electronic structure of the hydrogen-rich material YH3 changes in a specific way as pressure is applied. The results suggest a pressure-driven transition between different topological semimetal phases, which could inform how to tune electronic properties of materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unjustified neglect of spin-orbit coupling: Y is a 4d element, yet the BDI/AI classification and the nodal ring vanish if SOC gaps the 5-meV 'pseudo' degeneracies.","rationale":"The reader's weakest assumption (SOC neglect) is the single most load-bearing concern. The paper's symmetry analysis uses the spinless realization of the AZ+I classification (Eq. 5: T=K with T^2=1), and the entire BDI labeling depends on this. Yttrium is Z=39 with 4d valence orbitals; atomic 4d SOC splittings are typically tens of meV, far above the E_error=0.005 eV tolerance used to identify the pseudo nodal surface in Fig. 1d. If SOC is included, T^2=-1, so the localized T*I operator squares to -1, the symmetry class changes (e.g., AII), and the nodal ring may be gapped. The paper never performs that check, and its stated reason ('Y and H are light elements') is incorrect for Y. A single VASP calculation with LSORBIT=.TRUE. would settle whether the zero-pressure ring survives; without it, the central claim is not established. I also considered the structural phase transition at 21 GPa, but because the paper explicitly studies the hexagonal phase and even notes that 31 GPa exceeds the transition pressure, that issue is secondary: it qualifies the pressure range but does not invalidate the zero-pressure classification. The reader's CONDITIONAL verdict is appropriate; my concern does not change the verdict but reinforces the condition.","tokens_in":8892,"tokens_out":9280,"duration_ms":92033,"concrete_test":"Repeat the VASP calculations with spin-orbit coupling (e.g., LSORBIT=.TRUE. and a full-relativistic PAW potential for Y, using the same PBE functional and 25x25x25 k-mesh) at zero pressure and at P=28, 31, 32 GPa. Evaluate the minimum direct gap between the LUCB and HOVB on the k-space surface/ring identified in Fig. 1d. If the minimum gap along the ring exceeds 0.005 eV (or the integrated charge of the ring changes), the BDI/AI classification and the pressure-driven trivialization claim are invalidated. If the ring persists, report the new symmetry class (e.g., AII) and the corresponding Z2 invariant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that YH3 at zero pressure is a pseudo nodal surface semimetal in class BDI, and that the nodal ring annihilates around 31 GPa, is computed entirely without spin-orbit coupling (SOC). The authors justify this by asserting 'both Y and H are light elements' (Introduction, p.2). This is factually incorrect for yttrium (Z=39, a 4d transition metal). Atomic 4d spin-orbit splittings are typically tens of meV (e.g., ~50 meV for Y 4d states), which is an order of magnitude larger than the E_error=0.005 eV tolerance used to define the pseudo nodal surface in Fig. 1d. The entire symmetry classification is spinless: the operators in Eq. (5) use T=K, T^2=1, which is the spinless realization. With SOC, T^2=-1, so the BDI/AI classes of the AZ+I tenfold table no longer apply; the localized T*I operator squares to -1, and the nodal ring may be gapped. Nowhere does the paper compute the band structure with SOC to check whether the LUCB-HOVB gap remains closed on the purported nodal ring. If SOC opens a gap of tens of meV at the zero-pressure 'pseudo nodal surface,' the material is a trivial insulator at all pressures, and the central claim is an artifact of the SOC-free approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a first-principles study of hexagonal P-3c1 YH3 under hydrostatic pressure. Using PBE DFT with Wannier interpolation, the authors identify two bands near the Fermi level whose crossings form, within a 5 meV tolerance, a 'pseudo nodal surface' around the Γ point. They assign the zero-pressure crossings to class BDI of the AZ+I classification of Bzdusek and Sigrist, based on an approximate particle-hole symmetry; under pressure, they argue that this symmetry is gradually broken, leaving an AI nodal ring that shrinks and closes at about 31 GPa, after which YH3 becomes a trivial insulator. The paper presents no fitted parameters, uses standard computational tools, and gives a concrete pressure-evolution scenario.","tokens_in":9191,"tokens_out":8723,"duration_ms":94949,"significance":"If substantiated, the paper would add YH3 to the short list of concrete material realizations of the centrosymmetric AZ+I classification and would provide a pressure-driven route from a nodal-surface-like semimetal to a trivial insulator. The DFT setup is standard, the symmetry analysis with Bilbao irreps is detailed, and the WannierTools-based node search is reproducible in principle. The paper also makes a falsifiable prediction: the nodal ring annihilates near 31 GPa. However, the central classification is currently conditional on three load-bearing approximations: spin-orbit coupling is neglected without quantitative justification, the particle-hole symmetry is identified from the very band structure it is used to classify, and the claimed trivial Z2 charge is inferred rather than computed. These are testable and fixable, but they need to be addressed before the central claims can be accepted.","major_comments":[{"comment":"The entire AZ+I classification is carried out in the spinless realization, with T=K and T^2=1 in Eq. (5), justified by the statement that 'both Y and H are light elements.' This justification is inaccurate for yttrium (Z=39, a 4d transition metal), and the pseudo nodal surface is defined with an energy tolerance E_error=0.005 eV in Fig. 1d, which is an order of magnitude smaller than typical 4d spin-orbit splittings (tens of meV). If spin-orbit coupling opens a gap of that size on the purported nodal surface or nodal ring, the BDI/AI assignment and the pressure-driven trivialization would not apply as described. The manuscript contains no SOC-included calculation. Please provide band structures with SOC at the relevant pressures (0, 28, 31, and 32 GPa) and report the gap on the purported nodal ring, or otherwise quantitatively demonstrate that SOC effects are below the 0.005 eV tolerance used to define the nodes.","section":"Introduction, Eq. (5), Fig. 1d"},{"comment":"The existence of particle-hole symmetry P is inferred from the approximate mirror symmetry of the two bands about the Fermi level and is then used to assign the BDI class. This is partly circular, and no operator P is constructed or tested in the DFT/Wannier basis. Because P is only approximate, the exact band crossings are nodal lines rather than a true nodal surface, so the 'class BDI' label is a heuristic classification of near-degeneracies, not a symmetry-protected topological statement. Please quantify the P-breaking, for example by giving the norm of the anticommutator {P,H} or the maximal deviation of the two-band spectrum from particle-hole symmetry over the Brillouin zone, specify the action of P on the basis states used for the Wannier interpolation, and show how this measure evolves from 0 to 32 GPa. Without this, the zero-pressure 'class BDI' claim is not established at the quantitative level needed for the paper's title claim.","section":"pp. 8-9, Eqs. (3)-(6)"},{"comment":"The paper concludes that the nodal ring is topologically trivial and that the Berry phase for all occupied bands is quantized to 0, based solely on the observation that the ring shrinks continuously and gaps out without sudden changes. A continuous annihilation is not a proof of zero Berry phase; a nontrivial ring could also disappear through pair annihilation with another ring or through gap-closing events elsewhere in the Brillouin zone. Please compute the Berry phase, or the Z2 invariant of Ref. 13, on a loop enclosing the nodal ring at several pressures, or provide an independent symmetry-based argument. This is load-bearing for the claimed 'topological phase transformation' and for the statement that YH3 becomes a trivial insulator.","section":"p. 10, paragraph beginning 'To figure out whether...'"},{"comment":"All pressure-dependent calculations are performed in the P-3c1 structure up to 32 GPa, although the cited Ref. 21 places a structural transition to a cubic phase at 21 GPa. The claims for 28-32 GPa therefore describe a metastable or hypothetical phase rather than the ground-state material. Please either demonstrate that P-3c1 remains (meta)stable in this pressure range, or restrict the central claims to pressures below 21 GPa, or extend the calculation to the high-pressure phase. The abstract's 'above 31 GPa' statement needs to be qualified accordingly.","section":"Abstract and Fig. 2"}],"minor_comments":[{"comment":"The manuscript contains numerous typos and grammatical errors, including 'A lots of progress', 'can be understand', 'classification sheme', 'oftenly', the corrupted character sequence 'BerryâAZs', and 'zhe AZ+I'. A thorough language edit is needed.","section":"Throughout"},{"comment":"The definition of nodes by E_LUCB - E_HOVB < 0.005 eV should be accompanied by a convergence check of the Wannier interpolation at this energy scale, since 5 meV is close to typical interpolation errors. Please report the interpolation accuracy and how the node count in Fig. 1d depends on the tolerance near this value.","section":"Fig. 1d"},{"comment":"The relationship between the non-spatial-symmetry protection proposed here and the crystalline-symmetry protection proposed in Refs. 22 and 23 is asserted but not developed. Since those papers attribute the same or similar crossings to glide-plane or mirror symmetries, the manuscript should explicitly state whether the non-spatial protection is compatible with, or an alternative to, the crystalline-symmetry protection, and how the two descriptions can be distinguished in the band structure.","section":"Introduction, Refs. 22-23"},{"comment":"The two-band Hamiltonian in Eq. (6) omits a σy term. The text should state explicitly that this follows from the spinless time-reversal and inversion symmetries, which make the Hamiltonian real in the chosen basis, so that the reader can follow the symmetry constraints without additional derivation.","section":"Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The main risk to the paper is the neglected spin-orbit coupling: yttrium is a 4d transition metal, and the 5 meV node-search tolerance is small compared with typical 4d SOC splittings. If a SOC-included calculation shows that the nodes survive, the paper could be publishable after the requested Berry-phase computation and a more careful discussion of the approximate particle-hole symmetry. I would also ask the authors to reconcile their non-spatial-symmetry picture with the earlier glide-plane and mirror-symmetry interpretations cited in Refs. 22 and 23."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper gives a concrete scenario for pressure-driven evolution of a nodal ring in YH3, but the whole classification is done in the spinless limit, and for yttrium—a 4d element—that limit is not safe at the 5 meV scale the nodes are defined on. The authors say \"both Y and H are light elements\" (Introduction), which is simply not true for Y. Until they run the same calculation with SOC and show the nodes survive, the BDI/AI label and the trivialization story should be read as a hypothetical, not a property of YH3.\n\nWhat is genuinely useful: the symmetry analysis is careful. They use BANDREP to track the e- and h-band irreps, show the crossings arise from two-band overlap, and they make a clean pedagogical point about how approximate particle-hole symmetry upgrades a nodal ring (class AI) to a pseudo nodal surface (class BDI) in the AZ+I periodic table. The pressure dependence of the gap, with the ring shrinking and closing around 31 GPa, is directly visible in the band structures. That part of the paper is fine and well presented.\n\nThe soft spots, in order of severity. First, the SOC issue above; it is load-bearing because with T^2=-1 the T*I operator squares to -1 and the entire classification changes, and the pseudo-gap tolerance of 5 meV is small compared to expected 4d SOC of tens of meV. Second, the Z2 invariant is not computed; it is inferred because the ring shrinks and gaps out without a sudden change, which is suggestive but not a proof. Third, the threshold E_error = 0.005 eV is hand-set; the 'fuzzy surface' away from the exact nodal line is tolerance-dependent, and the BDI label only exists within that tolerance. Fourth, they study pressures up to 32 GPa even though the structure transitions at 21 GPa (their own Ref. 21); they don't address whether the P-3c1 phase is metastable there.\n\nThe paper is worth a serious referee because the question—can a realistic material exhibit an approximate BDI nodal surface, and does pressure trivialize it?—is legitimate and the authors provide a testable path. But the referee should insist on an SOC-inclusive calculation, a direct computation of the Berry phase / Z2 invariant, and a discussion of the structural stability range. With those, the story could be solid; without them, it remains a spinless model.","headline":"A plausible spinless story about pressure-driven trivialization in YH3, but the SOC neglect for a 4d element leaves the central claim unverified.","tokens_in":9681,"tokens_out":3642,"would_cite":false,"duration_ms":35996,"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 paper argues that hexagonal YH3 is a pseudo nodal surface semimetal in class BDI at zero pressure and that pressure above 31 GPa turns it into a trivial insulator.","keywords":["topological semimetal","nodal surface","nodal ring","yttrium trihydride","YH3","AZ+I classification","pressure-induced transition","first-principles calculation"],"falsifier":"A relativistic density-functional band calculation of $P\\bar{3}c1$ YH3, including spin-orbit coupling, at ambient pressure and at 28–32 GPa, resolving energies around the nodal ring; if the ring acquires a gap at any pressure below 31 GPa, the claimed pressure-driven trivialization fails, whereas survival of the ring would support the pseudo nodal surface picture.","tokens_in":8700,"feed_emoji":"🔬","tokens_out":6441,"duration_ms":63823,"temperature":0.7,"pith_summary":"The paper argues that the band crossings near the Fermi level in hexagonal YH3 are caused by the overlap of an electron-like and a hole-like band, and that their protection comes mainly from non-spatial symmetries rather than crystalline symmetries. Because time reversal, inversion, and an approximate particle-hole symmetry are present, the crossings can be viewed as a pseudo nodal surface in class BDI of the AZ+I classification of gapless topological matter. The word “pseudo” matters: away from a closed nodal line a small gap of about 0.005 eV opens, so the strictly protected object is a nodal ring in class AI. The paper then claims that hydrostatic pressure gradually breaks the particle-hole symmetry, shrinks the nodal ring, and at about 31 GPa collapses it to a point; above 31 GPa all crossings are gapped and YH3 becomes a trivial insulator. If correct, this gives a concrete material in which a pressure-driven electronic topology transition can be followed step by step in first-principles band structure.","feed_headline":"Pressure drives YH3 from topological semimetal to trivial insulator","feed_subtitle":"YH3's electron and hole bands touch at zero pressure; by 31 GPa, calculation says all crossings gap out.","key_machinery":"The load-bearing object is the effective two-band Hamiltonian $H(k)=f(k)\\sigma_z+g(k)\\sigma_x$ with symmetry operators $T=K$, $P=\\sigma_z K$, and $C=\\sigma_z$. When $f(k)=0$, particle-hole symmetry is exact and the nodes form the surface $g(k)=0$, which belongs to class BDI; the small term $f(k)$ breaks particle-hole symmetry and leaves only the nodal ring satisfying $f(k)=0$ and $g(k)=0$, which belongs to class AI. The paper locates this ring numerically by searching for k-points where the lowest unoccupied conduction band and the highest occupied valence band agree to within 0.005 eV, and it tracks the ring radius as a function of pressure.","core_discovery":"The core claim is that at zero pressure YH3 hosts accidental band crossings, not symmetry-enforced crystalline degeneracies: a lowest unoccupied conduction band and a highest occupied valence band overlap and mirror each other approximately across the Fermi level. With the approximate particle-hole symmetry taken as exact, the crossings form a nodal surface belonging to class BDI; once the small symmetry-breaking term is included, the protected object is a closed nodal ring belonging to class AI with a $Z_2$ Berry phase. The paper traces this ring under hydrostatic pressure and finds that it shrinks continuously from 28 GPa to 31 GPa, becomes a point at about 31 GPa, and is fully gapped above 32 GPa. Because the gapping happens smoothly without a sudden level crossing, the paper concludes that the nodal ring carries a trivial $Z_2$ invariant, so the pressure evolution is a topological phase transformation from a semimetal to a trivial insulator.","pith_inferences":["The paper leaves spin-orbit coupling untested; a relativistic calculation would reveal whether the 0.005 eV-scale crossings and the nodal ring survive, and this is the most direct check of the classification.","A similar electron-hole overlap mechanism might occur in other rare-earth hydrides, so the strategy of looking for bands mirrored around the Fermi level could identify additional BDI/AI nodal materials.","Because the particle-hole symmetry is only approximate, the robust experimental prediction is the class-AI nodal ring rather than the class-BDI surface; measurements should target the ring itself.","One could test the pressure evolution experimentally by measuring quantum oscillations or the Berry phase of the occupied bands around the ring as pressure crosses 31 GPa."],"forward_implications":["If the claim is correct, YH3 provides a realistic material realization of a gapless phase from the AZ+I classification, with a pressure knob that tunes it toward a trivial insulator.","The continuous shrinkage and disappearance of the nodal ring without a sudden level crossing establishes that the ring's $Z_2$ Berry phase is 0, meaning this nodal ring is not topologically protected.","The predicted electronic transition near 31 GPa is distinct from the reported structural phase transition near 21 GPa, so the two transitions can be separated and studied independently.","Above about 32 GPa YH3 should behave as a fully gapped trivial insulator, a statement that can be checked by pressure-dependent transport or optical measurements."],"supporting_citations":[{"why":"Supplies the tenfold AZ+I classification table that assigns the BDI and AI classes and gives the dimensionality and topological charges of the nodes.","marker":"[1]"},{"why":"Defines the $Z_2$ Berry phase invariant for nodal lines, which the paper uses to conclude that the YH3 nodal ring is trivial.","marker":"[13]"},{"why":"Provides the low-pressure $P\\bar{3}c1$ structure and the reported structural phase transition pressure of 21 GPa, the reference point for the electronic transition.","marker":"[21]"},{"why":"Earlier claim that YH3 is a nodal-line semimetal protected by glide-plane symmetry, the interpretation this paper revises.","marker":"[22]"},{"why":"Introduced the notion of a pseudo nodal surface with approximate degeneracies, which this paper adapts to describe the class-BDI limit.","marker":"[23]"},{"why":"Supplies the initial crystal structure of YH3 used as the starting point for the first-principles calculations.","marker":"[24]"},{"why":"Provides the Wannier-function construction used to build the tight-binding model from which the nodes are located.","marker":"[40]"},{"why":"Provides the node-finding tool used to identify k-points where the two bands agree within the 0.005 eV error tolerance.","marker":"[41]"}],"fun_headline_variants":["YH3 nodal ring collapses to a point at 31 GPa","Pressure turns YH3 semimetal into trivial insulator","Band crossings in YH3 vanish under pressure","From BDI to AI: YH3's topological phase shift","YH3's zero-pressure crossings gone by 31 GPa"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classification assumes spin-orbit coupling is negligible because yttrium and hydrogen are light elements, but yttrium is a 4d transition metal; if spin-orbit coupling gaps the nodes, the BDI/AI classification and the pressure-driven trivialization would not hold as described.","fun_headline_variants_meta":{"raw":{"variants":["YH3 nodal ring collapses to a point at 31 GPa","Pressure turns YH3 semimetal into trivial insulator","Band crossings in YH3 vanish under pressure","From BDI to AI: YH3's topological phase shift","YH3's zero-pressure crossings gone by 31 GPa"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000964,"raw_usage":{"total_tokens":4116,"prompt_tokens":967,"completion_tokens":3149,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":3066}},"tokens_in":583,"tokens_out":3149,"duration_ms":20004,"temperature":1.0,"reasoning_tokens":3066,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:43:22.795662+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A relativistic density-functional band calculation of $P\\bar{3}c1$ YH3, including spin-orbit coupling, at ambient pressure and at 28–32 GPa, resolving energies around the nodal ring; if the ring acquires a gap at any pressure below 31 GPa, the claimed pressure-driven trivialization fails, whereas survival of the ring would support the pseudo nodal surface picture.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the low-pressure $P\\bar{3}c1$ structure and the reported structural phase transition pressure of 21 GPa, the reference point for the electronic transition."},{"cited_title":"Shao , author T","cited_arxiv_id":null,"evidence_quote":"Earlier claim that YH3 is a nodal-line semimetal protected by glide-plane symmetry, the interpretation this paper revises."},{"cited_title":"Wang , author Y","cited_arxiv_id":null,"evidence_quote":"Introduced the notion of a pseudo nodal surface with approximate degeneracies, which this paper adapts to describe the class-BDI limit."},{"cited_title":"Persson ,\\ 10.17188/1199674 title Materials Data on YH 3 ( SG :165) by Materials Project \\ ( year 2016 ) NoStop","cited_arxiv_id":null,"evidence_quote":"Supplies the initial crystal structure of YH3 used as the starting point for the first-principles calculations."}],"review_version":1}