The robustness of composite pulses elucidated by classical mechanics. II. The role of initial state imperfection

Here is an explanation of the paper using simple language, analogies, and metaphors.

The Big Picture: Herding Cats with a Magnet

Imagine you are trying to herd a group of cats (these are the atoms in your experiment) into a specific room (a state called "population inversion"). In the world of Nuclear Magnetic Resonance (NMR), scientists use magnetic "whips" called pulses to flip these cats from one side of the room to the other.

Ideally, you want every single cat to flip perfectly at the exact same time. But in the real world, things go wrong in two ways:

The Whip is Wobbly: The magnetic field isn't perfectly even (some cats get a stronger whack than others).
The Cats are Scattered: The cats don't all start in the exact same spot; they are scattered around a little bit.

For decades, scientists have built special "Composite Pulses" (sequences of whips) to fix the wobbly whip problem. They figured out how to make the sequence robust against the uneven field.

This paper asks a new question: What happens if the cats are scattered to begin with? Does the special whip sequence still work if the herd isn't perfectly lined up?

The Setup: The "Levitt" Sequence

The authors focus on a famous, clever sequence of pulses invented by a scientist named Levitt. It's like a specific dance routine: Step Left, Spin Around, Step Left again.

Previous Work: They already proved this dance works great if the field is wobbly, assuming all cats start in a perfect line.
Current Work: They asked, "What if the cats are spread out in a circle before we start dancing? Does the dance still keep them together?"

The Tools: Measuring the "Blob"

To answer this, the authors needed a way to measure how "messy" the group of cats gets. They used two main tools:

1. The Macroscopic View: The "Shadow Area"

Imagine shining a light on your group of cats from above. You see a shadow on the floor.

The Goal: You want the shadow to stay the same size. If the shadow gets huge, the cats have spread out too much, and the pulse failed to keep them synchronized.
The Problem: In physics, you can't actually shrink the total volume of a system (like a balloon that can't get smaller). However, the shadow (the 2D projection) can get bigger if the balloon gets squashed and stretched in weird ways.
The Finding: They found that Levitt's dance keeps the shadow relatively small, even with scattered cats. It's a very good dance.

2. The Microscopic View: The "Shear Coefficients"

This is like looking at the individual cats to see how they are being twisted.

Imagine a deck of cards. If you push the top of the deck sideways, the cards slide over each other. This is called shear.
The authors calculated how much the "cards" (the atoms) were sliding past each other.
The Connection: They found a mathematical rule (the "Co-area formula") that links the sliding of individual cards to the size of the shadow. If the cards slide too much, the shadow gets bigger.

The Results: How Good is Levitt's Dance?

The authors tested Levitt's sequence against two types of chaos:

Scenario A: The Wobbly Field (RF Inhomogeneity)

What happened: The cats were scattered, and the magnetic field was uneven.
Result: Levitt's dance worked surprisingly well! The "shadow" only grew by about 20%. The cats stayed mostly together.
The Twist: The authors found they could tweak the dance slightly (by tilting the direction of the second spin) to make the shadow even smaller. They found a "super-version" of Levitt's dance that is even more robust.

Scenario B: The Wrong Tune (Resonance Offset)

What happened: The cats were scattered, and the magnetic field was slightly off-tune (like trying to dance to music that is slightly too fast or slow).
Result: Levitt's dance was still good, but the cats got a bit more scattered than in the first scenario. The shadow grew by about 8%.
The Verdict: Interestingly, for this specific problem, Levitt's original dance is already so good that the authors couldn't find a "better" version. It was already near-perfect.

The Conclusion: Why This Matters

New Error Source: This paper highlights that "starting in the wrong place" is a systematic error that scientists need to worry about, not just "bad equipment."
Levitt is a Hero: The famous Levitt pulse sequence is robust even when the starting conditions are messy. It's a very efficient dance.
A New Toolkit: The authors created a new way to measure "robustness" by looking at both the big picture (shadow size) and the small picture (sliding cards). This toolkit can be used to design better pulses for quantum computers and medical imaging in the future.

The Takeaway Metaphor

Think of Levitt's pulse sequence as a highly skilled conductor leading an orchestra.

Old View: We knew the conductor could handle it if the musicians played slightly out of tune (field inhomogeneity).
New View: This paper checked if the conductor could handle it if the musicians were sitting in the wrong seats (initial state imperfection).
Result: The conductor is amazing! The orchestra stays in sync even with the wrong seats. However, the authors found that if the conductor adjusts his baton angle just a tiny bit, he can get the orchestra to play even more perfectly.

In short: The paper proves that a classic magnetic pulse technique is incredibly resilient, but with a little bit of math and optimization, we can make it even better for real-world, messy experiments.

Here is a detailed technical summary of the paper "The robustness of composite pulses elucidated by classical mechanics. II. The role of initial state imperfection" by Jonathan Berkheim and David J. Tannor.

1. Problem Statement

Composite Pulses (CPs) are standard tools in Nuclear Magnetic Resonance (NMR) and quantum control used to correct for systematic pulse imperfections, such as RF field inhomogeneity ( $\Omega_1$ ) and resonance offsets ( $\Delta$ ). While previous literature extensively addresses these Hamiltonian imperfections, the imperfection of the initial state (a spread of initial conditions, or ICs, on the Bloch sphere) has not been treated as a distinct systematic error.

In many physical scenarios (e.g., solid-state NMR with lattice distortions or thermal ensembles), the system does not start from a single pure state but from a 2D distribution of states. The authors aim to:

