{"id":"c9447b89-19e1-43db-bb89-44302f475197","arxiv_id":"1908.03726","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A numerical Riemann-Hilbert solution maps the γ-generalized two-body log-gas onto a Muttalib-Borodin ensemble with an effective potential, interpolating the hard-edge exponent from -1/3 at γ=1 to -1/2 at γ=0.","lead":"This paper introduces a random matrix model with an extra two-body interaction of tunable strength γ and computes its eigenvalue density by solving a Riemann-Hilbert problem numerically. The authors show the density can be recast as a standard Muttalib-Borodin ensemble with a γ-dependent effective potential, with the hard-edge divergence exponent crossing from -1/3 to -1/2 as γ decreases from 1 to 0.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"One-cut support and the standard RH reduction are assumed for γ<1 without proof; a direct Euler-Lagrange consistency check is needed before the effective-potential mapping can be accepted.","rationale":"I read the paper as proposing a numerical Riemann-Hilbert method for the γ-ensemble and claiming that the effect of γ can be represented by a γ-dependent effective single-particle potential. For that equivalence to hold, the equilibrium measure must be one-cut on [0,b] with the hard edge at 0, and the CR contour construction must remain valid when γ is reduced from 1. The reader's weakest assumption identifies exactly this one-cut condition. My concern is a sharper version: the entire scalar reduction leading to equation (25) assumes not only one-cut support but also that the non-standard jump relation (12), which contains the unknown measure μ([x,b]), can be handled by the same fixed Joukowsky transform as in the γ=1 case. If the support geometry changes, f is not well-defined and equation (25) describes a different problem. I did not find an internal algebraic contradiction in the derivation, and the match with the known γ=1 and γ=0 endpoints is genuine supporting evidence. However, the numerical method is under-specified and no independent verification against the original Euler-Lagrange equation is shown, so the central claim currently rests on an unverified assumption rather than on proof. This does not warrant rejection because the construction is plausible, self-consistent, and reproduces known limits; it does warrant keeping the CONDITIONAL verdict and requesting the direct consistency test described above.","tokens_in":16,"tokens_out":16139,"duration_ms":244958,"concrete_test":"Directly test the ansatz by solving the original equilibrium problem for θ=2, V=2x: discretize probability measures on [0,L] and minimize the energy functional (5), or solve the Euler-Lagrange equation (6) for the density σ, at γ=0.1, 0.2, 0.4, 0.6, 0.8, and 1.0. Then check (i) the recovered support is one interval [0,b] with no gaps, and (ii) the density obtained from equations (25)-(29) with the same γ satisfies the original Euler-Lagrange equation on the support to within a small tolerance, e.g., residual below 1% of V, and has total mass 1. If either check fails at any γ, the effective-potential mapping and the exponent curve in Figure 7 are not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, that the γ-ensemble is equivalent to an MB ensemble with an effective potential Veff(x;γ), rests on the unproven assumption that the equilibrium measure of (4) is supported on a single interval [0,b] for every 0<γ≤1. Section 3 carries the whole γ<1 construction over from Claeys-Romano: equation (12) generalizes the jump condition with the unknown measure μ([x,b]) on the right-hand side, and the derivation of (16), (19), and the scalar integral equation (25) uses the fixed Joukowsky contour ν and the two-preimage relation Jc(s1)=Jc(s2)=x. This presupposes the one-cut hard-edge geometry. If the support splits or develops an interior gap for some γ, the contour ν no longer maps D to Hθ\\[0,b], the function f=N_++N_- is not single-valued along ν, and equation (25) solves a different model from the original γ-ensemble. No proof and no numerical check of the one-cut condition for γ<1 is provided; matching the endpoints γ=1 and γ=0 is necessary but not sufficient evidence. Because the headline claim is exactly an equivalence of equilibrium measures, this gap is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a log-gas model (4), the 'γ-ensemble', whose joint density contains the additional two-body factor |x_i^θ - x_j^θ|^γ on the non-negative half-line. The authors adapt the Claeys-Romano (CR) Riemann-Hilbert method for the Muttalib-Borodin ensemble to derive an integral equation (25) for an auxiliary function f, define an effective potential Veff by Veff' = f/x, and claim that the γ-ensemble is equivalent at the level of the global equilibrium density to an MB ensemble with potential Veff. For V(x)=2x and θ=2, they compute the hard-edge density and report a crossover of the edge exponent from -1/3 at γ=1 to -1/2 at γ=0 (Figure 7). A second example with r(x)=e^x and V=x^2/2 is treated in Section 6 and the Appendix. The main difficulties are that the one-cut support assumption is unproven for γ<1, the derivation of the integral equation contains a visible typo in Eq. (24), the numerical solution is undocumented, and the density formula (29) contains a factor 1/γ that appears inconsistent with the claimed MB mapping.","tokens_in":11204,"tokens_out":33946,"duration_ms":299207,"significance":"If correct, the effective-potential mapping is a useful simplification: it reduces the two-body interaction with strength γ to a one-body confining potential, and the hard-edge exponent crossover is a concrete quantitative prediction that can be tested by independent numerical simulation of the log-gas (4). The paper builds on a rigorous framework (CR) and treats both hard-edge and soft-edge examples, with the γ=0 and γ=1 limits matching known ensembles, which is a good sanity check. However, the central equivalence is not yet established: the one-cut assumption is load-bearing and unchecked, the derivation has gaps and typos, and the numerical evidence lacks reproducibility details. The approach is promising, but the manuscript needs substantial revision before the claims can be accepted.","major_comments":[{"comment":"The entire construction for 0<γ<1 presupposes that the equilibrium measure of (4) is supported on a single interval [0,b] with a hard edge at 0 and no interior gaps, but this one-cut condition is only stated for γ=1 in Section 2 and is neither proved nor numerically checked for γ<1. The Joukowsky transform (10), the contour ν, and the two-preimage relation Jc(s1)=Jc(s2)=x all depend on this geometry; if the support develops a gap or splits as γ varies, equation (25) solves a different model from the original γ-ensemble. I do not object to the generalized jump condition (12), which does follow from (13) and the Euler-Lagrange equation (6), but the unproven one-cut assumption is load-bearing for the claimed equivalence and needs a proof or at least a numerical test that the computed density satisfies (6) on [0,b] and the variational inequality outside.","section":"Section 3"},{"comment":"Equation (29) contains an explicit factor 1/γ in front of the integral, while the CR formula for the MB ensemble (γ=1) has no such factor. As written, the γ-ensemble density is (1/γ) times the MB density with effective potential Veff, so the statement that 'the density for γ<1 has the same expression as that for γ=1, except that V is replaced by Veff' is not literally correct; moreover, if both densities are normalized to 1, such a constant factor cannot be right. Please reconcile the factor with the effective-potential mapping and show that the raw output of (29) integrates to 1 without an ad hoc normalization.","section":"Section 4, Eq. (29)"},{"comment":"The two formulas displayed for N+(s1) and N-(s2) in Eq. (24) are identical, which must be a typo; the kernel φ in Eq. (26) indicates that the second line should involve \\bar h(y) in the denominators. As printed, the derivation of the integral equation (25) cannot be followed. Please correct Eq. (24) and include the intermediate steps that lead to (25), including the contour orientations and the branch choices for h and \\bar h.","section":"Section 3, Eq. (24)"},{"comment":"The numerical solution is not documented: the paper states that (25) and (20) are solved self-consistently but gives no discretization scheme, quadrature rule, iteration procedure, convergence criterion, or error analysis. Figure 7 reports 'uncertainties in the numerical estimates' but the source of these uncertainties is not explained. Without this information, the central numerical prediction (the exponent crossover) is not reproducible and the quoted exponents cannot be assessed.","section":"Section 5"},{"comment":"The soft-edge integral equation (39) differs from the hard-edge equation (25) by the absence of the constant term, and the density formula (44) has no 1/γ factor whereas (29) does. If both cases are meant to follow from the same effective-potential construction, the structural difference should be explained; otherwise the presence/absence of the 1/γ factor looks like an inconsistency rather than a feature of the different boundary conditions.","section":"Appendix, Eq. (39)"}],"minor_comments":[{"comment":"The y-axis label 'iy(×10^{10})' is unclear; the text says the tiny imaginary parts are numerical noise. Please relabel the axis as Im(Jc(s)) and state whether the plotted curve is the raw numerical map or a smoothed curve.","section":"Figures 2 and 11"},{"comment":"The caption says 'Normalized density'; please specify how the normalization is performed and whether Eq. (29) was already normalized or was rescaled for display.","section":"Figure 6 caption"},{"comment":"The inverse map h is double-valued; the paper uses h(x) and \\bar h(x) without defining which preimage is h(x) and without stating that \\bar h is the complex conjugate branch. Please define the branches precisely.","section":"Equations (23)-(26)"},{"comment":"The abstract says the density is computed by 'numerically solving the Riemann-Hilbert problem'; in fact the paper solves a scalar integral equation obtained after assuming an RH reduction. Please rephrase to match the actual method.","section":"Abstract and Section 1"},{"comment":"The Euler-Lagrange equation is written with V(x) on the left and the constant l on the right; in Eq. (9a) the constant appears as V(x)+l. Please unify the sign convention.","section":"Equation (6)"}],"recommendation":"major_revision","confidential_remarks":"The main mathematical risk is the one-cut support assumption; if the journal can accept a physics-level treatment, a numerical consistency check against direct sampling of (4) would go a long way. I would suggest inviting a referee familiar with Claeys-Romano's RH reduction to evaluate the corrected equations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fair warning before you read: this is a numerical extension of Claeys–Romano to a γ-parameter family of log-gases, and the headline claim—that the γ-ensemble is an MB ensemble with an effective potential—is plausible but rests on an unproven one-cut assumption. I don't think it's wrong, and it deserves a real referee, but it needs work before I'd trust the central equivalence.\n\nWhat's genuinely new: the model in Eq. (4) is a natural generalization of Muttalib–Borodin with a tunable exponent γ on the second interaction, motivated by the quasi-1D to 3D crossover in disordered conductors. The adaptation of the CR Riemann–Hilbert machinery to 0<γ<1 is real, and the numerical method gives a concrete way to compute densities for this family. The endpoint checks are honest: γ=1 recovers MB, and the hard-edge exponent goes from −1/3 to −1/2 as γ→0, matching the known Laguerre limit. That's a nontrivial sanity check, and I believe it. The extension to r(x)=e^x is a nice bonus that shows the technique isn't tied to the x^θ interaction. Citations look appropriate: CR, CW, Borodin, Forrester, and the Muttalib–Gopar–Wölfle line are the relevant anchors.\n\nSoft spots, in order of severity. First, the one-cut assumption. The whole construction uses a single contour ν and a two-preimage relation on [0,b]. For γ<1 there is no proof—and as far as I can tell no numerical check—that the equilibrium measure stays supported on one interval. If the support splits, the integral equation (25) solves a different model from the original γ-ensemble. This is load-bearing, because the effective-potential equivalence is exactly an equality of equilibrium measures. Second, equation (12) is asserted rather than derived; it's the generalized jump condition that drives everything. A direct Euler–Lagrange consistency check—plug the computed density back into (6)—would catch a wrong jump condition, and I'd want to see that in the paper. Third, the numerics are under-described: no discretization scheme, no convergence test, no error analysis, and no code. Figure 7 has error bars, but I can't reproduce the method from the text as it stands. Finally, the abstract and summary overclaim: \"any two-body interaction\" with a known Joukowsky transform is too broad, and the effective-potential picture is more a reformulation than a predictive mapping. That's not circularity, but it's not an independent check either.\n\nWho gets value: people working on log-gas/RH methods and mesoscopic transport toy models. This is a serious paper, not a desk reject. I'd send it to peer review, but with a strong request for a one-cut proof or numerical evidence, a derivation or consistency verification of (12), and proper numerical details.","headline":"A plausible numerical extension of the Claeys–Romano method to a γ-deformed Muttalib–Borodin ensemble, but the effective-potential equivalence rests on an unproven one-cut assumption and needs a consistency check before I'd fully trust it.","tokens_in":11818,"tokens_out":3034,"would_cite":false,"duration_ms":33709,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60B20","82B44"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper shows that a log-gas with an added two-body interaction is equivalent, for its equilibrium density, to a standard ensemble with a gamma-dependent effective potential.","keywords":["gamma-ensemble","log-gas model","random matrix ensemble","Riemann-Hilbert problem","Joukowsky transform","effective confining potential","hard-edge density exponent","disordered conductors"],"falsifier":"Compute the equilibrium density for $\\theta=2$, $V=2x$, and a representative $\\gamma$ such as $0.5$ by an independent numerical method (direct minimization of the energy functional or Monte Carlo sampling of the log-gas) and compare the support and density profile with the integral-equation solution. If the support has more than one component, or the near-zero density exponent differs from the computed value, the effective-potential claim fails for that $\\gamma$.","tokens_in":22,"feed_emoji":"🧮","tokens_out":9807,"duration_ms":121966,"temperature":0.7,"pith_summary":"The paper introduces a log-gas model, the $\\gamma$-ensemble, where a parameter $\\gamma$ controls the strength of a second two-body interaction on top of the usual logarithmic repulsion. Its central claim is that, despite this extra interaction, the equilibrium density is exactly that of the known $\\gamma=1$ ensemble, provided the confining potential is replaced by a $\\gamma$-dependent effective potential. The effective potential is obtained by solving a self-consistent linear integral equation, and the authors implement this numerically for a linear confining potential and for an exponential interaction. The main quantitative result is that the power-law divergence of the density at the hard edge changes continuously from $-1/3$ at $\\gamma=1$ to $-1/2$ at $\\gamma=0$. A reader would care because this provides a solvable interpolation between two known random-matrix limits and a concrete tool for interaction-strength crossovers.","feed_headline":"Hard-edge exponent slides from -1/3 to -1/2 as gamma weakens","feed_subtitle":"A gamma-dependent effective potential reduces the generalized random-matrix ensemble to a known solvable form.","key_machinery":"The central objects are the vector-valued Riemann-Hilbert problem for the transforms $g$ and $\\tilde g$ and the Joukowsky transform $J_c(s)=c(s+1)(s+1/s)^{1/\\theta}$, which maps a contour $\\nu$ in the $s$-plane onto the cut plane of the density. The technical step is to encode the $\\gamma$ modification in a function $f(x;\\gamma)$, defined so that $V'_{\\mathrm{eff}}=f/x$, and to convert the Riemann-Hilbert data into the linear integral equation (25) for $f$. This $f$ fixes the endpoint parameter $c$ through Eq. (20) and the density through Eq. (29), reducing the $\\gamma$-ensemble to an effective one-body problem.","core_discovery":"For $0<\\gamma\\le 1$, the joint density of Eq. (4) has an equilibrium measure whose support is taken to be a single interval $[0,b]$ with a hard edge. Within this one-cut regime, the vector-valued Riemann-Hilbert problem of the $\\gamma=1$ case still governs the problem, with the $\\gamma$ dependence entering through a function $f(x;\\gamma)$. Defining the effective potential by $V'_{\\mathrm{eff}}(x;\\gamma)=f(x;\\gamma)/x$, the function $f$ obeys the linear integral equation (25), whose kernel is built from the inverse of the Joukowsky transform. Solving this equation self-consistently with the endpoint condition (20) yields the equilibrium density through formula (29). For the worked example $V=2x$, $\\theta=2$, the computed density shows the hard-edge exponent crossing continuously from $-1/3$ at $\\gamma=1$ to $-1/2$ at $\\gamma=0$; the same scheme, with the appropriate Joukowsky transform, is applied to an exponential interaction on the whole real line, where the density has two soft edges and stays non-diverging while the effective potential varies with $\\gamma$.","pith_inferences":["If the effective-potential mapping extends beyond the one-point density to correlation kernels, the $\\gamma$-ensemble would inherit the integrable structure of the $\\gamma=1$ ensemble; this is a testable extension the paper leaves open.","The one-cut assumption may fail for some $\\gamma$ near $0$ or for other confining potentials, producing a multi-cut support; a numerical scan of the support as a function of $\\gamma$ would delimit where the construction is valid.","The same integral-equation method could be applied to other monotone functions $r(x)$ to test whether the $\\gamma=0$ hard-edge exponent is always $-1/2$ or depends on the form of the interaction.","The crossover curve of the exponent in Fig. 7 invites a scaling-law fit, which could expose a universal interpolation between the two known limits."],"forward_implications":["For any $\\gamma$ in $(0,1]$, the equilibrium density of the $\\gamma$-ensemble can be computed by solving a single linear integral equation, once the inverse Joukowsky map is known.","The hard-edge exponent is not fixed by the exponent $\\theta$ alone; it varies continuously with the interaction strength $\\gamma$, interpolating between the $\\gamma=1$ and $\\gamma=0$ limits.","The effective-potential picture means existing techniques for the $\\gamma=1$ ensemble, including correlation kernels, can in principle be adapted to the $\\gamma$-ensemble.","The method is not specific to power-law interactions: the same construction works for an exponential interaction $r(x)=e^x$, giving non-diverging densities with $\\gamma$-dependent effective potentials.","The $\\gamma=0$ limit connects to the known orthogonal-Laguerre behavior with hard-edge exponent $-1/2$."],"supporting_citations":[{"why":"introduced the original one-parameter log-gas with an added power-law interaction that this paper generalizes.","marker":"[5]"},{"why":"provided the exact solution and detailed analysis of the gamma=1 ensemble that serves as the baseline.","marker":"[6]"},{"why":"supplies the potential-theoretic existence and uniqueness of the equilibrium measure used in the derivation.","marker":"[16]"},{"why":"supplied the Riemann-Hilbert method, the Joukowsky transform, and the integral-equation structure that the paper adapts for gamma<1.","marker":"[17]"},{"why":"derived the Joukowsky transform for the exponential-interaction model used as the second example.","marker":"[18]"}],"fun_headline_variants":["Gamma tunes hard-edge exponent from -1/3 to -1/2","Effective potential explains gamma-dependent hard edge","Hard-edge exponent interpolates as gamma varies","Generalized RMT: edge exponent shifts with gamma","Gamma-dependent effective potential shifts hard-edge exponent"],"cache_read_input_tokens":13952,"weakest_assumption_plain":"The construction assumes the equilibrium measure is supported on a single interval $[0,b]$ for every $0<\\gamma\\le 1$; if the support splits into multiple intervals, the single-contour Joukowsky transform and the integral equation no longer apply.","fun_headline_variants_meta":{"raw":{"variants":["Gamma tunes hard-edge exponent from -1/3 to -1/2","Effective potential explains gamma-dependent hard edge","Hard-edge exponent interpolates as gamma varies","Generalized RMT: edge exponent shifts with gamma","Gamma-dependent effective potential shifts hard-edge exponent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1641,"prompt_tokens":849,"completion_tokens":792,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":465,"completion_tokens_details":{"reasoning_tokens":718}},"tokens_in":465,"tokens_out":792,"duration_ms":7513,"temperature":1.0,"reasoning_tokens":718,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:03:55.764847+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the equilibrium density for $\\theta=2$, $V=2x$, and a representative $\\gamma$ such as $0.5$ by an independent numerical method (direct minimization of the energy functional or Monte Carlo sampling of the log-gas) and compare the support and density profile with the integral-equation solution. If the support has more than one component, or the near-zero density exponent differs from the computed value, the effective-potential claim fails for that $\\gamma$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduced the original one-parameter log-gas with an added power-law interaction that this paper generalizes."},{"cited_title":"B 536 704–732 ISSN 0550-3213 URL http://dx.doi.org/10","cited_arxiv_id":null,"evidence_quote":"provided the exact solution and detailed analysis of the gamma=1 ensemble that serves as the baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the potential-theoretic existence and uniqueness of the equilibrium measure used in the derivation."},{"cited_title":"org/10.1088/0951-7715/27/10/2419","cited_arxiv_id":null,"evidence_quote":"supplied the Riemann-Hilbert method, the Joukowsky transform, and the integral-equation structure that the paper adapts for gamma<1."}],"review_version":1}