Quantum Gate Generation in Two-Level Open Quantum Systems by Coherent and Incoherent Photons Found with Gradient Search

TL;DR

This paper develops a gradient-based optimization method to generate quantum gates in open two-level systems by combining coherent and incoherent controls, leveraging spectral density modulation of incoherent photons.

Contribution

It introduces a novel gradient calculation for control of open quantum systems with both coherent and incoherent controls, enabling efficient optimization of quantum gate generation.

Findings

Optimized control protocols for single-qubit gates in open systems.

Gradient expressions derived for piecewise constant controls.

Method applicable to non-unitary processes and multi-level systems in future work.

Abstract

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 $\gamma_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 minimization 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.

Equations

$$\frac{d\rho (t)}{dt}=-i[H_0+H_c(t),\rho (t)]+\sum _k{\gamma }_k(t){\mathcal{D}}_k(\rho (t))$$

$$\frac{d\rho }{dt}=-i[H_0+Vu(t),\rho ]+\gamma {\mathcal{L}}_{n(t)}(\rho ),\rho (0)={\rho }_0.$$

$$\begin{gathered} {\mathcal{L}}_{n(t)}({\rho }_t)=n(t)\left({\sigma }^{+}{\rho }_t{\sigma }^{-}+{\sigma }^{-}{\rho }_t{\sigma }^{+}-\frac{1}{2}\{{\sigma }^{-}{\sigma }^{+}+{\sigma }^{+}{\sigma }^{-},{\rho }_t\}\right) \\ +\left({\sigma }^{+}{\rho }_t{\sigma }^{-}-\frac{1}{2}\{{\sigma }^{-}{\sigma }^{+},{\rho }_t\}\right), \end{gathered}$$

$$\rho=\dfrac{1}{2}\Bigg{(}\mathbb{I}+\sum_{\alpha\in\{x,y,z\}}r_{\alpha}\sigma_{\alpha}\Bigg{)},\quad r_{\alpha}=\mathrm{Tr}\rho\sigma_{\alpha},\quad\alpha\in\{x,y,z\},$$

$$\frac{d\mathbf{r}}{dt}=A(t)\mathbf{r}+\mathbf{b}=\left(B+B^{u}u(t)+B^{n}n(t)\right)\mathbf{r}+\mathbf{b},\mathbf{r}(0)={\mathbf{r}}_0,$$

$$B=\left(\begin{array}{ccc}-\frac{\gamma }{2}& \omega & 0\\ -\omega & -\frac{\gamma }{2}& 0\\ 0& 0& -\gamma \end{array}\right),B^u=2\mu \left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -1\\ 0& 1& 0\end{array}\right),B^n=-\gamma \left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 2\end{array}\right),\mathbf{b}=\left(\begin{array}{c}0\\ 0\\ \gamma \end{array}\right).$$

$$n(\omega )=\frac{1}{{\pi }^2}\frac{{\omega }^3}{\exp (\beta \omega )-1},$$

$$\underset{0}{\overset{\infty }{\int }}{n}_{\omega }(t)d\omega <\infty \forall t\ge 0,$$

$$\rho (t,u,n)=\Phi (t,u,n){\rho }_0.$$

$$J(\rho ,{\rho }_{\text{target}})=\Vert \rho -{\rho }_{\text{target}}{\Vert }^2=\mathrm{Tr}(\rho -{\rho }_{\text{target}}{)}^2,$$

$$F(u,n)=J\left(\Phi (T,u,n){\rho }_0,{\rho }_{\text{target}}\right)=\Vert \Phi (T,u,n){\rho }_0-{\rho }_{\text{target}}{\Vert }^2\to \inf .$$

$$D_{{\rho }_0^{(1)},\dots ,{\rho }_0^{(N)}}(\mathcal{W},\mathcal{W}{}^{\prime })=\frac{1}{N}\sum _{j=1}^{N}J\left(\mathcal{W}{\rho }_0^{(j)},\mathcal{W}{}^{\prime }{\rho }_0^{(j)}\right).$$

$$\begin{gathered} F_{U,N}\left(u,n;{\rho }_0^{(1)},\dots ,{\rho }_0^{(N)}\right)=D_{{\rho }_0^{(1)},\dots ,{\rho }_0^{(N)}}(\Phi (T,u,n),U\cdot U^{†})= \\ \frac{1}{N}\sum _{j=1}^{N}J\left(\Phi (T,u,n){\rho }_0^{(j)},U{\rho }_0^{(j)}{U}^{†}\right)=\frac{1}{N}\sum _{j=1}^N{\Vert \rho (T,u,n)-U{\rho }_0^{(j)}{U}^{†}\Vert }^2\to \inf . \end{gathered}$$

$$\rho (t,u,n)$$

$$U{\rho }_0^{(j)}{U}^{†}$$

$$F_{U,N}\left(u,n;{\rho }_0^{(1)},\dots ,{\rho }_0^{(N)}\right)=\frac{1}{2N}\sum _{j=1}^N{\left|{\mathbf{r}}^{(j)}(T,u,n)-{\mathbf{r}}_{\text{target}}^{(j)}\right|}^2\to \inf .$$

