{"id":"cc23663f-1070-4de0-9807-6260303f83a0","arxiv_id":"2412.19065","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"This study predicts the C1s X-ray absorption spectra and vibronic fine structure of CN+, CN, and CN-, assigning first peaks to sigma* (cation, neutral) and pi* (anion) transitions, with 0-0 energies of 280.7, 279.6, and 285.8 eV.","lead":"The paper predicts the carbon K-edge X-ray absorption spectra of the interstellar molecules CN+, CN, and CN- using quantum chemistry simulations. It assigns the peaks, explains why the anion absorbs more strongly, and calculates the vibrational fine structure, giving reference energies for astronomy and laboratory spectroscopy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported 0-0 energies (280.7/279.6/285.8 eV) are anchored to the BLYP spectrum via the Fig. 4 alignment, so the multiconfigurational wavepacket calculations do not independently validate the absolute positions; only the relative vibronic shapes are method-independent.","rationale":"This is a solid computational spectroscopy paper: the methods are standard, the CO+ benchmark reproduces the experimental peak separation well, and the qualitative physics—peak assignments, the degeneracy-based intensity argument for CN-, and core-hole bond-length contraction—are robust. The reader's weakest assumption correctly identifies δ transfer from CO+ as the main risk. My concern is sharper: the paper's own Fig. 4 caption states that all non-BLYP methods were aligned to the BLYP spectrum by adjusting the 0-0 peak. Therefore the apparent agreement among the three multiconfigurational methods in Fig. 4 is partly enforced, not emergent, and cannot be used as evidence for the absolute 0-0 energies. The abstract and conclusions quote 280.7, 279.6, and 285.8 eV as if they were the result of wavepacket calculations on multiconfigurational PECs; in the presented analysis, those numbers are anchored to BLYP plus the CO+-derived shift. This does not invalidate the paper's relative vibronic profiles or assignments, but it does mean the headline absolute energies are less independently supported than they appear. I would keep the CONDITIONAL verdict: the predictions are useful reference data, but the absolute energy scale should be treated as benchmark-dependent pending an independent calibration. The proposed test—reporting unaligned 0-0 energies for the three multiconfigurational methods—would settle how much of the reported agreement is real versus imposed by the alignment step.","tokens_in":20143,"tokens_out":9309,"duration_ms":92223,"concrete_test":"Recompute the RASSCF (Molpro), RASSCF (OpenMolcas), and RASPT2 vibronic 0-0 energies for CN+, CN, and CN- without the Fig. 4 alignment, applying only the method-specific δ calibrated on CO+ from Section II B. If the three unaligned 0-0 positions for a given species spread by more than about 0.3 eV, then the headline 280.7/279.6/285.8 eV values are set by the BLYP anchor rather than independently predicted by the multiconfigurational methods.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the set of 0-0 absorption energies. Section II B determines a method-specific shift δ by aligning CO+ to experiment and assumes it transfers unchanged to CN+, CN, and CN-. Section III E and the Fig. 4 caption make an additional, stronger assumption: \"all other methods were aligned to the BLYP spectrum by adjusting the 0-0 peak.\" Consequently, the RASSCF and RASPT2 curves in Fig. 4 agree with each other and with BLYP only in relative shape (0-1/0-2 spacings and intensities), not in absolute energy. The 0-0 positions quoted in the abstract and Section III E are therefore effectively BLYP-FCH positions shifted by the CO+-derived δ; they are not independently \"revealed\" by the anharmonic wavepacket calculations on multiconfigurational PECs. This matters because CN- has very different unoccupied orbitals and required a basis-set compromise for its wide-range spectrum, making the transfer of δ especially questionable for the anion. Without an independent calibration or an uncertainty estimate, the claimed accuracy of the three absolute energies rests on a single DFT functional plus one isoelectronic benchmark.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper predicts C1s X-ray absorption spectra for the astrochemically relevant species CN+, CN, and CN- using BLYP-FCH DFT, SA-RASSCF (in both Molpro and OpenMolcas), and MS-RASPT2. Each method is calibrated by a single ad hoc shift delta derived from the CO+ experimental spectrum, under the assumption that delta transfers to the carbon-nitrogen species. The authors assign the first absorption peak to C1s -> sigma* for CN+ and CN and to C1s -> pi* for CN-, explain the larger intensity in CN- by the twofold degeneracy of pi*, and compute vibronic fine structure with a quantum wavepacket method on anharmonic Morse-fitted multiconfigurational PECs. They report 0-0 absorption energies of 280.7 eV (CN+), 279.6 eV (CN), and 285.8 eV (CN-), together with core-hole-induced bond shortening and vibrational frequency increases.","tokens_in":20381,"tokens_out":7760,"duration_ms":66651,"significance":"If the absolute transition energies were independently validated, this would be a valuable reference set for X-ray astronomy and laboratory spectroscopy of interstellar CN species. The paper has clear strengths: it cross-checks two independent RASSCF implementations and RASPT2, reproduces the CO+ experimental peak separation at the multiconfigurational level, uses anharmonic wavepacket dynamics rather than only harmonic Franck-Condon factors, and supplies raw PEC data in the Supplementary Material. However, the headline absolute energies are not parameter-free: they depend on a single fitted shift and, in the final spectra, on an alignment of all methods to the BLYP spectrum. The relative vibronic profiles and the peak assignments are considerably more robust than the absolute energy calibration.","major_comments":[{"comment":"The reported 0-0 energies of 280.7, 279.6, and 285.8 eV are not independently determined by the multiconfigurational wavepacket calculations. The Fig. 4 caption states that \"all other methods were aligned to the BLYP spectrum by adjusting the 0-0 peak,\" and Section II B specifies that the BLYP spectrum is itself shifted by a CO+-derived delta. Consequently, the RASSCF and RASPT2 curves in Fig. 4 agree with BLYP only in relative vibronic shape, not in absolute energy. The abstract and Section III E present the three values as if they were revealed by the anharmonic wavepacket simulations on multiconfigurational PECs. Because Table I shows raw vertical excitation energies that differ by more than 1 eV between methods (for example, CN+ peak i at 282.22 eV for RASSCF/MOLPRO versus 281.01 eV for BLYP in the calibrated table), this alignment is not a negligible detail. I request either an independent calibration of delta for at least one additional species or an explicit uncertainty estimate for the 0-0 positions based on the spread of raw method-dependent energies; the current presentation overstates the evidence for the absolute values.","section":"Section III E and Fig. 4 caption"},{"comment":"The transferability of the method-specific shift delta from CO+ to CN+, CN, and CN- is assumed without quantitative support. CO+ is isoelectronic with CN, but not with CN+ or CN-, and Section III A explicitly notes that CN- has \"notably different unoccupied orbitals\" compared with the other systems. Core-hole relaxation, and hence the systematic error in the absolute core-excitation energy, may differ substantially for the anion. The authors should provide a test of transferability, for example by calibrating delta against an additional benchmark molecule or ion with a known C1s spectrum and by comparing unshifted RASPT2 vertical energies with BLYP across the three target species. Without such a test, the claimed accuracy of the absolute 0-0 energies, especially the CN- value of 285.8 eV, rests on a single fitted constant.","section":"Section II B"},{"comment":"For CN-, the wide-range electronic XAS spectrum was obtained with the smaller cc-pVDZ basis set and calibrated to the main aug-cc-pVTZ peak, because the aug-cc-pVTZ calculation covered only a narrow energy window. This compromise is acknowledged, but its effect on the relative intensities of the pi* first peak and the higher-energy C-np peaks is not quantified. Since the twofold degeneracy of pi* is used to explain the pronounced first absorption of CN-, it is important to show that the two-fold intensity enhancement is not an artifact of the smaller basis or of the calibration procedure. I ask the authors to state explicitly how the basis-set compromise affects the oscillator strengths, the relative peak intensities, and the value of delta for CN-, and to show any available aug-cc-pVTZ data over the full displayed range or to provide a conservative uncertainty estimate for the affected spectral region.","section":"Section II B and Fig. S1"}],"minor_comments":[{"comment":"The manuscript's species list is inconsistent: Table I and Fig. 2 include CO and CO+, while the abstract and main text focus on CN+, CN, and CN-; please clarify in the captions and text whether CO is used only as an additional reference or is an intended target.","section":"Table I and Fig. 2"},{"comment":"The statement that the ground-state PECs \"decrease in the order of CN+, CN, and CN-\" is not directly supported by the Te values in Table III (for example, RASPT2 Te is 0.082 eV for both CN+ and CN); please specify the reference energy for Te so that the ordering can be checked.","section":"Section III D and Table III"},{"comment":"The label \"non-physical shake-off processes\" for the negative IPs of the core-excited CN- is confusing: a negative vertical electron-detachment energy for a core-excited anion can be a physically meaningful autoionization threshold; please rephrase and clarify the interpretation.","section":"Table IV"},{"comment":"The paragraph beginning \"Multiconfigurational PECs for CN+/CN/CN- were simulated at varying bond distances...\" appears twice with only minor wording differences; please remove the duplicate.","section":"Section II C"},{"comment":"Please justify the choice of 0.05 eV HWHM lifetime broadening for the vibronic spectra and 0.4 eV HWHM for the electronic spectra, and specify whether the 0.4 eV includes instrumental broadening or only lifetime effects.","section":"Section II B and Section II D"}],"recommendation":"major_revision","confidential_remarks":"The core spectroscopy and assignments are likely sound, and the authors are transparent about the delta assumption, which is a point in their favor. The main problem is that the abstract and Section III E present the absolute 0-0 energies as multiconfigurational predictions when they are in fact anchored to the BLYP spectrum through the Fig. 4 alignment. This is fixable by rephrasing the claims, adding an uncertainty estimate, and testing the transferability of delta on one more benchmark. I see no indication of citation manipulation; self-citations are to the group's own methodological software and prior benchmarks and are appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first vibrationally resolved C1s XAS study for CN+, CN, and CN-, with consistent assignments and a nice explanation for the anion's intense first peak. The headline 0-0 energies (280.7, 279.6, 285.8 eV) are less independent than the abstract suggests: the multiconfigurational curves in Fig. 4 were aligned to the BLYP spectrum, so the absolute positions ultimately rest on BLYP plus the CO+-derived shift δ. The relative vibronic structure is method-independent and trustworthy.\n\nWhat's genuinely new: the specific spectra, peak assignments, anharmonic PECs, and vibronic fine structure for these three species. The degeneracy argument for CN- (two π* orbitals vs one σ*) is clean and explains the factor-of-two intensity difference. The CO+ benchmark reproduces the experimental peak separation well for all three multiconfigurational approaches, and the RASSCF results are consistent across two independent codes. Raw PECs are in the SI, which is good practice.\n\nSoft spots, in proportion: the largest is the transferability of δ. It is fitted on CO+ and assumed to hold for CN+, CN, and CN-, including the anion with its very different unoccupied orbital manifold and a basis-set compromise for the wide-range spectrum. No uncertainty is attached to the three absolute energies. That makes them conditional predictions, not definitive reference values. The vibronic intensity ratios and relative spacings are on much firmer ground. Minor: the captions are transparent about the alignment, but the abstract presents 0-0 energies without that qualification, so an unwary reader might over-credit the multiconfigurational methods for the absolute positions. The self-citations to the group's own codes are appropriate given the methods used.\n\nBottom line: this is a solid, honest computational spectroscopy paper that fills a real gap. It deserves serious peer review. I would ask the authors to state the provenance of the absolute energies more explicitly and, if possible, add an error estimate by comparing a second calibration species or by reporting the raw method-dependent shifts. The paper will be useful to anyone interpreting future X-ray observations or planning ion spectroscopy experiments.","headline":"Solid reference calculations for CN+/CN/CN- C1s XAS, but the headline 0-0 energies are anchored to BLYP via a CO+ shift rather than independently validated by the multiconfigurational wavepacket work.","tokens_in":20928,"tokens_out":3740,"would_cite":true,"duration_ms":35240,"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":"Predicted carbon K-edge spectra for CN+, CN, and CN- place the first vibronic lines at 280.7, 279.6, and 285.8 eV, with the CN- absorption assigned to a doubly degenerate π* transition.","keywords":["X-ray absorption spectroscopy","vibronic fine structure","core-hole effect","CN radical","interstellar molecules","multiconfigurational methods","quantum wavepacket method","charge-state effects"],"falsifier":"A gas-phase carbon K-edge X-ray absorption measurement on mass-selected CN+, CN, and CN- would settle the claim by checking the 0-0 peak positions against 280.7, 279.6, and 285.8 eV and by checking that the first CN- peak is about twice as strong as the first peaks of CN+ and CN. Disagreement beyond the calibrated method accuracy, or a wrong ordering of the first peaks, would falsify the calibrated predictions.","tokens_in":19947,"feed_emoji":"⚛️","tokens_out":6466,"duration_ms":62555,"temperature":0.7,"pith_summary":"The paper tries to establish what high-resolution carbon K-edge X-ray absorption spectra of CN+, CN, and CN- should look like, since no experimental spectra exist for these interstellar species. It argues that the first electronic peak is a C1s→σ* transition in CN+ and CN but a C1s→π* transition in CN-, and that the two-fold degeneracy of the π* orbitals makes CN-'s first band roughly twice as intense. Using multiconfigurational potential energy curves and a quantum wavepacket treatment, it predicts the 0-0 vibronic lines at 280.7, 279.6, and 285.8 eV, each profile dominated by the 0-0 peak. These numbers are anchored by calibrating each method to the measured CO+ spectrum, because CO+ is isoelectronic with CN.","feed_headline":"First X-ray peaks predicted for CN+, CN, and CN-","feed_subtitle":"Calculations place the 0-0 absorptions at 280.7, 279.6, and 285.8 eV and assign CN- to a degenerate π* transition.","key_machinery":"The carrying object is a set of anharmonic potential energy curves for the ground and lowest C1s core-excited states of each species, computed at the multiconfigurational level (SA-RASSCF and MS-RASPT2) and fitted to Morse functions; the spectra are obtained by propagating a quantum wavepacket on these curves. A second load-bearing piece is the method-specific ad hoc shift δ, determined by aligning the calculated CO+ spectrum to its measured first peak and then applied, unchanged, to CN+, CN, and CN-. The Morse fits supply equilibrium distances, well depths, and force constants, and the wavepacket propagation converts those curves into the 0-0, 0-1, and 0-2 vibronic stick intensities.","core_discovery":"On the paper's own terms, the central discovery is a set of quantitative predictions, not a measurement: vibrationally resolved C1s XAS profiles for the three charge states of CN, with assignments and peak energies. The first electronic absorption of CN+ and CN is C1s→σ*, that of CN- is C1s→π*, and the degeneracy of the π*xz and π*yz channels doubles the anion's first-peak intensity; the 0-0 energies are 280.7 eV (CN+), 279.6 eV (CN), and 285.8 eV (CN-). The calculations also show that creating a C1s core hole shortens the C-N bond by roughly 4-6 pm (6-7 pm at the DFT harmonic level) and raises the vibrational frequency by 300-400 cm−1 in all three species. The authors take the agreement among RASSCF, RASPT2, and DFT profiles, plus the calibration to CO+, as evidence that these are accurate predictions for astronomy and laboratory work.","pith_inferences":["The transferability of the CO+ calibration shift is the natural thing to test: because CN- has diffuse orbitals and strong relaxation, its absolute energies may carry a larger systematic error than the 0-1 relative spacings; a high-level calculation or measurement of the anion alone would settle this.","If the predicted lines are confirmed, the X-ray absorption bands could be used as diagnostics of charge state in astrophysical environments where CN+, CN, and CN- coexist, paralleling how their microwave emission is already used.","The same wavepacket-plus-Morse machinery could be applied to other isoelectronic triads to test whether the non-monotonic energy ordering seen here is a general feature of changing charge state."],"forward_implications":["If the predictions hold, the first carbon K-edge absorption of CN- should appear about 6 eV above those of CN+ and CN, giving observers a spectroscopic handle to separate the anion from the neutral and cation.","The assignment of the first peak to C1s→π* in CN- means its absorption should be roughly twice as intense as the corresponding C1s→σ* peaks, so intensity ratios can be used alongside energies for identification.","Each vibronic profile is dominated by the 0-0 line, so future high-resolution X-ray spectra should show a single prominent line at the predicted energy rather than a long vibrational progression.","The predicted bond shortening and frequency increase upon core excitation imply that vibrationally resolved spectra of all three species will shift intensity from higher vibrational levels toward the 0-0 line."],"supporting_citations":[{"why":"Provides the experimental CO+ carbon K-edge NEXAFS spectrum whose first peak fixes the ad hoc shift δ for every method.","marker":"[15]"},{"why":"Earlier validated simulation protocol for vibrationally resolved XAS of diatomics; establishes the anharmonic wavepacket treatment and the sensitivity of such spectra to anharmonicity.","marker":"[11]"},{"why":"Defines the full core hole approximation used for the DFT XAS calculations.","marker":"[57]"},{"why":"Defines the SA-RASSCF multiconfigurational method used to compute ground and core-excited potential energy curves.","marker":"[45–47]"},{"why":"Defines the MS-RASPT2 method used to add dynamic correlation to the potential energy curves and spectra.","marker":"[48, 49]"},{"why":"Supplies the Morse potential form used to fit all potential energy curves, producing the equilibrium distances, well depths, and width parameters that drive the wavepacket spectra.","marker":"[71]"}],"fun_headline_variants":["Charge state shifts CN X-ray absorption peaks","Predicted XAS reveals σ* vs π* in CN species","Degenerate π* doubles CN- first X-ray peak","Three CN charge states: distinct X-ray spectra","Quantum chemistry pins CN, CN± X-ray lines"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole absolute-energy scale rests on assuming that the ad hoc shift δ fitted to CO+ transfers unchanged to CN+, CN, and CN-; if the C1s core hole relaxes differently in these charge states, all predicted 0-0 energies would be off by roughly the same systematic amount while the spectral shapes would remain useful.","fun_headline_variants_meta":{"raw":{"variants":["Charge state shifts CN X-ray absorption peaks","Predicted XAS reveals σ* vs π* in CN species","Degenerate π* doubles CN- first X-ray peak","Three CN charge states: distinct X-ray spectra","Quantum chemistry pins CN, CN± X-ray lines"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000695,"raw_usage":{"total_tokens":3227,"prompt_tokens":1116,"completion_tokens":2111,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":732,"completion_tokens_details":{"reasoning_tokens":2034}},"tokens_in":732,"tokens_out":2111,"duration_ms":16180,"temperature":1.0,"reasoning_tokens":2034,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:57:38.626897+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A gas-phase carbon K-edge X-ray absorption measurement on mass-selected CN+, CN, and CN- would settle the claim by checking the 0-0 peak positions against 280.7, 279.6, and 285.8 eV and by checking that the first CN- peak is about twice as strong as the first peaks of CN+ and CN. Disagreement beyond the calibrated method accuracy, or a wrong ordering of the first peaks, would falsify the calibrated predictions.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental CO+ carbon K-edge NEXAFS spectrum whose first peak fixes the ad hoc shift δ for every method."},{"cited_title":"Zhang, M","cited_arxiv_id":null,"evidence_quote":"Earlier validated simulation protocol for vibrationally resolved XAS of diatomics; establishes the anharmonic wavepacket treatment and the sensitivity of such spectra to anharmonicity."},{"cited_title":"Triguero, L","cited_arxiv_id":null,"evidence_quote":"Defines the full core hole approximation used for the DFT XAS calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Morse potential form used to fit all potential energy curves, producing the equilibrium distances, well depths, and width parameters that drive the wavepacket spectra."}],"review_version":1}