{"id":"b906804a-286e-47dd-8cc0-060cf895cc27","arxiv_id":"1908.04984","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Charged obstacle size and shape set a threshold and control amplitude, velocity, and number of precursor solitons in a supersonic dusty plasma flow.","lead":"This dusty plasma experiment shows that the height, width, and slope of a charged obstacle's potential hill control whether precursor solitons form and how fast and large they are. A sharp potential step emits solitons ahead of the flow, a gentle slope only leaves wakes, and no solitons form below a threshold hill size.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (1) as printed is a forced linear KdV equation (the nonlinear advection term is missing), so the numerical soliton solutions and the claimed experimental–f-KdV agreement cannot be reproduced from the stated model.","rationale":"The reader identified the indirect dust-microprobe measurement as the weakest assumption, and that concern is legitimate: the reported height–width relation and threshold depend on an energy-conservation reconstruction that is not independently verified, and the triangular-object source width is hand-set. However, the printed f-KdV equation is a more specific and more easily checkable problem. As written, Eq. (1) has no nonlinear term, so it cannot support KdV solitons; this directly undermines the numerical-agreement component of the central claim. The experimental observations could in principle stand on their own, and the missing term may be a typographical error, but the manuscript as submitted does not allow a reader to verify the simulations. I therefore retain the reader's CONDITIONAL verdict: the paper should not be fully accepted until the model equation is corrected and the numerical results are reproduced with the correct equation, and until the potential-hill measurement is clarified with uncertainties or a direct calibration. This does not amount to rejection because the empirical trends are internally consistent and benchmarked against prior work, and the requested re-run would likely settle the issue without invalidating the observations if the typo is confirmed.","tokens_in":11792,"tokens_out":6562,"duration_ms":72478,"concrete_test":"Independently implement Eq. (1) exactly as printed (linear forced dispersive equation) with the Gaussian source function, A=4, and the supersonic source speed used in Sec. III A, and compare the output to Fig. 4(b). If no soliton-like upstream structures appear—which is expected for the linear equation—then ask authors to supply the corrected f-KdV equation (including the missing nonlinear term) and to re-run Figs. 8, 11, and 13 with it; alternatively, compare Eq. (1) with Ref. 40 to identify whether the nonlinear term was dropped in typesetting.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (1) reads ∂nd1/∂t + A∂nd1/∂ξ + (1/2)∂^3nd1/∂ξ^3 = (1/2)∂S2/∂ξ. Written this way, it is a forced linear dispersive equation: there is no term of the form nd1∂nd1/∂ξ, which is the nonlinearity that creates KdV solitons and their amplitude–width relation. A linear equation cannot produce the solitary structures shown in Fig. 4(b), the amplitude/width scalings in Fig. 8, or the threshold behaviour used to support the experimental trends in Secs. III B and III C. The text does not provide an alternative nonlinear term or an explanation of how 'solitons' arise from Eq. (1). This matters because the abstract and Sec. IV rest part of the central claim on 'qualitative agreement' with f-KdV simulations. If this is merely a typesetting omission, the paper needs to print the correct equation and re-verify the simulations; if it is the equation actually solved, then the numerical support is invalid as stated. This is more directly load-bearing than the already valid concern about indirect potential-hill measurement, because it affects the interpretability of every figure that claims f-KdV support.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports experiments in a Π-shaped DPEx device in which a supersonic dusty-plasma flow is driven over a charged copper wire or a triangular object, and the excitation of precursor solitons and wakes is studied as a function of the size and shape of the electrostatic potential hill. By varying the voltage drop across a resistance in series with the wire, the authors change the hill height and width, observing that decreasing hill height reduces the amplitude, velocity, and number of precursor solitons while increasing their width, with a threshold at Vwg=135 V below which no solitons appear. Replacing the wire with a triangular object, they find that a linearly rising slope produces only wakes, whereas a sharp edge produces both solitons and wakes. All experimental trends are compared qualitatively with numerical solutions of a forced Korteweg-de Vries (f-KdV) model using Gaussian or half-Gaussian source terms.","tokens_in":11988,"tokens_out":2704,"duration_ms":27767,"significance":"If the experimental trends hold, the paper would extend the recently established precursor-soliton phenomenon in complex plasmas by showing systematic control through obstacle size and shape, with potential relevance to space-debris and ionospheric interactions. The work has concrete strengths: it benchmarks against the earlier experiment of Jaiswal et al. (Ref. 29), uses PIV for flow-speed measurement, reports a threshold observation, and provides a physically motivated source-function modeling approach. However, the numerical support is currently compromised by the printed model equation, and the key potential-hill diagnostics rest on an indirect technique without error characterization. The central qualitative claims are plausible but require correction of the equation and more careful uncertainty reporting before the agreement with f-KdV can be considered established.","major_comments":[{"comment":"Eq. (1) as printed is a forced linear advection-dispersion equation: ∂nd1/∂t + A∂nd1/∂ξ + (1/2)∂^3nd1/∂ξ^3 = (1/2)∂S2/∂ξ. It contains no nonlinear term such as nd1∂nd1/∂ξ, so it cannot produce the solitary structures shown in Fig. 4(b), the amplitude-width scalings in Fig. 8, or the threshold behavior discussed in Sec. III B. The numerical solutions therefore cannot be reproduced from the model as stated. If the nonlinear term was omitted in typesetting, the correct f-KdV equation must be printed and the simulations re-verified; if Eq. (1) is the equation actually solved, the claimed qualitative agreement with f-KdV is unsupported.","section":"Sec. III A, Eq. (1)"},{"comment":"The height and width of the potential hill are inferred from the dust-microprobe energy-conservation technique of Ref. 41, with no independent calibration or cross-check, and Fig. 5(b) reports a power-law fit h ∼ w^−3.13 without error bars or goodness-of-fit statistics. The threshold Vwg=135 V is a single observation. Since these inferred values become the source amplitudes and widths for the f-KdV simulations, the quantitative scalings in Figs. 7 and 8 rest on uncharacterized measurement uncertainty.","section":"Sec. III B, Fig. 5"},{"comment":"The half-Gaussian source width is chosen as W/2 = 6 based on the quoted sheath size rather than on direct measurement of the triangular-object potential profile, and the wake-only versus soliton-plus-wake outcome is largely encoded in the orientation of the source function. The authors should demonstrate that the wake-only result for the rising slope is robust to a range of half-Gaussian widths and slope values, and should provide quantitative comparisons (e.g., amplitudes and wavelengths) rather than only qualitatively similar images.","section":"Sec. III C, Figs. 9-13"}],"minor_comments":[{"comment":"The axis label 'Voltage acorss Wire' contains a typo; it should read 'Voltage across Wire'.","section":"Fig. 5(a)"},{"comment":"References to 'Eq. I' should be 'Eq. (1)' to match the equation numbering used in the text.","section":"Throughout"},{"comment":"The paper states that the structures were confirmed to be solitons by checking the constancy of amplitude × width^2, but no quantitative data or plot is provided; please include this evidence.","section":"Sec. III A"},{"comment":"The plots in Figs. 6 and 7 show no error bars or point-to-point scatter; the authors should state the number of repeated experiments and the run-to-run variability.","section":"Sec. III B, Figs. 6-7"},{"comment":"The numerical method for solving the f-KdV equation (discretization, time step, boundary conditions) is not described; a brief statement of the scheme would improve reproducibility.","section":"Sec. II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports an interesting experimental extension of precursor soliton studies, but the printed model equation is linear and cannot generate solitons, which is a load-bearing issue for the numerical comparison. The authors should also strengthen the uncertainty analysis of the potential-hill measurement. If the equation is corrected and the simulations re-verified, the paper could be acceptable after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper reports a genuine extension of the group's 2016 precursor-soliton result: for the first time, it systematically varies the potential hill's height and width (via a variable resistance) and the obstacle's shape (wire vs. triangular wedge), and finds a threshold below which no solitons appear, plus a clean shape-dependent dichotomy where a linearly rising slope gives only wakes while a sharp edge gives solitons plus wakes. The images are consistent, the benchmark against prior work is sensible, and the qualitative trends in amplitude, velocity, width, and soliton number are internally coherent. That is a real contribution for the dusty-plasma and space-debris analog community.\n\nThe soft spots are real but mostly fixable. The biggest issue is Eq. (1): as printed, it is a forced linear dispersive equation with no nonlinear term, so it cannot produce the soliton solutions shown in Figs. 4(b), 8, 11(b), or 13(b). This is likely a typesetting omission rather than what the authors actually solved, but it is load-bearing because the abstract and Sec. IV rest part of the claim on \"qualitative agreement\" with f-KdV simulations. The paper must print the correct equation and confirm the simulations use it. If the equation shown is actually what was solved, the numerical support is invalid and the claim collapses to experiment-only.\n\nOther concerns are less severe but worth naming. The potential-hill height and width are measured with the dust-microprobe technique of Ref. 41, which assumes energy conservation; an independent or calibrated measurement would strengthen things. There are no error bars or replicate statistics anywhere, the threshold at Vwg = 135 V is a single point, and the soliton criterion (amplitude × width^2 constant) is asserted but not shown. On the simulation side, the source terms are partly hand-chosen: the half-Gaussian width for the triangular object is set from a sheath estimate, and the wake-vs-soliton outcome is substantially encoded in that choice. No data or code are provided.\n\nNone of this kills the experimental result, which stands on its own images and internal consistency. But the printed equation error must be corrected, and the authors should provide uncertainty estimates, clarify the measured-to-model mapping, and ideally release raw data or a direct potential measurement.\n\nWho is this for? Dusty-plasma experimentalists and anyone working on charged-object-plasma interactions. It deserves a serious referee, but the referee should demand the fixes above before acceptance.","headline":"A systematic experimental study of how obstacle hill size and shape control precursor solitons in a dusty plasma, but the printed f-KdV equation is missing its nonlinear term, so the numerical support needs correction before the paper can be fully trusted.","tokens_in":12590,"tokens_out":1942,"would_cite":false,"duration_ms":21935,"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":"The height and shape of a charged obstacle's electrostatic hill control the precursor solitons emitted ahead of a supersonic dusty-plasma flow.","keywords":["precursor solitons","dusty plasma","forced Korteweg-de Vries equation","potential hill","dust acoustic waves","wakes","supersonic flow","solitary waves"],"falsifier":"Repeat the resistor scan while independently measuring the potential profile around the wire; if precursor solitons appear for $V_{wg}$ above 135 V, or if the measured height-width relation departs from $h\\sim w^{-3.13}$, the threshold and size-scaling claims would be contradicted.","tokens_in":11497,"feed_emoji":"🌊","tokens_out":8534,"duration_ms":85133,"temperature":0.7,"pith_summary":"This paper claims that when a supersonic dusty plasma flows over a charged obstacle, the obstacle's electrostatic potential hill determines how many precursor solitons are emitted and how fast, large, and wide they are. Lowering the hill reduces the amplitude, velocity, and number of the solitons and widens them; below a threshold no solitons are excited at all. The shape of the hill matters as well: a linearly rising potential slope produces only downstream wakes, while a sharp potential rise produces both upstream solitons and downstream wakes. The experimental trends are reproduced qualitatively by forced Korteweg-de Vries simulations with Gaussian and half-Gaussian source terms. This matters because precursor solitons are a nonlinear wave signature of a moving disturbance, and the result shows how the disturbance's geometry sets whether that signature exists.","feed_headline":"Steep, tall barriers emit precursor solitons; low hills do not","feed_subtitle":"A dusty-plasma experiment shows soliton size, speed, count, and existence hinge on the obstacle's potential hill.","key_machinery":"The central object is the forced Korteweg-de Vries (f-KdV) equation for the dust density perturbation, with a moving source term $S_2$ that represents the electrostatic potential hill of the charged object. The source is taken as a Gaussian for the cylindrical wire and as a half-Gaussian for the triangular obstacle; supersonic motion of this source produces upstream solitons and downstream wakes, while reducing the source amplitude or broadening it weakens the solitons and eventually leaves only wakes. A companion measurement technique uses the dust particles themselves as microprobes: tracking their trajectories over the hill and applying energy conservation gives the hill's height and width, linking the measured voltage drop $V_{wg}$ to the source parameters used in the simulations.","core_discovery":"For a supersonic flow of dust fluid over a charged wire in a dusty plasma, the height and width of the electrostatic potential hill act as tunable controls for the emitted forced dust-acoustic solitary waves. As the voltage drop across a series resistor increases, the potential hill lowers and widens; in response, the precursor solitons become smaller in amplitude, slower, fewer, and wider, and at a threshold voltage of about 135 V the excitation stops entirely. Replacing the wire with a triangular object shows that a linearly rising potential slope excites only wakes, whereas the same object reversed, presenting a sharp potential jump, excites both upstream precursor solitons and downstream wakes. These observations are qualitatively matched by numerical solutions of the forced Korteweg-de Vries equation when the source term is a Gaussian for the wire and a half-Gaussian for the triangular object, with source amplitude and width chosen from the measured height-width relation of the potential hill.","pith_inferences":["A natural extension would be to sweep the ramp angle of the triangular object continuously and map the boundary between wake-only and soliton-plus-wake regimes, a scan the paper does not report.","If the sharp-gradient condition carries over to other fluids, ship-generated precursor solitons in shallow water should weaken or vanish as hull draft or blockage is reduced, making the dusty-plasma experiment a small-scale analogue of the coastal problem.","The measured height-width power law ($h\\sim w^{-3.13}$) for the wire's sheath could be tested at other discharge pressures and densities; if it holds, source parameters for the f-KdV model could be predicted without rerunning the microprobe measurement."],"forward_implications":["Because lowering the potential hill reduces soliton amplitude, velocity, and number while increasing width, the obstacle's height acts as a control parameter for the emitted nonlinear wave train.","The existence of a threshold hill height means precursor-soliton emission is not automatic for supersonic flow; below the threshold the fluid simply flows over the hill.","A linearly rising potential slope is insufficient to excite upstream solitons; only a sharp potential rise produces both precursor solitons and downstream wakes.","The f-KdV model with a Gaussian or half-Gaussian source term reproduces the experimental trends, so source shape must be included alongside flow speed in predictions of precursor excitation.","The results support interpreting natural observations, such as disturbances from objects moving through ionospheric or space plasmas, as size- and shape-dependent precursor-soliton events."],"supporting_citations":[{"why":"Provides the first experimental observation of precursor solitons in a flowing dusty plasma, which the present work uses as its benchmark.","marker":"[29]"},{"why":"Supplies the dust-particle microprobe technique used to measure the height and width of the potential hill.","marker":"[41]"},{"why":"Gives the forced KdV model equation and its solution method, from which the numerical source-term results are obtained.","marker":"[40]"},{"why":"Reports full fluid simulations confirming precursor soliton and wake excitation by a moving charged object.","marker":"[30]"},{"why":"Reports molecular-dynamic simulations confirming fore-wake excitations from moving charged objects.","marker":"[31]"},{"why":"Describes the DPEx device and its diagnostics, the experimental platform for all measurements.","marker":"[32]"},{"why":"Establishes the hydrodynamic precedent that upstream solitons from moving disturbances obey a forced KdV description and are Galilean invariant.","marker":"[4]"}],"fun_headline_variants":["Potential hill shape and height tune precursor soliton emission","Barrier steepness dictates soliton generation in dusty plasma","Precursor solitons vanish below obstacle potential threshold","Obstacle shape flips soliton wakes in dusty-plasma flow","Soliton emission tracks obstacle height, width, and shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The microprobe method's energy-conservation conversion from tracked particle motion to the height and width of the potential hill must faithfully represent the true electrostatic profile, because those numbers become the source terms that the numerical comparison depends on.","fun_headline_variants_meta":{"raw":{"variants":["Potential hill shape and height tune precursor soliton emission","Barrier steepness dictates soliton generation in dusty plasma","Precursor solitons vanish below obstacle potential threshold","Obstacle shape flips soliton wakes in dusty-plasma flow","Soliton emission tracks obstacle height, width, and shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3512,"prompt_tokens":1021,"completion_tokens":2491,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":2408}},"tokens_in":637,"tokens_out":2491,"duration_ms":17390,"temperature":1.0,"reasoning_tokens":2408,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:27:35.438186+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the resistor scan while independently measuring the potential profile around the wire; if precursor solitons appear for $V_{wg}$ above 135 V, or if the measured height-width relation departs from $h\\sim w^{-3.13}$, the threshold and size-scaling claims would be contradicted.","supporting_citations":[{"cited_title":"Jaiswal , author P","cited_arxiv_id":null,"evidence_quote":"Provides the first experimental observation of precursor solitons in a flowing dusty plasma, which the present work uses as its benchmark."},{"cited_title":"Arora , author P","cited_arxiv_id":null,"evidence_quote":"Supplies the dust-particle microprobe technique used to measure the height and width of the potential hill."},{"cited_title":"Sen , author S","cited_arxiv_id":null,"evidence_quote":"Gives the forced KdV model equation and its solution method, from which the numerical source-term results are obtained."},{"cited_title":"Kumar Tiwari \\ and\\ author A","cited_arxiv_id":null,"evidence_quote":"Reports full fluid simulations confirming precursor soliton and wake excitation by a moving charged object."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports molecular-dynamic simulations confirming fore-wake excitations from moving charged objects."},{"cited_title":"Jaiswal , author P","cited_arxiv_id":null,"evidence_quote":"Describes the DPEx device and its diagnostics, the experimental platform for all measurements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the hydrodynamic precedent that upstream solitons from moving disturbances obey a forced KdV description and are Galilean invariant."}],"review_version":1}