Quantum Sensing of Noisy and Complex Systems under Dynamical Control

We review our unified optimized approach to the dynamical control of quantum-probe interactions with noisy and complex systems viewed as thermal baths. We show that this control, in conjunction with tools of quantum estimation theory, may be used for inferring the spectral and spatial characteristics of such baths with high precision. This approach constitutes a new avenue in quantum sensing, dubbed quantum noise spectroscopy.


Introduction
Nanoscale sensing may be defined as the gathering of information concerning the structure, properties or dynamics of atomic, molecular, or condensed-matter systems by probes capable of nanometric resolution.Such probes are amenable to quantum-mechanical description even though their quantumness, i.e., quantum coherence or entanglement, may or may not play a role in the sensing process.If the goal is to acquire quantum information (QI), i.e., read out the state of the probed system, then measurements of quantum correlations within or among probes become indispensable [1].We, however, focus here on the much more robust and readily obtainable classical information on the system by simple quantum probes: either single or correlated qubits that are engineered to possess advantageous properties for the sensing process in question: Thus, superconducting qubits are excellent detectors of microwave or far-infrared photons [2]; Rydberg-atom qubits are highly sensitive probes of dipolar forces; and nitrogen-vacancy centers (NVC) in diamond are the most sensitive magnetometers or electrometers to date [3][4][5][6][7].We shall not be concerned with the probing of mechanical degrees of freedom (phonon modes) for which the most appropriate probes are harmonic oscillators, as in optomechanical setups [8].The peculiarity of the approach we present here is that rather than try to maximize the signal-to-noise ratio (SNR) obtainable by the probe, we set out to analyze and extract information from the noise exerted by the system on the probe [9].Namely, the spectrum of fluctuations may be a viable source of information concerning a noisy, complex system which we shall treat as a "bath", e.g., a large molecule, an aggregate of interacting spins, a liquid, or stray fields at diverse frequencies.This emerging area of research, nicknamed "noise spectroscopy" [9][10][11][12][13], requires a detailed consideration of the probing process, in which we schematically discern the following stages: change of the probe state under the influence of the noise and control; transfer and storage of this state or at least its partial information content to either a storage device or an information processor for processing or improving the sensing process; readout and analysis of the extracted information.

Bath-Optimized Task-Oriented Control (BOTOC)
Our general framework, dubbed BOTOC, allows for controlling the interaction of a quantum probe with its (usually multiple) baths, and is aimed at optimally executing a chosen operation, whether unitary or non-unitary.Examples are: state-transfer or storage with the highest fidelity possible, state purification, entanglement distribution, or energy transfer.To this end, we design the temporal control of the Hamiltonian that governs the system, subject to variational minimization or maximization (as the case may be) of a functional that quantifies the success probability of the operation.The control must be faster than the bath correlation (memory) time in order to counter or change bath-induced decoherence.For certain tasks, our control (BOTOC) may also benefit from the probe-bath interactions when the operation is non-unitary, e.g., state purification [34,35].BOTOC is not limited to pulsed control, as opposed to DD [36][37][38][39][40][41][42], and thus may drastically reduce the energy and control errors if a smoothly varying field is employed instead of pulses [34,35].These considerations remain true for our present purpose: the use of BOTOC as a tool of bath (quantum noise) diagnostics by quantum probes.
The bath response function is usually associated with a characteristic correlation or memory, τ, which separates the non-Markovian (t τ) and the Markovian (t τ) regimes, as in the case of an exponentially decaying bath-memory functional Φ(t) ∝ e − t τ .Within the bath memory time, the bath modes oscillate in unison, and preserve the "memory" of their prior interaction with the system, whereas after the correlation time has passed, the modes lose their mutual coherence and "forget" their previous interactions [50][51][52][53].

