We present the moving picture of a molecular bond, in phase-space, in real-time, at resolution limited by quantum uncertainty. The images are tomographically reconstructed Wigner distribution functions (WDF) obtained from four-wave mixing measurements on Br_{2}-doped ice. The WDF completely characterizes the dissipative quantum evolution of the system, which despite coupling to the environment retains quantum coherence, as evidenced by its persistent negative Wigner hole. The spectral decomposition of the WDF allows a direct visualization of wavefunctions and spatiotemporal coherences of the system and the system-bath interaction. The measurements vividly illustrate nonclassical wave mechanics in a many-body system, in ordinary condensed matter.

## I. INTRODUCTION

Dynamical systems are most completely described by their evolution in phase space, which can be represented and visualized by Wigner's quasi-probability distribution function (WDF) along conjugate momentum *p* and position *q* coordinates, *W*(*p*, *q*; *t*).^{1–3} The quasi-probability designation of the WDF underscores that unlike classical distributions it can locally assume negative values, and the negative volume of the WDF serves as a measure of non-classicality.^{4} Given a continuous record of the marginal (shadow^{5}) of the WDF, such as the observable probability density along coordinate space ρ(*q*; *t*) = ∫*W*(*q*, *p*; *t*)*dp*, the WDF can be reconstructed through the inverse Radon transform.^{6} This tomographic procedure is well established in quantum optics,^{7,8} and has been implemented to reconstruct the quantum state of a photon,^{9–12} of ions in an ion trap,^{13} of atoms in an atom beam,^{14} and of diatomic sodium molecules in the gas phase.^{15} Strictly positive WDFs arise in classically coherent collective dynamics.^{16} Superpositions that interfere in phase space exhibit negative Wigner holes are manifestly non-classical; and as such, are broadly designated as Schrodinger “cat” states.^{4} W-holes have been demonstrated in Bose condensates of photons^{11,12} and atoms.^{14,17} Here, we present the tomographic reconstruction of the dissipative motion of a molecule in ordinary matter: bromine in ice. The essential vibronic dynamics of this system has been previously characterized through electronically resonant transient grating (TG) measurements.^{18} The spectrally resolved extension of the same method (SRTG), which we implement here, allows imaging of the marginal ρ(*q*, *t*) = ψ^{*}(*q*, *t*)ψ(*q*, *t*) of the molecular bond after setting it into motion on its electronically excited state. Although measurement effects cannot be entirely eliminated, a faithful reconstruction of the WDF is possible and its negative hole can be demonstrated. We capture and present the complete wave mechanical motion of a bond in real time, with spatial resolution limited by quantum uncertainty, in moving images that vividly illustrate wave-particle duality in molecular matter.

The observations are remarkable because they are made in ordinary condensed matter, in which quantum states are expected to rapidly collapse through decoherence,^{19} effected by energy transfer (dissipation) and loss of deterministic phase evolution (dephasing) as a system interacts with its environment. The observed quantum coherent motion of the Br_{2} bond develops after electronic dephasing and after extensive vibrational dissipation. Despite coupling to the surrounding, the environmentally induced quantum coherence persists at the measurement temperature of T = 120 K, as evidenced by the persistence of its Wigner hole. Since the WDF provides a complete description of the state of the system, it can be spectrally decomposed to directly observe wavefunctions and the spatiotemporal coherence of the system-bath interaction. We visualize the entanglement |*s*, *b*⟩ + |*s*′, *b*′⟩ that appears as the inseparable superposition of the Br_{2} bond and the breathing motion of the cage of water molecules within which it enclathrates.

## II. THE SYSTEM — ORDINARY MATTER

The experiments are conducted on a thin solid film of Br_{2}/H_{2}O obtained by condensing a premixed vapor (1/500 mole ratio) on a sapphire plate held at 120 K. Under these conditions, the solid is expected to be low-density amorphous ice.^{20} Indeed, the visible absorption spectrum of Br_{2}/ice is indistinguishable from that of Br_{2} dissolved in water, as shown in Figure 1. In both cases, the transitions that give bromine its brown color are shifted and broadened. As is typical to ordinary condensed matter, the combination of structural disorder and fast electronic dephasing render the spectra featureless. There is no evidence of a discrete level structure, which is a prerequisite for quantum evolution.