$$0\le F_{U,N}\left(u,n,{\rho }_0^{(1)},\dots ,{\rho }_0^{(N)}\right)\le 2.$$

$$\begin{array}{rl}{\rho }_0^{(1)}& =|0⟩⟨0|=\left(\begin{array}{cc}1& 0\\ 0& 0\end{array}\right),\\ {\rho }_0^{(2)}& =|1⟩⟨1|=\left(\begin{array}{cc}0& 0\\ 0& 1\end{array}\right),\\ {\rho }_0^{(3)}& =|+⟩⟨+|=\frac{1}{2}\left(\begin{array}{cc}1& 1\\ 1& 1\end{array}\right),\\ {\rho }_0^{(4)}& =|i⟩⟨i|=\frac{1}{2}\left(\begin{array}{cc}1& -i\\ i& 1\end{array}\right).\end{array}$$

$$\begin{gathered} {\rho }_0=\left(\begin{array}{cc}{x}_1& x_3+i\text{leavevmode}{x}_4\\ x_3-i{x}_4& x_2\end{array}\right)=(x_1+x_4-x_3){\rho }_0^{(1)}+ \\ +(x_2+x_4-x_3){\rho }_0^{(2)}+2x_3{\rho }_0^{(3)}-2x_4{\rho }_0^{(4)}. \end{gathered}$$

$$\begin{array}{r}{\mathbf{r}}_0^{(1)}=(0,0,1),{\mathbf{r}}_0^{(2)}=(0,0,-1),{\mathbf{r}}_0^{(3)}=(1,0,0),{\mathbf{r}}_0^{(4)}=(0,1,0).\end{array}$$

$$\sqrt{Y}|0⟩⟨0|{\sqrt{Y}}^{†}$$

$$\sqrt{Y}|1⟩⟨1|{\sqrt{Y}}^{†}$$

$$n(t)=w(t{)}^2,$$

$$u(t)$$

$$w(t)$$

$$A_k=B+B^u{u}_k+B^n{w}_k^2,k=1,\dots ,M.$$

$$\begin{gathered} {\mathbf{r}}_k\equiv \mathbf{r}(t_k)=e^{A_k\Delta t_k}{\mathbf{r}}_{k-1}+{\mathbf{g}}_k=e^{A_k\Delta t_k}\cdots e^{A_1\Delta t_1}{\mathbf{r}}_0 \\ +e^{A_k\Delta t_k}\cdots e^{A_2\Delta t_2}{\mathbf{g}}_1+\cdots +e^{A_k\Delta t_k}{\mathbf{g}}_{k-1}+{\mathbf{g}}_k, \end{gathered}$$

$${\mathbf{g}}_k=(e^{A_k\Delta t_k}-\mathbb{I})A_k^{-1}\mathbf{b},\Delta t_k=t_k-t_{k-1},k=1,\dots ,M.$$

$$\dfrac{\partial\mathbf{r}(T)}{\partial v^{m}_{k}}=e^{A_{N}\Delta t_{N}}\dots e^{A_{k+1}\Delta t_{k+1}}\bigg{[}\dfrac{\partial}{\partial v_{k}^{m}}\left(e^{A_{k}\Delta t_{k}}\right)\mathbf{r}_{k-1}+\dfrac{\partial\mathbf{g}_{k}}{\partial(u_{k},w_{k})}\bigg{]},\\ k=1,\dots,M,\quad m=1,2.$$

$$\frac{\partial {\mathbf{g}}_k}{\partial v_k^1}=\frac{\partial {\mathbf{g}}_k}{\partial u_k}=\left(\frac{\partial }{\partial u_k}{e}^{A_k\Delta t_k}-(e^{A_k\Delta t_k}-\mathbb{I})A_k^{-1}{B}^u\right)A_k^{-1}\mathbf{b},$$

$$\frac{\partial }{\partial v_k^1}{e}^{A_k\Delta t_k}=\frac{\partial }{\partial u_k}{e}^{A_k\Delta t_k}=\underset{0}{\overset{\Delta t_k}{\int }}{e}^{A_{k}t}{B}^u{e}^{A_k(\Delta t_k-t)}\mathrm{d}t,$$

$$\frac{\partial {\mathbf{g}}_k}{\partial v_k^2}=\frac{\partial {\mathbf{g}}_k}{\partial w_k}=\left(\frac{\partial }{\partial w_k}{e}^{A_k\Delta t_k}-2w_k(e^{A_k\Delta t_k}-\mathbb{I})A_k^{-1}{B}^n\right)A_k^{-1}\mathbf{b},$$

$$\frac{\partial }{\partial v_k^2}{e}^{A_k\Delta t_k}=\frac{\partial }{\partial w_k}{e}^{A_k\Delta t_k}=2w_k\underset{0}{\overset{\Delta t_k}{\int }}{e}^{A_{k}t}{B}^n{e}^{A_k(\Delta t_k-t)}\mathrm{d}t,k=1,\dots ,M.$$

