We investigate the time-optimal solution of the selective control of two uncoupled spin 1/2 particles. Using the Pontryagin Maximum Principle, we derive the global time-optimal pulses for two spins with different offsets. We show that the Pontryagin Hamiltonian can be written as a one-dimensional effective Hamiltonian. The optimal fields can be expressed analytically in terms of elliptic integrals. The time-optimal control problem is solved for the selective inversion and excitation processes. A bifurcation in the structure of the control fields occurs for a specific offset threshold. In particular, we show that for small offsets, the optimal solution is the concatenation of regular and singular extremals.
1The PMP and the geometric techniques have been recently applied with success to a series of fundamental problems in quantum control, such as, to cite a few, the state to state transfer [33,34,35,36,37,38], the implementation of unitary gates [40,39,41,42,43], the simultaneous control of different systems [44,45,46,47,48] and the control of two-level quantum systems or spin systems in presence of relaxation [49,50,51,52,53,54,55,58,57,56,59].In magnetic resonance, a benchmark example for the selective control of spins is given by an inhomogeneous ensemble of spin 1/2 particles with different offsets [6,8]. In this paper, we propose to investigate the simplest selectivity problem, that is the simultaneous time-optimal control of two uncoupled spins by means of magnetic fields. The two spins are assumed to be initially at the thermal equilibrium state, i.e. the north pole of the Bloch sphere. We consider in this work both the selective excitation and inversion processes for which the goal is to steer one of the two spins towards the equator or the south pole, while bringing back the other to the initial state. We derive the global time-optimal solution with a constraint on the maximum available field intensity. For a large offset difference, the optimal pulse is regular of maximum intensity. We show the existence of a bifurcation for a specific offset threshold. For smaller offset difference, the optimal solution is the concatenation of regular and singular arcs, the singular solution corresponding to zero field.The article is organized as follows. In Sec. 2, we define the model system and we show how to apply the PMP in this case. Section 3 is dedicated to the presentation of the results. We derive the time-optimal solutions for the selective excitation and inversion of spins. We discuss how this minimum time varies as a function of the offset difference. A comparison with a direct numerical optimization is made in Sec. 5. A summary of the different results obtained and prospective views are presented in Sec. 6. Technical computations are reported in the appendices A, B and C.