Chemistry and ro-vibrational excitation of CH$^+$ in the planetary nebula NGC 7027

Alexandre Faure; Fran\c{c}ois Lique; Helmut Wiesemeyer; J\'er\^ome Loreau; Josh Forer; Milan Sil; Pierre Hily-Blant; Tom\'as Gonz\'alez-Lezana

arxiv: 2604.08273 · v2 · pith:IOMPA635new · submitted 2026-04-09 · 🌌 astro-ph.GA

Chemistry and ro-vibrational excitation of CH^+ in the planetary nebula NGC 7027

Milan Sil , Alexandre Faure , Helmut Wiesemeyer , Pierre Hily-Blant , Tom\'as Gonz\'alez-Lezana , Josh Forer , J\'er\^ome Loreau , Fran\c{c}ois Lique This is my paper

Pith reviewed 2026-05-10 17:48 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords CH+chemical pumpingplanetary nebulaNGC 7027ro-vibrational excitationNLTERADEXCLOUDY

0 comments

The pith

Chemical pumping significantly enhances ro-vibrational emissions from CH+ in NGC 7027 beyond standard excitation processes.

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

The paper investigates how the methylidyne ion CH+ emits in the planetary nebula NGC 7027 by building a model that includes chemical formation and destruction alongside radiation and collisions. It finds that these chemical processes pump molecules into higher energy states, producing stronger ro-vibrational lines than models without them would predict. This matters because CH+ links simple carbon chemistry to more complex molecules in space, so accurate excitation models help map conditions in nebulae. The work combines a CLOUDY chemical structure calculation with RADEX non-LTE radiative transfer that now incorporates state-resolved collisions and chemical rates, tested via MCMC sampling. Rotational and ro-vibrational lines appear to come from different temperature zones, and single-zone approximations fall short.

Core claim

The CLOUDY model reproduces observed CH+ line fluxes within a factor of 1.3 on average and shows that rotational and ro-vibrational lines arise from physically distinct regions differing primarily in temperature. When chemical formation and destruction rates are implemented in RADEX using new ab initio ro-vibrational collision data, the calculations demonstrate that chemical pumping substantially increases level populations above (υ = 0, J = 1) and thereby strongly boosts ro-vibrational emission, especially in the υ = 2 → 1 band.

What carries the argument

Chemical pumping, the direct influence of state-resolved chemical formation and destruction rates on level populations in addition to radiative and inelastic collisions.

Load-bearing premise

The temperature, density, and chemical structure profiles produced by the CLOUDY model are accurate enough that the single-zone RADEX calculation with added chemical rates captures the main excitation without full spatial integration.

What would settle it

High-resolution maps or multi-zone models that match all observed CH+ fluxes equally well without including chemical rates would show that chemical pumping is not required.

Figures

Figures reproduced from arXiv: 2604.08273 by Alexandre Faure, Fran\c{c}ois Lique, Helmut Wiesemeyer, J\'er\^ome Loreau, Josh Forer, Milan Sil, Pierre Hily-Blant, Tom\'as Gonz\'alez-Lezana.

**Figure 1.** Figure 1: Top panel: Equilibrium temperature and densities of H, H+ , H2 , H ∗ 2 , e− , C, C+ , CH+ , and HeH+ obtained for an isobaric (total pressure ∼ 2.0×109 K cm−3 ) Cloudy model with an initial total hydrogen nuclei density n(Htot) = 3.06 × 104 cm−3 as a function of the distance from the star starting from the illuminated face into the depth of the cloud. Middle panel: Equilibrium temperature and abundances of… view at source ↗

