{"id":"06da13b9-199f-4a71-9d7f-348397a56c6a","arxiv_id":"2505.01264","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A lumped-parameter cardiovascular model suggests that after lung resection, both increased pulmonary afterload and reduced right ventricular contractility can produce similar volume changes but opposite changes in right-sided pressures.","lead":"The authors built a computer model of the heart and blood vessels to simulate what happens after lung removal surgery, and compared two suspected causes of right-heart dysfunction: increased blood pressure in the lungs and weakened right ventricle contraction. The model shows the two causes produce opposite changes in specific pressures, which could help doctors tell which mechanism is at play.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted PASP/RVSP increase that distinguishes afterload from contractility loss depends on the unvalidated assumption in Eq. (21) that proximal impedance scales as 1/(1−α). A γ-sweep of the scaling exponent would show whether the differentiator survives.","rationale":"The reader's weakest_assumption is well placed: the afterload-scaling block, and Eq. (21) in particular, is the quantitative hinge for the paper's headline pressure signal. I considered two other candidate concerns. First, the six Windkessel parameters are fitted by exhaustive search to reference ranges without a formal identifiability criterion; this affects baseline fidelity but not necessarily the directional sensitivities or the sign of the afterload-versus-contractility pressure difference. Second, the paper concludes that post-operative RV dysfunction 'may result from a combination' of the two mechanisms without simulating the combined case; this is a real inferential gap, but the authors explicitly hedge with 'may' and state that the relative extent of the two mechanisms is unknown, so the central claim is a hypothesis rather than a demonstrated quantitative result. By contrast, Eq. (21) is a modeling choice that directly sets the magnitude of the afterload-induced PASP increase, the very signal proposed as the clinical differentiator. The gamma-sweep test isolates this dependency and would show whether the differentiator is robust or an artifact of the parallel-segment assumption. Because the authors already flag the need for validation and the conclusions are appropriately cautious, a conditional acceptance remains the right verdict.","tokens_in":17702,"tokens_out":7797,"duration_ms":88584,"concrete_test":"Re-run the afterload scenario with Z_a,post = Z_a,pre × (1−α)^(−γ) for γ ∈ {0, 0.25, 0.5, 0.75, 1}, keeping Eqs. (20) and (22) for R and C unchanged, and recompute the α = 50% changes in RVSP, PASP, PADP, and RVEF. If for γ ≤ 0.5 the afterload-induced PASP increase drops below the contractility-loss decrease, or becomes comparable to it, the proposed pressure differentiator is an artifact of the unvalidated Eq. (21). This single computational sweep directly tests the load-bearing assumption without requiring new clinical data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the afterload scaling in Eqs. (20)–(22), specifically Eq. (21) for proximal impedance. The model represents the pulmonary arterial tree as N identical parallel segments, so removing a fraction α scales Z_a, R, and C by 1/(1−α), 1/(1−α), and (1−α), respectively. This is plausible for distal resistance and compliance, which are distributed across the microvascular bed, but it is much less secure for Z_a, the characteristic impedance of the proximal, extra-parenchymal pulmonary arteries. A pneumonectomy removes one main branch, yet the remaining main and lobar arteries are still present, and their characteristic impedance depends on the geometry and stiffness of the remaining tree, not simply on the number of remaining parallel segments. Table 10 shows PASP is most sensitive to Z_a, so the predicted ~49% PASP increase at α = 50% — the key signal that is supposed to distinguish afterload increase from contractility loss — rests on an unvalidated scaling exponent. The authors acknowledge in Section 5.4 that \"the assumption that pulmonary vascular resistance proportionally relates to lung volume requires validation,\" but the same caveat applies even more strongly to Z_a. If a more realistic scaling gives a smaller PASP increase, the proposed pressure differentiator could become too small to detect or even reverse.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reviews clinical evidence on right ventricular (RV) dysfunction after lung resection and presents a closed-loop 0D lumped-parameter model of the cardiovascular system with four heart chambers, four valves, and systemic and pulmonary three-element Windkessel circulations. The circulation parameters are calibrated by an exhaustive search, and the model is used to simulate two isolated mechanisms: afterload increase, modeled by scaling pulmonary resistance, impedance, and compliance with the resected lung fraction (Eqs. 20-22), and RV contractility loss, modeled by reducing E_RVmax. Local and global sensitivity analyses are performed. The central result is that most volume and pressure indices change in the same direction under both mechanisms, but RVSP, PASP, and PADP increase with afterload and decrease with contractility loss, which the authors propose as a potential clinical differentiator. The paper concludes that postoperative RV dysfunction may result from a combination of both mechanisms and claims to be the first computational model of lung resection effects on the cardiovascular system.","tokens_in":17998,"tokens_out":7407,"duration_ms":71384,"significance":"If the modeling assumptions are accepted, the paper makes a useful contribution as a first mechanistic framework for interpreting postoperative RV dysfunction after lung resection. The model is transparent, the baseline outputs lie within physiological reference ranges, and the sensitivity analyses are thorough (about 3.67 million samples for Sobol indices). The proposal that three pulmonary pressures may help distinguish afterload from contractility mechanisms is clinically interesting and falsifiable. However, the central differentiator rests heavily on an unvalidated scaling assumption for proximal impedance, and the baseline 'validation' is partly circular because the same reference ranges are used for calibration and for the comparison in Table 9. These issues are addressable and do not undermine the potential value of the framework.","major_comments":[{"comment":"The proposed differentiator—PASP increasing with afterload but decreasing with contractility loss—rests on the parallel-segment scaling for proximal impedance Z_a, Eq. (21). Table 10 shows PASP is highly sensitive to Z_a, and Section 4.3 reports a 49.1% PASP increase at alpha=50%, the largest pressure change in the afterload scenario. The scaling Z_a ∝ 1/(1−alpha) is plausible for distal resistance but much less secure for the characteristic impedance of proximal, extra-parenchymal pulmonary arteries, which depend on the geometry and stiffness of the remaining large vessels rather than simply on the number of parallel segments. The authors acknowledge in Section 5.4 that 'the assumption that pulmonary vascular resistance proportionally relates to lung volume requires validation,' but they do not extend this caveat to Z_a. Since the sign and magnitude of the PASP/RVSP difference is the central claim, the manuscript should include a robustness analysis in which the scaling exponent of Z_a is varied (e.g., Z_a ∝ (1−alpha)^{−γ} for γ in [0,2]) and should report whether the opposite trends between the two mechanisms survive, or justify the scaling with anatomic/imaging data.","section":"Section 3.3.1, Eq. (21), Table 10, Section 4.3"},{"comment":"The baseline circulation parameters are obtained by an exhaustive search that selects the sample best matching physiological reference values, and the same reference values are then used in Table 9 to show that the baseline outputs are 'reasonably close.' This is a calibration loop rather than an independent validation, so Table 9 should not be presented as evidence of external validity. I recommend stating the objective function used to select the best sample, reporting the number of near-optimal parameter sets and the spread of their outputs, and reframing Table 9 as a check of calibration consistency rather than as validation. This does not invalidate the mechanism comparison, but it affects the quantitative confidence in the baseline state.","section":"Section 3.1 'Parameter estimation', Table 9"},{"comment":"The conclusion that 'post-op RV dysfunction may be caused by both post-op afterload increase and RV contractility loss simultaneously' is not directly tested: the simulations compare the two mechanisms in isolation. A combined simulation—for example, afterload increase together with a range of E_RVmax reductions—would test whether the combined state reproduces the clinical volume pattern and whether the proposed pressure differentiators remain separable when both mechanisms act together. Please either add such simulations or soften the conclusion to state that the data are consistent with a combination but do not directly test it.","section":"Section 5.5, Section 4.3"}],"minor_comments":[{"comment":"The elastance formula is typeset ambiguously: 'E(t) = Emax−E min / maxt(...) H1 H2 + E min' should read '(Emax − Emin)/max_t(H1 H2) × H1 H2 + Emin.' Please add the missing parentheses.","section":"Eq. (2)"},{"comment":"The statement that 'the capillary resistance (Rpul) variations ... had a highly negative impact on RVSP, PASP, PADP, and RVESV' appears to contradict Table 10, which reports positive local sensitivities for PADP (0.54), PASP (0.19), and RVESV (0.3) with respect to Rpul. The afterload-increase scenario, which increases Rpul, also raises PASP, PADP, and RVESV. Please correct either the sentence or the sensitivity table.","section":"Section 5.2, Table 10, Figure 5"},{"comment":"The notation 'α=100% n/N' is confusing because (1−α) appears in the denominators of Eqs. (20)-(22). Since α is treated as a fraction, please write α = n/N (or define α as a percentage consistently and use (1−α/100) in the formulas).","section":"Section 3.3.1, Eqs. (20)-(22)"},{"comment":"The limitations listed in Section 5.4 are appropriate, but the most consequential one—the scaling of Z_a with resected volume—should be elevated from a future-work remark to a first-class limitation with a concrete sensitivity analysis, as described in Major Comment 1.","section":"Section 5.4"},{"comment":"Minor typographical issues: 'systematic' should be 'systemic' in Section 5.1, and 'In deed' in Section 5.4 should be 'Indeed.'","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid first modeling step rather than a definitive mechanistic proof. The central differentiator is plausible but depends on the unvalidated Z_a scaling (Eq. 21), so the requested gamma-sweep is the key gate for the paper's headline claim. The baseline-calibration circularity is also worth addressing, but it is secondary. I would be willing to see a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. This is the first lumped-parameter model aimed at cardiovascular changes after lung resection, and it does something useful: it isolates afterload increase from RV contractility loss in a closed-loop heart/circulation model and compares both to clinical observations. The headline result—RVSP, PASP, and PADP rise with afterload and fall with contractility loss—is a genuine model output, but it is close to a structural consequence of the equations. Increasing pulmonary resistance and impedance directly raises proximal pulmonary pressure; reducing E_RVmax directly lowers RV-generated pressure. That does not make the paper worthless, but it means the proposed clinical differentiator is a hypothesis derived from the model, not an empirical finding from it.\n\nWhat the paper does well: the clinical literature review is organized and honest about contradictory pressure findings; the model description is transparent; baseline outputs sit in physiological ranges; and the sensitivity analysis is thorough, with local and Sobol indices agreeing on R_pul and E_RVmax as the dominant parameters. The authors also explicitly flag several limitations, including the need to validate the PVR scaling with lung volume.\n\nThe soft spots are real. The most load-bearing is the scaling of the pulmonary circulation after resection. Equations (20)–(22) assume all segments are identical parallel branches, so R and Za scale as 1/(1−alpha) and C as (1−alpha). That is defensible for distal resistance and compliance, but it is much weaker for Za, the characteristic impedance of the proximal pulmonary arteries; those vessels are not distributed uniformly through lung parenchyma. Since PASP is most sensitive to Za (Table 10), the predicted ~49% PASP increase at 50% resection rests on an unvalidated exponent. The stress-test suggestion is right: a simple gamma-sweep on the Za scaling would show whether the differentiator survives under a less favorable but still plausible scaling. The authors only flag the resistance scaling in Section 5.4, not Za.\n\nThe baseline parameter search also has an element of circularity: the six circulation parameters are chosen by matching the same reference ranges used later as the qualitative standard. This is disclosed and not fatal, but it means the match with physiology is partly by construction. Independent validation, ideally against postoperative patient data, is the missing step.\n\nThis paper is for people working on lumped-parameter circulatory models and clinicians interested in post-resection RV dysfunction. It deserves a serious referee: it is a careful first attempt at a clinically relevant question, the methods are reproducible in principle, and the weaknesses are identifiable and fixable. Ask the authors for code, a robustness sweep on the Za exponent, and a clearer separation between fitted and independently verified outputs.","headline":"First lumped-parameter model of lung resection effects; the pressure differentiator is real but structurally baked in, and the Za scaling assumption needs a robustness check before the clinical claim can carry weight.","tokens_in":18520,"tokens_out":3366,"would_cite":true,"duration_ms":35267,"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":"This paper claims that right-ventricular dysfunction after lung resection probably involves both increased afterload and reduced contractility, and that three right-heart pressures can separate the two mechanisms.","keywords":["lung resection","right ventricular dysfunction","lumped parameter modelling","afterload","contractility","pulmonary circulation","time-varying elastance","sensitivity analysis"],"falsifier":"A clinical study measuring RVEF, RVSP, PASP and PADP before and after lobectomy or pneumonectomy would settle the claim: if patients whose RVEF drops show right-heart pressures that stay flat or rise, the model's predicted opposite-pressure signature of contractility loss is contradicted, whereas a mix of falling and rising pressures across patients would support it.","tokens_in":17477,"feed_emoji":"🫀","tokens_out":8050,"duration_ms":78021,"temperature":0.7,"pith_summary":"After lung resection, patients often develop right-ventricular (RV) dysfunction, but clinicians cannot tell whether the cause is extra load on the right heart or a weakened right ventricle. This paper builds a computer model of the whole circulation and runs the two suspected mechanisms separately. It finds that almost every volume and pressure index changes in the same direction under both mechanisms, except three pressures on the right side: RV systolic pressure and pulmonary artery systolic and diastolic pressure rise when afterload increases and fall when contractility drops. The paper argues that if this holds clinically, those three pressures could identify which mechanism is at work in a given patient. The broader point is that postoperative RV dysfunction may be a mixture of both mechanisms, not afterload alone.","feed_headline":"Lung surgery: model separates two causes of right-heart strain","feed_subtitle":"In simulations, RVSP, PASP and PADP rise with afterload but fall with contractility loss, offering a way to tell the mechanisms apart.","key_machinery":"The load-bearing object is a closed-loop lumped-parameter (0D) circulation model with 14 state variables: four heart chambers driven by time-varying elastance, four dynamic valves, and systemic and pulmonary three-element Windkessel afterloads. Lung resection is imposed by rescaling pulmonary impedance and resistance as $1/(1-\\alpha)$ and compliance as $1-\\alpha$, where $\\alpha$ is the fraction of lung removed; contractility loss is imposed by lowering maximum RV elastance $E_{\\max}^{\\mathrm{RV}}$. The discriminating result is that volume indices respond similarly to both perturbations, while the right-heart pressures are pushed in opposite directions by downstream resistance and upstream contraction strength.","core_discovery":"The central claim is that after lung resection, right-ventricular dysfunction can be produced by either increased pulmonary afterload or decreased RV contractility, and that these two mechanisms are distinguishable by the sign of three pressure changes. In the simulations, removing up to 50% of lung segments raises RVSP by 19.5%, PASP by 49.1% and PADP by 28.9% when afterload alone rises; cutting RV contractility to half lowers RVSP by 21.7%, PASP by 13.5% and PADP by 7.12%. Volumes such as RVEDV, RVESV and RVEF move in the same direction under both mechanisms, which is why routine imaging cannot separate them. The authors therefore propose RVSP, PASP and PADP as candidate discriminating indices.","pith_inferences":["Editorial inference: testing the model's pressure-sign prediction would require a small prospective cohort with paired echocardiographic RVEF and right-heart pressure measurements before and after resection; most existing studies report pressures and ejection fraction separately, which may be why the distinction has not been noticed.","Editorial inference: the model omits ventricular interdependence and autonomic compensation; adding these could shift the magnitude of the predicted pressure splits and may reduce how cleanly the two mechanisms separate in real patients.","Editorial inference: if vascular recruitment or remodeling invalidates the proportional afterload scaling, the modeled load increase would overestimate the real postoperative load, shifting the balance of evidence toward contractility loss as the primary driver.","Editorial inference: the same parallel-segment rescaling logic could be applied to other settings that remove vascular beds, where a similar dissociation between volume-based and pressure-based indices might appear."],"forward_implications":["If the model is correct, clinicians could use RVSP, PASP and PADP together with RVEF to infer which mechanism dominates in a given patient after lung resection.","The model reproduces the clinical pattern that RVEDV and RVESV rise and RVEF falls after resection under either mechanism, so volume measurements alone cannot settle the mechanism debate.","The simulation predicts only small LV changes under both mechanisms, consistent with the clinical observation that LV function is largely preserved after lung resection.","Because afterload alone produces the same volume trends as contractility loss, the paper's results support treating postoperative RV dysfunction as a mixed afterload-plus-contractility problem rather than purely a loading problem."],"supporting_citations":[{"why":"Supplies the time-varying elastance model used for all four heart chambers.","marker":"[30, 31]"},{"why":"Supplies the dynamic valve model with opening and closing states used for all four valves.","marker":"[32]"},{"why":"Supplies the three-element Windkessel representation of arterial afterload for both circulations.","marker":"[39]"},{"why":"Provides reference parameter ranges used to search for the pulmonary and systemic Windkessel values.","marker":"[40]"},{"why":"Supplies physiological reference ranges for ventricular volumes and ejection fraction used to select the baseline parameters.","marker":"[42]"},{"why":"Supplies normal ranges for right-heart pressures and ejection fraction used to check the baseline simulation.","marker":"[43]"},{"why":"Documents the clinical uncertainty over afterload versus contractility as the cause of postoperative RV dysfunction, motivating the two simulated mechanisms.","marker":"[12]"},{"why":"Provides the clinical RV volume and ejection-fraction changes after lung resection that the model's trends are compared against.","marker":"[4]"},{"why":"Provides the largest clinical cohort reporting postoperative pulmonary pressure and PVR increases used in the comparison.","marker":"[8]"}],"fun_headline_variants":["RV pressure directions discriminate afterload vs contractility loss","After lung resection, three pressures separate two RV stress mechanisms","Lung surgery: load and pump injuries give opposite RV pressure signs","Opposite RV pressure trends tell if lung surgery harmed load or pump"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes pulmonary vessels behave as identical parallel segments, so removing a fraction $\\alpha$ of the lung scales resistance and impedance as $1/(1-\\alpha)$ and compliance as $1-\\alpha$; if the real vasculature recruits, collapses, or remodels rather than scaling in proportion, the modeled afterload increase does not match the postoperative circulation.","fun_headline_variants_meta":{"raw":{"variants":["RV pressure directions discriminate afterload vs contractility loss","After lung resection, three pressures separate two RV stress mechanisms","Lung surgery: load and pump injuries give opposite RV pressure signs","Opposite RV pressure trends tell if lung surgery harmed load or pump"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001312,"raw_usage":{"total_tokens":5341,"prompt_tokens":934,"completion_tokens":4407,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":4349}},"tokens_in":550,"tokens_out":4407,"duration_ms":28692,"temperature":1.0,"reasoning_tokens":4349,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:21:35.634830+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clinical study measuring RVEF, RVSP, PASP and PADP before and after lobectomy or pneumonectomy would settle the claim: if patients whose RVEF drops show right-heart pressures that stay flat or rise, the model's predicted opposite-pressure signature of contractility loss is contradicted, whereas a mix of falling and rising pressures across patients would support it.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the dynamic valve model with opening and closing states used for all four valves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies physiological reference ranges for ventricular volumes and ejection fraction used to select the baseline parameters."},{"cited_title":"Closed-loop real-time simulation model of hemodynamics and oxygen transport in the cardiovascu- lar system.BioMedical Engineering OnLine, 12(1):69, July 2013","cited_arxiv_id":null,"evidence_quote":"Supplies normal ranges for right-heart pressures and ejection fraction used to check the baseline simulation."},{"cited_title":"Perioperative cardiovascular pathophysiology in patients undergoing lung resection surgery: a narrative review.British Journal of Anaesthesia, 130(1):e66–e79, January 2023","cited_arxiv_id":null,"evidence_quote":"Documents the clinical uncertainty over afterload versus contractility as the cause of postoperative RV dysfunction, motivating the two simulated mechanisms."},{"cited_title":"McCall, Alex Arthur, Adam Glass, David S","cited_arxiv_id":null,"evidence_quote":"Provides the clinical RV volume and ejection-fraction changes after lung resection that the model's trends are compared against."},{"cited_title":"Elrakhawy, Mohamed A","cited_arxiv_id":null,"evidence_quote":"Provides the largest clinical cohort reporting postoperative pulmonary pressure and PVR increases used in the comparison."}],"review_version":1}