Full text

Keywords: quantum control; coherent control; incoherent control; open quantum systems; single-qubit gate; Hadamard gate.

1 Introduction

Optimal control of quantum systems has various applications ranging from quantum computing to optimal molecule discrimination **[1, 2, 3, 4, 5]**. Often controlled quantum systems are interacting with the environment. Hence, they have to be considered as open quantum systems. In some situations, the environment prevents from controlling the system. In other cases, it can be used as a useful resource. An example is via incoherent control proposed in **[6, 7]**, where time dependent decoherence rates induced by non-equilibrium spectral density of incoherent photons are generally used as control jointly with coherent control by lasers. This incoherent control uses the engineered environment as a useful resource. Such incoherent control was applied to general multilevel systems in the dynamical picture (i.e., using master equations for the description of the dynamics)  **[8]**, as well as to specific problems for single-qubit and two-qubit models, e.g., in **[9, 10, 11, 12]**, etc. Controllability of open quantum systems in the kinematic picture (using universally optimal Kraus maps for description of evolution of the system) was studied in **[13]**. The first experimental realization of Kraus maps studied in **[13]** for an open single qubit was done in **[14]**.

Engineered environments were used for manipulation of quantum systems in various contexts. Preparation of mixed states and generation of non-unitary evolutions are necessary for quantum computing with mixed states and non-unitary circuits **[15, 16]**. Engineered environments were suggested for improving quantum state engineering **[17]**, cooling of translational motion without reliance on internal degrees of freedom **[18]**, preparing entangled states **[19, 20]**, improving quantum computation **[17]**, inducing multiparticle entanglement dynamics **[21]**, and making robust quantum memories **[22]**. A scheme for steady-state preparation of a harmonic oscillator in the first excited state using dissipative dynamics was proposed **[23]**. Quantum harmonic oscillator state preparation by reservoir engineering was considered **[24]**. Optimal control for cooling a quantum harmonic oscillator by controlling its frequency was proposed **[25]**. Reachable states for a qubit driven by coherent and incoherent controls were analytically described using geometric control theory **[10]**. Controllability analysis of quantum systems immersed within an engineered environment was performed **[26]**. Coherent control and simplest switchable noise on a single qubit were shown to be capable of transferring between arbitrary $n$-qubit quantum states **[27]**. Quantum state preparation by controlled dissipation in finite time was proposed **[28]**. Precision limits to quantum metrology of open dynamical systems through environment control were studied **[29]**.

Various optimization methods are used to find controls for quantum systems including Krotov type approaches (**[30, 31, 32]**, Zhu–Rabitz **[33]** method, GRadient Ascent Pulse Engineering (GRAPE) **[34]**, Hessian based optimization in the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm and combined approaches **[35, 36]**, speed gradient method **[37]**, gradient free Chopped RAndom Basis (CRAB) optimization **[38]**, genetic algorithms **[39, 6]**, dual annealing **[11]**, machine learning **[40]**, etc. Quantum speed limit for a controlled qubit moving inside a leaky cavity has been studied **[41]**. Preparation of single qubit and two-qubit gates has been studied, e.g., in **[42, 43, 44, 32, 45, 46, 47, 48, 49]**, etc.

In this work, the incoherent control in the dissipator is implemented using the environment formed by incoherent photons and is used as a resource, together with the coherent control in the Hamiltonian, to control a single qubit system, based on the approach in **[6]**. The incoherent control induces time-dependent decoherence rates $\gamma_{k}(t)$ via time-dependent spectral density of the environment $n(t)$, in addition to the coherent control. The following master equation for the system density matrix $\rho(t)$ evolution was proposed **[6]**:

[TABLE]

Here, $H_{0}$ is a free system Hamiltonian, $H_{c}(t)$ is a Hamiltonian describing interaction of the system with coherent control, $k$ denotes all possible different pairs of energy levels in the system, and ${\cal D}_{k}$ is a Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) dissipator (we set Planck constant as $\hbar=1$). In addition to general consideration, two physical classes of the environment were exploited in **[6]** — incoherent photons and quantum gas, with two explicit forms of ${\cal D}_{k}$ derived in the weak coupling limit (describing atom interacting with photons) and low density limit (describing collisional type decoherence), respectively. In **[8]**, it was shown that for the master equation (1) with ${\cal D}_{k}$ derived in the weak coupling limit, a generic $N$-level quantum system becomes approximately completely controllable in the set of density matrices.

As particular control example, we consider here generation of the single-qubit gate $X$ (Pauli matrix $X$) and Hadamard gate $H$ for a two-level open quantum system, evolving according to the Gorini–Kossakowski–Sudarshan–Lindblad master equation with time-dependent coefficients determined by these coherent and incoherent controls. The control problem is formulated as the minimization of the objective functional, which is the sum of Hilbert–Schmidt distances between density matrices evolving under controls from some basis states for some time $T$ and target density matrices. The target density matrices are produced from the basis states using the ideal quantum gate. Note that using three specific mixed states is enough in the general case of a $N$-level quantum system **[32]**. We obtain an exact expression for the gradient of the objective functional with respect to piecewise constant controls and perform the optimization using a gradient type algorithm with an adaptive step size. The use of the adaptive step size leads to oscillating behavior of the norm of the gradient vs number of iterations, which is in contrast to Figure 2 in **[50]**. In that figure, the norm of the gradient of the objective first monotonically increases, and then, monotonically decreases with number of the iterations. We visualize the results by plotting in the Bloch ball evolution of the initial states under obtained controls which are in a good agreement with the ideal case.

