pith. sign in

arxiv: 2604.15930 · v1 · submitted 2026-04-17 · 🌌 astro-ph.GA

The Reaction between Atomic Carbon and Molecular Nitrogen as a Source of Cyanamide and Carbodiimide on Interstellar Ices

Kevin M. Hickson , Jean-Christophe Loison , Audrey Coutens This is my paper

Pith reviewed 2026-05-10 08:23 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords cyanamidecarbodiimideinterstellar icesatomic carbonmolecular nitrogensurface chemistryprotostellar environmentsastrochemical modeling
0
0 comments X

The pith

The ice-surface reaction of atomic carbon with molecular nitrogen is the dominant source of cyanamide and carbodiimide in dense clouds and protostellar environments.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines the reaction of ground-state atomic carbon with molecular nitrogen both in the gas phase and on amorphous solid water clusters to determine whether it can supply the missing formation routes for cyanamide and carbodiimide. Quantum calculations show that C atoms landing on N2-covered ice form CNN without a barrier, after which exothermic hydrogenation steps break the N-N bond through low-barrier cyclic intermediates to yield the N-C-N backbone of the target molecules. When these steps are added to a three-phase astrochemical model of low-mass protostar evolution, the surface channel accounts for the large majority of NH2CN and HNCNH production in both dense clouds and protostellar regions, resolving underpredictions in existing networks.

Core claim

The reaction of ground-state atomic carbon with molecular nitrogen already present on the surface of amorphous solid water produces CNN in a barrierless process. Successive exothermic hydrogenation of the resulting C-N-N intermediates allows N-N bond cleavage via cyclic structures with low barriers, generating molecules with N-C-N backbones. Insertion of this sequence into an updated three-phase astrochemical model shows that the ice-surface C + N2 route constitutes by far the dominant formation pathway for NH2CN and HNCNH in both dense interstellar clouds and low-mass protostellar environments.

What carries the argument

Barrierless formation of CNN on ASW followed by hydrogenation, N-N bond breaking, and low-barrier cyclic intermediates that rearrange to N-C-N structures.

If this is right

  • Astrochemical models that omit the C + N2 ice reaction will continue to underpredict NH2CN and HNCNH abundances in dense regions.
  • Surface production of N-C-N backbones from simple atomic carbon addition becomes the leading channel once N2 is present on the ice.
  • The cyclic intermediates identified in the calculations provide a general template for converting linear C-N-N species into branched nitrogen organics on ice.
  • Destruction rates of the hydrogenated C-N-N intermediates must remain low enough relative to rearrangement for the net yield to stay high.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar atomic-carbon additions to other simple molecules already adsorbed on ice could supply additional underpredicted organics in the same environments.
  • Direct infrared searches for CNN or its hydrogenated forms on interstellar ices would provide an independent test of the proposed sequence.
  • If the surface route dominates nitrogen-carbon coupling, the total budget of N-C-N species available for delivery to forming planets may be higher than gas-phase-only models suggest.

Load-bearing premise

The calculated barriers and exothermicities for CNN formation, hydrogenation, and N-N cleavage on small water clusters remain valid on real interstellar ice mantles and are not overturned by competing destruction channels or model uncertainties.

What would settle it

A laboratory measurement on realistic water-ice surfaces at 10 K that finds either a significant barrier to C + N2 addition or rapid destruction of the CNN intermediate before hydrogenation would falsify the dominance claim.

Figures

Figures reproduced from arXiv: 2604.15930 by Audrey Coutens, Jean-Christophe Loison, Kevin M. Hickson.

Figure 4
Figure 4. Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
read the original abstract

