Case Studies of Using the Partial-Structure R1 and the Single-Atom R1 to Assemble Small-Molecule Crystal Structures

Xiaodong Zhang

REVIEW 4 major objections 5 minor 1 references

Case Studies of Using the Partial-Structure R1 and the Single-Atom R1 to Assemble Small-Molecule Crystal Structures

T0 review · 4 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper shows that small-molecule crystal structures can be assembled faster and more accurately when the single-atom R1 search is guided by pre-known connectivity and fragments, with partial-structure R1 reserved for low-resolution…

desk verdict A clear, honest case study of guided sR1/pR1 assembly with a genuinely useful connectivity-guided completion step; the general strategy claims outrun the evidence, but the paper deserves peer review with conditions. read the letter →

arxiv 2412.15284 v1 pith:D55JG56K submitted 2024-12-18 cond-mat.mtrl-sci physics.comp-phphysics.data-an

Xiaodong Zhang This is my paper · ORCID

classification cond-mat.mtrl-sciphysics.comp-phphysics.data-an

keywords partialstructureR1single-atommolecularreplacementlowdataresolutionconnectivity-guidedsearchsmall-moleculecrystalfragmentassemblyheavy-atomsubstructure

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper identifies practical strategies for assembling small-molecule crystal structures using the R1 search family. For light-atom-only crystals, the recommended route is normal single-atom R1 to build a framework, then connectivity-guided single-atom R1 to finish the model; only when data resolution is low should partial-structure R1 place known fragments first. For heavy-atom-containing crystals, the paper recommends normal single-atom R1 for the heavy-atom substructure followed by connectivity-guided single-atom R1. On three test structures the completed models place essentially all atoms within 0.5 Å of the correct positions, and the heavy-atom example runs in 1310 seconds with 5 misplaced atoms versus 3840 seconds and 18 misplaced atoms for the unguided route. The practical upshot is that pre-knowledge turns a blind search into an orderly, planned assembly.

What carries the argument

The load-bearing object is the connectivity-guided single-atom R1 search. It works by taking a known atom A and a target atom Q expected to bond to A at distance r, then testing only grid points inside a spherical shell of radii r−0.5 Å to r+0.5 Å around A before running the usual R1 calculation; this converts a whole-cell search into a one-bond extension. The companion pR1 method searches for a known fragment by finding the deepest hole of a partial-structure R1 map in a 6-dimensional orientation-location space, split into separate 3-dimensional orientation and location searches to save time. The paper also feeds sharpened intensities (Fo2 multiplied by exp(2Bs2), with B from a Wilson plot) into both methods, which reduces the number of ghost atoms. These pieces turn pre-known connectivity, fragments, and bond lengths into constraints that make each search step small and targeted.

What would settle it

On a previously unsolved light-atom structure with two identical known fragments and 1.2 Å data, place the first fragment by the deepest pR1 hole, define the pR1 for the second fragment, and compare the deepest hole's predicted position with the true position from a later full refinement; if the true position is not the deepest hole, the central heuristic is falsified.

Watch

Extended reading notes

Core claim

The paper's central claim is that the two R1 search modes should be combined according to crystal type: start with the normal single-atom R1 (sR1) to establish a framework or heavy-atom substructure, then switch to the connectivity-guided sR1, in which the search for a new atom is restricted to a spherical shell around a known bonding partner at an approximate bond length. The partial-structure R1 (pR1), which orients and positions entire known fragments by locating the deepest hole of an R1 map in a six-dimensional orientation-location space, is reserved for low-resolution cases where sR1 alone fails. The paper shows on sample 1 that when the resolution is truncated to 1.2 Å the normal sR1 can no longer solve the structure, while pR1-placed benzene-star fragments let sR1 complete it; on sample 3, the recommended strategy finishes in 1310 seconds with only 5 misplaced atoms, against 3840 seconds and 18 misplaced atoms for the unguided normal sR1 run. These results are offered as a workflow, not as a claim that pR1 is generally preferable.

Load-bearing premise

The deepest hole in the pR1 map over orientation-location space marks the true position and orientation of a missing fragment, and the split into separate three-dimensional searches preserves that signal; if either part fails on structures beyond these three cases, the recommended strategies collapse.

Editorial extensions

If this is right

For light-atom-only structures at ordinary resolution, normal sR1 followed by connectivity-guided sR1 completes the model faster than normal sR1 alone for the finishing step (446 s vs 710 s on sample 1) and can target specific missing atoms.
At low resolution (sample 1 truncated to 1.2 Å), the normal sR1 alone can no longer solve the structure, but pR1 placement of known fragments makes completion possible.
For heavy-atom-containing structures, the normal-sR1-then-connectivity-guided-sR1 strategy cut total time from 3840 s to 1310 s and misplaced atoms from 18 to 5 on sample 3.
pR1 should be avoided whenever sR1 can do the job, because orienting and positioning fragments with pR1 took 780 s and 3700 s respectively for sample 1, making the full pR1-assisted route slower than unguided sR1 at full resolution.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

