Coherent control of a nanomechanical twolevel system
The Bloch sphere is a generic picture describing a coupled twolevel system and the coherent dynamics of its superposition states under control of electromagnetic fields Nielsen2000 . It is commonly employed to visualise a broad variety of phenomena ranging from spin ensembles bib:Vandersypen2005 and atoms bib:Haroche2006 to quantum dots bib:Hanson2008 and superconducting circuits You2011 . The underlying Bloch equations bib:Bloch1946a describe the state evolution of the twolevel system and allow characterising both energy and phase relaxation processes in a simple yet powerful manner bib:Vandersypen2005 ; bib:Yafet1963 ; bib:Burkard2004 .
Here we demonstrate the realisation of a nanomechanical twolevel system which is driven by radio frequency signals. It allows to extend the above Bloch sphere formalism to nanoelectromechanical systems. Our realisation is based on the two orthogonal fundamental flexural modes of a high quality factor nanostring resonator which are strongly coupled by a dielectric gradient field Faust2012a . Full Bloch sphere control is demonstrated via Rabi bib:Rabi1937 , Ramsey bib:Ramsey1950 and Hahn echo bib:Hahn1950 experiments. This allows manipulating the classical superposition state of the coupled modes in amplitude and phase and enables deep insight into the decoherence mechanisms of nanomechanical systems. We have determined the energy relaxation time and phase relaxation times and , and find them all to be equal. This not only indicates that energy relaxation is the dominating source of decoherence, but also demonstrates that reversible dephasing processes are negligible in such collective mechanical modes. We thus conclude that not only but also can be increased by engineering larger mechanical quality factors. After a series of groundbreaking experiments on ground state cooling and nonclassical signatures of nanomechanical resonators in recent years OConnell2010 ; Teufel2011 ; Chan2011 ; PhysRevLett.108.033602 ; 2012arXiv1206.5562P , this is of particular interest in the context of quantum information processing Nielsen2000 ; bib:Ladd2012 employing nanomechanical resonators bib:Stannigel2010 ; 2012arXiv1211.4456R .
While the dynamics of a twolevel system under the influence of a pulsed external electromagnetic field was observed in atomic and nuclear spin physics decades ago, a mechanical analogon to such a system remained elusive for a long time. Only recently, coherent exchange of energy quanta between a mechanical and an electrical mode was achieved: In 2010, O’Connell et al. OConnell2010 managed to control the swapping of a single quantum of energy between a qubit and a mechanical resonator, while Palomaki et al. 2012arXiv1206.5562P demonstrated the temporary storage of itinerant microwave photons in a mechanical resonator in 2012. At the same time, several approaches were employed to achieve purely mechanical resonant coupling either between separate resonators perisanu:063110 ; APEX.2.062202 ; PhysRevB.79.165309 or different modes of the same resonator kozinsky:253101 ; Faust2012a in the classical regime. So far, the pulsed coherent control of the system was prevented by weak coupling, low quality factors or the lack of a sufficiently strong and fast tuning mechanism.
We present the successful implementation of a purely mechanical twolevel system with coherent timedomain control (see also the experiments independently performed at NTT using parametric coupling Okamoto2012 ). To this end, we use a 250 nm wide and 100 nm thick, strongly stressed PhysRevLett.105.027205 silicon nitride beam resonator with a length of 50 µm dielectrically coupled to a pair of electrodes used for detection Faust2012 as well as actuation and tuning rieger:103110 . The two fundamental flexural modes of the mechanical resonator oscillating in the outofplane and inplane direction (see Fig. 1) are coupled by crossderivatives of the strong inhomogeneous electric field generated between the electrodes Faust2012a . A constant dc voltage of 15 V is used to dielectrically tune the system close to the resulting avoided crossing, while the signals generated by an arbitrary waveform generator (AWG) enable timeresolved control vicinal to the anticrossing (see Fig. 1b,c). Both voltages are added and combined with the rf actuation of the beam via a biastee and applied to one electrode. The other electrode is connected to a 3.6 GHz microstrip cavity, enabling heterodyne detection of the beam deflection Faust2012 after addition of a microwave bypass capacitor at the first electrode rieger:103110 . These components as well as the mechanical resonator are placed in a vacuum of mbar and cooled to K to improve the temperature stability as well as cavity quality factor. The microwave cavity is interfaced to the readout with a single coaxial cable and a circulator.
When the system is driven by an external white noise source and the AWG output voltage is swept, the avoided crossing of the two modes shown in Fig. 1c can be mapped out, exhibiting a frequency splitting Hz. With a quality factor and a linewidth of Hz at the resonance frequency , the system is clearly in the strong coupling regime of . For all measurements discussed in the following, an rf drive of 59 dBm at 7.539 MHz, resonantly actuating the beam at an AWG voltage of 0 V, is applied, which initialises the system in its inplane mode (see black circle in Fig. 1c). A 1 ms long, adiabatic voltage ramp up to 2.82 V brings the state to the point of minimal frequency splitting between the coupled modes. Here, the system dynamics is described by two hybrid modes formed by the inphase and outofphase combinations of the fundamental flexural modes. The adiabatic ramp thus transforms all the energy of the inplane mode into the lower hybrid state, such that the twolevel system, consisting of the two hybrid modes, is prepared in its lower state. As the drive frequency remains constant (dashed line in Fig. 1c), the beam is no longer actuated and its energy is slowly decaying.
Now, the application of a continuous pump tone with frequency will start Rabi oscillations bib:Vandersypen2005 ; bib:Rabi1937 between the lower and upper state, as shown in Fig. 2. They can be measured directly by monitoring the time evolution of the output power spectrum at the frequency of one of the hybrid modes, here shown for the upper state at 7.6028 MHz, and measured with a bandwidth of 10 kHz. All timeresolved measurements are averaged over 20 (Rabi oscillations and T1 measurement) or 10 pulse sequences (Ramsey fringes and Hahn echo). For a drive amplitude of 100 mV (half peaktopeak) we find a Rabi frequency of 8.3 kHz (see section II.A of the Supplementary Information for the frequency dependence of the Rabi oscillations). In principle, the decay of these oscillations is governed by both energy relaxation, characterised by a rate , and phase decoherence, characterised by or , where includes reversible processes caused by slow fluctuations or spatial inhomogeneity of the coupling. For clarity, we use these wellknown phenomenological constants in the same way as, e. g., in spin systems bib:Vandersypen2005 , as discussed in more detail in the Supplementary Information section I.
The exponential decay of a state’s energy defines . The corresponding measurement is shown in Fig. 3 for both lower and upper state: The system is once again prepared in the lower hybrid state. To reach the upper state, a subsequent pulse is applied, thus performing one half of a Rabi cycle which transfers the system to the upper state (see Supplementary Information section II for details on the frequency and amplitude calibration of the applied pulses). The exponential decay is then measured directly with a spectrum analyser at a bandwidth of 3 kHz, exhibiting different relaxation times ms and ms for the lower and upper mode, respectively. They correspond to the spectrally measured quality factors. Previously, it has been shown that, at maximum coupling, the two hybrid modes should have the same quality factor and thus time Faust2012a . However, both modes are affected by dielectric damping rieger:103110 , leading to the observed difference.
To measure the time, a pulse is used after the preparation in the lower state to bring the system into a superposition state between lower and upper hybrid mode. The frequency of the pulse is detuned to Hz, leading to a slow precession of the state vector around the zaxis of the Bloch sphere bib:Vandersypen2005 ; bib:Ramsey1950 . As a result, a second pulse after time does not always bring the system into the upper state, but a slow oscillation, the socalled Ramsey fringes, is observed when the delay between the two pulses is varied and the zprojection of the state vector is measured after the second pulse, as shown in Fig. 4. The decay constant of this oscillation is , while the decay of the mean value corresponds to an effective of both modes. The fit in Fig. 4b results in ms and ms. The energy relaxation time of the superposition state is identical to the reciprocal rate average of the two hybrid modes
as the mechanical energy oscillates between the two modes with frequency (see Supplemental Video).
By including an additional pulse at into the Ramsey pulse scheme and replacing the final pulse by an pulse to once again rotate to the upper state (see Fig. 5), the time can be measured in a Hahn echo experiment bib:Vandersypen2005 ; bib:Hahn1950 . The rotation flips the state vector in the xyplane of the Bloch sphere, thus reversing the effects of a fluctuating or inhomogeneous coupling strength in the second delay interval of and thereby canceling their contribution. The frequency of the pulses is once again exactly , as all three pulses need to be applied exactly around the same axis. The resulting decay curve represents , for which a value of ms can be extracted from the fit in Fig. 5b.
The good agreement between and clearly shows that reversible elastic dephasing, e. g. caused by temporal and spatial enviromental fluctuations or spatial inhomogeneities, does not noticeably increase decoherence. Although the experiment is performed with billions of phonons, they all reside in the same collective mechanical mode and thus all experience an identical environment. This strongly constrasts the behaviour found e. g. in spin qubits, where the hyperfine interaction with energetically degenerate nuclear spins causes (see Supplemental Information section III).
It is more surprising that the phase coherence time is equal to the average energy relaxation time . This indicates the absence of measurable elastic phase relaxation processes in the nanomechanical system, such that the observed loss of coherence is essentially caused by the energy decay of the mechanical oscillation (see also Supplemental Information section III). Earlier research PhysRevLett.105.027205 suggests that the dominant relaxation mechanism in silicon nitride strings is mediated by localised defect states of the amorphous resonator material, described as twolevel systems at low temperature. They facilitate energy relaxation by providing the momentum required to transform a resonator phonon into a bulk phonon. For this process to lead to elastic phase relaxation, an excited defect state would have to reemit the phonon back into the resonator mode, which is extremely unlikely due to the weak coupling between the two. In conclusion, we demonstrate coherent electrical control of a strongly coupled () nanomechanical twolevel system, employing the pulse techniques wellknown from coherent spin dynamics in the field of nanomechanics. Each superposition state of the two hybrid modes on the Bloch sphere can be addressed by a sequence of the described pulses. The presented system stands out by the finding that the elastic phase relaxation rate is negligible compared to the energy decay rate , leaving room for improvement of the coherence via increased quality factors.
In light of the recent breakthrough in groundstate cooling of nanomechanical resonators OConnell2010 ; Teufel2011 ; Chan2011 ; PhysRevLett.108.033602 , the coherent manipulation schemes presented here open new applications for nanomechanical systems in quantum information. Not only can they be used as efficient interfaces for quantum state transfers in hybrid quantum systems bib:Stannigel2010 ; bib:Meystre2012 , but by creating coupled, quantised resonators Brown2011 quantum computations can be carried out directly using nanoelectromechanical twolevel systems 2012arXiv1211.4456R .
Acknowledgements
Financial support by the Deutsche Forschungsgemeinschaft via Project No. Ko 416/18, the German Excellence Initiative via the Nanosystems Initiative Munich (NIM) and LMUexcellent, as well as the European Commission under the FETOpen project QNEMS (233992) is gratefully acknowledged. We thank G. Burkard for his comments on decoherence in a threelevel system and H. Okamoto, I. Mahboob and H. Yamaguchi for critically reading the manuscript.
Competing Interests
The authors declare that they have no competing financial interests.
Author Contributions
J.R. and M.J.S. designed and fabricated the sample, T.F. conducted the measurements and analysed the data. T.F., J.P.K. and E.M.W. wrote the paper with input from the other authors, the results were discussed by all authors.
Correspondence
Correspondence and requests for materials should be addressed to E.M.W. (email: weig@lmu.de).
References
 (1) Nielsen, M. A. & Chuang, I. L. Quantum Computation and Quantum Information (Cambridge University Press, 2000).
 (2) Vandersypen, L. M. K. & Chuang, I. L. NMR techniques for quantum control and computation. Rev. Mod. Phys. 76, 1037–1069 (2005).
 (3) Haroche, S. & Raimond, J.M. Exploring the Quantum: Atoms, Cavities, and Photons (Oxford University Press, USA, 2006).
 (4) Hanson, R. & Awschalom, D. D. Coherent manipulation of single spins in semiconductors. Nature 453, 1043–1049 (2008).
 (5) You, J. Q. & Nori, F. Atomic physics and quantum optics using superconducting circuits. Nature 474, 589–597 (2011).
 (6) Bloch, F. Nuclear Induction. Phys. Rev. 70, 460–474 (1946).
 (7) Yafet, Y. g Factors and SpinLattice Relaxation of Conduction Electrons. Solid State Physics 14, 1–98 (1963).
 (8) Burkard, G., Koch, R. H. & DiVincenzo, D. P. Multilevel quantum description of decoherence in superconducting qubits. Phys. Rev. B 69, 064503 (2004).
 (9) Faust, T. et al. Nonadiabatic Dynamics of Two Strongly Coupled Nanomechanical Resonator Modes. Phys. Rev. Lett. 109, 037205 (2012).
 (10) Rabi, I. I. Space Quantization in a Gyrating Magnetic Field. Phys. Rev. 51, 652–654 (1937).
 (11) Ramsey, N. F. A Molecular Beam Resonance Method with Separated Oscillating Fields. Phys. Rev. 78, 695–699 (1950).
 (12) Hahn, E. L. Spin Echoes. Phys. Rev. 80, 580–594 (1950).
 (13) O’Connell, A. D. et al. Quantum ground state and singlephonon control of a mechanical resonator. Nature 464, 697–703 (2010).
 (14) Teufel, J. D. et al. Sideband cooling of micromechanical motion to the quantum ground state. Nature 475, 359–363 (2011).
 (15) Chan, J. et al. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 478, 89–92 (2011).
 (16) SafaviNaeini, A. H. et al. Observation of Quantum Motion of a Nanomechanical Resonator. Phys. Rev. Lett. 108, 033602 (2012).
 (17) Palomaki, T. A., Harlow, J. W., Teufel, J. D., Simmonds, R. W. & Lehnert, K. W. State Transfer Between a Mechanical Oscillator and Microwave Fields in the Quantum Regime. ArXiv eprints (2012). eprint 1206.5562.
 (18) Ladd, T. D. et al. Quantum computers. Nature 464, 45–53 (2012).
 (19) Stannigel et al. Optomechanical transducers for longdistance quantum communication. Phys. Rev. Lett. 105, 220501 (2010).
 (20) Rips, S. & Hartmann, M. J. Quantum Information Processing with Nanomechanical Qubits. ArXiv eprints (2012). eprint 1211.4456.
 (21) Perisanu, S. et al. The mechanical resonances of electrostatically coupled nanocantilevers. Applied Physics Letters 98, 063110 (2011).
 (22) Okamoto, H., Kamada, T., Onomitsu, K., Mahboob, I. & Yamaguchi, H. Optical Tuning of Coupled Micromechanical Resonators. Applied Physics Express 2, 062202 (2009).
 (23) Karabalin, R. B., Cross, M. C. & Roukes, M. L. Nonlinear dynamics and chaos in two coupled nanomechanical resonators. Phys. Rev. B 79, 165309 (2009).
 (24) Kozinsky, I., Postma, H. W. C., Bargatin, I. & Roukes, M. L. Tuning nonlinearity, dynamic range, and frequency of nanomechanical resonators. Applied Physics Letters 88, 253101 (2006).
 (25) Okamoto, H. et al. Coherent phonon manipulation in coupled mechanical resonators. submitted to Nature (2012).
 (26) Unterreithmeier, Q. P., Faust, T. & Kotthaus, J. P. Damping of Nanomechanical Resonators. Phys. Rev. Lett. 105, 027205 (2010).
 (27) Faust, T., Krenn, P., Manus, S., Kotthaus, J. P. & Weig, E. M. Microwave cavityenhanced transduction for plug and play nanomechanics at room temperature. Nat Commun 3, 728 (2012).
 (28) Rieger, J., Faust, T., Seitner, M. J., Kotthaus, J. P. & Weig, E. M. Frequency and Q factor control of nanomechanical resonators. Applied Physics Letters 101, 103110 (2012).
 (29) Meystre, P. A short walk through quantum optomechanics. arXiv:1210.3619 (2012).
 (30) Brown, K. R. et al. Coupled quantized mechanical oscillators. Nature 471, 196–199 (2011).