We utilize the gradient descent method for optimization. In the context of quantum control, it is known as GRAPE. The gradient descent is a local method and does not necessarily converge to a global minimum. It can converge to a local minimum if such a local minimum exists. Therefore, for objective functionals with traps, i.e., local but not global minima, using slower global stochastic methods might be preferable. The analysis of traps for quantum systems was initiated in **[51]** and studied for various models in **[50, 52, 53, 54, 55, 56, 57, 58, 59, 60]**, etc.

The analysis of this work is related to the approach developed in **[61, 62]** to control of open quantum systems, where Kraus maps are described by points of a complex Stiefel manifold (strictly speaking, factor of the Stiefel manifold over some equivalence relation). This picture is called a kinematic picture for description of the dynamics. In these works, the gradient and Hessian based optimization methods were developed for optimization of control objectives for open quantum systems. In the present work, we consider the dynamical picture, where evolution of an open qubit is described by a master equation. The relation between the two approaches is that this master equation induces some time evolved Kraus map, which corresponds to some trajectory on a suitable factor of a Stiefel manifold. The gradient of the objective in the dynamic picture, which is used to find optimal controls, can be computed using chain rule for derivatives and gradient of the objective computed in the kinematic picture.

The structure of this paper is the following. In Section 2, the master equation for a qubit driven by coherent and incoherent photons is provided. In Section 3, incoherent control by photons is discussed in more details. In Section 4, the objective functional which describes the problem of generation of single-qubit quantum gates is formulated. In Section 5, the gradient of the objective is computed for piecewise constant controls. In Section 6, numerical optimization results are provided, including plots for optimal coherent and incoherent controls, and for the evolution of the qubit basis states in the Bloch ball under the action of the optimal controls. In Section 7, a connection of the obtained results with gradient optimization over Stiefel manifolds **[61, 62]** is discussed. The conclusions in Section 8 summarizes the results.

2 Master Equation for a Qubit Driven by Coherent and Incoherent Controls

We consider a master equation for an open two-level quantum system (qubit) driven by coherent and incoherent controls with time-dependent decoherence rate determined by $\gamma(t)=\gamma[n(t)+1/2]$, where $n(t)$ is the spectral density of incoherent photons surrounding the qubit at the time $t$:

[TABLE]

Here, $H_{0}=\omega\begin{pmatrix}0&0\\ 0&1\end{pmatrix}$ is the free Hamiltonian, $V=\mu\sigma_{x}=\mu\begin{pmatrix}0&1\\ 1&0\end{pmatrix}$ is the interaction Hamiltonian, $\omega$ is the transition frequency, $\mu>0$ is the dipole moment, $u(t)$ is coherent control (real-valued function), and $\gamma>0$ is the decoherence rate coefficient. The dissipative superoperator is

[TABLE]

where $n(t)\geq 0$ is incoherent control (non-negative real-valued function), matrices $\sigma^{\pm}=\dfrac{1}{2}(\sigma_{x}\pm i\sigma_{y})$, $\sigma^{+}=\begin{pmatrix}0&1\\ 0&0\end{pmatrix}$, $\sigma^{-}=\begin{pmatrix}0&0\\ 1&0\end{pmatrix}$, and $\sigma_{x}$, $\sigma_{y}$, $\sigma_{z}$ are the Pauli matrices.

We use convenient parameterization of the density matrix $\rho$ by Bloch vector $\mathbf{r}=(r_{x},r_{y},r_{z})\in\mathbb{R}^{3}$:

[TABLE]

where $\mathbb{I}$ is the $2\times 2$ identity matrix. In this representation, we obtain the following equation for the dynamics of the Bloch vector $\mathbf{r}(t)$:

[TABLE]

where $r_{0\alpha}=\mathrm{Tr}\rho_{0}\sigma_{\alpha}$, $\alpha\in\{x,y,z\},$ and

[TABLE]

3 Incoherent Control

Incoherent control as proposed in **[6]** is realized by using the environment as a resource to manipulate quantum systems. Physically incoherent control can be realized by using incoherent light and tailoring its distribution in frequency. In this case, control is density of photons $n_{\omega}(t)$ with frequency $\omega=|{\bf k}|c$ (${\bf k}$* is photon momentum and $c$ is the speed of light) and time $t$. In general, photons have polarization $\alpha$, and their proper distribution should include dependence on polarization, if any, so that incoherent control becomes a function of $\omega$, $\alpha$, and $t$. However, we neglect this dependence in the present work. For thermal light or black-body radiation, distribution of its photons has the famous Planck form*

