Papers
Topics
Authors
Recent
Search
2000 character limit reached

Q-PulsePol: Robust Finite-Pulse Quantum Control

Updated 5 July 2026
  • Q-PulsePol is a pulse-sequence-based quantum control scheme designed to enable efficient electron-nuclear polarization transfer by applying a phase-corrected inversion pulse under realistic pulse durations.
  • It employs a bimodal Floquet framework to isolate desired double-quantum channels while suppressing competing pathways through restored quadrature and time-reversal symmetries.
  • Robust against pulse imperfections and power errors, Q-PulsePol achieves near-ideal polarization transfer efficiency, offering significant advantages for bulk hyperpolarization in solid-state systems.

Q-PulsePol is a pulse-sequence-based quantum control scheme for polarization transfer between electron and nuclear spins, introduced as a finite-pulse-corrected variant of PulsePol for settings such as nitrogen-vacancy centers. Its defining modification is a phase adjustment of the central inversion pulse that restores quadrature symmetry in the interaction-frame spin Hamiltonian when microwave pulses have finite duration rather than being near-ideal and instantaneous. In the formulation of "Quadrature-Symmetric PulsePol for Robust Quantum Control Beyond the Ideal Pulse Approximation," this restoration suppresses the unwanted quantum pathway, preserves the intended single-mode transfer channel, and improves polarization-transfer efficiency under realistic pulse constraints (Jhamnani et al., 6 Apr 2026).

1. Driven-spin model and the finite-pulse problem

The starting point is the driven-spin Hamiltonian, written in the rotating frame for the electron,

H(t)=ΩSz+AzSzIz+AxSzIx+ω0nIz+Hμw(t),H(t)=\Omega S_z + A_z S_z I_z + A_x S_z I_x + \omega_{0n} I_z + H_{\mu w}(t),

with

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].

Here, Ω\Omega is microwave detuning from the electron resonance, (Az,Ax)(A_z,A_x) are secular/pseudo-secular hyperfine couplings, ω0n\omega_{0n} is the nuclear Larmor frequency, and ω1(t)\omega_1(t) describes the finite-duration pulse envelope (Jhamnani et al., 6 Apr 2026).

Finite-duration pulses enter through a piecewise ω1(t)\omega_1(t) with rise time trt_r, plateau ω1\omega_1, and fall time tft_f. During these finite ramps, the internal Hamiltonian

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].0

is not negligible. This is the mechanism by which the ideal “instantaneous”-pulse approximation fails.

Passing to the microwave interaction frame with

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].1

and to the nuclear free-precession frame with

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].2

the remaining interaction-frame Hamiltonian is

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].3

where

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].4

This formulation makes explicit that finite-pulse imperfections are encoded in the time dependence of Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].5 rather than being a perturbative afterthought. A plausible implication is that sequence performance is governed less by nominal flip angles alone than by the symmetry properties of the full interaction-frame trajectory.

2. Bimodal Floquet structure and transfer channels

Because the sequence is periodic with cycle time Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].6 and modulation frequency Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].7, Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].8 is expanded as

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].9

Substitution into Ω\Omega0 gives a bimodal Floquet Hamiltonian,

Ω\Omega1

The dynamic content relevant to dynamic nuclear polarization is carried by the Ω\Omega2 terms (Jhamnani et al., 6 Apr 2026).

For these harmonics,

Ω\Omega3

The Ω\Omega4 channel corresponds to double-quantum (DQ) transfer and the Ω\Omega5 channel to zero-quantum (ZQ) transfer when combined with Ω\Omega6. Resonance, or the Hartmann–Hahn condition, occurs when

Ω\Omega7

The effective strengths of the DQ and ZQ processes at harmonic Ω\Omega8 are described by

Ω\Omega9

and

(Az,Ax)(A_z,A_x)0

These scaling factors are the central quantitative diagnostics: they identify whether the periodic control enforces a clean uni-modal transfer or leaks weight into the unwanted channel.

3. Symmetry conditions, channel selectivity, and failure of standard PulsePol

Pure DQ transfer requires

(Az,Ax)(A_z,A_x)1

and

(Az,Ax)(A_z,A_x)2

Pure ZQ transfer requires the same equations with one sign flipped. In both cases, (Az,Ax)(A_z,A_x)3 and (Az,Ax)(A_z,A_x)4 must be in exact (Az,Ax)(A_z,A_x)5 quadrature (Jhamnani et al., 6 Apr 2026).

Two Floquet symmetries enforce these relations at all (Az,Ax)(A_z,A_x)6. The first is quadrature symmetry,

(Az,Ax)(A_z,A_x)7

which enforces

(Az,Ax)(A_z,A_x)8

thereby eliminating one quantum pathway. The second is XY–time-reversal symmetry,

(Az,Ax)(A_z,A_x)9

which further enforces

ω0n\omega_{0n}0

The first symmetry provides channel selectivity; the second provides maximal strength.

The finite-pulse pathology of standard PulsePol is traced to a specific symmetry-breaking event. In standard PulsePol, the finite rise and fall of the central inversion pulse at phase ω0n\omega_{0n}1 (ω0n\omega_{0n}2) break quadrature symmetry: the ω0n\omega_{0n}3 trajectory is no longer a quarter-cycle shifted copy of ω0n\omega_{0n}4. This is the mechanism identified by the bimodal Floquet analysis for the deterioration in fidelity under realistic pulse shaping. The result is not merely a reduction in desired coupling, but the simultaneous growth of the competing ZQ pathway.

4. Sequence construction and the central phase flip

Q-PulsePol is obtained by a single, minimal modification: the phase of the central ω0n\omega_{0n}5 pulse is flipped from ω0n\omega_{0n}6 to ω0n\omega_{0n}7. With this change, and for arbitrary pulse duration, the sequence recovers

