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

Abstract

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.

Accepted for publication in Phys Rev B

References in corpus (12)

• The density-matrix renormalization group in the age of matrix product states
• Matrix Product Density Operators: Simulation of finite-T and dissipative systems
• Matrix product states represent ground states faithfully
• A spin-adapted Density Matrix Renormalization Group algorithm for quantum chemistry
• Simulating Strongly Correlated Quantum Systems with Tree Tensor Networks
• CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry
• Applying matrix product operators to model systems with long-range interactions
• Density matrix renormalization group approach for many-body open quantum systems
• Efficient classical simulation of the semi-classical Quantum Fourier Transform
• Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations
• A quantum mechanical NMR simulation algorithm for protein-scale spin systems
• Exploiting Matrix Symmetries and Physical Symmetries in Matrix Product States and Tensor Trains

Cited by in corpus (14)

• Training Schrödinger's cat: quantum optimal control
• The density matrix renormalization group for ab initio quantum chemistry
• Tensor product methods and entanglement optimization for ab initio quantum chemistry
• Tensor Networks for Dimensionality Reduction and Large-Scale Optimizations. Part 2 Applications and Future Perspectives
• Parallel cross interpolation for high-precision calculation of high-dimensional integrals
• Imaginary-time matrix product state impurity solver for dynamical mean-field theory
• Fokker-Planck formalism in magnetic resonance simulations
• Tensor Networks for Big Data Analytics and Large-Scale Optimization Problems
• Excited-state DMRG made simple with FEAST
• Digital quantum simulation of NMR experiments
• Hybrid quantum-classical method for simulating high-temperature dynamics of nuclear spins in solids
• Free induction decays in nuclear spin-1/2 lattices with small number of interacting neighbors: the cases of silicon and fluorapatite
• Tensor product approach to quantum control
• Parallel algorithms for computing the tensor-train decomposition