Bromine is known to form clathrate hydrates — a ubiquitous class of crystalline solids in which the guest molecule occupies regular polyhedra made of hydrogen-bonded water molecules.^{21} In crystalline hydrates, the Br_{2} absorption bands undergo moderate shifts^{22} and long Raman progressions are observed, which are characteristic to the various cages the molecule is known to occupy.^{23} In water, the Raman spectrum is quenched.^{24} These differences are explained by local structure.^{25} In water and disordered ice, H_{2}O may bind to Br_{2} via its lone pair of electrons, as in the H_{2}O-Br_{2} complex.^{26} This geometry is not accessible in the clathrate cage where O atoms are fully H-bonded (Fig. 1). That a sub-ensemble of molecules enclathrates in ice was inferred through nonlinear time-resolved spectroscopy.^{18} Independent of initial excitation energy, at wavelengths ranging from 495 nm–565 nm, a coherent vibrational superposition consisting of *v* = 0–2 was observed to form on the electronic B-state, after the loss of 3000–4000 cm^{−1} in energy (Fig. 2(b)).^{18} Measurements in single clathrate crystals show that the wavepacket is prepared during rebound, after the stretching molecule collides with its immediate surrounding, and that the curve crossing located near *v* = 3 limits the highest vibrational level that can be trapped.^{27} It is this indirect, environmentally induced vibrational coherence that we characterize here.

## III. THE MEASUREMENT

As a spectrally resolved, electronically resonant, coherent four-wave mixing process, SRTG probes the phase-space.^{28} The scheme is illustrated in Fig. 2. Two identical non-collinear laser pulses with wavevectors *k*_{1} and *k*_{2} are used in coincidence to prepare a superposition of two vibrational wavepackets |B, *k*_{1}⟩ + |B′, *k*_{2}⟩ on the electronic B-state, and a delayed pulse resonant with the T ← B charge-transfer transition is used to induce coherent scattering between them, *P*^{(3)} = ⟨B′, *k*_{2}|μ(*k*_{4})|T⟩⟨T|μE(*k*_{3})|B, *k*_{1}⟩. This third order polarization over a volume of ∼10^{3} μm^{3} forms the signal which propagates along *k*_{4} = *k*_{1} *− k*_{2} *+ k*_{3}. The measured material response is understood with the aid of the time-circuit diagram (Fig. 2(c)). During *t*_{1}, the excitation evolves as a vibrational density |$\rho _{B,B^\prime } (t_1)$|$\rho B,B\u2032(t1)$. The probe pulse creates the electronic coherence, *ρ*_{B,T}(*t*_{2}; *t*_{1}), which evolves during *t*_{2} to close the measurement loop through radiation *P*^{(3)}(*t*) = Tr[*μρ*_{B,T}(*t*_{2}; *t*_{1})]. The two-time (*t*_{2}; *t*_{1})-correlation that determines *P*^{(3)} maps on the time-frequency domain of conjugate variables (*ω*, *t*_{1}):

The measurement consists of recording the power spectrum of this time correlation along *t*_{2} at a given delay *t*_{1} between writing and reading of the grating.

## IV. SPATIOTEMPORAL COHERENCE

The SRTG data are shown in Fig. 3(a). The oscillatory signal is unmistakably the wavepacket of a damped oscillator, albeit mapped along the spectral coordinate. Slices of the signal at early time show spectral sub-structure due to recursion in the signal. This can be seen in the (*t*_{2}, *t*_{1})-correlation map of Fig. 3(b), which is obtained by Fourier transforming the signal along the spectral coordinate:

