1
|
De Lorenzi F, Weinmann T, Bruderer S, Heitmann B, Henrici A, Stingelin S. Bayesian analysis of 1D 1H-NMR spectra. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2024; 364:107723. [PMID: 38936240 DOI: 10.1016/j.jmr.2024.107723] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/25/2024] [Revised: 06/14/2024] [Accepted: 06/14/2024] [Indexed: 06/29/2024]
Abstract
Extracting spin system parameters from 1D high resolution 1H-NMR spectra can be an intricate task requiring sophisticate methods. With a few exceptions methods to perform such a total line shape analysis commonly rely on local optimization techniques which for increasing complexity of the underlying spin system tend to reveal local solutions. In this work we propose a full Bayesian modeling approach based on a quantum mechanical model of the spin system. The Bayesian formalism provides a global optimization strategy which allows to efficiently include prior knowledge about the spin system or to incorporate additional constraints concerning the parameters of interest. The proposed algorithm has been tested on synthetic and real 1D 1H-NMR data for various spin systems with increasing complexity. The results show that the Bayesian algorithm provides accurate estimates even for complex spectra with many overlapping regions, and that it can cope with symmetry induced local minima. By providing an unbiased estimate of the model evidence the proposed algorithm furthermore offers a way to discriminate between different spin system candidates.
Collapse
Affiliation(s)
- Flavio De Lorenzi
- Institute of Applied Mathematics and Physics, Zurich University of Applied Sciences, Technikumstr. 71, 8400 Winterthur, Switzerland.
| | - Tom Weinmann
- Institute of Applied Mathematics and Physics, Zurich University of Applied Sciences, Technikumstr. 71, 8400 Winterthur, Switzerland
| | - Simon Bruderer
- Bruker Switzerland AG, Industriestr. 26, 8117 Fällanden, Switzerland
| | - Björn Heitmann
- Bruker Switzerland AG, Industriestr. 26, 8117 Fällanden, Switzerland
| | - Andreas Henrici
- Institute of Applied Mathematics and Physics, Zurich University of Applied Sciences, Technikumstr. 71, 8400 Winterthur, Switzerland
| | - Simon Stingelin
- Institute of Applied Mathematics and Physics, Zurich University of Applied Sciences, Technikumstr. 71, 8400 Winterthur, Switzerland
| |
Collapse
|
2
|
Vinding MS, Goodwin DL, Kuprov I, Lund TE. Optimal control gradient precision trade-offs: Application to fast generation of DeepControl libraries for MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2021; 333:107094. [PMID: 34794089 DOI: 10.1016/j.jmr.2021.107094] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2021] [Revised: 09/21/2021] [Accepted: 10/19/2021] [Indexed: 06/13/2023]
Abstract
We have recently demonstrated supervised deep learning methods for rapid generation of radiofrequency pulses in magnetic resonance imaging (https://doi.org/10.1002/mrm.27740, https://doi.org/10.1002/mrm.28667). Unlike the previous iterative optimization approaches, deep learning methods generate a pulse using a fixed number of floating-point operations - this is important in MRI, where patient-specific pulses preferably must be produced in real time. However, deep learning requires vast training libraries, which must be generated using the traditional methods, e.g., iterative quantum optimal control methods. Those methods are usually variations of gradient descent, and the calculation of the gradient of the performance metric with respect to the pulse waveform can be the most numerically intensive step. In this communication, we explore various ways in which the calculation of gradients in quantum optimal control theory may be accelerated. Four optimization avenues are explored: truncated commutator series expansions at zeroth and first order, a novel midpoint truncation scheme at first order, and the exact complex-step method. For the spin systems relevant to MRI, the first-order midpoint truncation is found to be sufficiently accurate, but also significantly faster than the machine precision gradient. This makes the generation of training databases for the machine learning methods considerably more realistic.
Collapse
Affiliation(s)
- Mads Sloth Vinding
- Center of Functionally Integrative Neuroscience (CFIN), Department of Clinical Medicine, Faculty of Health, Aarhus University, Denmark.
| | - David L Goodwin
- Institute for Biological Interfaces 4 - Magnetic Resonance, Karlsruhe Institute for Technology (KIT), Karlsruhe, Germany; Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford, UK.
| | - Ilya Kuprov
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, UK
| | - Torben Ellegaard Lund
- Center of Functionally Integrative Neuroscience (CFIN), Department of Clinical Medicine, Faculty of Health, Aarhus University, Denmark
| |
Collapse
|
3
|
Machnes S, Assémat E, Tannor D, Wilhelm FK. Tunable, Flexible, and Efficient Optimization of Control Pulses for Practical Qubits. PHYSICAL REVIEW LETTERS 2018; 120:150401. [PMID: 29756895 DOI: 10.1103/physrevlett.120.150401] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Revised: 12/05/2017] [Indexed: 06/08/2023]
Abstract
Quantum computation places very stringent demands on gate fidelities, and experimental implementations require both the controls and the resultant dynamics to conform to hardware-specific constraints. Superconducting qubits present the additional requirement that pulses must have simple parameterizations, so they can be further calibrated in the experiment, to compensate for uncertainties in system parameters. Other quantum technologies, such as sensing, require extremely high fidelities. We present a novel, conceptually simple and easy-to-implement gradient-based optimal control technique named gradient optimization of analytic controls (GOAT), which satisfies all the above requirements, unlike previous approaches. To demonstrate GOAT's capabilities, with emphasis on flexibility and ease of subsequent calibration, we optimize fast coherence-limited pulses for two leading superconducting qubits architectures-flux-tunable transmons and fixed-frequency transmons with tunable couplers.
Collapse
Affiliation(s)
- Shai Machnes
- Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
- Weizmann Institute of Science, 76100 Rehovot, Israel
| | - Elie Assémat
- Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
| | - David Tannor
- Weizmann Institute of Science, 76100 Rehovot, Israel
| | - Frank K Wilhelm
- Theoretical Physics, Saarland University, 66123 Saarbrücken, Germany
| |
Collapse
|
4
|
Majewski K. Rotation relaxation splitting for optimizing parallel RF excitation pulses with T 1- and T 2-relaxations in MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2018; 288:43-57. [PMID: 29414063 DOI: 10.1016/j.jmr.2018.01.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2017] [Revised: 12/18/2017] [Accepted: 01/09/2018] [Indexed: 06/08/2023]
Abstract
Exact solutions of the Bloch equations with T1- and T2-relaxation terms for piecewise constant magnetic fields are numerically challenging. We therefore investigate an approximation for the achieved magnetization in which rotations and relaxations are split into separate operations. We develop an estimate for its accuracy and explicit first and second order derivatives with respect to the complex excitation radio frequency voltages. In practice, the deviation between an exact solution of the Bloch equations and this rotation relaxation splitting approximation seems negligible. Its computation times are similar to exact solutions without relaxation terms. We apply the developed theory to numerically optimize radio frequency excitation waveforms with T1- and T2-relaxations in several examples.
Collapse
Affiliation(s)
- Kurt Majewski
- Siemens AG, CT RDA BAM ORD-DE, 80200 Munich, Germany.
| |
Collapse
|
5
|
Majewski K, Ritter D. First and second order derivatives for optimizing parallel RF excitation waveforms. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2015; 258:65-80. [PMID: 26232364 DOI: 10.1016/j.jmr.2015.06.010] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/20/2015] [Revised: 06/12/2015] [Accepted: 06/14/2015] [Indexed: 06/04/2023]
Abstract
For piecewise constant magnetic fields, the Bloch equations (without relaxation terms) can be solved explicitly. This way the magnetization created by an excitation pulse can be written as a concatenation of rotations applied to the initial magnetization. For fixed gradient trajectories, the problem of finding parallel RF waveforms, which minimize the difference between achieved and desired magnetization on a number of voxels, can thus be represented as a finite-dimensional minimization problem. We use quaternion calculus to formulate this optimization problem in the magnitude least squares variant and specify first and second order derivatives of the objective function. We obtain a small tip angle approximation as first order Taylor development from the first order derivatives and also develop algorithms for first and second order derivatives for this small tip angle approximation. All algorithms are accompanied by precise floating point operation counts to assess and compare the computational efforts. We have implemented these algorithms as callback functions of an interior-point solver. We have applied this numerical optimization method to example problems from the literature and report key observations.
Collapse
Affiliation(s)
- Kurt Majewski
- Siemens AG, CT RTC BAM ORD-DE, 80200 Munich, Germany.
| | - Dieter Ritter
- Siemens AG, HC IM MR R&D SYS PHYS, Post Box 32 60, 91050 Erlangen, Germany.
| |
Collapse
|
6
|
Maximov II, Vinding MS, Tse DHY, Nielsen NC, Shah NJ. Real-time 2D spatially selective MRI experiments: Comparative analysis of optimal control design methods. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2015; 254:110-120. [PMID: 25863895 DOI: 10.1016/j.jmr.2015.03.003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2014] [Revised: 03/12/2015] [Accepted: 03/13/2015] [Indexed: 06/04/2023]
Abstract
There is an increasing need for development of advanced radio-frequency (RF) pulse techniques in modern magnetic resonance imaging (MRI) systems driven by recent advancements in ultra-high magnetic field systems, new parallel transmit/receive coil designs, and accessible powerful computational facilities. 2D spatially selective RF pulses are an example of advanced pulses that have many applications of clinical relevance, e.g., reduced field of view imaging, and MR spectroscopy. The 2D spatially selective RF pulses are mostly generated and optimised with numerical methods that can handle vast controls and multiple constraints. With this study we aim at demonstrating that numerical, optimal control (OC) algorithms are efficient for the design of 2D spatially selective MRI experiments, when robustness towards e.g. field inhomogeneity is in focus. We have chosen three popular OC algorithms; two which are gradient-based, concurrent methods using first- and second-order derivatives, respectively; and a third that belongs to the sequential, monotonically convergent family. We used two experimental models: a water phantom, and an in vivo human head. Taking into consideration the challenging experimental setup, our analysis suggests the use of the sequential, monotonic approach and the second-order gradient-based approach as computational speed, experimental robustness, and image quality is key. All algorithms used in this work were implemented in the MATLAB environment and are freely available to the MRI community.
Collapse
Affiliation(s)
- Ivan I Maximov
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.
| | - Mads S Vinding
- Center for Insoluble Protein Structures (inSPIN), Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark.
| | - Desmond H Y Tse
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Niels Chr Nielsen
- Center for Insoluble Protein Structures (inSPIN), Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| | - N Jon Shah
- Institute of Neuroscience and Medicine 4, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; Department of Neurology, Faculty of Medicine, RWTH Aachen University, JARA, 52074 Aachen, Germany
| |
Collapse
|
7
|
Ehni S, Luy B. BEBE(tr) and BUBI: J-compensated concurrent shaped pulses for 1H-13C experiments. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2013; 232:7-17. [PMID: 23673080 DOI: 10.1016/j.jmr.2013.04.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2013] [Revised: 03/22/2013] [Accepted: 04/12/2013] [Indexed: 06/02/2023]
Abstract
Shaped pulses designed for broadband excitation, inversion and refocusing are important tools in modern NMR spectroscopy to achieve robust pulse sequences especially in heteronuclear correlation experiments. A large variety of mostly computer-optimized pulse shapes exist for different desired bandwidths, available rf-field strengths, and tolerance to B1-inhomogeneity. They are usually derived for a single spin 1/2, neglecting evolution due to J-couplings. While pulses with constant resulting phase are selfcompensated for heteronuclear coupling evolution as long as they are applied exclusively on a single nucleus, the situation changes for concurrently applied pulse shapes. Using the example of a (1)H,(13)C two spin system, two J-compensated pulse pairs for the application in INEPT-type transfer elements were optimized: a point-to-point pulse sandwich called BEBE(tr), consisting of a broadband excitation and time-reversed excitation pulse, and a combined universal rotation and point-to-point pulse pair called BUBI, which acts as a refocusing pulse on (1)H and a corresponding inversion pulse on (13)C. After a derivation of quality factors and optimization protocols, a theoretical and experimental comparison with conventionally derived BEBOP, BIBOP, and BURBOP-180° pulses is given. While the overall transfer efficiency of a single pulse pair is only reduced by approximately 0.1%, resulting transfer to undesired coherences is reduced by several percent. In experiments this can lead to undesired phase distortions for pairs of uncompensated pulse shapes and even differences in signal intensities of 5-10% in HSQC and up to 68% in more complex COB-HSQC experiments.
Collapse
Affiliation(s)
- Sebastian Ehni
- Institute of Organic Chemistry and Institute for Biological Interfaces, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 6, 76131 Karlsruhe, Germany
| | | |
Collapse
|
8
|
Schilling F, Glaser SJ. Tailored real-time scaling of heteronuclear couplings. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2012; 223:207-218. [PMID: 22982752 DOI: 10.1016/j.jmr.2012.07.024] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2012] [Revised: 07/21/2012] [Accepted: 07/31/2012] [Indexed: 06/01/2023]
Abstract
Heteronuclear couplings are a valuable source of molecular information, which is measured from the multiplet splittings of an NMR spectrum. Radiofrequency irradiation on one coupled nuclear spin allows to modify the effective coupling constant, scaling down the multiplet splittings in the spectrum observed at the resonance frequency of the other nuclear spin. Such decoupling sequences are often used to collapse a multiplet into a singlet and can therefore simplify NMR spectra significantly. Continuous-wave (cw) decoupling has an intrinsic non-linear offset dependence of the scaling of the effective J-coupling constant. Using optimal control pulse optimization, we show that virtually arbitrary off-resonance scaling of the J-coupling constant can be achieved. The new class of tailored decoupling pulses is named SHOT (Scaling of Heteronuclear couplings by Optimal Tracking). Complementing cw irradiation, SHOT pulses offer an alternative approach of encoding chemical shift information indirectly through off-resonance decoupling, which however makes it possible for the first time to achieve linear J scaling as a function of offset frequency. For a simple mixture of eight aromatic compounds, it is demonstrated experimentally that a 1D-SHOT {(1)H}-(13)C experiment yields comparable information to a 2D-HSQC and can give full assignment of all coupled spins.
Collapse
Affiliation(s)
- Franz Schilling
- Department of Chemistry, Technische Universität München, 85747 Garching, Germany
| | | |
Collapse
|
9
|
Momot KI, Takegoshi K. Sensitivity of the NMR density matrix to pulse sequence parameters: a simplified analytic approach. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2012; 221:57-68. [PMID: 22750252 DOI: 10.1016/j.jmr.2012.05.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2012] [Revised: 04/25/2012] [Accepted: 05/05/2012] [Indexed: 06/01/2023]
Abstract
We present a formalism for the analysis of sensitivity of nuclear magnetic resonance pulse sequences to variations of pulse sequence parameters, such as radiofrequency pulses, gradient pulses or evolution delays. The formalism enables the calculation of compact, analytic expressions for the derivatives of the density matrix and the observed signal with respect to the parameters varied. The analysis is based on two constructs computed in the course of modified density-matrix simulations: the error interrogation operators and error commutators. The approach presented is consequently named the Error Commutator Formalism (ECF). It is used to evaluate the sensitivity of the density matrix to parameter variation based on the simulations carried out for the ideal parameters, obviating the need for finite-difference calculations of signal errors. The ECF analysis therefore carries a computational cost comparable to a single density-matrix or product-operator simulation. Its application is illustrated using a number of examples from basic NMR spectroscopy. We show that the strength of the ECF is its ability to provide analytic insights into the propagation of errors through pulse sequences and the behaviour of signal errors under phase cycling. Furthermore, the approach is algorithmic and easily amenable to implementation in the form of a programming code. It is envisaged that it could be incorporated into standard NMR product-operator simulation packages.
Collapse
Affiliation(s)
- Konstantin I Momot
- School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, GPO Box 2434, Brisbane, Qld 4001, Australia.
| | | |
Collapse
|
10
|
Anand CK, Bain AD, Curtis AT, Nie Z. Designing optimal universal pulses using second-order, large-scale, non-linear optimization. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2012; 219:61-74. [PMID: 22617160 DOI: 10.1016/j.jmr.2012.04.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2011] [Revised: 03/26/2012] [Accepted: 04/03/2012] [Indexed: 06/01/2023]
Abstract
Recently, RF pulse design using first-order and quasi-second-order pulses has been actively investigated. We present a full second-order design method capable of incorporating relaxation, inhomogeneity in B(0) and B(1). Our model is formulated as a generic optimization problem making it easy to incorporate diverse pulse sequence features. To tame the computational cost, we present a method of calculating second derivatives in at most a constant multiple of the first derivative calculation time, this is further accelerated by using symbolic solutions of the Bloch equations. We illustrate the relative merits and performance of quasi-Newton and full second-order optimization with a series of examples, showing that even a pulse already optimized using other methods can be visibly improved. To be useful in CPMG experiments, a universal refocusing pulse should be independent of the delay time and insensitive of the relaxation time and RF inhomogeneity. We design such a pulse and show that, using it, we can obtain reliable R(2) measurements for offsets within ±γB(1). Finally, we compare our optimal refocusing pulse with other published refocusing pulses by doing CPMG experiments.
Collapse
|
11
|
de Fouquieres P. Implementing quantum gates by optimal control with doubly exponential convergence. PHYSICAL REVIEW LETTERS 2012; 108:110504. [PMID: 22540447 DOI: 10.1103/physrevlett.108.110504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2011] [Indexed: 05/31/2023]
Abstract
We introduce a novel algorithm for the task of coherently controlling a quantum mechanical system to implement any chosen unitary dynamics. It performs faster than existing state of the art methods by 1 to 3 orders of magnitude (depending on which one we compare to), particularly for quantum information processing purposes. This substantially enhances the ability to both study the control capabilities of physical systems within their coherence times, and constrain solutions for control tasks to lie within experimentally feasible regions. Natural extensions of the algorithm are also discussed.
Collapse
Affiliation(s)
- Pierre de Fouquieres
- Centre for Quantum Information and Foundations, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
12
|
Charnock GTP, Krzystyniak M, Kuprov I. Molecular structure refinement by direct fitting of atomic coordinates to experimental ESR spectra. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2012; 216:62-68. [PMID: 22300803 DOI: 10.1016/j.jmr.2012.01.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2011] [Revised: 01/04/2012] [Accepted: 01/06/2012] [Indexed: 05/31/2023]
Abstract
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 vibrational excursions from the equilibrium geometry.
Collapse
Affiliation(s)
- G T P Charnock
- Oxford e-Research Centre, University of Oxford, 7 Keble Road, Oxford OX1 3QG, UK
| | | | | |
Collapse
|
13
|
de Fouquieres P, Schirmer SG, Glaser SJ, Kuprov I. Second order gradient ascent pulse engineering. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2011; 212:412-417. [PMID: 21885306 DOI: 10.1016/j.jmr.2011.07.023] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2011] [Revised: 07/22/2011] [Accepted: 07/28/2011] [Indexed: 05/31/2023]
Abstract
We report some improvements to the gradient ascent pulse engineering (GRAPE) algorithm for optimal control of spin ensembles and other quantum systems. These include more accurate gradients, convergence acceleration using the Broyden-Fletcher-Goldfarb-Shanno (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.
Collapse
Affiliation(s)
- P de Fouquieres
- Centre for Quantum Computation, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
| | | | | | | |
Collapse
|
14
|
Castillo AM, Patiny L, Wist J. Fast and accurate algorithm for the simulation of NMR spectra of large spin systems. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2011; 209:123-130. [PMID: 21316274 DOI: 10.1016/j.jmr.2010.12.008] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2010] [Revised: 12/21/2010] [Accepted: 12/23/2010] [Indexed: 05/30/2023]
Abstract
The computational cost for the simulation of NMR spectra grows exponentially with the number of nuclei. Today, the memory available to store the Hamiltonian limits the size of the system that can be studied. Modern computers enable to tackle systems containing up to 13 spins [1], which obviously does not allow to study most molecules of interest in research. This issue can be addressed by identifying groups of spins or fragments that are not or only weakly interacting together, i.e., that only share weakly coupled spin pairs. Such a fragmentation is only permitted in the weak coupling regime, i.e., when the coupling interaction is weak compared to the difference in chemical shift of the coupled spins. Here, we propose a procedure that removes weak coupling interactions in order to split the spin system efficiently and to correct a posteriori for the effect of the neglected couplings. This approach yields accurate spectra when the adequate interactions are removed, i.e., between spins only involved in weak coupling interactions, but fails otherwise. As a result, the computational time for the simulation of 1D spectra grows linearly with the size of the spin system.
Collapse
Affiliation(s)
- Andrés M Castillo
- Departamento de Química, Universidad Nacional de Colombia, Bogotá D.C., Colombia
| | | | | |
Collapse
|
15
|
Hogben HJ, Krzystyniak M, Charnock GTP, Hore PJ, Kuprov I. Spinach--a software library for simulation of spin dynamics in large spin systems. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2011; 208:179-94. [PMID: 21169043 DOI: 10.1016/j.jmr.2010.11.008] [Citation(s) in RCA: 270] [Impact Index Per Article: 20.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2010] [Revised: 11/07/2010] [Accepted: 11/11/2010] [Indexed: 05/12/2023]
Abstract
We introduce a software library incorporating our recent research into efficient simulation algorithms for large spin systems. Liouville space simulations (including symmetry, relaxation and chemical kinetics) of most liquid-state NMR experiments on 40+ spin systems can now be performed without effort on a desktop workstation. Much progress has also been made with improving the efficiency of ESR, solid state NMR and Spin Chemistry simulations. Spinach is available for download at http://spindynamics.org.
Collapse
Affiliation(s)
- H J Hogben
- Physical and Theoretical Chemistry Laboratory, Chemistry Department, University of Oxford, South Parks Road, Oxford OX1 3QZ, UK
| | | | | | | | | |
Collapse
|