= Solution
The terminal term $\mathcal A$ is the objective reward. With the <Hilbert-Schmidt inner product> $\langle\!\langle A|\rho\rangle\!\rangle=\operatorname{Tr}(A^\dagger\rho)$, a Hermitian target observable gives its final <expected value>. The vector $\rho_v$ is the <density operator> represented as a vector in operator space; it retains the same physical content as the <density matrix>, not a wavefunction of the system.
The term $\mathcal D$ imposes the dynamical equation by a time-dependent <costate> $A_v$. It is a Lagrange-multiplier operator in the same operator space, propagated backward from a terminal condition. It is not required to be a positive, trace-one <density matrix>. Define
$$
K[f]= -\frac i\hbar\mathcal L_{\rm tot}[f].
$$
Then the enforced state equation is $\dot\rho_v=K[f]\rho_v$. The uncontrolled coherent part is $\mathcal L_0$, the dissipative contribution is $\mathcal L_D$, and each $f_m\mathcal L_m$ describes a control coupling. For <Quantum Hamiltonian> control, $\mathcal L_m\rho=[H_m,\rho]$.
The term $\mathcal C$ is a quadratic field-energy or fluence penalty. For the positive scale $p_0$ and positive weights $\lambda_m$, larger $\lambda_m$ make a given control amplitude more costly. Thus \b[maximizing $J$ trades terminal performance against control fluence while enforcing the dynamics]. The multiplier term vanishes on any dynamically feasible trajectory. The weights do not represent <decoherence> rates; they specify control-resource preferences.
Back to article page