The slices along the *t*_{2}-axis provide a detailed dissection of the electronic dephasing process of the probe coherence. They establish that the |B⟩⟨T| coherence is limited to inertial motion (Gaussian decay), with dephasing time equal to the probe pulse width of 60 fs. The recursion at *t*_{2} = 180 fs characterizes the curvature of the T-state by its harmonic frequency: *ω*_{e}(*T*) = 185 cm^{−1}. Frequency analysis along *t*_{1} characterizes the B-state potential by its harmonic frequency *ω*_{e}(*B*) = 173.2 cm^{−1} and anharmonicity *ω*_{e}*x*_{e}(*B*) = 1.75 cm^{−1}. These values are modified in ice in comparison to those of the free molecule: *ω*_{e} = 167.6 cm^{−1} and *ω*_{e}*x*_{e} = 1.08 cm^{−1}.^{29} The similarity of frequencies in B and T states indicates that the molecule retains its bond-order during excitation, consistent with the assignment of the probe resonance to Br_{2}^{−}:H_{2}O_{(s)}^{+}(T) ← Br_{2}(B) charge-transfer transition. The similarity of curvatures in Br_{2}(B) and Br_{2}^{−} potential minima^{30} implies that the difference potential *ΔV*(*q*) = *V*_{T} − *V*_{B} is a linear function of the bond length, *q*, and the Franck-Condon factors are nearly diagonal |$|\langle v({\rm B})|v^\prime (T)\rangle |^2 \sim \delta_{v,v^\prime }$|$|\u27e8v(B)|v\u2032(T)\u27e9|2\u223c\delta v,v\u2032$. Given the linear relation *ω* = *ΔV*(*q*) and the signal is trivially transformed from (*ω, t*_{1})*-*plane to (*q*, *t*_{1})-plane as implemented in Fig. 3(c). As long as ∂Δ*V*(*q*)/∂*q* ≠ 0, the color of the radiation defines the instantaneous bond length. Both turning points of the packet are seen, therefore the full amplitude of motion *Δq* remains within the probe bandwidth Δω = 350 cm^{−1} given by the convolution of pump and probe laser spectra. The requirement for capturing the complete marginal of the evolving density, (∂Δ*V*(*q*)/∂*q*)Δ*q* ⩽ Δω, is fulfilled. Finally, the slope of the difference can be directly read-off the terminal width of the signal, which is given by the thermal width of the nearly stationary state.

To eliminate the interference due to evolution past the probe pulse, we truncate *t*_{2} at 100 fs and back-transform the data:

The filtered signal, plotted against the coordinate axis in Fig. 3(c), shows the motion of the bond with amplitude that decays as it relaxes to form the stationary Gaussian *v* = 0 state. Inspection shows that the oscillating distribution is skewed, and this can be directly traced to the finite time of measurement (see below). The signal in this imposed instantaneous scattering limit is treated as the experimental marginal after correcting for bandwidth, by dividing it with the spectrum at *t* = 0:^{31}

The WDF is constructed through the inverse Radon transform of *ρ*^{e}(*q*, *t*), using the filtered back-transform algorithm over continuous 100-fs time segments that span a half-period of motion.^{32}

## V. THE WIGNER DISTRIBUTION FUNCTION AND ITS PERSISTENT HOLE

Snapshots of the continuously normalized WDF are shown in Fig. 4, and the complete animation is provided in Video S1. As highlighted by the zero-contours, we clearly see a persistent negative hole. A variety of metrics have been defined to quantify non-classicality, of which a convenient measure is the negative volume fraction:^{4}

which is plotted in Fig. 4. The measured W-hole is a periodic function of time, with amplitude that decays from an initial value of 0.3 to 0. Contrary to isolated systems, the interrogated density does not undergo unitary time evolution. It evolves and decays due to contact with the environment. The decay reflects relaxation of the system toward *v* = 0 — a state that does not sustain a hole: |$\delta \_(v = 0) = 0$|$\delta _(v=0)=0$. The periodic modulation is an artifact of the measurement. We illustrate the measurement effect in Fig. 5 by considering pure states of the isolated molecule for which exact calculations are possible. The origin of the modulation is clearest when considering the measurement in W-representation, where the signal is formed by the overlap of the WDF of the B-state and its copy on the T-state. Evolution under the pulse corresponds to rotation of the respective WDFs with relative phase velocity determined by detuning Δ = ℏω − *T*_{e}. As such, the measured distribution is subject to angular convolution under the probe pulse of width *τ*,

