{"id":"2fd4c582-2c5b-494a-be6a-73e2089e6a5e","arxiv_id":"2501.06851","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A harmonic crystal monolayer on a decagonal quasicrystal is predicted to adopt a nonzero twist angle and a striped distortion pattern, confirmed by simulation.","lead":"A crystalline layer on a quasicrystalline surface is predicted to twist to a nonzero angle rather than align with the surface's fivefold axes. This creates a striped distortion pattern with strongly directional sliding and friction properties.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Constrained boundary conditions and a single length ratio leave the claimed generic nonzero twist unproven; free-boundary or strain-relaxed tests are needed.","rationale":"The core second-order phonon calculation leading to Eqs. (10) and (11) is standard, and the parameter-free numerical comparison at weak coupling is genuine evidence that the constrained calculation is internally consistent. The reader's weakest-assumption identification is broadly correct: the system is constrained to a uniform, rigidly rotating harmonic lattice by the variational ansatz and by the fixed outer ring in the numerics. I sharpen this into a concrete, testable concern by noting that homogeneous strain is also excluded, and that the 'generally' in the abstract is supported by only one length ratio in the reviewed text. These are not accusations of error; they are statements of what evidence would be needed to move from a conditional to an unconditional acceptance. The proposed free-boundary simulation and strained variational minimization would settle whether the predicted nonzero misfit angle and striped state survive when the layer is allowed the relaxation channels available in a real unconstrained monolayer. Until such a check is reported, the conditional verdict remains appropriate.","tokens_in":7692,"tokens_out":12228,"duration_ms":150209,"concrete_test":"Run LAMMPS simulations of the same model with open/free boundaries (no fixed ring), starting from several initial orientations and coupling strengths, and compare the spontaneously selected global orientation and the dominant Fourier component of the relaxed displacement field with Eq. (11) and θopt ≈ 5.31°. Independently, augment the variational ansatz of Eq. (10) with a homogeneous strain tensor ε, minimizing E1-ph over both θ and ε with q = G − R(θ)(1+ε)τ. If free-boundary relaxation or strained minimization selects θ ≈ 5.31° and a single q pair, the concern is resolved; if a different angle, a multiple-q pattern, or a domain state emerges, the constrained boundary is load-bearing and the generality claim is weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the equilibrium layer 'will generally exhibit a nonzero misfit angle' rests on two load-bearing supports that are not fully demonstrated. First, Eq. (10) is minimized only over the global rotation angle θ of the overlayer reciprocal lattice τ; the reference lattice itself is not allowed to undergo a homogeneous strain ε that would change the matching condition to q = G − R(θ)(1+ε)τ. In crystal-on-crystal Novaco–McTague epitaxy, such homogeneous strain is part of the mechanism, so excluding it may shift the optimum. Second, the numerical simulations fix the outermost ring at perfect-lattice positions, which imposes θ and suppresses bulk uniform dilation or shear, and the main text presents only one length ratio apot/acoll = 5.4 µm/5.8 µm. The Supplemental Material is said to contain another ratio, but it was not part of the reviewed version. The observed excellent agreement at g = 10⁻⁴ validates Eq. (10) for the constrained geometry, not the unconstrained equilibrium orientation. If a free monolayer can relax by uniform strain or spatially varying rotation, the predicted θopt ≈ 5.31° and the clean stripe pattern may be artifacts of the boundary constraint rather than the generic equilibrium state claimed in the abstract.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the Novaco-McTague weak-corrugation theory of orientational epitaxy from crystalline substrates to a decagonal quasicrystalline substrate. A harmonic hexagonal monolayer is coupled to a quasiperiodic potential generated by p=5 interfering beams. Using a coherent-state variational ansatz and a one-phonon approximation, the authors derive a closed-form energy lowering, Eq. (10), and the equilibrium displacement field, Eq. (11). Minimizing Eq. (10) over the global rotation angle θ for one length ratio a_pot/a_coll = 5.4/5.8 yields θ_opt ≈ 5.31°, with the distortion dominated by a single pair of shortest reciprocal vectors q_min, hence a striped moiré pattern. LAMMPS simulations of a circular sample with a fixed outermost ring at perfect-lattice positions reproduce the energy curve quantitatively for g = 10^-4 and qualitatively for larger g. The paper claims that a nonzero misfit angle is generic for crystal-on-quasicrystal epitaxy.","tokens_in":7989,"tokens_out":8616,"duration_ms":95358,"significance":"If the generic claim were established, this would be a natural and important extension of Novaco-McTague physics to quasicrystalline substrates, with a concrete, falsifiable striped state and anisotropic tribological consequences. The paper has real strengths: Eq. (10) is a parameter-free analytic formula, the small-g simulation comparison is quantitative, and the stripe mechanism is clearly identified through the dominance of a single q_min pair. The main limitations are that the numerical validation is a self-consistency check of the same model, the boundary conditions impose the global rotation and suppress homogeneous strain, and only one length ratio is examined. These limitations do not invalidate the formalism, but they leave the 'generally' claim of the abstract under-supported.","major_comments":[{"comment":"Eq. (10) is minimized only over the global rotation θ of the reference lattice, with the matching condition q = G − R(θ)τ. A homogeneous strain ε of the overlayer would change this condition to q = G − R(θ)(1+ε)τ and can lower the energy in epitaxial systems; the manuscript neither includes this degree of freedom nor proves it is negligible. Since the numerical simulations fix the outermost ring at perfect-lattice positions, they explicitly suppress uniform dilation and shear. The claim that the free equilibrium is θ_opt ≈ 5.31° therefore requires either an analytic argument excluding homogeneous strain or a free-boundary/affine-relaxation simulation.","section":"Variational approach, Eq. (10)"},{"comment":"The generic statement in the abstract is supported by a single length ratio a_pot/a_coll = 5.4/5.8. The text calls this a 'generic lattice-incommensurate situation', but no analytic proof or multi-ratio scan is provided; the second ratio is only mentioned in the Supplemental Material and is not part of the reviewed manuscript. Because the position of q_min(θ) and hence the sign and size of θ_opt can depend on the ratio, a reader cannot exclude that the nonzero θ_opt is specific to the chosen ratio. Please provide an argument covering generic irrational ratios or a systematic scan over many ratios.","section":"Numerical, Fig. 3"},{"comment":"The fixed outermost ring imposes the twist angle and prevents spatially nonuniform relaxation such as domain walls, dislocations, or a bulk rotation that differs from the boundary rotation. The displacement field in Eq. (11) is a coherent sinusoidal modulation of the whole lattice, and the stripe pattern in Fig. 4 is obtained under this constraint. Since the abstract claims the equilibrium configuration generically has nonzero misfit with stripes, free-boundary or periodic simulations that allow the orientation and strain to relax are needed to show that the pattern is not an artifact of the boundary control.","section":"Numerical, Fig. 4"}],"minor_comments":[{"comment":"The phrase 'quite different from the the ordinary moiré pattern' contains a duplicated definite article.","section":"Introduction, p. 2"},{"comment":"The author list 'D. Shechtman, I. Blech, and D. Gratias, J. W. Cahn Phys. Rev. Lett.' is missing a comma between Gratias and Cahn, making the attribution unclear.","section":"References, Ref. [33]"},{"comment":"The vertical dashed line is labeled 'θ_opt ≃ 5.31°'; please state explicitly whether this value is the analytic minimum of Eq. (10) or the numerical minimum, since the two are not necessarily identical at larger g.","section":"Fig. 3, p. 3"},{"comment":"The statement that the theory 'perfectly agrees' with numerical optimization is precise only in the weak-coupling regime; Fig. 3 shows quantitative agreement at g = 10^-4 but only qualitative agreement at g = 0.01-0.05. Please qualify the claim to avoid overstating the range of quantitative validity.","section":"Abstract and Discussion"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know that this paper gives the first clean extension of Novaco–McTague orientational epitaxy to a quasicrystalline substrate. The genuinely new result is that the optimal distortion is dominated by a single inversion-related pair of reciprocal vectors, producing stripes rather than the hexagonal moiré of crystal-on-crystal epitaxy. That is a real prediction, not a routine application, and the authors back it with a parameter-free comparison to LAMMPS simulations at weak coupling. The agreement at g = 10⁻⁴ in Fig. 3 is quantitatively excellent and the displacement patterns in Fig. 4 match Eq. (11) convincingly.\n\nThe derivation is standard but carefully executed: coherent-state variational energy, one-phonon approximation, momentum conservation giving Eq. (10). The identification of the shortest q-vector pair as the main contributor is physically transparent and well supported. Citation practice looks honest—the relevant NM, Shiba, and colloidal experimental literature is cited.\n\nNow the soft spots, in proportion. The abstract says the layer \"will generally exhibit a nonzero misfit angle,\" but the evidence is one length ratio apot/acoll in the main text and a second in the Supplemental Material, which was not part of the reviewed version. That makes the \"generally\" claim stronger than the shown proof. More substantively, the variational space only includes global rotation, not homogeneous strain. In NM theory for crystals, strain is part of the mechanism, so omitting it could shift the optimal angle. The simulation fixes the outermost ring at perfect-lattice positions, which imposes θ and suppresses uniform dilation or shear. So the excellent numerical agreement validates Eq. (10) for the constrained geometry, not necessarily an unconstrained equilibrium. A free monolayer that can relax via strain or spatially varying rotation might find a slightly different optimum; the clean stripe pattern could be smeared by domain walls or dislocations at larger couplings or different ratios. These are addressable concerns, not fatal flaws. The weak-coupling regime is where the prediction lives, and there the uniform-distortion ansatz is plausible.\n\nOne minor quibble: \"perfectly agrees\" is a bit strong for the high-g cases, where the agreement is qualitative. But that does not undercut the main result.\n\nWho is this for? Researchers working on 2D epitaxy, moiré materials, and colloids on optical lattices. It deserves a serious referee. The referee should ask for the SM, a second length ratio, and a discussion (or short simulation) addressing homogeneous strain and free boundaries. With those, the generality claim could become solid.\n\nI would send this to peer review.","headline":"First clean extension of Novaco-McTague theory to quasicrystal substrates, predicting a striped twisted state; the model-level result is convincing but the 'generic' claim is not fully established.","tokens_in":8516,"tokens_out":2103,"would_cite":true,"duration_ms":23268,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["68.35.Af","68.08.De","62.10.+s","62.20.Qp"],"model":"deepseek-v4-flash","headline":"The paper predicts that a crystalline monolayer on a quasicrystalline substrate will generally sit at a nonzero misfit angle, forming a striped, not hexagonal, moiré pattern.","keywords":["quasicrystal epitaxy","twisted moiré","striped pattern","colloidal monolayer","one-phonon theory","misfit angle","orientational epitaxy","decagonal optical lattice"],"falsifier":"An experimental or numerical check with a free (not ring-constrained) colloid monolayer on a decagonal optical lattice at the paper's length ratio: if the equilibrium angle is zero, or if the displacement field shows a sixfold moiré rather than parallel stripes, the central prediction fails.","tokens_in":7523,"feed_emoji":"🌀","tokens_out":3884,"duration_ms":33082,"temperature":0.7,"pith_summary":"The paper extends the long-known twisted-epitaxy result of Novaco and McTague from crystalline to quasicrystalline substrates, predicting that a crystalline monolayer on a quasicrystal will generally adopt a nonzero misfit angle. The central claim is that the resulting equilibrium distortion is not a hexagonal moiré but a striped pattern set by a single pair of reciprocal wavevectors. This follows analytically from a variational one-phonon calculation and matches numerical relaxation of a colloid monolayer on a decagonal optical lattice without fitting parameters. A sympathetic reader should care because the striped state is highly anisotropic, which should show up in direction-dependent friction and other interface properties.","feed_headline":"Monolayers on quasicrystals twist into striped moirés","feed_subtitle":"Analytic one-phonon theory and colloid simulations predict a striped, anisotropic moiré state instead of the usual hexagonal one.","key_machinery":"The argument uses a coherent-state variational phonon state for the harmonic monolayer, expanded to one-phonon order in the quasiperiodic substrate potential $V(\\mathbf{x})$. Momentum conservation forces each substrate Fourier component $\\mathbf{G}$ to combine with a monolayer reciprocal vector $\\boldsymbol{\\tau}$, leaving only the $\\mathbf{q}=\\mathbf{G}-\\boldsymbol{\\tau}$ points in the first Brillouin zone that satisfy Eq. (9). The resulting energy is Eq. (10), a sum over these $\\mathbf{q}$ points weighted by $|\\mathbf{G}\\cdot\\boldsymbol{\\epsilon}_{\\mathbf{q},s}|^2/\\omega_{\\mathbf{q},s}^2$, and the displacement field is Eq. (11). Because the decagonal potential has $p(p-1)=20$ $\\mathbf{G}$ vectors that appear only in inversion pairs, the shortest $\\mathbf{q}$ vector and its opposite dominate, producing stripes.","core_discovery":"The equilibrium orientation of a two-dimensional crystal on a quasiperiodic substrate is generically twisted: the optimal energy per particle is given by the one-phonon expression of Eq. (10) and the displacement field by Eq. (11). In the weak-coupling regime this analytic result reproduces numerical relaxation quantitatively, and the equilibrium pattern is dominated by a single pair of shortest wavevectors, producing parallel stripes rather than the hexagonal moiré of crystal-on-crystal epitaxy. For the studied length ratio $a_{\\rm pot}/a_{\\rm coll}=5.4\\,\\mu\\mathrm{m}/5.8\\,\\mu\\mathrm{m}$, the optimal twist angle is $\\theta_{\\rm opt}\\simeq 5.31^\\circ$.","pith_inferences":["If the rigid-rotation constraint is relaxed, real layers could break into domains with different twist directions, so the clean striped state may coexist with domain walls; this is an extension the paper does not address.","Varying the $a_{\\rm pot}/a_{\\rm coll}$ ratio should tune both the optimal twist angle and the stripe spacing, a testable prediction that follows directly from Eq. (10).","At stronger corrugation the one-phonon approximation degrades, and the displacement pattern is expected to acquire contributions from additional phonon modes, as the paper's own Fig. 4 suggests.","The same $\\mathbf{q}=\\mathbf{G}-\\boldsymbol{\\tau}$ logic applied to other quasiperiodic substrates with different rotational symmetry should produce stripe patterns of characteristic wavelength and orientation, which could be checked in optical-lattice colloid experiments."],"forward_implications":["Crystal-on-quasicrystal interfaces should generically be misaligned, so quasicrystalline substrates can impose twist on an adsorbed crystal without any external rotation.","The equilibrium moiré is striped, so any property controlled by the moiré, such as adhesion, electronic modulation, or chemical reactivity, will be strongly direction dependent.","In the weak-coupling regime the interface is superlubric, but kinetic friction should preferentially excite the shortest-$\"\"\"$q$ phonon and therefore be anisotropic.","The one-phonon formula provides a parameter-free benchmark for simulations and future experiments at small corrugation, and the optimal angle for the studied ratio is $\"\"\"$\\theta_{\\rm opt}\\simeq 5.31^\\circ$.","The striped pattern should persist for other combinations of incompatible symmetry, such as hexagonal-on-square or square-on-decagonal interfaces."],"supporting_citations":[{"why":"Establishes the original Novaco-McTague twisted-epitaxy theory for crystal-on-crystal that this paper extends to quasicrystals.","marker":"[4]"},{"why":"Provides the detailed weak-coupling one-phonon derivation on which Eq. (10) is modeled.","marker":"[5]"},{"why":"Reports colloidal monolayer experiments on decagonal optical lattices that motivate the model potential and the proposed test system.","marker":"[12]"},{"why":"Shows quasicrystalline and Archimedean-tiling conformations of colloids on quasicrystal substrates, the baseline against which the twisted state is new.","marker":"[14]"},{"why":"Demonstrates twisted epitaxy in colloidal monolayers on crystals, the closest prior system for the numerical setup.","marker":"[10]"},{"why":"LAMMPS is the simulation engine used for the numerical relaxation compared against Eq. (10).","marker":"[27]"}],"fun_headline_variants":["Quasicrystal monolayers twist into striped states","Twisted stripes on quasicrystal epitaxy","Monolayers on quasicrystals form striped twists","Unexpected striped twist in quasicrystal epitaxy","Quasicrystal layers twist into anisotropic stripes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the overlayer stays a uniform harmonic crystal that rotates as a rigid whole, with a global displacement field and no domain walls, dislocations, or spatially varying rotation.","fun_headline_variants_meta":{"raw":{"variants":["Quasicrystal monolayers twist into striped states","Twisted stripes on quasicrystal epitaxy","Monolayers on quasicrystals form striped twists","Unexpected striped twist in quasicrystal epitaxy","Quasicrystal layers twist into anisotropic stripes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000135,"raw_usage":{"total_tokens":1060,"prompt_tokens":778,"completion_tokens":282,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":394,"completion_tokens_details":{"reasoning_tokens":210}},"tokens_in":394,"tokens_out":282,"duration_ms":3103,"temperature":1.0,"reasoning_tokens":210,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:49:52.249384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An experimental or numerical check with a free (not ring-constrained) colloid monolayer on a decagonal optical lattice at the paper's length ratio: if the equilibrium angle is zero, or if the displacement field shows a sixfold moiré rather than parallel stripes, the central prediction fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the original Novaco-McTague twisted-epitaxy theory for crystal-on-crystal that this paper extends to quasicrystals."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the detailed weak-coupling one-phonon derivation on which Eq. (10) is modeled."},{"cited_title":"Mikhael, M","cited_arxiv_id":null,"evidence_quote":"Reports colloidal monolayer experiments on decagonal optical lattices that motivate the model potential and the proposed test system."},{"cited_title":"Mikhael, J","cited_arxiv_id":null,"evidence_quote":"Shows quasicrystalline and Archimedean-tiling conformations of colloids on quasicrystal substrates, the baseline against which the twisted state is new."},{"cited_title":"Mandelli, A","cited_arxiv_id":null,"evidence_quote":"Demonstrates twisted epitaxy in colloidal monolayers on crystals, the closest prior system for the numerical setup."}],"review_version":1}