[TABLE]

where $\beta$ is the inverse temperature (in the Planck system of units, where reduced Planck constant $\hbar$, speed of light $c$, and Boltzmann constant $k_{\textrm{B}}$ equal one). However, it can generally be non-thermal with other shape $n_{\omega}(t)$, which can evolve in time. One can obtain non-thermal shapes, e.g., by filtering of black-body radiation, using light-emission diodes, etc. A plot with Planck distribution and some of its filterings are shown in Figure 1.

Since physical meaning of incoherent control is density of particles of the environment, mathematically, it is a non-negative quantity, $n_{\omega}(t)\geq 0$. In general, one can require that

[TABLE]

meaning that total density of photons at any time moment is finite. Decoherence rates for a quantum system, e.g., an atom, interacting with incoherent light, are determined by strength of interaction between the atom and light, and by this density of photons. Such incoherent control naturally makes decoherence rates time dependent. The magnitude of the decoherence rates affects the speed of decay of off-diagonal elements of the system density matrix and determines the values of the diagonal elements towards which the diagonal part of the system density matrix evolves. This property can be used, for example, for approximate generation of various density matrices **[8]**.

In our case of a single qubit, the qubit only has one transition frequency $\omega_{0}$. In the first approximation, the spectral density of incoherent light at only this frequency, $n(t)=n_{\omega_{0}}(t)$, affects the system. Thus, in this case, without loss of generality, one can use only thermal radiation with variable temperature. For multilevel systems, this is not the case and thermal distribution with variable temperature does not provide the full generality for incoherent control at all system’s transition frequencies.

4 Objective Functional for Single-Qubit Gate Generation

Let us introduce the evolution superoperator $\Phi(t,u,n)$, also called as dynamical map, which is a completely positive and trace preserving superoperator. The evolution superoperator $\Phi(t,u,n)$ acts on the initial density matrix $\rho_{0}$ and maps it to the solution of the system (2) evolving under controls $u$ and $n$ for $t\geq 0$:

[TABLE]

To find a control that steers the initial state $\rho_{0}$ to the target state $\rho_{\textrm{target}}$ as close as possible, one can use the Hilbert–Schmidt distance

[TABLE]

and solve the optimal control problem of finding controls $u$ and $n$ which minimize the objective functional

[TABLE]

It should be noted that there might be many different controls that steer the system to the given target state.

Generation of a unitary quantum gate $U$ means that we are looking for controls $u$ and $n$ which produce the evolution operator $\Phi(T,u,n)$ that is close enough to the ”target” dynamic map $U\cdot U^{\dagger}$. Hence, we need to use a functional which will describe similarity of two dynamic maps. Let $\rho^{(j)}$, $j=1,\dotsc,N$ be the basis density matrices in the space of $N\times N$ density matrices. Two linear maps coincide if they act identically on the basis. Therefore, the following functional can be used to compare two dynamic maps $\mathcal{W}$ and $\mathcal{W\,^{\prime}}$:

[TABLE]

Hence, for gate generation, we consider the objective functional

[TABLE]

Now, we formulate this problem using the Bloch parameterization. Make the correspondence between density matrices and vectors in the Bloch ball as

[TABLE]

Then, the problem (9) is formulated as

[TABLE]

It can be seen that values of the objective functional (12) are bounded, both below and above, as

[TABLE]

Finding exact maximum and minimum values for a fixed time $T$, basis $\rho_{0}^{(j)}$ and gate $U$ is not a trivial problem.

We use $N=4$ initial pure states $\rho_{0}^{(j)}$, $j=1,2,3,4$, which form a basis in the linear space of $2\times 2$ Hermitian matrices, to define action of the unitary gate $U$:

[TABLE]

Here, $\rho_{0}^{(j)}$ are pure state density matrices formed by eigenvectors of Pauli matrices $\sigma_{z}$ (vectors $|0\rangle$ and $|1\rangle$), $\sigma_{x}$ (vector $|+\rangle=(|0\rangle+|1\rangle)/\sqrt{2}$), and $\sigma_{y}$ (vector $|i\rangle=(|0\rangle+i|1\rangle)/\sqrt{2}$). Note that, as was mentioned in the Introduction, three special mixed states are enough in the general case of a quantum system of an arbitrary dimension **[32]**.

Any single-qubit initial density matrix $\rho_{0}$ can be uniquely decomposed in the basis (13):

[TABLE]

Therefore, this number of states is sufficient to define any linear operator: particularly, any unitary transformation $\rho\mapsto U\rho U^{\dagger}$. The corresponding Bloch vectors $\mathbf{r}_{0}^{(j)}$ are

[TABLE]

We also consider the objective functional $F_{U,2}$ with only two states $\rho_{0}^{(0)}=|0\rangle\langle 0|$ and $\rho_{0}^{(1)}=|1\rangle\langle 1|$ to show that this system is insufficient for generating quantum gates. Indeed, the quantum gate $\sqrt{Y}\sim\exp(-i\pi/4Y)$ (rotation along y-axis by $\pi/2$ angle in the Bloch ball) acts exactly as the Hadamard gate $H$ on the states $\rho_{0}^{(0)}$ and $\rho_{0}^{(1)}$ :