Analyze how CPs perform when acting on a 2D ensemble of initial conditions rather than a single point.
Determine if the robustness of established sequences (specifically Levitt's $90(x)180(y)90(x)$) is maintained when the initial state is distributed.
Develop a framework to quantify the "robustness" of a pulse sequence in terms of preserving the projected area of the ensemble while achieving population inversion.

2. Methodology

The authors extend their previous classical canonical framework to analyze the evolution of 2D distributions of initial conditions under pulse imperfections.

A. Theoretical Framework

Canonical Coordinates: The system is analyzed using canonical coordinates $(\phi, \eta) = (\tan^{-1}(y/x), z)$ on the Bloch sphere.
Liouville's Theorem: The authors note that while the extended phase space volume (including the imperfection parameter $w$ as a coordinate) is preserved under Hamiltonian flow (symplectic transformation), the projected 2D area of the ensemble on the $(\phi, \eta)$ plane is generally not preserved. It can expand or contract depending on the transformation.
Dual-Scale Evaluation: Robustness is assessed using two complementary metrics:
1. Macroscopic Measure (Projected Area): The total area $A$ covered by the union of distributions resulting from all possible imperfection values. The goal is to minimize the expansion ratio $R_{fi} = A_f / A_i$ .
2. Microscopic Measure (Shear Coefficients): Based on the stability matrix and the Co-area formula. The authors calculate the Jacobian of the transformation from the 3D space $(\phi_i, \eta_i, w)$ $(ϕ_{i}, η_{i}, w)$ to the 2D final space $(\phi_f, \eta_f)$ $(ϕ_{f}, η_{f})$ . The shear coefficient $G_{fi} = \sqrt{\det(J J^T)}$ $G_{f i} = det (J J^{T})$ quantifies local volume scaling.
  - $G_{fi} \gg 1$ : Expansion.
  - $G_{fi} \approx 1$ : Preservation.
  - $G_{fi} < 1$ (in specific regimes with multiplicity): Contraction.

B. Numerical Optimization

Case Study: Levitt's $90(x)180(y)90(x)$ pulse sequence.
Parameters: The study scans a 3D parameter space ( $\Delta\phi_2, \Delta\vartheta_2, \tau_1$ ) to find variants of Levitt's sequence that minimize the global area expansion ratio ( $R_{30}$ ) while maintaining high population inversion ( $\bar{\eta}_3 \to -1$ ).
Constraints: The optimization assumes symmetric 3-segment pulses and fixes the first rotation axis to the x-axis.

3. Key Results

A. RF Field Inhomogeneity Ensemble

Levitt's Performance: The sequence is robust in the $\eta$ $η$ -direction but causes significant expansion in the $\phi$ $ϕ$ -direction.
- Global Expansion: The projected area expands by 20% ( $R_{30} \approx 1.20$ ).
- Inversion Efficiency: High, with mean terminal population $\bar{\eta}_3 = -0.94$ (compared to $-0.99$ for a single IC).
- Mechanism: The shear coefficients show expansion in the first segment and preservation in the subsequent segments.
Optimization: By introducing a z-component to the rotation axes (tilting the pulses), the authors found variants that reduce the global expansion to $R_{30} \approx 1.08$ while maintaining comparable inversion efficiency. This suggests Levitt's sequence can be improved for inhomogeneous fields by adjusting the rotation axes out of the equatorial plane.

B. Resonance Offset Ensemble

Levitt's Performance: The sequence exhibits alternating expansion and contraction.
- Global Expansion: The area expands by only 8% ( $R_{30} \approx 1.08$ ).
- Inversion Efficiency: Moderate, with $\bar{\eta}_3 = -0.79$ .
- Mechanism: The sequence is robust in the $\phi$ -direction but less so in the $\eta$ -direction. The shear coefficients show sharp concentration near unity (contraction) in the final segments, which counteracts earlier expansion.
Optimization: Numerical scans revealed no significant improvement over Levitt's original sequence for this specific error type. The authors conclude that Levitt's sequence is already near-optimal for resonance offset ensembles.

C. Validation on Other Sequences

The method was successfully applied to Tycko's $180(0)180(120)180(0) $sequence, which was found to outperform Levitt's sequence for RF inhomogeneity ($ R_{30} = 1.05 $,$ \bar{\eta}_3 = -0.95$), validating the framework's utility for comparing different CP designs.

4. Key Contributions

Extension of Robustness Analysis: The paper formally treats the initial state distribution as a systematic error, extending stability analysis from 0D (single point) to 2D distributions on the Bloch sphere.
Dual-Scale Framework: Introduction of a unified evaluation method combining macroscopic area ratios and microscopic shear coefficients, linked theoretically by the Co-area formula. This provides a geometric explanation for why area expansion occurs even when Liouville's theorem holds in the extended space.
Geometric Insight into Levitt's Pulse: The study provides a geometric justification for the robustness of Levitt's pulse sequence, showing it maintains robustness even with IC spreads, though it is not perfect.
Optimization Strategy: Demonstration that adding a z-component to rotation axes (tilting pulses) can significantly improve robustness against RF inhomogeneity for ensembles, a modification not previously highlighted for this specific error type.

5. Significance

This work bridges the gap between theoretical stability analysis and practical experimental constraints where initial state inhomogeneity is unavoidable. By redefining robustness as the preservation of projected area rather than just point-wise fidelity, the authors provide a more rigorous metric for evaluating composite pulses in real-world ensembles. The findings suggest that while Levitt's sequence is remarkably robust, specific geometric modifications can yield further improvements, particularly for RF field inhomogeneity. The proposed methodology is general and applicable to any pulse sequence, offering a powerful tool for the design of high-fidelity quantum control protocols in NMR and quantum information processing.