At full resolution, unguided sR1 already solved sample 1 faster than the pR1-assisted workflow, so the practical value of pre-knowledge may lie mainly in low-resolution rescue and in targeting stubborn atoms rather than in speed at high resolution.
The 0.5 Å model-comparison metric ignores atom types, so reported success can coexist with type misassignments (sample 3 required manual correction of four I/Mo swaps); an automated pipeline would need a type-assignment step.
If the deepest-hole heuristic transfers, the same two-step pR1-to-connectivity-guided-sR1 plan could be tested on electron-diffraction data, where low resolution is common and fragment placement may become the standard first move.
A direct stress test would vary the assumed bond length r and shell width Δr to measure how sensitive the connectivity-guided sR1 is to imperfect pre-knowledge.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Here's my take on Zhang's sR1/pR1 case studies. The genuinely new pieces are the connectivity-guided sR1 completion step, the Wilson-sharpening tweak that cuts ghost atoms, and an explicit decision rule for when to use pR1 versus sR1. As a demonstration, the paper works: all three structures assemble correctly, the guided runs are faster than unguided sR1 on sample 3 (1310 s vs 3840 s, 5 vs 18 misplaced atoms), and the author is honest that pR1 is slow and should be avoided when sR1 can do the job. The writing is clear, and the SI gives enough algorithmic detail to reproduce the method.

The soft spots are real but proportionate. The central strategy recommendations rest on only three structures, and the low-resolution branch rests on a single unreported truncation test: the paper says that at 1.2 Å sample 1's sR1 fails but pR1-assisted sR1 succeeds, with no R1 values, orientation rankings, or timings shown. That's a one-line anecdote carrying a general claim. The 'deepest hole' heuristic from the author's prior paper is also unproven, and the whole method is validated against known structures rather than blind tests. There's no comparison with SHELX or other standard small-molecule software, so the practical wins over existing tools are not established. The 0.5 Å match tolerance is coarse—it treats a model as correct if every atom is within half an angstrom, which is a weak standard for a structure solution.

That said, these are not fatal flaws for a methods demonstration. The paper doesn't claim a new formalism; it claims a smarter way to use one. The author cites his own prior work appropriately, and the experimental setup is transparent about what was done and what was tweaked. The main problem is overgeneralization from a small sample, not sloppiness or circularity.

Who is this for? Crystallographers who work with difficult small-molecule structures where standard direct methods struggle. I'd send it to peer review, but conditionally: require release of the data and code, report the truncated-data pR1 results with proper R1-map rankings and failure statistics, and ideally add a couple more cases or a blind test. Without those, the 'necessary' and 'correct' strategy language should be softened to 'suggested'.

I wouldn't cite it in my own work yet, but I'd put it in front of a reading group for discussion.

Referee Report

4 major / 5 minor

Summary. The paper reports three case studies in which pre-known chemical information (fragment identity, connectivity, approximate bond lengths) is used to guide partial-structure R1 (pR1) and single-atom R1 (sR1) calculations for solving small-molecule crystal structures. Sample 1 is a light-atom-only structure assembled by first orienting and positioning four benzene-star fragments with pR1 and then completing the model with connectivity-guided sR1. Sample 2 uses the same pR1-plus-connectivity-guided-sR1 approach for a second light-atom structure. Sample 3 is a heavy-atom-containing structure assembled by first determining the heavy-atom substructure with normal sR1 and then completing the remaining side chains with connectivity-guided sR1. The paper concludes with two workflow recommendations: for light-atom structures, use normal sR1 for the framework and connectivity-guided sR1 to complete it, except at low resolution where pR1 is said to be necessary; for heavy-atom structures, use normal sR1 for the heavy-atom substructure and then connectivity-guided sR1.

Significance. If the proposed strategies are reliable, the paper offers practical guidance for using pR1/sR1 in small-molecule structure assembly: it demonstrates that pre-known fragments and connectivity can be exploited in a stepwise, planned manner, and it provides timing comparisons suggesting that connectivity-guided sR1 is faster than normal sR1 in the cases shown. The manuscript is transparent about its implementation choices, such as intensity sharpening and the specific grid parameters, and it explicitly names the central algorithmic assumption (deepest-hole hypothesis in SI S4). The main limitations are that the evidence consists of three single successful runs without deposited raw data or code, the success metric is purely positional, and the low-resolution branch of the central recommendation rests on a single unreported truncation experiment. These limitations currently make the generalized conclusions stronger than the evidence supports.