[TABLE]

However, these two gates are different.

5 Gradient of the Objective Functional

The main features of quantum control landscape are local (if they exist) and global extrema of the objective functional. Since extrema are critical points where gradient of the objective is zero, an important step for the analysis of quantum control landscapes is the derivation of the expression for gradient of the objective. Below, we derive the gradient for piecewise constant control. Incoherent control is constrained: $n(t)\geq 0$; therefore, we make a variable change:

[TABLE]

so that control becomes a pair of $u$ and $w$. The controls, $u(t)$ and $w(t)$, as piecewise constant functions, are defined by $M$-dimensional real vectors $u_{k}$ and $w_{k}$:

[TABLE]

where $0<t_{0}<t_{1}<\dots<t_{M}=T,$ and $\chi_{[t_{k-1},t_{k})}$ is the characteristic function of the half-open interval $[t_{k-1},t_{k})$. For the sake of brevity, we denote this pair as $v=\left(v^{1},v^{2}\right)=(u,w)$ and its components as $v_{k}=\left(v^{1}_{k},v^{2}_{k}\right)=\left(u_{k},w_{k}\right)$, $k=1,\dotsc,M$.

The matrix $A(t)$ in the right-hand side (r.h.s.) of the Equation (4) also takes only $M$ different values $A_{k}$:

[TABLE]

The solution of the equation with piecewise constant r.h.s can be expressed using the following recurrent formula:

[TABLE]

where

[TABLE]

*Differentiating the final state $\mathbf{r}_{M}=\mathbf{r}(T)$ (21) with respect to the *$k$th component of piecewise constant controls $u$ (19) and $w$ (20) gives

[TABLE]

The derivatives are given by

[TABLE]

The expressions (25) and (27) were found using the integral formula **[63]**:

[TABLE]

Finally, the derivatives of the functional (12) with respect to the components $u_{k}$ and $w_{k}$ of the controls $u$ and $w$, respectively, are

[TABLE]

6 Numerical Optimization for Generation of $X$ and $H$ Gates

In this section, we use the above derived gradient (29) of the objective function to implement a gradient optimization method for generating single-qubit quantum gates. We consider two single-qubit gates: the $X$ gate, which is $\sigma_{x}$ Pauli matrix, and the Hadamard gate $H$, i.e., $U\in\{X,H\}$:

[TABLE]

To find controls $u$ (19) and $w$ (20) that minimize the objective functional (12), we utilize the gradient descent method known in the context of quantum control as GRAPE. This is an iterative algorithm, where each update of the control $v=(u,w)$ describes moving against the gradient of the objective functional:

[TABLE]

Here, for brevity, we denote

[TABLE]

*where the basis matrices $\left\{\rho^{(j)}\right\}_{j=1}^{4}$ are fixed and chosen as (13), their Bloch representations are given by (15); $\mathrm{grad}_{v}$ denotes a $2M$-dimensional vector operator, each of the components acts according to formula (23); $v^{(l)}=\left(u^{(l)},w^{(l)}\right)$ is the control on the *$l$*th iteration; and $h^{(l)}$ is the size of the *$l$th step.

*The stopping criterion is the condition for the fidelity to be small enough, such that fidelity value on the *$l$th step is less than some precision parameter $\varepsilon_{1}>0$:

[TABLE]

This stopping criterion is satisfied when we find an optimal control $v=(u,w)$ that steers the initial basis states $\rho_{0}^{(j)}$ to the final states $\Phi\left(t,u,w^{2}\right)\rho_{0}^{(j)}$ that are different from the action of the unitary gate $U$ on them, $U\rho_{0}^{(j)}U^{\dagger}$, less than $\varepsilon_{1}$ in the sense of mean square of the Hilbert–Schmidt distances (9). Since each control $v=(u,w)$ defines some dynamical map $\Phi\left(t,u,w^{2}\right)$ on the space of density matrices, we attempt to construct a dynamical map which is as close as possible to the desired map $U\cdot U^{\dagger}$ produced by the target quantum gate $U$.

*Performance of gradient descent (31) largely depends on the choice of a step size sequence $h^{(l)}$. Often with a predefined sequence $h^{(l)}$, gradient descent method may produce a fidelity value sequence $G_{U}^{(l)}$ which is not monotonically decreasing. We propose the adaptive choice of step size and modification of the control update rule (31). In each *$(l+1)$*th iteration, we compute the potential next updated value of fidelity G_{U}\left(v^{(l)}-h^{(l)}\mathrm{grad}_{v}G\big{(}v^{(l)}\big{)}\right) and compare it with the last value on the *$l$th iteration $G_{U}^{(l)}$; if it is lower, we update the control; otherwise, the control and the fidelity value remain the same:

[TABLE]

Increase of the fidelity means that the step size is too large; therefore, we decrease the next step in the second case. We also increase the step size in the first case to boost up the gradient descent:

[TABLE]

