Source author record

D. V. Savostyanov

D. V. Savostyanov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA physics.chem-ph quant-ph Data Structures and Algorithms math.CA Numerical Analysis

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2014arXiv

A quantum mechanical NMR simulation algorithm for protein-scale spin systems

Nuclear magnetic resonance spectroscopy is one of the few remaining areas of physical chemistry for which polynomially scaling simulation methods have not so far been available. Here, we report such a method and illustrate its performance by simulating common 2D and 3D liquid state NMR experiments (including accurate description of spin relaxation processes) on isotopically enriched human ubiquitin - a protein containing over a thousand nuclear spins forming an irregular polycyclic three-dimensional coupling lattice. The algorithm uses careful tailoring of the density operator space to only include nuclear spin states that are populated to a significant extent. The reduced state space is generated by analyzing spin connectivity and decoherence properties: rapidly relaxing states as well as correlations between topologically remote spins are dropped from the basis set. In the examples provided, the resulting reduction in the quantum mechanical simulation time is by many orders of magnitude.

preprint2014arXiv

Exact NMR simulation of protein-size spin systems using tensor train formalism

We introduce a new method, based on alternating optimization, for compact representation of spin Hamiltonians and solution of linear systems of algebraic equations in the tensor train format. We demonstrate the method's utility by simulating, without approximations, a 15N NMR spectrum of ubiquitin --- a protein containing several hundred interacting nuclear spins. Existing simulation algorithms for the spin system and the NMR experiment in question either require significant approximations or scale exponentially with the spin system size. We compare the proposed method to the Spinach package that uses heuristic restricted state space techniques to achieve polynomial complexity scaling. When the spin system topology is close to a linear chain (e.g. for the backbone of a protein), the tensor train representation is more compact and can be computed faster than the sparse representation using restricted state spaces.

preprint2013arXiv

Grid-free powder averages: on the applications of the Fokker-Planck equation to solid state NMR

We demonstrate that Fokker-Planck equations in which spatial coordinates are treated on the same conceptual level as spin coordinates yield a convenient formalism for treating magic angle spinning NMR experiments. In particular, time dependence disappears from the background Hamiltonian (sample spinning is treated as an interaction), spherical quadrature grids are avoided completely (coordinate distributions are a part of the formalism) and relaxation theory with any linear diffusion operator is easily adopted from the Stochastic Liouville Equation theory. The proposed formalism contains Floquet theory as a special case. The elimination of the spherical averaging grid comes at the cost of increased matrix dimensions, but we show that this can be mitigated by the use of state space restriction and tensor train techniques. It is also demonstrated that low correlation order basis sets apparently give accurate answers in powder-averaged MAS simulations, meaning that polynomially scaling simulation algorithms do exist for a large class of solid state NMR experiments.

preprint2013arXiv

On the decay of elements of inverse triangular Toeplitz matrix

We consider half-infinite triangular Toeplitz matrices with slow decay of the elements and prove under a monotonicity condition that elements of the inverse matrix, as well as elements of the fundamental matrix, decay to zero. We also provide a quantitative description of the decay of the fundamental matrix in terms of p-norms. Finally, we prove that for matrices with slow log-convex decay the inverse matrix has fast decay, i.e. is bounded. The results are compared with the classical results of Jaffard and Veccio and illustrated by numerical example.

preprint2010arXiv

Wedderburn rank reduction and Krylov subspace method for tensor approximation. Part 1: Tucker case

New algorithms are proposed for the Tucker approximation of a 3-tensor, that access it using only the tensor-by-vector-by-vector multiplication subroutine. In the matrix case, Krylov methods are methods of choice to approximate the dominant column and row subspaces of a sparse or structured matrix given through the matrix-by-vector multiplication subroutine. Using the Wedderburn rank reduction formula, we propose an algorithm of matrix approximation that computes Krylov subspaces and allows generalization to the tensor case. Several variants of proposed tensor algorithms differ by pivoting strategies, overall cost and quality of approximation. By convincing numerical experiments we show that the proposed methods are faster and more accurate than the minimal Krylov recursion, proposed recently by Elden and Savas.