Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 Aug 29;11(35):eadx7094.
doi: 10.1126/sciadv.adx7094. Epub 2025 Aug 29.

Robust bilinear rotations

Affiliations

Robust bilinear rotations

Yannik T Woordes et al. Sci Adv. .

Abstract

Bilinear rotation elements allow interaction-dependent manipulations of spins in quantum technologies and, particularly, in spectroscopy. We used bilinear rotation in nuclear magnetic resonance (NMR) spectroscopy and derived several ways of introducing robustness into the filter element. We distinguished two performance levels: the coupling dependent inversion of polarization and the full bilinear π-rotation capability. In addition, all four essential variants of bilinear π rotations-BIRDd, BIRDr, BIRDd,X, and BIRDr,X filters-were given for all cases studied. In the first step, adiabatic CHIRP-type and BUBI/BUBU pulses lead to improved robustness with respect to offset/detuning effects and variations in B1-field strengths. In the second step, we optimized time-optimal coupling-compensated BIRD elements. Together with correspondingly optimized pulse shapes, we established fully coupling, offset, and B1-compensated bilinear π-rotation (COB-BIRD) elements and characterized them in theory and experiment. Overall, this demonstrated the use of the robust bilinear rotation capability on a partially aligned sample in a homodecoupled two-dimensional NMR experiment.

PubMed Disclaimer

Figures