Thus, the finite time of measurement leads to modulation, smoothing, and addition of wings to the WDF (see Fig. 5). The extent of distortion depends on details of potentials, laser bandwidth and detuning. However, independent of specifics, the visibility of the hole can only be reduced upon measurement since the full convolution of two W-functions leads to a strictly positive distribution.^{2} In effect, the measured δ_{−} is a lower bound of the true negative volume fraction. We may conclude that the evolving density is non-classical. It retains this character throughout the 4 ps of evolution. Moreover, a weak coherence transfer from the environment can be noticed in the slow growth of the hole near 2.5 ps.

## VI. TWO-SLIT EXPERIMENT IN MANY-BODY WAVE MECHANICS

The nonclassical nature of the motion of the bond is clearly seen. Consider for example the bimodal density along the outer ring in the first frame of Fig. 4. It represents a bond that is simultaneously stretched and compressed, with negative probability amplitude for being at the classical equilibrium geometry. The distribution imitates the Schrodinger cat state of Fig. 4(a), which is constructed from a superposition of two displaced Gaussians:^{4}

The comparison is to a “running cat,” with non-zero momentum: *q*_{0} = 0.8, *p*_{0} = 1.8 and reduced visibility *γ* = 0.5.^{33} This is the initial superposition state in which the bond is prepared with two laser pulses |B(*k*_{1}, *t*)⟩ + |B′(*k*_{2}, *t*)⟩, after rebounding from the cage of water molecules.^{27} The bimodal preparation of the density may be understood to arise from a trajectory ensemble that splits along two paths, e.g., one-bounce and two-bounce paths, to fall below the curve crossing shown in Fig. 2. In the process, the surrounding is driven into motion, and system and bath are entangled |*s*, *b*⟩ + |*s*′, *b*′⟩.^{34} To the extent that the process is sudden, classical coherence subject to phase noise would be expected. This is familiar in both solid^{35} and liquid^{36} phase pump-probe measurements where wavepacket motion is observed, which for the most part can be reproduced using classical trajectory ensembles.^{37} A bimodal classical distribution centered at *q*_{+} and *q*_{−} would be given by the first two terms in (7). The quantum interference between the two paths *W*_{int}, which gives the negative hole, is the nontrivial observation. As in any two-slit experiment, the interference identifies wave mechanics. Here, the interference of the entangled superposition implies phase coherent wave mechanics in a many-body environment. Since the measurement is limited to the Br-Br coordinate, the observable is the wave-mechanical motion of the Br_{2} bond, which can be recognized by comparison to the reference states of the quantum harmonic oscillator provided in Fig. 4(b)–4(d). The early evolution is in a superposition of *v* = 0, 1, 2, and relative amplitudes can be selected to emulate individual frames. As *v* = 2 decays, the quantum equivalent of an oscillating bond appears as the *c*_{0}|0⟩ *+ c*_{1}|1⟩ superposition shown in Fig. 4(c). As *v* = 1 decays, so does the hole (Fig. 4(d)). The spiral distributions seen in Fig. 4 (e.g., at 825 fs and 1375 fs) illustrate the distortion due to angular convolution (6), as in the simulations in Fig. 5. Note that the reference distributions are for pure states. In effect, after a period of evolution during which the system loses ∼0.5 eV of energy, in what must be regarded as classical evolution since it occurs in the continuum of states seen in the absorption spectrum of Fig. 1, a quantum coherent superposition of the system is prepared. Remarkably, while the prepared state dissipates, its quantum interference persists even though at 120 K both system and bath degrees of freedom are thermally accessible. The spectral decomposition, which we consider next, allows a dissection of the underlying dynamical contributions.

## VII. VISUALIZING WAVEFUNCTIONS AND THE SYSTEM-BATH ENTANGLEMENT

