{"id":"8cef0d40-9d16-4941-8286-31818f379102","arxiv_id":"2509.09998","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Solvent molecules inside DBTTF:HATCN cocrystals act as the main electron acceptors in the computed charge distribution, even though Raman data give a solvent-independent donor-acceptor transfer of about 0.11 e.","lead":"Six new crystal forms of the organic pair DBTTF and HATCN show that solvent molecules trapped between the molecules can take up electron density themselves, sometimes outpacing the intended electron acceptor. If confirmed, trapped solvent becomes a design knob for redirecting charge in optoelectronic cocrystals, not just an unwanted impurity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hirshfeld-only charge assignment leaves solvent-as-acceptor claim unvalidated; a Bader/Hirshfeld-I cross-check on D:A:THF2 would settle it.","rationale":"The reader's weakest assumption (Hirshfeld dependence) is exactly the load-bearing point. The paper's central novelty is that solvent participates as primary electron acceptor; if that is wrong, the abstract's headline claim is unsupported, though the structural and optical data remain valuable. The Raman calibration is also fragile (single reference, no error bars), but even a robust 0.11 e DCT would not establish solvent acceptance. Therefore the decisive test is a partitioning-independent or experimental verification of solvent charge. I do not see an internal inconsistency in the DFT calculations themselves; the issue is external validity of the charge assignment. The paper should be accepted conditionally pending this check, matching the reader's verdict.","tokens_in":104,"tokens_out":3139,"duration_ms":48658,"concrete_test":"Recompute the partial charges for D:A:THF2 from the same PBE0 density using two independent schemes: Bader analysis (e.g., via VASP/PAW or FHI-aims) and Hirshfeld-I (iterative). If the combined charge on the two THF molecules is not at least -0.20 e and/or HATCN does not become net positive, the 'primary electron acceptor' claim is an artifact of the chosen partitioning. A useful ancillary check is to compute C≡N frequency shifts for HATCN with explicit THF molecules at the experimental geometry to see whether the observed Raman shift is compatible with the Hirshfeld-derived charge distribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Table 2 is the sole quantitative support for the central claim that solvent molecules are primary electron acceptors. These values are Hirshfeld charges from PBE0 calculations on PBE-relaxed experimental geometries. Hirshfeld partitioning is a non-unique density decomposition: it compares the molecular density to a promolecule built from neutral atoms, so it does not measure a physical charge. In D:A:THF2 the two THF molecules intercalate with contacts as short as 2.54 Å to DBTTF and 2.77 Å to HATCN; overlapping densities are exactly where partitioning ambiguities are largest. No alternative population analysis (Bader, DDEC, Hirshfeld-I, electrostatic-potential fit) or experimental observable (core-level shift, vibrational mode of solvent, conductivity) is reported to confirm that THF/ACE are negatively charged. Moreover, the paper's own electronic-structure results are in tension with the claim: PDOS (Fig. 5) shows solvent states only in the deeper valence region and the CT1 excitation is assigned to DBTTF→HATCN; if THF accepted ~0.35 e, the frontier acceptor states would be expected to shift to the solvent. The Raman-derived 0.11 e is also only a net D→A estimate and cannot discriminate where the excess electron resides. Thus the headline mechanism rests on a single unverified partitioning scheme.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports six new DBTTF:HATCN cocrystals: five solvated 1:1 polymorphs grown from THF, acetone, or acetonitrile solution, and one solvent-free 3:2 polymorph grown by horizontal vapor deposition. The authors combine single-crystal X-ray diffraction, Raman spectroscopy, optical absorption, and DFT (PBE/PBE0, FHI-aims) to examine how included solvent molecules affect structure, charge transfer, and electronic/optical properties. They find that the optical gap and the main charge-transfer excitations are largely insensitive to solvent inclusion. From the Raman C≡N shift of HATCN, they extract a degree of charge transfer (DCT) of about 0.11 e via a two-part decomposition: a constant DCT-induced shift of 4.81 cm⁻¹ plus a dipole-induced shift linear in the solvent gas-phase dipole moment (Eq. 2). DFT Hirshfeld charges indicate that in solvated structures the solvent molecules carry substantial negative charge, up to -0.35 e total in D:A:THF2, and in D:A:THF2 and D:A:ACE2 HATCN itself becomes a net electron donor. This leads to the central claim that solvent molecules act as primary electron acceptors rather than passive templates, while the overall electronic structure remains robust.","tokens_in":18681,"tokens_out":4549,"duration_ms":56564,"significance":"If the solvent-as-acceptor mechanism is correct, this is a conceptually interesting result that would expand the design space for ternary organic cocrystals. The paper has clear strengths: six new crystal structures with deposited CCDC numbers, transparent experimental methods, a combined experimental/computational workflow, and a public data repository. The Raman and optical data are reported in a reproducible way. However, the headline claim rests on a single population-analysis scheme (Hirshfeld) and on a two-point linear decomposition whose constancy is partly built into the model. These issues are load-bearing because the abstract and conclusions generalize from the Hirshfeld charges to a 'primary electron acceptor' role for solvents. With additional charge-partitioning cross-checks or an experimental probe of the solvent charge state, the claim could become solid; as it stands, it is not yet established.","major_comments":[{"comment":"The central claim that solvents are primary electron acceptors rests entirely on Hirshfeld partial charges from PBE0 calculations on PBE-relaxed experimental geometries. Hirshfeld partitioning is promolecule- and basis-set-dependent, and the short solvent contacts (e.g., d_D:sol = 2.54 Å and d_A:sol = 2.77 Å in D:A:THF2, Table 1) place overlapping densities precisely in the region where the partitioning is most ambiguous. The -0.16/-0.19 e on THF and +0.09 e on HATCN in D:A:THF2 are therefore not robust without a second independent scheme. Please provide Bader, DDEC, Hirshfeld-I, or electrostatic-potential-fitted charges, or an experimental probe (solvent vibrational mode, core-level shift, or transport/conductivity signature). If an alternative partitioning changes the sign or magnitude of the solvent charge, the headline mechanism must be revised.","section":"Table 2 and 'Partial charge analysis from DFT'"},{"comment":"The solvent-as-acceptor picture is in tension with the computed electronic structure. The PDOS in Fig. 5 shows THF states only in the deeper valence region, and CT1 is assigned to DBTTF→HATCN transitions. The authors state that DFT 'points to an excess of charge accumulated on the solvent molecule rather than a state-selective charge-transfer,' but Hirshfeld charges do not identify the orbital character of the transferred charge. If THF were truly a primary acceptor with -0.35 e, one would expect acceptor-like states near the frontier; none is shown. This does not disprove the claim, but the paper should explicitly reconcile the PDOS/transition assignment with the large Hirshfeld solvent charges, or soften the 'primary electron acceptor' formulation to 'net electron accumulation according to Hirshfeld partitioning.'","section":"Optical Absorption and Electronic Structure (Fig. 5)"},{"comment":"Equation (2) decomposes the observed C≡N shift as Δν = 4.81 cm⁻¹ + 1.06 cm⁻¹/D·p, with the intercept taken as the DCT-induced shift and the slope fixed by two points (D3:A2 and D:A:MeCN). This assumes (i) the DCT-induced shift is exactly the same for the solvent-free 3:2 polymorph and all solvated 1:1 polymorphs, and (ii) the dipole-induced shift is strictly linear in the gas-phase dipole moment with zero intercept for zero dipole. Neither assumption is tested; the resulting constancy of the inferred DCT (0.11 e) is largely a consequence of the model structure. Please report uncertainties, use more than two data points to determine the dipole slope, and test the linearity assumption (e.g., with a nonpolar solvent or a wider range of dipole moments).","section":"Vibrational analysis and DCT (Eq. 2)"},{"comment":"The conversion of the 4.81 cm⁻¹ DCT shift to 0.11 e uses a single literature reference point (≈45 cm⁻¹ per electron from a crystal-violet HATCN anion salt). The raw C≡N shifts vary from 4.8 cm⁻¹ (D3:A2) to 9.0 cm⁻¹ (D:A:MeCN), and the measured values depend on peak fitting and on potential contributions from pristine phases (the authors note pristine DBTTF/HATCN in the MeCN-grown sample). An error estimate on the extracted DCT should be given, and the sensitivity of the result to the choice of reference calibration and to the two-point fit should be discussed.","section":"Raman calibration (Ref. 35 and Fig. 3)"}],"minor_comments":[{"comment":"The units in Eq. (2) should be written explicitly: Δν = 4.81 cm⁻¹ + 1.06 cm⁻¹·Debye⁻¹ × p.","section":"Eq. (2)"},{"comment":"The dashed red line labeled 'DCT-induced shift' is drawn as a constant; consider labeling it 'constant DCT contribution (4.81 cm⁻¹)' to avoid implying it is a fitted line.","section":"Fig. 3"},{"comment":"Typo in the diffraction software name: 'CrysAlisP rosoftware' should be 'CrysAlisPro software'.","section":"Experimental Methods"},{"comment":"The phrase 'non-invasive knob' in the conclusions is vague; consider replacing with a more precise description of the tunable structural/electronic parameter.","section":"Abstract/Conclusions"},{"comment":"The label 'DBTTF:HATCN:soln' in panel (a) is ambiguous; spell out 'THF2' or use the same notation as Table 1.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically competent and the experimental dataset is valuable. My main concern is that the abstract and conclusions state the solvent-as-primary-acceptor mechanism as a finding, while the only quantitative support is a single Hirshfeld partitioning. I would be satisfied if the authors either (a) add corroborating population analyses and reconcile the frontier-state picture, or (b) substantially reframe the claim as a partitioning-dependent observation that motivates further study. The Raman DCT value is reasonable but the decomposition needs error bars and more calibration points. This is a major revision rather than rejection because the claim is potentially correct and the underlying data are reproducible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: you should know this paper for the six new cocrystal structures and the careful structural characterization. The headline mechanism—solvent as primary electron acceptor—is not established by the evidence as presented; it rests on one Hirshfeld partition and sits in tension with the paper's own electronic structure results.\n\nWhat's genuinely new: six novel DBTTF:HATCN cocrystals, two stoichiometries (1:1 solvated with THF, ACE, MeCN; 3:2 solvent-free by horizontal vapor deposition), all in P2_1/c, with CIFs deposited (CCDC 2473715, 2475667–2475671) and data on Zenodo. The structural analysis is thorough: the correlations between solvent size, D:A separation, and HATCN dihedral distortions are documented. The optical result that CT1 sits near 1.4–1.5 eV across all forms is clean. The Raman fitting of the C≡N band is careful as far as it goes.\n\nThe soft spots. First, the quantitative DCT: Eq. (2) is a line through two points. The 'constant' 4.81 cm^-1 is just the zero-dipole point (the solvent-free 3:2), and the slope is fixed by MeCN. So the statement that the DCT is consistently 0.11 e is partly a consequence of the model, not the data; the raw shifts range from 4.8 to 9.0 cm^-1. The calibration to 45 cm^-1 per electron comes from a single crystal-violet HATCN anion reference. No error bars are given. This is a real weakness but not fatal, since the order of magnitude is plausible.\n\nSecond, and more important: Table 2 is the entire basis for the solvent-acceptor claim. Hirshfeld partitioning is scheme-dependent, and with contacts as short as 2.54 Å the overlapping densities are exactly where that ambiguity is largest. The paper's own PDOS (Fig. 5) shows solvent states only in the deeper valence region and CT1 assigned to DBTTF→HATCN. If THF truly held −0.35 e collectively, the frontier acceptor weight should move onto the solvent. The authors' distinction between charge accumulation and state-selective charge transfer does not resolve that contradiction. A Bader or Hirshfeld-I cross-check on D:A:THF2 would settle it. Without that, the headline mechanism is underdetermined.\n\nAlso minor: the D:A:MeCN sample contains pristine material, complicating that data point; and the paper itself calls the Raman shift an upper limit.\n\nWho this is for: people working on charge-transfer cocrystals or solvent inclusion in molecular crystals. They will want the structures and the optical data. The charge-transfer claims should be read as a testable hypothesis. A serious referee should engage with this paper—the structural work deserves it and the mechanism question is worth hashing out—but the authors should be asked to tone down the 'primary electron acceptor' language or add a second population analysis.","headline":"Six new cocrystal structures with a careful structural story; the solvent-as-acceptor mechanism is a hypothesis that outruns the present evidence.","tokens_in":19177,"tokens_out":4355,"would_cite":true,"duration_ms":47178,"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":"In solvated DBTTF:HATCN cocrystals, solvent molecules act as the primary electron acceptors, accumulating up to -0.35 e and even reversing HATCN's role to a net donor.","keywords":["charge transfer","cocrystals","DBTTF","HATCN","solvent inclusion","Raman spectroscopy","density functional theory","organic optoelectronics"],"falsifier":"Compute partial charges for the D:A:THF2 structure with an alternative partitioning method, such as Bader or electrostatic-potential fitting, and check whether THF still accumulates about -0.35 e; alternatively, measure the solvent's core-level binding energy by X-ray photoelectron spectroscopy and see whether it shifts as expected for negatively charged THF.","tokens_in":1252,"feed_emoji":"🧪","tokens_out":8118,"duration_ms":103973,"temperature":0.7,"pith_summary":"This paper asks whether residual solvent molecules trapped in donor-acceptor cocrystals are passive structural guests or active electronic participants. In six new DBTTF:HATCN cocrystals, the authors combine Raman spectroscopy with density-functional theory to show that overall charge transfer stays near 0.11 electrons across all structures, yet the partial-charge analysis reveals that in solvated crystals the solvent, not HATCN, is the main electron acceptor. In the THF-heavy structure, two THF molecules together hold about 0.35 electrons and HATCN actually donates charge. The claim matters because it turns solvent from a processing artifact into a design variable and complicates the standard donor-to-acceptor charge-transfer picture.","feed_headline":"Solvent molecules act as primary electron acceptors in cocrystals","feed_subtitle":"Trapped THF and acetone store up to 0.35 electrons, overturning the donor-acceptor charge picture.","key_machinery":"The analysis rests on two complementary probes. First, the Raman shift of HATCN's cyano (C≡N) stretching mode is used as a charge-transfer meter, decomposed into a charge-transfer contribution (4.81 cm^-1) and a solvent-dipole contribution (1.06 cm^-1 per Debye), calibrated to about 45 cm^-1 per electron from a reference HATCN anion salt to yield 0.11 e. Second, partial-charge partitioning on the experimental crystal structures assigns net charges to DBTTF, HATCN, and each solvent molecule, revealing where the electrons actually sit and showing that the solvent can become the dominant acceptor.","core_discovery":"The central discovery is that in solvated DBTTF:HATCN cocrystals the solvent molecules actively participate in charge transfer as primary electron acceptors. Raman measurements give a consistent degree of charge transfer of 0.11 e in every polymorph, and density-functional theory on solvent-free structures reproduces that value. But when the solvent is included in the simulation, the picture changes: in D:A:THF2 the two THF molecules accumulate -0.35 e total, DBTTF loses 0.26 e, and HATCN loses 0.09 e, making HATCN a net donor. A similar pattern appears with acetone. The authors conclude that charge distribution in these ternary crystals involves the whole atomic network, not just the cyano","pith_inferences":["If solvents genuinely store significant negative charge, then 'degree of charge transfer' measured from the acceptor's vibrational modes is an incomplete metric: it reports the acceptor's net change, not the destination of the electrons, so doping models for solvent-containing films may need revision.","A testable extension would be to vary the solvent's electron affinity, for example by using fluorinated or nitrile-containing analogues of THF and acetone, and watch whether the Raman shift and band gap respond as the solvent becomes a stronger or weaker acceptor.","The same reasoning suggests that other 'inert' intercalants, including water or halogenated solvents, might also participate electronically; residual-solvent characterization should become a standard check in donor-acceptor cocrystal studies."],"forward_implications":["Solution-grown 1:1 cocrystals preserve the intrinsic donor-acceptor electronic character: band gaps and absorption onsets stay near 1.4-1.5 eV regardless of solvent inclusion.","The degree of charge transfer measured by Raman is essentially constant, about 0.11 e, across solvent-free and solvated polymorphs, so the donor-to-acceptor transfer is consistent across processing routes.","In solvated crystals, charge can flow to the solvent from both donor and acceptor, so HATCN's role is not fixed; in D:A:THF2 it becomes a net electron donor.","The solvent dipole moment shifts the cyano vibration linearly, allowing the polarity effect and the charge-transfer effect to be separated in Raman data.","Increasing the number of solvent molecules widens the band gap by roughly 0.1 eV, indicating a subtle but systematic handle for tuning electronic properties."],"fun_headline_variants":["Solvent molecules become the primary electron acceptors in cocrystals","In solvated cocrystals, THF and acetone siphon electrons, not HATCN","Cocrystal solvents steal 0.35 e, overturning donor-acceptor roles","Solvent inclusion flips charge transfer: solvent accepts, donor loses","Trapped solvent acts as electron acceptor in DBTTF:HATCN cocrystals"],"cache_read_input_tokens":20480,"weakest_assumption_plain":"The load-bearing premise is that the computational partial-charge scheme correctly assigns electron density between the solvent molecules and the donor-acceptor stack; if that assignment is wrong, the claim that the solvent is the primary electron acceptor fails.","fun_headline_variants_meta":{"raw":{"variants":["Solvent molecules become the primary electron acceptors in cocrystals","In solvated cocrystals, THF and acetone siphon electrons, not HATCN","Cocrystal solvents steal 0.35 e, overturning donor-acceptor roles","Solvent inclusion flips charge transfer: solvent accepts, donor loses","Trapped solvent acts as electron acceptor in DBTTF:HATCN cocrystals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000823,"raw_usage":{"total_tokens":3464,"prompt_tokens":801,"completion_tokens":2663,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":2560}},"tokens_in":545,"tokens_out":2663,"duration_ms":22663,"temperature":1.0,"reasoning_tokens":2560,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:18:51.121722+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute partial charges for the D:A:THF2 structure with an alternative partitioning method, such as Bader or electrostatic-potential fitting, and check whether THF still accumulates about -0.35 e; alternatively, measure the solvent's core-level binding energy by X-ray photoelectron spectroscopy and see whether it shifts as expected for negatively charged THF.","supporting_citations":[],"review_version":1}