{"id":161,"date":"2024-06-10T13:54:20","date_gmt":"2024-06-10T13:54:20","guid":{"rendered":"https:\/\/sabba.me\/nmr\/?p=161"},"modified":"2024-06-10T13:58:41","modified_gmt":"2024-06-10T13:58:41","slug":"evolution-constrained-bounds-on-polarization-transfer-in-i_ns-spin-systems-part-i-refocused-inept","status":"publish","type":"post","link":"https:\/\/sabba.me\/nmr\/2024\/06\/10\/evolution-constrained-bounds-on-polarization-transfer-in-i_ns-spin-systems-part-i-refocused-inept\/","title":{"rendered":"Evolution-Constrained Bounds on Polarization Transfer in I_{N}S Spin Systems Part I: Refocused INEPT"},"content":{"rendered":"\n<p>It is seldom appreciated that the bounds on spin-order transfer are <em>evolution-constrained<\/em>. That is to say, a particular <em>type<\/em> of evolution (corresponding to some pulse sequence, effective Hamiltonian, or equivalent definition) &#8211; designed to convert some well-defined initial spin order configuration (be it a population, coherence, or arbitrary mixture thereof) into another well-defined target &#8211; may be associated with a distinct upper bound on the maximum achievable spin order transfer, an upper value that is in general <em>lower <\/em>than the <em>absolute <\/em>upper bound, irrespective of whether the absolute upper bound happens to be symmetry-constrained [<a href=\"https:\/\/www.sciencedirect.com\/science\/article\/pii\/S1090780715001950\" data-type=\"link\" data-id=\"https:\/\/www.sciencedirect.com\/science\/article\/pii\/S1090780715001950\">1<\/a>] or not [<a href=\"https:\/\/www.sciencedirect.com\/science\/article\/abs\/pii\/002223649090278H\" data-type=\"link\" data-id=\"https:\/\/www.sciencedirect.com\/science\/article\/abs\/pii\/002223649090278H\">2<\/a>].<\/p>\n\n\n\n<p>Oof, that was a mouthful. Here&#8217;s what I&#8217;m talking about (a table from a paper I&#8217;ve been writing for a few forevers):<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"803\" height=\"747\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-14.png\" alt=\"\" class=\"wp-image-179\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-14.png 803w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-14-300x279.png 300w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-14-768x714.png 768w\" sizes=\"auto, (max-width: 803px) 100vw, 803px\" \/><\/figure>\n\n\n\n<p>For now, let&#8217;s talk about just one of the columns from my table. The refocused INEPT pulse sequence is perhaps the most commonly used polarization transfer experiment in solution-state NMR spectroscopy. The experiment was described by Burum and Ernst <a href=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(80)90168-7\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(80)90168-7\">[3]<\/a> as a sequel to the famous INEPT sequence devised by Morris and Freeman [<a href=\"https:\/\/pubs.acs.org\/doi\/pdf\/10.1021\/ja00497a058\" data-type=\"link\" data-id=\"https:\/\/pubs.acs.org\/doi\/pdf\/10.1021\/ja00497a058\">4<\/a>].<\/p>\n\n\n\n<p>To my surprise, most NMR spectroscopists are unaware of the simple general expressions for both the optimal refocusing duration and the ensuing maximum enhancement factor, despite the fact that they were derived by David Doddrell and coworkers in the land <a href=\"https:\/\/www.youtube.com\/watch?v=XfR9iY5y94s&amp;ab_channel=MenAtWorkVEVO\" data-type=\"link\" data-id=\"https:\/\/www.youtube.com\/watch?v=XfR9iY5y94s&amp;ab_channel=MenAtWorkVEVO\">Down Under<\/a> in the 1980s [<a href=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(81)90186-4\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(81)90186-4\">5<\/a>, <a href=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1063\/1.444188\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1063\/1.444188\">6<\/a>, <a href=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(83)90303-7\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/10.1016\/0022-2364(83)90303-7\">7<\/a>, <a href=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1016\/0022-2364(85)90076-9\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1016\/0022-2364(85)90076-9\">8<\/a>, <a href=\"https:\/\/sci-hub.st\/10.1016\/S0065-3055(08)60611-5\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/10.1016\/S0065-3055(08)60611-5\">9<\/a>, <a href=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1016\/0022-2364(83)90259-7\" data-type=\"link\" data-id=\"https:\/\/sci-hub.st\/https:\/\/doi.org\/10.1016\/0022-2364(83)90259-7\">10<\/a>]. To save the reader time ruffling through references, these expressions can be found in a very brief JACS paper [<a href=\"https:\/\/pubs.acs.org\/doi\/pdf\/10.1021\/ja00393a067\" data-type=\"link\" data-id=\"https:\/\/pubs.acs.org\/doi\/pdf\/10.1021\/ja00393a067\">11<\/a>]. <\/p>\n\n\n\n<p>Let&#8217;s start with the INEPT sequence, which has a very simple theory in I<sub>N<\/sub>S spin systems. In a nutshell, INEPT accomplishes the conversion I<sub>z<\/sub> \u2192 [I<sub>x<\/sub> \u2192 2I<sub>y<\/sub>S<sub>z<\/sub> \u2192 2I<sub>z<\/sub>S<sub>z<\/sub> \u2192] 2I<sub>z<\/sub>S<sub>y<\/sub> (from pure I-spin longitudinal order to maximally I-correlated S-spin order, via a zz- or two-spin longitudinal- order intermediate), via the following steps:<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>The I-spins are tipped to the transverse plane by a 90<sub>y<\/sub> pulse i.e. something like I<sub>z<\/sub> \u2192 I<sub>x<\/sub>. [Alternatively, from another frame of reference, the <em>Hamiltonian <\/em>is tipped from 2\u03c0J\u00d7I<sub>z<\/sub>S<sub>z<\/sub> \u2192 2\u03c0J\u00d7I<sub>x<\/sub>S<sub>z<\/sub>].<\/li>\n\n\n\n<li class=\"has-medium-font-size\">Evolution for a total time period of \u03c4 = 1\/(2J) leads to a total rotation \u03c4\u00d72\u03c0J = \u03c0 about the I<sub>x<\/sub>S<sub>z<\/sub> axis (more properly \u03c0\/2 about the 2I<sub>x<\/sub>S<sub>z<\/sub> axis). This basically accomplishes the transformation I<sub>z<\/sub> \u2192 2I<sub>y<\/sub>S<sub>z<\/sub>.<\/li>\n\n\n\n<li class=\"has-medium-font-size\">The I spins are tipped to the longitudinal plane by a 90<sub>x<\/sub> pulse; i.e.  2I<sub>y<\/sub>S<sub>z<\/sub> \u2192 2I<sub>z<\/sub>S<sub>z<\/sub><\/li>\n\n\n\n<li>The S spins are tipped to the transverse plane for detection by 90<sub>-x<\/sub> pulse;  i.e.  2I<sub>z<\/sub>S<sub>z<\/sub> \u2192 2I<sub>z<\/sub>S<sub>y<\/sub>, leading to an antiphase multiplet most people associate with INEPT, that you can read more about in my paper [<a href=\"https:\/\/arxiv.org\/abs\/2404.03560\" data-type=\"link\" data-id=\"https:\/\/arxiv.org\/abs\/2404.03560\">12<\/a>].  <\/li>\n<\/ol>\n\n\n\n<p>[Note: there are some sign changes involved depending on the sign of the J-coupling, whether a 180 degree pulse is used, its phase, etc. which I will not waste our time worrying about.]<\/p>\n\n\n\n<p>The key thing to note about the INEPT sequence is that the I-spins are evolving under the action of the S-spins, of which &#8211; <em>by definition<\/em> &#8211; there is only one in I<sub>N<\/sub>S spin systems. Hence, the optimal duration of the INEPT sequence is always going to be 1\/(2J). It&#8217;s a no-brainer. <\/p>\n\n\n\n<p>Now, let&#8217;s consider the eponymous refocusing block present in refocused INEPT. In a nutshell, the refocusing block tries to accomplish the transformation 2I<sub>z<\/sub>S<sub>z<\/sub> [\u2192 2I<sub>z<\/sub>S<sub>y<\/sub> \u2192 S<sub>x<\/sub>] \u2192 S<sub>z<\/sub>. Note the symmetry with the INEPT sequence; it&#8217;s mirrored. Due to the fact that the S-spin is evolving under the action of (in general) multiple I-spins, things get a bit tricky, but it&#8217;s nothing we can&#8217;t handle. The problem essentially amounts to solving:<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"444\" height=\"109\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-2.png\" alt=\"\" class=\"wp-image-166\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-2.png 444w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-2-300x74.png 300w\" sizes=\"auto, (max-width: 444px) 100vw, 444px\" \/><\/figure>\n\n\n\n<p>Where:<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"47\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3-1024x47.png\" alt=\"\" class=\"wp-image-167\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3-1024x47.png 1024w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3-300x14.png 300w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3-768x35.png 768w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3-1536x70.png 1536w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-3.png 1690w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"753\" height=\"52\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-4.png\" alt=\"\" class=\"wp-image-168\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-4.png 753w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-4-300x21.png 300w\" sizes=\"auto, (max-width: 753px) 100vw, 753px\" \/><\/figure>\n\n\n\n<p>All we have to do is evaluate the first few solutions. We get:<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"707\" height=\"292\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-7.png\" alt=\"\" class=\"wp-image-172\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-7.png 707w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-7-300x124.png 300w\" sizes=\"auto, (max-width: 707px) 100vw, 707px\" \/><\/figure>\n\n\n\n<p>Now, the solution is staring us right in the face. It is easy to see by mathematical induction that we have a general solution for the transformation amplitude of longitudinal 2I<sub>z<\/sub>S<sub>z<\/sub> to S<sub>z<\/sub>:<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"747\" height=\"57\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-8.png\" alt=\"\" class=\"wp-image-173\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-8.png 747w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-8-300x23.png 300w\" sizes=\"auto, (max-width: 747px) 100vw, 747px\" \/><\/figure>\n\n\n\n<p>Obtaining the optimal duration of the refocusing tau (aka the time at which the maximum occurs) is a simple matter of high-school calculus:<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"312\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-10-1024x312.png\" alt=\"\" class=\"wp-image-175\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-10-1024x312.png 1024w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-10-300x91.png 300w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-10-768x234.png 768w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-10.png 1163w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<p>Leading to:<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"874\" height=\"222\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-11.png\" alt=\"\" class=\"wp-image-176\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-11.png 874w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-11-300x76.png 300w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-11-768x195.png 768w\" sizes=\"auto, (max-width: 874px) 100vw, 874px\" \/><\/figure>\n\n\n\n<p>In the thermodynamic limit, we encounter yet another high-school calculus problem:<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"479\" height=\"240\" src=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-13.png\" alt=\"\" class=\"wp-image-178\" srcset=\"https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-13.png 479w, https:\/\/sabba.me\/nmr\/wp-content\/uploads\/2024\/06\/image-13-300x150.png 300w\" sizes=\"auto, (max-width: 479px) 100vw, 479px\" \/><\/figure>\n\n\n\n<p>The reader is encouraged to compare these values with the so-called &#8220;entropy bound&#8221; (a term I despise with a passion) and the thermodynamic limit of Sorensen&#8217;s bound [<a href=\"https:\/\/doi.org\/10.1016\/0022-2364(91)90095-B\" data-type=\"link\" data-id=\"https:\/\/doi.org\/10.1016\/0022-2364(91)90095-B\">13<\/a>].<\/p>\n","protected":false},"excerpt":{"rendered":"<p>It is seldom appreciated that the bounds on spin-order transfer are evolution-constrained. That is to say, a particular type of evolution (corresponding to some pulse sequence, effective Hamiltonian, or equivalent definition) &#8211; designed to convert some well-defined initial spin order configuration (be it a population, coherence, or arbitrary mixture thereof) into another well-defined target &#8211; [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[7],"tags":[],"class_list":["post-161","post","type-post","status-publish","format-standard","hentry","category-theory"],"_links":{"self":[{"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/posts\/161","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/comments?post=161"}],"version-history":[{"count":4,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/posts\/161\/revisions"}],"predecessor-version":[{"id":180,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/posts\/161\/revisions\/180"}],"wp:attachment":[{"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/media?parent=161"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/categories?post=161"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sabba.me\/nmr\/wp-json\/wp\/v2\/tags?post=161"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}