{"id":"e48d38cb-e8b9-40d5-b776-b2d4400b7c71","arxiv_id":"1908.04058","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"At H2 columns above about 10^25 cm^-2, trapped infrared radiation heats inner dust and makes HCN vibrational lines bright and optically thick, explaining observations of buried galactic nuclei.","lead":"This paper models how thick dust cocoons around galactic nuclei trap infrared radiation and get much hotter inside, like a greenhouse. The models explain bright HCN vibrational lines seen in several buried nuclei and predict ring-like images and bright millimeter continuum.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Greenhouse amplification and the inferred HCN calibration are only as robust as the isotropic, smooth-column assumption that §4 itself flags; a 3D clumpy/disk test is needed.","rationale":"The reader and I identify the same soft spot, and the authors themselves flag it: §4 says spherical symmetry may overestimate Tdust. The central claim is not merely that backwarming exists—it is that the calibrated models reproduce observed HCN vibrational brightnesses with specific ΣIR and X_HCN. That calibration inherits the smooth-sphere assumption. I do not think this rises to rejection: the DUSTY benchmark, energy conservation, the plausible luminosities, and the predicted ring-like morphology in IC 860 are real supporting evidence. But it does justify the conditional verdict. Since the reader already reached CONDITIONAL and my concern is the same, no verdict change is needed; a 3D test would be the decisive next step.","tokens_in":34327,"tokens_out":6904,"duration_ms":88329,"concrete_test":"Use a 3D Monte Carlo radiative transfer code (e.g., RADMC-3D) to recompute the fiducial NGC 4418 model (ΣIR=2.2×10^8 L_sun/pc^2, NH2=10^25 cm^-2, q=1) with the same total dust mass and luminosity but distributed as clumps with volume filling factor f=0.1 and characteristic clump optical depth tau_20um≈30; compare the volume-averaged 14 μm mean intensity and the HCN nu2=1 f J=3−2 flux with the smooth spherical model. If either drops by more than a factor 2, the Fig. 12 calibration and the inferred ΣIR and X_HCN are geometry-limited; if the change is <30%, the spherical approximation is adequate for the claimed luminosities.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The argument's quantitative bridge is Fig. 2d: at NH2=10^25 cm^-2, tau_20um≈300 makes the volume-averaged 14 μm intensity ~10× larger, and §3.2 uses this field to populate HCN nu2=1 to optical thickness. That enhancement is computed for a smooth spherical rho∝r^-q density law with no clumping, so the high column covers every direction. The authors concede in §4 that this 'assumes isotropic column densities from the center and no clumpiness' and 'may overestimate Tdust' in real systems; clumps or a disk open low-column escape paths, reducing the mean mid-IR intensity and the nu2 population. Because the source comparison in Fig. 12 and Table 2 reads observed FHCN/ΔΩ as ΣIR (with X_HCN≈10^-6), even a factor 2 change in line flux—plausible for volume filling factors of a few tenths—shifts the inferred ΣIR by roughly a factor 4 given the saturated-area scaling of eq. (5). The greenhouse effect itself is physically plausible, and the DUSTY benchmark and IC 860 ring morphology are independent support; the weak point is specifically using a smooth, isotropic cocoon to convert observed brightness into physical parameters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents spherically symmetric radiative-transfer models of buried galactic nuclei (BGNs) with high H2 columns and high luminosity surface densities, for both AGN-like central heating and distributed starburst heating. The authors compute dust temperature profiles and emergent SEDs, and post-process them with an HCN excitation model including the nu2=1 bending state. They find that for NH2 about 10^25 cm^-2 or more, trapping of infrared radiation enhances inner dust temperatures and the mean mid-infrared intensity by more than an order of magnitude, pumping HCN vibrational states so that the nu2=1 J=3-2 and 4-3 lines become optically thick. They use the model grid to interpret observed HCN vibrational fluxes in NGC 4418, Arp 220W/E, Zw 049.057, IC 860, and Mrk 231, inferring Sigma_IR around (0.5-2)e8 Lsun pc^-2 with X_HCN about 1e-6 and luminosities consistent with independent estimates. The paper also makes predictions for line ratios, spatial profile shapes, and millimeter continuum brightness.","tokens_in":34536,"tokens_out":6948,"duration_ms":73182,"significance":"If correct, the paper provides a quantitative physical basis for the bright HCN vibrational emission in buried nuclei: the greenhouse effect of infrared trapping makes the nu2=1 lines a natural luminosity-surface-density diagnostic. The modeling is carefully done and the paper gives credit where due: the continuum code conserves energy to better than 1%, is benchmarked against DUSTY (Appendix A.3), and the temperature profiles are provided as analytic fits (Tables 3-4). The paper also gives falsifiable predictions, including saturated line ratios near 1.8, a central brightness drop or ring morphology as observed in IC 860, and absorption of millimeter continuum by the lines. The main risk is that the core quantitative result is derived for smooth, isotropic, spherical density distributions, an assumption the authors acknowledge may overestimate dust temperatures in real clumpy or disk systems; this affects the derived Sigma_IR calibration. Nonetheless, the paper is a significant advance in modeling BGNs and provides a framework that can be tested against higher-resolution observations.","major_comments":[{"comment":"The authors state in Section 4 that the spherical symmetry 'assumes isotropic column densities from the center and no clumpiness' and that this 'oversimplified smoothed density structure may overestimate the dust temperature as compared with real systems.' This caveat bears directly on the central claim: the factor of at least 10 enhancement of the 14 micron mean intensity at NH2=10^25 cm^-2 (Figure 2d) and the resulting optically thick HCN nu2=1 lines are computed for a smooth, isotropic cocoon. In a clumpy or disk-like medium, radiation escapes along low-column sightlines, reducing the mean mid-infrared intensity that pumps nu2=1. Because the calibration in Figure 12 and Table 2 maps observed F_HCN/Delta_Omega to Sigma_IR through exactly this mechanism (equation 5), the inferred Sigma_IR and luminosities are upper limits unless a filling-factor or three-dimensional geometry test is performed. I request a quantitative sensitivity study, for example clumpy or disk models with conservative filling factors, or at minimum an explicit statement that all derived Sigma_IR are upper limits under the smooth-sphere assumption.","section":"Section 4 (Discussion) and Section 3.2.2 / Figure 2d"},{"comment":"The source comparison uses fiducial values X_HCN=10^-6 and Delta_V=67 km/s that are themselves partly motivated by previous analyses of the same sources (for example NGC 4418 and Arp 220), and the observed HCN fluxes are then used to infer Sigma_IR. The agreement in Figure 12 is therefore a demonstration of consistency rather than an independent inversion. The degeneracy is substantial: with the saturated-area scaling of equation (5), a factor of about 2 uncertainty in X_HCN, which the authors assign in Section 4, translates into a factor of more than 2 change in the inferred Sigma_IR for the same line flux. The paper should provide a joint constraint plot, such as the allowed Sigma_IR-X_HCN locus per source, or otherwise quantify how the inferred physical parameters depend on the assumed abundance and velocity dispersion.","section":"Section 3.1, Section 3.2.5, Figure 12, Table 2"},{"comment":"The comparison with Arp 220W's 2.6 mm brightness relies on an extrapolation: for NH2 above 10^25 cm^-2 the authors 'simply assumed that the Tdust profile remains the same as for NH2=10^25 cm^-2' and argue that the inferred brightness temperatures are lower limits because Tdust increases with NH2. This monotonic increase is plausible but is not demonstrated at these columns, and the increasing optical depth at 2.6 mm could instead saturate the brightness. A self-consistent radiative-transfer calculation for NH2 around 10^26 cm^-2, or an explicit argument for why the profile is unchanged, is needed to support the favorability claim for an AGN in Arp 220W based on the high 2.6 mm brightness.","section":"Section 3.2.3 / Figure 11b"}],"minor_comments":[{"comment":"The phrase 'spectral enery distribution' is a typo and should read 'spectral energy distribution'.","section":"Section 2.1"},{"comment":"The text contains 'an squematic approach' and the Figure 15 caption contains 'opticallt thin'; these should be corrected to 'a schematic approach' and 'optically thin'.","section":"Appendix A.1 and Figure 15 caption"},{"comment":"The reference 'Dekel, & Burkert 2014' is incomplete, lacking a journal or preprint identifier, and the in-text citation 'Förster Schreiber et al. 2003' appears in the reference list with the year 1993; these should be reconciled.","section":"References"},{"comment":"Equation (5) is introduced with the remark that it is 'only valid for NHCN = 10^19 cm^-2'; the text should clarify whether this restriction also limits the applicability of the proportionality argument used in Section 3.2.4 for scaling line fluxes to other sources.","section":"Section 3.2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid modeling study with a clear central caveat that the authors themselves state. The revision should either add a geometry-robustness test or frame the derived Sigma_IR and HCN abundance as upper limits under the smooth-sphere assumption; this is feasible within the manuscript's scope. I do not see grounds for rejection, but the calibration claim needs to be hardened before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper does what a good modeling paper should do. It takes a known effect (backwarming in optically thick dusty environments), works out its consequences for buried galactic nuclei in enough detail to explain bright HCN nu2=1 lines, cold SEDs, and bright millimeter continuum with a single set of parameters, and then states the key caveats plainly. I think the broad picture is right, and the limitations are real but not fatal.\n\nWhat is genuinely new is the quantitative link. Rowan-Robinson discussed backwarming in 1982, Rolffs et al. applied it to hot cores, but this is the first systematic spherical radiative-transfer treatment aimed at extragalactic BGNs. It yields a simple calibration between HCN vibrational brightness and Sigma_IR, tested against six galaxies with luminosities that match independent far-IR absorption estimates. The continuum code is benchmarked against DUSTY, conserves energy to better than 1%, and the analytic fits to Tdust profiles make reimplementation feasible even without shipped code. The predicted ring-like morphology for saturated HCN lines is a nice, checkable outcome, and it matches the IC 860 maps.\n\nThe soft spot is exactly what the stress-test flags. The greenhouse boost that multiplies the 14 micron intensity by ~10 relies on columns of 1e25 cm^-2 covering essentially every line of sight. A smooth spherical r^-q density law ensures that; real nuclei with clumps or a disk will have low-column escape paths, reducing the mean mid-IR intensity and the nu2 population. The authors concede this in Section 4, saying the smooth structure may overestimate Tdust. Their counterargument (beam-averaged columns are high, radiation pressure may drive spheroidality) is plausible but not quantified. Because the source comparison reads observed HCN fluxes as Sigma_IR at fixed X_HCN, a factor 2 change in line flux can shift Sigma_IR by roughly a factor 4 through the saturated-area scaling of eq. (5). So the absolute values in Table 2 carry a systematic uncertainty larger than the listed ranges.\n\nThat said, the mechanism does not collapse. At 1e25 cm^-2 the trapping is strong enough that some enhancement survives moderate clumpiness, and the same calibration yields luminosities consistent with independent data across the whole sample. The 4-3/3-2 ratio mismatch in Arp 220W is minor. The authors are not hiding anything; they list the limitations themselves.\n\nWho is this for: observers and modelers working on obscured galactic nuclei, vibrational lines, or dusty radiative transfer. It deserves a serious referee. My recommendation: send it to peer review, with a request for a major revision where the clumping/geometry effect is either treated with a two-phase model or explicitly bounded.","headline":"A credible quantitative case for greenhouse-driven HCN vibrational emission in buried nuclei, with an honest and accurate statement of the spherical-symmetry limit that should be tested with 3D models.","tokens_in":35134,"tokens_out":2564,"would_cite":true,"duration_ms":31272,"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":"When hydrogen columns in buried galactic nuclei exceed about $10^{25}$ cm$^{-2}$, infrared photons are trapped so efficiently that the inner dust heats strongly and the HCN bending-mode lines become optically thick, explaining observed…","keywords":["buried galactic nuclei","HCN vibrational emission","dust radiative trapping","greenhouse effect","infrared luminous galaxies","active galactic nuclei","starburst","radiative transfer"],"falsifier":"A spatially resolved map of the HCN $\\nu_2=1$ $J=3-2$ line toward a nucleus with $N_{\\rm H_2}\\approx 10^{25}$ cm$^{-2}$ should show a central dip in line brightness (ring-like morphology), because the saturated line absorbs the bright 1.1 mm continuum; detecting a centrally peaked line instead would contradict the optically thick, greenhouse-driven picture.","tokens_in":34062,"feed_emoji":"🔥","tokens_out":6113,"duration_ms":63724,"temperature":0.7,"pith_summary":"This paper argues that the bright vibrational lines of HCN seen toward buried galactic nuclei are a direct consequence of a greenhouse effect in the dust: when the hydrogen column density exceeds roughly $10^{25}$ cm$^{-2}$, infrared photons emitted by dust cannot escape, so they heat the inner regions far above the temperature expected from the luminosity alone. The trapped light boosts the mid-infrared intensity by more than an order of magnitude, populating the $\\nu_2=1$ bending state of HCN and making the $J=3-2$ and $4-3$ vibrational lines optically thick. The authors show that a single set of parameters—HCN abundance near $10^{-6}$ relative to H$_2$ and surface brightness $\\Sigma_{\\rm IR}\\sim (0.5-2)\\times 10^8\\ L_\\odot\\ {\\rm pc}^{-2}$—reproduces the observed line fluxes in NGC 4418, Arp 220, IC 860, Zw 049.057, and Mrk 231, while also matching far-infrared photosphere temperatures and bright millimeter continuum. A sympathetic reading is that this is the first quantitative case that the same radiation trapping which heats the dust also powers the cyanopolyne vibrational emission, turning a qualitative idea into calibrated predictions.","feed_headline":"Dust cocoons trap infrared, heating nuclei and lighting up HCN","feed_subtitle":"A greenhouse effect in dense nuclear dust explains bright HCN vibrational lines and the galaxy luminosities behind them.","key_machinery":"The load-bearing object is the dust-temperature profile obtained from a spherically symmetric continuum radiative-transfer calculation in which radiation is carried by parallel rays through a power-law density cocoon ($\\rho\\propto r^{-q}$, $q=1$ or 1.5), heating shells by local absorption and re-emission. A compact blackbody at 1300 K represents the AGN case; a distributed energy deposition proportional to the dust mass and density represents the starburst case. The greenhouse effect appears in the equilibrated $T_{\\rm dust}(r)$: at high columns, inward (backwarming) fluxes almost cancel outward fluxes, so $\\Upsilon_{\\rm IR}=4\\pi r^2\\sigma T_{\\rm dust}(r)^4$ is not conserved and the inner shells are far hotter than the optically thin solution. These temperatures are then fed into a model of HCN with 25 rotational levels in the ground vibrational state and up to 48 levels in the $\\nu_2=1$ bending state, treating gas and dust as thermally coupled and including line–dust extinction and ro-vibrational overlaps. The single most important relation is that the HCN $\\nu_2=1\\ f\\ J=3-2$ line is in LTE at the local dust temperature wherever the 14 $\\mu$m continuum is optically thick, so the emergent line flux is set by the solid angle of the region where $T_{\\rm dust}\\approx 200$ K.","core_discovery":"The central claim is that radiative trapping, not a hotter central engine, is what makes buried galactic nuclei appear to have warm interiors. In the models, once $N_{\\rm H_2}\\gtrsim 10^{25}$ cm$^{-2}$ the optical depth at 20 $\\mu$m reaches hundreds, so any photon emitted in the inner shells is absorbed and re-emitted many times before it escapes. This backwarming raises the dust temperature in the inner third of the source to roughly 200–500 K, even when the externally observed spectral energy distribution looks cold, and raises the mean mid-infrared intensity inside the cocoon by more than a factor of ten. At those temperatures and columns, the HCN $\\nu_2=1$ state is populated so efficiently that the $J=3-2$ and $4-3$ vibrational lines saturate ($\\tau\\gtrsim 1$) over a substantial fraction of the source, with flux ratios near the optically thick value $(\\nu_{4-3}/\\nu_{3-2})^2\\approx 1.8$. The same greenhouse that traps the continuum therefore dictates the line luminosity, and the authors calibrate this to match the observed brightnesses in five galaxies with one fiducial abundance and a narrow range of surface brightness.","pith_inferences":["Editorial inference: If real buried nuclei are clumpy or disk-like, the greenhouse heating is likely weaker than the spherical models predict, so the inferred surface brightness values should be read as upper bounds; reproducing the observed lines in a clumpy medium would require even higher columns or luminosities.","Editorial inference: The same trapped-radiation mechanism should boost vibrational lines of other cyanopolynes, such as HC$_3$N $\\nu_7$ and $\\nu_6$ and HNC, which the authors list as future work; the ratio of HC$_3$N to HCN vibrational lines could serve as a cleaner thermometer of the inner cocoon.","Editorial inference: Because the photon-diffusion timescale ($\\sim 10^4$ yr) is comparable to AGN flickering timescales, a faded AGN could leave a fossil greenhouse cocoon that still emits HCN vibrational lines and shines at millimeter wavelengths, making a buried AGN resemble a starburst.","Editorial inference: The predicted steep rise of HCN vibrational luminosity with $\\Sigma_{\\rm IR}$ and its saturation at high columns can be tested by stacking unresolved galaxies that have measured compact dust masses; sources below a threshold of roughly $10^7\\ L_\\odot$ pc$^{-2}$ should show much weaker lines."],"forward_implications":["The HCN $\\nu_2=1$ $J=3-2$ and $4-3$ lines should be optically thick in buried galactic nuclei with $N_{\\rm H_2}\\gtrsim 10^{25}$ cm$^{-2}$ and $\\Sigma_{\\rm IR}\\gtrsim 10^7\\ L_\\odot$ pc$^{-2}$, with a flux ratio close to 1.8.","The same models predict bright, compact (sub)millimeter continuum with brightness temperatures of several hundred Kelvin, especially for AGN-heated cocoons, providing a way to spot buried active nuclei.","Far-infrared photosphere temperatures of 80–160 K emerge naturally from the greenhouse models, matching the temperatures inferred from high-lying molecular absorption lines in these galaxies.","A central dip or ring-like morphology of the HCN vibrational line is expected, since the optically thick line absorbs the bright millimeter continuum near the center.","For the diagnostics considered, AGN and starburst models give nearly identical HCN vibrational line fluxes, so the lines alone cannot distinguish the heating source; millimeter continuum peaks and central mass estimates are needed."],"supporting_citations":[{"why":"First extragalactic detection of the HCN $\\nu_2=1$ lines in NGC 4418, defining the observed phenomenon and supplying the line flux used for calibration.","marker":"Sakamoto et al. 2010"},{"why":"Provides HCN vibrational line fluxes for IC 860 and Zw 049.057, anchoring the comparison sample.","marker":"Aalto et al. 2015b"},{"why":"High-resolution 3 mm continuum observations of Arp 220W that give the central brightness and column density constraints used to favor AGN models.","marker":"Sakamoto et al. 2017"},{"why":"Far-infrared molecular absorption and compact source-size estimates for NGC 4418, used to set the source geometry, luminosity, and photospheric temperature.","marker":"González-Alfonso et al. 2012"},{"why":"Benchmark radiative-transfer calculation used to validate the continuum code's temperature profiles and emergent SEDs.","marker":"Ivezić & Elitzur 1997"},{"why":"Introduced the backwarming or greenhouse concept for dust in optically thick regions, the physical mechanism on which this paper builds.","marker":"Rowan-Robinson 1982"},{"why":"Analysis of the 14 $\\mu$m HCN band toward buried nuclei, providing an independent estimate of $X_{\\rm HCN}\\sim 10^{-6}$.","marker":"Lahuis et al. 2007"},{"why":"Provides the radiative-transfer method for molecular excitation with dust-mixed gas and line–dust extinction used in the HCN models.","marker":"González-Alfonso & Cernicharo 1999"}],"fun_headline_variants":["Greenhouse effect in galactic cores ignites HCN lines","Trapped IR heats buried nuclei, saturates HCN vibrational lines","Galactic greenhouses: dust traps IR, fires up HCN","Dust trapping turns on HCN glow in galaxy cores"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The models assume smooth spherical symmetry with isotropic column densities and no clumping; the authors acknowledge this may overestimate the dust temperature, because in real systems radiation can escape along low-column sightlines and weaken the greenhouse heating.","fun_headline_variants_meta":{"raw":{"variants":["Greenhouse effect in galactic cores ignites HCN lines","Trapped IR heats buried nuclei, saturates HCN vibrational lines","Galactic greenhouses: dust traps IR, fires up HCN","Dust trapping turns on HCN glow in galaxy cores"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1565,"prompt_tokens":1171,"completion_tokens":394,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":787,"completion_tokens_details":{"reasoning_tokens":322}},"tokens_in":787,"tokens_out":394,"duration_ms":4719,"temperature":1.0,"reasoning_tokens":322,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:53:14.757309+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spatially resolved map of the HCN $\\nu_2=1$ $J=3-2$ line toward a nucleus with $N_{\\rm H_2}\\approx 10^{25}$ cm$^{-2}$ should show a central dip in line brightness (ring-like morphology), because the saturated line absorbs the bright 1.1 mm continuum; detecting a centrally peaked line instead would contradict the optically thick, greenhouse-driven picture.","supporting_citations":[{"cited_title":"2017, ApJ, 849, 14","cited_arxiv_id":null,"evidence_quote":"High-resolution 3 mm continuum observations of Arp 220W that give the central brightness and column density constraints used to favor AGN models."},{"cited_title":"2012, A&A, 541, A4","cited_arxiv_id":null,"evidence_quote":"Far-infrared molecular absorption and compact source-size estimates for NGC 4418, used to set the source geometry, luminosity, and photospheric temperature."},{"cited_title":"Phillips, P","cited_arxiv_id":null,"evidence_quote":"Introduced the backwarming or greenhouse concept for dust in optically thick regions, the physical mechanism on which this paper builds."},{"cited_title":"1999, ApJ, 525, 845","cited_arxiv_id":null,"evidence_quote":"Provides the radiative-transfer method for molecular excitation with dust-mixed gas and line–dust extinction used in the HCN models."}],"review_version":1}