The general equation of δ direct methods and the novel <i>SMAR</i> algorithm residuals using the absolute value of ρ and the zero conversion of negative ripples

The general equation δ<sub>M</sub>(r) = ρ(r) + g(r) of the δ direct methods (δ-GEQ) is established which, when expressed in the form δ<sub>M</sub>(r) - ρ(r) = g(r), is used in the SMAR phasing algorithm [Rius (2020). Acta Cryst A76, 489-493]. It is shown that SMAR is based on the alternating minimization of the two residuals R<sub>ρ</sub>(χ) = ∫<sub>V</sub> [ρ(χ) - ρ(Φ)s<sub>ρ</sub>]<sup>2</sup> dV and R<sub>δ</sub>(Φ) = ∫<sub>V</sub> m<sub>ρ</sub>[δ<sub>M</sub>(χ) - ρ(Φ)s<sub>ρ</sub>]<sup>2</sup> dV in each iteration of the algorithm by maximizing the respective S<sub>ρ</sub>(Φ) and S<sub>δ</sub>(Φ) sum functions. While R<sub>ρ</sub>(χ) converges to zero, R<sub>δ</sub>(Φ) converges, as predicted by the theory, to a positive quantity. These two independent residuals combine δ<sub>M</sub> and ρ each with |ρ| while keeping the same unknowns, leading to overdetermination for diffraction data extending to atomic resolution. At the beginning of a SMAR phase refinement, the zero part of the m<sub>ρ</sub> mask [resulting from the zero conversion of the slightly negative ρ(Φ) values] occupies ∼50% of the unit-cell volume and increases by ∼5% when convergence is reached. The effects on the residuals of the two SMAR phase refinement modes, i.e. only using density functions (slow mode) supplemented by atomic constraints (fast mode), are discussed in detail. Due to its architecture, the SMAR algorithm is particularly well suited for Deep Learning. Another way of using δ-GEQ is by solving it in the form ρ(r) = δ<sub>M</sub>(r) - g(r), which provides a simple new derivation of the already known δ<sub>M</sub> tangent formula, the core of the δ recycling phasing algorithm [Rius (2012). Acta Cryst. A68, 399-400]. The nomenclature used here is: (i) Φ is the set of φ structure factor phases of ρ to be refined; (ii) δ<sub>M</sub>(χ) = FT<sup>-1</sup>{c(|E| - 〈|E|〉)×exp(iα)} with χ = {α}, the set of phases of |ρ| and c = scaling constant; (iii) m<sub>ρ</sub> = mask, being either 0 or 1; s<sub>ρ</sub> is 1 or -1 depending on whether ρ(Φ) is positive or negative.

Authors

Institutions

Publication Details

Journal
Acta Crystallographica Section A Foundations and Advances
Published
2024-11-19
DOI
https://doi.org/10.1107/s2053273324009628
Primary Topic
X-ray Diffraction in Crystallography
Type
article
Field-Weighted Citation Impact
0.00

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

The general equation of δ direct methods and the novel SMAR algorithm residuals using the absolute value of ρ and the zero conversion of negative ripples

Jordi Rius
Acta Crystallographica Section A Foundations and Advances
X-ray Diffraction in Crystallography
article

The general equation of δ direct methods and the novel SMAR algorithm residuals using the absolute value of ρ and the zero conversion of negative ripples

Jordi Rius
article en

Abstract

The general equation δ<sub>M</sub>(r) = ρ(r) + g(r) of the δ direct methods (δ-GEQ) is established which, when expressed in the form δ<sub>M</sub>(r) - ρ(r) = g(r), is used in the SMAR phasing algorithm [Rius (2020). Acta Cryst A76, 489-493]. It is shown that SMAR is based on the alternating minimization of the two residuals R<sub>ρ</sub>(χ) = ∫<sub>V</sub> [ρ(χ) - ρ(Φ)s<sub>ρ</sub>]<sup>2</sup> dV and R<sub>δ</sub>(Φ) = ∫<sub>V</sub> m<sub>ρ</sub>[δ<sub>M</sub>(χ) - ρ(Φ)s<sub>ρ</sub>]<sup>2</sup> dV in each iteration of the algorithm by maximizing the respective S<sub>ρ</sub>(Φ) and S<sub>δ</sub>(Φ) sum functions. While R<sub>ρ</sub>(χ) converges to zero, R<sub>δ</sub>(Φ) converges, as predicted by the theory, to a positive quantity. These two independent residuals combine δ<sub>M</sub> and ρ each with |ρ| while keeping the same unknowns, leading to overdetermination for diffraction data extending to atomic resolution. At the beginning of a SMAR phase refinement, the zero part of the m<sub>ρ</sub> mask [resulting from the zero conversion of the slightly negative ρ(Φ) values] occupies ∼50% of the unit-cell volume and increases by ∼5% when convergence is reached. The effects on the residuals of the two SMAR phase refinement modes, i.e. only using density functions (slow mode) supplemented by atomic constraints (fast mode), are discussed in detail. Due to its architecture, the SMAR algorithm is particularly well suited for Deep Learning. Another way of using δ-GEQ is by solving it in the form ρ(r) = δ<sub>M</sub>(r) - g(r), which provides a simple new derivation of the already known δ<sub>M</sub> tangent formula, the core of the δ recycling phasing algorithm [Rius (2012). Acta Cryst. A68, 399-400]. The nomenclature used here is: (i) Φ is the set of φ structure factor phases of ρ to be refined; (ii) δ<sub>M</sub>(χ) = FT<sup>-1</sup>{c(|E| - 〈|E|〉)×exp(iα)} with χ = {α}, the set of phases of |ρ| and c = scaling constant; (iii) m<sub>ρ</sub> = mask, being either 0 or 1; s<sub>ρ</sub> is 1 or -1 depending on whether ρ(Φ) is positive or negative.

Acta Crystallographica Section A Foundations and AdvancesVol. 81(1)
Institut de Ciència de Materials de Barcelona (ES)
Consejo Superior de Investigaciones Científicas, Agencia Estatal de Investigación
Openalex Percentile: Top 87%
X-ray Diffraction in Crystallography
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.