{"id":"7b2bec88-567a-4950-9648-0e83ac849705","arxiv_id":"2506.07819","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Honeycomb lattices of phenalenyl-tessellation molecules realize effective Dirac models with hopping parameters in positive integer ratios given by the number of inter-molecule connections.","lead":"A honeycomb lattice made from phenalenyl-tessellation molecules is shown to host low-energy Dirac bands that match a simple two-site effective model. The effective hopping strengths are set by integer connection counts between molecules, which tunes band gaps and Fermi velocities while allowing coexisting vacancy-localized zero modes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gapped-case effective model relies on a fitted tau, so the abstract's 'determined only by connections' holds only for hopping ratios, not the Hamiltonian's overall scale; quantitative predictions in Sec. III.C are not first-principles.","rationale":"The paper's central claim is that the low-energy Dirac bands of H-PTMs are captured by an effective honeycomb model whose parameters are set by the number of inter-PTM connections. The supporting derivation is exact at gap-closing points, where the sqrt3 zero-mode ansatz is valid, and the numerical comparisons for isotropic and critical cases are convincing, including the analytically derived tau values in the Appendix. The soft spot is the gapped case: there is no zero-energy eigenstate to normalize, so tau is fitted rather than derived. This matters because the abstract and conclusion assert that the model is 'determined only by the connections between neighboring effective atoms.' If that were true, tau would be fixed by the PTM connectivity alone; but the paper itself shows tau changes between structures with identical connection ratios (Fig. 4 vs Fig. 5), so the overall energy scale is a property of the individual PTM, not just the connections. The integer-ratio design principle, and hence the gap-opening/closing criterion, survives because tau cancels in the triangle inequality; the concrete test proposed would establish whether the gapped effective model can be made predictive without an ad hoc fit. This is essentially the reader's weakest assumption, so I agree with the CONDITIONAL verdict and recommend no change.","tokens_in":13672,"tokens_out":7567,"duration_ms":93131,"concrete_test":"Use the exact TB eigenvector at the band-edge k-point of the 3:1:1 H-PTM (where |E| is minimal) to compute tau without fitting: project the full low-energy subspace onto the sqrt3 zero-mode basis and evaluate the effective hopping matrix element. Compare the resulting tau and the predicted band gap/bandwidth with the fitted tau = 0.0172t and with the original TB bands. If the derived tau differs by more than a few percent or the effective band fails to match the original within ~5% across the Brillouin zone, the gapped-case effective model is not determined by connections alone and the abstract should be amended.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation of Eq. (8) assumes the low-energy eigenstates are the k-independent sqrt3 x sqrt3 zero modes of each PTM, so tau = 2|phi_alpha||phi_beta|t is obtained from zero-mode normalization at a gap-closing point (K for isotropic, M for critical). In the gap-opening case (Sec. III.C, ND0:ND1:ND2 = 3:1:1), no zero-energy state exists; the authors state tau = 0.0172t is 'numerically determined by averaging |phi_alpha||phi_beta| across all DZCs.' Thus the effective Hamiltonian is not determined only by connections: its overall scale, and therefore the gap magnitude and bandwidth, depend on a fitted parameter that varies with PTM internal structure. For example, Fig. 4 gives tau = 2/11 t while Fig. 5 gives tau = 6/sqrt(4843) t for the same 1:1:1 connection ratios. The integer-ratio prediction of gap opening/closing is robust because tau cancels in the triangle inequality, but the abstract's stronger claim that the model is 'determined only by the connections' and any quantitative use of the gapped effective model are unsupported without a first-principles derivation of tau at E != 0. The derivation in Sec. II.D also contains an apparent index error in the displayed expansion (the inter-PTM term is written with psi_sqrt3_B1alpha instead of psi_sqrt3_B1beta), further obscuring the argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the low-energy electronic bands of a periodic honeycomb arrangement of phenalenyl-tessellation molecules (PTMs) can be described by a two-site effective honeycomb tight-binding model whose hopping ratios are given by the integers N_D0:N_D1:N_D2, the numbers of inter-PTM double-zigzag-corner connections in the three bond directions. The authors derive this effective model from the √3×√3 zero modes of the constituent PTMs, compare it with full tight-binding calculations for isotropic, critical, and gap-opening connection patterns, discuss Fermi-velocity control as a function of arm length, and show that vacancies introduce localized zero modes that coexist with the Dirac bands.","tokens_in":13972,"tokens_out":8048,"duration_ms":83193,"significance":"If the central claim holds, the paper offers a concrete bottom-up design principle for graphene-based lattices: varying the connection counts between molecular building blocks quantizes the effective hopping ratios, enabling controlled gap opening and Fermi-velocity modulation without large strain. The zero-mode normalization is derived analytically for several isotropic and critical structures, and the qualitative gap-opening/closing criterion is robust because it depends only on the integer ratios, not on the overall scale τ. The coexistence of Dirac bands with vacancy-localized spins extends the potential for correlated-electron devices. However, the claim that the effective model is 'determined only by the connections' is too strong, since the overall hopping scale τ=2|φ_α||φ_β|t depends on the internal zero-mode structure, and in the gap-opening case τ is obtained by numerical fitting rather than analytically. These caveats do not undermine the integer-ratio result itself, which is the paper's central and most valuable contribution.","major_comments":[{"comment":"The statement that the effective model is 'determined only by the connections between neighboring effective atoms' is not supported by the derivation. Equation (9) defines τ=2|φ_α||φ_β|t, and the normalization factors depend on the internal structure of the PTMs: for the same 1:1:1 connection counts, Section III.A reports τ=2/11t for Fig. 4 but τ=6/√4843t for Fig. 5. Only the ratios τ_0:τ_1:τ_2=N_D0:N_D1:N_D2 are fixed by the connection counts; the overall energy scale is not. The abstract and Section V should be reworded to state that the hopping ratios (not the full model) are determined by the connections.","section":"Abstract, Section II.D, Eq. (8)"},{"comment":"In the gap-opening case, the √3×√3 zero-mode assumption fails because the eigenenergy is not zero, and the value τ=0.0172t is 'numerically determined by averaging the values of |φ_α||φ_β| across all DZCs.' Consequently the quantitative band structure and gap magnitude of the effective model in this section rest on a fitted parameter rather than a first-principles derivation. The qualitative gap-opening prediction is robust because the triangle inequality involves only the integer ratios, but the paper should state this limitation explicitly and temper the claim of a parameter-free validation.","section":"Section III.C and Fig. 7"},{"comment":"The replacement of the four terms by the √3-mode expression contains an index error: the first term is written as ψ√3*_A1α H_AB ψ√3_B1α and the third as ψ*_A1β H_AB ψ√3_B1α, whereas the preceding expansion shows that the inter-PTM term is ψ*_A1α H_AB ψ_B1β and the intra-β term is ψ*_A1β H_AB ψ_B1β. The correct simplified expression should read ψ√3*_A1α H_AB ψ√3_B1β + ψ√3*_A1α H_AB ψ_B1α + ψ*_A1β H_AB ψ√3_B1β. As printed, the derivation of Eq. (8) is obscured.","section":"Section II.D, displayed expansion after Eq. (7)"}],"minor_comments":[{"comment":"The value τ=6/√4863t should read 6/√4843t; with |φ_α|=3/√167 and |φ_β|=1/√29, the product is 6/√(167·29)=6/√4843. The main text uses the correct value.","section":"Appendix"},{"comment":"The fractions are typeset ambiguously (e.g., 'τ= 2 8NPU + 3t' should be τ=2t/(8N_PU+3)); please use clear fraction notation throughout Section IV.A.","section":"Equations (11)-(19)"},{"comment":"Reference [35] appears to be a duplicate of [34] or an erratum with changed pagination; please verify the citation and list it as an erratum if that is the case.","section":"References [34] and [35]"},{"comment":"In Fig. 10(c), the red-dashed effective-model lines are described as unchanged from Fig. 5(e), but the figure appears to show only the original model; please indicate the effective-model bands more clearly or state that they are omitted from the panel.","section":"Figure 10"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central design idea is interesting. The main concern is the overstatement in the abstract; the authors should also double-check the index error in Sec. II.D, which I believe is a typo rather than a conceptual flaw. The fitted τ in Sec. III.C should be presented as a limitation. No concerns about novelty or citation integrity."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Morishita et al. show that in honeycomb lattices of phenalenyl-tessellation molecules, the effective hopping parameters between neighboring PTMs are positive integer multiples of a single scale tau, with the integers equal to the number of inter-PTM corner connections. That is a clean design rule, and the paper backs it up with analytic derivations of tau for the isotropic and critical cases (2/11 t, 9/419 t, 6/sqrt(4843) t), all matching full tight-binding bands. The Fermi-velocity formulas for arm-extended PTMs, with limits of 1/2 and anisotropic 0/1 relative to graphene, are a useful bonus. They also posted the tight-binding calculation files, which makes the numerical work checkable.\n\nThe main weakness is the gap-opening case in Sec. III.C. There the wavefunction is no longer the sqrt3 x sqrt3 zero mode, so tau = 0.0172t is fitted numerically from the full model. That means the effective Hamiltonian is not fully predictive in the gapped regime—only the hopping ratios are set by connections, while the overall scale (hence gap size and bandwidth) still depends on molecular details. The abstract's 'determined only by the connections' is too strong; it should say 'determined only by connections up to an overall structure-dependent scale.' The qualitative gap opening is robust because tau cancels in the triangle inequality, so the design principle survives, but the quantitative band structure in gapped systems is not a first-principles prediction.\n\nThere is also an apparent index error in the derivation in Sec. II.D: in the expansion after Eq. (7), the third term shows psi_A1beta H_AB psi_sqrt3_B1alpha, which looks like a typo for B1beta. Worth fixing.\n\nNone of this is fatal. The central argument—that connection counts give integer hopping ratios and control gap closure—is sound and well tested for the gapless cases. The paper deserves a serious referee and will probably need minor revision to fix the overstatement and the gapped-case caveat.","headline":"Nice integer-ratio effective model for PTM honeycombs, but the abstract overstates 'connections-only' and the gapped case relies on a fitted tau; worth a serious look.","tokens_in":14493,"tokens_out":3609,"would_cite":true,"duration_ms":38327,"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":"This paper establishes a structural design rule: in a honeycomb lattice of phenalenyl-tessellation molecules, the low-energy Dirac bands are determined solely by the numbers of corner connections between neighboring molecules, so the…","keywords":["graphene","Dirac cones","phenalenyl-tessellation molecules","tight-binding model","zero modes","Fermi velocity","nanographene","band gap engineering"],"falsifier":"Run a full tight-binding calculation for an H-PTM with connection counts 2:1:1, where the effective model predicts a gapless Dirac band, and for 3:1:1, where it predicts a gap; the central claim fails if the exact spectrum shows a gap in the first case or no gap in the second, or if the gap-closing point shifts away from the effective model's predicted k-point (the M point for 2:1:1).","tokens_in":13489,"feed_emoji":"🧩","tokens_out":9274,"duration_ms":102543,"temperature":0.7,"pith_summary":"This paper tries to establish a design rule for a family of nanographene lattices: in a honeycomb arrangement of phenalenyl-tessellation molecules (PTMs), the low-energy Dirac bands are governed entirely by how many double-zigzag corners connect each molecule to its neighbors. That rule makes the effective hopping parameters appear as positive integers in fixed ratios, so by choosing the connection counts a designer can decide whether the system is a gapless Dirac metal or a gapped semiconductor. It would matter because it turns a molecular-geometry choice into quantitative electronic design, including tuning the Fermi velocity from near zero up to values comparable to graphene, and it allows vacancy-localized zero modes to coexist with the Dirac electrons. A sympathetic reader would care because this is a parameter-free route, up to one overall scale, to modulated Dirac physics in a carbon-based platform.","feed_headline":"Nanographene Dirac bands tuned by integer link counts","feed_subtitle":"Each molecule acts as an atom; corner-link counts set hopping ratios, gap, and Fermi velocity.","key_machinery":"The carrying object is the honeycomb PTM lattice (H-PTM), in which each phenalenyl-tessellation molecule acts as an effective atom and neighboring molecules are joined through double-zigzag corners (DZCs). The identity that does the work is Eq. (8), which maps the many-site bipartite lattice onto the textbook two-site honeycomb tight-binding model with effective hoppings $\\tau_i = N_{Di}\\tau$. The derivation rests on the zero-sum rule for bipartite zero modes: for the $\\sqrt{3}\\times\\sqrt{3}$ zero-mode wavefunctions on $\\alpha$- and $\\beta$-PTMs, all inter-PTM coupling terms vanish except those that touch the two sublattices through DZCs, leaving only the connection numbers $N_{Di}$. That is what converts molecular geometry into counting, and it is also what fixes $\\tau$ analytically for simple structures through the normalization of the zero-mode amplitudes.","core_discovery":"The central claim is that a two-site effective honeycomb model describes the low-energy Dirac bands of an H-PTM lattice exactly enough to decide gap opening and closing. The effective Hamiltonian is $H_{\\mathrm{eff}} = \\tau$ times the standard honeycomb matrix with entries $N_{D0} + N_{D1} e^{-i\\mathbf{k}\\cdot\\mathbf{a}'_1} + N_{D2} e^{-i\\mathbf{k}\\cdot\\mathbf{a}'_2}$, where $N_{D0}$, $N_{D1}$, $N_{D2}$ are the numbers of double-zigzag-corner connections in the three bond directions and $\\tau = 2|\\varphi_\\alpha||\\varphi_\\beta|t$ is set by the amplitudes of the $\\sqrt{3}\\times\\sqrt{3}$ zero modes on the two molecular sublattices. Since $\\tau$ is just a common scale, the three effective hoppings have positive integer ratios $N_{D0} : N_{D1} : N_{D2}$, and the triangular inequality for those integers decides whether the bands touch in Dirac cones, merge at a critical point, or open a gap. The paper verifies this against full tight-binding spectra for isotropic (1:1:1), critical (2:1:1), and gapped (3:1:1) cases, and shows that vacancies can introduce flat zero modes at $E=0$ without removing the Dirac cone.","pith_inferences":["Editorial inference: the integer-count rule suggests a combinatorial design chart, since every triple $(N_{D0}, N_{D1}, N_{D2})$ realizes a distinct anisotropic honeycomb model, including Dirac-point-merge and gap-opening phases, so one could in principle enumerate all achievable band structures.","Editorial inference: the coexistence of vacancy zero modes with a tunable Dirac band points toward systems where localized spins and itinerant electrons interact magnetically, a direction the paper mentions but does not compute.","Editorial inference: a testable extension would be to break the symmetry between the alpha and beta molecules, for example by removing atoms from only one type, to see whether the integer-hopping effective model predicts an energy gap at the Dirac point while the zero-sum derivation still holds.","Editorial inference: because the effective model depends only on connection counts, small energy shifts on interior carbon atoms that preserve the molecular tiling should leave the Dirac physics unchanged; this insensitivity could be checked by adding random small site-energy variations in tight-binding calculations."],"forward_implications":["The low-energy band character of an H-PTM is fixed by three integers, $N_{D0}$, $N_{D1}$, $N_{D2}$, and one overall scale $\\tau$; no other molecular detail matters for whether a gap opens.","Choosing integers that violate the triangle inequality opens a band gap without strain, and choosing equality places the Dirac-point merge at a controlled k-point such as the M point.","The Fermi velocity can be engineered from near zero (uniaxial lattices with long arms) up to a value comparable to graphene's (short isotropic arms), with explicit $1/L_{\\mathrm{arm}}$ scaling laws.","Vacancies inside PTM regions create localized zero modes at the Fermi level that coexist with the Dirac cone, giving a built-in coexistence of localized spins and itinerant electrons.","Because the hoppings are integer multiples of a common $\\tau$, the effective model is parameter-free apart from $\\tau$, making the design rules directly checkable by tight-binding calculation."],"supporting_citations":[{"why":"Defines the PTM construction and gives the topology-based counting of zero modes that underlies the effective-atom picture.","marker":"[20]"},{"why":"Establishes the sqrt3 x sqrt3 zero-mode shape and the zero-sum rule on PTM networks used to derive Eq. (8).","marker":"[21]"},{"why":"Supplies the tight-binding zero-mode analysis of honeycomb lattices and the triangular-inequality condition inherited by the effective model.","marker":"[28]"},{"why":"Shows how anisotropic hoppings on generalized honeycomb lattices control Dirac cones and gap opening, the behavior the effective model reproduces.","marker":"[30]"},{"why":"Provides the zero-sum rule for zero-energy eigenvectors used to force the uniform sqrt3 x sqrt3 amplitude pattern.","marker":"[38]"},{"why":"Fixes graphene's tight-binding parameter t and Dirac-cone form that serve as the comparison baseline for Fermi-velocity ratios.","marker":"[37]"},{"why":"Shows that vacancies on PTMs give localized zero modes with spin character, supporting the coexistence claim.","marker":"[22]"},{"why":"Supplies the multi-vacancy N=M=3 case in which six localized zero modes appear without disturbing the effective Dirac model.","marker":"[36]"}],"fun_headline_variants":["Integer bond ratios tune Dirac gap and velocity","Molecule-as-atom honeycomb: integer hops set Dirac bands","Nanographene Dirac bands from integer connection counts","Effective honeycomb model: integer ratios control Dirac cones","Phenalenyl tessellations: integer links dictate Dirac dispersion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each molecule's low-energy state keeps the uniform $\\sqrt{3}\\times\\sqrt{3}$ zero-mode shape for every wavevector that matters, which is exact only at the gap-closing point; once a gap opens the shape is modulated and the effective hopping must be fitted numerically rather than derived.","fun_headline_variants_meta":{"raw":{"variants":["Integer bond ratios tune Dirac gap and velocity","Molecule-as-atom honeycomb: integer hops set Dirac bands","Nanographene Dirac bands from integer connection counts","Effective honeycomb model: integer ratios control Dirac cones","Phenalenyl tessellations: integer links dictate Dirac dispersion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000358,"raw_usage":{"total_tokens":2003,"prompt_tokens":1069,"completion_tokens":934,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":856}},"tokens_in":685,"tokens_out":934,"duration_ms":10484,"temperature":1.0,"reasoning_tokens":856,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:24:01.531824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full tight-binding calculation for an H-PTM with connection counts 2:1:1, where the effective model predicts a gapless Dirac band, and for 3:1:1, where it predicts a gap; the central claim fails if the exact spectrum shows a gap in the first case or no gap in the second, or if the gap-closing point shifts away from the effective model's predicted k-point (the M point for 2:1:1).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the PTM construction and gives the topology-based counting of zero modes that underlies the effective-atom picture."},{"cited_title":"Ziatdinov, S","cited_arxiv_id":null,"evidence_quote":"Establishes the sqrt3 x sqrt3 zero-mode shape and the zero-sum rule on PTM networks used to derive Eq. (8)."},{"cited_title":"Hasegawa, R","cited_arxiv_id":null,"evidence_quote":"Supplies the tight-binding zero-mode analysis of honeycomb lattices and the triangular-inequality condition inherited by the effective model."},{"cited_title":"Ortiz, G","cited_arxiv_id":null,"evidence_quote":"Shows how anisotropic hoppings on generalized honeycomb lattices control Dirac cones and gap opening, the behavior the effective model reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the zero-sum rule for zero-energy eigenvectors used to force the uniform sqrt3 x sqrt3 amplitude pattern."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fixes graphene's tight-binding parameter t and Dirac-cone form that serve as the comparison baseline for Fermi-velocity ratios."},{"cited_title":"Morishita, Y","cited_arxiv_id":null,"evidence_quote":"Shows that vacancies on PTMs give localized zero modes with spin character, supporting the coexistence claim."}],"review_version":1}