{"id":"e7174eb6-21cb-4066-b0d8-6d87d05ece50","arxiv_id":"2508.12625","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Combining 1D image profile fitting along the outflow axis with SED fitting constrains the viewing angle of Cepheus A to about 34 degrees, but does not improve constraints for the poorly resolved G35.20-0.74N.","lead":"The authors develop a fitting method that combines spectral energy distribution fitting with 1D brightness profiles extracted along the outflow axes of massive protostars, using the ZT18 radiative transfer model grid. Applied to Cepheus A, the method tightens the inferred viewing angle; applied to the more distant G35.20-0.74N, SOFIA's spatial resolution limits its usefulness.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inclination constraint may be an artifact of the fixed outflow position angle and the linearly interpolated background; rotating the strip PA and varying the background model would settle this.","rationale":"The reader's weakest assumption (background subtraction) is real and I partially agree, but I see the more load-bearing issue as the joint dependence of the theta_view measurement on the fixed strip position angle, the post hoc strip width, and the background model. The paper is transparent about limitations and the external agreement with kinematic estimates of theta_view (26-28 deg) provides independent support, so CONDITIONAL is appropriate rather than REJECT. The proposed test directly checks whether the headline constraint survives variation of the two most important user-defined parameters. If it survives, the paper's claim is strengthened; if it fails, the headline inclination is an artifact of the extraction geometry.","tokens_in":23254,"tokens_out":1300,"duration_ms":14084,"concrete_test":"Re-run the IMPRO and combined fitting for Cep A with the strip PA rotated from 40 to 60 degrees in 2-degree steps, and separately replace the linear background with a constant background estimated at the outer edge of the observed strip. If the resulting theta_view distribution shifts by more than the quoted +5.1/-4.9 degree uncertainty, or widens beyond about 10 degrees, the claimed inclination constraint is not robust to the two largest user-defined analysis choices.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that IMPRO+SED fitting tightens theta_view for Cep A rests on comparing the observed 1D brightness profile to synthetic profiles from a single-axis symmetric ZT18 model, with several user-chosen quantities held fixed: the strip PA (50 deg), the strip width (20 arcsec, chosen post hoc to give the narrowest theta distribution), and the background, which is estimated at the core radius and linearly interpolated across the strip (Section 2, Eq. 2). The paper explicitly notes that the best combined fit is poor (chi^2=6.58) and that the 19.7 um profile shape is not reproduced (Section 3.1.4). If the true PA of the mid-infrared emission differs by even a few degrees from 50 deg, or if the background varies nonlinearly across the strip, the fitted A_V and the synthetic profile shape will change, and the narrow theta_view distribution could shift or widen. The model also assumes the profile asymmetry arises only from viewing geometry of a symmetric bipolar cavity, while Cep A shows a precessing jet (Cunningham et al. 2009) and nearby sources such as HW3c, so the observed asymmetry may have an additional origin. This does not invalidate the method, but it means the headline numeric value (33.9 deg) is not yet shown to be robust to the analysis choices the method itself leaves free.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops Image Profile (IMPRO) fitting: a method that extracts 1D mid-infrared brightness profiles along the outflow axis of massive protostars and combines them with SED fitting against the Zhang & Tan (2018) radiative-transfer model grid. The pipeline is described in detail, including strip extraction, beam-convolved rebinning, foreground extinction, and a linear background prescription. The method is applied to two SOFIA/FORCAST SOMA sources: Cepheus A, where the authors claim the viewing angle becomes tightly constrained (theta_view = 33.9 +5.1/-4.9 deg) compared to the almost unconstrained SED-only range, and G35.20-0.74N, where limited spatial resolution prevents significant improvement. The central claim is that MIR brightness profiles can break SED degeneracies in source inclination for well-resolved massive protostars, with JWST expected to extend the method to more distant sources.","tokens_in":23522,"tokens_out":4865,"duration_ms":55290,"significance":"If the central claim holds, the paper makes a useful methodological contribution: it provides a transparent, quantitative way to incorporate spatial brightness information into ZT18-grid fitting and demonstrates on Cep A that the inclination degeneracy can be broken. The pipeline is described with enough detail to be reproduced, and the inclination result is cross-checked against independent disk/outflow inclination measurements. The honest treatment of the G35.2N resolution limit is a strength, as is the explicit acknowledgement that the ZT18 grid does not simultaneously reproduce the SED and profile data for Cep A. However, the quantitative reliability of the headline theta_view constraint is not yet demonstrated, because the fit quality is poor and several user-chosen analysis choices (strip position angle, background model, good-model threshold) are not varied or tested.","major_comments":[{"comment":"The quoted inclination constraint theta_view = 33.9 +5.1/-4.9 deg is the mean and full range of 'good' models selected by chi2_combined <= chi2_min,combined + 2. Since the minimum reduced combined chi-square is 6.58 (stated in §3.1.3 and Figure 8), the absolute fit is unacceptable, and the +2 threshold is not a statistically meaningful confidence interval. The authors should either base uncertainties on a proper Delta-chi2 criterion or add a model-discrepancy term; as written, the 'tight constraint' claim is only a spread over grid models that are all systematically inconsistent with the data.","section":"§3.1.4, Table 1"},{"comment":"The background intensity I_bg is set to the observed profile value at the model core radius and linearly interpolated across the strip, with no uncertainty and no alternative background model. Because Eq. (2) couples I_bg and the foreground extinction A_V through f_AV, any error in the background changes the fitted A_V and the shape of the synthetic profile. For Cep A, the clump environment and the nearby HW3c source make the linear-background assumption particularly risky. A robustness test—for example, a polynomial background, a background measured from a parallel strip, or an additional background amplitude parameter—is needed to show that theta_view and A_V are stable.","section":"§2, Eq. (2)"},{"comment":"The strip position angle is fixed at 50 deg based on the CO outflow orientation, and only the strip width and the band combination are tested. The MIR morphology traces the warm outflow cavity, which may not be exactly aligned with the CO axis given the precessing jet in Cep A (Cunningham et al. 2009), and the synthetic profiles assume the strip is aligned with the model's symmetry axis. The authors should test the sensitivity of the theta_view distribution to the strip PA (e.g., 40-60 deg) and report whether the 33.9 deg result persists; without such a test, the headline constraint may be an artifact of the assumed strip orientation.","section":"§3.1.2, Fig. 1"},{"comment":"The combined fit has minimum reduced chi-square 6.58 and explicitly does not reproduce the 19.7 um profile shape; the best SED model and best IMPRO model imply core radii that differ by a factor of six (R_c = 0.3 pc versus 0.05 pc). This indicates that the ZT18 grid cannot simultaneously describe the SED and the MIR morphology under the adopted assumptions. The abstract and conclusions should state this limitation prominently, and the narrower good-model range should not be presented as a measurement uncertainty without adding a model-discrepancy component.","section":"§3.1.3-3.1.4, Fig. 9"}],"minor_comments":[{"comment":"Cunningham et al. (2009) appears twice in the reference list with identical entries; one duplicate should be removed.","section":"References"},{"comment":"'planer motion' should be 'planar motion'.","section":"§3.1.1"},{"comment":"The phrase 'stepsA_V' is missing a space; it should read 'steps A_V'.","section":"§2"},{"comment":"The parenthetical '20 ′′(see the discussion' is missing a space before the parenthesis.","section":"§2"},{"comment":"The notes state that uncertainties are the full range of the good-model distributions; it would be helpful to add a sentence clarifying that these are not confidence intervals, since the underlying model grid is sparsely sampled.","section":"Tables 1 and 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is closely tied to the authors' own ZT18 grid and SEDcreator package, and the novelty relative to Fedriani et al. (2023a) and Zhang et al. (2013b) should be made more explicit. The independent checks by Patel et al. (2005) and Sanna et al. (2017) are reassuring, but the poor absolute chi-square and untested analysis choices currently prevent the headline theta_view value from being treated as a robust measurement. With the requested robustness tests and more cautious wording, the paper would be a solid methodological contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the headline result, theta_view = 33.9 +5.1/-4.9 for Cep A instead of 29-90 degrees from SED-only fitting, is plausible and lines up with independent kinematics (26-28 degrees from maser/disk measurements). Second thing: the authors are upfront that the combined best fit is bad (reduced chi^2 = 6.58) and that the 19.7 micron profile shape is not reproduced. So read the narrow uncertainty as the spread of good models under their selection rule, not as a statement that the model family actually describes the data simultaneously.\n\nWhat is new is the systematic multi-band IMPRO pipeline: extracting 1D profiles from three FORCAST bands, marginalizing over foreground extinction, combining the profile chi^2 with the SED chi^2, and testing the whole thing on two sources. Zhang et al. (2013b) fitted one profile for G35.2N with an earlier model, so the novelty is incremental rather than groundbreaking, but this is a proper pipeline paper with a useful negative result. The G35.2N case honestly shows the method fails when the source is barely resolved, which is exactly the kind of calibration the community needs before JWST applications.\n\nThe soft spots are real but not fatal. First, the strip PA is fixed at 50 degrees and never rotated; Cep A's jet precesses and HW3c sits nearby, so the asymmetric profile that drives the inclination constraint could be affected by a few degrees of misalignment. Second, the background is linearly interpolated from the value at the core radius; if the clump emission is not linear, the fitted A_V and the synthetic profile comparison will shift. Third, the 20 arcsec strip width was chosen because it gives the narrowest theta distribution; the authors acknowledge it is user-defined, but that is still post hoc. Fourth, the SED and IMPRO constraints are not actually being satisfied together: the combined best model sits close to the IMPRO-only best model and under-predicts the 20-40 micron SED. Finally, no code or model grid is released, which makes adoption harder.\n\nI agree with the stress-test note: the central inclination result is externally corroborated, so I would not call it an artifact, but the paper has not shown it is robust to the analysis choices the method itself leaves free. Rotating the PA, varying the background model, and releasing the code would settle that.\n\nWho is this for? People doing SED fitting of massive protostars with grid models, and anyone planning MIR imaging of such sources with JWST. The method section is the contribution; the Cep A number is a test case. A serious referee should engage with this. I would send it out and ask for a revision that stress-tests the PA and background assumptions.","headline":"A clearly described SED+image-profile fitting pipeline that tightens the inclination constraint for Cepheus A, but the quoted precision comes with a poor joint chi^2 and several user-chosen analysis settings that are not fully stress-tested.","tokens_in":24101,"tokens_out":2400,"would_cite":true,"duration_ms":27620,"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":"A 1D mid-infrared strip along the outflow axis breaks the inclination degeneracy of a massive protostar, pinning Cepheus A at 34 degrees.","keywords":["massive protostars","SED fitting","mid-infrared image profiles","outflow cavities","viewing angle","radiative transfer models","SOFIA/FORCAST","Cepheus A"],"falsifier":"Re-run the Cepheus A fit with the background estimated from the median intensity just outside the model core radius instead of the linear interpolation; if the best-fitting $\\theta_{\\rm view}$ moves by more than the quoted $+5.1/-4.9$ degrees, the background assumption is the weak link. A complementary check is a resolved mid-infrared image showing outflow cavity edges that no ZT18 model with $\\theta_{\\rm view}\\approx 34^\\circ$ can reproduce, which would falsify the central claim directly.","tokens_in":23030,"feed_emoji":"🌌","tokens_out":9961,"duration_ms":98748,"temperature":0.7,"pith_summary":"This paper develops IMPRO, a fitting method that adds one-dimensional mid-infrared brightness profiles along a protostar's outflow axis to standard spectral energy distribution (SED) fitting. The aim is to break a known degeneracy in which the SED of a dense envelope seen through its outflow cavities can mimic a less dense envelope viewed closer to the midplane. The paper argues that the brightness asymmetry between the near-facing and far-facing outflow lobes carries independent information about the viewing angle. Applied to Cepheus A with SOFIA/FORCAST images at 19.7, 31.5, and 37.0 µm, the combined fit narrows the viewing angle from a nearly unconstrained 29–90 degrees to $33.9^{+5.1}_{-4.9}$ degrees, consistent with independent disk-based estimates. For the more distant G35.20-0.74N, the same pipeline does not improve on SED fitting, which the paper attributes to the three-times coarser linear resolution of the observations.","feed_headline":"Outflow-axis profiles pin a massive protostar's tilt","feed_subtitle":"Adding 1D mid-infrared brightness profiles to SED fits narrows Cepheus A's viewing angle to 34°.","key_machinery":"The carrying object is the one-dimensional brightness profile $I_{\\rm obs}$: the average intensity along a strip of width 20 arcseconds centered on the protostar and aligned with the outflow axis, rebinned to the observing beam. The same strip extraction is applied to synthetic images from the ZT18 grid, a grid of radiative transfer models of massive protostars built on the Turbulent Core Accretion scenario; the grid spans core mass, clump surface density, protostellar mass, and viewing angle sampled in $\\cos\\theta_{\\rm view}$, giving 8640 models before foreground extinction. Each synthetic profile is extinguished by a foreground $A_V$ using the grid's adopted extinction law, added to a linearly interpolated background determined at the model core radius, and compared with the observed profile through a reduced chi-square; that profile chi-square is averaged with the SED chi-square to rank models. The diagnostic power comes from the brightness asymmetry between the outflow lobes, which is strongest at 19.7 µm because extinction separates the near and far cavities.","core_discovery":"The central claim is that the radial brightness profile along the outflow axis, extracted from a 20-arcsecond strip in multi-band mid-infrared images and compared with the ZT18 turbulent-core radiative transfer grid, constrains the viewing angle $\\theta_{\\rm view}$ to about five degrees once the SED has fixed the broad physical parameters. For Cepheus A, SED-only fitting leaves $\\theta_{\\rm view}$ ranging from 29 to 90 degrees, while the combined SED+IMPRO fit gives $33.9^{+5.1}_{-4.9}$ degrees, agreeing with the 26–28 degree values inferred from disk and maser kinematics. The paper also shows that IMPRO-only and SED-only fits prefer different envelope parameters, and that the joint minimum reduced chi-square of 6.58 is worse than either method alone, which it reads as evidence that the model grid does not fully capture the sources. For G35.20-0.74N the profile fits do not tighten the parameter distributions, and the paper identifies the physical resolution of roughly 7700 au, against model peak offsets mostly below 2000 au, as the limiting factor.","pith_inferences":["As an extension, the same strip statistic applied to JWST NIRCam and MIRI images, with sub-arcsecond beams, would resolve about 70 au at Cepheus A's distance and directly test whether the 20-arcsecond strip washes out cavity-opening-angle information that the current pipeline cannot use.","If few-degree inclinations become routine, a sample-level comparison of outflow $\\theta_{\\rm view}$ with disk-plane orientations from ALMA could test whether massive protostellar outflows are perpendicular to their disks, a prediction of monolithic turbulent-core collapse that competitive accretion does not require.","The single-source, symmetric-outflow assumption is the first thing to relax; repeating the strip-width tests on synthetic images with a second embedded source would predict a systematic drift in fitted $\\theta_{\\rm view}$ with strip width, which could be checked against the multiple cores in G35.20-0.74N.","A direct test of the background prescription would be to fit synthetic images with known clump gradients; if the linear-background model biases $\\theta_{\\rm view}$ by more than the reported five-degree uncertainty, the quoted Cepheus A uncertainty is optimistic."],"forward_implications":["For nearby massive protostars whose outflow cavities are resolved, the viewing angle becomes a measured quantity with few-degree uncertainty, removing the largest source of degeneracy in envelope-model fitting.","A model must now reproduce both the SED and the spatial brightness distribution; the two can disagree, as in Cepheus A where the SED-only best model peaks at the source position while the observed profile peaks at an offset, so the combined fit exposes models that accidentally match one but not the other.","The method's usefulness is distance-limited: model peak offsets mostly sit below 2000 au, so sources beyond roughly 2 kpc require sub-arcsecond resolution before the profile adds information.","If the Cepheus A result is correct, the preferred model is a more compact, denser core with a lower protostellar mass (about 8 $M_\\odot$) than SED-only fitting suggested, but the elevated combined chi-square leaves room for missing physics such as cavity shape."],"supporting_citations":[{"why":"Supplies the five-parameter radiative transfer model grid whose synthetic images and SEDs are fitted by the pipeline.","marker":"Zhang & Tan 2018"},{"why":"Provides the revised SED fitting pipeline that treats foreground extinction as an independent parameter, used for the SED chi-square.","marker":"Fedriani et al. 2023a"},{"why":"Provides the SOFIA/FORCAST observations of Cepheus A and G35.20-0.74N and the 10 percent uncertainty convention adopted in the fits.","marker":"De Buizer et al. 2017"},{"why":"Supplies the extinction law applied to the synthetic profiles through foreground $A_V$ in Equation 2.","marker":"Kim et al. 1994"},{"why":"Presents an earlier brightness-profile fit along the G35.2N outflow with a similar model framework, used as a benchmark for the new pipeline.","marker":"Zhang et al. 2013b"},{"why":"Gives an independent Cepheus A disk inclination whose implied outflow inclination is compared with the fitted $\\theta_{\\rm view}$.","marker":"Patel et al. 2005"},{"why":"Gives methanol maser kinematics implying an outflow inclination of about 26 degrees, another independent comparison for the Cepheus A result.","marker":"Sanna et al. 2017"},{"why":"Provides the updated SED fitting treatment for multiple and crowded sources, including the Cepheus A luminosity values cited in the analysis.","marker":"Telkamp et al. 2025"}],"fun_headline_variants":["IMPRO fit narrows massive protostar's tilt to 34°","Image-profile fitting tightens protostar inclination","SED plus image profiles fixes protostar tilt angle","IMPRO narrows Cepheus A's viewing angle to 34°","One-dimensional profiles pin down protostar orientation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the observed brightness outside the model core radius is a smooth, linearly interpolated background and that a single protostar with a symmetric bipolar outflow dominates the strip; if the surrounding clump is clumpy, the strip contains a second source, or the outflow cavities are asymmetric, the fitted $A_V$ and $\\theta_{\\rm view}$ will be biased.","fun_headline_variants_meta":{"raw":{"variants":["IMPRO fit narrows massive protostar's tilt to 34°","Image-profile fitting tightens protostar inclination","SED plus image profiles fixes protostar tilt angle","IMPRO narrows Cepheus A's viewing angle to 34°","One-dimensional profiles pin down protostar orientation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000662,"raw_usage":{"total_tokens":3116,"prompt_tokens":1129,"completion_tokens":1987,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":745,"completion_tokens_details":{"reasoning_tokens":1901}},"tokens_in":745,"tokens_out":1987,"duration_ms":15452,"temperature":1.0,"reasoning_tokens":1901,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:20:59.858516+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the Cepheus A fit with the background estimated from the median intensity just outside the model core radius instead of the linear interpolation; if the best-fitting $\\theta_{\\rm view}$ moves by more than the quoted $+5.1/-4.9$ degrees, the background assumption is the weak link. A complementary check is a resolved mid-infrared image showing outflow cavity edges that no ZT18 model with $\\theta_{\\rm view}\\approx 34^\\circ$ can reproduce, which would falsify the central claim directly.","supporting_citations":[],"review_version":1}