{"id":"c739197a-d60e-4663-a36c-05a3bd435076","arxiv_id":"1908.04004","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors show that a GGA+U calculation with U above 1 eV turns insulating DFT UTe2 metallic, and the resulting Fermi surfaces imply time-reversal-invariant topological superconductivity with Majorana surface states for all four odd-parity pairing channels.","lead":"This paper uses a first-principles method with an adjustable Coulomb parameter to predict that UTe2 becomes metallic once electron correlations are strong enough, and that its Fermi surfaces then make it a topological superconductor for every odd-parity pairing symmetry allowed by the crystal.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that UTe2 is topological for all odd-parity pairings rests on GGA+U Fermi surfaces at U>1 eV; since U is undetermined and a competing calculation at U=0.51 eV gives an insulator, the predicted invariants are not pinned down.","rationale":"The reader's weakest assumption identifies the same load-bearing premise: the GGA+U method with an adjustable U, J=0, and AMF double counting determines the Fermi surfaces, and the physical U lies above 1 eV. My independent reading agrees. The topological step itself is well founded: given odd-parity, time-reversal-invariant pairing and the stated Fermi surfaces, the winding-number and Z2 formulas follow from published work, and the gap-node classification is drawn from established EAZ classifications. The paper also honestly flags the U ambiguity and suggests experimental checks. Therefore this is a correctness risk about the input Fermi surfaces, not an internal contradiction. The conditional verdict is the right one: the claim is well formed and plausible, but it is not a parameter-free derivation, and the absence of deposited input data prevents independent numerical verification. No adjustment to the reader's verdict is needed.","tokens_in":12624,"tokens_out":12189,"duration_ms":138007,"concrete_test":"Independently recompute the UTe2 GGA+U calculation with the same structure and WIEN2k settings but using the self-interaction-corrected (SIC/Liechtenstein) double-counting formula at U=1.0 and 1.1 eV with J=0.7 eV, then extract the occupation numbers at the eight TRIM in the folded BZ. If the resulting Fermi-surface topology is not regions (i)/(ii), or if any of the eight n(Ki) values in Table II change, then the predicted (omega, nu1, nu2, nu3)=(1,1,1,1) for all odd-parity pairings is not robust to the undetermined U and double-counting choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is not internally inconsistent; it is a parameter-dependent prediction. Table II's topological invariants are read off from GGA+U Fermi surfaces that exist only for U>1.0 eV with J=0 and the around-mean-field double-counting correction. The paper states explicitly: 'we cannot determine the value of U in the framework of the GGA+U method.' A competing DFT+U calculation at U=J=0.51 eV cited in the paper itself yields an insulator, so below the threshold the normal state has no Fermi surface and the entire analysis does not apply. Because U=1.0 eV sits at the insulator-metal transition, the region-(i) Fermi surfaces are tiny, and the n(Ki) values in Table II are sensitive to the double-counting scheme, k-mesh, and other numerical choices. No input files or raw Fermi-surface data are deposited, so the Table II entries cannot be independently recomputed. Thus the strongest claim holds only if the physical low-temperature 5f state is described by U above 1 eV with this particular GGA+U prescription — a premise the paper does not establish.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports GGA+U electronic structure calculations for the heavy-fermion superconductor UTe2. It finds an insulator-metal transition driven by the Coulomb U parameter, with metallic Fermi surfaces appearing for U > 1.0 eV, and classifies the resulting Fermi-surface topologies into three regions (i)-(iii). Assuming odd-parity pairing and time-reversal symmetry, the authors use Fermi-surface formulas to compute the 3D winding number and three 2D Z2 invariants in the folded Brillouin zone, concluding that for regions (i) and (ii) all four odd-parity pairing symmetries of the Immm space group give topological superconductivity with Majorana surface states on (100), (010), and (001). They also present a symmetry-based classification of gap nodes using effective Altland-Zirnbauer classes and propose experimental tests, including thermal-conductivity anisotropy, to determine the pairing symmetry. A schematic field-temperature phase diagram under magnetic fields along the b axis is also proposed.","tokens_in":12796,"tokens_out":5540,"duration_ms":59521,"significance":"If the metallic Fermi surfaces obtained at U > 1 eV describe the actual low-temperature state of UTe2, the paper provides a concrete and falsifiable prediction of an intrinsic time-reversal-invariant (class DIII) topological superconductor with Majorana surface states, a rare and sought-after situation. The use of well-established Fermi-surface formulas for topological invariants and the systematic EAZ classification of gap nodes are rigorous tools that give the paper a clear logical structure. The predictions for gap structures and surface states are testable by STM, ARPES, and thermal transport. However, the central result is conditional on the choice of the Coulomb parameter U, whose physical value is not determined in the paper, and the competing DFT+U calculation at lower U gives an insulating state. The strength of the conclusion therefore depends on an unestablished premise.","major_comments":[{"comment":"The entire topological prediction rests on the insulator-metal transition at U = 1.0 eV, yet the paper explicitly states (p.2) that \"we cannot determine the value of U in the framework of the GGA+U method.\" The competing DFT+U calculation at U = J = 0.51 eV cited as Ref. [42] yields an insulating state, and the bare DFT result is also insulating. Since Table II's invariants (1,1,1,1) for regions (i) and (ii) require the metallic Fermi surfaces, and region (iii) is trivial for most pairings, the claim that UTe2 is a topological superconductor for all odd-parity pairings is not established for the physical material unless an independent constraint on U is provided. The authors should either derive such a constraint (e.g., from specific heat, quantum oscillations, or ARPES) or explicitly frame the result as a conditional prediction over the unknown U window and identify experimental measurements that would discriminate between the regions.","section":null},{"comment":"The occupation numbers n(Ki) in Table II, which determine the parity of the topological invariants, are read off from Fermi surfaces that are tiny near the transition (region (i), U = 1.0 eV). The paper provides no convergence tests with respect to the k-mesh, no comparison of the two double-counting schemes beyond a statement in the Supplemental Material that results were crosschecked, and no deposited input files or raw Fermi-surface data. Since the parity of n(Ki) can change if a small pocket is missed or spuriously created by numerical noise, the authors should demonstrate stability of the Table II entries against numerical parameters (e.g., k-mesh density, RMTKmax, and double-counting prescription) and, ideally, provide the calculated occupancies at all TRIM for each U value.","section":null},{"comment":"The gap-node classification in Table S2 is derived under the assumption that the high-symmetry planes and axes intersect the normal-state Fermi surfaces (stated in the Supplemental Material). The paper does not explicitly verify this condition for the Fermi surfaces in Figs. 2(b)-(d). For instance, a point node on the kz axis (Λ) for the B1u state only affects the excitation spectrum if the Fermi surface actually crosses that axis. If the small electron and hole pockets in regions (i) and (ii) do not intersect the relevant symmetry lines, the predicted point nodes would be absent, and the associated surface Majorana states would not appear as claimed. The authors should check the intersection of each symmetry line and plane with the calculated Fermi surfaces and update Table III accordingly.","section":null}],"minor_comments":[{"comment":"The abstract states \"topological superconductivity at an intermediate U for all the odd-parity pairing symmetry,\" but the main text limits this claim to regions (i) and (ii); for region (iii) only the Au state may be topological with an even winding number. This wording should be made more precise to avoid overstatement.","section":null},{"comment":"The explanation of the folded Brillouin zone and the correspondence in Table I is terse. A short derivation or a sketch showing how the doubled unit cell is compatible with the (100), (010), (001) surfaces would help readers trust the use of the folded BZ for the FS formulas.","section":null},{"comment":"The caption says \"The values n(Ki)-180 are shown below,\" but the meaning of the constant 180 and why it is subtracted are not explained in the text. Please clarify.","section":null},{"comment":"The phase diagram in Fig. 3 is schematic and based on symmetry considerations plus qualitative energy arguments. The paper should state explicitly that the proposed two-phase structure is a conjecture, not a result of a microscopic calculation.","section":null},{"comment":"The paragraph on magnetism suggests nesting-induced finite-q fluctuations coexisting with ferromagnetic fluctuations. This is an interesting speculation, but it is not developed quantitatively; please indicate whether this is a suggestion or a derived result.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper presents a clearly structured and technically sound calculation of topological invariants once the Fermi surfaces are accepted, but the central prediction hinges on an undetermined Coulomb U. The authors should be encouraged to deposit their band-structure inputs and computed occupations so that the Table II values can be independently checked, and to engage more directly with the competing insulating result at lower U. The gap-node verification against the actual Fermi-surface location is a substantive technical issue that may alter the experimental signatures."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, clearly written theory paper that does what it says: it shows a GGA+U calculation for UTe2 that produces metallic Fermi surfaces for U above 1 eV, and then, assuming odd-parity pairing, it derives the topological invariants and gap-node structures for all four odd-parity representations of Immm. The topological classification is standard and the Fermi-surface formulas are applied carefully, including the fold into the double-sized cell, which matters here. The gap-node table is systematically derived from the effective Altland-Zirnbauer classification. The paper is also honest: it states explicitly that GGA+U cannot fix U, and it cites the competing calculation at U=J=0.51 eV that finds an insulator. It frames the results as predictions to be tested by ARPES, quantum oscillations, thermal conductivity, and STM, which is appropriate.\n\nThe main soft spot is exactly the one the reader flagged: the whole edifice rests on the U parameter. The metallic FSs exist only for U>1 eV with J=0 and the around-mean-field double-counting. That is a narrow window, and the competing calculation at lower U gives no FSs at all, so the topological invariants in Table II are not robust to that uncertainty. Region (i) at U=1.0 eV sits right at the transition with tiny pockets, so the occupation numbers there are sensitive to numerical details. No input files or raw FS data are deposited, so the Table II entries cannot be independently recomputed. This is a genuine limitation, but I don't think it is fatal: the paper is transparent about the parameter dependence, and the claim is conditional, not absolute. It is a prediction, not a derivation. The field is better served by having this on the table, with the caveats stated, than by waiting for a DMFT calculation that might not come soon.\n\nA couple of minor quibbles: the phase diagram under field is schematic, and the thermal conductivity anisotropy table should be read as rough guidance. The nesting discussion is speculative but interesting. None of this undercuts the central message.\n\nWho should read this? Anyone working on UTe2, heavy-fermion superconductivity, or topological superconductivity in uranium compounds. It is the kind of paper that will be cited frequently, even if some of the specific predictions turn out to be wrong.\n\nMy recommendation: definitely send it to peer review. It deserves a serious referee, and with a careful reading it should be published, perhaps conditionally, with the authors asked to deposit the DFT input/output data and to discuss the U-dependence and double-counting sensitivity in more detail. I would accept this paper for publication.","headline":"A timely, honest prediction paper mapping UTe2's odd-parity pairing onto topological invariants, but the metallic Fermi surfaces it depends on only exist for a Coulomb U the paper cannot pin down.","tokens_in":13360,"tokens_out":4313,"would_cite":true,"duration_ms":40653,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper predicts that electron correlations turn UTe2 into a time-reversal-invariant topological superconductor for every odd-parity pairing symmetry allowed by its crystal structure.","keywords":["UTe2","topological superconductor","odd-parity pairing","GGA+U","Majorana surface states","Fermi surface","insulator-metal transition","gap nodes"],"falsifier":"Measure the Fermi surface of UTe2 by quantum oscillations or ARPES: if the observed sheets match the small Fermi surfaces of the plain DFT calculation, or if the material is insulating at low temperatures, the predicted winding numbers and Majorana surface states do not follow. A complementary check is low-temperature thermal conductivity: for regions (i) and (ii), B2u pairing requires $\\kappa_b > \\kappa_{a,c}$ and B3u requires $\\kappa_a > \\kappa_{b,c}$, so observing different anisotropies would rule out the predicted gap-node structure.","tokens_in":12386,"feed_emoji":"⚛️","tokens_out":6380,"duration_ms":58554,"temperature":0.7,"pith_summary":"The paper argues that once electron correlations are included through a GGA+U correction, the uranium superconductor UTe2 is a metal with large Fermi surfaces rather than the insulator predicted by plain density-functional theory, and that this metallic state makes the superconductivity topologically nontrivial for every odd-parity pairing symmetry allowed by its Immm crystal structure. This matters because UTe2 is one of the strongest known candidate spin-triplet superconductors, and intrinsic time-reversal-invariant topological superconductors with Majorana surface states are rare. The paper translates the topology into concrete gap structures and measurable signatures, including a thermal-conductivity anisotropy that can distinguish the four pairing channels.","feed_headline":"Coulomb repulsion turns UTe2 into a topological superconductor","feed_subtitle":"All four odd-parity pairing symmetries become topologically nontrivial, with Majorana states on the surface.","key_machinery":"The load-bearing identity is the Fermi-surface formula for the winding number, $\\omega = \\tfrac{1}{2}\\sum_{\\mathbf{K}_i} n(\\mathbf{K}_i) \\pmod 2$, together with analogous formulas for the $\\mathbb{Z}_2$ invariants $\\nu_1$, $\\nu_2$, and $\\nu_3$ on the $k_x=0$, $k_y=0$, and $k_z=0$ time-reversal-invariant planes, evaluated in the folded Brillouin zone because the doubled unit cell is compatible with the (100), (010), and (001) surfaces. The occupation numbers come from GGA+U band structures, so the chain runs from the correlation-strength parameter U to Fermi-surface topology to superconducting topology. Gap structures are then classified by the effective Altland-Zirnbauer symmetry class on high-symmetry mirror planes and rotational axes.","core_discovery":"The central discovery is that a Coulomb interaction strength U above about 1.0 eV drives an insulator-metal transition in UTe2, producing Fermi surfaces whose occupation numbers at the eight time-reversal-invariant momenta make the three-dimensional winding number and all three two-dimensional Z2 invariants nontrivial for the moderate-U Fermi surfaces labeled (i) and (ii). Consequently, for each of the four odd-parity irreducible representations Au, B1u, B2u, and B3u, the gapped superconducting bulk is topologically nontrivial, and Majorana surface states appear on the (100), (010), and (001) surfaces. The gap structure depends on pairing symmetry: Au is fully gapped, B1u has point nodes on the kz axis, B2u on the ky axes, and B3u on the kx axes for regions (i) and (ii), while in region (iii) only Au remains topological with an even winding number and B1u becomes fully gapped.","pith_inferences":["If confirmed, UTe2 would become a material platform for exploring non-Abelian Majorana quasiparticles in a three-dimensional superconductor; the paper itself does not discuss quantum-computation applications.","The results suggest a practical rule of thumb for actinide heavy-fermion compounds where plain DFT gives a small gap: a GGA+U Fermi surface at intermediate U may be enough for a topological classification. A direct dynamical-mean-field calculation would test whether the U-dependence of these Fermi surfaces survives a more rigorous treatment of correlations.","Because the topology holds for all four odd-parity pairings at moderate U, the paper implies that the topological superconductivity is controlled by the Fermi-surface geometry rather than by any one pairing mechanism. A testable extension is to feed experimentally measured Fermi surfaces, from quantum oscillations or ARPES, into the same winding-number and Z2 formulas without assuming a specific U"],"forward_implications":["If the central claim is right, UTe2 is an intrinsic time-reversal-invariant class DIII topological superconductor, and clean (100), (010), and (001) surfaces should show zero-energy Majorana states detectable by tunneling spectroscopy or ARPES.","The predicted pairing-symmetry dependence of gap nodes gives distinct thermal-conductivity anisotropies: B1u gives $\\kappa_c > \\kappa_{a,b}$, B2u gives $\\kappa_b > \\kappa_{a,c}$, and B3u gives $\\kappa_a > \\kappa_{b,c}$ for regions (i) and (ii), so a single bulk measurement can select the pairing channel.","Under a magnetic field along the b axis, the theory predicts two superconducting phases, a line-nodal Bu state at low field and a point-nodal Au state at high field, offering an explanation for the re-entrant superconductivity and a testable field dependence of specific heat and Knight shift.","The nesting of the calculated Fermi surfaces implies that UTe2 sits near multiple magnetic instabilities, which is consistent with the absence of magnetic order despite strong ferromagnetic fluctuations."],"supporting_citations":[{"why":"Supplies the existing DFT band structure with a small gap that the paper's correlation-corrected calculation overturns.","marker":"[5]"},{"why":"Provides a competing first-principles calculation with small Fermi surfaces and ARPES data that the paper argues are incompatible with transport.","marker":"[22]"},{"why":"Gives transport evidence for a large carrier density that favors the large Fermi surfaces obtained at intermediate U.","marker":"[23]"},{"why":"The first-principles band-structure code used for the GGA and GGA+U calculations.","marker":"[24]"},{"why":"Introduces the GGA+U method that the paper uses to incorporate electron correlations.","marker":"[25]"},{"why":"Supplies the Fermi-surface formula for the parity of the three-dimensional winding number.","marker":"[11]"},{"why":"Provides the Fermi-surface formulas for the two-dimensional Z2 invariants used to identify Majorana surface states.","marker":"[13]"},{"why":"Gives the point-group classification of ordering parameters used to label the odd-parity pairing channels.","marker":"[17]"},{"why":"Provides the effective Altland-Zirnbauer classification used to determine symmetry-protected gap nodes.","marker":"[21]"},{"why":"Reports the NMR Knight-shift decrease below Tc that the paper uses to constrain the d-vector orientation under magnetic fields.","marker":"[39]"}],"fun_headline_variants":["Coulomb U drives UTe2 into topological superconductivity","UTe2 topological superconductivity enabled by metal transition","All odd-parity pairings in UTe2 turn topological at moderate U","First-principles UTe2: insulator-metal transition yields Majorana states","UTe2: Coulomb-induced metal state hosts topological superconductivity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that a corrected density-functional calculation with an adjustable Coulomb repulsion U above about 1 eV captures the real low-temperature electronic state of the uranium 5f electrons, and the paper itself states that the GGA+U method cannot determine U.","fun_headline_variants_meta":{"raw":{"variants":["Coulomb U drives UTe2 into topological superconductivity","UTe2 topological superconductivity enabled by metal transition","All odd-parity pairings in UTe2 turn topological at moderate U","First-principles UTe2: insulator-metal transition yields Majorana states","UTe2: Coulomb-induced metal state hosts topological superconductivity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000231,"raw_usage":{"total_tokens":1490,"prompt_tokens":953,"completion_tokens":537,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":449}},"tokens_in":569,"tokens_out":537,"duration_ms":6140,"temperature":1.0,"reasoning_tokens":449,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:54:47.050782+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Fermi surface of UTe2 by quantum oscillations or ARPES: if the observed sheets match the small Fermi surfaces of the plain DFT calculation, or if the material is insulating at low temperatures, the predicted winding numbers and Majorana surface states do not follow. A complementary check is low-temperature thermal conductivity: for regions (i) and (ii), B2u pairing requires $\\kappa_b > \\kappa_{a,c}$ and B3u requires $\\kappa_a > \\kappa_{b,c}$, so observing different anisotropies would rule out the predicted gap-node structure.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the existing DFT band structure with a small gap that the paper's correlation-corrected calculation overturns."},{"cited_title":"Fujimori, I","cited_arxiv_id":null,"evidence_quote":"Provides a competing first-principles calculation with small Fermi surfaces and ARPES data that the paper argues are incompatible with transport."},{"cited_title":"Blaha, K","cited_arxiv_id":null,"evidence_quote":"The first-principles band-structure code used for the GGA and GGA+U calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the GGA+U method that the paper uses to incorporate electron correlations."},{"cited_title":"Sato, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the Fermi-surface formula for the parity of the three-dimensional winding number."},{"cited_title":"Sumita, T","cited_arxiv_id":null,"evidence_quote":"Provides the effective Altland-Zirnbauer classification used to determine symmetry-protected gap nodes."},{"cited_title":"Nakamine, S","cited_arxiv_id":null,"evidence_quote":"Reports the NMR Knight-shift decrease below Tc that the paper uses to constrain the d-vector orientation under magnetic fields."}],"review_version":1}