Supplemental Material to ”Coherent control of a nanomechanical twolevel system”
I Relaxation times
In the main text, the relaxation constants and are introduced phenomenologically to define the exponential decay times extracted from the energy relaxation and Hahn echo experiments. But different from the twostate spin systems associated with experiments on the Bloch sphere, the mechanical system investigated here actually has three states: the two coupled hybrid modes, i. e. the lower and upper state as well as the ground state (the thermally occupied phonon bath) into which a phonon can relax from either mode, see Fig. S1.
For one, this makes it necessary to prepare the lower state prior to any measurement, as none of the states of the nanomechanical twolevel system is automatically populated (except for the comparatively weak thermal excitation). Furthermore, it introduces additional terms to the Bloch equations: Assuming two independent decay rates and for the upper and lower state and following reference , the rotatingframe Bloch equations can be written as
(S1)  
(S2)  
(S3) 
Here, is the detuning between drive frequency and coupling strength and reflects the drive strength and corresponds to the frequency of Rabi oscillations. is the average decay rate, the relaxation time in the equator plane of the Bloch sphere, the relaxation time along the z direction and , and denote the respective components of the state vector. The phenomenological parameters and can be identified as the sum of the respective rates.
The measured values of are consistent with the quality factors of the corresponding modes, which are only limited by and . This implies a negligibly small . As the measurements show that is equal to , must also be negligible. Thus, the two coherence times and are solely limited by the average mechanical damping of the two resonator modes.
Reference 31 also introduces an additional term in equation (S3) taking into account the difference between the two decay rates and . It leads to a slow tilt of the state vector of a superposition state away from the equator plane towards the state with the smaller decay rate. However, this effect plays no role in the classical system presented here: In a classical superposition state, the energy is distributed between both modes. The system performs an oscillation between its two fundamental modes, changing the direction of rotation with the precession frequency. Thus, the time spent in each mode is equal, and the state just experiences the average decay constant. A plot of the phase space trajectory for one precession cycle is shown in Fig. S2, using a 10 times exaggerated coupling strength. The horizontal and vertical axes correspond to an inplane and outofplane oscillation, while a diagonal motion is associated with the lower and upper hybrid mode. An animated version of this plot is available as an ancillary file on arXiv.org.
Ii Pulse calibration
The first step in characterising the system is to measure an avoided crossing as shown in Figure 1 of the main text. From its fit, the approximate frequency splitting and the AWG voltage required to adiabatically tune to the lower state is extracted .
ii.1 Pulse frequency
To precisely determine the correct pump tone, the fitted frequency splitting is not accurate enough. Instead, the frequency of Rabi oscillations is monitored while sweeping the pump frequency. The quadratic dependence for small detunings , as shown in Fig. S3, allows to fit the measured points and extract the lowest Rabi frequency and thus the exact pump frequency corresponding to zero detuning.
ii.2 Pulse length and spacing
The pulses applied to the system should be as short as possible to allow fast control sequences. The lower limit of the pulse length is one period of the pump signal, as abrupt voltage jumps of a chopped up sine wave will disturb the system. The shortest applied pulse (a pulse) is thus set to a duration of 1/(pump frequency) so that it consists of exactly one sine wave. The other pulses are correspondingly longer, as shown in the pulse schemes in Figures 35 of the main text.
The delay time between two pulses in Figure 4 and 5 also needs to be a multiple of 1/(pump frequency), as otherwise the second rotation will not be carried out along the same axis. This is apparent from the Ramsey fringe experiment, where the observed beating pattern, and thus the effective rotation angle, is caused by an intentional detuning of 500 Hz.
ii.3 Pulse amplitude
After both length and frequency of a pulse are fixed as described above, its amplitude has to be adjusted to define the rotation angle achieved with each pulse. To this end, four sine periods (i. e. a pulse at the desired amplitude) are applied to the system, varying the amplitude. The population of the upper state is measured. At zero amplitude, the pulse has no effect and all energy remains in the lower state. With increasing amplitude, the achieved rotation angle increases, and the first minimum corresponds to the desired character of the pulse. The population of the upper state can be described as
(S4) 
where A is the pulse amplitude and the amplitude corresponding to a rotation of . A fit to the measured data, as shown in Fig. S4, can be used to extract the correct pulse amplitude. This measurement has to be repeated with the detuned pump frequency used for the Ramsey fringes experiment, as the slightly different frequency also leads to a small shift in the amplitude of the pulse.
Iii Phase relaxation
To define the elastic phase relaxation rate two situations have to be distinguished: In systems where energy relaxation occurs only from the upper to the lower level (e.g. references ) the phase relaxation rate is defined as
(S5) 
In contrast to that, in a system dominated by spontaneous decay to a third state and no relaxation between the upper and lower level (i. e. and finite , see equation (S3)), the phase relaxation rate , in this case equal to (defined in equation (S1) and equation (S2)), is given as
(S6) 
This is the case for the system presented here. As and are equal within the measurement acurracy, can not be determined from the experiment. This shows that the measured phase decoherence is solely caused by energy relaxation. Processes changing the phase but preserving the state’s energy seem to play no role.
To compare this mechanical twolevel system to other coherent systems, it helps to take a general look at the possible decay processes: Inelastic processes in which energy is transfered to a thermal bath are irreversible. They are directly represented in the time and also pose a limit to via the two above equations S5 or S6. Irreversible elastic interactions lead to a nonzero and thus reduce , whereas reversible phase decay processes can be measured and controlled e. g. by a Hahn echo experiment and only decrease .
In most coherent nanoscale solidstate systems, the coupling to a fluctuating thermal bath of phonons, photons, twolevel systems or (nuclear) spins leads to one or more of the above processes. For example, superconducting qubits suffer from flux, charge and photon noise. In gatedefined spin qubits e. g. nuclear spin and phonon interactions limit the performance, whereas NV centres in diamond couple to the surrounding nuclear spin bath. In the amorphous dielectric system presented here, relevant loss mechanisms occur via defect states with a broad energy spectrum often associated with twolevel systems , and the phononic environment, as long as no additional electronic noise is introduced via the measurement devices and tuning voltages. As the mechanical modes under investigation are situated within the suspended beam, they effectively reside inside a phonon cavity and couple extremely weakly to the phonon bath of the bulk sample via the narrow clamping points . The exchange of energy of the discrete longwavelength resonator modes and the continous shorterwavelength phonon spectrum of the beam is found to be mediated by the defects. As only higher energy phonons with small wavelengths can effectively transmit energy through the clamps to the bath, scattering of thermally excited higher energy phonons at defect states in a three particle interaction is the most likely process. These inelastic processes destroy energy as well as phase and likely explain why we find .
Supplementary References

