{"id":"05051d82-bdad-4cb8-8d55-1a90e6895a12","arxiv_id":"2412.14034","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The deepest hole of a partial-structure R1 map is used to orient and locate missing molecular fragments, demonstrated on four small-molecule structures.","lead":"A crystallography method uses the deepest holes of a partial-structure R1 map to find the orientation and position of missing molecular fragments. The paper reports successful placements on four small-molecule crystal structures, with several manual corrections and tuning steps.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sample 4 contradicts the central claim: the lowest-R1 'deepest hole' for ring 25 is wrong, and the correct solution is chosen manually using chemical knowledge (Section 10).","rationale":"The reader's weakest_assumption is the core approximation in Eq. (1), specifically the omitted cross terms between the modeled block and the unmodeled tail. That is a real and important concern, but it is not the most load-bearing one for the paper's central claim. The paper itself provides a direct counterexample to the claim that the deepest pR1 hole determines the correct solution: in Section 10, the lowest-R1 orientation for ring 25 is wrong, and the authors adopt the second-lowest orientation based on chemical bonding considerations. In Section 9, ring 0 must be deleted and re-discovered after visual inspection. Therefore, even if Eq. (1) is accepted as an adequate approximation, the scoring function alone does not select the correct orientation and location in the reported tests. The successes are achieved with human pattern recognition and chemical judgment as an essential part of the pipeline. This is not an external consensus disagreement; it is an internal inconsistency with the strong claim in the abstract. The paper has genuine strengths: it clearly explains the two-step orientation/location search, explicitly retains the unmodeled-atom baseline in the pR1 target, and reports four structures that end within 0.5 Angstrom of the known answers after manual curation. Because the paper is a demonstration rather than a formal proof, a conditional verdict is appropriate, but the central claim should be restated as 'pR1 provides useful candidate solutions that require chemical filtering', not 'the deepest hole determines the orientation and location'. The reader's conditional verdict already captures the need for caution, so I do not recommend changing the verdict; rather, the concern sharpens the reason for the condition.","tokens_in":9938,"tokens_out":6055,"duration_ms":58546,"concrete_test":"Re-run the sample 4 pipeline exactly as described through ring 25, with no human override: after placing rings 0-24 (using the same manual deletions already reported), evaluate all candidate orientations for ring 25 and select the one with lowest R1 that passes the ghost/triangular-bond rules. If the selected ring is the isolated one from orientation 16, the central claim is falsified under the paper's own scoring. If a modified automatic rule (e.g., requiring a chemically sensible connection when alternatives are within ~0.001 in R1) is proposed, test that rule on all four samples and report how many manual interventions remain.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central hypothesis is that the deepest hole of a pR1 map determines the orientation and location of a missing fragment. The paper's own Section 10 provides a counterexample. After rings 0-24 are placed, ring 25 is determined by trying all candidate orientations: orientation 16 gives the lowest R1 (0.5827) but places the ring isolated, with no chemically meaningful bonds; orientation 2 gives the second-lowest R1 (0.5836) and is the correct solution. The authors explicitly override the lowest-R1 choice ('the result of this orientation should be adopted') because chemical knowledge says so. Thus, under the algorithm's own scoring rule, the deepest hole does not determine the correct orientation and location in one of the four test cases. Similarly, in Section 9, benzene ring 0 is deleted and re-discovered after visual inspection, so the 'strict calculation' reported in Section 8 as 'a complete success' actually required manual correction. These internal observations are more decisive than the Eq. (1) approximation concern: even granting Eq. (1), the R1 ranking is not a reliable selector of the true solution. The tiny R1 difference (0.0009) suggests the target is nearly flat near the correct solution. The paper would need to weaken the claim to 'pR1 generates candidate solutions that, after chemical filtering, can complete structures', and provide an automatic decision rule that reproduces the reported successes without human judgment.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Zhang presents pR1 as a generalization of sR1, hypothesizing that the deepest hole of a pR1 map determines the orientation and location of a missing fragment. The paper describes a two-step protocol: first detect candidate orientations from holes in the pR1 map of a free-standing fragment in 3D rotation space; then, for each trial orientation, locate the fragment by the deepest hole (or global minimum) of a pR1 map in translation space, ranking orientation-location candidates by R1. The approach is demonstrated on four small-molecule crystal structures, and final models match the known structures within 0.5 Å. A spherical 'completely disoriented fragment' model is introduced to predict possible locations and speed up the search. The central claim is that pR1 is a successful MR target.","tokens_in":10261,"tokens_out":5062,"duration_ms":41489,"significance":"If the central claim held, the paper would challenge the longstanding view that R1-type targets are unsuitable for molecular replacement, and it would offer a conceptually simple alternative to maximum-likelihood targets. The strengths are the clear falsifiable formulation and the use of externally known crystal structures as ground truth for validation. However, the evidence is limited to four small-molecule examples, and the paper explicitly reports manual interventions and an ad hoc n=1000 scaling. The deepest-hole criterion is contradicted in Section 10 for ring 25. The demonstrated pipeline is therefore not yet an automatic MR method, but it provides candidate solutions that, after chemical filtering, complete the tested structures. The significance is moderate and conditional: the results are promising enough to warrant revision, but the present claims outrun the evidence.","major_comments":[{"comment":"The central hypothesis that the deepest hole (lowest R1) determines the missing fragment is contradicted by the paper's own sample 4 result. Section 10 reports that for ring 25, orientation 16 yields R1=0.5827 and an isolated ring, while orientation 2, with R1=0.5836, is correct and is adopted manually. Thus the lowest-R1 selector fails in one of four test cases; the claim must be weakened to 'pR1 generates candidate placements that may require chemical filtering,' and an automatic decision rule must be supplied.","section":"Abstract, Section 3, Section 10"},{"comment":"The reported successes rely on human judgment at several steps: ring 0 is placed at (0.3,0.3,0.3) by hand in Section 9 (and similarly molecule 0 in Section 6 and the first Cu atom in Section 7), and Section 9 states that benzene ring 0 'is deleted and re-discovered' after visual inspection, contradicting the 'complete success' claimed in Section 8 for the strict calculation. Section 10 likewise requires visual identification and deletion of 12 suspected bad rings. Since no automated decision rule is provided for these steps, the method as demonstrated is not a complete MR algorithm, and the 0.5 Å validation applies to the final manually adjusted model.","section":"Sections 6, 7, 8, 9, 10"},{"comment":"The core approximation in Eq. (1) drops all cross terms between the known atoms/current fragment (groups 1 and 2) and the tail of other missing atoms (group 3). This approximation is inherited from Zhang & Donahue (2024) and is not rederived or tested here. Because every pR1 hole and every placement in Sections 6-10 is computed under this approximation, its adequacy is load-bearing: if the omitted cross terms are substantial, the deepest hole need not coincide with the true fragment position. The authors should at least provide a numerical test of the approximation, e.g., comparing Eq. (1) with exact |F|² for the known structures.","section":"Section 2, Eq. (1)"},{"comment":"The 'completely disoriented fragment' model introduces a parameter n that is first defined as n≡1 (Eq. 3) and later arbitrarily set to 1000 for two of the fragment types 'to resume prediction power.' This is a free parameter tuned to the examples, with no physical or statistical justification. The claim that the model 'can predict the locations' of missing fragments is accordingly not robust; the authors should either derive n from first principles or abandon the disoriented-model prediction as evidence.","section":"Section 8, Eqs. (2)-(3)"}],"minor_comments":[{"comment":"The symbol R1 is used without a formal definition; the reader is not told how R1 is computed from Fc² (e.g., R1 = Σ||Fo|-|Fc||/Σ|Fo|, or a variant). Please define R1 and pR1 explicitly.","section":"Section 2"},{"comment":"The procedure for identifying 'holes' in a pR1 map (e.g., local minima, threshold criteria) is not specified; the reported counts (14340, 780, 4176, 432 holes) cannot be reproduced without this information.","section":"Sections 6-10"},{"comment":"The relationship between Eq. (2), Eq. (3) and the nomenclature 'completely disoriented model' is not fully explained; in particular, the derivation of the sinc term G(4πsr_i) should be sketched or referenced.","section":"Section 8"},{"comment":"The 'clustering ghost atoms' and 'triangular bonding' rules are referenced only to the previous paper; a brief summary is needed for the present paper to be self-contained.","section":"Section 3"},{"comment":"The rotation and coordinate-conversion matrices would benefit from a worked example and from verification of the sign conventions; the current presentation is hard to check.","section":"Supporting information S3/S4"},{"comment":"'General MR calculations' is overstated given the small-molecule scope; please qualify the claim to the demonstrated regime.","section":"Section 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is heavily self-referential, relying on the author's prior sR1 paper for the core equations and filtering rules; the four examples all come from the author's own crystallography, and no code or data are released. Given the journal's scope, the paper is more of a research report than a finished methods paper; the counterexample in Section 10 and the manual interventions should be addressed before publication. Also note that the 'n=1000' parameter is not justified and appears data-dependent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is an honest, clearly written proof-of-concept for a partial-structure R1 target in molecular replacement, but the paper's central claim — that the deepest hole of a pR1 map selects the correct orientation and location — is contradicted by its own Section 10. For ring 25 in sample 4, orientation 16 gives R1 = 0.5827 and a chemically meaningless placement; the correct orientation 2 has R1 = 0.5836 and is adopted only because the chemist sees it makes a proper bond. That is a 0.0009 R1 difference, so the target is nearly flat near the correct answer. Similarly, sample 3 required deleting and re-discovering ring 0 after visual inspection, so the 'complete success' of the strict calculation is not fully automatic.\n\nWhat the paper does well: the idea of retaining the unmodeled-atom baseline term (the tail fj+1^2+...+fN^2) in an R1 target is conceptually interesting and clearly explained in Section 12. The two-step search (orientation from a free-standing fragment, then location from a pR1 map) is a sensible decomposition and is described with enough detail to reimplement. All four final models match the known structures to 0.5 Å, which is a real result, albeit on small molecules.\n\nThe soft spots are not minor. The n=1000 scaling in the disoriented-fragment model is ad hoc and unexplained; a spherical-average scattering factor times 1000 has no physical basis as stated. The manual first placements at (0.3,0.3,0.3) and the lucky guesses about which orientation labels correspond to which fragments in samples 1 and 2 mean the method's success on those samples is partly handcrafted. And Equation (1), the core approximation, drops cross terms between the modeled fragment and the tail; this is inherited from the earlier sR1 paper without sensitivity analysis here.\n\nThe honest reading is that pR1 produces candidate placements that, after chemical filtering, can lead to a correct structure. That is a weaker claim than the abstract's, and it is the claim the evidence actually supports. The paper would be improved by an automatic decision rule that reproduces the successes without human judgment.\n\nThis is a paper for crystallographers who think about MR targets, not for users looking for a turnkey tool. It deserves a serious referee: the concept is new, the writing is clear, and the counterexample is instructive rather than hidden. I would send it to review with a request to soften the claim and to address the manual steps explicitly.","headline":"The deepest-hole hypothesis is not borne out by the paper's own Section 10, where the lowest-R1 solution is discarded for chemical reasons; still, the pR1 baseline idea is worth a referee's time.","tokens_in":10774,"tokens_out":2678,"would_cite":false,"duration_ms":22115,"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 crystallographic target called pR1 locates missing molecular fragments by the deepest hole of its map, completing four test structures to within 0.5 Å.","keywords":["partial structure R1","molecular replacement","single atom R1","pR1 map","completely disoriented fragment model","R1 target","crystallography","structure solution"],"falsifier":"For a known structure with one fragment removed, compute the pR1 map from Equation (1) and check whether the deepest hole coincides with the missing fragment; finding a structure where the deepest hole consistently marks a different position while a shallower hole marks the true fragment would falsify the deepest-hole hypothesis.","tokens_in":9705,"feed_emoji":"🔬","tokens_out":7077,"duration_ms":56605,"temperature":0.7,"pith_summary":"This paper proposes that a crystallographic target called partial structure R1 (pR1) can solve molecular replacement: the deepest hole in a pR1 map identifies the orientation and location of a missing molecular fragment. The idea generalizes the earlier single-atom R1 method by treating any known geometrical fragment as the current model, and it splits the six-dimensional search into separate orientation and location stages. The paper tests the method on four small-molecule structures, including up to 26 benzene rings, and reports that every final model matches the correct structure to within 0.5 Å. The paper argues that pR1 works where the traditional R1 target fails because it retains a baseline term for all atoms still outside the model.","feed_headline":"Deepest hole of a pR1 map finds missing molecular fragments","feed_subtitle":"A partial-structure R1 map places missing fragments in test crystals to within 0.5 Å.","key_machinery":"The central object is the pR1 map, a function that gives the R1 agreement value for each trial orientation and position of a missing fragment. It is computed from Equation (1), the core approximation: $F_c^2(hkl)$ equals the squared cosine/sine sums over atoms 1 to j (the known atoms plus the trial fragment) plus the tail $f_{j+1}^2(hkl)+\\dots+f_N^2(hkl)$ of squared scattering factors for the atoms j+1 to N that are not yet placed. The load-bearing mechanism is the deepest-hole hypothesis: the orientation and location that make the pR1 map deepest are the true orientation and location of the missing fragment. To keep the calculation practical, the paper divides the search into a free-standing-fragment orientation scan, where holes in a 3-D orientation space provide candidate orientations, and a location scan in a 3-D location space, with the global minimum found by coarse grid and five-step refinement.","core_discovery":"The central claim is that the deepest hole of a pR1 map determines the orientation and location of a missing fragment. pR1 is defined by the same core approximation used for sR1: the expected squared structure factor is the squared sum of cosine and sine contributions from the known atoms and the trial fragment, plus a tail of squared scattering factors for the remaining unmodeled atoms. According to the paper, keeping this tail accounts for the baseline effect of the other missing atoms, and dropping it is what makes the ordinary R1 target weak; maximum-likelihood targets such as LLGI mitigate the same problem by taking a gain. With this target, the paper determines two S2O2C12 molecules, four N4C9 and four PF6 fragments in sample 2, and seven or twenty-six benzene rings in samples 3 and 4, reporting final models with all atoms within 0.5 Å of the correct coordinates.","pith_inferences":["If the deepest-hole hypothesis is general, pR1 maps could also serve as an incomplete-model diagnostic: the deepest hole should point at whatever part of the structure is missing, not only in molecular replacement.","The paper's comparison with maximum-likelihood targets is left qualitative, so a natural testable extension is a head-to-head pR1 versus LLGI search on the same data to see where the baseline term most changes the ranking.","The complete-disorientation model's failure for S2O2C12 and N4C9 at n=1 and its recovery at n=1000 suggests an unexplored tuning rule for this predictor; understanding when large n is needed could extend the method to more rigid fragments.","Because the method is demonstrated only on four small structures, macromolecular application is open; the orientation-equivalence test that ignores atom types would likely need type awareness for proteins."],"forward_implications":["pR1 can be used as a molecular-replacement target for small-molecule crystals: in all four test structures, the final model is within 0.5 Å of the correct structure.","Retaining the tail term in Equation (1) is the paper's explanation for why pR1 succeeds where the traditional R1 target fails.","A completely disoriented, spherically averaged fragment model can predict the locations of missing fragments, shortening the search enough to handle structures with many fragments such as 26 benzene rings.","Errors in molecular replacement can be found and fixed by deleting fragments that bond incorrectly and resuming the calculation, as done for samples 3 and 4."],"supporting_citations":[{"why":"Introduces the sR1 concept and the core approximation, Equation (1), that pR1 inherits.","marker":"Zhang & Donahue, 2024"},{"why":"Defines the molecular-replacement problem and the rotation-function/translation-function framework that pR1 targets replace.","marker":"Rossmann & Blow, 1962"},{"why":"States the general belief that R1 is not a good MR target, the baseline the paper argues pR1 overturns.","marker":"McCoy, 2017"},{"why":"Describes the LLGI likelihood target that the paper compares with pR1 in explaining why pR1 works.","marker":"Read & McCoy, 2016"},{"why":"An early use of the traditional R1 target in translation searches, providing the contrast case for pR1's retained tail.","marker":"Eventoff et al., 1975"}],"fun_headline_variants":["pR1 map holes reveal missing fragment orientation","Deepest pR1 hole determines fragment placement","Partial R1 map finds missing molecular pieces","Molecular replacement with partial R1 maps","pR1 holes locate fragments to 0.5 Å"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole method rests on the core approximation of Equation (1), which assumes that the modeled atoms and the unmodeled tail contribute to the squared structure factor as a squared sum plus a sum of squared scattering factors, with all cross terms between the two groups omitted.","fun_headline_variants_meta":{"raw":{"variants":["pR1 map holes reveal missing fragment orientation","Deepest pR1 hole determines fragment placement","Partial R1 map finds missing molecular pieces","Molecular replacement with partial R1 maps","pR1 holes locate fragments to 0.5 Å"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000275,"raw_usage":{"total_tokens":1683,"prompt_tokens":1022,"completion_tokens":661,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":590}},"tokens_in":638,"tokens_out":661,"duration_ms":5832,"temperature":1.0,"reasoning_tokens":590,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:32:51.854531+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a known structure with one fragment removed, compute the pR1 map from Equation (1) and check whether the deepest hole coincides with the missing fragment; finding a structure where the deepest hole consistently marks a different position while a shallower hole marks the true fragment would falsify the deepest-hole hypothesis.","supporting_citations":[{"cited_title":"& Donahue, J","cited_arxiv_id":null,"evidence_quote":"Introduces the sR1 concept and the core approximation, Equation (1), that pR1 inherits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the molecular-replacement problem and the rotation-function/translation-function framework that pR1 targets replace."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the LLGI likelihood target that the paper compares with pR1 in explaining why pR1 works."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"An early use of the traditional R1 target in translation searches, providing the contrast case for pR1's retained tail."}],"review_version":1}