{"id":"c39e37bd-e12b-4b91-b2a7-0747049622d0","arxiv_id":"2412.15103","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A two-component superconducting order parameter following the lower eigenvector of the Landau free-energy matrix acquires a Berry phase around a degeneracy, reversing the Josephson current.","lead":"This paper shows that a superconducting order parameter with two components can pick up a geometric Berry phase when interaction parameters are cycled in a loop. It connects Landau's classic theory of phase transitions to topological physics and predicts a measurable current reversal in a Josephson junction.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Adiabatic locking to the lower eigenvector is the load-bearing assumption; Eq. (11) is also internally inconsistent as written, so Eq. (13) needs a direct TDGL check before the Josephson signature is accepted.","rationale":"The paper's central mechanism is adiabatic following of the lower eigenvector of the quadratic Landau matrix. If that following fails, the Berry-phase integral in Eq. (13) is no longer the phase acquired by the physical order parameter. This is precisely the condition the paper states but does not verify quantitatively for the loops used in the Josephson calculation. I could not find a more fundamental flaw: the identification of A-eigenvectors with order-parameter directions, the Dirac/Weyl point geometry, and the Josephson sign reversal are all internally consistent under the stated hierarchy. The main concrete weakness is that the derivation of the phase equation is abbreviated; Eq. (11) as printed is not a correct equation of motion for the λ-dependent equilibrium magnitude, which makes the adiabatic elimination argument less transparent than it should be. A direct numerical integration of the TDGL would settle whether the Berry-phase formula holds for experimentally plausible parameters and whether the Josephson-current reversal is robust. The reader's weakest assumption identifies the same physical edge of the argument; I would keep the CONDITIONAL verdict but require the numerical check, or an explicit quantitative bound on the loop parameters, before accepting the experimental signature.","tokens_in":16460,"tokens_out":21708,"duration_ms":208615,"concrete_test":"Integrate the full TDGL equation (Eq. 9 with the explicit free energy of SM I) along the red and green loops of Fig. 2(b) at finite ramp speed, using the same parameters (Λ = 10, N(0)V0 = 0.4, ϕ = 2.7) and a temperature slightly below T_c(r). Extract the phase of Δ after one period and repeat for ramp times tau spanning tau_- to 10^4 tau_- and for several loop radii that cross the locking condition a_+ ~ |B| |Δ_0|^2. Check that the phase tends to π as tau/tau_- → ∞ and quantify the deviation when tau_- ≳ tau or a_+ ≲ |B| |Δ_0|^2. If the deviation is not controlled, the predicted Josephson-current reversal needs a quantitative validity condition.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Eq. (13), depends entirely on the order parameter staying on the instantaneous lower-eigenvector ray Δ̂_-(lambda) of A(lambda), which requires the hierarchy tau_+ << tau_- << tau (before Eq. 10) and a_+ >> |B_{αβγδ}| |Δ_0|^2 (after Eq. 8). These conditions are necessary: if the quartic term reorients the order parameter, or if tau_- is not negligible relative to the ramp time tau, the phase evolution is not the Berry connection of A. The paper states these conditions but never checks them quantitatively for the loops in Fig. 2(b) or for the Josephson current calculation. A related internal issue is that substituting the ansatz (10) into the TDGL equation cannot give dΔ_0/dt = 0 as printed: Δ_0(λ(t)) is time-dependent through λ, and even in the scalar limit the adiabatic solution has Δ̇_0 = Δ_0'(λ)λ̇, not zero. The real-part projection should describe the small lag that slaves the magnitude to the instantaneous minimum; as written, Eq. (11) is inconsistent with Eq. (10) unless the loop is traversed infinitely slowly and subleading corrections are discarded. This does not by itself falsify Eq. (12), but it shows the adiabatic elimination is not fully carried out and strengthens the need for a direct numerical check of the accumulated phase.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an extension of Landau theory, called \"topological Landau theory,\" for multicomponent order parameters transforming under the same irreducible representation. In this framework, the quadratic Landau matrix A(λ) plays the role of a Hamiltonian, the stabilized order parameter is its lowest eigenvector Δ̂_-(λ), and a cyclic adiabatic variation of the thermodynamic parameters λ produces a Berry phase given by Eq. (13). The authors derive a BCS mean-field free energy and a time-dependent Ginzburg-Landau (TDGL) equation for a tetragonal superconductor with two attractive partial waves in the trivial representation. They analyze two models, one time-reversal symmetric with a thermodynamic Dirac point and one time-reversal breaking with a Weyl point, and compute a Josephson current signature of the Berry phase.","tokens_in":16776,"tokens_out":25391,"duration_ms":226229,"significance":"If the central claim holds, this is a genuinely useful conceptual extension of Landau theory: it identifies the quadratic Landau matrix as a parameter-space Hamiltonian and the order-parameter direction as a geometric object with Berry-phase content. The paper has real strengths: the Landau free energy and TDGL equation are derived from a microscopic BCS Hamiltonian; the Berry phase in Eq. (13) follows from the eigenvector of A without fitted parameters; the static prediction that the ground state is close to Δ̂_- near T_c is checked numerically by minimizing the full free energy (M≈1 near T_c); and the Josephson current reversal in Fig. 3 is a concrete, falsifiable signature. The main weakness is that the adiabatic elimination leading to Eq. (13) is not carried out rigorously, and the adiabatic conditions are not quantitatively verified for the loops used in the Josephson calculation.","major_comments":[{"comment":"The derivation of the central result is not internally consistent. Substituting the quasistatic ansatz (10) into the TDGL equation (9) cannot produce Eq. (11), dΔ0/dt = 0, because Δ0(λ(t)) depends on time through λ and its derivative is generally nonzero. More seriously, if the state were exactly on the instantaneous minimum, the right-hand side of Eq. (9) would vanish (using Δ0^2 = -a_-/b_-), forcing dΔ/dt = 0, which is incompatible with the λ-dependence of Δ̂_-(λ) unless Δ0 is constant along the loop. The correct adiabatic elimination must keep the small lag of the magnitude behind the instantaneous minimum; the imaginary part of the projected equation then reproduces the leading-order phase equation (12), while the real part gives the relaxation of that lag rather than Eq. (11). As written, the passage from Eqs. (10)-(12) to the Berry phase (13) is incomplete and should be redone with an explicit expansion in τ_+/τ and τ_-/τ.","section":"Adiabatic dynamics, Eqs. (10)-(12)"},{"comment":"The adiabatic hierarchy τ_+ ≪ τ_- ≪ τ and the condition a_+ ≫ |B_{αβγδ}| |Δ0|^2 are stated but never checked for the actual loops in Fig. 2(b) or for the Josephson current calculation in Fig. 3. Since τ_- diverges at the critical lines a_- = 0 and the topologically nontrivial loops encircle the metallic region at T > T_c^0, the paths may pass close to the critical surface, where the order parameter is not locked to the instantaneous lower eigenvector and Eq. (13) ceases to apply. Moreover, the current in Fig. 3 is computed from the instantaneous equilibrium gap at each parameter point, not from a solution of the TDGL equation, so it does not by itself validate the geometric-phase accumulation. I request a quantitative check of the hierarchy for the paths used, or a direct numerical integration of Eq. (9) along at least one nontrivial loop confirming the accumulated phase and the current reversal.","section":"Topological Josephson effect, Fig. 3"}],"minor_comments":[{"comment":"The statement that Δ̂_- has the lower critical temperature is reversed. From Eq. (7), the eigenvector of V with the larger attractive eigenvalue V0+r has the higher T_c, and that eigenvector is Δ̂_-, the lower eigenvector of A. The same mislabeling appears in the Weyl paragraph, where Δ̂_- is called the less attractive channel.","section":"Dirac and Weyl points, after Eqs. (14) and (15)"},{"comment":"In the time-reversal symmetric model the Berry connection vanishes everywhere, so the π phase is a holonomy due to the double-valuedness of the real eigenvector on the circle rather than a loop integral of A. The text notes this, but the discussion should explicitly distinguish this branch holonomy from the connection-based Berry phase of Eq. (13), because the two mechanisms are conceptually different.","section":"Dirac and Weyl points, TRS model"},{"comment":"The numerical parameters for Fig. 2 are given, but the temperature T at which M≈1 is evaluated, the radii of the red/green/blue loops, and the ramp time τ relative to τ_- are not specified. These values are needed to assess the adiabatic hierarchy and to make the Josephson prediction reproducible.","section":"Figs. 2 and 3"},{"comment":"In deriving the cubic contribution to the TDGL equation, the supplement assumes Δ_k(t1) ≃ Δ_k(t); the validity condition for this low-frequency approximation should be stated explicitly, since it is one of the assumptions behind the adiabatic dynamics.","section":"Supplementary Material II"}],"recommendation":"major_revision","confidential_remarks":"The paper makes a serious and promising contribution, and I do not see a fundamental obstruction to the central idea. The main issue is that the derivation of Eq. (13) is not carried out rigorously: Eq. (11) is not a correct adiabatic equation of motion, and the adiabatic hierarchy is not verified for the paths used in the Josephson calculation. Both issues are fixable within the manuscript's scope. There is also a clear mislabeling of which eigenvector has the higher critical temperature. I would support publication after a careful revision addressing these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core idea is genuinely nice: treat the quadratic Landau matrix A(λ) as a Hamiltonian, the stabilized order parameter as its lower eigenvector, and cyclic variations of λ as loops that pick up a Berry phase. The specific D4h two-partial-wave superconductor, the Dirac and Weyl points in parameter space, and the predicted Josephson current reversal are new relative to the cited momentum-space monopole literature. The microscopic derivation of the Landau free energy and the TDGL equation in the supplement is careful and standard, and the full free-energy minimization confirms M ≈ 1 near T_c, which lends real support to the central eigenvector-locking assumption.\n\nThe paper does have soft spots, and they are worth fixing before publication. The stress-test note is right: Eq. (11) is internally inconsistent as written. Δ₀(λ(t)) carries time dependence through λ, so dΔ₀/dt is not zero in an adiabatic process; it is a small lag term set by λ̇. Printing dΔ₀/dt = 0 as an exact equation is wrong, even if the intended meaning is that the magnitude slavishly follows the instantaneous minimum. This does not invalidate Eq. (12) or the Berry phase formula, but it shows the projection onto the amplitude direction was not carried out cleanly. A corrected adiabatic elimination, keeping the leading λ̇ correction, would remove the inconsistency. There is also a wording slip about critical temperatures: the more attractive channel has a higher T_c, yet the text says “lower transition temperature” in the Weyl case and labels the eigenvectors by “lower and higher critical temperatures” in a way that is easy to read as backwards. That should be clarified.\n\nThe bigger substantive request is a direct numerical check of the accumulated phase via the TDGL equation for the loops in Fig. 2(b) and, if possible, for the Josephson junction calculation. The paper states the adiabatic hierarchy τ₊ ≪ τ₋ ≪ τ but never verifies it quantitatively for the specific parameter values used in the Josephson simulation. The Josephson reversal is a derived consequence, not an input, so the logic is sound, but a numerical integration of the full TDGL for one non-contractible loop would close the loop. No code or data are provided, which makes independent verification slower than it needs to be.\n\nNone of this undermines the central argument. The framework is clear, the derivations are mostly reproducible, and the experimental signature is falsifiable. This is a solid theory Letter that deserves referee time. My recommendation: engage with it, send to a good referee, and expect a conditional acceptance after the adiabatic derivation is cleaned up and a TDGL check is added.","headline":"A genuinely new and mostly sound extension of Landau theory to Berry phases of multicomponent order parameters, with a sloppy adiabatic derivation that needs a fix and a numerical TDGL check.","tokens_in":17271,"tokens_out":3212,"would_cite":true,"duration_ms":21404,"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":"Multi-component Landau order parameters acquire a Berry phase when thermodynamic parameters are cycled adiabatically, and in the models presented this phase is π around a Dirac point and qΩ around a Weyl point, reversing a Josephson…","keywords":["topological Landau theory","Berry phase","order parameter topology","Dirac point","Weyl point","superconductivity","Ginzburg-Landau theory","Josephson effect"],"falsifier":"Numerically integrate the full time-dependent Ginzburg-Landau equation (Eq. 9) without the adiabatic ansatz for a closed loop encircling the Dirac point; if the phase change of the gap after one cycle differs from $\\pi$ by more than the numerical error, the geometric-phase prediction is wrong. A complementary experimental falsifier: a Josephson junction starting from $\\Delta^L_k = -i\\Delta^R_k$ that cycles the right superconductor's parameters around the Dirac point and returns the current to its initial direction rather than reversing it would refute the claimed topological phase.","tokens_in":16269,"feed_emoji":"⚡","tokens_out":17942,"duration_ms":128209,"temperature":0.7,"pith_summary":"Landau theory prescribes a free-energy expansion in powers of the order parameter, and the paper shows that when several components of the order parameter belong to the same irreducible representation, the quadratic coefficient matrix acts like a Hamiltonian: the stabilized order parameter is its lowest eigenvector, and that eigenvector can carry a geometric (Berry) phase as thermodynamic parameters are varied in a closed loop. Working through the time-dependent Ginzburg-Landau equation in the adiabatic limit, the paper finds that in a tetragonal superconductor with two pairing channels of the same symmetry, loops around thermodynamic analogs of Dirac and Weyl points give the order parameter a $\\pi$ phase in the time-reversal-preserving case and a $q\\Omega$ phase in the time-reversal-breaking case, with monopole charge $q = -1/2$. These phases are not dynamical details: they reverse the direction of a Josephson current through a junction, giving a measurable signature. If correct, the paper turns Landau theory into a topological theory, connecting phase diagrams to the geometry of eigenvectors.","feed_headline":"One loop in parameter space flips a supercurrent","feed_subtitle":"Cycling the superconductor's parameters once around the Dirac point gives it a π Berry phase and flips the current.","key_machinery":"The central object is the quadratic free-energy matrix $A(\\lambda)$, whose matrix elements $A_{\\alpha\\beta}$ are the second-order coefficients coupling order-parameter components that transform in the same irreducible representation. Because the free energy is $\\sum_{\\alpha\\beta}\\Delta_\\alpha^* A_{\\alpha\\beta}\\Delta_\\beta$ to this order, the stable ordered state is the normalized eigenvector $\\hat{\\Delta}_-(\\lambda)$ with the smallest eigenvalue, and the adiabatic dynamics of the time-dependent Ginzburg-Landau equation keeps the order parameter on this eigenvector. The Berry connection $A_{-,j} = i\\hat{\\Delta}_-^\\dagger \\partial_j \\hat{\\Delta}_-$ then encodes the geometric phase acquired over a loop, and the degeneracies of $A(\\lambda)$, points where its two eigenvalues coincide, act as Dirac or Weyl points that source Berry curvature in parameter space. The derivation requires the hierarchy of relaxation timescales $\\tau_+ \\ll \\tau_- \\ll \\tau$, where $\\tau_\\pm = \\hbar\\eta/|a_\\pm|$ are the relaxation times along the two eigenvectors and $\\tau$ is the timescale of parameter variation. This machinery converts a thermodynamic phase diagram into a synthetic band structure whose 'bands' are the eigenvectors of $A(\\lambda)$ and whose gaps close at the diabolical points.","core_discovery":"The central claim is that under cyclic adiabatic variation of the thermodynamic parameters $\\lambda$, the superconducting order parameter $\\boldsymbol{\\Delta}$ remains locked to the instantaneous lowest eigenvector $\\hat{\\Delta}_-(\\lambda)$ of the quadratic Landau matrix $A(\\lambda)$ and acquires the Berry phase $\\varphi = \\varphi_0 + \\oint_C i\\hat{\\Delta}_-^\\dagger d\\hat{\\Delta}_-$, the loop integral of the Berry connection (Eq. 13). In the time-reversal-symmetric model $V = V_0 I_2 + r(\\cos\\phi\\,\\sigma_z + \\sin\\phi\\,\\sigma_x)$, the phase diagram contains a Dirac point at $r = 0$; a loop encircling it yields a $\\pi$ Berry phase, equivalently a sign change of the order parameter, while a contractible loop yields zero. In the time-reversal-breaking model $V = V_0 I_2 + r(\\cos\\theta\\,\\sigma_z + \\sin\\theta(\\cos\\phi\\,\\sigma_x + \\sin\\phi\\,\\sigma_y))$, the degeneracy is a Weyl point of monopole charge $q = -1/2$, and a loop subtending solid angle $\\Omega$ yields the Berry phase $q\\Omega$. In both models the geometric phase appears in the Josephson current: after one cycle around the non-contractible loop, the current reverses direction, while a contractible path leaves it unchanged. The paper thereby extends Landau's theory of phase transitions to include the topology of the order parameter.","pith_inferences":["A testable extension: if a loop is traversed faster than the relaxation time $\\tau_-$, the order parameter will lag behind the eigenvector $\\hat{\\Delta}_-(\\lambda)$ and the phase after one cycle should deviate from the geometric value; measuring this deviation would probe the breakdown of the adiabatic locking that the paper's derivation assumes.","The same mechanism suggests that interaction-parameter loops could act as synthetic gauge-field sources in other tunable systems, such as ultracold atoms near a Feshbach resonance, to engineer arbitrary geometric phases on demand rather than only the specific $\\pi$ and $q\\Omega$ values computed here.","The paper does not analyze the critical behavior at the diabolical points themselves; a natural next step would be to compute critical exponents at the Dirac and Weyl points of $A(\\lambda)$, where two critical surfaces touch, to see whether they differ from ordinary Landau critical exponents.","The Josephson reversal is one possible readout; another would be to embed the same superconductor in an interference device where the Berry phase appears as a shift in the critical-current interference pattern, providing an independent experimental route."],"forward_implications":["A multi-component order parameter in Landau theory is generically topologically nontrivial: its cyclic adiabatic evolution produces a geometric phase fixed by the loop's winding around degeneracies, not by the drive rate.","The $\\pi$ Berry phase from the Dirac point means two successive loops return the order parameter to itself, and one loop reverses the direction of the Josephson current; in the Weyl case the current reversal is controlled by the solid angle subtended.","The mechanism generalizes beyond superconductivity: any continuous transition whose order parameter has several components transforming under the same irrep of the symmetry group will exhibit the same geometric phase under parameter cycling.","The geometric phase is fixed by topology: for any loop that does not cross a critical surface or a degeneracy, the phase is determined by whether the loop encircles the point, so perturbations of the path do not change the result."],"supporting_citations":[{"why":"Establishes the Landau free-energy expansion in powers of the order parameter that the paper extends.","marker":"[1]"},{"why":"Classifies superconducting states by pairing-symmetry irreps, motivating two same-irrep components.","marker":"[6]"},{"why":"Supplies the analogy mapping parameter space to momentum space and the quadratic matrix to a Bloch Hamiltonian.","marker":"[8]"},{"why":"Introduces the time-dependent relaxation dynamics of the order parameter near a phase transition.","marker":"[11]"},{"why":"Provides the time-dependent Ginzburg-Landau equation used to derive the adiabatic phase dynamics.","marker":"[12]"},{"why":"Supplies the Berry connection and geometric-phase formalism underlying Eq. (13).","marker":"[16]"},{"why":"Provides the Josephson current formula used to compute the current reversal signature.","marker":"[25]"}],"fun_headline_variants":["Topological Berry phase flips superconducting current","Cycling parameters reverses supercurrent via Berry phase","Adiabatic loop in parameter space gives π Berry phase","Topological Landau theory: Berry phase reverses supercurrent"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the order parameter remains locked to the instantaneous lowest eigenvector $\\hat{\\Delta}_-(\\lambda)$ of the quadratic matrix for the whole closed loop, which requires the high-energy eigenvector to stay negligible and the parameter evolution to be slow compared with the relaxation time $\\tau_-$; near the critical point $\\tau_-$ diverges, the quartic term can reorient the order parameter, and the Berry-phase formula (13) stops applying.","fun_headline_variants_meta":{"raw":{"variants":["Topological Berry phase flips superconducting current","Cycling parameters reverses supercurrent via Berry phase","Adiabatic loop in parameter space gives π Berry phase","Topological Landau theory: Berry phase reverses supercurrent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00141,"raw_usage":{"total_tokens":5732,"prompt_tokens":1014,"completion_tokens":4718,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":4656}},"tokens_in":630,"tokens_out":4718,"duration_ms":27052,"temperature":1.0,"reasoning_tokens":4656,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:38:57.458593+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the full time-dependent Ginzburg-Landau equation (Eq. 9) without the adiabatic ansatz for a closed loop encircling the Dirac point; if the phase change of the gap after one cycle differs from $\\pi$ by more than the numerical error, the geometric-phase prediction is wrong. A complementary experimental falsifier: a Josephson junction starting from $\\Delta^L_k = -i\\Delta^R_k$ that cycles the right superconductor's parameters around the Dirac point and returns the current to its initial direction rather than reversing it would refute the claimed topological phase.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Landau free-energy expansion in powers of the order parameter that the paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the time-dependent relaxation dynamics of the order parameter near a phase transition."},{"cited_title":"Schmid, A time dependent Ginzburg-Landau equation and its application to the problem of resistivity in the 6 mixed state, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the time-dependent Ginzburg-Landau equation used to derive the adiabatic phase dynamics."}],"review_version":1}