Source author record

Ilya Kuprov

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

23 published item(s)

preprint2022arXiv

Neural networks in pulsed dipolar spectroscopy: a practical guide

This is a methodological guide to the use of deep neural networks in the processing of pulsed dipolar spectroscopy (PDS) data encountered in structural biology, organic photovoltaics, photosynthesis research, and other domains featuring long-lived radical pairs and paramagnetic metal ions. PDS uses distance dependence of magnetic dipolar interactions; measuring a single well-defined distance is straightforward, but extracting distance distributions is a hard and mathematically ill-posed problem requiring careful regularisation and background fitting. Neural networks do this exceptionally well, but their "robust black box" reputation hides the complexity of their design and training - particularly when the training dataset is effectively infinite. The objective of this paper is to give insight into training against simulated databases, to discuss network architecture choices, to describe options for handling DEER (double electron-electron resonance) and RIDME (relaxation-induced dipolar modulation enhancement) experiments, and to provide a practical data processing flowchart.

preprint2021arXiv

J-driven Dynamic Nuclear Polarization for sensitizing high field solution state NMR

Dynamic nuclear polarization is widely used to enhance the sensitivity of nuclear magnetic resonance. It is presently the method of choice for enhancing the sensitivity of high-field solid state NMR experiments performed under cryogenic conditions. In liquids, however, its efficiency decays rapidly with magnetic field Bo and with the rotational correlation time, leading to negligible enhancements in mid- and high-field solution-state NMR experiments for all but exceptional cases. This study discusses a potential solution to this fundamental sensitivity problem, that relies on biradical species possessing inter-electron exchange couplings, that are on the order of the electron Larmor frequency.

preprint2020arXiv

Bi-selective pulses for large-area atom interferometry

We present designs for the augmentation 'mirror' pulses of large-momentum-transfer atom interferometers that maintain their fidelity as the wavepacket momentum difference is increased. These bi-selective pulses, tailored using optimal control methods to the evolving bi-modal momentum distribution, should allow greater interferometer areas and hence increased inertial measurement sensitivity, without requiring elevated Rabi frequencies or extended frequency chirps. Using an experimentally validated model, we have simulated the application of our pulse designs to large-momentum-transfer atom interferometry using stimulated Raman transitions in a laser-cooled atomic sample of $^{85}$Rb at 1 $μ$K. After the wavepackets have separated by 42 photon recoil momenta, our pulses maintain a fringe contrast of 90% whereas, for adiabatic rapid passage and conventional $π$ pulses, the contrast is less than 10%. Furthermore, we show how these pulses may be adapted to suppress the detrimental off-resonant excitation that limits other broadband pulse schemes.

preprint2020arXiv

Neural network interpretation using descrambler groups

The lack of interpretability and trust is a much-criticised feature of deep neural networks. In fully connected nets, the signalling between inner layers is scrambled because backpropagation training does not require perceptrons to be arranged in any particular order. The result is a black box; this problem is particularly severe in scientific computing and digital signal processing (DSP), where neutral nets perform abstract mathematical transformations that do not reduce to features or concepts. We present here a group-theoretical procedure that attempts to bring inner layer signalling into a human-readable form, the assumption being that this form exists and has identifiable and quantifiable features - for example, smoothness or locality. We applied the proposed method to DEERNet (a DSP network used in electron spin resonance) and managed to descramble it. We found considerable internal sophistication: the network spontaneously invents a bandpass filter, a notch filter, a frequency axis rectifier, frequency division multiplexing, group embedding, spectral filtering regularisation, and a Fourier-like transform that maps harmonic functions into Chebyshev polynomials - in ten minutes of unattended training from a random initial guess.

preprint2016arXiv

Fokker-Planck formalism in magnetic resonance simulations

This paper presents an overview of the Fokker-Planck formalism for non-biological magnetic resonance simulations, describes its existing applications and proposes some novel ones. The most attractive feature of Fokker-Planck theory compared to the commonly used Liouville - von Neumann equation is that, for all relevant types of spatial dynamics (spinning, diffusion, flow, etc.), the corresponding Fokker-Planck Hamiltonian is time-independent. Many difficult NMR, EPR and MRI simulation problems (multiple rotation NMR, ultrafast NMR, gradient-based zero-quantum filters, diffusion and flow NMR, off-resonance soft microwave pulses in EPR, spin-spin coupling effects in MRI, etc.) are simplified significantly in Fokker-Planck space. The paper also summarises the author's experiences with writing and using the corresponding modules of the Spinach library - the methods described below have enabled a large variety of simulations previously considered too complicated for routine practical use.