ω0n\omega_{0n}8

together with

ω0n\omega_{0n}9

In the paper’s formulation, these relations ensure perfect uni-modal DQ transfer with full scaling factor even for long pulses (Jhamnani et al., 6 Apr 2026).

For the dominant Hartmann–Hahn condition, the construction is given at ω1(t)\omega_1(t)0, with ω1(t)\omega_1(t)1. One full cycle consists of eight ω1(t)\omega_1(t)2 segments and one ω1(t)\omega_1(t)3 inversion in the middle.

Pulse Phase Operation
ω1(t)\omega_1(t)4 ω1(t)\omega_1(t)5 ω1(t)\omega_1(t)6, ω1(t)\omega_1(t)7
ω1(t)\omega_1(t)8 ω1(t)\omega_1(t)9 ω1(t)\omega_1(t)0, ω1(t)\omega_1(t)1
ω1(t)\omega_1(t)2 ω1(t)\omega_1(t)3 ω1(t)\omega_1(t)4, ω1(t)\omega_1(t)5
ω1(t)\omega_1(t)6 ω1(t)\omega_1(t)7 ω1(t)\omega_1(t)8, ω1(t)\omega_1(t)9
trt_r0 trt_r1 trt_r2, trt_r3
trt_r4 trt_r5 trt_r6, trt_r7
trt_r8 trt_r9 ω1\omega_10, ω1\omega_11
ω1\omega_12 ω1\omega_13 ω1\omega_14, ω1\omega_15
ω1\omega_16 ω1\omega_17 ω1\omega_18, ω1\omega_19

The pulse durations are

tft_f0

and the free delays tft_f1 are chosen such that the total cycle time tft_f2 satisfies

tft_f3

with the typical choice tft_f4. In operational terms, Q-PulsePol is therefore not a wholesale redesign of PulsePol, but a symmetry-restoring phase schedule.

5. Polarization-transfer dynamics

For a single tft_f5–tft_f6 pair under perfect uni-modal DQ transfer at harmonic tft_f7, the net flip-flop rotation per cycle is

tft_f8

After tft_f9 repeated cycles, the polarization-transfer efficiency is

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].00

In the small-Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].01 regime,

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].02

whereas at long contact times the transfer saturates at Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].03 (Jhamnani et al., 6 Apr 2026).

This expression ties the Floquet analysis directly to the experimentally relevant observable. The function of Q-PulsePol is to preserve Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].04 near its ideal value while keeping the ZQ channel negligible. In that sense, symmetry restoration is not an abstract property of the control frame; it sets the net rotation angle accumulated per cycle and therefore the attainable build-up of nuclear polarization.

The emphasis on single-mode transfer is particularly important for bulk hyperpolarization in solids. The abstract characterizes Q-PulsePol as a practical and reliable scheme for bulk hyperpolarization of nuclear spins in solids using a single-mode (zero-quantum or double-quantum) transfer. This suggests that channel purity, rather than only raw coupling amplitude, is part of the sequence’s practical utility.

6. Robustness, design rules, and high-field relevance

Pulse amplitude errors Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].05 and resonance offset errors Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].06 enter through distortions of the trajectories Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].07 and Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].08 and therefore modify Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].09 and Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].10. Under finite-pulse conditions, standard PulsePol is reported to degrade strongly when

Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].11

is close to Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].12–Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].13: finite-pulse effects suppress Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].14 by Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].15–Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].16 and allow Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].17 to grow. By contrast, Q-PulsePol maintains Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].18 within Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].19–Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].20 of the ideal value down to Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].21 and enforces Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].22 (Jhamnani et al., 6 Apr 2026).

The numerical comparison given in the source is specific. For Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].23 amplitude errors or Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].24 detuning around Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].25 and Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].26, Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].27 remains above Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].28 of its ideal value for Q-PulsePol, while standard PulsePol can drop below Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].29. These figures summarize the practical consequence of the symmetry analysis: realistic microwave-power limitations and pulse ramps no longer force a large trade-off between selectivity and efficiency.

The design rules stated for high-field NV centres are correspondingly direct:

  • Quadrature symmetry: choose phases Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].30 such that Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].31; the minimal central inversion-pulse flip Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].32 suffices.
  • Harmonic choice: work at Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].33 unless the NV–Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].34 coupling demands another harmonic; ensure Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].35.
  • Finite-power regime: use Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].36 to reduce cycle-to-cycle cross-talk, but Q-PulsePol remains robust even at Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].37.
  • Bulk hyperpolarization: for bulk hyperpolarization via spin diffusion one needs strict uni-modal transfer; Q-PulsePol provides that even under realistic pulse rise/fall times.
  • Inhomogeneity compensation: compensate static microwave-power gradients by composite Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].38 pulses or by adding small phase-alternating supercycles; the quadrature constraint must still be preserved.
  • High-field trade-off: at high Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].39 one can trade longer Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].40 (lower Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].41) for lower Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].42; Q-PulsePol will retain its selectivity.
  • Calibration: calibrate Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].43, Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].44, and Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].45 by measuring the offset-profile and build the DNP build-up curve; the maximal enhancement occurs when Hμw(t)=ω1(t)[Sxcosϕ(t)+Sysinϕ(t)].H_{\mu w}(t)=\omega_1(t)[S_x\cos\phi(t)+S_y\sin\phi(t)].46 is near its ideal value, checked by observing the width of the resonance.

In this formulation, Q-PulsePol bridges idealized quantum control with realistic pulse engineering and establishes design rules for spin-based quantum control protocols. A plausible implication is that the sequence’s broader significance lies in showing how finite-bandwidth control errors can be corrected by restoring exact Floquet symmetries, rather than only by increasing microwave power or shortening pulses.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (1)

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Q-PulsePol.