Motion in a quantum state is contained in the off-diagonal elements of the density matrix. Given the small number of contributing vibrational states, it is possible to decompose the spatiotemporal coherence into its principle spectral components, as a sum of system and system-bath contributions:

The task is simplified by the nearly harmonic level structure and the prevailing weak damping limit. Accordingly, the spectral components can be retrieved by adding and subtracting a time-shifted copy of ρ^{e}(*q*, *t*) from itself, with time-shift taken at half-period of the characteristic frequencies:^{38} |$\omega _{v,v^\prime } = \omega _s$|$\omega v,v\u2032=\omega s$, 2*ω*_{s}, *ω*_{b}, i.e., at fundamental, overtone, and bath frequencies. The residue yields the non-oscillatory diagonal density |$\sum\nolimits_v {p_v (t)\left| v \right\rangle \left\langle v \right|}$|$\u2211vpv(t)vv$, with slowly varying time-dependent population weights, *p*_{v}(*t*). We present the stepwise dissection of the spatiotemporal coherence and the WDF of components.

The spatiotemporal coherence at the fundamental frequency is isolated by using the fundamental period, *T*_{s} = 194 fs, in the shift-subtract operation:

The extracted map is shown in Figure 6(a), with the abscissa in cm^{−1} units of the raw data. We will refer to this axis as the spatial coordinate, since its linear *ω → q* transformation (which was implemented in Fig. 3(c)) is apparent by inspection. Its spatial sections (Fig. 6(b)) can be recognized as somewhat distorted products of familiar harmonic oscillator wavefunctions φ_{0}(*q*)φ_{1}(*q*) and φ_{1}(*q*)φ_{2}(*q*) of similar appearance. At early time the spatial profile along a horizontal cut is distorted due to evolution under the pulse, and this effect can be removed by considering the profile on a moving frame, i.e., along a skewed line (see Fig. 6(b)). The required slope, ∼100 fs, identifies the slippage to be due to motion of the overtone during the measurement. Consistent with this, the effect disappears as *v* = 2 decays (Fig. 6(b)). The temporal profile (Fig. 6(c)) is modulation at the fundamental beat, as clarified by the power spectrum (Fig. 6(d)). The spectrum establishes the purity of the isolated component. The fundamental beat decays with a time constant of *t* = 1.17 ps.

The shift-sum operation approximates the diagonal density, reinforces the overtone coherence, and adds a phase shift to components that evolve more slowly:

