{"id":"2cbb8803-6833-4018-abd8-3cae87e2b343","arxiv_id":"2504.16338","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A survey plus curated database of electron mobility and diffusion data for liquid argon and xenon, arguing that liquid-specific coherent scattering and screened potentials are required to match TPC measurements.","lead":"This review compiles experimental and theoretical results on how electrons drift and diffuse through liquid argon and xenon, and debuts an open database of past measurements for detector simulation. It argues that gas-phase scattering data alone fail to describe liquid behavior, and that structure-aware ab initio models are needed.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim of quantitative LAr transport agreement without empirical tuning is under-constrained: the characteristic-energy check rests on one dataset with unmeasured temperature, and no sensitivity analysis is given for key model choices such as the matching radius rm.","rationale":"The paper's central claim (Section 3.2.3) is that a parameter-free ab initio framework reproduces LAr transport coefficients, supported almost entirely by Figure 6, a reproduction from the authors' earlier papers. The two most load-bearing vulnerabilities are (i) the absence of any sensitivity analysis for the non-trivial model inputs, notably the matching radius rm, for which different physically motivated definitions exist in the literature (Section 3.2.2), and (ii) the fact that the characteristic-energy validation rests on a single liquid-argon measurement whose temperature is explicitly stated to be unmeasured. These are not internal inconsistencies, but they prevent the 'no empirical tuning' assertion from being verified to the required standard: a small change in rm or a plausible temperature shift could move the calculated curves outside the experimental scatter. The reader's weakest-assumption analysis focused on inelastic scaling; I agree that the inelastic treatment is also uncertain, but for the specific central claim about transport coefficients at low field, the elastic validation gap is more directly load-bearing. The suggested concrete test (varying rm within the range of published definitions) would directly probe whether the agreement is robust. Accordingly, the correct verdict remains CONDITIONAL: the review is valuable and scientifically sound, but the central quantitative claim needs supporting sensitivity and uncertainty analysis before it can be taken as established.","tokens_in":35720,"tokens_out":6785,"duration_ms":68221,"concrete_test":"Recompute the Liq+Coh drift velocity and characteristic energy in Figure 6 with the matching radius rm set to the Wigner-Seitz radius rm = (4πN/3)^{-1/3} ≈ 4.2 a0 instead of the first turning point of Ueff (≈ 4.3 a0), keeping all other inputs fixed. If either transport coefficient changes by more than the experimental uncertainty (or the model-experiment discrepancy) in Figure 6, the result is sensitive to a physically reasonable modeling choice and the 'no empirical tuning' claim is not robust. A secondary check is to repeat the characteristic-energy comparison using a density consistent with the plausible temperature range of the Shibamura dataset (e.g., 84–90 K) and verify whether the single liquid data point remains within the model band.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Section 3.2.3) is that the structure-modified, screened-potential framework reproduces LAr transport coefficients 'without the need for empirical tuning.' Two gaps prevent this claim from being established on the evidence presented. First, Figure 6 is reproduced from Boyle et al. (2015) with no in-paper verification and no sensitivity analysis for the model's key structural choices. In particular, the matching radius rm enters the scattering calculation (Section 3.2.2) and is defined by a specific criterion (first turning point of Ueff), but alternative criteria exist in the literature (e.g., the Wigner-Seitz radius used by Atrazhev and collaborators), yielding nearby but different values; the dependence of the converged transport coefficients on rm is not shown. If a ±5% change in rm shifts the computed drift velocity or characteristic energy by more than the experimental scatter, the 'no-empirical-tuning' claim would fail even if the fixed-rm result agrees. Second, the characteristic-energy panel of Figure 6 compares the liquid model (85 K) to only one liquid experimental dataset, Shibamura et al. (1979), whose caption notes an 'unmeasured liquid temperature'; the other characteristic-energy data shown are gas-phase at 77 K and 288 K. An unknown temperature implies unknown density and structure factor, so the agreement of that single point with the 85 K calculation is not a controlled validation. Together, the quantitative agreement claimed is under-constrained: one key observable is validated against an uncharacterized datum, and the elastic model's parameter sensitivity is not examined. Hence the central claim, as stated, is not yet demonstrated to the standard implied.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reviews the experimental and theoretical understanding of electron transport in liquid argon and xenon, with emphasis on time projection chamber applications. It introduces an open-access database of mobility and diffusion measurements, surveys historical and modern swarm experiments, and contrasts empirical models (NEST, BNL) with the authors' multi-term Boltzmann framework that incorporates structure-modified elastic cross-sections, screened polarization potentials, and bulk potential modifications. The central theoretical claim is that this framework reproduces measured transport coefficients in liquid argon quantitatively without empirical tuning (Section 3.2.3). The review also discusses inelastic scattering approximations, electron self-trapping, scintillation, doped mixtures, and gas-liquid interfaces.","tokens_in":36097,"tokens_out":4441,"duration_ms":40911,"significance":"If the central claim holds, the paper offers a parameter-light ab initio alternative to empirical transport models for noble-liquid detectors, and the open-access database is a useful community resource for benchmarking swarm data. The review is careful in flagging inconsistent legacy measurements and in distinguishing flux from bulk transport coefficients, and it explicitly acknowledges the limitations of current inelastic models, for which no validated liquid cross-sections exist. These strengths make the review informative. However, the quantitative validation of the central claim is imported from self-cited prior papers and is not accompanied by in-paper error or sensitivity analysis, which limits the weight that can be placed on the 'no empirical tuning' assertion.","major_comments":[{"comment":"The central claim that the structure-modified, screened-potential framework reproduces LAr transport 'without the need for empirical tuning' is under-constrained by the evidence shown. Figure 6 is reproduced from Boyle et al. (2015) without in-paper verification, and no sensitivity analysis is provided for the matching radius r_m defined in Section 3.2.2, which is set by the first turning point of U_eff. Because alternative criteria, such as the Wigner-Seitz radius used by Atrazhev and collaborators, give nearby but different values, the dependence of W and (3/2)D_T/µ on r_m should be shown; if a small change in r_m shifts the computed coefficients by more than the experimental scatter, the quantitative agreement is not robust. In addition, the characteristic-energy panel of Figure 6 compares the 85 K liquid calculation to only one liquid dataset, Shibamura et al. (1979), whose caption states the liquid temperature was unmeasured; an unknown temperature implies unknown density and structure factor, so this single comparison is not a controlled validation. I recommend either adding a sensitivity and uncertainty analysis or tempering the 'no empirical tuning' claim to what the reproduced figure actually supports.","section":"Section 3.2.3, Figure 6"},{"comment":"The inelastic-scattering treatment assumes that liquid-phase cross-sections can be approximated by gas-phase cross-sections shifted to the measured liquid band gap, with atomic excitation cross-sections used for excitons and n>1 Wannier excitons omitted. The paper itself concedes that 'the fraction of clusters that support exciton formation versus those contributing only perturbed atomic transitions remains unknown' and that the ionization data used for validation in Figure 7 are limited to low reduced fields. Since the rate coefficients shown in Figure 8 and the scintillation discussion in Section 4.1 depend on these inelastic cross-sections, the quantitative part of the theoretical narrative for inelastic processes and high-field behavior is not established. The limitation is acknowledged, but it should be presented as an open problem rather than as part of the validated framework described in the abstract.","section":"Section 3.2.4"},{"comment":"The new database is a central contribution, but many entries have missing density and pressure, and the text states that densities are estimated along the saturated liquid line from NIST when not reported. Since the reduced field E/N is obtained using these densities, an unquantified density estimate propagates directly into the quantities plotted in Figures 1-2. The repository should document, for each entry, whether density was reported or estimated, the estimation method, and an uncertainty estimate; without this, the transparent benchmarking goal is only partially met.","section":"Section 2.1 and Tables 1-2"}],"minor_comments":[{"comment":"The Introduction contains the duplicated phrase 'for for'; please correct it.","section":"Section 1"},{"comment":"The notation in Equation (18) should be checked: the screening function f(t) and the polarizability alpha_d(t) are used in a product inside the integral, but the distinction between the screening function and the dipole polarizability is not consistently annotated in the surrounding text.","section":"Section 3.2.2"},{"comment":"The legend labels contain typos ('Hal ern', 'Towse dBaile') and the reference 'Townsend and A., B. V.' is malformed; these should be corrected.","section":"Figure 6 caption"},{"comment":"The Author Contributions list contains 'IS' twice, with one entry appearing to be a duplicate; please remove the redundant entry.","section":"Author Contributions"},{"comment":"The text uses 'Wigner-Seitz diameter' in one place and 'Wigner-Seitz cell radii' earlier in the paper; please standardize the terminology.","section":"Section 3.2.2"}],"recommendation":"major_revision","confidential_remarks":"The review leans heavily on the authors' own prior publications for the central validation (Boyle et al. 2015, 2016; White et al. 2018; Simonović et al. 2019; Boyle et al. 2024). For a review article this is acceptable, but the editor may want to ensure that the 'quantitative agreement' claim is not overstated relative to the reproduced figures. The manuscript is within scope for physics.app-ph and would benefit from a shortened theory derivation or an explicit statement that Equations (13)-(22) are standard from the cited literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid, useful review. The genuinely new artifact is the GitHub database of measured mobilities and diffusion coefficients in LAr and LXe, with reduced quantities and temperature grouping. That will help detector developers and modelers who currently have to hunt through scattered old papers. The survey of experiments is careful and critical, particularly about legacy data quality and the empirical nature of NEST, and the inelastic-scattering section honestly concedes how speculative the liquid-phase cross-section scaling is. The theoretical synthesis of the authors' own multi-term Boltzmann work is coherent and well placed in context. So the paper earns its place as a reference review. The soft spots are real but not fatal. The central claim that the structure-modified, screened-potential model reproduces LAr transport 'without empirical tuning' is under-constrained. Figure 6 is imported from Boyle et al. (2015) with no in-paper verification, no sensitivity analysis on the matching radius rm, and the characteristic-energy panel leans on a single liquid data point whose temperature, and therefore density and structure factor, is unmeasured. I checked the figure caption in the text: Shibamura et al. is indeed listed as 'unmeasured liquid temperature.' That does not void the drift-velocity agreement, which has several well-characterized liquid datasets, but it does mean the diffusion comparison is not a controlled test. On the 'no empirical tuning' phrasing: the model does involve reasonable physical choices (matching radius criterion, band-gap shifts, density estimates), so the claim should be softened to something like 'no cross-section fitting.' The database tables also omit uncertainty values and versioning, which limits benchmarking value. Minor: no sensitivity to rm is shown, and the reader's concern about excimers and n>1 Wannier excitons being omitted is exactly what the paper itself flags, so that is an acknowledged limitation rather than an oversight. On balance, the review is honest, citation is thorough, and the central transport framework is not invalidated by these gaps. The audience is detector physicists and swarm theorists; a serious referee should engage with it. My recommendation: accept after revision, with a request to either add a brief sensitivity study for rm or explicitly qualify the quantitative-agreement claim, and to add uncertainty and versioning to the database.","headline":"A genuinely useful review and curated database, but the paper's flagship ab initio validation claim is weaker than the prose suggests.","tokens_in":36621,"tokens_out":1410,"would_cite":true,"duration_ms":17435,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper shows that a multi-term Boltzmann solver with structure-modified elastic cross-sections and screened, bulk-modified potentials reproduces measured electron transport coefficients in liquid argon without empirical tuning.","keywords":["electron transport","noble liquids","liquid argon","liquid xenon","Boltzmann equation","coherent scattering","swarm experiments","time projection chambers"],"falsifier":"Measure the drift velocity and longitudinal and transverse diffusion of electrons in ultrapure liquid argon at 85–90 K over $10^{-4}$ to $10^{-1}$ Td with independent density and impurity characterization. If the inferred momentum-transfer cross-section differs from the Liq+Coh prediction by more than the combined uncertainties, the screening and coherence model fails. Separately, a high-field ionization-coefficient measurement in liquid xenon above about 10 Td would discriminate between the four inelastic scenarios the paper compares, since the existing data stop at low reduced fields.","tokens_in":35539,"feed_emoji":"⚛️","tokens_out":8548,"duration_ms":80622,"temperature":0.7,"pith_summary":"This review makes a specific claim: low-energy electron drift and diffusion in noble liquids can be computed from the Boltzmann equation without empirical adjustment, provided the scattering input is made liquid-specific. The demonstration case is liquid argon, where gas-phase cross-sections fail completely, coherent scattering alone brings the drift velocity to the right order but overestimates the characteristic energy, and only the combination of coherent scattering, screened polarization, and the bulk potential reproduces both measured drift velocity and characteristic energy. The paper also compiles decades of liquid argon and xenon swarm measurements into a new open database with standardized reduced quantities, and it identifies inelastic processes—exciton formation, perturbed atomic excitations, interband transitions—as the unresolved frontier for scintillation and high-field behavior. A sympathetic reader would care because modern dark-matter and neutrino time projection chambers drift electrons over metre scales, and predictive transport coefficients directly affect signal reconstruction and detector design.","feed_headline":"Electron drift in liquid argon predicted without empirical tuning","feed_subtitle":"Once coherent scattering and screened potentials are both included, the model matches measured drift and diffusion.","key_machinery":"The load-bearing object is the structure-modified differential cross-section $\\Sigma(v,\\theta) = \\sigma(v,\\theta)\\,S(2m_e v/\\hbar \\sin(\\theta/2))$, which multiplies the binary electron–atom cross-section by the liquid's static structure factor $S(q)$ to encode coherent scattering. Around it sits an effective potential $U_{\\mathrm{eff}}(r) = U_1(r) + U_2(r)$, where $U_1$ is the focus atom's static potential plus its polarization potential screened by surrounding atoms (through Lekner's screening function $f(r)$, with the Lorentz screening factor at large $r$), and $U_2$ is the averaged contribution of the bulk atoms. Scattering phase shifts are evaluated at a matching radius $r_m$, avoiding artificial potential shifts and setting the conduction-band energy scale. These pieces feed a multi-term spherical-harmonic expansion of the Boltzmann equation, and the work they do is to make the low-energy elastic collision operator genuinely liquid rather than gas-like; the paper shows that dropping either the coherence factor or the potential modification breaks agreement with measured drift velocity and characteristic energy.","core_discovery":"The paper's central claim is that a multi-term solution of the Boltzmann equation, using structure-modified elastic cross-sections and ab initio potentials screened and shifted by the liquid environment, reproduces measured electron transport coefficients in liquid argon without tuning. In the elastic regime below a few Townsends, applying gas-phase cross-sections scaled by density is insufficient: the Ramsauer minimum is suppressed in the liquid, the momentum-transfer cross-section becomes nearly energy-independent at low energy, and both coherent scattering (encoded through the static structure factor) and the potential modifications must be included simultaneously for quantitative agreement. The same methodology is reported to extend to liquid xenon, liquid krypton, and positron transport in liquid helium. For inelastic channels the paper is deliberately cautious: liquid-phase excitation and ionization cross-sections are approximated by gas-phase cross-sections shifted to the measured band gap, exciton-forming clusters are not separated from perturbed atomic transitions, and the remaining validation data are limited to low reduced fields.","pith_inferences":["If the band-gap-shifted gas-phase cross-sections for inelastic scattering are validated against high-field drift, diffusion, and ionization data from large TPCs, the same scaling could be used to predict scintillation yields without invoking fitted light-yield parameters; this is an extension the paper leaves for future work.","The open database's standardized $E/N$ and $N\\mu$ entries could be used as training data for inverse-swarm machine-learning extraction of liquid-phase cross-sections, in parallel with the gas-phase deep-learning methods the paper cites, giving a testable route to close the inelastic gap.","The non-linear mixture behaviour at high packing fraction suggests that small dopant concentrations might be tuned to engineer mobility or diffusion, but only if the partial structure factors are known; that tuning knob is an implicit consequence of the review's mixture formalism."],"forward_implications":["Gas-phase scaling is not a valid short-cut at liquid densities: any transport model that ignores coherence and potential screening will mispredict low-field drift velocity and characteristic energy in liquid argon.","The same ab initio elastic machinery can be carried over to liquid xenon, liquid krypton, and positron systems, giving detector simulations a parameter-free alternative to empirically fitted mobility and diffusion.","In doped and mixed liquids, transport coefficients can lie outside the range bounded by the pure species, so mixture modeling requires partial static structure factors and cannot be interpolated from pure-fluid properties.","Modeling the gas–liquid interface as a smooth density gradient with a spatially varying conduction-band energy $V_0$ yields an effective field that can either inhibit or assist electron extraction in dual-phase time projection chambers.","Accurate inelastic rate coefficients for interband transitions and excitations, which become dominant above roughly $20$ Td in liquid xenon, are required inputs for scintillation and discharge modeling; they are not yet known from first principles."],"supporting_citations":[{"why":"Provides the liquid-argon transport calculation and the Gas, Gas+Coh, and Liq+Coh cross-section sets whose comparison with experiment carries the central claim.","marker":"Boyle et al. (2015)"},{"why":"Supplies the theory of hot electrons in dense media and the coherent-scattering modification of the elastic collision operator.","marker":"Cohen and Lekner (1967)"},{"why":"Introduces the screening function f(r), the Lorentz screening factor, the matching radius, and the V0 shift used to define liquid scattering potentials.","marker":"Lekner (1967)"},{"why":"The benchmark gas-phase electron–argon momentum-transfer cross-section used to validate the underlying interaction potential in the dilute limit.","marker":"Buckman et al. (2000)"},{"why":"Supplies the 85 K liquid-argon drift velocity data that the Liq+Coh calculation must match.","marker":"Miller et al. (1968)"},{"why":"Supplies additional 85 K liquid-argon drift velocity data used in the transport comparison.","marker":"Halpern et al. (1967)"},{"why":"Provides the cut-off atomic potential and Wigner–Seitz construction for electron scattering in atomic liquids, an alternative treatment the review builds on.","marker":"Atrazhev and Timoshkin (1996)"},{"why":"The empirical simulation framework whose fitted transport curves are compared against the ab initio results and against the compiled experimental data.","marker":"Szydagis et al. (2025)"},{"why":"Introduces the approach of shifting gas-phase cross-sections to the liquid band gap for inelastic channels in liquid argon.","marker":"Garland et al. (2018b)"}],"fun_headline_variants":["No-tuning model reproduces electron drift in liquid argon","Liquid argon electron drift explained without adjustable parameters","Ab initio model predicts electron drift in liquid argon","Liquid argon drift predicted with no empirical knobs","Coherent scattering and screened potentials fix liquid argon drift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the liquid can be represented by gas-phase electron–atom potentials plus classical structure-factor and screening corrections, and that inelastic channels can be approximated by gas-phase cross-sections shifted to the measured liquid band gap; if those proxies do not capture the actual liquid environment, the quantitative agreement and its extension to high fields would not follow.","fun_headline_variants_meta":{"raw":{"variants":["No-tuning model reproduces electron drift in liquid argon","Liquid argon electron drift explained without adjustable parameters","Ab initio model predicts electron drift in liquid argon","Liquid argon drift predicted with no empirical knobs","Coherent scattering and screened potentials fix liquid argon drift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000651,"raw_usage":{"total_tokens":2980,"prompt_tokens":934,"completion_tokens":2046,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":550,"completion_tokens_details":{"reasoning_tokens":1971}},"tokens_in":550,"tokens_out":2046,"duration_ms":15585,"temperature":1.0,"reasoning_tokens":1971,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:05:00.576959+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the drift velocity and longitudinal and transverse diffusion of electrons in ultrapure liquid argon at 85–90 K over $10^{-4}$ to $10^{-1}$ Td with independent density and impurity characterization. If the inferred momentum-transfer cross-section differs from the Liq+Coh prediction by more than the combined uncertainties, the screening and coherence model fails. Separately, a high-field ionization-coefficient measurement in liquid xenon above about 10 Td would discriminate between the four inelastic scenarios the paper compares, since the existing data stop at low reduced fields.","supporting_citations":[{"cited_title":"J., McEachran, R","cited_arxiv_id":null,"evidence_quote":"Provides the liquid-argon transport calculation and the Gas, Gas+Coh, and Liq+Coh cross-section sets whose comparison with experiment carries the central claim."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the theory of hot electrons in dense media and the coherent-scattering modification of the elastic collision operator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the screening function f(r), the Lorentz screening factor, the matching radius, and the V0 shift used to define liquid scattering potentials."},{"cited_title":"S., Howe, S., and Spear, W","cited_arxiv_id":null,"evidence_quote":"Supplies the 85 K liquid-argon drift velocity data that the Liq+Coh calculation must match."}],"review_version":1}