{"id":"2c62a484-0cf9-4e7c-9b11-ec5012b14316","arxiv_id":"1908.08354","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The Airy distribution is measured in a colloidal experiment for the first time, and new large-deviation results including conditional position distributions and two third-order phase transitions are derived.","lead":"This paper reports the first direct experimental measurement of the Airy distribution, the probability law for the area under a Brownian excursion, using dilute colloidal particles. It also derives new conditional position statistics and identifies two third-order dynamical phase transitions in the large deviation function.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. D5 is algebraically inconsistent: at the stated Xc2 it violates the area constraint, so the intermediate-branch rate function and the second third-order transition are unsupported as printed.","rationale":"The paper has genuinely strong components: the experimental histogram of the Airy distribution, the Vervaat-transform construction, the DV and OFM derivations that reproduce known tails, and the simulation checks in Figs. 2, 3, and 5. The central weak spot is the large-deviation calculation of the conditioned position distribution in Appendix D, which is also the source of the headline 'two third-order dynamical phase transitions.' The reader identified the tangent/vanishing-interval construction as the weakest assumption. My independent check found a more specific and concrete failure: Eq. D5, which determines the contact point τ of the intermediate solution, is algebraically inconsistent with the area constraint at the stated critical point Xc2. This is not a disagreement with consensus or a stylistic concern; it is an internal inconsistency in the printed derivation. It directly affects the quantitative prediction (D7) and the proof that the second transition is third order. The error is correctable, and the final claim might survive after correction, so the appropriate verdict remains conditional rather than rejection. I therefore leave the reader's verdict unchanged but would update the stated condition: the authors must correct Eq. D5 and re-verify the intermediate-branch rate function and the order of the transition at Xc2.","tokens_in":16720,"tokens_out":24738,"duration_ms":244729,"concrete_test":"Recompute the area of the trajectory defined by Eqs. D4 and D5 at Xc2 = 3A/(2t). If ∫ x dt ≠ A, the printed D5 is invalid. Then solve the area constraint 2l^3 − 3(2ã−1)l^2 − 1 = 0 for τ, replace D5, recompute the action difference (D7) at Xc2, and check whether the third derivative of the large-deviation function still jumps at Xc2. A direct numerical minimization of the action for a few X around Xc2 with, say, T=1, t=0.4, and large A would independently confirm which branch is optimal and whether the second transition is indeed third order.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim of two third-order dynamical phase transitions rests on the Appendix D variational construction, and the weak point is not the tangent construction in principle but its implementation. For the intermediate optimal trajectory (D4), the parameter τ is fixed by the area constraint. Integrating D4 and setting the area equal to A gives X [ (τ−t)/3 + t/2 − t^3/(6(τ−t)^2) ] = A. With l=(τ−t)/t and ã=A/(Xt), this is 2l^3 − 3(2ã−1)l^2 − 1 = 0. At the claimed second critical value Xc2 = 3A/(2t), one has ã = 2/3, and the only positive root is l = 1, so τ = 2t. The printed Eq. D5, however, gives τ/t ≈ 2.054 at ã = 2/3, i.e. l ≈ 1.054. Substituting this into D4 yields an excursion whose area is ≈ 1.052A, not A (for example, T=1, t=0.4, A=1, X=3.75). Thus the trajectory used to evaluate the intermediate-regime action difference (D7) does not satisfy the conditioning constraint. Consequently, the printed equations do not establish equality of the intermediate and supercritical branches at Xc2, and the third-order nature of the second transition does not follow as written. The correct τ is obtained from the cubic above; with that correction the transition may survive, but the manuscript currently contains an internal mathematical error in the core derivation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first experimental measurement of the Airy distribution (AD) using a dilute colloidal suspension, comparing the measured histogram of areas under Vervaat-transformed Brownian excursions with the exact scaling function in Eq. (3). The authors combine the Donsker-Varadhan (DV) formalism and the optimal fluctuation method (OFM) to derive, respectively, the small- and large-area tails of the AD and the single-time position distribution of a Brownian excursion conditioned on a specified area. For small areas they find that this conditional distribution coincides with the Ferrari-Spohn distribution and provide an exact mapping to a parabolic absorbing wall; for large areas they derive Gaussian fluctuations in the subcritical regime and non-Gaussian forms beyond critical values Xc1(t) and Xc2(t), which they interpret as two third-order dynamical phase transitions that merge into a single transition at t=T/2. An appendix derives the tails of the distribution of the area under the square of a Brownian excursion. I checked the algebraic concern raised about Eq. (D5): when the first term is read as a fraction rather than as a product, Eq. (D5) is the correct solution of the area-constraint cubic, and at X=Xc2 (namely ã=2/3) it gives τ/t=2, so the intermediate-regime trajectory satisfies the area constraint; the apparent inconsistency is therefore not present.","tokens_in":16999,"tokens_out":17389,"duration_ms":159642,"significance":"If the results hold, this is a significant paper: it provides the first direct laboratory verification of the Airy distribution, places the known asymptotic tails in a physically transparent large-deviation framework, explains the Ferrari--Spohn coincidence through an exact canonical/microcanonical mapping, and makes falsifiable predictions about third-order dynamical phase transitions in the conditional position distribution. The derivations in Appendices B and D are first-principles and reproduce known tails without adjustable parameters, and the experimental comparison in Fig. 1 is parameter-free. The algebraic error that initially appeared to undermine the second phase transition does not survive scrutiny; what remains are local presentation issues, most notably a cross-reference typo in Appendix D and a need for clearer typesetting of Eq. (D5).","major_comments":[],"minor_comments":[{"comment":"The sentence 'where τ is given by (D6)' should read 'where τ is given by (D5)'; Eq. (D6) is the supercritical trajectory and does not define τ.","section":"Appendix D, just below Eq. (D7)"},{"comment":"The first term of Eq. (D5) is line-broken in the present text and is easy to misread as a product. Please typeset it unambiguously as the fraction (2ã−1)^2 / [2(√((2ã−1)^3+1)+1)^{2/3}], or equivalently present τ/t in the explicitly symmetric Cardano form ã+1/2 + (1/2)[(√((2ã−1)^3+1)+1)^{2/3} + (√((2ã−1)^3+1)−1)^{2/3}], so that the area constraint at X=Xc2 is visibly satisfied.","section":"Appendix D, Eq. (D5)"},{"comment":"The claim that the second critical line Xc2(t) is a third-order transition is stated without an explicit check that the first two derivatives of the large-deviation function are continuous and the third derivative jumps there; adding a one-line verification would make the central claim easier to audit.","section":"Main text, after Eqs. (14)--(16)"},{"comment":"The abstract says there are two singularities for large areas, whereas for the symmetric observation time t=T/2 the two transitions merge into one; please add a qualifier such as 'for generic observation times' to avoid ambiguity.","section":"Abstract and Sec. II"},{"comment":"The experimental conditional averages are based on only 22 (large-area) and 200 (small-area) trajectories; adding error bars or confidence bands would make the agreement with Eqs. (8) and the constant-area prediction quantitative.","section":"Fig. 2(a)"},{"comment":"The broadening of the measured histogram near the maximum is attributed to particle polydispersity; a quantitative estimate, for example a convolution over a distribution of diffusivities, would strengthen this explanation.","section":"Fig. 1"}],"recommendation":"minor_revision","confidential_remarks":"The paper is well within the scope of the journal and the experimental novelty is supported by the data shown. I have no further confidential concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—this one is worth your time but needs a careful look at the appendix. The headline accomplishment is real: the authors have measured the area distribution of Brownian excursions in a colloidal experiment, and the histogram in Fig. 1 is a credible first direct measurement of the Airy distribution. The Vervaat construction from tracked trajectories is sensible, and the agreement with the closed-form f(ξ) is good, modulo a slight broadening near the peak that polydispersity can plausibly explain. That part should survive review.\n\nThe theory is a mixed bag. The small-area side is solid: the DV calculation reproduces the known tail without fitting, and the exact mapping to the Ferrari-Spohn model with a parabolic wall (Appendix C) is a neat, clean result. The subcritical Gaussian position distribution and the variance formula (18) check out against simulations.\n\nThe problem is in Appendix D, in the claim of two third-order transitions for t ≠ T/2. Eq. D5, which is supposed to fix τ via the area constraint, is algebraically inconsistent. At the stated Xc2 = 3A/(2t), the area constraint (integrating D4) gives τ/t = 2, while the printed D5 gives τ/t ≈ 2.054. Using the printed τ, the trajectory in D4 has area ≈ 1.052 A, not A. That means the intermediate-branch action difference D7 is evaluated on a trajectory that does not satisfy the conditioning constraint. And if you use the correct τ = 2t, D7 and D8 do not agree at Xc2—they differ by 5A^2/(D0 T^3). So the second third-order transition, as printed, is not established. The t = T/2 result in the main text is safe because the two critical lines merge there, so this is a localized but load-bearing error in the general case.\n\nMinor issues: the experimental histogram lacks error bars; the large-area conditioned window contains only 22 trajectories; and the phase transitions are tested only in simulation. None of these undercut the central experimental claim.\n\nVerdict: send it to a competent referee. The first measurement of the Airy distribution is a genuine step forward, and the DV/FS mapping is worth publishing. But the authors need to fix Appendix D and re-derive the intermediate branch before the two-transition statement can stand. I'd cite the experimental result, not the two-transition claim, until it is corrected.","headline":"First direct measurement of the Airy distribution is credible, but an algebraic error in Appendix D undermines the second third-order phase transition as printed.","tokens_in":17566,"tokens_out":19235,"would_cite":true,"duration_ms":156480,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60J65","60F10","82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"The Airy distribution—the probability law of the area under a Brownian excursion—is measured directly in a colloidal experiment, and its large-deviation function is shown to contain two third-order dynamical phase transitions.","keywords":["Airy distribution","Brownian excursion","large deviation theory","optimal fluctuation method","dynamical phase transition","conditional position distribution","colloidal experiment","area statistics"],"falsifier":"A concrete test: simulate large numbers of Brownian excursions conditioned on area $\\tilde A$, record the position at $t=T/2$, and compare $-\\ln \\tilde p_{\\tilde A}(z,1/2)$ with the predicted $g(z)$ from the paper (quadratic for $z<3$, cubic for $z>3$). The transition is sharp only as $\\tilde A\\to\\infty$, so the test is whether the third-derivative jump appears and its location tends to $z=3$ as $\\tilde A$ grows; a crossover whose location drifts with $\\tilde A$ or that never sharpens would falsify the transition claim. A similar check applies to the asymmetric critical lines $X_{c1}(t)$ and $X_{c2}(t)$.","tokens_in":16513,"feed_emoji":"🔬","tokens_out":6970,"duration_ms":65142,"temperature":0.7,"pith_summary":"This paper reports the first direct experimental measurement of the Airy distribution: the probability law of the area swept out by a Brownian excursion, a random path that starts and ends at zero and stays positive in between. The experiment tracks dilute colloidal particles, converts their trajectories into excursions, and records the histogram of areas, which matches the theoretical scaling curve. The paper also derives the position distribution of an excursion conditioned on having a given area, in both the small-area and large-area limits. For large areas this conditional distribution develops two singularities in its large-deviation function, interpreted as third-order dynamical phase transitions; for small areas it coincides with the distribution of a Brownian bridge pushed away from a moving wall, and the paper explains that coincidence through an exact mapping.","feed_headline":"First lab measurement of the Airy distribution of Brownian excursions","feed_subtitle":"Colloidal particle tracks confirm the area law and expose two sharp phase transitions in conditioned paths.","key_machinery":"The machinery has four main pieces. (1) The scaling form $f(\\xi)$ with its Laplace transform from Brownian excursion theory, giving the exact area distribution used for comparison. (2) The Wiener action $s[x(t)]=\\frac12\\int_0^T \\dot x^2(t)\\,dt$, whose minimization under the area constraint yields the optimal large-area trajectory $x_A^*(t)=(6At/T^2)(1-t/T)$. (3) The tilted-generator eigenvalue problem of the long-time large-deviation method, whose ground state is an Airy function and gives the small-area tail and the stationary conditional position distribution. (4) The tangent construction of the calculus of one-sided variations, which produces the correct constrained optimal trajectories when the unconstrained parabola would cross the origin; the switch between solution branches at $X_{c1}(t)$ and $X_{c2}(t)$ is what produces the third-order transitions. A biased-ensemble transformation connects the conditioned-excursion problem to the moving-wall problem and explains the small-area coincidence.","core_discovery":"The central claim is that the Airy distribution, previously a largely mathematical object, can be observed directly in a simple tabletop experiment and that its conditioned statistics conceal a nontrivial large-deviation structure. Using dilute colloidal suspensions and single-particle tracking, the authors construct Brownian excursions and measure the distribution of the area $A=\\int_0^T x(t)\\,dt$, finding agreement with the analytic scaling function $f(\\xi)$ in $P(A,T)=D_0^{-1/2}T^{-3/2}f(A/\\sqrt{D_0T^3})$. They then compute the conditional single-time position distribution of an excursion with prescribed area $A$. In the small-area limit this distribution is stationary in time and equals the wall-constrained bridge distribution; in the large-area limit it is non-stationary and its large-deviation function contains two third-order singularities, at critical position values $X_{c1}(t)$ and $X_{c2}(t)$, corresponding to the onset of zero-contact intervals in the optimal trajectory. The paper further shows that the two tails of the Airy distribution are governed by two different large-deviation techniques, and it derives the corresponding tails for the area under the square of an excursion.","pith_inferences":["If the third-order transitions survive finite-size corrections, the same kink structure should appear in other observables built from time-integrated powers of the excursion, not only the area and its square.","The correspondence between conditioned and biased ensembles suggests that the Airy-function ground state may be a universal signature of any one-dimensional process whose tilted generator has a linear potential, so the same distribution could appear in confined or active systems.","Because the measured histogram is slightly broadened near its maximum by particle polydispersity, higher-resolution measurements with more monodisperse particles should reduce the residual discrepancy and sharpen the observed tails."],"forward_implications":["The Airy distribution acquires direct experimental confirmation, grounding its many applications in inventory theory, data-storage algorithms, graph theory, and interface fluctuations.","The explicit formulas for the conditional position distribution give simulation and future single-particle experiments a precise benchmark to test against.","The predicted third-order dynamical phase transitions can be sought experimentally or numerically; they should sharpen as the dimensionless area $\\tilde A$ grows.","The exact mapping to the moving-wall problem transfers results between the two settings and indicates a shared universality class.","The derived tails for the distribution of the area under the square of an excursion extend the same two large-deviation methods to a new observable."],"supporting_citations":[{"why":"supplies the original probabilistic derivation of the Laplace transform of the Airy distribution.","marker":"[7]"},{"why":"independently provides the Laplace transform via excursion theory, establishing the object under study.","marker":"[8]"},{"why":"gives the closed analytic form of the scaling function used for comparison with the measured histogram.","marker":"[10]"},{"why":"provides the path-integral derivation of the Laplace transform, linking the Airy distribution to the action formalism used here.","marker":"[13]"},{"why":"supplies the tilted-generator eigenvalue method that yields the small-area tail and the stationary conditional position distribution.","marker":"[26–29]"},{"why":"provides the tangent construction of the calculus of one-sided variations that generates the supercritical optimal trajectories.","marker":"[50]"},{"why":"establishes the equivalence between biased and conditioned ensembles, used in the mapping to the wall-constrained distribution.","marker":"[45]"},{"why":"provides the wall-constrained bridge distribution that the small-area conditional distribution reproduces.","marker":"[44]"},{"why":"supplies the transformation used to construct Brownian excursions from experimental and simulated bridges.","marker":"[56]"},{"why":"supplies pre-exponential corrections to the tails and the normalization needed for quantitative comparison.","marker":"[21]"}],"fun_headline_variants":["First direct measurement of Airy distribution in lab","Colloid tracks confirm Airy area law, reveal phase jumps","Brownian excursion area distribution measured for first time","Lab test of Brownian area uncovers third-order transitions","Airy distribution from colloidal motion, with large-deviation twist"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument's load-bearing premise is that when the unconstrained optimal trajectory would cross the origin, the true constrained minimizer is exactly zero over a finite interval; if this tangent-construction ansatz is not exact, the two predicted third-order phase transitions could be artifacts of the variational method rather than real properties of the conditional distribution.","fun_headline_variants_meta":{"raw":{"variants":["First direct measurement of Airy distribution in lab","Colloid tracks confirm Airy area law, reveal phase jumps","Brownian excursion area distribution measured for first time","Lab test of Brownian area uncovers third-order transitions","Airy distribution from colloidal motion, with large-deviation twist"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000375,"raw_usage":{"total_tokens":1997,"prompt_tokens":938,"completion_tokens":1059,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":980}},"tokens_in":554,"tokens_out":1059,"duration_ms":9529,"temperature":1.0,"reasoning_tokens":980,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:41:22.731859+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: simulate large numbers of Brownian excursions conditioned on area $\\tilde A$, record the position at $t=T/2$, and compare $-\\ln \\tilde p_{\\tilde A}(z,1/2)$ with the predicted $g(z)$ from the paper (quadratic for $z<3$, cubic for $z>3$). The transition is sharp only as $\\tilde A\\to\\infty$, so the test is whether the third-derivative jump appears and its location tends to $z=3$ as $\\tilde A$ grows; a crossover whose location drifts with $\\tilde A$ or that never sharpens would falsify the transition claim. A similar check applies to the asymmetric critical lines $X_{c1}(t)$ and $X_{c2}(t)$.","supporting_citations":[{"cited_title":"von Smoluchowski, Annal","cited_arxiv_id":null,"evidence_quote":"supplies the original probabilistic derivation of the Laplace transform of the Airy distribution."},{"cited_title":"Langevin, Comp","cited_arxiv_id":null,"evidence_quote":"independently provides the Laplace transform via excursion theory, establishing the object under study."},{"cited_title":"Tak´ acs, Adv","cited_arxiv_id":null,"evidence_quote":"gives the closed analytic form of the scaling function used for comparison with the measured histogram."},{"cited_title":"Louchard, J","cited_arxiv_id":null,"evidence_quote":"provides the path-integral derivation of the Laplace transform, linking the Airy distribution to the action formalism used here."},{"cited_title":"Chetrite and H","cited_arxiv_id":null,"evidence_quote":"provides the tangent construction of the calculus of one-sided variations that generates the supercritical optimal trajectories."},{"cited_title":"Janson and G","cited_arxiv_id":null,"evidence_quote":"establishes the equivalence between biased and conditioned ensembles, used in the mapping to the wall-constrained distribution."},{"cited_title":"This unexpected coincidence can be ex- plained by an exact mapping that we found between the two systems","cited_arxiv_id":null,"evidence_quote":"provides the wall-constrained bridge distribution that the small-area conditional distribution reproduces."},{"cited_title":"Mazzolo, J","cited_arxiv_id":null,"evidence_quote":"supplies the transformation used to construct Brownian excursions from experimental and simulated bridges."},{"cited_title":"Medalion, E","cited_arxiv_id":null,"evidence_quote":"supplies pre-exponential corrections to the tails and the normalization needed for quantitative comparison."}],"review_version":1}