Fig. 1.
Fig. 1.. Theoretical comparison of the BIRD and JC-BIRD pulse sequences on their efficacy of spin inversion and universal rotation to a variation in J coupling.
Original BIRD (A) and JC-BIRD (B) pulse sequences (4, 5) and their robustness against J coupling variations (C and D). Solid and open rectangles correspond to 90° and 180° pulses, respectively; phases are x unless indicated otherwise, and transfer delays are calibrated to match Δ = 1/(2 × 1JCH). The phases ϕ1 to ϕ4 of the BIRD sequence (A) can be manipulated according to Table 1 to obtain the four possible effective bilinear rotations. If ϕ4 is not indicated, the dashed pulse is left out. The ΦPP (C) and ΦUR* (D) fidelities of both the BIRDd,X element (solid blue, ϕ3 = x, no dashed 180° pulse applied) and the JC-BIRDd,X element (dashed orange) are simulated using on-resonant hard pulses and Δ = 1/(2 × 185 Hz).
Fig. 2.
Fig. 2.. The use of offset-compensated adiabatic or optimal control–derived pulses in the BIRD and JC-BIRD to extend the robustness into offset variation.
(A) Alternative pulse sequence for BIRD elements based on a pair of adiabatic pulses to achieve offset-compensated optimal transfer of the BIRDd bilinear rotation. Trapezoids correspond to adiabatic pulses with their sweep direction indicated by the arrows. Universal rotation refocusing (UR-180°) and point-to-point inversion pulses [PP(z-z)] are presented as open rectangles with a dash or a curved diagonal, respectively. (B) A JC-BIRD sequence using a pair of BUBI and a single BUBU pulse sandwich for offset compensation is shown. Transfer delays are calibrated to match Δ = 1/(2 × 1JCH). (C, E, G, and I) The simulations of ΦPP (left column) and (D, F, H, and J) ΦUR* (right column) of the BIRD elements with respect to the offset frequency vS and J coupling for hard-pulse BIRDd [(C) and (D)] and JC-BIRDd [(E) and (F)], as well as the offset-compensated adiabatic BIRDd [(G) and (H)] and BUBI/BUBU-JC-BIRDd [(J) and (K)]. Hard pulses were simulated with an rf amplitude of 25 kHz, corresponding to a 10-μs 90° pulse. Adiabatic WURST40 pulses with T = 1 ms, Q = 5, and a sweep-width ΔvS = 40 kHz and BUBI/BUBU pulse sandwiches as specified in Table 2 were used in the simulations. Contour levels in all cases equal 0.8 (indigo), 0.9 (orange), 0.94 (ultramarine), and 0.97 (magenta).
Fig. 3.
Fig. 3.. TOP performance curves for optimal control–derived JC-BIRD elements.
Continuous shapes as well as hard-pulse delay sequences with, neglecting 180° pulses in refocused delays, (n + 1) proton pulses and (n) delays were evaluated using ΦUR and taking only the on-resonance performance with respect to J into account. Different TOP curves are represented with different colors as specified on the top of the graph. While the shape TOP curve marks roughly the physical limit for the bilinear rotation, hard-pulse delay sequences are needed for actual robust implementations. The performances of the original BIRD (2p1d due to the two 90° pulses flanking the refocused delay) and JC-BIRD (4p3d) sequences as well as the 4p3d sequence used subsequently for the COB-BIRD implementation are annotated correspondingly.
Fig. 4.
Fig. 4.. Comparison of experimental and simulated data to the universal π rotation of the BIRDd and the COB-BIRDd.
Experimental demonstration (black spectral series) of the BIRDd (A) and the COB-BIRDd (B) sequences, as well as their comparison to corresponding simulations (superimposed red lines). Experimental data were obtained as described in the main text, and simulations show ΦPP(IxI) fidelities on the highlighted elements using on-resonant hard pulses. The COB-BIRDd pulse sequence is given in Fig. 5.
Fig. 5.
Fig. 5.. The four essential COB-BIRD pulse sequences and their experimental fidelities for spin inversion and universal π rotation against coupling and offset variation.
(A) The general sequence for the COB-BIRD. For effective bilinear π rotations in the coupling range of 120 to 250 Hz, all delays δ1 = 2.583 ms. Universal rotation refocusing (UR-180°) and point-to-point inversion pulses [PP(z-z)] of pulse sandwiches are applied as indicated. By replacing the gray box with either (B) or (C) and setting the flip angle ϕ5 accordingly, all four different COB-BIRD bilinear rotations can be implemented. In order to design the COB-BIRDr,X or the COB-BIRDd,X, fragment (B) is used with ϕ5 set to −148° or 32°, respectively, and in order to design the COB-BIRDd or COB-BIRDr, fragment (C) is used with ϕ5 set to −148° or 32°, respectively. As an example, the ΦPP (D) and ΦUR* (E) fidelities for the COB-BIRDr,X variant are simulated with respect to vS and J. Contours are scaled the same way as in Fig. 2.
Fig. 6.
Fig. 6.. Experimental comparison of various refocusing elements applied in a 13C-detected, 1H-decoupled, and refocused J-INEPT experiment.
(A) The corresponding sequence with a gray box for different inversion elements and a COB3-INEPT transfer step is shown together with the structure and aromatic numbering of nicotine on top of the spectra. All pulses are applied along x unless indicated otherwise, and solid and open bars refer to 90° and 180° pulses, respectively, unless a different flip angle (136°) is annotated. For the refocusing elements, a simple 180° hard-pulse spin echo (B) is compared to the BIRDd,X (C), the JC-BIRDd,X (D), and the COB-BIRDd,X (E) elements. All four 2D spectra were applied to a partially aligned sample of (−)-nicotine as explained in the main text. For each experiment, slices of the aromatic signals are extracted over the F1 dimension and shown with the corresponding measured T coupling and signal-to-noise-ratio (S/N). The transfer delays in the BIRDd,X and JC-BIRDd,X are set to 125 Hz to accommodate the small couplings, and the COB-BIRDd,X is applied as presented in Fig. 5. Gradients of 1-ms duration and 200-μs recovery delay are applied in percent of the maximum gradient strength of nominally 53.5 G/cm according to G1 = 7%, G2 = 25%, G3 = −G4 = 10%, G5 = 60%, and G6 = 80%. Delays are δg = 1.2 ms, δ2 = 0.5401 ms, δ3 = 1.065 ms, and δ4 = 1.0702 ms. For a detailed description including phase cycling and the particular 180° element used, refer to the Supplementary Materials (fig. S14).

References

    1. Pfeiffer F., Werninghaus M., Schweizer C., Bruckmoser N., Koch L., Glaser N. J., Huber G. B. P., Bunch D., Haslbeck F. X., Knudsen M., Krylov G., Liegener K., Marx A., Richard L., Romeiro J. H., Roy F. A., Schirk J., Schneider C., Singh M., Södergren L., Tsitsilin I., Wallner F., Riofrío C. A., Filipp S., Efficient decoupling of a nonlinear qubit mode from its environment. Phys. Rev. X 14, 041007 (2024).
    1. Google Quantum AI and Collaborators , Quantum error correction below the surface code threshold. Nature 638, 920 (2025). - PMC - PubMed
    1. A.M. Stephens, A.G. Fowler, L.C.L. Hollenberg, Universal fault tolerant quantum computation on bilinear nearest neighbor arrays. arXiv:quant-ph/0702201 (2008).
    1. Garbow J. R., Weitekamp D. P., Pines A., Bilinear rotation decoupling of homonuclear scalar interactions. Chem. Phys. Lett. 93, 504–509 (1982).
    1. Uhrín D., Liptaj T., Kövér K. E., Modified BIRD pulses and design of heteronuclear pulse sequences. J. Magn. Reson. A 101, 41–46 (1993).