Erratum: Nonlinear spectroscopy of controllable many-body quantum systems (2014 New J. Phys. 16 092001)

We establish a novel approach to probing spatially resolved multitime correlation functions of interacting many-body systems, with scalable experimental overheads. Speci ﬁ cally, designing nonlinear measurement protocols for multidimensional spectra in a chain of trapped ions with single-site addressability enables us, for example, to distinguish coherent from incoherent transport processes, to quantify potential anharmonicities, and to identify decoherence-free subspaces.


Introduction
To unravel the spectral and dynamical properties of composite quantum systems of increasing complexity constitutes an indispensable prerequisite for robust control in diverse areas of modern quantum science, from quantum information processing [1][2][3] over photo-induced chemical reactions [4] to the primary processes of light-energy conversion in nature and technology [5,6]. A key ingredient of truly complex quantum systems are strong correlations over broadly distributed energy and time scales [7], which may lead to critical behavior or emergent phenomena and become manifest in such systems' dynamics [8,9].
For a precise characterization of the latter, higher-order correlation functions need to be probed, which poses a formidable experimental challenge [10][11][12]. Much progress has been accomplished in this regard in ultrafast, multidimensional spectroscopy of molecular aggregates, where experiments are often carried out under extremely challenging conditions characterized by short time scales and tight spatial confinement [12,[13][14][15]. In cold-matter experiments, however, where complexity is constructed by controlled assembly of individual components [8,[16][17][18][19][20], Ramsey-type pump-probe techniques are still the preferred experimental tool to probe many-body dynamics [21].
In a Ramsey experiment, the system evolves through a coherent superposition of groundand excited-states in a two-pulse sequence. Such experiments can resolve the energies of contributing excited states, as well as decoherence rates caused by environmental couplings; they are, however, restricted in their combined frequency-time resolution [12], and cannot resolve manifolds of higher excited states or the contribution of individual pathways. For this reason, in chemical physics, one typically employs four-pulse sequences combined with a phase-dependent selection of quantum pathways to resolve involved transport processes, exciton couplings, and conformational changes in large molecular aggregates, such as photosynthetic complexes.
In this paper we combine this systematic construction of measurement protocols of higherorder correlation functions with the extraordinary control and, in particular, the single-site addressability of cold-matter systems [1,22,23]. We elaborate how our general formalism provides unmatched possibilities for the direct experimental assessment of spatially resolved multitime correlation functions, in experiments that rely on current trapped-ion technology. In particular, second-and fourth-order signals will be defined to monitor intricate dynamical features in the vibrational and electronic degrees of freedom.

Nonlinear measurement protocols
Inspired by multidimensional optical spectroscopy, we develop a general formalism to systematically construct nonlinear measurement protocols in quantum many-body systems using a diagrammatic theory. To provide physical intuition, we consider the example of vibrational degrees freedom in a chain of N ions in a linear trap. However, the obtained theoretical formalism is completely general and independent of the system or degree of freedom.
The ion-trap potential can be characterized by trap frequencies ( , , ) x y z ν ν ν , with the confinement much stronger in the transverse than in the axial direction ( : 1 4 . Transverse vibrations along the x-direction can then be described by a tight-binding Hamiltonian of local phonons including a tunneling term [24], 4 We assume that the harmonic confinement along the two uncoupled transverse directions is of the same order of magnitude, i.e., x where a i † denotes the creation operator for a local phonon at site i, a |0 |1 i i † 〉 = 〉. The local trap frequencies and the coupling matrix can be microscopically derived as [24] u u ∑ with tunable strength U can be induced by a standing electromagnetic wave, effectively generating a Bose-Hubbard model with long-range couplings [24]. The basic building blocks for multidimensional spectroscopic protocols are excitation and readout schemes. A vibrational excitation of the form D e e a e a ( ) can be generated by suitable, perturbative pulse sequences [26,27], where 1 α ≪ , and ϕ can be controlled through the pulse parameters. The superoperator describes the corresponding change of the density matrix induced by this weak interaction. In the course of this paper, we consider the impulsive limit for the interaction , which means that the duration of the light pulses is much shorter than the time scale of the vibrational dynamics. This is justified since the characteristic time scale of the Hamiltonian evolution can be adjusted by controlling the distance between neighboring ions via the parameter β [24]. Excitations of local phonons can be read out via the observable (j labels the measured ion), which is accessible by mapping the vibrational population onto the ionʼs electronic state [28] (see appendix A).

( )
Using these ingredients, multidimensional spectra are obtained by scanning the time delays between a sequence of pulses and detecting the observable A j . After a series of m interactions on the ions i i ,..., m 1 , the signal is given by where the expectation value is taken with respect to the non-equilibrium density matrix created by a succession of short, impulsive interaction events and unperturbed time evolution, The latter is described by the Greenʼs function t t ( ) exp =   , where  denotes the Liouville superoperator [12,29], which we assume is time-independent: [29] can be simulated in trappedion systems [30][31][32].
These signals constitute a generalization of Ramsey spectroscopy, and provide the possibility to distinguish the contributions of individual quantum pathways. To see this, note that each time an excitation is created or destroyed due to the interaction ( ) , the phase ϕ k of the kth pulse is imprinted on the quantum state. Hence, different contributions to the total signal can be distinguished by their dependence on these phases. One can exploit this to post-select signals pertaining to certain pathways with distinct combinations of coherences and populations: To this end one employs the phase cycling protocol originally developed in nuclear magnetic resonance [10]. It is implemented by repeating the experiment for a small discrete set of phases followed by Fourier analysis [10,33]. Individual pathways can be represented by Feynman diagrams [12], which offer an intuitive interpretation of the signals (see also figure 1(a)). Further information on coherent pathways and their selection via phase cycling can be found in appendix B.
Let us stress that the way such signals are obtained in controlled quantum systems holds considerable advantages over related methods, such as phase-matched heterodyne ensemble spectroscopy of bulk materials [11,12] or phase-sensitive fluorescence measurements [34], e.g., of single molecules [33]. Specifically, we point out three key differences between these methods and the approach proposed here.
• The first consists of the ability to create localized excitations due to the micrometer separation of the ions, which induces superpositions of a multitude of eigenstates with a single interaction. This is not possible, for instance, when working on molecular aggregates, where interacting chromophores are separated only by few nanometers, two orders of magnitude below the diffraction limit of optical light. • Second, since artificial quantum systems, as considered in the present manuscript, are typically well isolated from the environment, the induced excitations do not decay naturally within the relevant experimental time scales [1]. Instead, the described fluorescence signal must be induced externally via coupling to a short-lived state, providing full control over the delay between readout and the final interaction. Direct field measurements of the Raman scattered light are not suitable for these single quantum systems due to the insufficient number of scattered photons. This is in stark contrast to fluorescence measurements in molecular systems, where fluorescence signals are created by spontaneous decay after a random time, and heterodyne measurements of stimulated emission, which is induced immediately following the last excitation. Therefore, methods that were developed for phase-matched heterodyne electronic spectroscopy cannot be implemented straightforwardly and instead must be adjusted to the current fluorescencebased situation. However, as we will show later in this manuscript, analogs of well-known measurement protocols from electronic two-dimensional (2D) spectroscopy can be defined and interpreted similarly.
• Third, aside from strong spatial confinement, experiments on molecular aggregates are also characterized by extremely challenging time scales. For instance, typical time scales for energy transport and coherence decay on photosynthetic complexes are on the order of picoseconds [5,12]. These parameters render precise experimental control of such molecular aggregates extremely demanding. In contrast, quantum optical experiments, such as trapped ions, allow the study of energy transport under well-controlled conditions on microsecond time scales with milliseconds of coherent dynamics [26,27].
Let us also briefly comment on the difference between single-shot experiments on an ensemble of systems in bulk and experiments on single quantum systems, which are repeated many times to obtain a statistically relevant sample of results. In the latter case, ensemble averaging over slightly fluctuating conditions leads to decoherence effects whose phenomenology is similar to that experienced in bulk materials. In ion-trap experiments, those effects are especially relevant in the context of quantum computation, and consequently have been analyzed thoroughly [35]. For the application of the proposed methods, these effects do not pose a problem, as demonstrated, for instance, by the successful implementation of related methods with single molecules [33].
Making use of the advantages discussed above, we can now design nonlinear signals to target specific dynamical and spectral features of engineered quantum matter systems. Current quantum optical experiments struggle with the efficient analysis of long-range coherences, which are responsible for notoriously elusive quantum critical phenomena [36], and whose role is also currently debated in quantum transport processes [15] which in turn may be studied in quantum simulations [20]. Moreover, in light of their increasing complexity, appropriate methods to certify and diagnose quantum computations are lacking-in particular, for the unambiguous identification of sources of error, such as anharmonic corrections to the trap potential and the precise characterization of detrimental effects induced by external environments [35]. In the following, we show how nonlinear spectroscopy is able to improve on this situation, by demonstrating a selection of applications for examples of second-and fourth order signals.

Coherent and incoherent phonon transport
As a first application, we scrutinize the role of quantum coherence for phonon transport in an ion chain via a second-order signal that monitors the space-and time-resolved spectral decomposition of a local excitation. To this end, we introduce the single quantum coherence (SQC) signal, which is extracted by phase cycling with respect to the phase 1 2 ϕ ϕ − ( figure 1(a)) from the total second-order signal S t t ( , ) i i j , ; (2) 1 2 1 2 . It contains time-correlation functions whose spatial resolution is controlled by the choice of excitation and readout pulses: the first tunable delay t 1 between the excitation pulses at ions i 1 and i 2 and a second tunable delay t 2 before readout at ion j are scanned. To unveil spectral properties, we Fourier-transform the signal with respect to one or both time delays, e.g., In figure 1, we simulate [37] SQC signals for a chain of five ions initialized in their individual ground states |0〉 with a confinement parameter of 0.1 β = . The first two pulses create vibrational excitations at the left-most ion, i.e., i i 1 1 2 = = , which are then probed, after interaction-induced propagation, at the central ion, j = 3. For the interpretation of the spectrum, it is most instructive to represent the local excitations in terms of single exciton states (energy eigenstates of the chain) e | j 〉, e.g., a c e |0 |1 | j j j 1 † 1 1 〉 = 〉 = ∑ 〉, which evolve with their respective eigenfrequencies ω j . We emphasize that the index j now refers to an eigenstate of the chain, i.e., we now work in the exciton basis-in contrast to equation (1), where the Hamiltonian was described on a local-site basis. Since the first pulse induces the coherence |1 0| 〉〈 , which is selected via its phase signature ϕ 1 (see figure 1(a)), the Fourier transform with respect to the first time delay t 1 reveals the single exciton frequencies ω j . During t 2 , the contributions with the phase label 1 2 ϕ ϕ − involve coherences between distinct excitons, such that the Fourier transform of t 2 reveals the energy differences ij i j ω ω ω = − ( figure 1(b)). As we can easily see, these determine the periods at which an excitation coherently tunnels between the ions: the probability for a phonon localized on ion a to tunnel to ion b at a time t is given by p t e c e ( ) | 1 | |1 | , We will therefore denote the frequency differences ω ij as tunnel frequencies. We remark here that the measurement protocol that reveals these frequencies only requires two localized excitation pulses and one readout pulse, independently of the length of the chain.
In the right panel of figure 1(b), we now add local dephasing, described by L a a i i i † γ = [30,31]. Coupling to the environment causes additional transitions, which induce incoherent transport. To understand this, we plot S t | ( , )| 1, 1; 3 (SQC) 1 2 Ω in figure 1(c), which displays the evolution of the initial stateʼs frequency components during t 2 and thereby allows one to monitor the time-resolved exciton transport. Without dephasing (top panel), there is no signal at ω 4 , since the associated breathing mode leaves the central ion immobile [25], and is therefore not detected by local readout. However, in the presence of local dephasing, the breathing mode is incoherently coupled to other modes, which do have a finite amplitude at the central ion. Thus, the bottom panel of figure 1(c) shows the emergence of a signal at ω 4 with increasing t 2 . We observe progressively incoherent (non-oscillatory) transport, leading to a pronounced 2D signal at ( , ) ( , 0) 1 2 4 Ω Ω ω = (peak A). The transient oscillations on top are due to short-lived coherences between the excitons e | 4 〉 and e | 2 〉, and give rise to another, weak 2D signal at ( , ) 4 42 ω ω (peak B). Thus, our local readout scheme provides full information on spatially resolved transport and on the coherently or incoherently populated transporting eigenmodes of the chain.

Detection of anharmonicities
The second order signal employed in our example above can only probe excitonic states within the single excitation manifold (due to the very structure of the impulsive interaction ( ) j ϕ  ), and is therefore unsuitable to detect any anharmoniticity of the phonon spectrum. For the latter purpose, we need to probe the double excitation manifold, and, hence, to find an appropriate fourth-order signal. Such an observable is given by the double quantum coherence (DQC), which has been successfully employed to probe electronic correlations in molecules [38] (see also [39]). For our present purposes, we define the analogous fluorescence-based signal, consisting of the four-pulse pathways with phase signature 1 figure 2(a)). The first two pulses create a double exciton (two-phonon) state f | i 〉, which is subsequently probed by the third and fourth pulses. During t 1 , the created coherence evolves with frequencies ω ω < < …, respectively. Let us consider the example N = 2, where all pulses are applied to the same ion and the times t 2 and t 4 are kept at zero. Note that, as pointed out in the introduction, and in contrast to other methods, this approach allows us to scan the time interval t 4 as well, which can be extremely useful for studies of, for instance, population decay. For the detection of anharmonicities, however, only the time intervals t 1 and t 3 need to be considered. Since for vanishing anharmonicity, U = 0, the energy of each double exciton state f | i 〉 is given by the sum of two single exciton energies, f j k i ω ω ω = + , transitions so that 11 ω ′ and 22 ω ′ coincide in  Ω ω = since only if the first pulse excites e | i 〉, the second pulse can promote the system into f | j 〉. For instance, we do not observe a peak at ( , ) 1 31 ω ω′ , since 2 When we add anharmonicity to the system (U 0 ≠ ), two effects can be observed [figure 2(c)]: First, degeneracies are lifted, and the signals at 11 ω ′ and 22 ω ′ can be resolved. Second, the anharmonicity U perturbs the symmetry of double excitons. This opens up new excitation pathways, e.g., f | 3 〉 can be accessed via e | 1 〉, creating a signal at ( , ) 1 31 ω ω′ . The strength of the anharmonicity U can be inferred by comparing the distance of the ij ω ′-transitions to the possible harmonic transitions, which in turn can be recovered from single-exciton frequencies along the Ω 1 -axis.

2D lineshapes
So far, we have demonstrated how to assess spectral properties of the vibrational degree of freedom of the ion chain. However, each ion also carries an electronic degree of freedom, which can be treated as a spin-1/2 system. The vibrational modes then mediate the coupling between the different spins. Applications of this type of interactions for quantum information purposes potentially suffer from the detrimental influence of environmental noise, what renders a precise characterization of decoherence mechanisms highly desirable in any such experiment. In our last example to illustrate the versatility of nonlinear spectroscopic techniques in combination with single-site resolution, we demonstrate how the line shape of a 2D signal can reveal the nature and strength of the environmental coupling and identify a decoherence-free subspace. We further show how this can be used to certify the fidelity of Bell-state generation.
The formalism for the construction of nonlinear signals introduced before can be readily adapted to other degrees of freedom by identifying appropriate excitation and readout schemes. For spins, we consider the interaction , corresponding to a measurement of the population at spin j [40].
The effective spin-spin interaction is described by the Mølmer-Sørensen Hamiltonian, which for two spins takes the form H The two ions are initialized in their electronic ground state |00〉, and all excitations and the readout are carried out on the same spin. One can readily predict four resonances at ( , , which are eigenstates of H MS with eigenvalues 2 Ω ± [41]. Without dissipation, all four peaks are delta peaks, yet the coupling to an environment can broaden them and gives rise to distinct 2D lineshapes. Figure 3 γ σ σ = ⊗ . Local dephasing broadens the resonances equally along both frequency axes, whereas collective dephasing only affects the width along Ω 1 . The selected quantum pathway evolves in the coherence | | Ψ Ψ 〉〈 ± ∓ during t 2 , which is part of a decoherence-free subspace [42]. For local dephasing, the peak width along Ω 1 allows us to infer the error probability of the corresponding gate U i H exp (  2 ) MS MS π Ω = − , since both of these quantities scale approximately linearly with the dephasing strength γ for 0.1 γ Ω < ( figure 3(b)). The error probability is quantified via the fidelity denotes the state at the outcome of the noisy gate.

Conclusion
We have presented a powerful method to probe spatially resolved multitime correlation functions for the investigation of complex non-equilibrium dynamics and discussed experimental realizations with trapped ions. Our methods generalize previous Ramsey-type techniques inasmuch as they allow for the construction of pulse sequences of arbitrary order and phase-coherently select the contribution of individual quantum pathways.
Our diagrammatic theory provides an intuitive description as well as a new language for cold matter experiments and can be readily extended to include recent developments from nonlinear spectroscopy, such as pulse-shaping and optimal control theory [43,44]. We envision that, beyond trapped ions, implementations with cold Rydberg gases [45,46] and spins of neutral atoms in optical lattices are within experimental reach [21,23]. Ω Ω as read off from two spins initially prepared in the state |00〉 and subject to the Mølmer-Sørensen interaction. The 2D lineshapes clearly distinguish between local and collective dephasing processes. Only under collective dephasing, the system evolves in a decoherence-free subspace during t 2 , as manifest in the delta-shaped signal along Ω 2 . (b) For weak local dephasing ( 0.1 γ Ω < ), the fullwidth-half-maximum (FWHM) along Ω 1 scales approximately linearly with the corresponding Mølmer-Sørensen gate error 1 −  , since both relate to the strength of the noise γ (inset).
We thereby provide a versatile toolbox for controlled quantum systems, opening up a wide range of possibilities to systematically study many-body effects in complex quantum systems, such as the survival of coherences under different environmental and internal couplings, as well as the role of excitonic states in quantum transport processes. The number of pulses required to probe a certain multipoint correlation function is independent of the system size, which renders this method scalable in the limit of increasing particles. This establishes an important step towards the experimental certification of quantum effects in large-scale quantum devices.
D e e a e a ( ) , the created coherences will carry the phase-shift e iϕ ± . All the contributions in a sequence of pulses can be represented by Feynman ladder diagrams, each one representing one excitation pathway with a characteristic phase signature (the combined phase shift of all pulses in the applied sequence), as for instance in figures 1(a) and 2(a). The diagrams are to be read as follows: • Time runs from bottom to top, the left vertical depicts the evolution of the ket, and the right one the evolution of the bra side of the density matrix. • Each excitation (de-excitation) is described by an arrow pointing towards (away from) the density matrix. Each de-excitation adds a factor ( 1) − to the overall sign of the diagram. • To yield a signal, the diagram has to end up in an excited state population when the fluorescence is collected.
The total signal is described by a coherent superposition of all pathways. A complete pedagogical introduction to this formalism is beyond the scope of the current paper, but can be found in a recent review [47], which is specifically written to appeal to quantum opticians, as well as to be appropriate for standard textbooks [11,12].
Phase cycling is a post-processing method that allows us to extract the contribution of subsets of pathways from the total signal, by exploiting their dependence on the phases [10]. To this end, the phases ϕ i of the individual pulses have to be scanned over a discrete set of values.
For example, consider a pulse sequence, where the phase signatures kΔϕ, k k 0, 1 ,..., max = ± ± occur, where in a perturbative treatment of the low-intensity pulses as considered above, k max can typically be assumed to be small (i.e., k 1 max = or 2). From the total signal, which decomposes into a sum of terms with different phase shifts, j max δϕ π = + . Thus, by scanning the phase shifts over a specifically chosen set, followed by an inverse discrete Fourier analysis of the obtained spectra, it is possible to experimentally select the contribution of individual pathways.