General Formulation of BOTOC
To lowest (i.e., second) order in the probe-bath interaction, ĤI , the probability change of an Hermitian operator P over time t can be shown to obey where denotes averaging with respect to the initial state of the probe in the interaction picture [34,35].
We assume a d dimensional Hilbert-space of the probe, and expand the probe-bath interaction Hamiltonian of a controlled (driven or modulated probe) as a sum of products of probe and bath operators, In the interaction picture we may expand in terms of (Hermitian, orthonormal) basis operators Ŝj .This expansion defines a (real orthogonal) rotation matrix ε(t) with elements in the probe's space These matrix elements are the correlation functions of the basis operators in the presence of a (driving) control.
The finite-time Fourier transform of this rotation matrix is denoted by the spectral functional We next define a bath autocorrelation matrix Φ(t) with elements and its (finite-time) Fourier transform, known as the bath spectral response (fluctuation spectrum) The inverse of the bandwith over which G(ω) varies considerably, defines the bath correlation (memory) time τ.
Finally, we define a Hermitian matrix Γ with elements Then we may cast the probability change of the operator in Equation (1) in the form of an overlap of two matrices In Equation (9) we have used the rotation matrix spectrum ε t (ω) and the matrix Γ to define the probe spectral matrix, dubbed the "filter" functional determined by the control field-intensity spectrum and its time duration t.Equation (10) demonstrates that the change P over a given time t is determined by the spectral overlap between the controlled-probe and the bath dynamics [34,35,[46][47][48][49].The bath-spectral matrix G(ω) must be positive semi-definite for all ω.
In general, F t (ω) is Hermitian but not necessarily positive semi-definite.This reflects the fact that the probability P can increase or decrease over t, according to the choice of Γ, depending on whether our control is meant to maximize or minimize P [34,35].
This formulation may be used to optimize the success probability P of any operation performed on a quantum probe whose weak interaction with thermal baths is dynamically controlled [34,35].It applies to any such control, be it coherent or projective (non-unitary), continuous or pulsed.Moreover, this optimization may be done under constraints on the total control energy E introduced by controlling the probe.
If P is the decoherence probability that should be minimized, this implies we must resort to the QZE-based control, namely, minimization of the integral product (overlap) of G(ω) and F t (ω).The filter functional is the control handle on our ability to optimally execute the task in the presence of a given bath or several baths [34,[46][47][48][49].
QZE-based control is required in operations related to QI storage and transmission [43,68], where bath effects are detrimental and must be suppressed.As shown here, it is likewise instrumental for extracting information on the bath characteristics [9,15].Regardless of the chosen form of control, the controlled-probe dynamics must then be Zeno-like, namely, result in suppressed probe-bath interaction.
Alternatively, we may implement anti-Zeno (AZE) control, namely, maximized overlap of G(ω) and F t (ω) (again, under control-energy constraints).AZE-control is instrumental for non-unitary operations that change the probe's entropy and benefit from efficient interaction with a bath for their execution [31,[51][52][53]70].If the underlying dynamics is anti-Zeno-like, probe-bath interaction will be enhanced and thereby facilitate these tasks.
Entanglement via the bath may require maximized bipartite coupling, but minimized single-partite coupling with the bath [72][73][74].This requires an interplay of QZE and AZE that may depend on the quantum statistics of the bath [65].

Dynamical Control of Qubit Dephasing (Decoherence)
A simple application of the foregoing framework is the control of qubit probe decoherence, or proper dephasing, that causes random frequency fluctuation characterized by the correlation time, τ.The dephasing rate, if no control is applied, is: where Φ(t) is dephasing autocorrelation function.
At times much longer than the correlation time, the decoherence rate approaches its asymptotic value R = lim t→∞ R(t) that satisfies R τ −1 .We assume that the control (driving) field has a Rabi frequency V(t) with amplitude This is in fact a specific example of Equations ( 9) and (10), that yields a dynamically modified decoherence rate according to the universal Kofman-Kurizki (KK) formula [44][45][46][47][48]50,55,56,60] where the filter function F t (ω) is given by Equations ( 5) and (11) with and It follows from Equation ( 16) that G(ω) is a symmetric function, The proper dephasing rate associated with an exponentially decaying correlation function is given by R = Aτ.
In the presence of a constant V 0 (continuous wave (CW) Hence, provided V(t) is sufficiently strong, the dephasing rate R can be reduced by the factor (V 0 τ) 2  1.This means that strong Rabi splitting can shift the qubit probe away from the noise bandwidth, or average out the noise effects.Periodic π-pulses may similarly suppress dephasing, if their period satisfies V 0 = π τ .

Modulation Forms
Under a general modulation with frequencies the long-time decoherence rate is (Figure 1d): Let us apply to the qubit probe an impulsive periodic phase-modulation, ε(t) = e i[t/τ]ϕ , consisting of phase jumps ϕ with period τ.This modulation can be a periodic AC-Stark shift to control decay or a periodic resonant field to control decoherence.The power spectrum (filter function) of this control consists of equidistant Fourier (harmonic) peaks with the following frequencies and weights: For small phase shifts, ϕ 1, the main contribution is that of the k = 0 harmonic, the higher harmonic weights scaling as (Figure 1b) As |ϕ| increases, the difference between the k = 0 and k = 1 weights diminishes until the phase jumps by ϕ = π Then, the filter function F t (ω) has two main symmetrically shifted peaks (Figure 1c) while higher harmonics |λ k | 2 decrease with k as (2k − 1) −2 .This modulation then coincides with dynamical decoupling (DD) pulses [36][37][38][39][40][41][42].

Optimal Decoherence Control of a Qubit
It is appropriate to compare the effect of dynamical decoupling (DD) [36][37][38][39][40][41][42] with that of BOTOC.DD does not change with the bath spectrum ( ).With an increasing pulse rate, DD shifts the weight of the system spectrum ( ) towards higher frequencies, until the overlap Equation ( 14) becomes sufficiently small.This DD procedure may, however, fail if the bath spectrum has a multipeak structure because then the overlap may increase rather than decrease with pulse rate whenever the vicinity of a peak is probed.By contrast, BOTOC is optimally tailored to any bath spectrum once it has been properly diagnosed by the probe.
Hence, DD sequences [36][37][38][39][40][41][42] may not be optimal for decoherence suppression, when they are insensitive to the bath spectral response (coupling spectrum).Instead, here we resort to the BOTOC approach outlined above [34,35], subject to minimal-energy constraints, so as to deduce the optimal, energy-constrained control of a decohering qubit probe by modulation, for any given bath-coupling spectrum.

Optimal Decoherence Control of a Qubit
It is appropriate to compare the effect of dynamical decoupling (DD) [36][37][38][39][40][41][42] with that of BOTOC.DD does not change with the bath spectrum G(ω).With an increasing pulse rate, DD shifts the weight of the system spectrum F t (ω) towards higher frequencies, until the overlap Equation ( 14) becomes sufficiently small.This DD procedure may, however, fail if the bath spectrum has a multipeak structure because then the overlap may increase rather than decrease with pulse rate whenever the vicinity of a peak is probed.By contrast, BOTOC is optimally tailored to any bath spectrum once it has been properly diagnosed by the probe.
Hence, DD sequences [36][37][38][39][40][41][42] may not be optimal for decoherence suppression, when they are insensitive to the bath spectral response (coupling spectrum).Instead, here we resort to the BOTOC approach outlined above [34,35], subject to minimal-energy constraints, so as to deduce the optimal, energy-constrained control of a decohering qubit probe by modulation, for any given bath-coupling spectrum.

Control for Bath Diagnostics
Our main purpose here is to dynamically control qubit-bath coupling not as a means of fighting decoherence but as a tool of bath-spectrum (noise) diagnostics [9][10][11], and thereby study the dynamics of decoherence.This diagnostics is again based on Equation ( 14), the KK universal formula [44][45][46]50,55,56,[60][61][62].Even if the ultimate goal of any modulation scheme is to reduce, or, if possible, eliminate the decoherence, one must first know the probe-bath coupling spectra so as to obtain the optimal modulation using BOTOC [34,35].This information is normally not available a-priori and thus one must use known, prescribed controls for revealing the unknown bath spectra that characterize the probe-bath coupling interaction, so then one can tailor optimal modulations for determined bath parameters.
The protocol consists in changing the filter function F t (ω), e.g., by varying the control-field Rabi frequency, measuring the corresponding decoherence rate R(t) and inferring G(ω) from the KK formula.The qubit probe is prepared at t = 0 in a superposition of its excited (| ↑ ) and ground (|↓ +| ↑ ), corresponding to a mean coherence that decays as R(t) being the dynamically controlled decoherence rate.
The probability of the qubit to occupy the symmetric or antisymmetric superposition of the upper and lower states is then given by One may infer R(t) from these probabilities that are obtainable from measurements.We have experimentally demonstrated the reconstruction of the bath-coupling spectrum G(ω) by measurements performed on a trapped cold-atom ensemble [10], and on nuclear spin systems [11].By monitoring the decoherence rate R(t) using control fields that generate narrow band filters F t (ω), like the ones described in Equation ( 24) and Figure 1a, we could scan the overlap integral in the KK formula of Equation ( 14) to infer the bath-coupling spectrum in the weak-coupling limit.

Maximized Information on Bath's Parameters by Dynamical Control
In order to further advance the noise spectroscopy approach and broaden its applicability, it is imperative to find the best general strategy for extracting environmental information from the qubit probe decoherence.Extracting the full noise spectrum is a very demanding process, and usually contains more information than required for diagnostic purposes.The strategy we adopt to improve the noise spectroscopy approach is to estimate key parameters that characterize the bath spectrum and make a faster and precise diagnostic.We aim at minimizing the error in estimating these bath key parameters x B by N m measurements performed on the qubit probe (Figure 1a), which obeys the bound derived by resorting to quantum estimation tools [75].
Here, the dependence of the controlled decoherence rate R(t) ≡ R(x B , t) on G(ω) ≡ G(x B , ω) is used in the expression for the quantum Fisher Information (QFI) [9,13,75] The relative error ξ(x B , t), Equation (31), may be minimized either per measurement N m or per total interrogation time T = N m t [9], both alternatives being of practical importance in experimental situations.Then, for the minimal relative error per measurement, the QFI has to be maximized: which amounts to measuring the qubit state at the optimal time and optimally controlling the qubit probe.
For the minimal ξ(x B , t) per total interrogation time, the QFI per unit of time needs to be maximized: Typical environments of interest for sensing are those undergoing a diffusion process characterized by a generalized Ornstein-Uhlenbeck spectral density with α ≥ 2 and the renormalization factor A α .We have shown that there is an ultimate accuracy bound for the estimation per measurement of the key parameters x B , such as the probe-bath coupling strength g and the correlation time τ, that may be achieved only under adequate control of the probe [9].This knowledge allows for the design of the appropriate dynamical control strategy for each key parameter.For example, the maximal accuracy on g is attained by stroboscopic measuring of the probe at the QZE regime [9], while τ is attained by a CPMG sequence control on the probe which generates a narrow spectral filter [9,13].
On other hand, if uncontrollable noise sources are also present, which are not of interest for the sensing process such as an intrinsic T 2 decoherence due to a white noise (Markovian) process or due to pulse imperfections, we have resorted to more elaborate controls to selectively address the targeted noise source and key parameters.We therefore nicknamed the method "selective dynamical recoupling" (SDR) [14,16].We have already implemented SDR to selectively monitor diffusion processes characterized by Equation ( 35) with α = 2 and A α = π −1 , so as to determine the probe-environment coupling strengths [18] and the environmental correlation time, as well as diffusion restriction lengths [16,17,19] in systems of biological and chemical interest.

Dynamical Control of Quantum Transfer as Means of Diagnosing a Bath in Hybrid Systems
Commonly, manipulations of QI can be schematically divided into three stages: "writing-in," "storage," and "reading-out" [20][21][22]76].Realistically, some systems are better suited for writing-in or reading-out than for storage, and vice versa.This has prompted the suggestion of hybrid, composite quantum systems [2]: quantum operations are rapidly performed and efficiently written in a qubit susceptible to decoherence, e.g., a Josephson (superconducting) qubit; then the QI is transferred (directly or via a link) to a storage qubit resilient to decoherence (encoded in an ensemble of, e.g., ultracold atoms); then, on demand, transferred back and read-out from the same fragile qubit.
Our aim is to examine a strategy for accomplishing the transfer with the highest fidelity and (or) speed to store or extend the life of the sensing information for improving the resolution of the sensing process.We also discuss relevant effects for extracting information in the bath characteristics from this transfer.This quantum (state)-transfer in hybrids will be affected from a noisy subspace that is fragile against decoherence to a quiet, robust subspace, by choosing an appropriate dynamical control [43] that optimizes this QI transfer.

Optimized Transfer from Noisy to Quiet Qubits
The transfer of a quantum state from a fragile (noisy) qubit probe to a robust (quiet) qubit may be used as a means of diagnosis, i.e., extracting information on the bath.Here we consider two resonant qubits at ω 0 with temporally controlled interaction characterized by V(t), the interaction strength that is adjustable by an external laser field [68], and descripted by the Hamiltonian Ĥc = V(t) σnoisy x ⊗ σquiet x where σi , i = x, y, z are the Pauli matrices.We choose the system-bath interaction Hamiltonian to be where B(t) is the bath operator B rotating with the free-bath ĤB .This model represents proper dephasing in the probe noisy-qubit due to the bath operator B, whereas the target quiet qubit is robust against decoherence.The accumulated phase is our control parameter.In the absence of decoherence, the state transfer from qubit 1 to qubit 2 can be completed if at the final time, T, the phase φ(T) satisfies φ(T) = π 2 .Then, any state of the (probe) noisy qubit is mapped onto that of quiet qubit, initially in the ground state, for any normalized γ, β.By contrast, under decoherence, since the source noisy qubit is affected by the bath, the longer the information stays there, the lower the fidelity of the transfer.On the other hand, if the transfer is very fast, it may result in population leakage from | ↓ noisy (| ↓ quiet ) into | ↑ noisy (| ↑ quiet ).Such leakage [77] signifies the violation of the rotating wave approximation (RWA) when the transfer rate is comparable to the level distance of the qubits, ω 0 .
We dwell here on bath effects in the RWA so that there is no leakage due to RWA breakdown, i.e., the transfer time being much slower than ω 0 −1 .
The general expression derived from Equations ( 9) and ( 10) for the averaged transfer fidelity at T is then: The energy constraint E, defined in Equation (28), determines the minimum possible transfer time Here we illustrate the transfer fidelity for a Lorentzian bath spectrum.Surprisingly, when the memory time τ of the bath is comparable to or larger than the characteristic transfer time τ T min , the best solution actually performs a full transfer, φ(t) = π 2 , well within the control time, but rather than stop at φ = π 2 , it overshoots the transfer so that φ(t) > π 2 , and then returns slowly to π 2 .During this interval, the information partially returns from the target (storage) qubit to the source (noisy) qubit probe.Then, similarly to the echo method, the noise acts to correct itself so that the non-Markovian bath effect is undone.This requires transfer times significantly larger than the minimal transfer time T min , up to 10T min or more, in order to achieve fidelity increases (up to 50% in Figure 2).than stop at = , it overshoots the transfer so that ( ) > , and then returns slowly to .During this interval, the information partially returns from the target (storage) qubit to the source (noisy) qubit probe.Then, similarly to the echo method, the noise acts to correct itself so that the non-Markovian bath effect is undone.This requires transfer times significantly larger than the minimal transfer time , up to 10 or more, in order to achieve fidelity increases (up to 50% in Figure 2).This method may be beneficial for the optimization of QI among different modules [2,66,67]: from superconducting qubit probes via a microwave resonator to ensembles of ultra-cold atoms or from NV-centers probe to nuclear spins [78].
To conclude, the state-transfer optimization within hybrid open systems, from a "noisy"qubit probe to its "quiet" counterpart, reveals the profound dependence of the transfer fidelity on the bath-memory time.Intriguingly, the fastest transfer is not optimal (for a given transfer energy).Hence, transfer fidelity maximization can disclose information on the bath non-Markovianity.Alternatively, one may use this method to significantly improve the storage and retrieval of QI in spectrally inhomogeneous spin ensembles [69] by preselecting an optimal portion of the ensemble via suitable control.

Dynamical Control of State Transfer via Noisy Quantum Channels and the Implications for Sensing Diagnostics
A similar approach may be used to exploit noisy quantum channels for state-transfer between the boundary (edge) qubits that are subject to temporal modulation and used for sensing.Alternatively, they can be exploited for diagnosis of the quantum channel.This modulation controls the coupling of the boundary system to the channels viewed as a fermionic bath.By means of this modulation one may design a spectral filter that suppresses the populations of noisy channel eigenmodes that impede successful transfer [43].We show that under optimal modulation, the fidelity and the speed of transfer can be improved by several orders of magnitude at the fastest possible transfer for a given fidelity.This filtering effect may be used to infer the noisy channel characteristics from the transfer fidelity and speed.

Boundary Qubit Probes and State Transfer through a Quantum Channel
The Hamiltonian of a chain of N + 2 spin-1/2 parcticles with XX nearest-neighbor interactions has the form Ĥ = Ĥch + Ĥqp (40) Here Ĥch and Ĥqp are the quantum channel and boundary-coupling probe Hamiltonians, respectively, and V i are the spin-exchange couplings.To simplify the treatment, we diagonalize Ĥch into the form Ĥch = , where b † k populates a fermionic eigenstate |ω k of energy ω k .We may further simplify the model upon noticing that there is a single, zero-energy fermionic mode in the channel, k = z = N+1 2 [43] if we assume mirror (inversion) symmetry of the couplings V i = V N−i for odd N. The boundary-qubit probes are resonantly coupled to ω z with an effective, controlled interaction strength V(t).This resonant interaction is described by an effective Hamiltonian of the form whereby the source and target (probe) mode operators are resonantly coupled to the zero-energy (k = 0) mode operator.These three fermionic modes constitute a "system" S that interacts with a "bath" B comprised of all other channel modes k (except for k = 0).On the basis of collective-mode operators associated with sums over the odd (even) bath-modes, the S-B interaction has the form where V k is the Fourier coefficient of the coupling to the kth mode.The dynamical control of this multipartite system then becomes similar to single-qubit level modulation [56,60-62] ( Figure 3).
Technologies 2016, 5, 1 12 of 22 whereby the source and target (probe) mode operators are resonantly coupled to the zero-energy (k = 0) mode operator.These three fermionic modes constitute a "system" S that interacts with a "bath" B comprised of all other channel modes k (except for k = 0).On the basis of collective-mode operators associated with sums over the odd (even) bath-modes, the S-B interaction has the form where is the Fourier coefficient of the coupling to the kth mode.The dynamical control of this multipartite system then becomes similar to single-qubit level modulation [56,[60][61][62] ( Figure 3).In frequency space, the probes are resonantly coupled by a fermionic-mode of the channel whose energy is defined to be the energy origin (zero).These 3 degrees of freedom constitute the effective system (green rectangular bars); while the remaining channel modes are taken to be the bath (red lines) with a gapped spectrum ± ( ).Under noise, the perturbed energies change/fluctuate bounded by the (colored gapped) Wigner-semicircle.In the gap, the optimal spectral-filters for the high fidelity statetransfer , ( ) correspond to dynamical boundary-control with ( ): = 0 (black dotted), = 2 (green); (b) Transfer infidelity as a function of the transfer time T for the gapped Wigner-semicircle spectral bath: static control = 0 (black dotted) and dynamic control = 2 (green); (c) Averaged transfer infidelity 1 − ( ) as a function of the perturbation strength of the noisy channel.Dynamical = 2 control (empty squares) strongly reduced the infidelity obtained under static control = 0 (empty circles).A fluctuating noisy channel results in less damage than the static-noisy channel → ∞, where its fidelity can even approach its unperturbed value in the Markovian limit → 0 ( = 0, green solid circles) .Here N = 29 channel-qubits and ≡ 1.
To evaluate the transfer of any state of the source to the target qubit probes, we resort to a second-order solution of our non-Markovian master equation for the system density matrix [43] to obtain the state-transfer fidelity averaged over all possible input states.Similarly to Section 5, complete state transfer through the mediated mode, in the absence of a bath, occurs when the accumulated phase over time , Equation (37), due to the modulation control ( ), satisfies ( ) = √ .In frequency space, the probes are resonantly coupled by a fermionic-mode of the channel whose energy is defined to be the energy origin (zero).These 3 degrees of freedom constitute the effective system (green rectangular bars); while the remaining channel modes are taken to be the bath (red lines) with a gapped spectrum G ± (ω).Under noise, the perturbed energies change/fluctuate bounded by the (colored gapped) Wigner-semicircle.In the gap, the optimal spectral-filters for the high fidelity state-transfer F T,− (ω) correspond to dynamical boundary-control with V(t): p = 0 (black dotted), p = 2 (green); (b) Transfer infidelity as a function of the transfer time T for the gapped Wigner-semicircle spectral bath: static control p = 0 (black dotted) and dynamic control p = 2 (green); (c) Averaged transfer infidelity 1 − F(T) as a function of the perturbation strength ξ V of the noisy channel.Dynamical p = 2 control (empty squares) strongly reduced the infidelity obtained under static control p = 0 (empty circles).A fluctuating noisy channel results in less damage than the static-noisy channel τ → ∞ , where its fidelity can even approach its unperturbed value in the Markovian limit τ → 0 (p = 0, green solid circles) .Here N = 29 channel-qubits and V i ≡ 1.

Optimized Filter Design for Channel Diagnostics
To evaluate the transfer of any state of the source to the target qubit probes, we resort to a second-order solution of our non-Markovian master equation for the system density matrix [43] to obtain the state-transfer fidelity averaged over all possible input states.Similarly to Section 5, complete state transfer through the mediated mode, in the absence of a bath, occurs when the accumulated phase over time T, Equation (37), due to the modulation control V(t), satisfies φ(T) = π √ 2 .

Optimized Filter Design for Channel Diagnostics
Modulation can control the transfer infidelity through the convoluted overlap of the bath and filter spectra governed by Here the ± indices are associated with the odd (even) parity modes that define the bath correlation functions φ ± (t) = ∑ k∈odd(even) V 2 k e −iω k t of Equations ( 6) and ( 7) and the rotations due to the modulation 4), ( 5) and (11).The optimal modulation control minimizes the overlap integrals of G ± (ω) and F T, ± (ω) for a given transfer time T by the variational Euler-Lagrange (E-L) method under energy constraint E, Equation ( 29).
Since we assume the channel to be symmetric with respect to the source and target qubits and the number of eigenvalues to be odd, the central mode results in being invariant to noise.Then, by considering this central mode as part of the system, the bath spectrum G ± (ω) has a central gap (Figure 3a).
Therefore, to minimize the transfer infidelity, we need a narrow bandpass filter centered on the gap that minimizes the overlap between G ± (ω) and F T,± (ω).Since G − (ω) has a narrower gap than G + (ω) , optimal analytical solutions for the control of the boundary probes can be found by the BOTOC method to be [79]: Here p = 0 signifies static control, while p = 2 stands for dynamical control during the time T and V 0 is T-dependent to satisfy φ(T) = π √ 2 .
For p = 0, the minimal-energy constraint E min = π 2 4T min , yields the narrowest filter around the bath-gap, but has many wiggles on the tails (Figure 3a) which overlap with bath-modes damaging the transfer.By contrast, for the same transfer time T, the p = 2 bandpass is wider and requires higher control energy, E 2 = 3 2 E min , but this filter is much lower (flat, without wiggles) at bath-frequencies.Hence, the bandpass filter width and its overlap with bath-energies as a function of T can optimize the modulations V(t).When there is a short T, p = 0 yields the highest fidelity, because the central peak of the filter is the narrowest.With a longer T, p = 2 will give rise to higher fidelity, because the filter tails crucially affect the overlap (Figure 3a,b), and the p = 2 filter can reduce the transfer fidelity by orders of magnitude in a noisy gapped bath.
The best tradeoff is obtainable for p = 2, as it yields a filter function without spectral tails.By chopping off the tails, one may reduce by several orders of magnitude the transfer infidelity (error) as well as the transfer time (Figure 3b).Thus, transfer guided by the optimized filter function can disclose the channel spectrum, particularly the effective gap width.

Diagnosing Noise Effects in Quantum Channels
Noise may affect the coupling strengths, causing: > 0 being the noise strength.We consider static (fluctuating) noise when ∆ i is time-independent (dependent).These noises affect the bath-energy modes, while the central (zero) energy ω z remains invariant, yet they are amenable to control: (i) Dynamical boundary-control can make the channel most robust against static noise, because it filters out the bath-energies that damage the transfer.In Figure 3c p = 2 control is shown to be advantageous compared with the static control case p = 0, at the expense of increasing the transfer time by only a factor of 2. (ii) Modulation is not helpful against Markovian noise.Remarkably, arbitrarily high fidelities can then be achieved by slowing down the transfer time, i.e., by decreasing V 0 , because in a Markovian bath, the fast coupling fluctuations suppress disorder-localization effects that hamper the transfer fidelity.(iii) For non-Markovian fluctuating noise V i + V i ∆ i (t), that randomly varies with correlation time τ (Figure 3b), in contrast to Markovian noise with τ → 0 , optimized dynamical control can strongly reduce the infidelity that lies between the static and Markovian limits, provided the bath-spectrum is gapped.
Hence, BOTOC dynamical control of the tradeoff between transfer speed and fidelity of qubit-state transfer through a quantum channel can disclose the characteristics of noise on qubit-qubit couplings.The fidelity and/or the transfer speed can be enhanced by orders of magnitude by appropriate control, while their robustness against noise on the qubit-qubit couplings remains intact.Remarkably, static noise is more detrimental but also more controllable than fluctuating noise.Hence, control allows for a clear spectral signature of the noise.
The detailed spectral distribution of the S−B coupling G(ω), is not required for BOTOC-based control design, only the spectral width and approximate mode spacing, which can be diagnosed by the methods of Section 4. Such channel noise diagnostics should suffice for designing the optimal tradeoff of the fidelity versus time transfer T by appropriate temporal modulation of the coupling.

Dynamical Control of Multipartite Probes for Bath Diagnostics
Decoherence control of a multiqubit probe coupled to a thermal bath, is much more challenging than that of a single-qubit probe, for a few reasons: (i) inter-qubit entanglement is more rapidly destroyed by the baths than single-qubit coherence [79]; (ii) cross-decoherence, namely mutual coupling of qubits via the baths, may complicate the control.We have nevertheless extended our decoherence control approach to multipartite scenarios [47,48,80].
The extension of our control formulation to decay and dephasing due to interaction with finite-temperature baths discloses additional information on the bath characteristics when one monitors the multipartite disentanglement in the time and spectral domains.
Namely, if one deals with a multipartite probe, one must apply multiple modulations that address the different qubits, on top of performing single-qubit modulations.This is essential in order to control the cross-coupling spectra of all possible qubit pairs [47,48], which are extremely important for reducing disentanglement and decoherence.
Our goal here is to exploit different modulations that effectively address different bath modes as a means of inferring the coupling strengths of these bath modes to different qubits comprising the probe.Subsequently, one may gather a wealth of information concerning the spatial and spectral structure of the bath.

Multipartite Decoherence Matrix
The decoherence rate (either due to decay or dephasing) of multiqubit probes weakly interacting with a thermal bath is an example of Equations ( 9) and (10), and similarly to Equation ( 14), it can be expressed in the spectral domain as a matrix with elements Here the diagonal elements of the decoherence matrix, R ii , are the single-qubit probe decoherence rates.At long times and without control, the single-qubit probe modulation power spectrum, F T,ii (ω) (see Equation (11), with Γ = I, where I is the identity operator), becomes a δ-function centered at the qubit energy separation, ω a , resulting in a coupling between the qubit probe and a particular mode of the bath, with coupling strength G(ω a ), in accordance with Fermi's Golden Rule.
The off-diagonal element, R ij , is comprised of the cross-relaxation rates of qubit i and j connected through the interaction with common bath modes, involving the overlap of the cross-coupling spectrum and the individual modulation spectra of the two qubits.Without modulations, the qubits comprising the probe that are coupled to different modes will not experience any cross-decoherence.Then, to induce cross-decoherence, the qubits should be properly controlled so as to couple them to the common modes with the same strength.Conversely, in order to eliminate cross-decoherence, one should locally modulate the qubits, at different rates, thereby coupling the two qubits to different modes.
What information can be inferred concerning the bath from such modulations?This issue will be discussed below.

Multiqubit Probe Modulations for Reconstructing Coupling Spectra
The multiqubit coupling spectra matrix, G ij (ω), which is usually unknown and is the objective of the sensing, can be reconstructed through the knowledge of: (i) the decoherence matrix, R ij (ω), which can be experimentally obtained via measurements of the decohering multiqubit probe; and (ii) the multiqubit modulation power-spectra matrix, F T,ij (ω), which can be calculated for each specific known modulation.
Although our approach applies to any number of excitations, simpler solutions are obtainable for dynamically modified relaxation of singly excited multiqubit states.
Both the coupling spectrum matrix, G ij (ω), and the dynamical modulation matrix, F T,ij (ω) are then complex and not real, as in the single-qubit scenario.This occurs because, as opposed to the single-qubit probe scenario, the off-diagonal terms of the multiqubit modulation matrix are determined by the phases between modulations of the different qubits comprising the probe.The diagonal elements of the decoherence matrix are easily separable into real decay rates and imaginary Lamb shifts.
We first consider a probe comprised of N identical qubits, i.e., set ω i ≡ ω a , subject to impulsive-phase modulations (see below Equation ( 23)), requiring AC Stark shifts multiple of 2π at any chosen time t = T.Then, modulations only affect the decoherence matrix R ij (t), but do not change the relative phases of the entangled qubits when their states are probed or manipulated at t = T.
Without any modulations, decoherence in this multiqubit probe has no inherent symmetry, yet one can symmetrize the decoherence by appropriate controls.To this end different, local, phase-locked modulations may be applied to the individual qubits comprising the probe and may be chosen to cause controlled interference between the coupling of these qubits to the bath.The F T,ij (ω) matrices can then satisfy 2N requirements at all times and be tailored to impose the symmetries detailed below.By contrast, a global modulation, i.e., when the same control is applied to all the qubit probes, characterized by F T,ij (ω) = |ε T,i (ω)| 2 is not guaranteed to satisfy N 1 symmetrizing requirements at all times (Figure 4b).
The most convenient symmetry is that of identical coupled particles (ICP), which would emerge if all the controlled qubits comprising the probe could be subject to the following NxN matrix fully symmetrized for their decoherence The result would then be a (N − 1)-dimensional single-excitation decoherence-free subspace (DFS).An initial multi-probe state in this DFS would be preserved for all times, i.e., it would lose neither its population nor its initial entanglement and the fidelity This symmetry requires N(N − 1)/2 conditions which cannot be ensured using N modulating fields.A more limited symmetry is that of independent identical particles (IIP), whereby the N different qubits comprising the probe experience the same single-qubit decoherence R(t) but no cross-decoherence (Figure 4a).This is achievable for ε T,i (ω) ≈ ε T,i δ(ω − ∆ i ).The spectral shifts ∆ i can be different, so as to couple each qubit to a different frequency range of bath modes to ensure that their cross-coupling vanishes: The N single-qubit decoherence rates can be made equal by an appropriate choice of N parameters {∆ i }: where δ i,j is Kronecker's delta.Here, results conveniently express the fidelity of the initial multi-probe state as F(t) = F p (t)F c (t), where F c (t) quantifies the correlation preservation, which is normalized with the population preservation F p (t) [47,48,80].
for dynamically modified relaxation of singly excited multiqubit states.
Both the coupling spectrum matrix, ( ), and the dynamical modulation matrix, , ( ) are then complex and not real, as in the single-qubit scenario.This occurs because, as opposed to the single-qubit probe scenario, the off-diagonal terms of the multiqubit modulation matrix are determined by the phases between modulations of the different qubits comprising the probe.The diagonal elements of the decoherence matrix are easily separable into real decay rates and imaginary Lamb shifts.We first consider a probe comprised of N identical qubits, i.e., set ≡ , subject to impulsivephase modulations (see below Equation ( 23)), requiring AC Stark shifts multiple of 2 at any chosen time t = T.Then, modulations only affect the decoherence matrix ( ), but do not change the relative phases of the entangled qubits when their states are probed or manipulated at t = T. Without any modulations, decoherence in this multiqubit probe has no inherent symmetry, yet one can symmetrize the decoherence by appropriate controls.To this end different, local, phase-locked modulations may be applied to the individual qubits comprising the probe and may be chosen to cause controlled interference between the coupling of these qubits to the bath.The , ( ) matrices can then satisfy 2N requirements at all times and be tailored to impose the symmetries detailed below.By contrast, a global modulation, i.e., when the same control is applied to all the qubit probes, characterized by , ( ) = | , ( )| is not guaranteed to satisfy N ≫ 1 symmetrizing requirements at all times (Figure 4b).(|↓↑ ±| ↑↓ ), coupled to a zero-temperature bath.The overall fidelity F(T) = F p (T)F c (T) is in terms of the correlation preservation F c (T), and the population preservation F p (T).Two qubits comprising a probe in a cavity.
Then, by considering two qubit probes in a cavity as exemplified in Figure 4c, the IIP symmetry, Equation (49), results in the correlation fidelity F c (t) = 1, yet this symmetry does not preclude and still permits population loss, so that the overall fidelity is then the same as the population fidelity, F(t) = F p (t) = e −Re{R(t)} .If the single-qubit decoherence rate R(t) can be suppressed, i.e., if G ii (ω i + ∆ i ) is sufficiently small, this overall fidelity F(t) will be kept close to 1 (Figure 4d).
Thus, the protocol outlined here is to use different modulations for each qubit, and measure the decoherence rate for each such modulation.The result is a vector of decoherence rates that depends on the modulation parameters.We can then find the coupling spectrum that best fits the calculated (F T,ij (ω)) and the measured (R ij (t)) quantities.
A novel feature of a multiqubit probe under control is the different modulation of each qubit, ε j (t) = ε j (t) (Figure 4a).The phase difference of the local modulations allows us to take advantage of interference between different decoherence channels, i.e., cross-decoherence.Different modulations may be designed to cause destructive interference and thus eliminate the cross-decoherence.Conversely, modulations of different qubits may give rise to a decoherence-free-subspace. Hence, control may bring out new distinctive characteristics of the bath cross-correlations between different bath modes or channels, and thus provide more bath-sensing information.

Conclusions
In this mini-review, we have surveyed our universal approach to the dynamical control of both single and multiple qubit probes subject to noise or decoherence by their environments viewed as baths for the purpose of diagnosing the noise spectra, thereby inferring the bath (environment) characteristics or extending the lifetime of the spin-probe state by transferring the information to a more isolated system for improving the resolution of the sensing process [78].Such probes may operate at nanoscales, as magnetometers, thermometers, sensors for imaging or monitoring chemical and biological processes [9,[14][15][16][17][18][19].
Our analysis of dynamically modified coupling between a qubit and a bath is based on the Kofman-Kurizki (KK) universal formula [44][45][46]50,55,56,[60][61][62].We have presented an approach (dubbed BOTOC-bath-optimized task-oriented control) towards finding the optimal time-dependence of the system Hamiltonian, such that a given observable attains a desired value [34,35,81].This approach adapts the spectrum of the system modulation to the bath-coupling spectrum.In general, BOTOC requires information on the peak positions of the bath spectrum.Then, one may optimize the pulse timings and/or continuous control-parameters with respect to a given bath spectrum.Hence, BOTOC encompasses both pulsed and continuous modulation of the system Hamiltonian as limiting cases of control, and thus allows for vastly more freedom compared to control that is restricted to stroboscopic pulses as in some dynamical decoupling (DD) sequences [36][37][38][39][40][41][42]: BOTOC does not require rapidly changing control fields that approximate brief π-pulses.These features allow efficient optimization under a minimal energy constraint.On the other hand, a standard pulse sequence may be preferable experimentally.
Our main point here is that the same control can also be applied to map out the bath spectrum by measuring the decoherence rate for a narrow-band modulation centered at different frequencies [10,11].Furthermore, this control may be optimized based on quantum estimation theory to infer the bath characteristics with the ultimate precision that this theory allows [9] and to selectively address key parameters by filtering out other sources of decoherence [14][15][16][17][18][19].
Our analysis of multiple, field-driven qubits that are coupled to partly correlated baths or undergo local random dephasing allows us to design both global and local control, based on different-but synchronous-local modulations of the different qubits [80].The merits of local vs.global modulations are essentially twofold:

•
Local modulation can effectively decorrelate the different dephasings of the multiple qubits, i.e., eliminate their cross-decoherence, resulting in their equal dephasing rates.For two qubits, the singlet and triplet Bell-states acquire the same dynamically modified decoherence rate.

•
For different couplings to a bath, one can better preserve any initial state by local modulation, which can reduce the mixing with other states, than by global modulation.Local modulation which eliminates the cross-decoherence terms, increases the fidelity more than the global modulation alternative.For two qubits, local modulation better preserves an initial Bell-state, whether a singlet or a triplet, compared to global π-phase "parity kicks." The conclusion from this analysis is that whenever qubits collectively couple to common bath modes (e.g., in ion traps [49] or cavities [47]) local modulations of individual qubits can eliminate cross-decoherence provided they are faster than the correlation (memory) time of the bath, but not necessarily than the much shorter multipartite coherence time.By observing the multipartite disentanglement, as a function of both local and global control, we may infer the cross-correlation between different bath modes (channels), both spatially and spectrally.Since dynamical decoupling (DD) techniques have also been recently employed in noise spectroscopy [82][83][84][85][86][87][88][89], in the wake of our pioneering work [10,11,[13][14][15][16][17][18][19] it is appropriate to compare the two approaches.
Our BOTOC (universal-formula) approach is based on a general framework for optimized design of the modulation, which is desirable in order to accomplish a given task.This optimization consists in adjusting the "score" that quantifies the success of the task, such as the targeted fidelity, purity, entropy, or energy transfer by dynamical modification of the system-bath coupling spectrum on demand [26,27,29].
Several major differences make BOTOC advantageous compared to the DD approach: (i) BOTOC relaxes the DD assumption that the control fields must be either very short or very strong.
In our formalism, the control fields are considered concurrently with the coupling to the bath, hence allowing a much wider variety of pulse sequences, ranging from continuous modulation all the way to DD sequences.(ii) Dynamical decoupling suggests using the same pulse sequence (be it periodic, optimized or concatenated), regardless of the shape of the bath spectrum.By contrast, BOTOC explicitly considers the bath spectrum and allows optimal tailoring of the modulation to a given bath spectrum.In many cases, the standard π-phase "bang-bang" is then found to be inadequate or non-optimal compared to dynamic control based on the optimization of the universal formula.
Whereas our BOTOC approach reduces to the DD method in the particular case of proper dephasing or decay via coupling to spectrally symmetric (e.g., Lorentzian or Gaussian) noise baths with limited spectral width, phase modulation advocated for the suppression of coupling to baths with frequency cutoff or other non-monotonic spectra is, however, drastically different from the DD method which may fail for multipeak spectra.(iii) Our BOTOC universal strategy has far broader applicability than DD: It can simultaneously control, unlike DD [10,11], both decay and decoherence (proper dephasing) by either pulsed or continuous wave (CW) modulation of the system-bath coupling governed by a simple universal formula.BOTOC has been generalized by us to finite temperatures and to qubits driven by an arbitrary time-dependent field, which may cause the failure of the rotating-wave approximation [11].It has also been extended to the analysis of multi-level systems, where quantum interference between the levels may either inhibit or accelerate the decay [19].(iv) Even if DD is adequate for independently decohering qubits, its extension to correlated multipartite systems is highly nontrivial.By contrast, BOTOC modulations with low energy decorrelate the different proper dephasings of the multiple two-level systems (TLS), resulting in equal dephasing rates for all states.For two TLS, we have shown that the singlet and triplet Bell-states acquire the same dynamically modified dephasing rate.This should be beneficial compared to standard DD based on global "bang-bang" (π-phase flips) if both the triplet and the singlet states are used (intermittently) for information transmission or storage.

Open Issues-Outlook
An open issue of the approach is the inclusion of higher orders in the system-bath coupling, which becomes important for strong or resonant system-bath coupling, so that the present perturbative approach cannot be applied.Concatenated DD [57] that treats increasingly higher order corrections of the noise, where each concatenation level of the control pulses reduces the previous level's induced errors, cannot be easily incorporated into our formalism since it goes beyond the second-order approximation used in our derivation of the universal formula.
Another standard assumption that may fail is that the initial state of the system and bath is factorized.Yet, if the system and the bath are in equilibrium, their state is entangled [55][56][57].
Hence, system-bath correlation is an important consideration that requires information on the corresponding parameters.It may be difficult to obtain such data with sufficient experimental precision.A way out is offered by a "closed" loop, based on actual measurements performed on the controlled system in real time rather than on prior knowledge of the bath properties.Such closed-loop control would allow efficient optimization, but at the cost of losing any insight into the physical mechanisms behind the result obtained.
To conclude, the surveyed results show that dynamical control of decoherence is a promising avenue not only for combatting multipartite decoherence, since their rate is not required to become faster as the number of qubits grows, but also for inferring bath spectra.This new avenue may advance the quantum-sensing approach referred to as "quantum noise spectroscopy".

Figure 1 .
Figure 1.(a) Probing bath-induced decoherence by a system probe that is subject to control.Key bath parameters that define the bath spectrum ( ), such as the correlation time τ and the bath-coupling coupling g, are estimated from measurements of the probe in the (anti-)symmetric superposition states.Spectral Filters ( ) designed by: (a) Free evolution (blue sinc filter function), and Car-Purcell-Meiboom-Gill (CPMG) sequence control (green narrow band pass filter); (b) Modulation with small phase shifts, φ ≪ 1; (c) π-pulses (φ = π) modulation; (d) Quasi-periodic modulation.

Figure 1 .
Figure 1.(a) Probing bath-induced decoherence by a system probe that is subject to control.Key bath parameters that define the bath spectrum G(ω), such as the correlation time τ and the bath-coupling coupling g, are estimated from measurements of the probe in the (anti-)symmetric superposition states.Spectral Filters F t (ω) designed by: (a) Free evolution (blue sinc filter function), and Car-Purcell-Meiboom-Gill (CPMG) sequence control (green narrow band pass filter); (b) Modulation with small phase shifts, ϕ 1; (c) π-pulses (ϕ = π) modulation; (d) Quasi-periodic modulation.

Figure 2 .
Figure 2. (a) Schematic view of the quantum state transfer (QST) process from a "noisy"qubit probe, due to the interaction with a bath, to the "quiet" qubit that may serve as robust storage of the information.Controlling the interaction V(t) allows optimized QST; (b) Lowest achievable averaged infidelity as a function of the transfer time T normalized to the infidelity at the fastest transfer time T min at a given transfer energy, plotted for (solid blue line) τ = 0 (Markovian) and (dashed red line) non-Markovian baths τ/T min = 10.For non-Markovian baths (τ T min ) two plateaux (regions of insensitivity to T) are found.The first plateau is independent of the memory time.The second plateau is lower for longer τ.

Figure 3 .
Figure 3. (a) Top: Controlled qubit probes coupled though a symmetric noisy quantum channel.In frequency space, the probes are resonantly coupled by a fermionic-mode of the channel whose energy is defined to be the energy origin (zero).These 3 degrees of freedom constitute the effective system (green rectangular bars); while the remaining channel modes are taken to be the bath (red lines) with a gapped spectrum ± ( ).Under noise, the perturbed energies change/fluctuate bounded by the (colored gapped) Wigner-semicircle.In the gap, the optimal spectral-filters for the high fidelity statetransfer , ( ) correspond to dynamical boundary-control with ( ): = 0 (black dotted), = 2 (green); (b) Transfer infidelity as a function of the transfer time T for the gapped Wigner-semicircle spectral bath: static control = 0 (black dotted) and dynamic control = 2 (green); (c) Averaged transfer infidelity 1 − ( ) as a function of the perturbation strength of the noisy channel.Dynamical = 2 control (empty squares) strongly reduced the infidelity obtained under static control = 0 (empty circles).A fluctuating noisy channel results in less damage than the static-noisy channel → ∞, where its fidelity can even approach its unperturbed value in the Markovian limit → 0 ( = 0, green solid circles) .Here N = 29 channel-qubits and ≡ 1.

Figure 3 .
Figure 3. (a) Top: Controlled qubit probes coupled though a symmetric noisy quantum channel.In frequency space, the probes are resonantly coupled by a fermionic-mode of the channel whose energy is defined to be the energy origin (zero).These 3 degrees of freedom constitute the effective system (green rectangular bars); while the remaining channel modes are taken to be the bath (red lines) with a gapped spectrum G ± (ω).Under noise, the perturbed energies change/fluctuate bounded by the (colored gapped) Wigner-semicircle.In the gap, the optimal spectral-filters for the high fidelity state-transfer F T,− (ω) correspond to dynamical boundary-control with V(t): p = 0 (black dotted), p = 2 (green); (b) Transfer infidelity as a function of the transfer time T for the gapped Wigner-semicircle spectral bath: static control p = 0 (black dotted) and dynamic control p = 2 (green); (c) Averaged transfer infidelity 1 − F(T) as a function of the perturbation strength ξ V of the noisy channel.Dynamical p = 2 control (empty squares) strongly reduced the infidelity obtained under static control p = 0 (empty circles).A fluctuating noisy channel results in less damage than the static-noisy channel τ → ∞ , where its fidelity can even approach its unperturbed value in the Markovian limit τ → 0 (p = 0, green solid circles) .Here N = 29 channel-qubits and V i ≡ 1.

Figure 4 .
Figure 4. Overlap of coupling spectrum ( ) , and modulation matrix elements , ( ) = , * ( ) , ( ) , resulting in modified decoherence matrix elements (overlap), for: (a) crossdecoherence elimination (IIP symmetry); and (b) global modulation (ICP symmetry); (c) Two qubits comprising a probe in a cavity.The qubits are coupled to the cavity modes (thin lines) and subject to local control fields (green lines); (d) Fidelity as a function of the probed time T of the IIP symmetry for the initial entangled two-qubit probe state | (0) = √ (|↓↑ ± |↑↓ ) , coupled to a zerotemperature bath.The overall fidelity ( ) = ( ) ( ) is in terms of the correlation preservation ( ), and the population preservation ( ).Two qubits comprising a probe in a cavity.

Figure 4 .
Figure 4. Overlap of coupling spectrum G ij (ω), and modulation matrix elements F T,ij (ω) = 1 t ε * T,i (ω)ε T,j (ω), resulting in modified decoherence matrix elements (overlap), for: (a) cross-decoherence elimination (IIP symmetry); and (b) global modulation (ICP symmetry); (c) Two qubits comprising a probe in a cavity.The qubits are coupled to the cavity modes (thin lines) and subject to local control fields (green lines); (d) Fidelity as a function of the probed time T of the IIP symmetry for the initial entangled two-qubit probe state |ψ(0) = 1√ 2