Reactions occurring on the ice-covered surfaces of interstellar dust grains are considered to be among the most important sources of complex species in the interstellar medium. Despite this, molecules such as cyanamide, NH2CN, are largely underpredicted by current astrochemical models suggesting that the network of reactions currently used to describe this species and its tautomer carbodiimide, HNCNH, are incomplete. Here, we performed a theoretical investigation of the reaction of ground state atomic carbon C(3P) with molecular nitrogen N2 in both the gas-phase and on the surface of amorphous solid water (ASW) clusters to examine its potential importance in the formation of NH2CN and HNCNH. We show that the reaction of gas-phase C-atoms with N2 molecules already present on the ASW surface results in the barrierless formation of CNN. Following exothermic hydrogenation reactions, the N-N bond of the C-N-N bearing intermediates is broken allowing the formation of molecules with N-C-N backbones through cyclic intermediates over low barriers. To test the importance of these processes to NH2CN and HNCNH formation, these reactions were included in a three-phase astrochemical model of low-mass protostellar evolution employing a reaction network that was updated to better describe the formation and destruction pathways of related small nitrogen bearing molecules. These simulations demonstrate that the ice surface reaction between C and N2 represents by far the dominant source of NH2CN and HNCNH in protostellar environments and in dense clouds.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper claims that the gas-phase C(3P) + N2 reaction on ASW ice surfaces proceeds barrierlessly to CNN, followed by exothermic hydrogenation, N-N bond cleavage via low-barrier cyclic intermediates, yielding NH2CN and HNCNH. When these surface reactions are added to an updated three-phase astrochemical network, the C + N2 ice pathway becomes by far the dominant source of both species in dense clouds and low-mass protostellar environments, addressing current model underpredictions.

Significance. If the barrierless character and dominance hold, the work supplies a previously missing efficient route to two N-bearing species whose abundances are systematically underpredicted, strengthening the case that grain-surface chemistry controls the inventory of small prebiotic molecules. The integration of explicit DFT cluster results into a full evolutionary model is a constructive step; the absence of free parameters in the core derivation is a positive feature.

major comments (2)
  1. [Computational Methods / Surface Reaction Results] Quantum chemical calculations on ASW clusters: the reported barrierless CNN formation and subsequent N-N cleavage rest on small clusters (typically <20 H2O). Edge effects and incomplete solvation can raise effective barriers by 1-5 kcal/mol relative to extended ice; without periodic-boundary or larger-cluster benchmarks, the assumption that these steps remain barrierless under realistic interstellar conditions is not yet load-bearing for the dominance claim.
  2. [Astrochemical Simulations] Astrochemical modeling section: the updated network is inserted into the three-phase code and the simulations then conclude that the new C + N2 surface route dominates. This introduces a circularity risk; without tabulated sensitivity runs that vary the rates of competing formation/destruction channels (e.g., CNN photodissociation or H2O reactions) or an explicit listing of all network changes, it is unclear whether dominance is an emergent result or an artifact of the chosen updates.
minor comments (3)
  1. [Abstract] The abstract states 'low-barrier' cyclic intermediates but supplies neither numerical barrier heights nor zero-point-corrected values; these should be reported explicitly with estimated uncertainties.
  2. [Computational Methods] No error bars or convergence tests are mentioned for the DFT energetics; adding a short table of basis-set and functional sensitivity would improve reproducibility.
  3. [Figures] Figure captions for the reaction schemes could explicitly label the ASW cluster size and the level of theory used for each stationary point.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives new reaction pathways via quantum chemical calculations on ASW clusters (barrierless CNN formation, exothermic hydrogenation, N-N cleavage via cyclic intermediates), then inserts the resulting rates into an updated three-phase astrochemical network and runs forward simulations of protostellar evolution to quantify relative contributions. This is a standard predictive modeling workflow: the dominance conclusion is an output of the simulation comparing production fluxes across all network pathways, not a re-expression of the input rates or a fit to the target abundances. No self-definitional steps, fitted parameters renamed as predictions, or load-bearing self-citations appear in the abstract or described chain; the network update is presented as an improvement to describe related N-bearing species, with the new C+N2 channel tested rather than presupposed.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on unverified quantum chemical accuracy for surface reactions and the assumption that the updated reaction network captures all relevant pathways without introducing bias toward the new source.