Vandersypen, L. M. K. & Chuang, I. L. NMR techniques for quantum control and computation. Rev. Mod. Phys. 76, 1037–1069 (2005).

Burkard, G., Koch, R. H. & DiVincenzo, D. P. Multilevel quantum description of decoherence in superconducting qubits. Phys. Rev. B 69, 064503 (2004).

Faust, T. et al. Nonadiabatic Dynamics of Two Strongly Coupled Nanomechanical Resonator Modes. Phys. Rev. Lett. 109, 037205 (2012).

Pöttinger, J. & Lendi, K. Generalized Bloch equations for decaying systems. Phys. Rev. A 31, 1299–1309 (1985).

Hu, X., de Sousa, R. & Sarma, S. D. Decoherence and dephasing in spinbased solid state quantum computers. Proceedings of the 7th International Symposium on Foundations of Quantum Mechanics in the Light of New Technology, eds.Yoshimasa A. Ono, K. Fujikawa und Kazuo Fujikawa, World Scientific (or: condmat/0108339) 3–11 (2002).

Drake, G. W. F. (ed.) Springer Handbook of Atomic, Molecular, and Optical Physics, page 1004 (Springer, 2006).

Burkard, G. Private communication (2012).

Houck, A., Koch, J., Devoret, M., Girvin, S. & Schoelkopf, R. Life after charge noise: recent results with transmon qubits. Quantum Information Processing 8, 105–115 (2009).