preprint2016arXiv

Modified Newton-Raphson GRAPE methods for optimal control of spin systems

Quadratic convergence throughout the active space is achieved for the gradient ascent pulse engineering (GRAPE) family of quantum optimal control algorithms. We demonstrate in this communication that the Hessian of the GRAPE fidelity functional is unusually cheap, having the same asymptotic complexity scaling as the functional itself. This leads to the possibility of using very efficient numerical optimization techniques. In particular, the Newton-Raphson method with a rational function optimization (RFO) regularized Hessian is shown in this work to require fewer system trajectory evaluations than any other algorithm in the GRAPE family. This communication describes algebraic and numerical implementation aspects (matrix exponential recycling, Hessian regularization, etc.) for the RFO Newton-Raphson version of GRAPE and reports benchmarks for common spin state control problems in magnetic resonance spectroscopy.

preprint2016arXiv

Pseudocontact shifts from mobile spin labels

This paper presents a detailed analysis of the pseudocontact shift (PCS) field induced by a mobile spin label that is viewed as a probability density distribution with an associated effective magnetic susceptibility anisotropy. It is demonstrated that non-spherically-symmetric density can lead to significant deviations from the commonly used point dipole approximation for PCS. Analytical and numerical solutions are presented for the general partial differential equation that describes the non-point case. It is also demonstrated that it is possible, with some reasonable approximations, to reconstruct paramagnetic centre probability distributions from the experimental PCS data.

preprint2015arXiv

Auxiliary matrix formalism for interaction representation transformations, optimal control and spin relaxation theories

Auxiliary matrix exponential method is used to derive simple and numerically efficient general expressions for the following, historically rather cumbersome and hard to compute, theoretical methods: (1) average Hamiltonian theory following interaction representation transformations; (2) Bloch-Redfield-Wangsness theory of nuclear and electron relaxation; (3) gradient ascent pulse engineering version of quantum optimal control theory. In the context of spin dynamics, the auxiliary matrix exponential method is more efficient than methods based on matrix factorizations and also exhibits more favourable complexity scaling with the dimension of the Hamiltonian matrix.

preprint2015arXiv

Training Schrödinger's cat: quantum optimal control

It is control that turns scientific knowledge into useful technology: in physics and engineering it provides a systematic way for driving a system from a given initial state into a desired target state with minimized expenditure of energy and resources -- as famously applied in the Apollo programme. As one of the cornerstones for enabling quantum technologies, optimal quantum control keeps evolving and expanding into areas as diverse as quantum-enhanced sensing, manipulation of single spins, photons, or atoms, optical spectroscopy, photochemistry, magnetic resonance (spectroscopy as well as medical imaging), quantum information processing and quantum simulation. --- Here state-of-the-art quantum control techniques are reviewed and put into perspective by a consortium uniting expertise in optimal control theory and applications to spectroscopy, imaging, quantum dynamics of closed and open systems. We address key challenges and sketch a roadmap to future developments.

preprint2014arXiv

A partial differential equation for pseudocontact shift

It is demonstrated that pseudocontact shift (PCS), viewed as a scalar or a tensor field in three dimensions, obeys an elliptic partial differential equation with a source term that depends on the Hessian of the unpaired electron probability density. The equation enables straightforward PCS prediction as well as analysis of experimental PCS data in systems with multiple and / or distributed unpaired electron centres.

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

A standard format and a graphical user interface for spin system specification

We introduce a simple and general XML format for spin system description that is the result of extensive consultations within Magnetic Resonance community and unifies under one roof all major existing spin interaction specification conventions. The format is human-readable, easy to edit and easy to parse using standard XML libraries. We also describe a graphical user interface that was designed to facilitate construction and visualization of complicated spin systems. The interface is capable of generating input files for several popular spin dynamics simulation packages.

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.

preprint2012arXiv

Spin system trajectory analysis under optimal control pulses

Several methods are proposed for the analysis, visualization and interpretation of high-dimensional spin system trajectories produced by quantum mechanical simulations. It is noted that expectation values of specific observables in large spin systems often feature fast, complicated and hard-to-interpret time dynamics and suggested that populations of carefully selected subspaces of states are much easier to analyze and interpret. As an illustration of the utility of the proposed methods, it is demonstrated that the apparent "noisy" appearance of many optimal control pulses in NMR and EPR spectroscopy is an illusion - the underlying spin dynamics is shown to be smooth, orderly and very tightly controlled.

preprint2011arXiv

Chemically Induced Dynamic Nuclear Polarization of 19F Nuclei