The spatiotemporal coherence at the overtone frequency can be isolated from (9) using *T*_{s}*/*4 time for the shift-subtract operation:

The result is shown in Figure 7. The modulation in time occurs at the center of the spatial distribution, as would be expected for the φ_{0}(*q*)φ_{2}(*q*) cross-term, identifiable by the spatial cuts in Fig. 7(b). The modulation decays with a time constant of 0.57 ps, at nearly half the fundamental lifetime of 1.17 ps. The power spectrum of the time slice (Fig. 7(d)) shows the overtone beat, 2*ω*_{s} = 330 cm^{−1}, flanked by weak sidebands that identify system-bath combination beats. Note that a 2:1 ratio of decay rates between *v* = 2 and *v* = 1 is to be expected for a harmonic oscillator subject to linear coupling to the environment, since γ_{v, v − 1}∝|⟨*v*′|*q*|*v*⟩|^{2} = |⟨*v*′|*a* + *a*†|*v*⟩|^{2} = *v*γ_{1, 0}.

The diagonal density is obtained by the repeated shift sum operation:

The resulting map is shown in Figure 8. It is void of any modulation at Br_{2} vibrational frequencies; however, a faint, slow modulation can be seen in its time profile (Fig. 9(c)). The population decays with a time constant of *t* = 1.9 ps. Its bimodal spatial profile at early time (Fig. 8(b)) identifies it as being dominated by |$\varphi _1^2 (q) + \varphi _2^2 (q)$|$\phi 12(q)+\phi 22(q)$. At *t* = 1.4 ps, the system can be seen to have decayed into |$\varphi _1^2 (q) + \varphi _0^2 (q)$|$\phi 12(q)+\phi 02(q)$. We see the vibrational dissipation through the wavefunctions.

The system-bath coherence is isolated by repeated shift-sum operations, by varying *T*_{b} and demanding maximum visibility:

The result is shown in Figure 9. The time-course of this coherence shows that after an initial swing, it rings for the duration of the population, with a characteristic frequency of *ω*_{b} = 38 cm^{−1}. A second frequency near the overtone of the bath, at 2*ω*_{b} ∼ 80 cm^{−1}, is apparent in the time profile (Fig. 9(b)), and can be seen as the combination band in both system-bath and overtone coherence (Figs. 9(d) and 7(d)). Such persistent lattice motion is usually reserved to zone-boundary phonons.^{39} Although the system is embedded in ice, contrary to phonons of amorphous ice,^{40} the cage motion appears as a nondispersive resonant mode. The observed frequencies are characteristic to breathing modes of the clathrate cage observed in neutron scattering.^{41–43} The observed dynamical localization of the cage modes distinguishes clathrate hydrates from ice and the large disparity in their thermal conductivities. Although an extended clathrate is not expected in this 1:500 guest:host composition, the persistence of the system-bath coherence suggests that we are interrogating a sub-ensemble of enclathrated molecules. We recognize that the local cage is effective in mechanically isolating the molecule from the rest of the environment. Indirectly, the molecule also appears to be electrically screened as if in a Faraday cage. The latter is inferred by noting that the measurements would not be possible if the probe transition were inhomogeneously broadened. Given the large dipole of the charge-transfer transition, its decoupling from dipolar disorder of amorphous ice requires effective screening; presumably, through the rapid mechanical reaction of the hydrogen-bonded cage and clathrate network to dielectric fluctuations in the medium. In effect, the primary bath coordinate that is directly coupled to the system can be identified as the cage-breathing mode, which in turn is poorly coupled to the rest of the solid.

The WDFs of the spectral components are illustrated by phase-space snapshots in Fig. 10 along with reference states of the harmonic oscillator. The WDF of the diagonal density is simply the weighted sum of wavefunctions in W-representation (see Figs. 10(a)–10(c) and Video S2)—once the WDF is reconstructed tomographically, or interferometrically,^{44} wavefunctions become observables.^{45} The comparison with the reference states highlights the fidelity with which the system wavefunctions are observed. The coherences can be seen as the time dependent quantum interference between members of the superposition. The coherence at the fundamental frequency (2nd column in Fig. 10) starts as highly non-classical *c*_{0},_{1}|0⟩⟨1| *+ c*_{1,2}|1⟩⟨2|, and only after dissipation of *v* = 2 does it approach the odd distribution |0⟩⟨1| (see Figs. 10(d) and 10(e), and Video S3), which conforms to the classical notion of a bond that alternatively stretches and compresses. The |2⟩⟨0| overtone coherence (3rd column in Fig. 10), which lasts for the duration of the *v* = 2 lifetime, does not have a classical equivalent — the bond simultaneously materializes at its stretched and compressed extremes (Fig. 6(f) and Video S4). Since this spectral component does not have any interferences, its evolution is a simple rotation of the |2⟩⟨0| phase space distribution (see Video S4). While we clearly see dissipation by the change in composition of the states, there is no evidence of dephasing or decoherence along the system coordinate. Clearly, the system is not isolated. However, the system coherence is isolated through the projective measurement |$\hat P_S$|$P\u0302S$ along the Br_{2} coordinate, enforced by the resonant interrogation of the density, which has the effect of tracing out the bath, *Tr*[|$\hat P_S$|$P\u0302S$*ρ*] = *Tr*_{b}[*ρ*]. To the extent that the total density |*s*_{i}, *b*_{k}⟩⟨*s*_{j}, *b*_{l}| can be decomposed as product states of system and bath, the coherence of the system can be isolated *Tr*_{b}[*ρ*] = *ρ*_{ij}* δ*_{k,l} in diagonal bath states and therefore time evolution strictly given by the system. Identifying the system states as the vibrations of Br_{2}, we obtain the first term in Eq. (8).

The inseparable part of the system-bath superposition, |*s*, *b*⟩ + |*s′*, *b′*⟩, namely the entanglement, is seen as the system-bath coherence in phase space (see last column in Fig. 10 and Video S5). After an initial period where a negative trough is observed, the entanglement appears as the phase portrait of a pure momentum state (as in the last frame of Fig. 10) that alternates sign at every cycle (see Video S5). The interaction appears as periodic momentum exchanges that last throughout the observation time of 4 ps. The coupled motion completely preserves phase coherence. Given the separation of time scales between system and bath, *ω*_{e}/*ω*_{b} ∼ 5, the inseparable motion can be regarded as evolution of the base states of the system, i.e., system states that adiabatically follow a cage that rings throughout the B-state lifetime. This in turn relegates the requirement for preservation of the W-hole to phase coherent evolution of the cage-breathing mode. The observed entanglement |*s*, *b*⟩ + |*s′*, *b′*⟩ is reduced to identifying *b* with vibrations of the immediate cage, which in turn are coupled to the rest of the bath. The extent to which the cage mode induces coherence in the lattice is interesting to consider, since this would be the mechanism for nucleating a “cat”-state that may be macroscopically distinguishable,^{46} e.g., through thermal conductivity. That through entanglement the driven system may induce a more extended coherence in the bath did not materialize in recent semi-classical simulations of a closely related system: I_{2} in a Kr cluster.^{47} This can be shown to occur in the more special case of a superfluid helium bath, where adiabatic evolution of rotational states of a guest molecule can be observed to follow the circulation they induce in the droplet.^{48} Here, it may be argued that the guest and its immediate cage are well isolated from the rest of the lattice, and that it is therefore not too surprising to see quantum coherent evolution of an isolated system of two coupled modes. Their reduction to the model of two bilinearly coupled qubits in a quantum harmonic bath allows exact solutions in the spin-boson framework.^{49} For a sudden disturbance away from thermal equilibrium, the concurrence due to bath coherence, which is closely related to entanglement, very nearly approximates the time profile of the coherence we see in Fig. 9(c).^{50} Rather than quantifying the extent of such system-bath entanglement, here we are satisfied by characterizing its phase-space portrait.

## VIII. CONCLUSIONS

We presented the complete quantum mechanical moving picture (see Video S1) of a molecular bond, captured at a rate of 1 frame every 5 fs, and at resolution limited by quantum uncertainty. The non-classical motion, as certified by the negative volume in its Wigner distribution function, was recorded in ordinary matter through a parametric four-wave mixing process in which a volume coherence on a scale of ∼10^{3} μm^{3} is prepared and interrogated. The quantum coherent response of the ensemble allows it to be viewed as identical copies of one state — a pure state limited in span to the phase space of three system vibrations. Its tomographic reconstruction yields the complete phase-space description of a system undergoing coherent dissipation. Through the spectrally decomposed Wigner distribution functions, we observe the wavefunctions in which the system develops and the spatiotemporal coherence of the system-bath entanglement. The generality of the observation of arrested decoherence is underscored by its realization in rather ordinary matter — impure ice. We show that hidden behind the dephasing window, buried deep in the thermally accessible density of states of the medium, quantum coherent mechanics survives and can be observed with properly tailored measurements. Indeed, such coherence is implied in realizations of interferometric control in condensed matter, the most remarkable of which are demonstrations in the liquid phase.^{51} Our measurements demonstrate them through the direct observation of wavefunctions.

## ACKNOWLEDGMENTS

This research was completed under support of the NSF Center for Chemistry at the Space-Time Limit (CHE-082913). This work has benefited from fruitful discussions with J. A. Cina and N. Schwentner over the years. We thank G. Kerenskaya for preparing Figures 1–3.

## REFERENCES

_{+}⟩ + b|G

_{−}⟩, where |G

_{+}⟩ = G(q+q

_{0},p-p

_{0}) and |G

_{−}⟩ = G(q-q

_{0}, p-p

_{0}) with p

_{0}≠ 0.