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

Abstract

The general equation δ_M(r) = ρ(r) + g(r) of the δ direct methods (δ-GEQ) is established which, when expressed in the form δ_M(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_ρ(χ) = ∫_V [ρ(χ) - ρ(Φ)s_ρ]² dV and R_δ(Φ) = ∫_V m_ρ[δ_M(χ) - ρ(Φ)s_ρ]² dV in each iteration of the algorithm by maximizing the respective S_ρ(Φ) and S_δ(Φ) sum functions. While R_ρ(χ) converges to zero, R_δ(Φ) converges, as predicted by the theory, to a positive quantity. These two independent residuals combine δ_M 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_ρ 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) = δ_M(r) - g(r), which provides a simple new derivation of the already known δ_M 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) δ_M(χ) = FT^-1{c(|E| - 〈|E|〉)×exp(iα)} with χ = {α}, the set of phases of |ρ| and c = scaling constant; (iii) m_ρ = mask, being either 0 or 1; s_ρ is 1 or -1 depending on whether ρ(Φ) is positive or negative.