{"id":"6d40c67a-26e4-47b7-b159-f766ddad3575","arxiv_id":"2501.11802","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Peropyrene's molecular plasmon is barely changed by nitrogen substitution, but oxygen substitution disrupts the pi-electron network, red-shifting and weakening the resonance, while charge doping tunes it.","lead":"Using computer simulations, this paper shows that adding nitrogen or oxygen atoms to a peropyrene molecule changes how its electrons oscillate, shifting and weakening its plasmonic resonances. A generalist might read it because it offers concrete design rules for building tiny optical antennas and color-tunable electrochromic devices from single molecules.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"S1's classification as a collective plasmon rests on PBE-based PI/GPI/TCM descriptors never benchmarked for peropyrene; if PBE inflation of Coulomb coupling overstates collectivity, the conjugation-controlled-plasmon narrative loses its object.","rationale":"The paper's central claim is that the longitudinal dipole plasmon of peropyrene is controlled by pi-conjugation, with nitrogen substitution barely perturbing it and oxygen substitution disrupting it, while charge doping tunes the resonance. The whole chain depends on identifying the dominant low-energy state (S1 in peropyrene, S3 in peropyrene-O) as a genuine collective plasmon. This identification rests on PI, GPI, and TCM, which are imported from prior work and applied without independent validation for peropyrene. No experimental absorption spectrum of peropyrene is compared, and no higher-level wavefunction method is used as a benchmark. The only cross-check is LR-TDDFT versus RT-TDDFT, but both use the PBE functional (Gaussian 16 and GPAW), so their agreement demonstrates internal consistency rather than physical correctness. PBE is known to overdelocalize pi-systems and may inflate the deexcitation coefficients in Eq. 2, increasing N and PI, and may shift TCM features in a way that mimics collective behavior. The manuscript's own limitation statement (Section 3.3, on electrochromic colors) acknowledges method dependence but does not apply the caveat to the plasmonic classification. The qualitative distinction between N and O doping could still survive a more accurate treatment, so the appropriate outcome is to preserve the reader's CONDITIONAL verdict, with the condition being an independent validation of the plasmonic descriptors. The absence of deposited data makes the required benchmark harder to reproduce, adding to the need for an explicit test.","tokens_in":14007,"tokens_out":5981,"duration_ms":69175,"concrete_test":"Recompute the S1 (peropyrene) and S3 (peropyrene-O) excitations with a range-separated hybrid (CAM-B3LYP) and with EOM-CCSD/cc-pVDZ on the same optimized geometry, then recalculate the number of excited electrons N (Eq. 2), PI, and the TCM halo displacement (red-halo versus dashed-line energy difference) at each level. If N drops below about 1.1, or if the primary TCM halo aligns with the excitation energy under any method, the collective assignment is a PBE artifact and the oxygen-doping explanation would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that substitution and doping tune a genuine dipolar molecular plasmon by altering pi-conjugation—stands or falls on the assignment of S1 (peropyrene) and S3 (peropyrene-O) as collective rather than single-particle excitations. The evidence for collectivity comes entirely from PBE-based LR-TDDFT with imported PI/GPI/TCM criteria (Eqs. 1–6), never benchmarked against experiment or a higher-level method for these molecules. The manuscript itself concedes in Section 3.3 that 'the absorption spectra are influenced by the computational methods employed,' but it applies that caveat only to predicted electrochromic colors, not to the plasmonic classification. PBE's known self-interaction error overdelocalizes pi-systems and can inflate the deexcitation coefficients (Y in Eq. 2), thereby raising the number of excited electrons N and PI, and it can shift the TCM halo relative to the excitation energy in a way that mimics strong electron-electron interaction. If S1 is actually a bright single-particle exciton, the downstream narrative—oxygen disrupting conjugation, reducing participating valence electrons, red-shifting and weakening the resonance—has no verified collective object. The LR-TDDFT versus RT-TDDFT agreement in Figure 1 is reassuring but uses the same PBE functional in both codes, so it tests internal consistency, not the validity of the plasmon assignment.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the optical excitations of peropyrene and two derivatives (peropyrene-N and peropyrene-O) using LR-TDDFT and RT-TDDFT with the PBE functional. The authors identify the main low-energy absorption peak as a longitudinal dipole molecular plasmon by applying the plasmonicity index (PI), the generalized plasmonicity index (GPI, here reported as E_plas), and transition contribution maps (TCMs). They argue that nitrogen substitution preserves the conjugated pi system and causes a slight blueshift, while oxygen substitution disrupts conjugation, reduces the number of participating valence electrons, and redshifts and weakens the resonance. They further show that charge doping tunes the resonance energy and that the trends differ between peropyrene/peropyrene-N and peropyrene-O, correlating with ELF-pi delocalization changes. The paper concludes that conjugated structure controls molecular plasmonic resonance and field enhancement.","tokens_in":14252,"tokens_out":3204,"duration_ms":35776,"significance":"If the plasmonic assignment is correct, the work provides a useful design rule for tuning molecular plasmons through heteroatom substitution and charge doping, and it demonstrates a multi-descriptor (PI, GPI, TCM, ELF-pi) workflow for small PAHs. The study has strengths: it combines two independent TDDFT implementations (Gaussian16 and GPAW), includes both LR and RT propagation, and reports qualitative agreement between them; it explicitly identifies the need for multi-dimensional analysis rather than relying on a single transition density. The main significance hinges on the reliability of the PI/GPI/TCM criteria to classify S1 (and analogous states) as collective plasmons rather than strongly correlated single-particle excitations; this point is not independently benchmarked for peropyrene, so the central claim carries a correctness risk that should be addressed before publication.","major_comments":[{"comment":"The classification of S1 (peropyrene), S3 (peropyrene-O), and related states as longitudinal dipole plasmons rests entirely on PBE-based PI, GPI, and TCM descriptors that are not benchmarked for peropyrene or its derivatives. PBE's known self-interaction error can over-delocalize pi electrons and may inflate the deexcitation coefficients Y in Eq. (2), increasing N and PI, and can also affect the TCM halo position. Since the central claim that conjugation controls a genuine collective plasmon depends on S1 being collective rather than a bright single-particle exciton, the manuscript should validate the assignment with a higher-level method (e.g., hybrid functional, BSE/GW, or coupled-cluster) or against experimental absorption spectra. The manuscript's own caveat in Section 3.3 that \"the absorption spectra are influenced by the computational methods employed\" is applied only to electrochromic colors, not to the plasmonic classification, and therefore does not address this concern.","section":"Section 3.1, Figure 2 and Eqs. (1)-(6)"},{"comment":"The replacement of GPI by E_plas is justified by assuming a constant damping energy Gamma for all excitations. This is a free parameter, and the paper does not test whether the qualitative conclusions (e.g., which states have 'high' or 'low' GPI) are robust to the value of Gamma. Since the paper explicitly states it reports E_plas \"referring to them as GPI values for brevity,\" the absolute criterion for calling a state plasmonic (GPI > 1 in the original formulation) is effectively bypassed. The authors should either specify the threshold they use for E_plas or show that the assignment of S1, S3, etc., as plasmonic is insensitive to the choice of Gamma within a reasonable range.","section":"Section 2.1, Eqs. (4)-(5)"},{"comment":"The explanation for the reduced absorption peak intensity in peropyrene-O attributes the decrease to a reduction in the number of valence electrons participating in the collective oscillation. However, the paper does not directly quantify this number: Eq. (2) defines N as the number of excited electrons, and PI is a monotonic function of N, but the manuscript never reports N or PI values for the longitudinal modes of the three neutral molecules. The ELF-pi maps in Figures 3c-3e show localization differences, but they are not converted into a quantitative measure of participating electrons. Without this quantitative link, the narrative that oxygen doping 'reduces the valence electrons participating in the plasmonic excitations' remains qualitative and should be supported by reporting N or a related integrated quantity for each molecule.","section":"Section 3.2, Figures 3a-3e"}],"minor_comments":[{"comment":"There is a typo in 'tuexcitations' where the word should be split as 'the excitations.'","section":"Introduction, paragraph 2"},{"comment":"The phrase 'the electron oscillation modes practicing in the plasmonic excitation' should read 'participating in the plasmonic excitation'.","section":"Section 3.2, paragraph after Figure 4"},{"comment":"The caption uses 'inserts' while the text uses 'insets'; please unify the terminology and ensure that the red and purple circle markings described in the caption are clearly visible in the figure.","section":"Section 3.1, Figure 1 caption"},{"comment":"The notation 'N𝑠𝑖𝑛𝑔𝑙𝑒' mixes a Unicode character with Latin text; please typeset it as N_single or a similar standard form.","section":"Section 2.1, Eq. (3)"},{"comment":"Reference [55] is incomplete in the reference list; the full citation for Multiwfn (Lu, T.; Chen, F. J. Comput. Chem. 2012, 33, 580-592) should be provided.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the computational workflow is clearly described. The main risk is the unvalidated transfer of PI/GPI/TCM criteria from prior work to peropyrene; a benchmark against a higher-level method or experiment would substantially strengthen the central claim. If the authors can provide such validation or at least a sensitivity analysis of the descriptors, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The worthwhile part of this paper is the object, not the method: peropyrene and its N/O derivatives have not been through the standard PI/GPI/TCM descriptor pipeline before, and the qualitative doping trends look believable. It is a competent application study, not a new principle.\n\nThe story — N substitution preserves pi-conjugation and slightly blue-shifts the longitudinal dipole mode; O substitution disrupts conjugation, cuts peak intensity, and red-shifts it; charge doping tunes the resonance, with ELF-pi bifurcation values tracking the delocalization changes — is coherent, and the evidence is internally consistent. LR- and RT-TDDFT (both PBE, different codes) give similar spectra, and PI, GPI, TCM, transition densities, and ELF-pi all point the same way. Credit where due: the authors flag in Section 3.3 that spectra and computed colors are method-dependent, and they treat mixed single-particle/plasmonic states explicitly rather than declaring every state a clean plasmon.\n\nThe real soft spot is the one the stress-test names: the collectivity label on S1 rests entirely on PBE-based descriptors imported without benchmarking against experiment or a higher-level method, and PBE's known overdelocalization could plausibly inflate PI. But the stress-test overshoots if it implies the paper falls if S1 is a bright exciton. The actual payload is the relative trends across molecules and doping levels at one fixed level of theory; those survive the absolute classification question. The LR/RT agreement tests numerical convergence, not the physical assignment, and the authors should say so.\n\nMinor: data are \"available upon request\" rather than deposited, and the Gaussian damping width is a free parameter — but since GPI is reported as Eplas with Gamma folded in, that is largely cosmetic.\n\nThis is for molecular plasmonics and electrochromic device modelers. It deserves a serious referee. My recommendation: conditional acceptance, with the main request being a one-paragraph statement of what the descriptors license — relative trends as the claim, not the absolute plasmon label — and ideally a single higher-level benchmark (BSE/GW or a range-separated functional) on the key excited states. I would not desk-reject it.","headline":"A competent application study where peropyrene is run through standard PI/GPI/TCM descriptors for the first time; the relative doping trends hold up, but the absolute plasmon label on S1 needs sharper benchmarking claims.","tokens_in":14788,"tokens_out":4197,"would_cite":true,"duration_ms":40980,"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":"The paper claims that the longitudinal dipole plasmonic resonance of peropyrene is controlled by its conjugated π-electron structure: nitrogen substitution keeps the conjugation intact and blueshifts the resonance, oxygen substitution…","keywords":["molecular plasmon","peropyrene","conjugated structure","time-dependent density functional theory","plasmonicity index","transition contribution map","charge doping","polycyclic aromatic hydrocarbons"],"falsifier":"A calculation of the three molecules' absorption spectra with a range-separated hybrid functional, or a measurement of those spectra in solution, would settle the claim: the predicted trend is a slight blueshift for nitrogen substitution, a clear redshift and intensity drop for oxygen substitution, and a monotonic redshift with charge doping for peropyrene and peropyrene-N; if the oxygen-doped peak is not red-shifted and weaker, the conjugation-control story fails.","tokens_in":13772,"feed_emoji":"🧪","tokens_out":4708,"duration_ms":43549,"temperature":0.7,"pith_summary":"This paper tries to establish that in peropyrene, a small polycyclic aromatic hydrocarbon, the longitudinal dipole plasmon resonance is controlled by the molecule's conjugated π-electron structure rather than by the atoms' elemental identity alone. It argues that nitrogen substitution at the ends preserves the conjugated framework and only slightly blueshifts the resonance, while co-doping with oxygen breaks π delocalization, lowering the number of electrons that oscillate and thereby red-shifting the resonance and cutting its intensity. It also argues that adding or removing charges modifies the conjugated structure and moves the resonance in ways that depend on which derivative is doped. A sympathetic reader would care because this gives a chemical design rule for molecular plasmonic devices: choose substituents that preserve or break conjugation to set the frequency, strength, and spatial field enhancement of the plasmon.","feed_headline":"Nitrogen blueshifts, oxygen red-shifts peropyrene's molecular plasmon","feed_subtitle":"Which atoms replace the carbon ends decides whether the resonance shifts blue or red and how much light the molecule concentrates.","key_machinery":"The central object is the longitudinal dipole plasmonic resonance of peropyrene, assigned by three descriptors: the plasmonicity index (PI), a count of how much the excitation involves deexcitation relative to a single-electron transition; the generalized plasmonicity index (GPI), the Coulomb self-energy of the transition density divided by a damping width; and transition contribution maps (TCMs), which show which occupied-to-unoccupied orbital transitions build an excited state and whether their energy differences match the excitation energy. The paper also uses the electron localization function restricted to π orbitals (ELF-π) to visualize whether conjugation is preserved or disrupted. These tools together identify which excited states are collective oscillations and connect the degree of π delocalization to the resonance peak's energy, intensity, and field distribution.","core_discovery":"On the paper's own terms, the central discovery is that the valence electrons taking part in the longitudinal dipole plasmonic oscillation of peropyrene are predominantly π-conjugated electrons, so anything that alters the π-conjugated structure directly alters the plasmon. Nitrogen substitution leaves the conjugated structure essentially intact, so peropyrene-N keeps a similar absorption peak with only a slight blueshift caused by a small geometric contraction. Oxygen substitution disrupts the conjugated structure, so peropyrene-O has a significantly weaker absorption peak and a red-shifted resonance energy because fewer valence electrons participate in the collective oscillation. Charge doping acts through the same channel: it changes the degree of π-electron delocalization, moving the resonance in different directions for peropyrene and peropyrene-N versus peropyrene-O. The paper further claims that these changes in conjugated structure change the electron oscillation mode and therefore where the local field enhancement is strongest (molecular center for peropyrene and peropyrene-N, molecular edge for peropyrene-O).","pith_inferences":["If the conjugation-preservation rule is general, then other polycyclic aromatic frameworks, such as larger peropyrene homologues or acenes, should show the same directionality: substituents that preserve delocalization shift the resonance only mildly, while conjugation-breaking substituents produce red-shifts with intensity loss.","The paper's field-enhancement maps suggest a testable design target: asymmetric doping could move the hot spot from the molecular center to the edges, which would matter for placing a reaction site or a second molecule at the maximum field.","A direct experimental check would be to measure the absorption spectra of peropyrene and its N/O-doped derivatives in an electrochromic gel; the predicted ordering (peropyrene-N slightly blue of peropyrene, peropyrene-O red and weak) should be visible if the computational assignment holds.","Because the PI, GPI, and TCM criteria were imported from larger systems, the same molecules could be re-examined with an independent plasmonicity diagnostic to see whether the $S_1$ assignment survives a change in method."],"forward_implications":["Substituent engineering that preserves π conjugation (like nitrogen at the ends) can shift a molecular plasmon's resonance without drastically weakening it, offering a mild blue-tuning knob.","Substituents that break π delocalization (like oxygen) provide a strong red-tuning knob at the cost of peak intensity and field-enhancement strength.","Charge doping of peropyrene and peropyrene-N gives a gradual red-shift of the longitudinal dipole resonance with increasing doping level, a broad tuning range useful for electrochromic devices.","The same resonance can be dark (zero oscillator strength) yet plasmonic, as with the $S_9$ mode, so identifying molecular plasmons requires multidimensional analysis rather than looking only at transition density or absorption strength."],"supporting_citations":[{"why":"Supplies the notion of plasmons in molecules and the electron-electron interaction criterion used to identify collective excitations.","marker":"[18]"},{"why":"Supplies the generalized plasmonicity index (GPI) method used to quantify the plasmonic character of each resonance.","marker":"[24]"},{"why":"Supplies the transition contribution map (TCM) tool used to decompose excited states into orbital transitions.","marker":"[25]"},{"why":"Supplies the plasmonicity index (PI) descriptor used to distinguish collective from single-particle excitation.","marker":"[26]"},{"why":"Provides the end and central plasmon resonance behavior in linear atomic chains that underpins the electron-number and intensity argument.","marker":"[35]"},{"why":"Supplies experimental evidence of molecular plasmons in polycyclic aromatic hydrocarbons that motivates the tunability study.","marker":"[38]"},{"why":"Provides the experimental electrochromic-device context and charge-tunable color behavior that the doping analysis builds toward.","marker":"[43]"},{"why":"Provides the earlier acene result with electron-hole doping trends that the paper contrasts when finding coexisting plasmonic and single-particle excitations.","marker":"[49]"},{"why":"Supplies the electron localization function (ELF) method that the paper adapts as ELF-π to visualize π delocalization and conjugation disruption.","marker":"[65]"}],"fun_headline_variants":["N and O steer peropyrene's molecular plasmon to blue and red","Peropyrene's plasmon turns blue with N, red with O","Oxygen adds red, nitrogen adds blue to peropyrene's plasmon","N substitution preserves peropyrene's blue plasmon; O breaks it red","Conjugation controls peropyrene's plasmon: N blue, O red"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central story assumes that the plasmonicity index, generalized plasmonicity index, and transition contribution maps correctly tell collective plasmon oscillations apart from strongly correlated single-particle transitions in small molecules like peropyrene; if those descriptors mislabel a single-particle transition as a plasmon (or vice versa), the claim that the conjugated structure controls a collective oscillation would not follow.","fun_headline_variants_meta":{"raw":{"variants":["N and O steer peropyrene's molecular plasmon to blue and red","Peropyrene's plasmon turns blue with N, red with O","Oxygen adds red, nitrogen adds blue to peropyrene's plasmon","N substitution preserves peropyrene's blue plasmon; O breaks it red","Conjugation controls peropyrene's plasmon: N blue, O red"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00142,"raw_usage":{"total_tokens":5727,"prompt_tokens":935,"completion_tokens":4792,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":4697}},"tokens_in":551,"tokens_out":4792,"duration_ms":36674,"temperature":1.0,"reasoning_tokens":4697,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:50:25.447044+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation of the three molecules' absorption spectra with a range-separated hybrid functional, or a measurement of those spectra in solution, would settle the claim: the predicted trend is a slight blueshift for nitrogen substitution, a clear redshift and intensity drop for oxygen substitution, and a monotonic redshift with charge doping for peropyrene and peropyrene-N; if the oxygen-doped peak is not red-shifted and weaker, the conjugation-control story fails.","supporting_citations":[{"cited_title":"and Jacob, C","cited_arxiv_id":null,"evidence_quote":"Supplies the notion of plasmons in molecules and the electron-electron interaction criterion used to identify collective excitations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the generalized plasmonicity index (GPI) method used to quantify the plasmonic character of each resonance."},{"cited_title":"P., Kuisma, M., Puska, M","cited_arxiv_id":null,"evidence_quote":"Supplies the transition contribution map (TCM) tool used to decompose excited states into orbital transitions."},{"cited_title":"and Mandado, M","cited_arxiv_id":null,"evidence_quote":"Supplies the plasmonicity index (PI) descriptor used to distinguish collective from single-particle excitation."},{"cited_title":"and Gao, S","cited_arxiv_id":null,"evidence_quote":"Provides the end and central plasmon resonance behavior in linear atomic chains that underpins the electron-number and intensity argument."},{"cited_title":"E., Manjavacas , A","cited_arxiv_id":null,"evidence_quote":"Supplies experimental evidence of molecular plasmons in polycyclic aromatic hydrocarbons that motivates the tunability study."},{"cited_title":"J., Lauchner, A., Cui, Y","cited_arxiv_id":null,"evidence_quote":"Provides the experimental electrochromic-device context and charge-tunable color behavior that the doping analysis builds toward."},{"cited_title":"B., Wang, H","cited_arxiv_id":null,"evidence_quote":"Provides the earlier acene result with electron-hole doping trends that the paper contrasts when finding coexisting plasmonic and single-particle excitations."}],"review_version":1}