{"id":"2cec077d-82a3-49f0-9421-380539404455","arxiv_id":"1908.05320","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A next-to-leading-order QCD calculation with kT broadening reproduces the LHCb bbar to J/psi J/psi correlation data for all measured observables and momentum cuts.","lead":"This paper compares LHCb measurements of bottom quark pair correlations, seen through J/psi pairs, to a next-to-leading-order QCD calculation with transverse momentum broadening. The calculation matches the data and shows that the J/psi decay products are insensitive to the broadening that affects the parent bottom quark pairs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's own Sec. III caveat that the Gaussian kT kick is not exactly equivalent between initial- and final-state implementations at NLO is the load-bearing assumption; if the pair-only kick distorts 2-to-3 kinematics, the claimed description of the LHCb |Delta phi| and pTp shapes could change.","rationale":"We read the central claim as a phenomenological demonstration: the exclusive HVQMNR NLO framework with a single Gaussian kT smearing and Peterson fragmentation can describe the normalized LHCb pair-shape data, and the parent bbar distributions are more responsive to kT than the accessible J/psi pair distributions. The reader's weakest assumption points to the kT implementation. This is well placed: the paper itself flags the LO/NLO inequivalence in Sec. III, and the affected observables (small-|Delta phi| tail, pTp shoulder, AT shape) are precisely the ones the model is said to reproduce. We do not see an internal inconsistency or a violation of conservation laws that would make the current implementation invalid, but the ambiguity is real. Our proposed check is deliberately direct: rerun the same code with the kick on the full final state; if the pair observables stay within the theoretical uncertainty band, the claim stands as stated; if they move beyond it, the agreement with data is partly an artifact of the chosen kick scheme. We did not choose the B-to-J/psi decay-model specification as the single concern because, although it is under-described, the J/psi pair distributions are directly compared to data and are less kT-sensitive; the kT implementation affects both the bbar and, indirectly, the interpreted sensitivity statement. The lack of a quantitative goodness-of-fit is a presentation issue, but this stress-test focuses on the argument's load-bearing point; the reader already made the CONDITIONAL verdict based on it. Hence no verdict change.","tokens_in":23805,"tokens_out":8617,"duration_ms":89989,"concrete_test":"Recompute Figs. 1-6 with the Gaussian kT of Eq. (3) applied to the full final state (including the light parton in 2-to-3 events) rather than boosting only the QQ pair; equivalently, apply it to the initial-state partons. Compare the normalized |Delta phi|, pTp, AT, and M shapes for pT>2 and 7 GeV with the curves shown and with the mass/scale bands of Figs. 8-9. If the bbar distributions shift by more than the uncertainty band or the LHCb error bars, the NLO equivalence caveat is load-bearing. Also verify that the relative kT sensitivity of bbar vs J/psi pairs (default vs <kT^2>=0) is unchanged.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that NLO HVQMNR plus Gaussian kT broadening and Peterson fragmentation describes all LHCb bbar to J/psi J/psi observables, and that parent bbar pairs are more kT-sensitive than J/psi pairs. The load-bearing assumption is the treatment of the kT kick. In Sec. III, the kick is applied to the QQ pair in its rest frame via Eq. (3), and the text explicitly notes that 'the two methods of introducing kT are equivalent if the calculation is LO but at NLO a light parton in the final state can make the correspondence inexact.' This is not a minor technicality: the NLO 2-to-3 configurations with a hard light parton are precisely what populate the small-|Delta phi| region and the shoulder in pTp that LHCb observes. If the kT should act on the full final state (or initial state) rather than only on the pair, the bbar |Delta phi|, pTp, and AT shapes could shift significantly. Because the same kT treatment underlies the comparison in Sec. V of <kT^2>=0 vs default, the conclusion that parent bbar pairs are more sensitive than J/psi pairs is also entangled with this implementation choice. The paper's agreement with data is asserted visually, so this caveat is not bracketed by a quantitative uncertainty; it is the weakest point in the argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents next-to-leading-order calculations of b-bbar pair production in the exclusive HVQMNR framework, with Peterson fragmentation and Gaussian intrinsic-kT broadening, and compares the resulting b-meson-pair and J/psi-pair observables with LHCb measurements from p+p collisions at 7 and 8 TeV for four minimum pT thresholds. The paper reports good agreement for all six pair observables, studies the sensitivity of the results to kT broadening, evaluates mass and scale uncertainties, examines the rapidity dependence, and extends the model to p+Pb and Pb+Pb collisions by increasing the broadening and modifying the fragmentation parameter in order to illustrate possible cold-nuclear-matter effects. The main physical conclusions are that the NLO-plus-kT model can describe the LHCb correlated observables and that parent b-bbar pairs are more sensitive to kT broadening than the J/psi pairs produced in their decays.","tokens_in":24156,"tokens_out":10923,"duration_ms":113418,"significance":"If the central claim holds, the paper provides a useful baseline: a fixed-parameter NLO calculation tuned to single-particle data reproduces correlated b-bbar observables, and the comparison is a genuine prediction test rather than a fit, since the parameters mb, muF/m, muR/m, epsilon_P, and <kT^2> were set in Ref. [14] from total bbar cross sections, FONLL B-meson pT, and Upsilon pT data. The paper covers six observables and four pT cuts, which is a broad and valuable test, and the explanation that the isotropic B decay decorrelates the J/psi directions is a clean, nontrivial insight. The nuclear-matter section is explicitly illustrative rather than a quantitative prediction. The main weaknesses are the unquantified NLO kT-kick scheme, the mismatched selection used for the 'bb' comparisons, and the absence of any statistical measure of the claimed agreement; these need to be addressed before the central claim can be considered established.","major_comments":[{"comment":"The calculated 'bb' distributions are obtained by applying the minimum pT cut to the parent B mesons, while the LHCb 'bb' data are selected by the daughter J/psi pT thresholds and the parent B pT is not measured. Because the B to J/psi decay is not collinear, a B meson with pT below the nominal threshold can produce a J/psi above it (the J/psi momentum in the B rest frame is about 1.7 GeV), and conversely a B meson above the threshold can produce a J/psi below it. The two phase-space regions are therefore different, and the comparison of |Delta phi*|, pTp, and AT in Figs. 1-6 may be biased. Please quantify this effect by repeating the 'bb' calculation with the actual acceptance, i.e. generating B pairs, decaying them, applying only the J/psi pT cuts, and constructing the vertex-derived direction variables, or otherwise justify that the B pT cut is a faithful proxy for the LHCb selection.","section":"Sec. IV, Figs. 1-6"},{"comment":"The kT kick is applied to the QQ pair in its rest frame, and the text explicitly states that this procedure is equivalent to an initial-state kick only at leading order, while at NLO a final-state light parton makes the correspondence inexact. This is the central approximation behind all the paper's results, including the comparison with <kT^2>=0 in Sec. V. Since the 2-to-3 configurations with a hard light parton are precisely what populate the small-|Delta phi| region and the shoulder in pTp, the claimed description of these shapes depends on an unquantified scheme choice. I request an estimate of the effect obtained by applying the kick to the full NLO final state instead of only to the pair, or a phase-space argument establishing that the difference is negligible for the observables under study.","section":"Sec. III, Eq. (3)"},{"comment":"The paper's central claim that the calculations 'reproduce the data very well' and are in 'good agreement' is supported only by visual inspection. No chi-square, pull, or other goodness-of-fit metric is reported for any observable at any pT threshold. Given that the LHCb data and uncertainties are public, a quantitative comparison (e.g., chi2 per degree of freedom computed with the correlated systematic uncertainties, or a table of pulls) is feasible and would turn the qualitative claim into a testable statement. If such a metric cannot be provided, the wording should be softened to state that the calculations are consistent with the data within the model uncertainties.","section":"Sec. IV, Figs. 1-6; Sec. IX"}],"minor_comments":[{"comment":"The displayed uncertainty bands are constructed from mass and scale variations only; they do not include the uncertainty in <kT^2> or the scheme ambiguity discussed in Sec. III. This should be stated explicitly in the captions of Figs. 8 and 9, since a reader may otherwise interpret the bands as the total theoretical uncertainty.","section":"Sec. VI, Eqs. (4)-(5)"},{"comment":"The nuclear-matter discussion uses only the central EPS09 set and models energy loss by changing epsilon_P to the e+e- value. The text is appropriately cautious, but the final paragraph of Sec. IX should state more explicitly that the discriminating power between broadening and energy loss is demonstrated only within this toy model, not as a quantitative prediction for data.","section":"Sec. VIII"},{"comment":"There are several typographical and grammatical errors that should be corrected: 'The bb pair distribution include' in Sec. IV; 'decrase' in Sec. VI; 'colliosions' in Ref. [4]; 'senstive' in Sec. IX; 'These contributions and summed together' in Sec. III; and 'The changes in the shadowing ratios is then small' in Sec. VIII.","section":"Throughout"},{"comment":"The caption uses 'J/psi' instead of the journal-style 'J/psi' used elsewhere; please make the notation consistent.","section":"Fig. 7"},{"comment":"For reproducibility, a table summarizing the central parameter values (mb, muF/m, muR/m, epsilon_P, Delta, and <kT^2> at the relevant energies) and their sources would be helpful; these values are currently scattered through the text and Ref. [14].","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The model parameters and much of the validation come from the author's own Ref. [14]; this is a legitimate reuse rather than circularity, because the LHCb pair observables studied here are new and were not used in fixing the parameters. The referee report focuses on three load-bearing issues: the mismatched parent-B selection for the 'bb' comparisons, the unquantified NLO kT-kick scheme, and the absence of quantitative goodness-of-fit measures. The nuclear-matter section is explicitly illustrative and should not be a deciding factor for publication. The paper is within the scope of a phenomenological QCD journal if the central p+p comparison is made quantitatively reliable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a follow-through application of Vogt's established HVQMNR+kt model to the LHCb bbar -> J/psi J/psi pair observables, and the interesting physics claim is that the parent bbar pairs are clearly more sensitive to the kT broadening than the J/psi pairs, because the B->J/psi decay randomizes the direction. That claim looks right and is useful: it tells the heavy-ion community to measure B-meson pair correlations rather than J/psi pairs if they want to isolate cold-nuclear-matter broadening. The paper deserves credit for being predictive: mb, scales, epsilon_p and <kT^2> were all fixed in Ref. [14], so the agreement with the LHCb pair shapes is not a fit to those data. The uncertainties are presented as envelopes from mass and scale variations, and the author is explicit that the nuclear section is illustrative, not definitive.\n\nThe soft spot is the one the author admits in Sec. III: the Gaussian kT kick is applied to the QQ pair in its rest frame, and the text says the initial-state and final-state implementations are only exactly equivalent at LO. The stress-test note is right that this is not a minor technicality. The NLO 2->3 configurations with a hard light parton are precisely what populate the small-|Delta phi| region and the shoulder in pTp, so if the kick should act on the full final state instead of only the pair, the shapes could shift. The paper does not quantify how much this approximation affects the comparison. That said, the caveat cuts both ways: the fact that the shapes come out close to the data with the pair-level kick is evidence that the approximation is not wildly wrong, at least for these observables. Still, without a goodness-of-fit metric, 'good agreement' is a visual statement.\n\nMinor points: no code or parameter files are shipped, so the comparison is not independently reproducible without reimplementing HVQMNR plus the kT treatment. The rapidity-dependence section is fine but not central. The nuclear modification discussion is explicitly illustrative and should not be read as a prediction.\n\nBottom line: this is a solid, honest phenomenological paper within an established program. It is not groundbreaking, but it gives a useful baseline and a clear experimental takeaway. I would send it to peer review and suggest the author add a quantitative comparison (chi2 or similar), address the kT equivalence caveat more quantitatively, and consider releasing the inputs. I would cite it if I were writing about heavy flavor correlations or cold nuclear matter.","headline":"A solid, honest NLO+kt comparison to LHCb bbar pair data; the useful result is that J/psi decay washes out kT sensitivity, but the pair-level kT kick remains an acknowledged approximation that deserves a quantitative check.","tokens_in":24712,"tokens_out":2739,"would_cite":true,"duration_ms":28440,"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":"One Gaussian transverse-momentum smearing reproduces all LHCb bottom-pair correlation measurements.","keywords":["bottom quark pair production","J/psi pair correlations","transverse momentum broadening","next-to-leading order QCD","HVQMNR","cold nuclear matter","LHCb forward rapidity","heavy flavor correlations"],"falsifier":"A direct measurement of $B$-hadron pair correlations with the same transverse-momentum cuts, without the $J/\\psi$ decay in between, would test the decay-decorrelation claim: if the parent $B$ pairs show no stronger dependence on broadening than the $J/\\psi$ pairs do, the central mechanism is wrong. Alternatively, a resummed or parton-shower-matched next-to-leading-order calculation with explicit all-order treatment of low pair transverse momentum should reproduce the same LHCb distributions; if it cannot match the data with the same single-kick prescription, the broadening model is inadequate.","tokens_in":23593,"feed_emoji":"⚛️","tokens_out":9365,"duration_ms":81806,"temperature":0.7,"pith_summary":"The paper argues that all six correlated observables of $b\\bar b$ pairs measured by LHCb through $B \\to J/\\psi$ decays can be reproduced by a next-to-leading-order QCD calculation once the pair is given a single Gaussian transverse-momentum kick. The kick, with an energy-dependent width, does the work that resummation and non-perturbative effects would otherwise have to do, so the model supplies a practical baseline for heavy-flavor pair production in proton-proton collisions. A second claim is that the parent $b\\bar b$ observables are more sensitive to this broadening than the $J/\\psi$ pair observables, because the $B \\to J/\\psi$ decay randomizes the $J/\\psi$ direction relative to its parent hadron. On this basis the paper identifies which correlated observables best separate cold nuclear matter effects: pair rapidity responds mainly to changes in fragmentation, while azimuthal separation responds mainly to additional broadening. A trustworthy baseline lets a single set of measurements in proton-nucleus and nucleus-nucleus collisions disentangle competing nuclear effects.","feed_headline":"Single Gaussian kick reproduces LHCb bottom-pair correlations","feed_subtitle":"NLO calculation with one transverse-momentum smearing matches every J/psi pair observable LHCb measured.","key_machinery":"The load-bearing object is the exclusive HVQMNR Monte Carlo for heavy-quark pair production at next-to-leading order, combined with a Gaussian transverse-momentum smearing, $g_p(k_T) = \\frac{1}{\\pi\\langle k_T^2\\rangle} e^{-k_T^2/\\langle k_T^2\\rangle}$, applied to the pair in the final state, and Peterson fragmentation with $\\epsilon_P = 0.0004$. The average kick grows logarithmically with energy, $\\langle k_T^2\\rangle = 1 + \\frac{\\Delta}{3}\\ln(\\sqrt{s}/20~\\mathrm{GeV})~\\mathrm{GeV}^2$, with $\\Delta=1$ as the default for proton-proton collisions, giving about $3\\mathrm{GeV}^2$ at 7 TeV. This single smearing stands in for the resummation that the fixed-order code does not include, softening the back-to-back peak in azimuth and filling the low pair-$p_T$ region; the paper notes that the initial-state and final-state implementations of the kick are exactly equivalent only at leading order.","core_discovery":"The central claim is that an exclusive next-to-leading-order treatment of $b\\bar b$ production with a single Gaussian intrinsic $k_T$ broadening and Peterson fragmentation describes every LHCb measurement of $b\\bar b \\to J/\\psi J/\\psi$ correlations at 7 and 8 TeV, for minimum $J/\\psi$ transverse momenta of 2, 3, 5 and 7 GeV in the forward rapidity range $2<y<4.5$. The same calculation reproduces the shapes of $|\\Delta\\phi|$, $|\\Delta y|$, the pair rapidity, the transverse momentum asymmetry $A_T$, the pair transverse momentum, and the pair mass. The paper's key interpretive result is that the parent $b\\bar b$ distributions retain a clear imprint of the $k_T$ broadening, while the $J/\\psi$ pair distributions are almost unaffected, because the decay randomizes the $J/\\psi$ momentum direction relative to the parent $B$ meson. Therefore measurements of $J/\\psi$ decay products alone are not a good probe of initial-state broadening; direct $B$-meson pair observables would be. The paper also shows that mass and scale variations change normalization more than shape, and that its nuclear-modification scenarios separate effects: pair rapidity tracks fragmentation, while azimuthal separation tracks broadening.","pith_inferences":["If the decay-decorrelation argument is correct, the same Gaussian-smearing model applied to charm pairs should produce much stronger azimuthal modifications at accessible LHCb transverse momenta, since for charm the kick width is comparable to the quark mass; this gives a sharper test of the mechanism than bottom pairs do.","The illustrative nuclear scenarios imply a quantitative strategy for heavy-ion data: a simultaneous fit of nuclear modification factors as functions of pair rapidity and azimuthal separation could separate energy loss from transverse-momentum broadening without relying on a full energy-loss calculation.","A dedicated resummed or parton-shower-matched calculation with explicit all-order treatment of low pair transverse momentum would show whether the single Gaussian kick is absorbing a real physical effect or merely providing a tuned stand-in."],"forward_implications":["The $p+p$ shapes of all six LHCb pair observables are stable between 7 and 8 TeV, so calculations at a single energy can be compared with the combined data.","Because the $J/\\psi$ decays decorrelate the final state from the parent pair, future measurements of reconstructed $B$-meson pairs would expose more of the intrinsic broadening signal that the decay channel hides.","Measured pair-rapidity and azimuthal-separation ratios in proton-nucleus and nucleus-nucleus collisions can be used together: the first is sensitive mainly to fragmentation and energy loss, the second mainly to $k_T$ broadening.","Mass and scale variations chiefly change the normalization of the pair distributions rather than their shapes, so comparisons of distribution shapes to data are robust against these theoretical uncertainties."],"supporting_citations":[{"why":"Supplies the exclusive next-to-leading-order heavy-quark pair Monte Carlo code (HVQMNR) in which all pair observables are computed.","marker":"[13]"},{"why":"Defines the model parameters adopted here: bottom mass, scale choices, the Peterson parameter, and the energy-dependent Gaussian broadening.","marker":"[14]"},{"why":"Provides the LHCb measurements of $b\\bar b \\to J/\\psi J/\\psi$ correlations at 7 and 8 TeV that all calculations are compared against.","marker":"[15]"},{"why":"A positive-weight next-to-leading-order Monte Carlo whose agreement with data LHCb interpreted as small NLO effects; this paper argues an NLO calculation plus broadening is the right test.","marker":"[21]"},{"why":"Supplies the Peterson fragmentation function used to hadronize bottom quarks into $B$ mesons, with the parameter later modified to model energy loss.","marker":"[26]"},{"why":"Supplies the Gaussian $k_T$ smearing function and the argument that adding the kick in the initial or final state is equivalent when the kick is not too large.","marker":"[31]"},{"why":"Supplies the $B \\to J/\\psi X$ branching ratio used to normalize the $J/\\psi$ pair yields.","marker":"[34]"},{"why":"Supplies the EPS09 nuclear parton distribution functions used to model shadowing in the proton-nucleus and nucleus-nucleus cold-matter scenarios.","marker":"[42]"}],"fun_headline_variants":["One Gaussian smearing matches all LHCb b-bbar pair data","Parent b-bbar keeps broadening signal; J/psi pairs wash it out","Single-kT model fits every LHCb bottom-pair correlation","NLO with one kT kick nails LHCb bottom-pair observables"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim rests on the assumption that one Gaussian transverse-momentum kick, applied with a single energy-dependent width to the produced pair, is an adequate stand-in for all non-perturbative and resummation effects on the pair's transverse momentum; if that equivalence fails at next-to-leading order, the calculated shapes can change.","fun_headline_variants_meta":{"raw":{"variants":["One Gaussian smearing matches all LHCb b-bbar pair data","Parent b-bbar keeps broadening signal; J/psi pairs wash it out","Single-kT model fits every LHCb bottom-pair correlation","NLO with one kT kick nails LHCb bottom-pair observables"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000732,"raw_usage":{"total_tokens":3409,"prompt_tokens":1211,"completion_tokens":2198,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":827,"completion_tokens_details":{"reasoning_tokens":2117}},"tokens_in":827,"tokens_out":2198,"duration_ms":16108,"temperature":1.0,"reasoning_tokens":2117,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:17:19.420091+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of $B$-hadron pair correlations with the same transverse-momentum cuts, without the $J/\\psi$ decay in between, would test the decay-decorrelation claim: if the parent $B$ pairs show no stronger dependence on broadening than the $J/\\psi$ pairs do, the central mechanism is wrong. Alternatively, a resummed or parton-shower-matched next-to-leading-order calculation with explicit all-order treatment of low pair transverse momentum should reproduce the same LHCb distributions; if it cannot match the data with the same single-kick prescription, the broadening model is inadequate.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exclusive next-to-leading-order heavy-quark pair Monte Carlo code (HVQMNR) in which all pair observables are computed."},{"cited_title":"Vogt, Heavy Flavor Azimuthal Correlations in Cold Nuclear Matter, Phys","cited_arxiv_id":null,"evidence_quote":"Defines the model parameters adopted here: bottom mass, scale choices, the Peterson parameter, and the energy-dependent Gaussian broadening."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Provides the LHCb measurements of $b\\bar b \\to J/\\psi J/\\psi$ correlations at 7 and 8 TeV that all calculations are compared against."},{"cited_title":"Norrbin and T","cited_arxiv_id":null,"evidence_quote":"Supplies the Peterson fragmentation function used to hadronize bottom quarks into $B$ mesons, with the parameter later modified to model energy loss."},{"cited_title":"Frixione, M","cited_arxiv_id":null,"evidence_quote":"Supplies the Gaussian $k_T$ smearing function and the argument that adding the kick in the initial or final state is equivalent when the kick is not too large."},{"cited_title":"Altarelli, G","cited_arxiv_id":null,"evidence_quote":"Supplies the $B \\to J/\\psi X$ branching ratio used to normalize the $J/\\psi$ pair yields."},{"cited_title":"As noted in Ref","cited_arxiv_id":null,"evidence_quote":"Supplies the EPS09 nuclear parton distribution functions used to model shadowing in the proton-nucleus and nucleus-nucleus cold-matter scenarios."}],"review_version":1}