{"id":"3a86f7c2-c202-43c3-a0d0-7daeb559f15b","arxiv_id":"2506.06241","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In situ GIFAD and SDRS during PbI2 MBE growth on graphene reveal a 1% tensile strain in the first monolayer that relaxes by 3-5 ML, and a 50 meV ARPES shift attributed to charge transfer from graphene to PbI2.","lead":"Researchers grew thin films of lead iodide on graphene, monitoring the growth in real time with two synchronized surface probes. They report that the first film layer is stretched by about 1 percent and that graphene's electronic bands shift about 50 meV, which they attribute to electron transfer from graphene to lead iodide.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1% monolayer strain claim is statistically marginal against the stated GIFAD precision (±0.01 Å^-1) and Fig. 5a has no error bars, so the structural pillar of the 'far exceeds vdW' conclusion is not yet established.","rationale":"The reader's weakest_assumption is exactly the most load-bearing concern. The manuscript's conclusion 'far exceeds that expected from vdW coupling' is primarily a statement about interface interaction strength, and the only quantitative structural evidence is the ~1% tensile strain. The text provides reciprocal vector precision of ±0.01 Å^-1 but shows no error bars in the strain-versus-thickness plot, leaving the statistical significance unresolved. The charge-transfer interpretation is secondary, and the paper itself hedges by offering two possible explanations ('much stronger than pure vdW' vs. 'retains vdW character but involves a charge transfer'). Thus, if the strain is not real, the structural pillar collapses, and the optical peak shift loses its proposed cause. A concrete re-analysis of the raw data can settle this. Since the concern is addressable without overturning the well-supported layer-by-layer growth and azimuthal alignment, the conditional verdict remains appropriate; no change to the reader's recommendation is needed.","tokens_in":13993,"tokens_out":9868,"duration_ms":95025,"concrete_test":"Re-analyze the raw GIFAD images used to construct Fig. 5a: fit the first-order PbI2 diffraction peaks at each thickness (especially 1 ML and >5 ML) with the same peak-fitting routine, propagate the fit uncertainty to the lattice parameter, and perform a two-sample t-test comparing the 1 ML value to the >5 ML value. If the difference is not significant at p<0.05 (or if the 1 ML lattice parameter does not exceed the relaxed value by more than the combined uncertainty), then the 1% strain claim is not supported and the central conclusion must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the PbI2–graphene interaction 'far exceeds that expected from vdW coupling' rests most heavily on the ~1% tensile strain measured in the first PbI2 monolayer (Fig. 5a). The paper quotes reciprocal vectors with an uncertainty of ±0.01 Å^-1 (about 0.6% relative). A 1% lattice expansion corresponds to an approximately 1% decrease in reciprocal vector, i.e., about 0.016 Å^-1 at g=1.59 Å^-1, which is only ~1.6 times the stated uncertainty on a single peak position. If the same precision applies to the strained 1 ML data, the strain is not significant at the 2σ level. Moreover, Fig. 5a is plotted without error bars, and the thickness at exactly 1 ML is assigned from the first GIFAD oscillation maximum, which may be affected by island coalescence, peak-fitting ambiguities, or the coexistence of exposed graphene and island-free regions. Because the optical correlation in Fig. 6c (the blue shift of the 2.93 eV interband peak) is interpreted as a consequence of this strain, the entire structural–optical narrative hinges on a measurement whose precision is not demonstrated. If the strain is an artifact, the structural and optical legs of the 'strong interaction' claim lose their quantitative support, leaving only the indirect and hedged charge-transfer inference.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports an in-situ correlated study of PbI2 growth on bilayer graphene/SiC(0001) using grazing-incidence fast atom diffraction (GIFAD) and surface differential reflectance spectroscopy (SDRS) during molecular beam epitaxy, supplemented by LEED, core-level photoemission, and ARPES. The authors find layer-by-layer growth, a unique azimuthal alignment (PbI2 armchair parallel to graphene zigzag), and report a ~1% tensile strain in the first PbI2 monolayer that relaxes by 3–5 ML, a correlated blue shift of the 2.93 eV interband transition for thin layers, and a ~50 meV shift of the graphene bands toward the Fermi level for one monolayer, interpreted as charge transfer from graphene to PbI2. The central conclusion is that the PbI2–graphene interface interaction \"far exceeds that expected from vdW coupling.\"","tokens_in":14296,"tokens_out":4848,"duration_ms":48352,"significance":"The paper's methodological contribution is strong: the simultaneous, real-time combination of GIFAD and SDRS during van der Waals epitaxy is unusual and enables a direct structural–optical correlation on the same growth sequence. The GIFAD and LEED data convincingly establish layer-by-layer growth and a well-defined azimuthal orientation, and the ARPES measurement is carefully discussed with respect to the mixed-layer graphene substrate. If the 1% strain claim is robust, the work would provide quantitative evidence that PbI2–graphene coupling is not merely weak vdW bonding, with implications for interface engineering in halide-based optoelectronics. However, the quantitative central claim rests heavily on a strain measurement whose statistical precision is not demonstrated in the paper as presented.","major_comments":[{"comment":"The claim of a ~1% tensile strain in the first PbI2 monolayer is not supported with error bars in Fig. 5a, and the stated precision of the GIFAD reciprocal-vector measurement, ±0.01 Å⁻¹ (about 0.6%), is comparable to the reported 1% expansion. For a reciprocal vector of 1.59 Å⁻¹, a 1% strain corresponds to a change of about 0.016 Å⁻¹, only about 1.6 times the quoted single-peak uncertainty. The authors must provide the uncertainty on each lattice-parameter data point, the number of independent measurements or growth runs, and a statistical test (e.g., a t-test or confidence interval) showing that the 1 ML value differs significantly from the relaxed value. Without this, the structural pillar of the \"far exceeds vdW\" conclusion is not established.","section":"Growth and structural properties, Fig. 5a"},{"comment":"The assignment of the first GIFAD oscillation maximum to exactly one complete monolayer is central to the strain measurement. In a layer-by-layer growth, the first reflectivity maximum may not correspond to a perfectly closed monolayer if island coalescence or partial second-layer nucleation occurs before the first layer is complete. The manuscript should justify the 1 ML calibration—for instance by comparing the GIFAD oscillation phase with ex-situ atomic force microscopy or with the evolution of the diffraction-peak intensities—and explain how the reciprocal-vector measurement at that coverage averages over any coexistence of bare graphene and PbI2 islands, which could bias the apparent lattice parameter.","section":"Growth and structural properties, Figs. 3b and 5a"},{"comment":"The blue shift of peak 2 (from 2.92 eV at 1 ML to 2.95 eV at 4 ML) is about 30 meV, and the attribution of this nonmonotonic shift to the lattice strain observed by GIFAD is speculative without quantitative support. The authors should report the SDRS spectral resolution, the peak-fitting uncertainty, and error bars on the peak positions, and ideally demonstrate a point-by-point correlation between the strain relaxation and the optical shift on the same sample. As written, the link between Fig. 5a and Fig. 6c is qualitative and does not independently corroborate the strain claim.","section":"Optical response, Fig. 6c"},{"comment":"The conclusion that the PbI2–graphene interaction \"far exceeds that expected from vdW coupling\" is stronger than the evidence presented. In the Results section, the authors themselves hedge: \"the interface interaction is either much stronger than pure vdW and involves chemical effects or, most probably, retains its vdW character but involves a charge transfer allowed by a favourable band alignment.\" To support the strong claim, the paper needs a quantitative benchmark—for example, comparison with known vdW epitaxy systems or with calculated binding energies and charge-transfer magnitudes—or the conclusion should be softened to state that the interaction involves measurable strain and charge transfer beyond the simplest physisorption picture.","section":"Summary and conclusions"}],"minor_comments":[{"comment":"The figure caption labels two panels as \"b\" (the surface reflectivity and the lattice-mismatch inset); the panels should be renumbered and referred to consistently in the text.","section":"Figure 3"},{"comment":"The sentence \"The observed strain can be explained by the lattice mismatch ... with (2g_PbI2 − g_Gr)/2g_PbI2 = 7.2%\" is confusing because a 7.2% mismatch does not by itself explain a 1% strain. The authors should clarify the proposed mechanism, e.g., partial strain accommodation by the flexible PbI2 layer, and state why only a fraction of the mismatch is accommodated.","section":"Growth and structural properties"},{"comment":"The text states that for PbI2 the armchair direction is the most corrugated and \"the other direction does show diffraction peaks from the hexagonal PbI2 lattice,\" but the assignments of armchair and zigzag directions in Fig. 4a and 4b should be made explicit and consistent with the earlier definitions for graphene.","section":"Results and discussion"},{"comment":"Reference [1] is incomplete; it lists only authors, volume, and page numbers (306, 666) without the article title or a complete journal citation.","section":"References"},{"comment":"The extraction of the optical band gap using Tauc plots is described only briefly; the authors should specify whether the direct or indirect Tauc formula was used for each thickness and how the absorption coefficient was derived from the differential reflectance data.","section":"Optical response, Fig. 6d"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of a surface-science or 2D-materials journal and the experimental approach is novel. The main risk is that the central quantitative claim (1% strain) may not survive scrutiny once error bars and statistics are added; if the strain cannot be established, the conclusion should be correspondingly weakened. I would encourage the editor to request the raw data for Fig. 5a and Fig. 6c as part of the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nWhat you should know: this is a genuinely new experimental study—real-time correlated GIFAD and SDRS during MBE growth of PbI2 on graphene—and it produces the first quantitative thickness-resolved picture of the interface: layer-by-layer growth, a unique azimuthal alignment, a 1% tensile strain in the first ML that relaxes by 3-5 ML, a non-monotonic shift of the 2.93 eV interband transition, and an ARPES shift of about 50 meV. The technique combination is the real novelty; few groups can do this.\n\nThe paper does several things well. The GIFAD oscillations cleanly show layer-by-layer growth. The diffraction patterns establish a well-defined orientation (PbI2 armchair parallel to graphene zigzag), consistent with Sinha. The optical data as a function of thickness show that peak 2 behaves differently from the others, which is a real observation worth understanding. The ARPES shows a systematic shift of all graphene features toward the Fermi level when 1 ML PbI2 is present—that's a real effect.\n\nThe soft spots are mostly about how the claims are framed. The headline 1% strain is the biggest one. The paper quotes reciprocal-vector precision of ±0.01 Å-1, which is ~0.6% relative. A 1% lattice change corresponds to about 1.6 times that uncertainty, and the strained 1-ML point in Fig. 5a has no error bars. As it stands, that measurement is not shown to be significant at the 2σ level. This is a fixable problem—report the spread over multiple runs or the uncertainty on the peak fitting—but until then the structural pillar of the \"far exceeds vdW\" conclusion is shaky. The optical blue shift of peak 2 is consistent with strain, but it's not a unique fingerprint; the authors should be clearer that it's a correlation, not proof. The ARPES shift is the most direct evidence for charge transfer, but it's still an inference; the text says so, the abstract doesn't. I'd also push back on the abstract's \"demonstrated by an energy shift\"—that's stronger than what the data can prove.\n\nWho is this for? Experimentalists working on vdW epitaxy, 2D halides, and in-situ growth monitoring. The GIFAD/SDRS combination will be of interest to the MBE and surface-science crowds. It deserves a serious referee: the dataset is real, the method is novel, and the weaknesses are addressable with added analysis and more careful wording. I'd send it to review, with the expectation that the strain claim needs to be backed up or softened.\n\nBest,\n[Your name]","headline":"A real methodological advance in correlating structure and optics during vdW epitaxy, but the headline 1% strain claim is statistically under-supported as presented.","tokens_in":14903,"tokens_out":3788,"would_cite":false,"duration_ms":37503,"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":"A single monolayer of PbI2 on graphene is stretched by about 1% and pulls charge out of graphene, evidence that this van der Waals interface is far stronger than expected.","keywords":["PbI2","graphene","van der Waals epitaxy","interface strain","charge transfer","GIFAD","surface differential reflectance spectroscopy","ARPES"],"falsifier":"Measure the in-plane lattice parameter of a single PbI2 monolayer on graphene by an independent real-space or diffraction method, such as scanning tunneling microscopy of the atomic lattice or low-energy electron diffraction spot positions with error bars, and check whether it is really about 4.61 Å rather than 4.56 Å; repeated grazing-incidence diffraction measurements at several azimuths and beam energies with propagated uncertainties would settle whether the 1% strain is real.","tokens_in":13808,"feed_emoji":"⚛️","tokens_out":8472,"duration_ms":69074,"temperature":0.7,"pith_summary":"This paper sets out to show that the interface between one monolayer of lead iodide (PbI2) and graphene is not the weak, passive contact usually assumed for van der Waals stacks. During layer-by-layer growth on graphene/SiC(0001), the first PbI2 monolayer is stretched by about 1% in plane, the strain relaxes only after 3–5 monolayers, the 2.93 eV interband absorption tracks that strain, and graphene's bands shift about 50 meV toward the Fermi level, which the authors read as charge transfer from graphene to PbI2. The evidence comes from synchronizing two real-time probes, helium-atom diffraction and surface reflectance spectroscopy, on the same growing film, plus photoemission on the finished monolayer. If the paper is right, heterostructure devices made from PbI2 and similar halide layers cannot treat the first interface as electronically inert.","feed_headline":"Graphene stretches the first PbI2 layer by 1 percent","feed_subtitle":"Real-time diffraction and optics show the strain fades by 5 monolayers and graphene bands shift 50 meV.","key_machinery":"The argument is carried by three synchronized or complementary probes. Grazing incidence fast atom diffraction (GIFAD) sends low-energy helium atoms along the surface; because helium is inert, the diffraction pattern reflects the electron-density corrugation of the last atomic plane, giving both reciprocal-lattice vectors (hence in-plane lattice parameter and strain) and relative diffraction-order intensities (hence changes in charge-density distribution). Surface differential reflectance spectroscopy (SDRS) records the ultraviolet/visible absorption of the growing film in real time, linking the 2.93 eV transition to the strained monolayer. Angle-resolved photoemission spectroscopy measures the band dispersion of graphene before and after one monolayer, providing the 50 meV shift attributed to charge transfer. The synchronization of GIFAD and SDRS during growth is what lets the authors correlate structural strain, optical peak position, and thickness on the same sample.","core_discovery":"On the authors' account, PbI2 grows on bilayer graphene/SiC(0001) in a layer-by-layer van der Waals-like mode with the armchair direction of PbI2 parallel to the zigzag direction of graphene. The first monolayer is under roughly 1% tensile strain, deduced from diffraction reciprocal vectors, and relaxes to the bulk lattice parameter over 3–5 monolayers, while the surface electron-density distribution keeps evolving until 9–10 monolayers. The 2.93 eV interband transition shows a non-monotonic shift that tracks this strain, first moving to higher energy and returning by 10 monolayers. Photoemission shows all graphene spectral features shift about 50 meV toward the Fermi level when one PbI2 monolayer is present, which the authors interpret as charge transfer from graphene to PbI2. Together these observations support their stated conclusion that the PbI2–graphene interaction far exceeds that expected from van der Waals coupling.","pith_inferences":["If the interaction is truly stronger than van der Waals, one testable extension is to grow the same PbI2 monolayer on an inert layered substrate such as hBN; a much smaller strain and no 50 meV shift would identify a graphene-specific chemical contribution rather than a generic overlayer effect.","The diffraction-intensity evolution up to 9–10 monolayers suggests the electronic perturbation outlives the lattice strain, which could be checked by thickness-dependent work-function or core-level measurements to see whether charge transfer continues beyond the first monolayer.","The correlation between the 2.93 eV peak shift and strain implies that SDRS could be used as a general in-situ strain monitor for other 2D halide heterostructures, provided a reference interband or excitonic transition is identified.","The photoemission substrate contains an admixture of one to three graphene layers, so part of the 50 meV shift could reflect charge redistribution among graphene layers rather than a pure graphene-to-PbI2 transfer; repeating the measurement on monolayer graphene would separate these contributions."],"forward_implications":["Ultrathin PbI2 devices must treat the first monolayer as strained, since its band structure, optical gap, and likely exciton properties differ from bulk-like PbI2 up to 3–5 monolayers.","The 2.93 eV interband transition can serve as an in-situ optical fingerprint of interface strain, because its energy shift tracks the diffraction-measured lattice relaxation.","The measured 50 meV graphene band shift implies that a PbI2 overlayer electronically dopes the graphene, so PbI2/graphene contacts are not electrically neutral interfaces.","Layer-by-layer growth with a fixed azimuthal alignment over large areas makes molecular beam epitaxy a viable route to thickness-controlled PbI2 films for photodetectors and perovskite-related devices."],"supporting_citations":[{"why":"Supplies single-layer PbI2 calculations showing structural flexibility and strain-induced band modulation, used to explain why a 1% tensile strain is plausible.","marker":"[12]"},{"why":"Provides thickness-dependent band-gap calculations that assign the indirect-to-direct transition and the optical gap trend with layer number.","marker":"[13]"},{"why":"Reports atomic structure and epitaxial alignment of monolayer PbI2 nanodisks on graphene, supporting the observed armchair-on-zigzag orientation.","marker":"[17]"},{"why":"Explains GIFAD's sensitivity to the electron density of the last atomic plane, the basis for interpreting the diffraction-order intensities as an electronic fingerprint.","marker":"[24]"},{"why":"Describes surface differential reflectance spectroscopy, the optical probe synchronized with GIFAD to correlate structure and optics during growth.","marker":"[25]"},{"why":"Defines van der Waals epitaxy, the weak-coupling growth framework whose expected interface strength the paper argues is exceeded.","marker":"[33]"},{"why":"Gives the bulk PbI2 lattice parameter used as the unstrained reference value for measuring the 1% expansion.","marker":"[34]"},{"why":"Supports the band-alignment and interfacial charge-transfer mechanism invoked to explain the strong electronic interaction.","marker":"[36]"},{"why":"Provides strain- and layer-dependent optical calculations for PbI2, used to attribute the 2.93 eV peak shift to lattice strain rather than thickness.","marker":"[40]"},{"why":"Reports charge transfer at a MoSe2/graphene interface, the comparison case for interpreting the photoemission shift as a proximity effect across a van der Waals gap.","marker":"[56]"}],"fun_headline_variants":["PbI2 on graphene shows 1% strain in the first layer, then relaxes","Graphene stretches the first PbI2 layer by 1%, then relaxes by 5 layers","Charge transfer from graphene to PbI2 shifts bands by 50 meV","Layer-by-layer PbI2 growth on graphene shows strain and charge transfer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on measuring a 1% stretch in the first monolayer with a diffraction technique whose stated precision is about 0.6%, and the one-monolayer data points have no error bars; if that apparent stretch is a fitting or thickness-calibration artifact, the main claim collapses.","fun_headline_variants_meta":{"raw":{"variants":["PbI2 on graphene shows 1% strain in the first layer, then relaxes","Graphene stretches the first PbI2 layer by 1%, then relaxes by 5 layers","Charge transfer from graphene to PbI2 shifts bands by 50 meV","Layer-by-layer PbI2 growth on graphene shows strain and charge transfer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00156,"raw_usage":{"total_tokens":6267,"prompt_tokens":1014,"completion_tokens":5253,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":5162}},"tokens_in":630,"tokens_out":5253,"duration_ms":33943,"temperature":1.0,"reasoning_tokens":5162,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:58:08.635960+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the in-plane lattice parameter of a single PbI2 monolayer on graphene by an independent real-space or diffraction method, such as scanning tunneling microscopy of the atomic lattice or low-energy electron diffraction spot positions with error bars, and check whether it is really about 4.61 Å rather than 4.56 Å; repeated grazing-incidence diffraction measurements at several azimuths and beam energies with propagated uncertainties would settle whether the 1% strain is real.","supporting_citations":[],"review_version":1}