axioms (2)
  • domain assumption Quantum chemical calculations on amorphous solid water clusters reliably predict reaction barriers and products for C + N2 and subsequent hydrogenation steps
    The paper bases barrierless formation and low-barrier paths on these computations.
  • domain assumption The three-phase astrochemical model with updated nitrogen network accurately represents interstellar conditions without missing major formation or destruction routes
    Dominance conclusion depends on this model completeness.

pith-pipeline@v0.9.0 · 5586 in / 1394 out tokens · 78401 ms · 2026-05-10T08:23:22.795770+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

  1. [1]

    NH + H2CN → NH2 + HCN → H2CNN + H -273 -116 3.0e-11 3.0e-11 0 0 0 0 Radical-radical reaction

    work page

  2. [2]

    NH + H2NC → NH2 + HCN → HNCNH + H -398 -241 3.0e-11 3.0e-11 0 0 0 0 Radical-radical reaction

    work page

  3. [3]

    NH2 + H2CN→ NH3 + HCN → H + HCNNH2 -335 +161 6.0e-11 0 0 0 Deduced from (Yelle et al. 2010)

    work page 2010

  4. [4]

    NH2 + H2NC→ NH3 + HNC → H + H2NCNH -408 +42 6.0e-11 0 0 0 Radical-radical reaction

    work page

  5. [5]

    CN + NH2 → NH + HCN → H + HNCN -145 -101 1.0e-10 1.0e-10 0 0 0 0 Radical-radical reaction

    work page

  6. [6]

    2009, Meads et al

    CN + NH3 → HCN + NH2 -83 2.77e-11 -1.14 0 (Blitz et al. 2009, Meads et al. 1993, Talbi & Smith 2009, Sims et al. 1994) 7. CN + CH3NH2 → HCN + CH3NH → HCN + CH2NH2 → NH2CN + CH3 → CH3NHCN + H -116 -146 -151 -76 4.0e-10 0 0 0 0 0 Global rate from (Sleiman et al. 2018). Their branching ratios are derived from theoretical calculations for pathways with a barr...

    work page 2009

  7. [7]

    H + CNN → HCNN → HNNC → c-HNNC → c-HCNN → HNCN → CN + NH → CH + N2 -308 -295 -171 -283 -452 +27 -171 0.3 0 0 0 0.7 0 0 0 0 CNN is a triplet state with C=N=N configuration

    H + NCN → HNCN -344 1 0 Radical-radical reaction 9. H + CNN → HCNN → HNNC → c-HNNC → c-HCNN → HNCN → CN + NH → CH + N2 -308 -295 -171 -283 -452 +27 -171 0.3 0 0 0 0.7 0 0 0 0 CNN is a triplet state with C=N=N configuration. Both HCNN and HNNC should be produced. Some HCNN should isomerize into HNCN and most of the HNNC (this work). 10. H + HNCN → NH2CN ...

    work page

  8. [8]

    H + HCNN → H2CNN → CH2(a1A1) + N2 -410 -237 0.0 1 0 0 No exit TS on the singlet surface leading to CH2(a1A1) + N 2. C H2(a1A1) + N 2 can give back to S2 → CH2(X3B1) + N2 → HCNNH → HCN + NH(X3-) → NH2CN → HNCNH -291 -294 -188 -543 -538 0.4 9 0 0 0.3 0.2 0 0 H2CNN or relaxed to C H2(X3B1) + N 2 or C H2(a1A1) can react with H2O species nearby. CH2(X3B1) + N...

    work page

  9. [9]

    2016, Lu et al

    C + N2 → CNN → N + CN -124 +196 1 0 No barrier (Hickson et al. 2016, Lu et al. 2023, Urzúa-Leiva & Denis-Alpizar 2021)

    work page 2016

  10. [10]

    1996, Le Picard & Canosa 1998, Berman et al

    CH + N2 → HCNN -137 1 0 (Brownsword et al. 1996, Le Picard & Canosa 1998, Berman et al. 2007)

    work page 1996

  11. [11]

    CH3 + NH2CNH → NH2CN + CH4 -342 1 0 Radical-radical reaction, we neglect the adduct (H2N-C(CH3)=NH) formation. 22. CH3 + NH2CHNH2 → NH2CHNH + CH4 → NH3 + HNC + CH4 -326 -280 1 0 0 Radical-radical reaction, we neglect the adduct (H2N-CH(CH3)-NH2) formation. The TS for NH2CHNH→ NH3 + HNC is located 343 kJ/mol above the NH2CHNH energy. 23. N + CN → NCN → 3CN...

    work page

  12. [12]

    N + CH3NH → CH3 + N2 + H -282 1 0 Radical-radical reaction. 27. N + CH2NH2 → NCH2NH2 → H2CN + NH2 → NH2CHN + H → NH2CNH + H → NH2CHNH → HCN + NH3 → NH2CN + H2 -282 -189 -171 -146 -569 -523 -486 0 0 0 0 0 1 0 0 Calculations at M06-2X/AVTZ level. We favor HCN + NH3, the TS for NH2CHNH → HCN + NH 3 being localized -306 kJ/mol below the N + CH2NH2 energy. 28....

    work page

  13. [13]

    N + H2CNN → H2CN + N2 → H + HCN + N2 -485 -380 0.2 0.8 0 Radical-radical reaction

    N + HCNN → HCN + N2 -790 1 0 The attack of the nitrogen atom on C or terminal N leads in both cases to HCN + N2 30. N + H2CNN → H2CN + N2 → H + HCN + N2 -485 -380 0.2 0.8 0 Radical-radical reaction. Most of the H2CN should dissociate into H + HCN considering the exothermicity of the reaction

    work page

  14. [14]

    NH + H2CN → NH2 + HCN → H2CNNH → H2CNN + H → HCNNH + H → CH2 + N2 + H -273 -265 -116 0 +2 0.5 0.5 0 0 0 0 0 H2CNNH is likely produced but is not considered in the network

    NH + CN → HNCN → NCN + H -479 -136 1 0 0 0 32. NH + H2CN → NH2 + HCN → H2CNNH → H2CNN + H → HCNNH + H → CH2 + N2 + H -273 -265 -116 0 +2 0.5 0.5 0 0 0 0 0 H2CNNH is likely produced but is not considered in the network

    work page

  15. [15]

    NH + HCNH → NH2 + HCN → HNCNH + H -305 -277 0.5 0.5 0 0

    work page

  16. [16]

    NH + CH2NH2 → NH2 + CH2NH → HNC + NH3 + H -221 -143 0.5 0.5 0 0

    work page

  17. [17]

    NH2 + CN → NH2CN → HNCNH → HNCN + H -500 -496 -101 0.8 0.2 0 0 0 0 Some NH 2CN may isomerize into HNCNH (this work)

    NH + CH3NH → NH2 + CH2NH → N2 + H2 + CH3 -251 -380 0.5 0.5 0 0 36. NH2 + CN → NH2CN → HNCNH → HNCN + H -500 -496 -101 0.8 0.2 0 0 0 0 Some NH 2CN may isomerize into HNCNH (this work). 37. NH2 + H2CN → H2CNNH2 → NH3 + HCN -258 -335 0.5 0.5 0 0 H2CNNH2 is likely produced but is not considered in the network. The TS (H2CNNH2→ HCN + NH3) = -10 kJ/mol below th...

    work page

  18. [18]

    2009, Meads et al

    CN + NH3 → HCN + NH2 -83 1 0 (Blitz et al. 2009, Meads et al. 1993, Talbi & Smith 2009, Sims et al. 1994) 40. HCO + CNN → HCNN + CO → HNNC + CO → HNCN + CO → CNNCHO → NNCCHO -248 -235 -392 -230 -240 0.5 0 0.5 0 0 0 0 CNN is a triplet state with C=N=N configuration

    work page 2009

  19. [19]

    HCO + HCNN → H2CNN + CO → CH2 + N2 + CO → HCNNCHO → HCN + NCHO → HCN + HNCO → HCOCHNN -350 -232 -271 -138 -496 -374 0.0 1 0.4 9 0 0 0.5 0 0 0 S4 0 43

    HCO + NCN → HNCN + CO → NCNCHO -283 -268 1 0 0 42. HCO + HCNN → H2CNN + CO → CH2 + N2 + CO → HCNNCHO → HCN + NCHO → HCN + HNCO → HCOCHNN -350 -232 -271 -138 -496 -374 0.0 1 0.4 9 0 0 0.5 0 0 0 S4 0 43. HCO + HNCN → NH2CN + CO → HNCNH + CO → NCN + H2CO → HCONHCN -339 -334 -18 -359 0.8 0.2 0 0 0 0

    work page

  20. [20]

    HCO + NH2CN → HNCN + H2CO +37 0 Endothermic reaction

    work page

  21. [21]

    HCO + HNCNH → HNCN + H2CO +33 0 Endothermic reaction

    work page

  22. [22]

    HCO + NH2CNH → NH2CN + H2CO → NH2CHNH + CO -275 -363 0.5 0.5 0 0

    work page

  23. [23]

    CH2OH + CNN → HCNN + H2CO → HNNC + H2CO → HNCN + H2CO → C2H2 + OH + N2 -181 -168 -325 -264 0.5 0 0.5 0.5 CNN is a triplet state with C=N=N configuration

    HCO + NH2CHNH2 → NH2CHNH + H2CO → NH2CH2NH2 + CO -259 -295 0.5 0.5 0 0 48. CH2OH + CNN → HCNN + H2CO → HNNC + H2CO → HNCN + H2CO → C2H2 + OH + N2 -181 -168 -325 -264 0.5 0 0.5 0.5 CNN is a triplet state with C=N=N configuration

    work page

  24. [24]

    CH2OH + NCN → HNCN + H2CO -217 1 0

    work page

  25. [25]

    CH2OH + HNCN → NH2CN + H2CO → HNCNH + H2CO -272 -268 0.8 0.2 0 0

    CH2OH + HCNN → H2CNN + H2CO → CH2 + N2 + H2CO -283 -165 0.2 0.8 0 0 51. CH2OH + HNCN → NH2CN + H2CO → HNCNH + H2CO -272 -268 0.8 0.2 0 0

    work page

  26. [26]

    CH2OH + NH2CN → HNCN + CH3O +6 0 Endothermic reaction. 53. CH3O + CNN → HCNN + H2CO → HNNC + H2CO → HNCN + H2CO → N2 + CO + CH3 -215 -202 -359 -529 0 0 0.5 0.5 CNN is a triplet state with C=N=N configuration

    work page

  27. [27]

    CH3O + NCN → HNCN + H2CO -251 1 0

    work page

  28. [28]

    CH3O + HCNN → H2CNN + H2CO → CH2 + N2 + H2CO -317 -199 0.1 0.9 0 0

    work page

  29. [29]

    CH3O + HNCN → NH2CN + H2CO → HNCNH + H2CO -306 -302 0.8 0.2 0 0

    work page

  30. [30]

    CH3O + NH2CN → HNCN + CH3OH -28 1 1760 M06-2X/AVTZ calculations (this work)

    work page

  31. [31]

    CH3O + HNCNH → HNCN + CH3OH -32 1 1580 M06-2X/AVTZ calculations (this work) Berman M.R., Tsuchiya T., Gregusova A., Perera S.A., Bartlett R.J., 2007, J. Phys. Chem A., 111, 6894 Blitz M.A., Seakins P.W., Smith I.W.M., 2009, Physical Chemistry Chemical Physics, 11, 10824 Brownsword R.A., Herbert L.B., Smith I.W.M., Stewart D.W.A., 1996, J. Chem. Soc. Farad...

    work page 2007