McDermott, R. Materials origins of decoherence in superconducting qubits. IEEE Transactions on Applied Superconductivity 19, 2 – 13 (2009).

Rigetti, C. et al. Superconducting qubit in a waveguide cavity with a coherence time approaching 0.1 ms. Phys. Rev. B 86, 100506 (2012).

Petta, J. R. et al. Coherent manipulation of coupled electron spins in semiconductor quantum dots. Science 309, 2180–2184 (2005).

Reilly, D. J. et al. Suppressing spin qubit dephasing by nuclear state preparation. Science 321, 817–821 (2008).

Bluhm, H. et al. Dephasing time of GaAs electronspin qubits coupled to a nuclear bath exceeding 200 µs. Nat Phys 7, 109–113 (2011).

Hanson, R., Kouwenhoven, L. P., Petta, J. R., Tarucha, S. & Vandersypen, L. M. K. Spins in fewelectron quantum dots. Rev. Mod. Phys. 79, 1217–1265 (2007).

Takahashi, S., Hanson, R., van Tol, J., Sherwin, M. S. & Awschalom, D. D. Quenching spin decoherence in diamond through spin bath polarization. Phys. Rev. Lett. 101, 047601 (2008).

Pohl, R. O., Liu, X., Thompson, E. Lowtemperature thermal conductivity and acoustic attenuation in amorphous solids. Rev. Mod. Phys. 74, 991–1013 (2002).

Cole, G. D., WilsonRae, I., Werbach, K., Vanner, M. R. & Aspelmeyer, M. Phonontunnelling dissipation in mechanical resonators. Nat Commun 2, 231 (2011).