{"id":"0d2b4843-fe2c-4435-98f1-df4d46274482","arxiv_id":"2412.14078","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"On Cu(111), dynamical screening from the substrate reduces or quenches the local spin moments of adsorbed transition-metal phthalocyanines in an orbital-filling-dependent way.","lead":"This paper models transition-metal phthalocyanine molecules on a copper surface using density functional theory plus a quantum impurity model. It finds that the metal contact dynamically screens the molecular spin moments, reducing or fully quenching them for most of the transition-metal series.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; the main caveat is the fixed U/J parameters, which the paper itself acknowledges as a limitation.","rationale":"The reader's weakest assumption correctly identifies the fixed U/J parameters as a potential soft spot. My stress-test analysis reached the same concern: the quantitative values of S_scr and the exact categorization of partial vs. full quenching depend on U and J, which are fixed across all TMs with no sensitivity study. This is a legitimate caveat for a CONDITIONAL rating. However, the paper's central qualitative trend—Cr retains the largest screened moment, Fe is strongly suppressed, Co/Ni are quenched—is supported by the orbital-resolved physics (d_xy nonbonding, d_R filling, charge/spin fluctuations), by the DFT+U comparisons, and by multiple experimental references for FePc, CoPc, and CuPc. The mechanism is internally consistent and not circular, since U and J are not fitted to the experimental moments. Therefore, while a sensitivity analysis on U/J would strengthen the paper, its absence does not invalidate the core claim. My recommendation is to keep the original CONDITIONAL verdict (i.e., 'UNCHANGED' relative to the reader's verdict). The stress-test did not identify a more severe or more fundamental objection, such as an internal inconsistency or a missing proof, beyond the acknowledged parameter uncertainty. I also note that the paper explicitly discusses the temperature dependence and time-scale dependence, which adds credibility. The suitable concrete test is a small set of additional CTQMC runs varying U and J, which directly addresses the reader's caveat and could upgrade the paper to ACCEPT if the trend is stable.","tokens_in":17888,"tokens_out":878,"duration_ms":9131,"concrete_test":"Repeat the CTQMC calculations for CrPc, FePc, and CoPc with U=3.5 eV and U=4.5 eV (keeping J=1.0 eV fixed), and also with J=0.8 eV and J=1.2 eV for fixed U=4.0 eV, re-evaluating S_scr and the d_R occupations. If the ordering (Cr > Mn > Fe > Co ≈ Ni in S_scr) and the qualitative classification (retained vs. suppressed vs. quenched) remain unchanged, the central claim is parameter-robust. Report the sensitivity of the FePc/CoPc boundary.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that dynamical screening, driven by orbital-selective hybridization and electron correlation, suppresses or quenches local spin moments in TMPc/Cu(111) interfaces, with a clear trend across the 3d series (Cr retains, Fe suppressed, Co/Ni quenched). The key assumption is that the three free parameters (U=4.0 eV, J=1.0 eV, and the DFT-derived hybridization) adequately capture the physics for all seven TMs. The paper explicitly states: 'In all of our calculations, we have used U = 4.0 eV and J = 1.0 eV.' No sensitivity analysis is performed. If the true U or J changes across the series (e.g., later 3d elements typically have larger U), the quantitative screening values—and especially the FePc vs. CoPc boundary—could shift. However, this is a standard parameter choice in the field, and the qualitative trend is robust to modest variations in U and J because the mechanism (d_xy being nonbonding and d_R filling) is driven mainly by hybridization and occupation. The paper is not circular: it uses U and J as fixed inputs, not fitted to experimental moments. The DFT+U comparison and experimental references provide independent support. Thus, the verification of the fixed U/J is a worthwhile sensitivity check, but it does not undermine the central qualitative conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates dynamical screening of local spin moments in transition-metal phthalocyanines (TiPc through NiPc, plus CuPc) adsorbed on Cu(111). Using DFT to construct multiorbital Anderson impurity models, solved with continuous-time quantum Monte Carlo, the authors compute the spin susceptibility chi(tau) and extract instantaneous (tau=0) and screened (tau=beta/2) effective moments. They find that orbital-dependent hybridization and 3d filling control the degree of screening: CrPc retains a large screened moment, FePc is strongly suppressed, (Co-Ni)Pc are quenched, and CuPc keeps a sizable moment. The mechanism is attributed to the strongly hybridized, screened d_{x^2-y^2} orbital, the weakly hybridized nonbonding d_xy orbital, and the metallic out-of-plane orbitals, with charge fluctuations away from d_R half-filling enhancing screening. The paper also compares the many-body results with DFT+U and discusses implications for probe-dependent experimental measurements.","tokens_in":18141,"tokens_out":7667,"duration_ms":67134,"significance":"If the results hold, the paper provides a systematic and physically transparent explanation of how metal-molecule contact reduces or quenches local spin moments, with a clear orbital-selective mechanism and a design rule based on d_R filling and hybridization. The systematic scan across the first-row transition-metal series, the orbital-resolved analysis of chi(tau), the use of a full rotationally invariant Coulomb tensor, and the explicit comparison with DFT+U and experiments are notable strengths. The central claim is not circular: U and J are fixed inputs, not fitted to the experimental moments being explained. The main caveats are the fixed Coulomb parameters across the series and the absence of reported Monte Carlo statistical errors, both of which affect the quantitative boundary between partial screening and full quenching.","major_comments":[{"comment":"The Coulomb parameters are fixed to U = 4.0 eV and J = 1.0 eV for all seven transition-metal species, with no sensitivity analysis. The quantitative screened moments in Figure 3 and, in particular, the boundary between strong suppression (FePc) and full quenching (CoPc, NiPc) depend on these choices, and U and J are known to vary across the 3d series. I request a sensitivity test for the borderline cases (e.g., U = 3.5-4.5 eV, J = 0.8-1.2 eV) or per-element parameters from constrained random-phase approximation, with a statement of whether the retention/suppression/quenching categorization of Figure 3 survives.","section":"Methods (many-body calculations)"},{"comment":"No Monte Carlo statistical errors are reported for the CTQMC results. Since the paper's central claim includes vanishing screened moments for CoPc and NiPc and a near-vanishing moment for FePc, the reader cannot judge whether these values are zero or merely small within statistical uncertainty. Please report error bars or an accuracy measure for chi(tau = beta/2) and the extracted S_scr values, at least for the systems with small screened moments.","section":"Results, Spin-Spin Correlation and Effective Local Moments"},{"comment":"The statement that 'quenching of spin moments of FePc on Cu surfaces has been seen in recent experiments' is supported by Refs. 48 and 53, but Ref. 48 reports an NO2-driven ferrous-to-ferric transition and Ref. 53 reports magnetic-anisotropy switching on oxidized Cu(110), neither of which directly demonstrates spin-moment quenching on Cu(111). Please either cite direct experimental evidence for FePc moment quenching on Cu surfaces or rephrase the claim to reflect the cited findings.","section":"Results, Quenching of Local Moments"}],"minor_comments":[{"comment":"In the sentence 'The screened moments show a reduction of and, in cases, even quenching of the local moments', the word 'of' appears to be missing its object; please rephrase, e.g., 'show a reduction and, in some cases, even quenching'.","section":"Results, Spin-Spin Correlation and Effective Local Moments"},{"comment":"The clause 'therefore, CuPc shows the largest screened moment among the TMPc molecules' appears to be a typo for CrPc, given the surrounding discussion and Figure 3; please correct and ensure consistency with the 'Complete Filling' paragraph.","section":"Results, At Half-Filling"},{"comment":"The displayed formulas for chi, M, and the charge/spin fluctuation correlators are garbled in the manuscript text; please provide clean typeset definitions, particularly Eq. (4) for M^2 and Eqs. (5)-(6) for the fluctuation correlators.","section":"Equations (2)-(6)"},{"comment":"The hybridization function symbol (Delta) is missing from the text/caption where 'imaginary part of the hybridization function' appears; please ensure the symbol renders correctly.","section":"Figure 2 caption"},{"comment":"Reference 39 and Reference 47 are the same paper (Arruda et al., Phys. Chem. Chem. Phys. 2020, 22, 12688-12696); please merge or renumber them.","section":"References"},{"comment":"Please correct 'Hellman Feynman forces' to 'Hellmann-Feynman forces'.","section":"Methods"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is the published ACS Nano article (2023, 17, 5974-5983) posted to arXiv. If this is being evaluated as a new journal submission, the editor should clarify the intended contribution relative to the published version; as a preprint it is a useful record of the published work, but the novelty for a new submission would need to be stated explicitly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nFor a paper on dynamical spin screening at metal-molecule interfaces, here is my read. It systematically scans Ti–Ni phthalocyanines on Cu(111) with DFT + Anderson impurity model + CTQMC, and the central finding is credible: instantaneous moments are atomic-like, but long-time screened moments vary strongly, with CrPc retaining the largest (S_scr ≈ 1.6), FePc strongly suppressed, and Co/Ni quenched. The mechanism is orbital-selective hybridization combined with the filling of the d_R subspace (d_xy plus out-of-plane orbitals); d_xy is weakly coupled and, when half-filled, preserves the moment.\n\nWhat is genuinely new: the full first-row scan under a common protocol, and the orbital-resolved susceptibility analysis that yields a simple filling rule. The work is careful in separating instantaneous (τ=0) from screened (τ=β/2) moments, and it supports the trend with computed occupations, charge/spin fluctuations, and several experimental comparisons (Fe quench, Co quench, CuPc persistence). The DFT+U comparison is also instructive: it explains where static mean-field fails.\n\nNow the soft spots. The fixed Hubbard U=4.0 eV and Hund's J=1.0 eV for all seven TMs is the main caveat. The paper states this in Methods. Real U varies by ~1 eV across the 3d series, and a sensitivity sweep would be needed to confirm the exact boundary between partial suppression and full quenching, especially for FePc. Also, no CTQMC statistical errors are reported, so we cannot gauge the noise on the S_scr values. These are real but minor; they do not undermine the qualitative ordering, which is driven by occupation and hybridization shape, not by fine-tuning.\n\nOne more point: the paper's connection to experimental probe time scales (fast probes like XAS vs. slow probes like INS) is plausible but speculative; it is an extrapolation, not a calculation. That is worth keeping in mind when citing.\n\nWho is this for? Researchers in molecular spintronics or strongly correlated interfaces. It is a solid reference and deserves peer review—the flaws are in the details, not the architecture. I would engage with it and recommend citing with a note on the fixed U/J.\n\nBest","headline":"Systematic Ti–Ni scan with fixed U/J gives a credible orbital-filling rule for screening at TMPc/Cu(111); the trend holds, but quantitative boundaries need a sensitivity check.","tokens_in":18681,"tokens_out":3531,"would_cite":true,"duration_ms":30992,"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":"At a copper contact, orbital-dependent dynamical screening reduces or quenches the local spin moments of most transition-metal phthalocyanines.","keywords":["molecular spintronics","transition-metal phthalocyanines","dynamical screening","local magnetic moment","Anderson impurity model","quantum fluctuations","Cu(111) surface","orbital-dependent hybridization"],"falsifier":"Recompute the screened moments with U scanned over a plausible range (for example, 3 to 5 eV) and J from 0.7 to 1.3 eV while holding all other inputs fixed. If CoPc or NiPc acquires a nonzero screened moment anywhere in that range, or if the ordering CrPc > FePc reverses, then the retention-versus-quenching categorization is an artifact of the fixed parameters rather than a robust prediction.","tokens_in":17694,"feed_emoji":"🧲","tokens_out":9623,"duration_ms":78384,"temperature":0.7,"pith_summary":"This paper tries to establish that the spin moment of a transition-metal ion in a phthalocyanine molecule is not a fixed property once the molecule touches a metal surface: the metallic contact dynamically screens the moment, so what a measurement sees depends on the probe's timescale. The paper builds a five-orbital Anderson impurity model from density-functional-theory data for Ti, V, Cr, Mn, Fe, Co, and Ni phthalocyanines on Cu(111) and compares the instantaneous moment (short times, near-atomic) with the long-time screened moment. It finds the screened moment is reduced in most cases and fully quenched for CoPc and NiPc, while CrPc keeps the largest screened moment because its four nonbonding $d_R$ orbitals sit near half-filling. The result matters because metal-molecule contact is unavoidable in spintronic devices, and the moment that persists after screening is the one a device could actually use.","feed_headline":"Copper contact quenches most molecular spin moments","feed_subtitle":"Dynamical screening shrinks or erases the local moments of most metal-phthalocyanine molecules on Cu(111).","key_machinery":"The central object is the imaginary-time spin susceptibility $\\chi(\\tau)$ of a five-orbital Anderson impurity model, a model of a localized correlated orbital coupled to a bath of itinerant electrons, solved with a continuous-time quantum Monte Carlo solver. $\\chi(0)$ gives the square of the instantaneous local moment, while $\\chi(\\tau = \\beta/2)$ gives the square of the dynamically screened moment after quantum fluctuations have acted, and comparing these two values quantifies the screening. The mechanism is driven by the orbital-resolved hybridization function: strong in-plane hybridization screens $d_{x^2-y^2}$, weak hybridization lets $d_{xy}$ retain spin when half-filled, and metallic out-of-plane hybridization screens $d_{xz}$, $d_{yz}$, and $d_{z^2}$ with an efficiency that grows with filling.","core_discovery":"The central claim is that orbital-dependent hybridization and Coulomb repulsion jointly produce charge and spin fluctuations that screen the local 3d moment of adsorbed transition-metal phthalocyanines, with the degree of screening set by how far the weakly bonded $d_R$ subspace (the $d_{xy}$, $d_{xz}$, $d_{yz}$, and $d_{z^2}$ orbitals) sits from half-filling. The strongly bonded $d_{x^2-y^2}$ orbital loses its moment almost entirely; the nearly nonbonding $d_{xy}$ orbital can preserve part of the moment when half-filled; and the out-of-plane orbitals hybridize with the copper surface and screen strongly, more so as their occupation grows. The resulting ordering is that CrPc retains the largest screened moment, FePc is strongly suppressed, CoPc and NiPc are quenched, and CuPc keeps a sizable moment carried by its $d_{x^2-y^2}$ hole. The paper states this as a general guide: for a large long-time moment, keep the $d_R$ subspace close to half-filling, and away from half-filling expect charge fluctuations to suppress or erase the moment.","pith_inferences":["The fixed Coulomb parameters U = 4.0 eV and J = 1.0 eV may shift the boundary between suppression and quenching if varied across the series; the qualitative orbital-filling rule is more robust than the exact screened-moment values.","The same half-filling rule might predict behavior on other metallic contacts: weaker hybridization (for example, on Au or Ag, or through a spacer layer) should push the quenching threshold toward heavier transition metals.","A direct test would compare spin moments from a sub-picosecond probe (XAS or XMCD) and a slow probe (neutron scattering or SQUID) on identical TMPc/Cu(111) samples; a systematic mismatch would confirm the timescale picture.","Ligand or substrate engineering that pins $d_{xy}$ near half-filling could preserve a screened moment even in late transition-metal phthalocyanines where it is currently quenched."],"forward_implications":["Fast probes such as X-ray absorption or emission spectroscopy should see near-atomic local moments on TMPc/Cu(111), while slower probes such as inelastic neutron scattering should see reduced or vanishing moments.","Static density-functional calculations with a Hubbard U will agree with the screened moment near half-filling (CrPc) but will overestimate the moment for MnPc, FePc, and TiPc, where charge fluctuations are strong.","The vanishing long-time moments computed for CoPc and NiPc provide a mechanism for experiments that report no persisting spin on these molecules on copper surfaces.","For device design, the practical route to a persistent molecular spin at a metal contact is to keep the weakly hybridized $d_R$ orbitals near half-filling, as in CrPc.","Lowering temperature increases screening, so the FePc moment should be further suppressed at low temperature, matching existing experiments."],"supporting_citations":[{"why":"Defines the screened local moment through the long-time limit of the spin susceptibility, the quantity the paper uses as the screened spin.","marker":"[42]"},{"why":"Shows dynamical screening of local moments in correlated metals, the bulk analogue of the molecule-surface effect studied here.","marker":"[43]"},{"why":"Supplies the timescale argument that fast probes see instantaneous moments while slow probes see screened moments.","marker":"[45]"},{"why":"Supports local-moment dynamics and screening in doped charge-transfer insulators, justifying the use of the long-time susceptibility value.","marker":"[44]"},{"why":"Establishes that the $d_R$ subspace of the 3d shell carries the magnetism of metal-organic molecules and motivates the many-body treatment.","marker":"[52]"},{"why":"Provides experimental evidence that FePc moments are suppressed on surfaces, a comparison point for the computed FePc quenching.","marker":"[48]"},{"why":"Reports adsorption-induced magnetic changes for FePc on copper, supporting the strong suppression found for Fe.","marker":"[53]"},{"why":"Shows CoPc on Cu(111) is highly susceptible to screening, matching the computed full quenching of the Co moment.","marker":"[55]"},{"why":"Reports a persisting local moment in CuPc layers, consistent with the paper's finding that CuPc retains a screened moment.","marker":"[50]"}],"fun_headline_variants":["Most spin moments vanish when molecules meet copper","Orbital filling decides which spin moments survive Cu contact","Half-filled d orbitals resist screening at Cu contacts","Screening at metal-molecule junction kills most spin moments"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative results assume one fixed Coulomb interaction strength (U = 4.0 eV) and one fixed Hund coupling (J = 1.0 eV) applied to all seven transition-metal ions; if these parameters should differ across the series, the point at which screening becomes full quenching would move.","fun_headline_variants_meta":{"raw":{"variants":["Most spin moments vanish when molecules meet copper","Orbital filling decides which spin moments survive Cu contact","Half-filled d orbitals resist screening at Cu contacts","Screening at metal-molecule junction kills most spin moments"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000995,"raw_usage":{"total_tokens":4227,"prompt_tokens":970,"completion_tokens":3257,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":3195}},"tokens_in":586,"tokens_out":3257,"duration_ms":22014,"temperature":1.0,"reasoning_tokens":3195,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:30:08.306151+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the screened moments with U scanned over a plausible range (for example, 3 to 5 eV) and J from 0.7 to 1.3 eV while holding all other inputs fixed. If CoPc or NiPc acquires a nonzero screened moment anywhere in that range, or if the ordering CrPc > FePc reverses, then the retention-versus-quenching categorization is an artifact of the fixed parameters rather than a robust prediction.","supporting_citations":[],"review_version":1}