where $c\geq 1$, $0<d<1.$

Thus, we produce a monotonically (not strictly) decreasing sequence of fidelity values $G_{U}^{(l)}$. This scheme allows us to adapt step size to the local change of the functional landscape. This approach may significantly reduce the number of iterations needed to reach the desired accuracy of optimization. The adaptive scheme (35) depends on the three parameters, which were chosen as: $h^{(0)}=1$, $c=1.1$, $d=0.5$.

The important question is whether the adaptive scheme (34) and (35) may stuck in a non-critical point $x$ (i.e., in a point with gradient $\nabla f(x)\neq 0$), i.e., may constantly reduce size of the step without decrease of the objective functional. If gradient is calculated with absolute precision, there would exist a step size $h$ which is sufficient to decrease the value of the objective functional. This intuitive statement can be proved by the following simple reasoning. Consider a differentiable functional $F:D\subset\mathbb{R}^{n}\to\mathbb{R}$. Then, in a vicinity of an internal point $x\in D$:

[TABLE]

Due to the definition of $o(h)$, for any $\epsilon>0$, there exists $h>0$ such that $o(h)<\epsilon h$. If $x$ is not a critical point (so that $\nabla F(x)\neq 0$), we can take $\epsilon=|\nabla F(x)|^{2}$ and $F(x-h\nabla F(x))<F(x)$. Thus, the answer to the posed question is ”no” — after a finite number of iterations (35), the step size will be small enough to decrease the objective functional value. However, practice of numerical experiments shows that the adaptive scheme sometimes does not provide the objective functional value decrease over quite large number of iterations. In order to deal with such cases, we utilize one more stop criterion. The second criterion is the following: if $G_{U}^{(l)}=G_{U}^{(r)},r=l-1,l-2,\dotsc,l-L_{\textrm{stuck}},l\geq L_{\textrm{stuck}}$, the iterations also are stopped. We set $L_{\textrm{stuck}}=20$ as the maximum number of ”stationary” iterations.

Numerical computations were performed by writing a Python program using the Numpy library for fast matrix operations and the SciPy library function scipy.linalg.expm for matrix exponentials computation using the Padé approximation. The following values of the system parameters were used: the transition frequency $\omega=1$, the dipole moment $\mu=0.1$, and the decoherence rate coefficient $\gamma=0.01$. We use regular partition of the time segment $[0,T]$ with $T=5$ into $M=10$ segments, such that each segment has length $\Delta t_{k}=T/M=0.5$. Integral formulae were approximated by the trapezoidal rule with the number of partitions $N_{\rm partition}=20$.

Initial guesses for a coherent control $u^{(0)}$ and for an incoherent control $w^{(0)}$ were chosen as discretization of the sine function and the normal distribution function (similar to physical realization of coherent and incoherent controls): $u^{(0)}_{k}=\sin{(2\pi t_{k-1}/T)}$, $n^{(0)}_{k}={w^{(0)}_{k}}^{2}=\exp{\left(-4(t_{k-1}/T-1/2)^{2}\right)}$. The precision parameter was $\varepsilon=0.001$.

Figure 2 illustrates optimized coherent (subplot a) and incoherent controls (subplot b) for generating the gate $X$ (blue line) and the Hadamard gate $H$ (purple line), convergence of the fidelity — the objective functional $G_{U}^{(l)},\;U\in\{X,H\}$, (33) (subplot c)—and norm of the gradient of the objective functional (subplot d) vs number $l$ of the optimization algorithm (31) iterations under the optimized coherent and incoherent controls. As shown on subplot (a), for generation of $H$ gate, the fidelity goes to a small value of $10^{-3}$ in about $30$ iterations. For this gate, the convergence is fast. For generation of $X$ gate, the convergence to such a low value of the fidelity is much slower; even after $120$ iterations, the fidelity is still slightly above this value. Thus, high-precision generation of $X$ gate is more difficult with the present gradient approach. For Hadamard gate, which is found to be easily generated with this approach, higher (compared to the case of $X$ gate) values of incoherent control were obtained, as shown on subplot (b). The difference of scales is large, since our target in this analysis is a unitary gate, and as confirmed by the plots, it is natural to expect that incoherent control, while non-trivial, will have much smaller values than coherent. These gates are good testing examples and the results are in the agreement with expectation. The next step will be to apply the method for generation of target non-unitary processes, where incoherent control component should be more significant, as was found before for the problem of mixed state generation **[11]**. Coherent controls, which were found here as the results of optimization, for $X$ and $H$ gates are to some degree symmetric around a constant value $u\approx 2$.

Figure 3 shows optimal trajectories of the basis states (15) in the Bloch ball for generation of the quantum gate $X$ (subplot a) and the Hadamard quantum gate $H$ (subplot b), where four initial states, corresponding to $\rho_{0}^{(j)}$ ($j=1,2,3,4$) are indicated by green dots, final states by red squares, and target states are centers of the black circles. Optimal trajectories are very close to the Bloch sphere (the longest distance to the sphere is up to $0.05$). However, they are located in the interior of the Bloch ball due to unavoidable non-unitary evolution of the system and the presence of incoherent control.

