{"id":"aff179fd-5af1-49b5-9a0c-bb6068551a67","arxiv_id":"1908.05092","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"All eight non-trivial Archimedean lattices host quantum spin Hall phases at some band fillings, together with flat bands, Dirac cones, and high-degeneracy points, according to tight-binding and DFT calculations.","lead":"This paper maps the electronic and topological properties of eight Archimedean lattices using tight-binding models, finding flat bands, Dirac fermions, and quantum spin Hall phases at certain fillings. A generalist should read it because it offers a systematic guide for designing 2D materials and photonic or phononic analogs with targeted band structures.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The intrinsic-SOC term in Eq. (3) is underspecified: the intermediate index k is neither summed nor defined, so the Fig. 4 Z2 map may depend on an arbitrary convention in lattices with multiple two-step paths.","rationale":"I read the paper as a systematic tight-binding survey whose central contribution is the claim that each non-regular Archimedean lattice can host a QSH phase at some occupation once intrinsic SOC is added. That claim is a model-level statement, so the most decisive weakness is not the absence of an ab initio lambda (which affects material realization) but whether the Hamiltonian generating the Z2 map is actually defined. Eq. (3) is the only definition of the SOC term, and it leaves the intermediate index k unresolved. In lattices with triangular motifs this is not a cosmetic issue: a two-step path is not unique, and the different choices can give opposite signs or cancel. The paper's DFT validation of pz bands without SOC cannot resolve this, because it never exercises H_SO. I therefore flag this as the load-bearing concern. I want to credit the authors for a standard Z2 calculation (Wannier charge centers) and for stable, phonon-checked DFT structures; those parts are solid. The concern is fixable by specifying the H_SO convention and, ideally, by showing the Z2 map for both conventions or a test lattice. With that clarification the model claim could stand; as written it is conditional. This is consistent with the reader's CONDITIONAL verdict, though for a more internal reason than the reader's emphasis on unspecified lambda and Rashba.","tokens_in":11978,"tokens_out":14646,"duration_ms":159650,"concrete_test":"Reimplement Eqs. (2)-(5) for the kagome lattice (3,6,3,6) with alpha=3.0/d_nn, t=1, lambda=1. For each second-neighbor pair (i,j), enumerate all intermediate sites k that are nearest neighbors of both i and j, and evaluate H_SO in two ways: (A) sum over all such k, and (B) choose exactly one k per pair following a crystallographically consistent orientation rule. Compute the band structures and the Z2 invariant at the occupations marked in Fig. 4 under both conventions. If convention (A) yields zero or different Z2 values, or if (A) and (B) disagree on any Z2=1 entry, the paper must disclose which convention was used and the central claim needs qualification; if both conventions reproduce the same Fig. 4 entries, the ambiguity is not numerically load-bearing.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Equation (3) defines H_SO = i \\sum_{i,j} \\lambda_{ij} c_i^\\dagger \\sigma \\cdot (d_{kj} \\times d_{ik}) c_j, but the index k is never defined or included in the sum. In the standard Kane-Mele construction k is the unique intermediate atom of a second-neighbor pair; this is unambiguous for honeycomb, but many of the eight studied lattices (e.g., kagome (3,6,3,6), (3^4,6), (3,4,6,4), (3^3,4^2)) contain triangular cycles in which a pair i,j has two common nearest-neighbor sites. If one sums over all such k, the cross-product contributions from a triangular plaquette have opposite signs and can cancel, possibly making H_SO vanish; if one instead selects one k per pair, the resulting Hamiltonian depends on an arbitrary orientation rule that is not stated and that may break the lattice symmetry. The central claim, that all studied lattices present a QSH phase, is computed from this Hamiltonian (Fig. 3 and Fig. 4 at alpha=3.0/d_nn), yet the paper neither specifies the k-sum convention nor gives a numerical lambda. As written, the Z2 map is therefore not uniquely defined; the model must state how H_SO is evaluated for multi-path pairs before the QSH claim can be assessed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies the electronic structure of eight non-trivial Archimedean lattices using a single-orbital tight-binding model with exponentially decaying hoppings and an intrinsic spin-orbit coupling (SOC) term. The authors classify the zero-SOC band structures, identifying Dirac points, flat bands, partially flat bands, and high-degeneracy points, and they compute the Z2 topological invariant by Wannier charge centers for each lattice as a function of occupied Kramers pairs of bands. The central claim is that every studied Archimedean lattice hosts a quantum spin Hall (QSH) phase at some particular occupation, signaled by spin-textured edge states in nanoribbon calculations. The paper also presents PBE-DFT relaxations, phonon dispersions, and pz-projected band structures for three planar carbon allotropes with the (3,12^2), (4,6,12), and (4,8^2) lattices as a material realization of the non-SOC tight-binding band features.","tokens_in":12270,"tokens_out":13226,"duration_ms":142017,"significance":"If the central claim is correct, the paper provides a useful systematic catalog of single-particle band features and topological phases across the Archimedean lattices, with potential applications to carbon allotropes, metal-organic frameworks, and photonic and phononic analogs. The non-SOC band-structure analysis and the DFT phonon-stability calculations for three carbon allotropes are careful and use standard, reproducible methods; the Wannier-charge-center approach for Z2 is appropriate. The significance is currently limited because the SOC Hamiltonian in Eq. (3) is not uniquely defined for lattices with multiple two-step paths and because no numerical SOC strength is specified, so the central QSH prediction is not yet a fully quantitative, falsifiable statement. These issues are local and can be addressed in revision.","major_comments":[{"comment":"The intrinsic SOC term is not uniquely defined as written. The intermediate site k in d_{kj} × d_{ik} is never introduced or summed over, and the double sum over i,j cannot determine k. In the honeycomb lattice each second-neighbor pair has a unique common nearest neighbor, which is the standard Kane-Mele convention, but the Archimedean lattices studied here contain square and other small cycles (e.g., (4,8^2), (3,4,6,4), (3^2,4,3,4), (3^3,4^2)) in which a pair of sites can be connected by two different two-step paths. If all intermediate k are summed, the two paths around a square give opposite cross-product contributions and may cancel; if one path is selected, the rule is unspecified and may break lattice symmetry. Because the Z2 map in Fig. 4 and the edge states in Fig. 3 are computed from this Hamiltonian, the central claim depends on an arbitrary convention. Please specify the convention, prove that the resulting Hamiltonian is Hermitian and symmetry-preserving, and verify that the Z2 assignments are independent of the choice.","section":"II, Eq. (3)"},{"comment":"The SOC strength λ is never assigned a numerical value and the normalization Nλ in Eq. (5) is not defined, so the gap sizes, edge-state dispersions, and the physical regime of the QSH prediction are unspecified. The Z2 invariant is insensitive to the magnitude of λ as long as the relevant gap remains open, but the statement that all studied lattices present a QSH phase in some particular occupation is computed at a single value of the hopping decay parameter, α = 3.0/d_nn, and at an unstated λ. Please state the λ/t ratio used for Figs. 3 and 4, define Nλ, and show that the Z2 occupations are stable over a reasonable range of α and λ, or explicitly limit the claim to the chosen parameter set.","section":"III.B, Figs. 3 and 4"}],"minor_comments":[{"comment":"The statement that the discussion is \"validated within density functional theory calculations\" is too strong: the DFT calculations in Section III.C and Fig. 5 are performed without SOC and validate only the pz-derived, non-SOC band features of three carbon allotropes. Please qualify the claim accordingly.","section":"Abstract and III.C"},{"comment":"The title advertises a \"topological flat band,\" but no Chern number or other topological characterization of the flat bands is computed; consider rephrasing to \"flat bands, Dirac fermions, and quantum spin Hall phase.\"","section":"Title"},{"comment":"Please specify the nanoribbon geometry used for the edge-state calculations, including the edge termination and width, so that the T1 and T2 edge states can be reproduced.","section":"III.B, Fig. 3"},{"comment":"There are several typographical and grammatical errors, including \"arises\" for \"arise\" (Abstract), \"presents\" for \"present\" (Section III.B), \"Bellow\" for \"Below\" (Section III.C), \"backscaterring\" for \"backscattering\" (Fig. 3 caption), and \"V ASP\" for \"VASP\" (Section III.C).","section":"Throughout"},{"comment":"References 11 and 40 are the same Kane-Mele paper and should be consolidated.","section":"References"},{"comment":"The caption states that the occupation number indicates \"the number of occupied Kramers pairs bands\"; this should read \"Kramers pairs of bands.\"","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a 2019 arXiv preprint, and the authors may need to update the reference list with subsequent work on Archimedean lattices. The main obstacle to publication is the underspecified SOC Hamiltonian in Eq. (3); once the convention is stated and the parameter values are provided, the work would be a useful catalog of flat-band, Dirac, and QSH physics in these lattices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious referee if the authors fix one ambiguity. The genuinely new thing is a single-orbital tight-binding survey of all eight non-trivial Archimedean lattices, with a Z2 invariant map across fillings, edge states, flat bands, and Dirac crossings. That systematic treatment is not in the prior literature, and it gives a practical design guide for photonic, phononic, and magnonic analogs. The DFT section is a solid proof of principle for three carbon allotropes: phonons are clean, and the pz-derived bands match the TB picture. The main problem is Eq. (3). The index k in the spin-orbit term is not defined and not summed over. In honeycomb, each second-neighbor pair has a unique intermediate atom, so Kane-Mele is unambiguous. In several of these lattices there are multiple two-step paths between a given pair of sites (e.g., diagonal pairs in squares or other cycles). Without a stated convention for choosing k, H_SO is not uniquely defined, and the Fig. 4 Z2 map may depend on that convention. This is not a cosmetic detail: summing over all k can give cancellation, while picking one k per pair introduces an arbitrary orientation rule that may break lattice symmetry. The authors need to state the convention explicitly and show that the Z2 results are robust to alternatives. Minor issues: the abstract says 'all Archimedean lattices' but only the eight non-trivial ones are treated, with the three standard lattices excluded. Lambda is never given a numerical value, though for a small intrinsic SOC the Z2 classification likely only requires it to be non-zero and small. The DFT checks are done without SOC, so they do not test the QSH claim directly; they validate the non-interacting band structure, which is fine for a proof of principle but leaves the material-level QSH phases as predictions rather than confirmations. Overall, the paper is a clean classification exercise using standard methods. If the SOC Hamiltonian is properly specified, it belongs in the literature as a reference for Archimedean lattice physics. My recommendation: send it to peer review, but require the authors to clarify Eq. (3) and to recompute the Z2 invariant for at least one multi-path lattice under both plausible conventions to demonstrate robustness.","headline":"Useful systematic tight-binding catalog of the eight non-trivial Archimedean lattices, but the spin-orbit Hamiltonian in Eq. (3) is underspecified and needs fixing before the Z2 map is fully defined.","tokens_in":635,"tokens_out":1673,"would_cite":true,"duration_ms":44881,"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":"Every non-trivial Archimedean lattice hosts a quantum spin Hall phase at some filling.","keywords":["Archimedean lattices","quantum spin Hall effect","tight-binding model","spin-orbit coupling","Dirac fermions","flat bands","Z2 topological invariant","two-dimensional carbon allotropes"],"falsifier":"Repeat the $\\mathbb{Z}_2$ calculation of Fig. 4 for the eight lattices with nearest-neighbour-only hopping ($\\alpha = 20/d_{nn}$) and check whether at least one occupation still gives $\\mathbb{Z}_2 = 1$; the paper only reports the $\\mathbb{Z}_2$ map at $\\alpha = 3/d_{nn}$, so a lattice that loses its topological phase at $\\alpha = 20/d_{nn}$ would show the result depends on the hopping range. Alternatively, compute the intrinsic spin-orbit coupling strength for the $(4,8^2)$ carbon allotrope from first principles and see whether the predicted edge states survive.","tokens_in":11782,"feed_emoji":"⚛️","tokens_out":7499,"duration_ms":68385,"temperature":0.7,"pith_summary":"The paper claims that each of the eight non-trivial Archimedean lattices—the uniform tilings of the plane beyond the triangular, square, and honeycomb cases—becomes a quantum spin Hall insulator at some electronic filling once intrinsic spin-orbit coupling is turned on. It maps which occupations of Kramers pairs give $\\mathbb{Z}_2 = 1$, showing backscattering-protected helical edge states in nanoribbon calculations. Along the way it catalogues type-I and type-II Dirac fermions, flat bands, and high-degeneracy points in these lattices, and validates the tight-binding picture with density functional theory on three stable carbon allotropes. If correct, any material that realizes one of these tilings is a candidate topological insulator at a predictable filling.","feed_headline":"Every Archimedean lattice hosts a quantum spin Hall phase","feed_subtitle":"Each of eight non-trivial tilings has a filling with Z2=1; carbon allotropes confirm the bands.","key_machinery":"The central object is the tight-binding Hamiltonian $H_{TB} = H_0 + H_{SO}$ with one orbital per site, nearest-neighbour hopping $t = 1$ as the energy unit, and hopping and spin-orbit amplitudes both decaying as $\\exp(-\\alpha d)$ with inter-site distance, normalized at the nearest-neighbour distance. Because the lattices are mirror symmetric in the plane, the Rashba term is set to zero, leaving intrinsic spin-orbit coupling as the gap-opening mechanism. The topological verdict is delivered by the $\\mathbb{Z}_2$ invariant computed from the evolution of Wannier charge centers, and the same Hamiltonian is interpreted as applying to photonic, phononic, and magnonic systems as well as electrons.","core_discovery":"Within a single-orbital tight-binding model with hoppings and intrinsic spin-orbit coupling that decay exponentially with distance, the paper finds that all eight non-trivial Archimedean lattices have at least one Kramers-pair band occupation with $\\mathbb{Z}_2 = 1$, meaning a quantum spin Hall phase with helical, backscattering-protected edge states. The topological phases appear when spin-orbit coupling opens gaps at Dirac crossings and at degenerate flat-band points; the paper computes the $\\mathbb{Z}_2$ invariant by tracking Wannier charge centers and maps it as a function of band filling. For the $(3,12^2)$, $(4,6,12)$, and $(4,8^2)$ carbon allotropes, the $p_z$-derived bands from density functional theory reproduce the Dirac crossings and flat bands of the model, supporting the claim that the model captures real materials.","pith_inferences":["The paper stops short of computing realistic spin-orbit coupling strengths; from first principles, intrinsic spin-orbit coupling in carbon is weak, so the quantum spin Hall phase may require proximity-induced spin-orbit coupling or heavier-element versions of these lattices to be observable.","Because the density functional theory check tests only $p_z$ bands without spin-orbit coupling, the model's zero-Rashba assumption may break in the low-symmetry oblique $(3^3,4^2)$ lattice or on substrates; a buckling calculation would settle this.","The accidental degeneracy in $(4,6,12)$ that lacks a defined $\\mathbb{Z}_2$ invariant could be tuned by strain or superlattice design to open a topological gap, turning a currently undefined point into a phase transition."],"forward_implications":["For each of the eight non-trivial Archimedean lattices, the $\\mathbb{Z}_2$ map identifies specific Kramers-pair occupations at which the system is a quantum spin Hall insulator, so a material with that lattice and filling should show helical edge states.","The same single-orbital model applies to photonic, phononic, and magnonic realizations, so the predicted gaps and edge modes can be probed in classical-wave lattices without electrons.","The three carbon allotropes $(3,12^2)$, $(4,6,12)$, and $(4,8^2)$ are dynamically stable and their $p_z$-projected bands reproduce the Dirac and flat-band features, making them concrete candidates for the predicted phases once spin-orbit coupling is introduced by proximity.","The type-II Dirac crossing in $(3^3,4^2)$ and the pseudospin-1 and pseudospin-2 degeneracies in $(4,8^2)$ and $(3^4,6)$ provide specific band features to search for in photoemission or transport experiments."],"supporting_citations":[{"why":"Defines the Archimedean tiling classification that selects the eight non-trivial lattices under study.","marker":"[25]"},{"why":"Provides the intrinsic spin-orbit Hamiltonian and the quantum spin Hall criterion used to model the phases.","marker":"[40]"},{"why":"Gives the Wannier charge-center expression for the $\\mathbb{Z}_2$ invariant used to map the phases.","marker":"[51]"},{"why":"Supplies the method to compute $\\mathbb{Z}_2$ without inversion symmetry, needed for the oblique lattice.","marker":"[52]"},{"why":"Proposes graphenylene, the $(4,6,12)$ carbon allotrope used as a material realization.","marker":"[29]"},{"why":"Reports the stability of planar carbon sheets, supporting the phonon validation of the density functional theory structures.","marker":"[61]"},{"why":"Provides the formation-energy scale of graphyne used to argue the carbon allotropes are experimentally feasible.","marker":"[59]"}],"fun_headline_variants":["Quantum spin Hall phase found in every Archimedean lattice","All Archimedean lattices predicted to host topological insulators","Carbon allotropes confirm spin Hall phase in all Archimedean tilings","Every Archimedean tiling yields a topological insulator","Universal quantum spin Hall effect in Archimedean lattices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole $\\mathbb{Z}_2$ map depends on the single-orbital tight-binding model with exponential hopping and intrinsic spin-orbit coupling and zero Rashba terms; if a real Archimedean material has significant multi-orbital mixing, substrate-induced Rashba, or a different hopping decay, the predicted topological gaps and edge states may not survive.","fun_headline_variants_meta":{"raw":{"variants":["Quantum spin Hall phase found in every Archimedean lattice","All Archimedean lattices predicted to host topological insulators","Carbon allotropes confirm spin Hall phase in all Archimedean tilings","Every Archimedean tiling yields a topological insulator","Universal quantum spin Hall effect in Archimedean lattices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000245,"raw_usage":{"total_tokens":1486,"prompt_tokens":849,"completion_tokens":637,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":465,"completion_tokens_details":{"reasoning_tokens":547}},"tokens_in":465,"tokens_out":637,"duration_ms":5412,"temperature":1.0,"reasoning_tokens":547,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:23:30.936354+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the $\\mathbb{Z}_2$ calculation of Fig. 4 for the eight lattices with nearest-neighbour-only hopping ($\\alpha = 20/d_{nn}$) and check whether at least one occupation still gives $\\mathbb{Z}_2 = 1$; the paper only reports the $\\mathbb{Z}_2$ map at $\\alpha = 3/d_{nn}$, so a lattice that loses its topological phase at $\\alpha = 20/d_{nn}$ would show the result depends on the hopping range. Alternatively, compute the intrinsic spin-orbit coupling strength for the $(4,8^2)$ carbon allotrope from first principles and see whether the predicted edge states survive.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Archimedean tiling classification that selects the eight non-trivial lattices under study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Wannier charge-center expression for the $\\mathbb{Z}_2$ invariant used to map the phases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proposes graphenylene, the $(4,6,12)$ carbon allotrope used as a material realization."}],"review_version":1}