{"id":"011d051c-0ad0-44c3-8531-2c606b64d7ab","arxiv_id":"2506.21644","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A new set of ODEs for the looping pendulum is presented with a MATLAB simulation and a four-ratio clay-bob experiment, but the simulation uses a constant angular speed and the experimental match is about 13 percent off.","lead":"Researchers derive differential equations for a looping pendulum, in which a light bob winds around a rod to hold a heavier bob, and compare simulations with clay-bob experiments. Generalist readers may care because the model claims to extend known capstan-based looping pendulum theory to strings with non-negligible mass, a step toward flexible-string robotics and nonlinear dynamics teaching.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MATLAB simulation uses a hard-coded omega=1.5 instead of the state variable dtheta, so the reported fit does not test Eq. (21).","rationale":"The paper's central claim has two necessary conditions: the ODEs are correct and the code solves them. The second condition is objectively false because the supplied code uses a constant omega=1.5 instead of the state-derived dtheta in the centripetal term of the heavy-bob equation. This is not a small correction: once the bob swings, dtheta is of order several rad/s, so the centripetal term changes by a large factor and the simulation trajectory cannot be regarded as a solution of Eq. (21). I agree with the reader's REJECT, but I put the emphasis on the code-equation mismatch rather than on the Eq. (17) kinematic identity. The Eq. (17) issue is real, yet under the reported parameters its omitted -r_r*d2theta term contributes only of order 10^-2 m/s^2 to d2y, far less than the error introduced by fixing omega=1.5. The code bug is therefore the decisive problem, and it is visible directly from the manuscript, so no external reproduction is required to see that the validation claim fails. Since the bottom-line verdict is unchanged, verdict_should_be is UNCHANGED.","tokens_in":15586,"tokens_out":12715,"duration_ms":152041,"concrete_test":"Take the shipped MATLAB code, change only the d2y expression in odefun from omega^2 to dtheta^2 (with dtheta=Y(2); the initial angular velocity is 0), keep all other parameters and the tspan values from §6.2, and recompute the four y_final values in the Sim. Data column. If any value changes by more than 10% relative to the published table, the published simulation is not a valid test of Eq. (21); if all values change by less than 10%, the hard-coded omega is numerically benign for this observable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing point is not the physical plausibility of the model but whether the submitted simulation solves the submitted equations. It does not. In the odefun in §2.6.1 and the Appendix, the d2y formula evaluates the centripetal term as Mu_d*m*l*omega^2*exp(...) with omega a scalar parameter set to 1.5 in the main script. The derivation defines omega as dtheta/dt (Section 2.2.2), so the code should use Y(2), the current angular velocity. Instead, the simulation injects a fixed 1.5 rad/s into the heavy-bob equation regardless of the state; at t=0 the initial angular velocity is dtheta0=0, so the code includes a centripetal acceleration that the initial condition says is absent. Once the bob swings, dtheta is not 1.5: for l about 0.5 m released from the horizontal, the angular speed near the bottom is of order sqrt(2g/l), roughly 6 rad/s, making omega^2 differ from the hard-coded value by a factor of order 15. The reported comparison in Figures 25/26 is therefore not a test of Eq. (21); it is a fit of a different, unspecified dynamical system. This alone invalidates the central claim as presented.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models the looping pendulum, in which a light bob wraps a massive string around a cylindrical rod and thereby supports a heavier bob that falls vertically. The author separates the dynamics into a rotating component and a non-rotating component, and derives three coupled equations: Eq. (6) for the angular acceleration of the light bob, Eq. (21) for the vertical acceleration of the heavy bob, and Eq. (22) for the time-varying free length of the string. A MATLAB ode45 simulation of these equations is compared with video-tracked experiments for four mass ratios (7.88, 12.92, 18.16, 29.03), and the paper claims close agreement between simulation and experiment, accrediting the derived theory.","tokens_in":15716,"tokens_out":4094,"duration_ms":43337,"significance":"If the derivation and validation were sound, the paper would offer a useful extension of existing massless-string looping-pendulum models to include string mass, packaged as a compact ODE system with a falsifiable prediction of fall distance versus mass ratio. The manuscript has strengths: the full MATLAB code is provided, the experimental procedure is described in detail, and the author acknowledges limitations such as low frame rate and a small number of trials. However, as detailed below, the central validation is undermined by a physically incorrect radial-acceleration assumption and by a simulation that does not solve the stated equations; the claimed agreement with experiment therefore does not support the paper's central claim.","major_comments":[{"comment":"The radial acceleration of the light bob is set to a_r = d2y/dt2 - l*omega^2, with the justification that m and M are connected by an inextensible string and therefore have the same translational acceleration. This is not correct: the free length l = L + y - rr(pi+theta) changes at the rate dl/dt = dy/dt - rr*dtheta/dt, so the radial coordinate of m is l, and its radial acceleration is d2l/dt2 - l*omega^2 = d2y/dt2 - rr*d2theta/dt2 - l*omega^2. The term -rr*d2theta/dt2 is omitted. This omission propagates through Eq. (18) into Eq. (21), so the principal dynamical equation for the heavy bob does not follow from the stated assumptions.","section":"§2.3.3.2, Eq. (17)"},{"comment":"The MATLAB function odefun evaluates the centripetal contribution to d2y as Mu_d*m*l*omega^2*exp(Mu_d*(pi+theta)), where omega is a scalar parameter passed from the main script and set to omega=1.5. In the derivation, omega denotes dtheta/dt (Section 2.2.2), so the code should use the state variable dtheta = Y(2), not a fixed constant. At the initial condition dtheta0=0 the hard-coded omega injects a nonzero centripetal acceleration, and during the motion dtheta differs substantially from 1.5 rad/s. The simulation therefore does not solve Eqs. (6), (21), and (22); it solves a different dynamical system, so Figures 25 and 26 are not a test of the derived theory.","section":"§2.6.1 and Appendix (odefun)"},{"comment":"For each trial, the simulation time span is set to the experimentally measured stopping time (0.3 s, 0.35 s, 0.4 s, 0.37 s), as stated in Section 6.2 and implemented in the main script. The reported predicted fall distance is y(end), evaluated at exactly the time when the experiment was observed to stop. Since the model has no stopping criterion, this is not a closed prediction of the total fall distance; it is an evaluation at a data-supplied time, so the comparison in Figure 25 is partly circular and cannot validate the model's predictive content.","section":"§6.2 and main script (tspan)"}],"minor_comments":[{"comment":"The torque from the string's weight is divided by two with the comment that the string mass is 'averaged at the center of the string'; this should be justified explicitly, since the torque distribution along a finite-length string is not self-evidently equivalent to a point mass at its midpoint.","section":"§2.2.2.1, Eq. (1)"},{"comment":"After dividing by e^{-mu_d phi}, the exponential in the constant term changes from e^{-mu_d phi} to e^{mu_d phi} without explanation; the algebra should be checked and the sign convention clarified.","section":"§2.3.2.4, Eq. (14)"},{"comment":"The text states the rod radius is 0.003 m, but the MATLAB code sets rr = 0.0003. This factor-of-ten discrepancy affects Eq. (22) and the wrapped-length correction; it should be resolved.","section":"§5.2 and main script"},{"comment":"The trajectory comparison for the heavy bob is described as having 'not much resemblance' between experiment and simulation, which is in tension with the abstract's claim of very close agreement; this discrepancy should be addressed.","section":"§6.3"}],"recommendation":"reject","confidential_remarks":"The manuscript appears to be a school research report. The identified issues—especially the incorrect radial-acceleration assumption and the hard-coded omega replacing the state variable—mean the reported numerical-experimental agreement does not validate the derived equations. While the topic is suitable for a pedagogical journal, the present form does not meet the standards of a research publication in physics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The best thing about this paper is the attempt to generalize the capstan-based looping pendulum model to a string with mass. Equations (6), (21), and (22) are a genuine extension of the massless-string treatment in reference (4), and the mass-ratio versus fall-distance scan is a reasonable experimental question. The paper is clearly written, includes full code and data tables, and is honestly transparent about its limitations: only four mass ratios, a 30 fps camera, and the string-bending problem. Those are real assets for a school-level research report.\n\nBut the central claim—that the experimental findings fit the simulation and thereby accredit the theory—does not survive contact with the submitted code. In the odefun, the centripetal term in the d2y equation uses omega^2 with omega fixed as a scalar 1.5, not the state variable dtheta. The paper derives omega as dtheta/dt, so the simulation is not a solution of equation (21). At t=0 the initial angular velocity is zero, yet the code injects a nonzero centripetal acceleration. Later, the angular speed can be several rad/s, so a fixed 1.5 is not a good proxy. This alone invalidates the reported fit as evidence for the model.\n\nThe derivation also has a kinematic error. The paper sets the radial acceleration of the light bob to d2y/dt2 - l*omega^2. That ignores the time variation of the unwrapped length l. Since l = L + y - rr(pi+theta), the second derivative of l involves theta'', not just y''. The approximation might be defensible in a limit, but the paper does not state or justify that limit, and it is load-bearing because it enters the expression for tension and hence equation (21).\n\nA third, smaller issue: the simulation time span is taken directly from the experimentally measured stopping time for each trial. So the final y(end) is not an independent prediction; it is partly a function of the data. The authors are upfront about this, but it further weakens the validation.\n\nWho is this paper for? It could be a useful educational example of the looping pendulum, and the experimental data might be worth secondary use. But as a research contribution it does not work. The code-equation mismatch and the kinematic gap are visible from the text alone, and the 13% systematic discrepancy is blamed on string bending without a quantitative check. I would not send this to peer review in its current form; it needs a corrected derivation and a code that actually solves the stated equations before it deserves referee time.","headline":"The paper's finite-string-mass generalization of the looping pendulum equations is a plausible extension, but the supplied MATLAB code does not solve those equations and the derivation has a real kinematic gap, so the reported validation does not hold up.","tokens_in":16369,"tokens_out":2587,"would_cite":false,"duration_ms":30288,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives three coupled differential equations for the looping pendulum and reports that simulations built on them closely match video-tracked measurements of the heavy bob's fall.","keywords":["looping pendulum","IYPT","rotational mechanics","torque","friction","capstan equation","nonlinear dynamics","differential equations"],"falsifier":"Track the heavy bob's fall $y(t)$ and the free hanging length $l(t)$ separately with high-speed video; if the second time derivative of $l$ differs from the second derivative of $y$, the assumption used to eliminate the string tension is false and the two-equation system is not a faithful model.","tokens_in":15250,"feed_emoji":"🪢","tokens_out":7319,"duration_ms":72516,"temperature":0.7,"pith_summary":"The paper attempts to establish that the looping pendulum, in which a light bob sweeps around a cylindrical rod while supporting a heavier bob, is governed by a system of three ordinary differential equations. The derivation splits the motion into a rotating part and a non-rotating part, and treats the tension in the wrapped string through the capstan equation for a string with non-negligible mass. The paper solves these equations numerically and compares the predicted vertical fall of the heavy bob with video-tracked experiments across four mass ratios. It reports that the experimental fall distances follow the same nonlinear rising trend as the simulation, and takes this agreement as evidence that the derived theory is valid.","feed_headline":"Three equations predict the looping pendulum's fall","feed_subtitle":"A capstan-based model with a massed string matches video-tracked fall distances across mass ratios.","key_machinery":"The load-bearing object is the capstan-equation tension profile along the string wrapped around the rod. The paper derives, for an infinitesimal string element, a differential relation between tension and wrap angle, integrates it to find an exponential factor $e^{\\mu_d(\\pi+\\theta)}$, and substitutes this into the heavy bob's vertical equation of motion. The other two equations are the angular equation obtained by writing torque as the time derivative of angular momentum with a time-varying moment of inertia, and the geometric relation $l = L + y - r_r(\\pi+\\theta)$ that couples the free length to the fall distance and the wrap angle.","core_discovery":"The paper's central claim is that the looping pendulum is described by three equations: an angular equation for the light bob, a vertical equation for the heavy bob whose acceleration contains an exponential capstan factor, and a geometric constraint for the free string length. The equations extend earlier work by giving the string a non-negligible mass, so the light bob's mass does not need to be accounted for separately. The paper reports that numerical solutions of these equations reproduce the experimentally observed trend, in which the heavy bob's vertical fall increases nonlinearly with the mass ratio, and it takes this agreement as validation of the derived theory.","pith_inferences":["The published numerical solver fixes $\\omega$ at 1.5 rad/s instead of taking it from the state variable $d\\theta/dt$; re-running with $\\omega$ extracted from the solution would show whether the reported agreement comes from the equations themselves or from that fixed parameter.","The derivation sets the light bob's radial acceleration equal to the heavy bob's vertical acceleration minus a centripetal term, which treats the free length $l$ as if it changed at the same rate as $y$; a corrected treatment with a separate $d^2l/dt^2$ term could alter the predicted fall distances.","The model's structure matches tendon-pulley systems in robotics, where a thin string wraps around a pulley and friction governs the tension; testing the same equations on a robotic tendon with measured parameters would be a direct transfer."],"forward_implications":["The heavy bob's fall distance rises nonlinearly with the mass ratio, with the increase flattening at larger ratios; this is the trend the paper's data and simulation both show.","A string with non-negligible mass can be modelled without separately tracking the mass of the light bob's string segment, because the capstan tension integration accounts for the string's linear density.","The same three-equation system yields both the spiral trajectory of the light bob and the vertical trajectory of the heavy bob, so a single simulation describes both components.","Because the tension scales as $e^{\\mu_d(\\pi+\\theta)}$, small changes in the dynamic friction coefficient or the wrap angle translate into large changes in the force supporting the heavy bob."],"supporting_citations":[{"why":"Introduces the looping pendulum problem from the 2019 IYPT that the paper sets out to model.","marker":"(1)"},{"why":"Previous massless-string derivation of the looping pendulum that this work generalises to a massed string.","marker":"(4)"},{"why":"States the capstan equation used to derive the exponential tension profile around the rod.","marker":"(5)"},{"why":"Supplies the dynamic friction coefficient of the cotton string used in the numerical simulation.","marker":"(11)"},{"why":"Earlier research on the looping pendulum whose trajectory plots the paper compares with its own.","marker":"(13)"}],"fun_headline_variants":["Three equations catch the looping pendulum's fall","Massed string model nails looping pendulum drop","Looping pendulum's nonlinear drop, now predictable","Capstan math predicts looping pendulum's fall"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the free hanging length of string changes exactly as fast as the heavy bob falls, which lets it equate the light bob's radial acceleration with the heavy bob's vertical acceleration minus the centripetal term.","fun_headline_variants_meta":{"raw":{"variants":["Three equations catch the looping pendulum's fall","Massed string model nails looping pendulum drop","Looping pendulum's nonlinear drop, now predictable","Capstan math predicts looping pendulum's fall"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000127,"raw_usage":{"total_tokens":1043,"prompt_tokens":803,"completion_tokens":240,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":184}},"tokens_in":419,"tokens_out":240,"duration_ms":3173,"temperature":1.0,"reasoning_tokens":184,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:35:02.743314+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track the heavy bob's fall $y(t)$ and the free hanging length $l(t)$ separately with high-speed video; if the second time derivative of $l$ differs from the second derivative of $y$, the assumption used to eliminate the string tension is false and the two-equation system is not a faithful model.","supporting_citations":[],"review_version":1}