{"id":"6575caba-f79e-4d99-a318-5465bdd45511","arxiv_id":"2412.12847","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A quantum chemistry study finds 44 alpha and 59 beta conformers of D-xylopyranose and proposes 10 per anomer covering >80% of the population between 298 and 1068 K.","lead":"This paper maps the shapes (conformers) of the sugar xylopyranose at temperatures used in biomass pyrolysis using high-level quantum chemistry. It identifies the most common conformers and their populations, giving kinetic modelers a shortlist of structures to use in pyrolysis simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 10-conformer, >80% population claim appears inconsistent with Table 6: at 1068 K, 22 α and 26 β conformers are needed to reach only ~90% of the full population, so the >80% figure is likely computed after renormalizing the 1%-cutoff-reduced space.","rationale":"I read the paper's central claim as the identification of 10 conformers per anomer that remain above 80% of the Boltzmann population over the pyrolysis temperature window, because that is the result the abstract, Section 3.3, and the conclusions all emphasize as the actionable outcome for future kinetic studies. The reader's weakest-assumption analysis focused on CREST/GFN2-xTB sampling completeness. That is a legitimate concern, but the more directly load-bearing problem is internal: the paper's own Table 6 says that many more than 10 conformers are needed to cover 85-92% of the full population at 678 and 1068 K. Section 3.2 explicitly describes renormalization after applying the 1% population cutoff, and Figure 9 appears to plot abundances within that reduced space. The 80% figure is therefore likely relative to the truncated set, not to the full conformational ensemble, which changes the meaning of the claim and the recommended conformer list. This is checkable from the SI without any new sampling or electronic-structure calculations, which is why it is the strongest single test. I do not recommend changing the reader's conditional verdict: the paper is plausible and valuable, but its headline deliverable should be verified and, if necessary, restated with the normalization made explicit once the SI is available.","tokens_in":18044,"tokens_out":3353,"duration_ms":31733,"concrete_test":"Using the SI population tables or by recomputing Boltzmann weights from the reported ΔG values, sum the full-space weights of the ten conformers named in Figure 8 for each anomer at 678 K and 1068 K without applying the 1% cutoff or any renormalization. If the summed share is below 80% for either anomer at either temperature, the central claim is an artifact of the reduced-space renormalization; if it is at or above 80%, the claim survives but must be restated with the normalization explicitly defined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central deliverable is the recommendation that 10 conformers per anomer capture more than 80% of the Boltzmann population across 298-1068 K, making them the kinetically relevant space for pyrolysis kinetics (Section 3.3, Figure 9, Conclusions). This claim is not supported by the paper's own full-population data. Section 3.2 and Table 6 report that at 1068 K, 92% of the α population requires 22 conformers (α-C1 to α-C22) and 85% of the β population requires 26 conformers (β-C1 to β-C26); at 678 K the corresponding numbers are 16 and 21 conformers. If those cumulative lists are correct, ten conformers cannot maintain >80% of the full Boltzmann population at the upper temperatures. The only way the >80% figure can hold is if the populations plotted in Figure 9 are renormalized within the reduced space obtained after discarding conformers with relative population below 1%; Section 3.2 states that this cutoff reduces the space to 23 α and 31 β conformers and that percentages are 'renormalized with respect to the reduced conformational space.' Eighty percent of that truncated, renormalized space is not equivalent to 80% of the full conformational population. The abstract and conclusions state the 10-conformer claim without this normalization caveat, so the paper's key practical recommendation is currently an overstatement by its own numerical tables. This concern is independent of whether the CREST search was complete: even if the 44/59 conformer sets are exhaustive, the size of the proposed reduced set is not justified by the reported full-space populations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a computational conformational analysis of alpha- and beta-D-xylopyranose in the gas phase over the temperature range 298–1068 K. Conformers are generated with the CREST/GFN2-xTB metadynamics workflow, reoptimized and characterized at the revDSD-PBEP86-D3(BJ)/jun-cc-pVTZ level, with single-point energies refined at the DLPNO-CCSD(T)-F12/cc-pVTZ-F12 level. Thermochemical properties are computed with the mRRHO approximation, and Boltzmann populations are derived from relative Gibbs free energies. The authors report 44 alpha and 59 beta conformers, identify a 1% population cutoff that reduces the space to 23 and 31 conformers, and recommend a subset of 10 conformers per anomer as \"kinetically relevant\" for future multiconformational kinetic studies. They also report an anomeric ratio of 61:39 favoring alpha at 298 K, consistent with some prior gas-phase studies.","tokens_in":18400,"tokens_out":3740,"duration_ms":31489,"significance":"If the results are validated, the paper offers a useful conformer library and thermochemical dataset for xylopyranose, a key hemicellulose model compound. The electronic-structure protocol is high in quality for this class of molecules: double-hybrid DFT geometries and frequencies with DLPNO-CCSD(T)-F12 single-point energies represent a defensible level of theory, and the Boltzmann analysis is standard. The comparison of the anomeric ratio with literature values is a welcome check. However, the central practical deliverable—the claim that 10 conformers per anomer capture more than 80% of the conformational population—is, as stated in the abstract and conclusions, not supported by the paper's own full-population data. This overstatement, together with the absence of the supplementary tables and coordinates that the text repeatedly references, makes the manuscript unsuitable for acceptance in its current form.","major_comments":[{"comment":"The claim that 10 conformers per anomer maintain more than 80% of the Boltzmann population across 298–1068 K is inconsistent with the paper's own Table 6. At 1068 K, the table states that 22 alpha conformers (alpha-C1 to alpha-C22) are needed to describe 92% of the alpha population, and 26 beta conformers (beta-C1 to beta-C26) are needed to describe 85% of the beta population. Because these are cumulative subsets built from the most populated conformers, ten conformers cannot contain more than 80% of the full Boltzmann population at that temperature. The only interpretation consistent with Figure 9 is that the percentages are renormalized within the reduced conformational space obtained after discarding conformers with relative population below 1%, as stated in Section 3.2. The abstract and Conclusions present the 10-conformer, >80% statement without this caveat, overstating the practical recommendation. Please either revise the claim to specify the reduced-space normalization explicitly or recompute the coverage with respect to the full 44/59 conformer sets and report the actual number of conformers needed to reach a given coverage.","section":"Abstract, Section 3.2, Section 3.3, Table 6, Figure 9"},{"comment":"The completeness of the GFN2-xTB/CREST conformational search is not independently verified. The paper notes that a comprehensive search is computationally prohibitive and relies on the iMTD-GC algorithm, but no convergence check is reported (for example, repetition with different random seeds, comparison with an independent sampling method, or a plot of the number of new conformers versus simulation time). Since every downstream quantity—the energy ordering, the Boltzmann populations, and the proposed 10-conformer subsets—depends on the identified conformer set being effectively complete within the energy window of interest, the authors should provide some evidence of convergence or explicitly characterize the sensitivity of the reported populations to possible missing conformers.","section":"Section 2.1 and Section 3.1"},{"comment":"The manuscript contains multiple unresolved cross-references to SI tables: \"Table ??\" in Section 3.1.1, \"Table ??, reported in SI\" in Section 3.1.2, \"The SI tables ??, ?? report the numerical values\" in Section 3.1.3, and \"the Cartesian coordinates of the optimized geometries are reported in the SI\" in Section 3.3. The SI is not included in the arXiv submission, so the numerical thermochemical tables and Cartesian coordinates on which the population analysis is based cannot be inspected. Without these data, the conformer assignments and the key numerical results (including the 10-conformer coverage) are not independently verifiable. Please provide the SI or include the essential numerical tables in the main text for review.","section":"Throughout (Sections 3.1.1, 3.1.2, 3.1.3, 2.2, Figure 8 caption)"}],"minor_comments":[{"comment":"The word \"neglectible\" should be \"negligible\" in the sentence describing the population deviation after applying the cutoff.","section":"Section 3.2"},{"comment":"The temperature is given as \"673K\" in the Conclusions paragraph (\"at 298 K, 673K and 1068 K\"), whereas 678 K is used throughout the rest of the manuscript; please unify.","section":"Section 4"},{"comment":"The sentence introducing Table 5 refers to \"the conformer structures most affected by entropic effects,\" but Table 5 is a list of conformer labels rather than structures; please rephrase for clarity.","section":"Section 3.1.3"},{"comment":"The vibrational frequency scale factor of 0.982 is said to be \"computed using the FREQ program.\" Please specify the underlying reference data and the exact functional/basis set for which this scale factor was derived, and cite the corresponding entry in the FREQ database.","section":"Section 2.2"},{"comment":"In the reproduced version of Figure 7, the horizontal lines indicating the population subsets are difficult to read; please ensure sufficient resolution and contrast in the final version.","section":"Figure 7"},{"comment":"The row for Schmidt et al. in the gas phase reports Delta-G = 0.05 kcal/mol favoring alpha (52/48), while the text in Section 1 states that \"experimental values slightly favored the beta form\"; this discrepancy between the table and the narrative should be clarified.","section":"Table 4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has real strengths—high-level electronic structure, a standard and transparent Boltzmann analysis, and a direct comparison of the anomeric ratio with earlier work. The main issue is the overstatement of the 10-conformer coverage claim in the abstract and conclusions, which is internally inconsistent with Table 6. This is fixable by rewording or reanalysis. The persistent 'Table ??' cross-references and the absence of the SI are also blocking for a proper evaluation of reproducibility. I would be willing to review a revised version that addresses these points. I do not see evidence of deliberate misrepresentation; rather, the paper appears to conflate the reduced-space renormalization with full-space population coverage."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, useful computational study, and the main deliverable—a DLPNO-quality conformer library with temperature-dependent populations for both xylopyranose anomers—is worth having. But the paper's headline convenience claim is currently overstated. The 10-conformer subsets that 'remain above 80%' are 80% of the truncated space left after removing conformers below 1% population, not 80% of the full Boltzmann population. By the paper's own Table 6, at 1068 K you need 22 α-conformers to reach 92% and 26 β-conformers to reach 85%. The text and conclusions don't consistently carry the 'renormalized' caveat, so a reader will walk away thinking ten conformers capture four-fifths of the real population at pyrolysis temperatures. They don't. That's a fixable framing problem, but it's the central practical claim, so it needs to be fixed before anyone builds a kinetic model on it.\n\nWhat is genuinely new: no earlier study gives both anomers at this level (rDSD geometries/frequencies, DLPNO-CCSD(T)/F12 single points, mRRHO thermochemistry) over 298–1068 K. The 44/59 conformer sets, the puckering labels, and the explicit population tables are a useful reference. The anomeric ratio check against experiment (61/39 in favor of α) is a reasonable validation, and the population calculation is definitional, not a fit, so there is no circularity problem.\n\nSoft spots beyond the normalization issue: the SI is missing from the version I saw, including Cartesian coordinates and the numerical thermochemical tables, so independent verification is impossible right now; there are broken cross-references and at least one temperature inconsistency (678 vs 673 K); and there is no uncertainty estimate on the populations, which matters because the mRRHO treatment and the 0.982 frequency scale factor carry real model dependence. The completeness of the CREST/GFN2-xTB search is the usual unknown; it would be good to see at least a statement about replicate sampling or a test of the iMTD-GC parameters, but I wouldn't call that a fatal flaw—it's inherent to this kind of conformational exploration.\n\nBottom line: this paper deserves referee time. The calculation is careful, the data are useful, and the main flaw is a reporting/normalization problem rather than a methodological one. For a reader building multipath pyrolysis kinetics, the work is valuable even after the caveat is made explicit—you just need to either soften the 10-conformer claim or present it honestly as '10 conformers capture >80% of the reduced space, and the reduced space itself retains roughly 94% of the full population at 1068 K.' I'd send it to review, ask for the SI, and ask for a rewrite of Sections 3.3 and the abstract to align with Table 6.","headline":"Useful DLPNO-quality conformer library for both xylopyranose anomers, but the 10-conformer '>80% population' claim is an artifact of renormalizing to the 1%-cutoff-reduced space and needs to be reframed.","tokens_in":18955,"tokens_out":3936,"would_cite":true,"duration_ms":33224,"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":"A subset of ten conformers per anomer retains more than 80% of the Boltzmann population of D-xylopyranose from 298 to 1068 K, defining the kinetically relevant conformational space for pyrolysis kinetics.","keywords":["xylopyranose","conformational analysis","hemicellulose pyrolysis","metadynamics","GFN2-xTB","Boltzmann population","mRRHO","anomeric ratio"],"falsifier":"Re-run the conformational search with different metadynamics seeds or a longer simulation, or scan the hydroxyl dihedral grid at the rDSD level, and check whether any new conformer appears within about 1 kcal mol−1 of α-C1 (or enters the top ten in free energy at 1068 K). Finding such a missed minimum would directly falsify the claim that the proposed ten conformers per anomer cover >80% of the equilibrium population across the pyrolysis temperature range.","tokens_in":17865,"feed_emoji":"🔥","tokens_out":9787,"duration_ms":80016,"temperature":0.7,"pith_summary":"This paper asks which three-dimensional shapes of D-xylopyranose, the sugar unit of hemicellulose, actually exist in the gas phase at temperatures where biomass undergoes fast pyrolysis (298-1068 K). Using a coarse semi-empirical search followed by high-level refinement, it finds 44 distinct conformers of the α anomer and 59 of the β anomer within about 10 kcal/mol of each anomer's lowest energy structure. The authors then use Boltzmann statistics on computed free energies to show that a much smaller set, ten conformers per anomer, always contains more than 80% of the equilibrium population across the whole pyrolysis temperature range. The practical point is that future ab initio kinetic studies of xylopyranose decomposition can restrict themselves to these twenty structures without losing the chemically relevant reactive states. The result also gives a gas-phase α/β anomeric ratio of 61:39 at 298 K, consistent with earlier reports of a slight preference for the α form in the absence of solvent.","feed_headline":"Ten conformers per anomer capture 80% of the reactive population","feed_subtitle":"Metadynamics plus high-level thermochemistry trims 44 α and 59 β conformers to the 20 that matter at 298–1068 K.","key_machinery":"The carrying mechanism is the Boltzmann population distribution built on Gibbs free energies rather than bare electronic energies. Populations are computed as $N_i/N = g_i e^{-\\Delta G_i/k_B T} / \\sum_j g_j e^{-\\Delta G_j/k_B T}$ using mRRHO free energies, where the modified rigid-rotor harmonic-oscillator treatment of Grimme supplies anharmonic corrections for hindered hydroxyl rotations and low-frequency ring modes. A population cut-off (discarding conformers below 1% relative abundance) reduces the 44/59 conformer lists to 23/31, and a further selection isolates the ten per anomer that maintain >80% of the population across all three reference temperatures; these selections are then proposed as the input manifold for ab initio kinetic calculations.","core_discovery":"The central discovery is a rigorously reduced conformational space for both anomers of D-xylopyranose under pyrolysis conditions. The metadynamics search with the GFN2-xTB Hamiltonian generated 57 α and 73 β candidate minima; after reoptimization and vibrational analysis at the revDSD-PBEP86-D3(BJ) level, and single-point energies at DLPNO-CCSD(T)/F12, 44 α and 59 β distinct conformers remained, spanning roughly 10 kcal mol−1 in zero-point-corrected electronic energy. Using modified rigid-rotor harmonic-oscillator free energies, the authors computed Boltzmann populations at 298, 678, and 1068 K and found that the five lowest-energy conformers hold 95% of the α population at 298 K (97% for β), with the population spreading to higher-energy conformers as temperature rises. From this analysis they define a selection of ten conformers per anomer whose combined relative abundance stays above 80% over the entire 298-1068 K range, and they argue that this twenty-conformer set is the kinetically relevant conformational space for multiconformational ab initio studies of pyrolysis reaction kinetics.","pith_inferences":["The same reduced-set workflow could be applied to xylofuranose or to glucopyranose to test whether ten conformers per monomer is a general rule for pyranose sugars at high temperature, or an accident of xylose's symmetry.","The population cut-off is thermodynamic, not kinetic; whether the ten-conformer set also carries the reaction flux depends on transition-state energies, so a natural next test is to compute key decomposition rate constants with the full and reduced sets and compare the resulting product yields.","Because the paper notes solvent reverses the α/β preference, adding implicit solvation or microsolvated water could shift the ten-conformer selection; a useful extension would map how the subset changes with water activity relevant to condensed-phase pyrolysis."],"forward_implications":["Kinetic models of hemicellulose pyrolysis can start with these twenty conformers instead of a single lowest-energy structure, so computed barriers and rate constants will reflect the actual reactive population.","At 1068 K roughly 35% of the population leaves the low-energy manifold; any single-conformer or few-conformer treatment would miss that fraction, so multiconformational rate calculations are necessary at pyrolysis temperatures.","The gas-phase α/β ratio of 61:39 at 298 K provides a clean benchmark for force-field and implicit-solvent models, which must reproduce this intrinsic preference before adding solvation.","The temperature dependence of relative free energies identifies which conformers change rank with temperature, guiding where explicit anharmonicity or hindered-rotor treatments matter most."],"supporting_citations":[{"why":"Provides the CREST program and iMTD-GC algorithm that performs the RMSD-metadynamics search generating the initial conformer pools.","marker":"30"},{"why":"Defines the GFN2-xTB tight-binding Hamiltonian used as the fast potential in the conformational search.","marker":"31"},{"why":"Describes the RMSD-based metadynamics biasing approach that drives exploration of the potential energy surface.","marker":"32"},{"why":"Supplies the revDSD-PBEP86 double-hybrid functional used for geometry reoptimization and harmonic vibrational analysis.","marker":"39"},{"why":"Provides the DLPNO-CCSD(T) method used to refine electronic energies of the DFT-optimized conformers.","marker":"44"},{"why":"Introduces the modified rigid-rotor harmonic-oscillator (mRRHO) approximation used to compute anharmonic free energies and Boltzmann populations.","marker":"56"},{"why":"Earlier MD/FEP study of the α/β anomeric equilibrium in gas and solution, providing the comparison value for the computed anomeric ratio.","marker":"4"},{"why":"Microwave-spectroscopy study of α/β-D-xylopyranose conformers whose hydrogen-bonding analysis is used to validate the stability ordering.","marker":"10"}],"fun_headline_variants":["20 conformers capture 80% of xylopyranose's reactive states","Metadynamics plus CCSD(T) pin down 20 key xylopyranose conformers","Xylopyranose pyrolysis: just 20 conformers dominate 80% of population","Reducing xylopyranose's conformational space to 20 for kinetics","Ten α and ten β conformers cover most pyrolysis-relevant states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole reduction is only as complete as the initial semi-empirical search: if CREST/GFN2-xTB metadynamics missed any low-energy minimum, then the energy ranking, the Boltzmann weights, and the final ten-conformer subsets could all be missing a kinetically important structure.","fun_headline_variants_meta":{"raw":{"variants":["20 conformers capture 80% of xylopyranose's reactive states","Metadynamics plus CCSD(T) pin down 20 key xylopyranose conformers","Xylopyranose pyrolysis: just 20 conformers dominate 80% of population","Reducing xylopyranose's conformational space to 20 for kinetics","Ten α and ten β conformers cover most pyrolysis-relevant states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1650,"prompt_tokens":1059,"completion_tokens":591,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":486}},"tokens_in":675,"tokens_out":591,"duration_ms":5638,"temperature":1.0,"reasoning_tokens":486,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:39:37.958619+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the conformational search with different metadynamics seeds or a longer simulation, or scan the hydroxyl dihedral grid at the rDSD level, and check whether any new conformer appears within about 1 kcal mol−1 of α-C1 (or enters the top ten in free energy at 1068 K). Finding such a missed minimum would directly falsify the claim that the proposed ten conformers per anomer cover >80% of the equilibrium population across the pyrolysis temperature range.","supporting_citations":[],"review_version":1}