This study explores, both theoretically and experimentally, the photochemically induced dynamic nuclear polarization (photo-CIDNP) of 19F nuclei, the associated spin relaxation, cross-relaxation and cross-correlation effects, as well as potential applications of 19F CIDNP to protein structure and folding problems. Keywords: CIDNP, relaxation, fluorine, NMR, NOE, TR-CIDNP, tc5b, Trp-cage, GFP

preprint2011arXiv

Molecular structure refinement by direct fitting of atomic coordinates to experimental ESR spectra

An attempt is made to bypass spectral analysis and fit internal coordinates of radicals directly to experimental liquid- and solid-state electron spin resonance (ESR) spectra. We take advantage of the recently introduced large-scale spin dynamics simulation algorithms and of the fact that the accuracy of quantum mechanical calculations of ESR parameters has improved to the point of quantitative correctness. Partial solutions are offered to the local minimum problem in spectral fitting and to the problem of spin interaction parameters (hyperfine couplings, chemical shifts, etc.) being very sensitive to distortions in molecular geometry.

preprint2011arXiv

On the accuracy of the state space restriction approximation for spin dynamics simulations

We present an algebraic foundation for the state space restriction approximation in spin dynamics simulations and derive applicability criteria as well as minimal basis set requirements for practically encountered simulation tasks. The results are illustrated with NMR, ESR, DNP and Spin Chemistry simulations. It is demonstrated that state space restriction yields accurate results in systems where the time scale of spin relaxation processes approximately matches the time scale of the experiment. Rigorous error bounds and basis set requirements are derived.

preprint2011arXiv

Parallel density matrix propagation in spin dynamics simulations

Several methods for density matrix propagation in distributed computing environments, such as clusters and graphics processing units, are proposed and evaluated. It is demonstrated that the large communication overhead associated with each propagation step (two-sided multiplication of the density matrix by an exponential propagator and its conjugate) may be avoided and the simulation recast in a form that requires virtually no inter-thread communication. Good scaling is demonstrated on a 128-core (16 nodes, 8 cores each) cluster.

preprint2011arXiv

Second order gradient ascent pulse engineering

We report some improvements to the gradient ascent pulse engineering (GRAPE) algorithm for optimal control of quantum systems. These include more accurate gradients, convergence acceleration using the BFGS quasi-Newton algorithm as well as faster control derivative calculation algorithms. In all test systems, the wall clock time and the convergence rates show a considerable improvement over the approximate gradient ascent.

preprint2010arXiv

Diagonalization-free implementation of spin relaxation theory for large spin systems

The Liouville space spin relaxation theory equations are reformulated in such a way as to avoid the computationally expensive Hamiltonian diagonalization step, replacing it by numerical evaluation of the integrals in the generalized cumulant expansion. The resulting algorithm is particularly useful in the cases where the static part of the Ha-miltonian is dominated by interactions other than Zeeman (e.g. in quadrupolar reson-ance, low-field EPR and Spin Chemistry). When used together with state space re-striction tools, the algorithm reported is capable of computing full relaxation supero-perators for NMR systems with more than 15 spins.

preprint2009arXiv

Derivatives of spin dynamics simulations

We report analytical equations for the derivatives of spin dynamics simulations with respect to pulse sequence and spin system parameters. The methods described are significantly faster, more accurate and more reliable than the finite difference approximations typically employed. The resulting derivatives may be used in fitting, optimization, performance evaluation and stability analysis of spin dynamics simulations and experiments. Keywords: NMR, EPR, simulation, analytical derivatives, optimal control, spin chemistry, radical pair.

preprint2007arXiv

Polynomially scaling spin dynamics simulation algorithm based on adaptive state space restriction

The conventional spin dynamics simulations are performed in direct products of state spaces of individual spins. In a general system of n spins, the total number of elements in the state basis is >4^n. A system propagation step requires an action by an operator on the state vector and thus requires >4^2n multiplications. It is obvious that with current computers there is no way beyond about ten spins, and the calculation complexity scales exponentially with the spin system size. We demonstrate that a polynomially scaling algorithm can be obtained if the state space is reduced by neglecting unimportant or unpopulated spin states. The class of such states is surprisingly wide. In particular, there are indications that very high multi-spin orders can be dropped completely, as can all the orders linking the spins that are remote on the interaction graph. The computational cost of the propagation step for a ktuples-restricted densely connected n-spin system with k<<n is O(n^2k). In cases of favourable interaction topologies (narrow graphs, e.g. in protein NMR) the asymptotic scaling is linear.