We note that on Figure 3, the norm of the gradient has oscillating behavior vs number of iterations. This is in contrast to the case of **[50]** (see Figure 2 in that paper), where the norm of the gradient of the objective first monotonically increases, and then, monotonically decreases with the number of iterations. This difference in the behavior is expected to be due to the adaptive step size used in the present work.

7 Discussion: Quantum Gate Generation as Optimization Over Complex Stiefel Manifolds

Our analysis in this work is directly related to the gradient optimization approach which was developed for open quantum systems using a proposed formulation of completely positive trace preserving dynamics of open quantum systems as points of complex Stiefel manifold (rigorously speaking, as points in quotient space of complex Stiefel manifolds over some equivalence relation). The corresponding theory of gradient optimization over complex Stiefel manifolds was developed for two-level **[61]** and general $N$–level quantum systems **[62]**.

Recall that the dynamics of density matrix of an $N$-level open quantum system can be described as a completely positive trace preserving map (Kraus map), which has an operator sum representation

[TABLE]

where $K_{i}$ are $N\times N$ complex matrices such that

[TABLE]

and $\mathbb{I}_{N}$ is the identity $N\times N$ operator. Denote $N^{3}\times N$ matrix

[TABLE]

The condition (36) is equivalent to

[TABLE]

The set of matrices which satisfy this condition forms a complex Stiefel manifold $V_{N}(\mathbb{C}^{N^{3}})$, which is the set of all orthonormal $N$-frames in $\mathbb{C}^{N^{3}}$. As is well known, a given physical evolution (Kraus map) has a non-unique Kraus map operator-sum representation. Thus, different sets of $\{K_{i}\}$ may describe the same physical evolution and one should take the appropriate equivalence class ${V_{N}(\mathbb{C}^{N^{3}})}{/\!\!\sim}$. This equivalence was taken into account in **[62]**. In particular, the Hadamard gate corresponds to the $8\times 2$ matrix

[TABLE]

which is a point of the complex Stiefel manifold $V_{2}(\mathbb{C}^{8})$.

Coherent and incoherent controls which act on our one-qubit system (and more generally, on any quantum system) induce some time evolving Kraus map $\Phi_{t}$, which induces some trajectory on the Stiefel manifold. Thus, optimization of gate generation can be formulated as optimization over complex Stiefel manifold. The corresponding theory, including gradient and Hessian-based optimization, in the kinematic picture was developed in **[61, 62]**. In our example, optimal controls plotted on Figure 2 induce some trajectory on the Stiefel manifold. Since the dimension of the manifold is $2N^{4}-N^{2}=28$, the corresponding evolution can not be drawn as one can draw for evolution of single-qubit states on the Bloch ball, but the Stiefel matrices for any moment in time can be computed.

8 Conclusions

In this work, we have exploited the coherent control and the environment formed by incoherent photons as a resource for generation of single-qubit gates in a two-level open quantum system evolving according to the Gorini–Kossakowski–Sudarshan–Lindblad master equation with time-dependent coefficients. The coherent control enters in the Hamiltonian, and the incoherent control enters in the dissipator and induces time-dependent decoherence rates of the system. The exact expression for the gradient of the objective functional with respect to piecewise constant controls has been obtained, and the optimization has been performed using the gradient type algorithm. The optimal coherent and incoherent controls and the basis states trajectories in the Bloch ball have been computed. The numerical optimization shows that generation of the $H$ gate is fast (a value of $10^{-3}$ is obtained in about $30$ iterations), giving low values of fidelity. For the $X$ gate optimization, it takes a much higher number of iterations ($120$* iterations were not enough). Thus, high-precision generation of $X$ gate is more difficult with the present approach, that might be due to the limits of controllability of the system [10]. For the Hadamard gate, higher values of incoherent control were obtained, compared to the case of the $X$ gate. The coherent controls, obtained as the results of optimization for $X$ and $H$ gates, are found to some degree to be symmetric around a constant value $u_{*}\approx 2$. In this work, we use the adaptive step size of the gradient search to analyze the landscape. This leads to oscillating behavior of the norm of the gradient vs number of iterations; this contrasts Figure 2 in [50], where the norm of the gradient of the objective firstly monotonically increases, and then, monotonically decreases with the number of iterations. While we use the particular model for the dissipative superoperator, the approach can be extended to other cases as well. Here, we develop the theoretical basis of the approach and apply it for target unitary gates as a testing example. In this case, the incoherent control magnitude is obviously small. The next step is to apply the method for generation of target non-unitary processes, where incoherent control component should be more significant, as was found before for the problem of mixed state generation [11].*

Acknowledgements

The new results in Sections 3–6 were obtained in Steklov Mathematical institute and supported by Russian Science Foundation under grant No. 22-11-00330, https://rscf.ru/en/project/22-11-00330/, and in Section 7 in the University of Science and Technology MISIS under the federal academic leadership program “Priority 2030” (MISIS Grant No. K2-2022-025)