Embedding with a rigid substructure

From MaRDI portal



Abstract: This paper presents a new distance geometry algorithm for calculating atomic coordinates from estimates of the interatomic distances, which maintains the positions of the atoms in a known rigid substructure. Given an Mimes3 matrix of coordinates for the rigid substructure mathbfX, this problem consists of finding the Nimes3 matrix mathbfY that yields of global minimum of the so-called STRAIN, i.e. [ min_{mathbf Y} left| �egin{bmatrix} mathbf{XX}^ op & mathbf{XY}^ op \ mathbf{YX}^ op & mathbf{YY}^ op end{bmatrix} ,-, �egin{bmatrix} mathbf A & mathbf B \ mathbf B^ op & mathbf C end{bmatrix} ight|_{mathsf F}^2 ~, ] where mathbfA=mathbfXXop , and mathbfB,mathbfC are matrices of inner products calculated from the estimated distances. The vanishing of the gradient of the STRAIN is shown to be equivalent to a system of only six nonlinear equations in six unknowns for the inertial tensor associated with the solution Y . The entire solution space is characterized in terms of the geometry of the intersection curves between the unit sphere and certain variable ellipsoids. Upon deriving tight bilateral bounds on the moments of inertia of any possible solution, we construct a search procedure that reliably locates the global minimum. The effectiveness of this method is demonstrated on realistic simulated and chemical test problems.











This page was built for publication: Embedding with a rigid substructure

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1381358)