major comments (4)

[Section 5] The low-resolution branch of conclusion (1) is supported only by the sentence: 'For sample 1, when data resolution is truncated to 1.2 Å, the sR1 method alone can no longer solve the structure, but with the pR1 method’s help to orient and position four benzene-star fragments, the sR1 method can then complete the model.' No truncated-data R1-map values, orientation rankings, location-search rankings, timings, or success metrics are reported, and no failure statistics or additional truncation levels are shown. Because the paper states that pR1 is 'necessary' in this regime, this claim is currently a single anecdote. The authors should report the quantitative truncated experiment, including the rank of the true orientation and location for each fragment and the final model accuracy, and ideally provide more than one truncated-data trial.
[SI S4] The entire pR1 procedure rests on the 'general hypothesis' that the deepest hole of a pR1 map in a 6-dimensional orientation-location space determines the missing fragment, together with the further assumption that the 6-D search can be split into two independent 3-D searches. The paper does not test these assumptions in the case studies; for each sample the text merely states that the pR1 method 'correctly oriented' or 'correctly positioned' the fragment. Since the recommended strategies inherit these assumptions, the paper should report, for each fragment, the rank of the true orientation among the candidate orientations and the rank of the true location among the candidate locations. Without this evidence, the workflow recommendation cannot be distinguished from a heuristic that worked in the demonstrated examples.
[Section 4 and S3] The reported success metric counts atoms whose positions are within 0.5 Å regardless of atom type, yet sample 3's initial normal-sR1 run produced ten type misassignments (4 I as Mo, 3 Mo as I, one Mo as S, one S as I, one N as S) that were then 'corrected' manually. Because the model comparison in S3 explicitly disregards atom types, the claim that the strategy 'results in better quality of a model' is not supported with respect to chemical identity. The authors should report type-aware match counts and specify whether the type corrections are part of the algorithm or a manual crystallographer intervention, since that distinction materially affects the reproducibility of the claimed strategy.
[Sections 2, 4, 5] All conclusions are derived from a single successful run for each of three structures under one implementation, with no raw data or code deposited and timings reported as single measurements on one device. The paper generalizes to a 'usual strategy' and a 'correct strategy,' which requires stronger evidence than three anecdotal successes. The authors should either provide additional runs and ideally further test structures, or explicitly scope the conclusions to the demonstrated cases and present the work as a proof-of-concept rather than a validated general protocol.

minor comments (5)

[Introduction] There is a typographical error: 'the pR1are two new model-searching techniques' should read 'the pR1 are two new model-searching techniques'.
[SI S1] The radiation wavelengths are given as λ = 0.71073 nm and λ = 1.54178 nm, but Mo Kα is 0.71073 Å and Cu Kα is 1.54178 Å; as written these values are off by a factor of ten.
[Abstract and Section 5] The phrase 'and then uses the connectivity-guided sR1 method' in the Abstract and in Section 5 should be 'and then use' to agree with the subject 'strategy'.
[S3] The description of the model-comparison algorithm says the shift-and-overlap procedure is repeated for all atoms in both models, but it is not explicitly stated whether the inversion of model B is tried in both the original and shifted frames; please clarify.
[Figures 1–3] The figure captions summarize the steps, but the actual figures are not included in the text supplied for review; the published version should ensure that the step numbering in the figures matches the step labels used in the text.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the R1 objective is data-derived, the pre-known fragments are external chemical inputs, and the deepest-hole heuristic is tested against correct models rather than assumed as a theorem.

full rationale

The paper does not derive pR1 or sR1 from first principles; it applies the implementation from Zhang & Donahue (2024) and reports three case studies. The R1 objective is computed against measured intensities (sharpened via a Wilson B parameter), and the pre-known fragments, connectivity, and approximate bond lengths are external chemical knowledge, not quantities fitted from the diffraction data. The load-bearing heuristic in SI S4, 'The general hypothesis is that the deepest hole of a pR1 map in a 6-dimensional orientation-location space determines the orientation and location of a missing fragment,' is explicitly labeled a hypothesis, and the paper's final models are compared against the correct structures within a 0.5 Å tolerance. That comparison makes the hypothesis externally checkable rather than assumed by construction: if the deepest-hole search had selected wrong orientations or positions, the reported agreement with the correct model would not hold. The Section 5 low-resolution claim about sample 1 at 1.2 Å truncation is underreported, with no R1 values, orientation rankings, or failure statistics, but this is an evidence-completeness weakness, not circularity: the statement is not equivalent to its inputs by definition. The self-citations to Zhang & Donahue (2024) identify the software implementation and the inherited search heuristic; they do not smuggle in the conclusion or invoke a uniqueness theorem to forbid alternatives. No fitted parameter is renamed as a prediction, and no known result is repackaged as a derivation. The central recommendation is a workflow suggestion grounded in timed case studies, and although the number of examples is small, the reasoning chain does not reduce to its own inputs.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

