pith. sign in

arxiv: 2604.09006 · v1 · submitted 2026-04-10 · 🌌 astro-ph.GA

Photoionization modelling of circumstellar nebulae using irregular grains

P. Jim\'enez-Hern\'andez , S.J. Arthur , D. Guirado , O. Mu\~noz , J. Martikainen , L. Sabin , W.J. Henney This is my paper

Pith reviewed 2026-05-10 18:14 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords photoionization modellingcircumstellar nebulaeirregular dust grainshexahedral grainsspherical grainsinfrared luminositydust massMRN size distribution
0
0 comments X

The pith

Irregular hexahedral dust grains produce up to 60 percent higher infrared luminosities than spherical grains in photoionization models of circumstellar nebulae.

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

The paper tests the impact of grain shape on models of ionized gas around evolved stars by replacing standard spherical dust particles with irregular hexahedral ones. Optical properties for the irregular grains come from a scattering database and are fed into a spectral synthesis code alongside the usual spherical properties. Both sets of grains follow the same power-law size distribution from 0.005 to 0.25 micrometers. The resulting model spectra diverge most at longer wavelengths and for larger grain sizes, with the infrared peak luminosity rising by as much as 60 percent when hexahedral grains are used. This shift implies that earlier calculations based on spherical grains may have required more total dust to match observed fluxes.

Core claim

When the same MRN size distribution is applied to both populations, the nebular continuum calculated with irregular hexahedral grains shows infrared peak luminosities up to 60 percent higher than the equivalent spherical-grain models, with the discrepancy largest for graphite and for the biggest grains. The authors conclude that the spherical-grain assumption therefore tends to overestimate the dust mass required to reproduce observed emission from photoionized circumstellar nebulae.

What carries the argument

Opacities derived from the scattering properties of irregular hexahedral grains in the TAMUdust2020 database, inserted into the cloudy code and compared directly to spherical-grain opacities over an MRN size distribution for graphite, amorphous carbon, and silicate.

If this is right

  • The difference between spherical and hexahedral models grows with grain size and is strongest for graphite.
  • Spherical-grain models systematically overestimate the dust mass needed to produce a given infrared luminosity.
  • The choice of grain shape affects the derived nebular continuum most at wavelengths near the infrared peak.
  • Amorphous carbon and silicate grains show smaller but still measurable differences compared with graphite.

Where Pith is reading between the lines

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

  • Dust-mass estimates derived from infrared observations of other dusty environments, such as supernova remnants or star-forming regions, may need similar shape corrections.
  • Radiative-transfer codes used for circumstellar disks or AGN tori could be updated with irregular-grain libraries to test consistency with nebula results.
  • Laboratory measurements of scattering from other non-spherical shapes could be added to the same database to map how sensitive the luminosity offset is to exact geometry.

Load-bearing premise

The optical properties measured for laboratory hexahedral grains accurately describe the actual dust particles present in circumstellar nebulae.

What would settle it

A side-by-side comparison of the predicted infrared spectral energy distributions from spherical versus hexahedral models against high-resolution infrared observations of a planetary nebula or similar circumstellar shell, checking whether the higher-luminosity hexahedral predictions reduce the required dust mass to match the data.

Figures

Figures reproduced from arXiv: 2604.09006 by D. Guirado, J. Martikainen, L. Sabin, O. Mu\~noz, P. Jim\'enez-Hern\'andez, S.J. Arthur, W.J. Henney.

