{"id":"1178309e-fa13-4d09-a5f8-5de696c6f936","arxiv_id":"2501.18952","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"RT-TDDFT simulations of L/D-Phenylglycinol on achiral and Na-doped chiral aluminum clusters show strong coupling, polaritonic mode formation, and plasmon-enhanced circular dichroism that grows with the number of coupled molecules.","lead":"The authors used atom-scale quantum simulations to show that chiral molecules strongly couple to aluminum nanoparticles, forming hybrid light-matter states and amplifying the molecules' circular dichroism signals. The work suggests design rules for highly sensitive chiral sensors and optoelectronic devices, though some headline claims about 'coupling dichroism' outrun the data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported CD enhancement factors rest on an unsupported additivity assumption: subtracting the isolated cluster's rotatory strength from the coupled system does not isolate molecular CD, and large-gap/denominator checks suggest it fails.","rationale":"The paper is a serious RT-TDDFT study; the strong-coupling signatures (LP/UP splitting, TCM sign flips, induced-density phase reversal, gap dependence) are internally consistent and independently evidenced. My concern is not about the existence of strong coupling, nor about the raw CD spectra, but about the specific quantity the paper calls the 'enhancement chirality factor.' The factor is the ratio of a subtracted quantity to an isolated-molecule peak. This is only meaningful if the cluster contribution is perfectly additive. The paper itself flags the definition as 'nebulous,' and the strong-coupling analysis shows delocalized transitions, so the subtraction cannot be assumed valid. The smoking gun is already in the submitted data: at 8 Å, where coupling is described as nearly gone, the factor is still ~1.8–3.5 instead of ~1; and for N=8 the denominator is near zero, producing 50× factors. A simple large-gap convergence check would settle whether the subtraction is valid. I therefore agree with the reader's conditional verdict: the qualitative conclusions may stand, but the quantitative enhancement factors need redefinition or error analysis before acceptance.","tokens_in":22504,"tokens_out":6125,"duration_ms":56774,"concrete_test":"Recompute ΔR(d)=R_total(d)−R_cluster_isolated and η=ΔR/R_mol_isolated for Al201@L-PG at d=8, 15, 25, and 50 Å using the same RT-TDDFT settings. If additivity holds, η(d)→1 as d→∞; if η(50 Å) remains above ~1.2, the subtraction does not isolate molecular CD and the enhancement factors in Tables S7–S10 are invalid as a measure of molecular CD. As a second check, re-evaluate the Al201@(L-PG)8 factor after shifting the selected peak by ±0.05 eV; if the factor changes by more than 20%, the reported 50.2× enhancement is an artifact of near-zero denominator.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—molecular CD enhanced 3–7× at 2 Å gap—is computed from Eq. (1.2): R_enh = R_total − R_cluster, divided by R_molecule. This assumes the cluster's rotatory strength is unchanged by the molecule and strictly additive. In a strongly coupled hybrid, the Kohn-Sham transitions are delocalized and the induced densities (Fig. 3) show cluster and molecule contributions reorganized; there is no justification for additivity. The paper itself calls the enhancement definition 'nebulous' and performs manual peak selection. Internal evidence already suggests the subtraction fails: at d=8 Å the authors state interaction is weakened/negligible, yet Tables S7–S10 give enhancement factors of 1.77–3.54, not ~1. If the subtraction correctly isolated molecular CD, factors should converge to 1 as coupling vanishes; no 15–50 Å factors are reported even though the absorption spectra fully overlap the isolated cluster there. The N=8 factors are additionally unstable: Table S11 gives Al201@(L-PG)8 factor 50.20 because the isolated-molecule denominator is −33.90 (near zero), so slight peak shifts or broadening changes the factor dramatically. Thus the 'enhanced molecular CD' numbers are not a robust observable; they mix cluster, interaction, and selection artifacts.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports real-time time-dependent density functional theory (RT-TDDFT) simulations of chiral phenylglycinol (L/D-PG) molecules coupled to Al201 and Na-doped chiral Al197Na4 clusters. The authors identify strong coupling with Rabi splitting into lower and upper polaritons, supported by transition contribution maps (TCMs) and induced-density analysis. They extract coupling strengths and linewidths by fitting a velocity-coupled harmonic oscillator model to the simulated absorption spectra. The central quantitative claim is that the molecular circular dichroism (CD) is enhanced by factors of roughly 3–7 at a 2 Å gap, with the enhancement depending on gap, cluster chirality, and the number of coupled molecules. Multi-molecule systems show collective effects, particularly an anomalously large enhancement factor at N=8.","tokens_in":22735,"tokens_out":4458,"duration_ms":41255,"significance":"If the quantitative claims were established, this would be a valuable first-principles account of plasmon-enhanced chirality and coupling dichroism at the atomic scale, extending previous TDDFT studies of CD enhancement to the strong-coupling regime. The work has notable strengths: it uses a validated, open-source RT-TDDFT implementation; it provides a systematic gap and molecule-number series; and the TCM/induced-density analysis gives a physically appealing picture of bonding and antibonding polaritonic modes. However, the central quantitative result—the CD enhancement factor—depends on an unsupported additivity assumption and is internally inconsistent with the large-gap limit. The paper also overstates chirality-dependent coupling and decay-rate differences that are within fitting uncertainty. These issues affect the main claims rather than only the presentation.","major_comments":[{"comment":"The enhancement factor in Eq. (1.2) is defined by subtracting the isolated cluster rotatory strength from the coupled-system rotatory strength, implicitly assuming that the cluster contribution is unchanged by the molecule and strictly additive. In a strongly coupled hybrid the Kohn-Sham transitions are delocalized, and the induced densities in Figure 3 show that cluster and molecule contributions are reorganized; the paper itself calls the enhancement definition 'nebulous' (page 10). This assumption fails an internal consistency test: at d=8 Å, where the text states the interaction is 'weakened or even negligible' (page 9), Tables S7–S10 report enhancement factors of 1.77–3.54, which should approach 1 if the subtraction isolated the molecular CD. No enhancement factors are reported for 15–50 Å, where the absorption spectra fully overlap the isolated cluster (Figures 2a and 2c). The quantitative claim of 3–7× molecular CD enhancement is therefore not established.","section":"Gap Dependent Enhanced Chirality, Eq. (1.2)"},{"comment":"The large N=8 enhancement factor (50.20 for Al201@(L-PG)8) is an artifact of a near-zero denominator: the isolated (L-PG)8 rotatory strength is -33.90, so small spectral shifts or broadening change the quotient dramatically. The conclusion of a distinct enhancement mechanism for N=8 (pages 19–20) rests on this unstable ratio. A robust observable should be insensitive to minor peak-selection variations; the present quantity is not.","section":"Table S11 and Figure 9d"},{"comment":"The abstract and conclusions claim that both the coupling factor and the decay rate are modulated by the chirality of the molecules and the cluster. However, the fitted g values for L-PG versus D-PG differ by at most 0.005 eV (e.g., 0.352 vs. 0.348 eV at 2 Å for Al201@L/D-PG, Tables S1 and S2), and the linewidth differences are ≤0.01 eV. The text itself acknowledges that coupling the same cluster to a chiral molecule or its enantiomer 'does not significantly affect the coupling strength' (page 10). Without uncertainty estimates from the MCMC fit, these differences are within likely fitting error, so the chirality-dichroism claim for g and decay rates is overstated.","section":"Tables S1–S4 and Figure 2f"}],"minor_comments":[{"comment":"The sentence 'the conventional criterion for the strong coupling state (i.e., )54' has a missing formula; please insert the explicit inequality (e.g., 2g > (γ_ex + γ_pl)/2).","section":"Page 9, strong-coupling criterion"},{"comment":"Since the emcee MCMC sampler is used, the authors should report credible intervals or standard errors for the fitted parameters (g, ω, γ) rather than only point values.","section":"Tables S1–S4, fitting"},{"comment":"The legend includes '9Å' although the main text sets the gap values to 2–8, 15, 25, and 50 Å; please clarify whether a 9 Å calculation was performed or whether this is a typographical error.","section":"Figure S3 legend"},{"comment":"The phrase 'achiral/chiral clusters induce significant spectral shifts' is imprecise: the cluster absorption peak remains near 7.7 eV, and the shifts appear in the polariton peaks. Please rephrase to distinguish cluster resonance from hybrid-mode energies.","section":"Abstract"},{"comment":"The phrase 'a time limiter was introduced' is vague; please specify what this limiter controls (e.g., maximum propagation time or a stopping criterion).","section":"Computational Details"}],"recommendation":"major_revision","confidential_remarks":"The RT-TDDFT simulations themselves appear carefully done and the qualitative polaritonic-mode analysis is interesting. The main problem is that the quantitative CD enhancement factors are built on an additivity assumption that the data themselves contradict at large gaps. I would ask the authors to either validate Eq. (1.2) by reporting enhancement factors at 15–50 Å and showing convergence to 1, or to reframe the paper as a qualitative study of polaritonic CD without claiming specific enhancement factors. The N=8 result should be reanalyzed with a stable definition of enhancement. The authors should also temper the coupling-dichroism claims for g and linewidths unless fitting uncertainties are reported. This is fixable within revision, hence major_revision rather than reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about this one. First, it's a genuinely new set of RT-TDDFT calculations: strong coupling of a chiral molecule (L/D-phenylglycinol) to an aluminum cluster, including a Na-doped chiral cluster, with systematic gap dependence and a multi-molecule N=2,4,8 study. The raw spectra, TCMs, and induced densities give a credible picture of bonding/antibonding polaritons and collective enhancement. Second, the headline claims about 'coupling dichroism' and the CD enhancement factors are not supported by the paper's own numbers.\n\nThe qualitative physics is solid. The methods are standard LCAO-RT-TDDFT, and the TCM analysis is a good way to show which transitions make up the polaritons. The gap-dependent coupling strengths and the N-dependent lower-polariton intensity are plausible and interesting. This is a real extension of earlier work by Härkönen et al. to chiral clusters and multi-molecule effects.\n\nThe soft spots are in the quantitative claims. The fitted coupling strengths for L vs D enantiomers differ by at most 0.005 eV (Tables S1–S4), well within typical fitting error for spectra like these; the authors themselves say molecular chirality does not significantly affect the coupling strength. So the title and abstract overstate the dichroism in the coupling. More importantly, the enhancement chirality factor R_enh = R_total − R_cluster is computed assuming the cluster's rotatory strength is additive and unchanged in the coupled system. Strong coupling delocalizes the transitions; the induced densities show the cluster and molecule reorganize. The paper's own data show the subtraction fails: at d=8 Å, where coupling is negligible, the enhancement factors are still 1.77–3.54 instead of converging to 1. And the N=8 factors (50.20 for Al201@(L-PG)8) come from a near-zero denominator in the isolated molecule. The authors even call the definition 'nebulous.' That is honest, but it means the numbers in Tables S7–S11 are not a robust observable.\n\nThis paper is worth engaging with—the simulations are first-principles and the qualitative story is new. But it needs major revision: report fitting uncertainties, define a better enhancement measure that doesn't rely on additivity, and adjust the claims accordingly. I would send it to peer review, but I would expect the referee to push hard on the enhancement factor.","headline":"Solid first-principles RT-TDDFT study of chiral molecule–Al cluster strong coupling, but the 'coupling dichroism' and CD enhancement factors are overstated and rest on a subtraction the paper's own data contradict.","tokens_in":23329,"tokens_out":3165,"would_cite":true,"duration_ms":29999,"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":"Atomistic simulations show that a single chiral molecule strongly couples to a plasmonic aluminum cluster at gaps up to 8 Å, forming bonding and antibonding polaritonic modes that amplify the molecule's circular dichroism signal by…","keywords":["Plasmonics","RT-TDDFT","Circular Dichroism","Strong Coupling","Induced Density","TCM","Electric-field Enhancement"],"falsifier":"Recompute the rotatory strength of the coupled system using a transition-density decomposition that isolates the molecular component without subtracting the cluster spectrum; if the inferred molecular CD peaks no longer show a 3- to 7-fold enhancement relative to the isolated molecule at a 2 Å gap, the reported enhancement factors are artifacts.","tokens_in":22203,"feed_emoji":"🌀","tokens_out":10480,"duration_ms":86909,"temperature":0.7,"pith_summary":"Strong coupling between a molecule and a plasmonic nanoparticle is usually studied for achiral systems; this paper asks whether the same coupling can amplify the weak circular dichroism (CD) of a chiral molecule and whether the nanoparticle's own chirality plays a role. The authors simulate, with real-time time-dependent density functional theory, single L- and D-Phenylglycinol molecules at controlled distances from an achiral aluminum cluster (Al201) and a chirally doped cluster (Al197Na4). They find that at gaps of 8 Å and below the systems enter the strong-coupling regime, producing lower (bonding) and upper (antibonding) polaritonic modes. The molecular CD signal is enhanced by factors of roughly 3 to 7 at a 2 Å gap, with the chiral cluster giving larger enhancement; the coupling strength and decay rate are modulated by the handedness of both components. Adding more molecules (up to eight) increases the lower polariton absorption and, at eight, reveals an additional enhancement channel from molecule-molecule interactions.","feed_headline":"Plasmonic coupling boosts molecular CD up to 7-fold","feed_subtitle":"At 2 Å gaps, a chiral molecule and an aluminum cluster form hybrid polaritons that amplify circular dichroism.","key_machinery":"The main mechanism is the formation of hybrid polaritonic modes in a strongly coupled molecule–cluster system. The authors use real-time time-dependent density functional theory (RT-TDDFT) to simulate the full electronic dynamics, and then analyze the results with transition contribution maps (TCMs), which decompose the optical absorption into individual electron–hole transitions, and with induced-density plots to visualize the phase relationship between the molecular and plasmonic dipoles. Coupling strengths and linewidths are extracted by fitting the absorption spectra to a velocity-coupled harmonic oscillator model. The circular-dichroism enhancement factor is defined as the ratio of the molecular rotatory strength in the coupled system (after subtracting the isolated cluster's contribution) to that of the isolated molecule at the same peak.","core_discovery":"The central discovery, stated in the paper's own terms, is that strong plasmon–molecule coupling, rather than merely a static field enhancement, is responsible for the observed CD amplification. The authors show that the coupled system's transition contribution maps and induced densities display two polaritonic modes: a lower-energy mode where the molecular and plasmonic dipoles are in phase, and an upper-energy mode where they are out of phase. This antibonding/bonding structure is the signature of strong coupling. In that regime, the molecular CD peaks are enhanced by factors between roughly 3 and 7 at a 2 Å gap, and the enhancement is larger when the cluster itself is chiral (Al197Na4) than when it is achiral (Al201). The paper further claims that the coupling strength g, the linewidths, and even the sign of the induced response depend on the handedness of the molecule and of the cluster, and that the number of coupled molecules tunes the lower polariton intensity in a collective manner.","pith_inferences":["The subtraction procedure used to isolate the molecular CD contribution assumes the cluster's rotatory strength is unchanged by coupling; if that additivity breaks down, the reported enhancement factors (3–7) may be overestimates. I would test this by recomputing the cluster's rotatory strength in the presence of the molecule's static charge distribution.","The paper finds little difference in coupling strength between L- and D-Phenylglycinol; this may be specific to the molecular orientation and the symmetric position of the molecule on the {100} facet. Rotating the molecule or using a different chiral molecule could reveal handedness-dependent coupling, which would be a useful test.","The collective enhancement at eight molecules is attributed to intermolecular interactions, but the paper does not separate the contribution of molecule–molecule electronic coupling from near-field or geometric effects. This could be disentangled by comparing the same eight-molecule system with the molecules placed far apart on a flat surface.","The simulations use a specific exchange-correlation functional and neglect thermal and solvent effects; real solution-phase measurements might show different enhancement magnitudes, so experimental verification is needed."],"forward_implications":["If strong coupling at gaps up to 8 Å is robust, single-molecule CD spectroscopy on plasmonic nanoparticles becomes feasible without the need for large ensembles, opening a route to label-free chiral sensing.","Doping a metal cluster to make it chiral (here with four Na atoms) provides an additional control knob for CD enhancement, suggesting that chiral plasmonic materials can be engineered to boost molecular signals.","The dependence of coupling strength and linewidth on molecular and cluster handedness implies that CD spectra of strongly coupled systems contain information about the absolute configuration of the molecule, not just its concentration.","The roughly linear growth of the lower polariton with the number of molecules up to four, and the deviation at eight, indicates that molecular coverage is a tunable parameter for optimizing collective strong-coupling effects."],"supporting_citations":[{"why":"Supplies the RT-TDDFT method and the Al201 cluster template for strong plasmon–molecule coupling.","marker":"[15]"},{"why":"Previous TDDFT study of molecular CD enhancement on silver nanoparticles, which this paper extends to strong coupling and chiral clusters.","marker":"[23]"},{"why":"Provides the LCAO-RT-TDDFT implementation for computing rotatory strengths (circular dichroism) in metal–organic clusters.","marker":"[35]"},{"why":"Introduces the velocity-coupled oscillator model used to fit the spectra and extract coupling strengths and linewidths.","marker":"[18]"},{"why":"Establishes the coupled-oscillator description of macroscopic observables in strongly coupled nanoparticle–molecule systems.","marker":"[55]"},{"why":"Defines the strong-coupling criterion used to justify the interpretation of the spectra.","marker":"[54]"}],"fun_headline_variants":["Chiral coupling dichroism yields 7-fold CD gain","Strong coupling between chiral molecule and plasmon boosts CD 7x","Hybrid polaritons amplify molecular chirality 7-fold","Coupling dichroism in chiral molecule-plasmon system amplifies CD","Chiral molecule-plasmon coupling yields 7x CD enhancement"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reported molecular CD enhancement assumes the cluster's rotatory strength is the same in the coupled system as in isolation, so that subtracting it from the total rotatory strength yields a purely molecular signal.","fun_headline_variants_meta":{"raw":{"variants":["Chiral coupling dichroism yields 7-fold CD gain","Strong coupling between chiral molecule and plasmon boosts CD 7x","Hybrid polaritons amplify molecular chirality 7-fold","Coupling dichroism in chiral molecule-plasmon system amplifies CD","Chiral molecule-plasmon coupling yields 7x CD enhancement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001642,"raw_usage":{"total_tokens":6565,"prompt_tokens":1025,"completion_tokens":5540,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":5450}},"tokens_in":641,"tokens_out":5540,"duration_ms":34229,"temperature":1.0,"reasoning_tokens":5450,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T21:49:51.405737+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the rotatory strength of the coupled system using a transition-density decomposition that isolates the molecular component without subtracting the cluster spectrum; if the inferred molecular CD peaks no longer show a 3- to 7-fold enhancement relative to the isolated molecule at a 2 Å gap, the reported enhancement factors are artifacts.","supporting_citations":[],"review_version":1}