The central claim rests on a handful of hand-chosen algorithmic parameters and on the inherited deepest-hole heuristic. No new physical entities are introduced.

free parameters (6)

Wilson B parameter = not reported; estimated from slope b of Wilson plot
Estimated from the intensity falloff by linear regression in SI S2 and used to sharpen intensities. It is a data-fitted parameter and affects ghost-atom behavior.
Translation grid step = 0.4 angstroms
Chosen by hand in SI S4 for the coarse orientation-location search. It sets the resolution of the search and therefore influences success and runtime.
Rotation angle step = 5 degrees
Chosen by hand in SI S4 for the coarse rotational search. It affects whether the true orientation is found and how long the search takes.
Connectivity shell tolerance = 0.5 angstroms (values 0.3, 0.4, or 0.5 allowed)
Used in the connectivity-guided sR1 method to limit the search for a bonded atom to a spherical shell. The choice affects which atoms can be found.
Model match tolerance = 0.5 angstroms
Used in SI S3 to decide whether an atom in one model matches an atom in another. The reported 'all atoms within 0.5 angstroms' success statement depends on this threshold.
Approximate bond lengths = 1.39 angstroms C-C, 2.0 angstroms S-S, 1.8 angstroms S-C
Supplied as pre-known chemical inputs for the connectivity-guided sR1 searches. They are not fitted to the target data, but the search shells are centered on them.

assumptions (5)

domain assumption The deepest hole of a pR1 map in orientation-location space determines the true fragment orientation and location.
Invoked in SI S4 and in Zhang and Donahue (2024). This unproven heuristic underlies all pR1 and sR1 placements.
domain assumption The normal sR1 method, with ghost-atom clustering and triangular-bonding exclusion rules, correctly locates single atoms.
Carried over from the prior method and used in all three samples. No proof is given in this paper.
domain assumption Wilson sharpening with a U=0 atomic model matches the data better and reduces ghost atoms.
Section 1 states this as a finding without independent demonstration, and it motivates the sharpening step used throughout.
domain assumption The known expected structures and fragment geometries used as pre-knowledge are correct.
Sections 2 to 4 rely on this. If the pre-knowledge is wrong, the assembled model inherits the error. This is an application condition rather than a mathematical axiom.
standard math Crystallographic standard transformations and linear algebra, such as cell-to-Cartesian conversion and rotation matrices, are correct.
SI S4 uses these standard formulas without proof.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Case Studies of Using the Partial-Structure R1 and the Single-Atom R1 to Assemble Small-Molecule Crystal Structures." pith.science (2026). https://pith.science/paper/D55JG56K

@misc{pith2026241215284,
  author       = {Pith},
  title        = {Pith review of: Case Studies of Using the Partial-Structure R1 and the Single-Atom R1 to Assemble Small-Molecule Crystal Structures},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/D55JG56K}},
  note         = {Machine review of arXiv:2412.15284}
}

read the original abstract

This paper demonstrates how pre-knowledge of a crystal structure, including the constituent fragments, the atomic connectivity, and the approximate bond lengths, etc., can be utilized in the partial-structure R1 (pR1) and the single-atom R1 (sR1) calculations. It has been shown that taking advantage of pre-known information the calculations can proceed in an orderly and well-planned manner. Furthermore, in the case of the sR1 calculation, computer time can also be greatly reduced. Because the pR1 calculation is more time-consuming than the sR1 calculation, when there is a choice between the pR1 and the sR1, the former should be avoided. So, the usual strategy of assembling a light-atom-only structure should start with the normal sR1 method to determine a basic framework of the structure, and then uses the connectivity-guided sR1 method to complete the model. However, when the data resolution is low, the first step is necessary to use the pR1 method to assemble the known fragments. For a heavy-atom-containing structure, the correct strategy starts with the normal sR1 method to determine the heavy-atom substructure, and then uses the connectivity-guided sR1 to complete the model.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

1 extracted references

[1]

& Donahue, J

Zhang, X. & Donahue, J. P. (2024). Acta Cryst. A80, 237-248

2024

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.

[1] [1]

& Donahue, J

Zhang, X. & Donahue, J. P. (2024). Acta Cryst. A80, 237-248

2024