OPTESIM, a Versatile Toolbox for Numerical Simulation of Electron Spin Echo Envelope Modulation (ESEEM) that Features Hybrid Optimization and Statistical Assessment of Parameters

• Type of Material: Article
• Authors: Li Sun, Emory University; Jessica Hernandez-Guzman, Emory University; Kurt Warncke, Emory University
• Language: English
• Date: 2009-09
• Publisher: Elsevier
• Publication Version: Author Accepted Manuscript (After Peer Review)
• Copyright Statement: © 2009 Elsevier Inc. All rights reserved.
• License: Creative Commons Attribution-NonCommerical-NoDerivs 3.0 Unported
• Final Published Version (URL): http://dx.doi.org/10.1016/j.jmr.2009.05.012
• Title of Journal or Parent Work: Journal of Magnetic Resonance
• ISSN: 1090-7807
• Volume: 200
• Issue: 1
• Start Page: 21
• End Page: 28
• Research Categories: Physics, General

Abstract

Electron spin echo envelope modulation (ESEEM) is a technique of pulsed-electron paramagnetic resonance (EPR) spectroscopy. The analyis of ESEEM data to extract information about the nuclear and electronic structure of a disordered (powder) paramagnetic system requires accurate and efficient numerical simulations. A single coupled nucleus of known nuclear g value (gN) and spin I=1 can have up to eight adjustable parameters in the nuclear part of the spin Hamiltonian. We have developed OPTESIM, an ESEEM simulation toolbox, for automated numerical simulation of powder two- and three-pulse one-dimensional ESEEM for arbitrary number (N) and type (I, gN) of coupled nuclei, and arbitrary mutual orientations of the hyperfine tensor principal axis systems for N>1. OPTESIM is based in the Matlab environment, and includes the following features: (1) a fast algorithm for translation of the spin Hamiltonian into simulated ESEEM, (2) different optimization methods that can be hybridized to achieve an efficient coarse-to-fine grained search of the parameter space and convergence to a global minimum, (3) statistical analysis of the simulation parameters, which allows the identification of simultaneous confidence regions at specific confidence levels. OPTESIM also includes a geometry-preserving spherical averaging algorithm as default for N>1, and global optimization over multiple experimental conditions, such as the dephasing time ( ) for three-pulse ESEEM, and external magnetic field values. Application examples for simulation of 14N coupling (N=1, N=2) in biological and chemical model paramagnets are included. Automated, optimized simulations by using OPTESIM lead to a convergence on dramatically shorter time scales, relative to manual simulations.