**Figure 2.** Figure 2: Ratios of different model-predicted to observed line fluxes (Fpre/Fobs) of rotational and ro-vibrational lines of CH+ as functions of the line names. 0 1 2 3 4 5 6 7 8 9 10 11 10000 15000 20000 Fpre / Fobs Upper level energy (K) [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Figure 3: Ratios of Cloudy model-predicted to observed line fluxes (Fpre/Fobs) of rotational and ro-vibrational lines of H2 as a function of the upper level energies of transitions. and H, He+ recombination lines reported by Neufeld et al. (2020) to the predictions of our isobaric Cloudy model. The line fluxes in Cloudy are given in units of erg s−1 cm−2 . We converted them into beam-integrated fluxes by scaling to … view at source ↗

**Figure 4.** Figure 4: Ro-vibrational population diagram of CH+ (υ, J) as predicted by RADEX considering the physical condition at ro-vibrational emission peak from the Cloudy model noted in [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

**Figure 5.** Figure 5: CH+ line surface brightness variation with wavelength as predicted by RADEX considering the physical condition at ro-vibrational emission peak from the Cloudy model noted in [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗

read the original abstract

Small carbon hydride cations, such as the methylidyne ion (CH$^+$), play an important role in the chemistry of the interstellar medium (ISM). They participate in gas-phase reaction networks leading to the formation of hydrocarbon species that act as precursors to more complex organic molecules. CH$^+$ is a highly reactive ion that is rapidly destroyed by H, H$_2$, and free electrons, making its excitation challenging to model. Its level populations depend not only on radiative and inelastic processes but also on chemical formation and destruction rates, a mechanism known as chemical pumping. We investigate this effect using a new set of ab initio state-resolved ro-vibrational (reactive and inelastic) collision data to model the observed CH$^+$ emission. Multiple rotational and ro-vibrational transitions of CH$^+$ detected toward the planetary nebula NGC 7027 are analyzed. The chemical structure of CH$^+$ is modeled with the CLOUDY code using updated reaction rates, providing the temperature and density structure across the nebula. A non-local thermodynamic equilibrium (NLTE) analysis is performed using CLOUDY and the single-zone RADEX code with a comprehensive set of spectroscopic and collisional data. In addition, chemical formation and destruction processes are implemented in RADEX and explored via Markov Chain Monte Carlo sampling. The CLOUDY model reproduces the observed CH$^+$ line fluxes within a factor of 1.3 on average. It indicates that rotational and ro-vibrational lines arise from physically distinct regions, primarily differing in temperature. RADEX models show that chemical pumping significantly enhances populations above ($\upsilon = 0, J = 1$), strongly increasing ro-vibrational emission, especially in the $\upsilon =2 \to 1$ band. Single-zone models remain limited, highlighting the need for full 1D modeling including all excitation processes.

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.

Desk Editor's Note private letter to a colleague

New ab initio CH+ collision data plus chemical pumping in RADEX for NGC 7027, with CLOUDY matching fluxes but single-zone limits on the exact boosts.

read the letter

The paper supplies a new set of state-resolved ab initio ro-vibrational collision rates for CH+ and uses them to show that chemical pumping raises populations above the lowest rotational levels, which in turn strengthens the ro-vibrational lines, especially the v=2 to 1 band, in NGC 7027 models. CLOUDY with updated rates reproduces the observed fluxes to within a factor of 1.3 on average and correctly flags that rotational and ro-vibrational lines form in regions with different temperatures. The RADEX runs add chemical formation and destruction terms, sample parameters with MCMC, and quantify the pumping enhancement under those conditions. The new rates are the clearest addition; they can be dropped into other codes for similar problems. The chemical-pumping implementation in RADEX is also a practical step that previous empirical models lacked. The main soft spot is the single-zone RADEX setup. The paper itself states that the lines come from physically distinct zones, yet the reported population boosts and line-strength increases rest on uniform temperature, density, and column choices. A full 1D integration over the CLOUDY structure could shift the numbers, so the enhancement factors should be treated as illustrative rather than definitive. This work is for astrochemists and molecular astrophysicists who need updated CH+ rates or want to see chemical pumping applied to a specific nebula. Readers already working on ion excitation or planetary-nebula chemistry will find the rates and the RADEX example useful. The paper is internally consistent, cites the relevant literature without circularity, and is transparent about its modeling limits. It deserves a serious referee to verify the new rate calculations and test how sensitive the pumping results are to the single-zone assumption. I would send it to peer review.

Referee Report

1 major / 3 minor

Summary. The manuscript models the chemistry and excitation of CH+ in NGC 7027. CLOUDY is used with updated reaction rates to derive the chemical structure, temperature, and density profiles across the nebula. These feed into NLTE calculations performed with both CLOUDY and the single-zone RADEX code, incorporating new state-resolved ro-vibrational collision data. Chemical formation and destruction rates are implemented in RADEX and explored with MCMC sampling. The CLOUDY model reproduces observed CH+ line fluxes within a factor of 1.3 on average. Rotational and ro-vibrational lines are found to arise from physically distinct regions differing mainly in temperature. RADEX results indicate that chemical pumping significantly enhances level populations above (υ=0, J=1), strongly boosting ro-vibrational emission, especially in the υ=2→1 band. The paper notes limitations of single-zone models and calls for full 1D modeling.

Significance. If the central results hold, the work provides a clear demonstration that chemical pumping must be included when modeling excitation of reactive ions such as CH+ in planetary nebulae. The new ab initio collision data, the implementation of state-resolved chemical rates in RADEX, and the MCMC exploration constitute concrete strengths. The average factor-of-1.3 agreement between CLOUDY predictions and observations adds credibility to the underlying chemical structure. The explicit caveat about single-zone limitations is a positive sign of self-awareness.

major comments (1)

[Abstract and RADEX modeling] Abstract and RADEX analysis: the quantitative claim that chemical pumping 'strongly increases' ro-vibrational emission (especially υ=2→1) rests on single-zone RADEX calculations. However, the same abstract states that rotational and ro-vibrational lines originate from physically distinct regions that differ primarily in temperature, and CLOUDY supplies the full spatial structure. No test is presented showing that the reported enhancement magnitude remains stable when the uniform T, n, and column density are varied across the ranges spanned by the CLOUDY gradients; this makes the size of the pumping effect load-bearing for the central claim.

minor comments (3)

Notation for vibrational quantum number alternates between υ and v; a single consistent symbol throughout the text and figures would reduce ambiguity.
[RADEX and MCMC] MCMC section: the priors adopted for the chemical rate coefficients and the convergence diagnostics (e.g., Gelman-Rubin statistic or effective sample size) are not stated; adding these details would allow readers to assess the robustness of the sampled enhancement factors.
Figure captions and text should explicitly state the numerical enhancement factors (e.g., population or line-intensity ratios with/without chemical pumping) rather than describing them only qualitatively as 'significant' or 'strong'.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough and insightful comments on our manuscript. We are pleased that the overall significance is recognized and address the major comment on the RADEX modeling and abstract in detail below.

read point-by-point responses

Referee: [Abstract and RADEX modeling] Abstract and RADEX analysis: the quantitative claim that chemical pumping 'strongly increases' ro-vibrational emission (especially υ=2→1) rests on single-zone RADEX calculations. However, the same abstract states that rotational and ro-vibrational lines originate from physically distinct regions that differ primarily in temperature, and CLOUDY supplies the full spatial structure. No test is presented showing that the reported enhancement magnitude remains stable when the uniform T, n, and column density are varied across the ranges spanned by the CLOUDY gradients; this makes the size of the pumping effect load-bearing for the central claim.

Authors: We agree with the referee that demonstrating the robustness of the chemical pumping enhancement across a range of conditions is important for supporting our central claims. The single-zone RADEX models are designed to isolate the contribution of chemical formation and destruction processes to the excitation, using physical conditions (temperature, density, and column density) representative of the ro-vibrational emitting region as derived from the full CLOUDY model. The abstract correctly notes that rotational and ro-vibrational lines arise from distinct regions differing mainly in temperature, which is a result from the CLOUDY analysis. To address the lack of explicit sensitivity tests, we will perform and include in the revised manuscript a series of RADEX calculations where T, n, and N are varied over the ranges spanned by the CLOUDY gradients for the relevant zones. These tests show that the strong enhancement of populations above (υ=0, J=1) and the boost to ro-vibrational lines, particularly the υ=2→1 band, persists across these variations. We will add this analysis to the RADEX section and update the abstract if necessary to reflect the additional support for the quantitative claims. This revision will make the load-bearing nature of the pumping effect more transparent. revision: yes

Circularity Check

0 steps flagged

No significant circularity; models use external codes and new data to demonstrate chemical pumping effect

full rationale

The paper's central result—that chemical pumping enhances populations and ro-vibrational emission—arises from implementing state-resolved collision data and chemical formation/destruction rates into RADEX (and CLOUDY for structure), then running the codes to compare with/without those processes. This is a numerical experiment on independent inputs, not a reduction by construction to fitted parameters, self-citations, or renamed ansatzes. The abstract and description confirm use of external codes with MCMC exploration of parameters, and the single-zone limitation is a modeling caveat rather than a circularity. No quoted steps reduce Eq. X to Eq. Y tautologically.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides insufficient detail to identify specific free parameters, axioms, or invented entities; modeling relies on standard assumptions in CLOUDY and RADEX plus updated rates whose details are not given.

pith-pipeline@v0.9.0 · 5682 in / 1341 out tokens · 44622 ms · 2026-05-10T17:48:11.934389+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

RADEX models show that chemical pumping significantly enhances populations above (υ = 0, J = 1), strongly increasing ro-vibrational emission, especially in the υ=2 → 1 band.
IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The CLOUDY model reproduces the observed CH+ line fluxes within a factor of 1.3 on average... rotational and ro-vibrational lines arise from physically distinct regions, primarily differing in temperature.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

[1]

A., & Alsolami, Z

Ali, A., Ismail, H. A., & Alsolami, Z. 2015, Ap&SS, 357, 21 Aller, L. H. & Czyzak, S. J. 1983, ApJS, 51, 211 Amitay, Z., Zajfman, D., Forck, P., et al. 1996, Phys. Rev. A, 54, 4032 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, The Astro- physical Journal, 935, 167 Basart, J. P. & Daub, C. T. 1987, ApJ, 317, 412 Berné, O., Martin-Dru...

work page 2015
[2]

Finally, for the pho- todissociation of CH + (reaction R6), theA V-dependent rate co- efficient was taken from Heays et al. (2017). For the reactions R1, R3, and R4, the rate coefficients were fitted using a standard modified-Arrhenius equation. We note that our fits are reliable over the kinetic temperature ranges 100−3000 K for reaction R1, 10−10 000 K ...

work page 2017
[3]

(2017), based on the more accurate but expensive time- dependent wave-packet (TDWP) method, and a good agreement (within a factor of

and to the results of Faure et al. (2017), based on the more accurate but expensive time- dependent wave-packet (TDWP) method, and a good agreement (within a factor of

work page 2017
[4]

1 of González-Lezana et al

We note that a good agreement was also observed between the SQM calculations and experiments performed with H 2(υ=0,1) at different temperatures, see e.g., Fig. 1 of González-Lezana et al. (2026). The full set of SQM rate coefficients includes 1194 state- to-state transitions and covers the kinetic temperature range from 30 to 3000 K. We note that the ori...

work page 2026
[5]

The formation tempera- tureT form was taken as 2/3×(E H2(v,J)−∆H o 0), as observed and recommended by Faure et al

andE CH+(υ′,J ′) andE H2(υ,J) are the ro-vibrational energies of CH + and H 2, respectively. The formation tempera- tureT form was taken as 2/3×(E H2(v,J)−∆H o 0), as observed and recommended by Faure et al. (2017). We note that in the mea- surements of the deuterated variants of the H 2 +H + 2 reaction, about two-thirds of the reaction exothermicity was ...

work page 2017
[6]

(2017), typically within a factor of

was found to be in good agree- ment with the more accurate TIQM calculations of Faure et al. (2017), typically within a factor of

work page 2017
[7]

The SQM thermal rate co- efficient thus matches well with the thermal TIQM value and also with various measurements, except below∼50 K where a strong decrease of the rate coefficient was observed in the ion- trap experiment of Plasil et al. (2011). This surprising result was interpreted by these authors as a dramatic loss of reactiv- Article number, page ...

work page 2011
[8]

and that the corresponding rate coefficients do not signif- 10-13 10-12 10-11 10-10 10-9 10-8 10-7 10-6 10-5 101 102 103 CH+(v=0, J=0) + e- → CH+(v', J' ) + e- Experiment (v'=0, J'=1) Theory (v'=0, J'=1) Theory (v'=0, J'=2) Theory (v'=0, J'=3) Theory (v'=1, J'=0) Theory (v'=1, J'=1) Theory (v'=1, J'=2) Theory (v'=1, J'=3) Rate coefﬁcient (cm3s-1) Kinetic ...

work page 2022
[9]

D.1.Emissivity profile of CH + pure rotational lines and ro- vibrational lines as a function of distance

CH+ ro-vibrational lines (v = 1-0) Fig. D.1.Emissivity profile of CH + pure rotational lines and ro- vibrational lines as a function of distance. 10-26 10-25 10-24 10-23 10-22 10-21 10-20 10-19 10-18 10-17 0.014 0.016 0.018 0.02 0.022 0.024 0.026 0.028 H2 lines emissivity (erg cm-3 s-1) Distance (pc) v = 0-0 v = 1-0 v = 2-1 Fig. D.2.Emissivity profile of ...

work page 2026
[10]

CH+ lines Wavelength Upper state Surface Brightness (erg cm −2 s−1 sr−1) (υ′,J ′)→(υ,J) (µm) energy/k B (K) Without pumping With pumping υ=1→0 (1,0)→(0,1) [P(1)] 3.6876 3942 3.14×10 −5 4.22×10 −5 (1,1)→(0,0) [R(0)] 3.6146 3980 1.18×10 −5 2.04×10 −5 (1,1)→(0,2) [P(2)] 3.7272 3980 4.36×10 −5 7.52×10 −5 (1,2)→(0,3) [P(3)] 3.7689 4058 6.27×10 −5 1.10×10 −4 (1...

work page 1952

[1] [1]

A., & Alsolami, Z

Ali, A., Ismail, H. A., & Alsolami, Z. 2015, Ap&SS, 357, 21 Aller, L. H. & Czyzak, S. J. 1983, ApJS, 51, 211 Amitay, Z., Zajfman, D., Forck, P., et al. 1996, Phys. Rev. A, 54, 4032 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, The Astro- physical Journal, 935, 167 Basart, J. P. & Daub, C. T. 1987, ApJ, 317, 412 Berné, O., Martin-Dru...

work page 2015

[2] [2]

Finally, for the pho- todissociation of CH + (reaction R6), theA V-dependent rate co- efficient was taken from Heays et al. (2017). For the reactions R1, R3, and R4, the rate coefficients were fitted using a standard modified-Arrhenius equation. We note that our fits are reliable over the kinetic temperature ranges 100−3000 K for reaction R1, 10−10 000 K ...

work page 2017

[3] [3]

(2017), based on the more accurate but expensive time- dependent wave-packet (TDWP) method, and a good agreement (within a factor of

and to the results of Faure et al. (2017), based on the more accurate but expensive time- dependent wave-packet (TDWP) method, and a good agreement (within a factor of

work page 2017

[4] [4]

1 of González-Lezana et al

We note that a good agreement was also observed between the SQM calculations and experiments performed with H 2(υ=0,1) at different temperatures, see e.g., Fig. 1 of González-Lezana et al. (2026). The full set of SQM rate coefficients includes 1194 state- to-state transitions and covers the kinetic temperature range from 30 to 3000 K. We note that the ori...

work page 2026

[5] [5]

The formation tempera- tureT form was taken as 2/3×(E H2(v,J)−∆H o 0), as observed and recommended by Faure et al

andE CH+(υ′,J ′) andE H2(υ,J) are the ro-vibrational energies of CH + and H 2, respectively. The formation tempera- tureT form was taken as 2/3×(E H2(v,J)−∆H o 0), as observed and recommended by Faure et al. (2017). We note that in the mea- surements of the deuterated variants of the H 2 +H + 2 reaction, about two-thirds of the reaction exothermicity was ...

work page 2017

[6] [6]

(2017), typically within a factor of

was found to be in good agree- ment with the more accurate TIQM calculations of Faure et al. (2017), typically within a factor of

work page 2017

[7] [7]

The SQM thermal rate co- efficient thus matches well with the thermal TIQM value and also with various measurements, except below∼50 K where a strong decrease of the rate coefficient was observed in the ion- trap experiment of Plasil et al. (2011). This surprising result was interpreted by these authors as a dramatic loss of reactiv- Article number, page ...

work page 2011

[8] [8]

and that the corresponding rate coefficients do not signif- 10-13 10-12 10-11 10-10 10-9 10-8 10-7 10-6 10-5 101 102 103 CH+(v=0, J=0) + e- → CH+(v', J' ) + e- Experiment (v'=0, J'=1) Theory (v'=0, J'=1) Theory (v'=0, J'=2) Theory (v'=0, J'=3) Theory (v'=1, J'=0) Theory (v'=1, J'=1) Theory (v'=1, J'=2) Theory (v'=1, J'=3) Rate coefﬁcient (cm3s-1) Kinetic ...

work page 2022

[9] [9]

D.1.Emissivity profile of CH + pure rotational lines and ro- vibrational lines as a function of distance

CH+ ro-vibrational lines (v = 1-0) Fig. D.1.Emissivity profile of CH + pure rotational lines and ro- vibrational lines as a function of distance. 10-26 10-25 10-24 10-23 10-22 10-21 10-20 10-19 10-18 10-17 0.014 0.016 0.018 0.02 0.022 0.024 0.026 0.028 H2 lines emissivity (erg cm-3 s-1) Distance (pc) v = 0-0 v = 1-0 v = 2-1 Fig. D.2.Emissivity profile of ...

work page 2026

[10] [10]

CH+ lines Wavelength Upper state Surface Brightness (erg cm −2 s−1 sr−1) (υ′,J ′)→(υ,J) (µm) energy/k B (K) Without pumping With pumping υ=1→0 (1,0)→(0,1) [P(1)] 3.6876 3942 3.14×10 −5 4.22×10 −5 (1,1)→(0,0) [R(0)] 3.6146 3980 1.18×10 −5 2.04×10 −5 (1,1)→(0,2) [P(2)] 3.7272 3980 4.36×10 −5 7.52×10 −5 (1,2)→(0,3) [P(3)] 3.7689 4058 6.27×10 −5 1.10×10 −4 (1...

work page 1952