Christiana Sabett, Ariel Hafftka, Kyle Sexton, Richard G. Spencer
Determination of the distribution of magnetic resonance (MR) transverse relaxation times is emerging as an important method for materials characterization, including assessment of tissue pathology in biomedicine. These distributions are obtained from the inverse Laplace transform (ILT) of multiexponential decay data. Stabilization of this classically ill-posed problem is most commonly attempted using Tikhonov regularization with an L2 penalty term. However, with the availability of convex optimization algorithms and recognition of the importance of sparsity in model reconstruction, there has been increasing interest in alternative penalties. The L1 penalty enforces a greater degree of sparsity than L2, and so may be suitable for highly localized relaxation time distributions. In addition, Lp penalties, 1 < p < 2, and the elastic net (EN) penalty, defined as a linear combination of L1 and L2 penalties, may be appropriate for distributions consisting of both narrow and broad components. We evaluate the L1, L2, Lp, and EN penalties for model relaxation time distributions consisting of two Gaussian peaks. For distributions with narrow Gaussian peaks, the L1 penalty works well to maintain sparsity and promote resolution, while the conventional L2 penalty performs best for distributions with broader peaks. Finally, the Lp and EN penalties do in fact outperform the L1 and L2 penalties for distributions with components of unequal widths. These findings serve as indicators of appropriate regularization in the typical situation in which the experimentalist has a priori knowledge of the general characteristics of the underlying relaxation time distribution. Our findings can be applied to both the recovery of T2 distributions from spin echo decay data as well as distributions of other MR parameters, such as apparent diffusion constant, from their multiexponential decay signals.
{"title":"L1, Lp, L2, and elastic net penalties for regularization of Gaussian component distributions in magnetic resonance relaxometry","authors":"Christiana Sabett, Ariel Hafftka, Kyle Sexton, Richard G. Spencer","doi":"10.1002/cmr.a.21427","DOIUrl":"10.1002/cmr.a.21427","url":null,"abstract":"<p>Determination of the distribution of magnetic resonance (MR) transverse relaxation times is emerging as an important method for materials characterization, including assessment of tissue pathology in biomedicine. These distributions are obtained from the inverse Laplace transform (ILT) of multiexponential decay data. Stabilization of this classically ill-posed problem is most commonly attempted using Tikhonov regularization with an L<sub>2</sub> penalty term. However, with the availability of convex optimization algorithms and recognition of the importance of sparsity in model reconstruction, there has been increasing interest in alternative penalties. The L<sub>1</sub> penalty enforces a greater degree of sparsity than L<sub>2</sub>, and so may be suitable for highly localized relaxation time distributions. In addition, L<sub><i>p</i></sub> penalties, 1 < <i>p </i>< 2, and the elastic net (EN) penalty, defined as a linear combination of L<sub>1</sub> and L<sub>2</sub> penalties, may be appropriate for distributions consisting of both narrow and broad components. We evaluate the L<sub>1</sub>, L<sub>2</sub>, L<sub><i>p</i></sub>, and EN penalties for model relaxation time distributions consisting of two Gaussian peaks. For distributions with narrow Gaussian peaks, the L<sub>1</sub> penalty works well to maintain sparsity and promote resolution, while the conventional L<sub>2</sub> penalty performs best for distributions with broader peaks. Finally, the L<sub><i>p</i></sub> and EN penalties do in fact outperform the L<sub>1</sub> and L<sub>2</sub> penalties for distributions with components of unequal widths. These findings serve as indicators of appropriate regularization in the typical situation in which the experimentalist has a priori knowledge of the general characteristics of the underlying relaxation time distribution. Our findings can be applied to both the recovery of T<sub>2</sub> distributions from spin echo decay data as well as distributions of other MR parameters, such as apparent diffusion constant, from their multiexponential decay signals.</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 2","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21427","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90054162","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nonuniform sampling (NUS) offers NMR spectroscopists a means of accelerating data collection and increasing spectral quality in multidimensional (nD) experiments. The data from NUS experiments are incomplete by design, and must be reconstructed prior to use. While most existing reconstruction techniques compute point estimates of the true signal, Bayesian statistics offers a means of estimating posterior distributions over the signal, which enable more rigorous quantitation and uncertainty estimation. In this article, we describe the variational approach to approximating Bayesian posterior distributions, and illustrate how it can be applied to extend existing results from Bayesian spectrum analysis and compressed sensing. The new NUS reconstruction algorithms resulting from variational Bayes are computationally efficient, and offer new insights into the concepts of spectral sparsity and optimal sampling in NMR experiments.
{"title":"Variational Bayesian analysis of nonuniformly sampled NMR data","authors":"Bradley Worley","doi":"10.1002/cmr.a.21428","DOIUrl":"10.1002/cmr.a.21428","url":null,"abstract":"<p>Nonuniform sampling (NUS) offers NMR spectroscopists a means of accelerating data collection and increasing spectral quality in multidimensional (<i>n</i>D) experiments. The data from NUS experiments are incomplete by design, and must be reconstructed prior to use. While most existing reconstruction techniques compute point estimates of the true signal, Bayesian statistics offers a means of estimating <i>posterior distributions</i> over the signal, which enable more rigorous quantitation and uncertainty estimation. In this article, we describe the <i>variational</i> approach to approximating Bayesian posterior distributions, and illustrate how it can be applied to extend existing results from Bayesian spectrum analysis and compressed sensing. The new NUS reconstruction algorithms resulting from variational Bayes are computationally efficient, and offer new insights into the concepts of spectral sparsity and optimal sampling in NMR experiments.</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 2","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21428","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80497629","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Advances in alternative sampling and processing","authors":"David Rovnyak, Adam D. Schuyler","doi":"10.1002/cmr.a.21458","DOIUrl":"10.1002/cmr.a.21458","url":null,"abstract":"","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 2","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21458","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85443133","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The starting point of all NMR experiments is a spin polarization which develops when we place the sample in static magnetic field B0. There are excess of spins aligned along B0 (spin up with lower energy) than spins aligned opposite (spin down with higher energy) to the field B0. A natural question is what is the source of this excess spin polarization because relaxation mechanisms can flip a up spin to a down spin and vice-versa. The answer lies in the density of states. When a molecule with spin down flips to spin up it loses energy. This energy goes into increasing the kinetic energy of the molecule in the gas/solution phase. At this increased kinetic energy, there are more rotational-translational states accessible to the molecule than at lower energy. This increases the probability the molecule will spend in spin up state (higher kinetic energy state). This is the source of excess polarization. In this article, we use an argument based on equipartition of energy to explicitly count the excess states that become accessible to the molecule when its spin is flipped from down to up. Using this counting, we derive the familiar Boltzmann distribution of the ratio of up vs down spins. Although prima facie, there is nothing new in this article, we find the mode counting argument for excess states interesting. Furthermore, the article stresses the fact that spin polarization arises from higher density of states at increased kinetic energy of molecules.
{"title":"Conservation of energy, density of states, and spin lattice relaxation","authors":"Navin Khaneja","doi":"10.1002/cmr.a.21457","DOIUrl":"https://doi.org/10.1002/cmr.a.21457","url":null,"abstract":"<p>The starting point of all NMR experiments is a spin polarization which develops when we place the sample in static magnetic field <i>B</i><sub>0</sub>. There are excess of spins aligned along <i>B</i><sub>0</sub> (spin up with lower energy) than spins aligned opposite (spin down with higher energy) to the field <i>B</i><sub>0</sub>. A natural question is what is the source of this excess spin polarization because relaxation mechanisms can flip a up spin to a down spin and vice-versa. The answer lies in the density of states. When a molecule with spin down flips to spin up it loses energy. This energy goes into increasing the kinetic energy of the molecule in the gas/solution phase. At this increased kinetic energy, there are more rotational-translational states accessible to the molecule than at lower energy. This increases the probability the molecule will spend in spin up state (higher kinetic energy state). This is the source of excess polarization. In this article, we use an argument based on equipartition of energy to explicitly count the excess states that become accessible to the molecule when its spin is flipped from down to up. Using this counting, we derive the familiar Boltzmann distribution of the ratio of up vs down spins. Although prima facie, there is nothing new in this article, we find the mode counting argument for excess states interesting. Furthermore, the article stresses the fact that spin polarization arises from higher density of states at increased kinetic energy of molecules.</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 3","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21457","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91838675","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"NMR Concepts","authors":"","doi":"10.1002/cmr.a.21436","DOIUrl":"https://doi.org/10.1002/cmr.a.21436","url":null,"abstract":"","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 1","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21436","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"109174259","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"DVD Review","authors":"","doi":"10.1002/cmr.a.21368","DOIUrl":"https://doi.org/10.1002/cmr.a.21368","url":null,"abstract":"","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 1","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21368","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"109174258","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A simple approach how to calculate analytical expressions for unbalanced steady-state free precession (ubSSFP) signals of arbitrary dephasing order is presented. Dephasing order is the number of effective gradient dephasing cycles that magnetization immediately after an RF-pulse has experienced during the ubSSFP sequence. Based on the obtained equations, which are in accordance with existing literature, the sensitivity of ubSSFP signals to T2∗ is derived under the assumption of a Lorentzian frequency distribution resulting from static field inhomogeneities. Further, the phases of all ubSSFP signals are calculated and a general expression of how to use them for B0-fieldmapping is given. The derivation is supported by the extended phase graph (EPG) model, and as such the work also acts as a comprehensive introduction to the formal description of SSFP. In addition, the balanced SSFP (bSSFP) sequence is explored. The connection of bSSFP to ubSSFP is shown, and potential T2∗-sensitivity of bSSFP is examined based on numerical simulations. It is shown that ubSSFP signals with negative dephasing order have a refocusing character and behave similar to spin-echo signals. Conversely, ubSSFP signals with zero or positive dephasing order can be regarded as T2∗-weighted. The behavior of bSSFP depends largely on the exact distribution of frequencies. For instance, for a narrow spherical distribution, bSSFP acts like a spin-echo sequence, while for a Lorentzian distribution a refocusing behavior does not occur.
{"title":"Steady-state free precession signals of arbitrary dephasing order and their sensitivity to T2∗","authors":"Jochen Leupold","doi":"10.1002/cmr.a.21435","DOIUrl":"10.1002/cmr.a.21435","url":null,"abstract":"<p>A simple approach how to calculate analytical expressions for unbalanced steady-state free precession (ubSSFP) signals of arbitrary dephasing order is presented. Dephasing order is the number of effective gradient dephasing cycles that magnetization immediately after an RF-pulse has experienced during the ubSSFP sequence. Based on the obtained equations, which are in accordance with existing literature, the sensitivity of ubSSFP signals to <i>T</i><sub>2</sub><sup>∗</sup> is derived under the assumption of a Lorentzian frequency distribution resulting from static field inhomogeneities. Further, the phases of all ubSSFP signals are calculated and a general expression of how to use them for B<sub>0</sub>-fieldmapping is given. The derivation is supported by the extended phase graph (EPG) model, and as such the work also acts as a comprehensive introduction to the formal description of SSFP. In addition, the balanced SSFP (bSSFP) sequence is explored. The connection of bSSFP to ubSSFP is shown, and potential <i>T</i><sub>2</sub><sup>∗</sup>-sensitivity of bSSFP is examined based on numerical simulations. It is shown that ubSSFP signals with negative dephasing order have a refocusing character and behave similar to spin-echo signals. Conversely, ubSSFP signals with zero or positive dephasing order can be regarded as <i>T</i><sub>2</sub><sup>∗</sup>-weighted. The behavior of bSSFP depends largely on the exact distribution of frequencies. For instance, for a narrow spherical distribution, bSSFP acts like a spin-echo sequence, while for a Lorentzian distribution a refocusing behavior does not occur.</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"46A 1","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21435","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87722707","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sandra S. Eaton, Lukas B. Woodcock, Gareth R. Eaton
Nitroxide biradicals have been prepared with electron-electron spin-spin exchange interaction, J, ranging from weak to very strong. EPR spectra of these biradicals in fluid solution depend on the ratio of J to the nitrogen hyperfine coupling, AN, and the rates of interconversion between conformations with different values of J. For relatively rigid biradicals EPR spectra can be simulated as the superposition of AB splitting patterns arising from different combinations of nitrogen nuclear spin states. For more flexible biradicals spectra can be simulated with a Liouville representation of the dynamics that interconvert conformations with different values of J on the EPR timescale. Analysis of spectra, factors that impact J, and examples of applications to chemical and biophysical problems are discussed.
{"title":"Continuous wave electron paramagnetic resonance of nitroxide biradicals in fluid solution","authors":"Sandra S. Eaton, Lukas B. Woodcock, Gareth R. Eaton","doi":"10.1002/cmr.a.21426","DOIUrl":"10.1002/cmr.a.21426","url":null,"abstract":"<p>Nitroxide biradicals have been prepared with electron-electron spin-spin exchange interaction, <i>J</i>, ranging from weak to very strong. EPR spectra of these biradicals in fluid solution depend on the ratio of <i>J</i> to the nitrogen hyperfine coupling, <i>A</i><sub>N</sub>, and the rates of interconversion between conformations with different values of <i>J</i>. For relatively rigid biradicals EPR spectra can be simulated as the superposition of <i>AB</i> splitting patterns arising from different combinations of nitrogen nuclear spin states. For more flexible biradicals spectra can be simulated with a Liouville representation of the dynamics that interconvert conformations with different values of <i>J</i> on the EPR timescale. Analysis of spectra, factors that impact <i>J</i>, and examples of applications to chemical and biophysical problems are discussed.</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"47A 2","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-05-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21426","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45562301","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"DVD Review","authors":"","doi":"10.1002/cmr.a.21431","DOIUrl":"https://doi.org/10.1002/cmr.a.21431","url":null,"abstract":"","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"45A 6","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21431","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"109175464","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Christopher Anand, Bob Berno, Stephen Boulton, Michael A. Brook, Richard Epand, Timothy R. Field, Gillian R. Goward, Paul Hazendonk, Giuseppe Melacini
In this tribute to our friend, mentor and colleague Alexander Davidson Bain we have collectively recapitulated the milestones of his career at McMaster university. We start from Alex's scientific and educational achievements and continue with his accomplishments as a community and infrastructure builder. We attempt to provide a sense of the breadth and depth of his seemingly endless scientific contributions while at McMaster by briefly summarizing selected representative examples from his body of work. Following Alex's lead, the scientific account is mixed with anecdotes and “bits of wisdom” we fondly remember from our interactions and collaborations with him. We also touch upon his brilliant and nurturing educational style and his “aggregator” role within the McMaster and wider NMR communities. We conclude with a more personal picture of Alex D. Bain, in which his scientific excellence and profound intellect are inextricably tied to his kind, nurturing and optimistic character and to his uniquely wry humor. He was not just a “good guy,” he was the epitome of the “good guy.”
{"title":"A tribute to Alexander Davidson Bain: An NMR pioneer and mentor at McMaster University","authors":"Christopher Anand, Bob Berno, Stephen Boulton, Michael A. Brook, Richard Epand, Timothy R. Field, Gillian R. Goward, Paul Hazendonk, Giuseppe Melacini","doi":"10.1002/cmr.a.21418","DOIUrl":"10.1002/cmr.a.21418","url":null,"abstract":"<p>In this tribute to our friend, mentor and colleague Alexander Davidson Bain we have collectively recapitulated the milestones of his career at McMaster university. We start from Alex's scientific and educational achievements and continue with his accomplishments as a community and infrastructure builder. We attempt to provide a sense of the breadth and depth of his seemingly endless scientific contributions while at McMaster by briefly summarizing selected representative examples from his body of work. Following Alex's lead, the scientific account is mixed with anecdotes and “bits of wisdom” we fondly remember from our interactions and collaborations with him. We also touch upon his brilliant and nurturing educational style and his “aggregator” role within the McMaster and wider NMR communities. We conclude with a more personal picture of Alex D. Bain, in which his scientific excellence and profound intellect are inextricably tied to his kind, nurturing and optimistic character and to his uniquely wry humor. He was not just a “good guy,” he was the epitome of the “good guy.”</p>","PeriodicalId":55216,"journal":{"name":"Concepts in Magnetic Resonance Part A","volume":"45A 6","pages":""},"PeriodicalIF":0.6,"publicationDate":"2018-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1002/cmr.a.21418","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81401056","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}