Figure 1
Figure 1. Figure 1: Opacity parameters for silicate grains computed for two grain shapes: hexahedral (solid lines) and spherical (dotted lines), across the 10 bins of our size distribution (see [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Spectra from cloudy photoionization models with silicate grains, comparing spherical (blue line) and hexahedral-shaped (orange line) grains. The first two panels show models using only the grain sizes corresponding to bins 1 and 10 from [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Opacity parameters for graphite grains computed for two grain shapes: hexahedral (solid lines) and spherical (dotted lines), across the 10 bins of our size distribution (see [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Spectra from cloudy photoionization models with graphite grains, comparing spherical (blue line) and hexahedral-shaped (orange line) grains. The first two panels show models using only the grain sizes corresponding to bins 1 and 10 from [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Opacity parameters for AC amorphous carbon for two grain shapes: hexahedral (solid lines) and spherical (dotted lines), across the 10 bins of our size distribution (see [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Spectra from cloudy photoionization models with AC amorphous carbon grains, comparing spherical (blue line) and hexahedral-shaped (orange line) grains. The first two panels show models using only the grain sizes corresponding to bins 1 and 10 from [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Opacity parameters for BE amorphous carbon for two grain shapes: hexahedral (solid lines) and spherical (dotted lines), across the 10 bins of our size distribution (see [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Spectra from cloudy photoionization models with BE amorphous carbon grains, comparing spherical (blue line) and hexahedral-shaped (orange line) grains. The first two panels show models using only the grain sizes corresponding to bins 1 and 10 from [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The extinction efficiency (Qext), single-scattering albedo (𝜔), and asymmetry factor (𝑔) obtained from the TAMUdust2020 database, together with the imaginary part of the refractive index for each dust species. The opaque region of the lines indicates the values used to calculate the opacities. The different colours refer to the bin size used for each calculation. enhance far-IR emission without requiring l… view at source ↗
Figure 10
Figure 10. Figure 10: Dust grain temperature profiles from our photoionization models, computed using the MRN size distribution divided into 10 size bins for silicate, graphite and amorphous carbon grains. Solid lines represent models with hexahedral grains, while dotted lines correspond to spherical grain models. Future improvements will require more comprehensive measure￾ments of dust optical properties across the full range… view at source ↗
read the original abstract

We study the effects of using the optical properties of irregular hexahedral grains in photoionization models of circumstellar nebulae around evolved stars. Dust opacities for the irregular grains were obtained from the scattering properties available in the TAMUdust2020 database and these were implemented in the spectral synthesis code cloudy. A sample of photoionization models that use opacities from both spherical and irregular hexahedral grains across a standard MRN size distribution (0.005 to 0.25 um) was produced. We consider the optical properties of graphite, amorphous carbon and silicate dust grains and find that differences between the model nebula continua calculated using spherical and irregular dust grains increase with the grain size, especially for graphite. In particular, we find that the luminosities at the infrared peak for the hexahedral grain models can be up to 60% higher than those from the equivalent spherical grain models for the largest grains. This result suggests that traditional spherical grain assumptions may lead to an overestimate of the dust mass in photoionized nebulae.

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

3 major / 2 minor

Summary. The manuscript implements optical properties of irregular hexahedral grains from the TAMUdust2020 database into the Cloudy photoionization code for circumstellar nebulae around evolved stars. Using the MRN size distribution (0.005–0.25 μm) for graphite, amorphous carbon, and silicate grains, it compares the resulting model continua to those from equivalent spherical grain models and reports that IR peak luminosities can be up to 60% higher for hexahedral grains (especially larger graphite grains). This leads to the suggestion that spherical grain assumptions may cause dust mass overestimates in photoionized nebulae.

Significance. If the models are normalized to identical total dust mass, the reported shape-dependent opacity effect would be a useful benchmark for dust modeling in nebulae, highlighting how non-spherical grains can alter IR emission predictions in standard codes like Cloudy. The direct use of an external scattering database provides a concrete, reproducible implementation that could be adopted by others. However, the quantitative implication for mass estimates depends on validation of the mass equivalence and applicability to real circumstellar environments.

major comments (3)
  1. [Abstract and §2 (model setup)] Abstract and model description: The central claim that hexahedral models produce up to 60% higher IR luminosities (and thus that spherical models overestimate dust mass) requires that both grain populations are compared at identical total dust mass. The manuscript does not specify how the effective radius or size parameter from TAMUdust2020 hexahedra maps to grain volume/mass relative to the spherical case under the same MRN limits. If the 'size' definitions yield different masses, the luminosity difference partly reflects mass normalization rather than pure shape effects, directly undermining the mass-overestimate conclusion.
  2. [§2 (model setup) and §4 (discussion)] Model setup and discussion: The MRN size distribution is applied equally to both grain shapes without justification or sensitivity tests for its suitability to circumstellar ejecta (as opposed to ISM). Circumstellar dust formation may produce different size distributions, so the reported 60% difference, while internally consistent for the chosen distribution, does not securely support conclusions about real nebulae without additional checks or references to observed circumstellar grain properties.
  3. [§3 (results) and §5 (conclusions)] Results and conclusions: No direct comparison of the modeled continua to observed spectra of circumstellar nebulae is presented to test whether hexahedral models improve fits or alter inferred parameters. Without such validation, the suggestion that spherical assumptions lead to mass overestimates remains a modeling difference rather than a demonstrated improvement for observational applications.
minor comments (2)
  1. [Abstract] The abstract states that 'a sample of photoionization models' was produced but provides no details on the number of models, the range of stellar parameters, or nebular densities explored. Adding a brief summary or table of the parameter space would clarify the generality of the 60% figure.
  2. [§2 (model setup)] Notation for grain size (e.g., how 'largest grains' are defined quantitatively) and the exact definition of 'infrared peak luminosity' (e.g., wavelength range or peak value) could be stated more explicitly in the methods to aid reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive comments. We respond to each major comment below, indicating revisions where appropriate.

read point-by-point responses

  1. Referee: Abstract and §2 (model setup): The central claim that hexahedral models produce up to 60% higher IR luminosities requires identical total dust mass. The manuscript does not specify how the effective radius or size parameter from TAMUdust2020 hexahedra maps to grain volume/mass relative to the spherical case under the same MRN limits.

    Authors: We agree clarification is needed. The TAMUdust2020 size parameter for hexahedra is the equivalent-volume radius, which we matched to the spherical case when applying the MRN limits. This ensures equal grain volumes and masses (assuming identical densities) across corresponding sizes. We will revise §2 to explicitly state this volume equivalence and confirm the mass normalization. revision: yes

  2. Referee: §2 (model setup) and §4 (discussion): The MRN size distribution is applied equally without justification or sensitivity tests for its suitability to circumstellar ejecta (as opposed to ISM).

    Authors: MRN is a standard assumption in nebular modeling when specific circumstellar distributions are unknown, and it has been employed in prior studies of evolved-star nebulae. We will add supporting references in §2 and a note in §4 that the 60% difference is for this distribution, with exploration of alternatives reserved for future work. revision: partial

  3. Referee: §3 (results) and §5 (conclusions): No direct comparison of the modeled continua to observed spectra of circumstellar nebulae is presented to test whether hexahedral models improve fits or alter inferred parameters.

    Authors: The manuscript quantifies shape-induced differences in controlled photoionization models rather than claiming improved observational fits. The mass-overestimate implication follows directly from the modeled luminosity contrast at fixed mass. We will revise §5 to clarify the scope and identify observational validation as future work. revision: partial

Circularity Check

0 steps flagged

No circularity: direct model comparison using external opacities

full rationale

The paper implements scattering properties from the external TAMUdust2020 database into the Cloudy code, applies the standard MRN size distribution uniformly to both spherical and hexahedral grains, and computes the resulting nebular continua. The reported up-to-60% IR luminosity difference is a direct numerical output of these models rather than any fitted parameter, self-defined quantity, or self-citation chain. No load-bearing step reduces to its own inputs by construction; the comparison is self-contained against the chosen external data and code.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim depends on the accuracy of the TAMUdust2020 optical properties for hexahedral grains and on the assumption that the MRN size distribution can be applied unchanged to irregular grains.

free parameters (1)
  • MRN size distribution limits
    The 0.005 to 0.25 um range is adopted as standard without new fitting, but the choice still constitutes a modeling parameter.
axioms (2)
  • domain assumption Optical properties from TAMUdust2020 accurately represent real irregular grains in nebular conditions
    The database values are inserted directly into cloudy without additional validation against nebular observations.
  • domain assumption MRN size distribution applies equally to spherical and hexahedral grains
    The same distribution is used for both grain shapes without adjustment for shape-dependent coagulation or destruction.

pith-pipeline@v0.9.0 · 5507 in / 1455 out tokens · 35170 ms · 2026-05-10T18:14:07.339986+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

3 extracted references · 3 canonical work pages · 1 internal anchor

  1. [1]

    F., Huffman D

    Astropy Collaboration et al., 2013, A&A, 558, A33 Astropy Collaboration et al., 2018, AJ, 156, 123 Biermann P., Harwit M., 1980, ApJ, 241, L105 Bohren C. F., Huffman D. R., 1998, Absorption and Scattering of Light by Small Particles. Wiley-VCH Bussoletti E., Colangeli L., Borghesi A., Orofino V., 1987, A&AS, 70, 257 Chatzikos M., et al., 2023, Rev. Mex. A...

    work page internal anchor Pith review doi:10.48550/arxiv.astro-ph/0107183 2013

  2. [2]

    The shaded grey region indicates wavelengths where opacities were derived from TAMUdust2020 output data

    The contribution of the stellar spectrum is shown in black, and green diamonds mark the IR continuum peak of each model. The shaded grey region indicates wavelengths where opacities were derived from TAMUdust2020 output data. MNRAS000, 1–14 (2026) 16Jiménez-Hernández et al. 101 102 103 [ m] 1033 1034 1035 1036 1037 1038 1039 L [erg s 1] Bin 2 101 102 103 ...

    work page 2026

  3. [3]

    The shaded grey region indicates wavelengths where opacities were derived from TAMUdust2020 output data

    The contribution of the stellar spectrum is shown in black, and green diamonds mark the IR continuum peak of each model. The shaded grey region indicates wavelengths where opacities were derived from TAMUdust2020 output data. MNRAS000, 1–14 (2026)

    work page 2026