In this work, we consider a model of two qubits driven by coherent and incoherent time-dependent controls. The dynamics of the system is governed by a Gorini–Kossakowski–Sudarshan–Lindblad master equation, where coherent control enters into the Hamiltonian and incoherent control enters into both the Hamiltonian (via Lamb shift) and the dissipative superoperator. We consider two physically different classes of interaction with coherent control and study the optimal control problem of state preparation formulated as minimization of the Hilbert–Schmidt distance’s square between the final density matrix and a given target density matrix at some fixed target time. Taking into account that incoherent control by its physical meaning is a non-negative function of time, we derive an analytical expression for the gradient of the objective and develop optimization approaches based on adaptation for this problem of GRadient Ascent Pulse Engineering (GRAPE). We study evolution of the von Neumann entropy, purity, and one-qubit reduced density matrices under optimized controls and observe a significantly different behavior of GRAPE optimization for the two classes of interaction with coherent control in the Hamiltonian.

The GRadient Ascent Pulse Engineering (GRAPE) method is widely used for optimization in quantum control. GRAPE is gradient search method based on exact expressions for gradient of the control objective. It has been applied to various coherently controlled closed and open quantum systems. In this work, we adopt the GRAPE method for optimizing objective functionals in open quantum systems driven by both coherent and incoherent controls. In our case, a tailored or engineered environment acts on the controlled system as control via time-dependent decoherence rates γ i ( t ) or, equivalently, via spectral density n ω ( t ) of the environment. To develop the GRAPE approach for this problem, we compute gradient of various objectives for general N-level open quantum systems for the piecewise constant class of control. The case of a single qubit is considered in details and solved analytically. For this case, an explicit analytical expression for evolution and objective gradient is obtained via diagonalization of a 3 × 3 matrix determining the system’s dynamics in the Bloch ball. The diagonalization is obtained by solving a cubic equation via Cardano’s method. The efficiency of the algorithm is demonstrated through numerical simulations for the state-to-state transition problem and its complexity is estimated. Robustness of the optimal controls is also studied.

In this work, we consider an environment formed by incoherent photons as a resource for controlling open quantum systems via an incoherent control. We exploit a coherent control in the Hamiltonian and an incoherent control in the dissipator which induces the time-dependent decoherence rates γk(t) (via time-dependent spectral density of incoherent photons) for generation of single-qubit gates for a two-level open quantum system which evolves according to the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) master equation with time-dependent coefficients determined by these coherent and incoherent controls. The control problem is formulated as maximization of the objective functional, which is the sum of Hilbert-Schmidt norms between four fixed basis states evolved under the GKSL master equation with controls and the same four states evolved under the ideal gate transformation. The exact expression for the gradient of the objective functional with respect to piecewise constant controls is obtained. Subsequent optimization is performed using a gradient type algorithm with an adaptive step size that leads to oscillating behaviour of the gradient norm vs iterations. Optimal trajectories in the Bloch ball for various initial states are computed. A relation of quantum gate generation with optimization on complex Stiefel manifolds is discussed. We develop methodology and apply it here for unitary gates as a testing example. The next step is to apply the method for generation of non-unitary processes and to multi-level quantum systems.