Projecting an ultra-strongly-coupled system in a non-energy-eigenbasis with a driven nonlinear resonator.

We explore the problem of projecting the ground-state of an ultra-strong-coupled circuit-QED system into a non-energy-eigenstate. As a measurement apparatus we consider a nonlinear driven resonator. We find that the post-measurement state of the nonlinear resonator exhibits a large correlation with the post-measurement state of the ultra-strongly coupled system even when the coupling between measurement device and system is much smaller than the energy scales of the system itself. While the projection is imperfect, we argue that because of the strong nonlinear response of the resonator it works in a practical regime where a linear measurement apparatus would fail.

Introduction

A quantum measurement typically projects the quantum state of a system into an eigenstate of a measured observable . In quantum measurement theory, the measurement apparatus interacts with the target system due to an interaction Hamiltonian , where denotes an operator of the apparatus and denotes a coupling strength[1,2]. After allowing the system and apparatus to interact for some time, they become strongly correlated. A subsequent measurement on the apparatus itself implements the desired projection of the target system, and the readout of the apparatus is associated with the eigenvalues of the system observable . The relationship between the operator and the natural Hamiltonian describing the system is typically important. For example, to realize a quantum non-demolition (QND) measurement, a minimal condition is that the observable should commute with the target system Hamiltonian[2-6].

Alternatively, as a means of control, in addition to measurement, it is interesting to consider strongly non-QND measurements[7-9]. Surprisingly, such non-energy-eigenbasis projective measurements are sometimes not straightforward. Ideally, if the system observable to be measured does not commute with the system Hamiltonian, must have matrix elements to induce transitions between the energy eigenstates. Importantly, however, if the coupling strength is much smaller than the energy-level separation of the system, such transition matrix elements disappear under a rotating-wave approximation[10] (see Appendix A for details), and we cannot project the system into the eigenbasis of ; the system remains only minimally perturbed by the measurement apparatus, and as such stays in its energy eigenbasis.

On the other hand, in the ideal case, if the coupling between the system and apparatus is much larger than the system energy, one can perform a projective measurement much faster than the typical time scale of the system, hence realizing the desired non-energy eigenbasis measurements[11-13]. However, if these energy scales are comparable, the dynamics, and the subsequent quantum measurement process, become much more complicated than the cases described above. Understanding the interaction between the apparatus and system, and their dynamics, is important not only for explaining the mechanism of quantum projective measurements away from ideal parameter regimes, but also to achieve a higher level of control over quantum systems.

The ultra-strong coupling (USC) regime between atoms and light[14,15] is especially attractive and practical areas in which to explore the possibility of non-energy eigenbasis measurements. This is because the ground state of this system exhibits non-trivial entanglement between the atom and photons, and virtual excitations, which are difficult to probe with energy eigenbasis measurements alone. If one cannot dynamically control the coupling strength between light and matter due to experimental limitations, being able to measure and control the non-eigenstate structure of such systems would allow one to more carefully validate the properties of the system, choose between competing models of the same physical phenomena[16], and manipulate them for quantum information purposes.

In addition, when the coupling strength between light and matter becomes extremely strong (also known as ultra-strong coupling (USC) limit), so that it is comparable with the cavity resonance frequency, it is predicted that a new ground state will emerge[17-35]. Such a regime was recently experimentally demonstrated[36-38], and has the potential to be realized in a large variety of devices[14,15], including hybrid systems with artificial enhancement of coupling and collective effects[39-47]. Interestingly, non-eigenbasis measurements on a ground state of a system in the USC regime could potentially be used to induce an optical cat state, which is itself a resource for quantum information processing[22,23].

Given such potential benefits, and open problems to be solved, this ultra-strongly-coupled system is attractive as an example with which to investigate the problem of non-energy eigenbasis measurements. Although there are several previous works studying how to detect virtual-excitations in the ground state in an ultra-strongly-coupled system[17,19-21,25,48,49], here we focus only on how to perform non-energy eigenbasis measurements on the ground state of such a system. This would allow us, for example, to gain some information of the weights of the different components in the ground-state, as well as projecting the system into superspositions of eigenstates.

In this paper, we specifically analyze the full dynamics of an ultra-strongly coupled system interacting with a nonlinear measurement apparatus, in the above described situation (where the measured system observable does not commute with the system Hamiltonian). We evaluate the dynamics of the measurement apparatus during the interaction period, the back-action of the measurements on the system, and the correlations that build up between the system and the apparatus.

Such properties are typically studied when one tries to examine in detail a quantum measurement process[50-53]. Although there exist theoretical proposals to use a detector that continuously monitors the system[12], here we consider a binary-outcome measurement performed on the measurement apparatus after the measurement apparatus and system have been allowed to interact. Such a binary-outcome measurement is understood to induce a strong correlation with the system[54,55].

While linear resonators are used as a standard method for quantum measurement in cavity quantum electrodynamics and circuit quantum electrodynamics, in some cases a nonlinearity has been employed to improve qubit readout[4-6,54,56,57]. Due to the bifurcation effect, the state of the nonlinear resonator becomes highly sensitive to the state of the system, which enables one to implement a high-visibility readout.

Here, full numerical modelling and a low-energy approximation allow us to understand how a driven nonlinear resonator interacting with the ultra-strongly-coupled system operates, albeit with limited fidelity, in a regime where one might expect it not to work at all. This is surprising because, although the coupling between the nonlinear measurement device and the ultra-strongly-coupled system is weak compared to system energy scales, we find that, on reasonable time-scales, a strong correlation between them builds up, leading to a partial non-energy-eigenbasis projection of the USC system. We find that such a non-trivial strong correlation comes from the AC Stark effect induced by the driven nonlinear resonator.

To understand better this build up of correlations, we evaluate how quantum correlations, such as entanglement and quantum discord, are generated between system and measurement device during their interaction period. These additional results reveal the conditions neccessary for a non-energy eigenbasis measurement to be realized with this system.

The remainder of this paper is organized as follows. First, we introduce the USC system and its ground state. Second, we discuss the nature of interaction between the nonlinear resonator and the USC system, and we introduce a coarse-graining measurement of the nonlinear resonator itself. Third, we present numerical results to show how a strong correlation arises, even in a parameter regime where the coupling strength may be incorrectly considered to be negligible. A two-level approximation of the USC system allows us to understand how to gain information on the ground-state structure, what it means for the USC system to project onto superpositions of eigenstates, and the relationship between the measurement process and the AC-stark effect.

Ultra-strong Coupling Between Light and Matter

The Hamiltonian of light in a single-mode cavity ultra-strongly coupled to matter (where the matter is well described by a two-level system) is, in its simplest form, given by the Rabi model[58]where () is an annihilation (creation) operator for the single-mode cavity/resonator, () denotes the qubit (resonator) frequency, and is the coupling strength between resonator (light) and qubit (matter).

Recall that, when the matter is in the form of a superconducting flux qubit, as in the recent ultra-strong coupling experiments in[36-38], is diagonal in the persistent-current basis of and of the superconducting flux qubit. Here, we assumed that the flux qubit is operated at the symmetric point.

Throughout this paper we assume that the qubit frequency is much smaller than the resonator frequency, allowing us later to use an adiabatic approximation. In this case, in the limit , we can approximately write the ground state of this system as[17]where is the ratio of the coupling strength and resonator energy . As an example, using parameters close to those used in[36], we plot the function of the reduced density matrix of where the atom is traced out in Fig. 1. The definition of the function for a state is , where is a coherent state for a complex number . We plot the real part of in the axis while we plot the imaginary part of in the axis.

Figure 1

The Q function of the reduced density matrix of the cavity in the ground state , where . Here , , .

In this regime, it is straightforward to understand that spectroscopic measurements can excite transitions between eigenstates, and give information on the energy-level structure. However, it is difficult to, for example, probe the relative weights of the and states in the ground-state, or project the entire system into a state which is not in the eigenbasis, because the energy scales of the system are so large compared with the coupling strength between the nonlinear resonator and the flux qubit.

Using a Nonlinear Resonator As a Measurement Device

Here, as a measurement apparatus, we consider a driven nonlinear resonator whose energy depends on the state of the qubit. It is well understood that a nonlinear resonator can exhibit bistability[50,52,59-61], which makes such a device sensitive to small changes in external fields. In addition, the nonlinearity induces a rapid change in the photon number under driving[50], compared to the linear case. When used as a measurement device, the fast evolution and the sensitivity of the steady-state to weak fields result in a strong and fast correlation of the nonlinear resonator state with the qubit being measured, potentially giving a means to implement a rapid projective measurement. One should note that, typically, the state of the nonlinear resonator is itself measured by standard homodyne techniques[58], and this measurement provides the information about the qubit state.

It is worth mentioning that there are some theoretical proposals[11] to treat such a measurement device as a two level system when the measurement outcomes are binary. However, since such a simplification cannot easily quantify the strength of the correlation between the target qubit and measurement apparatus during the measurement process, we now model the measurement apparatus with a proper Hamiltonian as we will describe below.

The total system, composed of the ultra-strongly-coupled light-matter system, and the nonlinear resonator measurement device, can be described by the Hamiltonian in the laboratory frame[17,20,21,50,52,60,61]where is the annihilation operator of the nonlinear system, denotes the detuning between the nonlinear resonator energy and driving frequency, is the nonlinearity strength, denotes the driving strength of the nonlinear resonator, and is the driving frequency of the nonlinear resonator. In addition, is the coupling between the qubit and the nonlinear resonator, which is not derived from the dispersive approximation to a dipole coupling, but is intrinsic (see Appendix B for details.). Depending on the state of the qubit, the frequency of the nonlinear resonator changes. We set the parameters of the system such that when the qubit is in , the nonlinear resonator experiences the bifurcation effect to generate the high amplitude state. On the other hand, when the state is in , the state of the nonlinear resonator stays in the low amplitude state.

It is possible to activate the effective interaction between the measurement apparatus and ultra-strongly coupled system by starting driving the nonlinear resonator, because the vacuum state of the nonlinear resonator without driving makes the effective interaction negligible. In the rotating frame defined by  , and by applying the rotating-wave approximation, we have

In order to include the loss of photons from the nonlinear resonator, we adopt the following Lindblad master equation, valid when the coupling between nonlinear resonator and its environment is weak, and when the coupling between the nonlinear resonator and the qubit is weak[50,52,60,61]where denotes the photon leakage rate from the nonlinear cavity. The potential losses from the ultra-strongly coupled system are described later.

Coarse-graining of the measurement outcome

After the qubit and the measurement apparatus have been allowed to interact for some time (see Fig. 2), we need to perform a measurement on the measurement apparatus itself. Ideally, one could apply a projection operator on the nonlinear resonator, where is an eigenvector of the quadrature operator . However, due to imperfections in the measurement setup, one cannot resolve arbitrarily small differences in the state of the resonator. Normally, to describe more realistically the measurement process, one takes this into account by considering the integrated signal-to-noise[62], where the noise can include contributions from vacuum fluctuations and noise in the measurement apparatus itself. Here, instead we employ a “coarse graining” approximation described by the following operator with Gaussian noisewhere is the width of the error of the measurement process, and the post-measurement state is described by . Similar coarse graining approaches have been made in refs. [63,64]. This approach allows us to consider the transition from small to large noise situations without being specific about the source of the noise.

Figure 2

A schematic of the measurement process. The USC system and the nonlinear resonator are allowed to interact for some time t, after which the resonator is projected onto two different quadratures. We then analyze the correlation between the conditioned USC system state and the resonator.

Correlations between the nonlinear resonator and the qubit should occur after they have interacted for some time, and, for the parameter regime we use in this work, typically the nonlinear resonator state with () corresponds to an outcome where the qubit was initially in its excited (ground) state. We can describe the post-measurement state of the ultra-strongly-coupled (USC) system as (see Appendix C for details)where erfc is the complementary error function and and N′ are normalization factors.

In the limit when , we obtain . In this case, the measurement results do not contain any information of the post-measurement state of the qubit. On the other hand, we obtainandin the limit , which corresponds to an ideal projective measurement that can perfectly distinguish or .

Full Dynamics of the USC System and Nonlinear Measurement Device

Using parameters from[36], we numerically[65,66] solve Eq. (2), with the following measurement protocol (see Fig. 2): We assume the USC system is in the initial state , and that the nonlinear measurement apparatus is in its uncoupled and undriven ground state. We then allow them (the USC system and the measurement apparatus) to interact for a time up to . We then apply an instantaneous coarse-grained measurement to the measurement system.

To begin, in Fig. 3, we show how the function of the resonator part of the USC system depends on the resolution of the coarse-graining measurement. One sees that the change in the state of the USC resonator is much stronger when the measurement resolution is higher (corresponding to a decrease of the coarse-graining variance ).

Figure 3

The Q function of the reduced density matrix of the resonator ultra-strongly-coupled to the qubit after the coarse graining-measurement. We consider in (a) a projection into at a time for , and (b) the same projection for . These examples confirm that, as we increase the value of , the change of the function induced by the measurement becomes smaller. We set , , , , , , , , and .

In addition, in Fig. 4 we plot the post-measurement observable for the state of the qubit, for the different post-measurement outcomes and .

Figure 4

after the coarse-graining measurements that projects the state into or depending on the measurement results. We plot for the nonlinear resonator (a) and for the linear resonator (b) as the measurement apparatus. Here, we set the coarse-graining value as . For the other parameters, we use the same as those in Fig. 3.

For comparison, we consider both a linear resonator () and a nonlinear resonator () as the measurement devices. In addition, in Fig. 5, we show the average photon number inside the measurement resonator, for the case of a nonlinear and a linear measurement device. We can see that the number of photons at time ns is almost the same. In all figures, for the nonlinear measurement resonator, when we set the coarse-graining value as , the post-measurement states of the USC cavity and the qubit change significantly, depending on the measurement outcome, and so we observe a clear measurement backaction on the ultra-strongly coupled system.

Figure 5

The average number of photons in the nonlinear resonator (dashed green curve) and in the linear resonator (blue continous curve). The parameters used are the same as those in Fig. 4.

The state of the qubit (Fig. 4) depends strongly on the outcome of the coarse-grained measurement in the case of the nonlinear resonator, indicating a clear measurement backaction on the ultra-strongly coupled system. On the other hand, we cannot observe the measurement backaction in the case of a linear resonator. Interestingly, although we set the coupling strength as approximately 300 times smaller than the qubit energy in these examples, the nonlinear resonator can still affect the state of the system on reasonable time scales.

An ideal result would return perfect post-selection correlations, such that and , indicating the fidelity of our result is imperfect. However, a naive application of the rotating-wave approximation to the system and measurement device coupling term suggests that the influence of system and measurement apparatus on each other should be entirely negligible, and that there should be no correlation at all.

It is worth mentioning that such an approximation should also take into account the norm of the operator in the interaction term, which, for the driven nonlinear resonator, can be large. In our case, the number of photons in the nonlinear resonator is around 20. Surprisingly, even if we consider the norm of the operator in the interaction term (which corresponds to the number of the photons in the nonlinear resonator), the effective coupling strength ( multiplied by the norm) is still approximately 15 times smaller than the qubit energy in these non-energy eigenbasis measurements. So we can conclude that a clear difference between the case with a finite and the case without in our numerical simulations cannot be explained by a simple application of the rotating wave approximation.

In addition, a comparison to a linear measurement device for the same regime shows almost no correlation at all, so while the nonlinear device cannot perfectly project onto eigenstates of , it is capable of performing partial projection. In the next sections we will show that a better result can of course be reached by decreasing or increasing , and we will explain why the nonlinear resonator works, albeit imperfectly, in the difficult parameter regime when is much smaller than .

Comparison to the QND limit

To further compare our non-energy eigenbasis measurements with an ideal quantum non-demolition (QND) measurement, it is instructive to look at the behavior of the function of the nonlinear resonator, as shown in Fig. 6. Here, we consider our non-energy eigenbasis measurement, quantum non-demoliton measurements for the limit (which makes the measurement satisfy the QND condition ), and null measurements with .

Figure 6

The Q functions of the nonlinear resonator for several conditions: (a) Numerical simulation of the full Hamiltonian described in Eq. (1). (b) The two-level system approximation described by the Hamiltonian in Eq. (3). (c) Ideal QND measurement, which is possible in the limit . (d) When the nonlinear resonator does not couple at all with the qubit. We use the same parameters as those in Fig. 4.

We observe a clear difference between our non-energy eigenbasis measurements and measurements in the QND limit. In particular, the probability to obtain the high-amplitude state in the nonlinear resonator becomes much larger for QND measurements than that for the non-energy eigenbasis measurement case.

Also, the figures show that, roughly speaking, the probability to obtain the high-amplitude state of the resonator for the non-energy eigenbasis measurements lies between the case of the QND measurements and null measurements.

Low-energy two-level approximation

To give an intuitive explanation for why the nonlinear resonator measurement apparatus can become strongly correlated with the USC system, even when the coupling between measurement apparatus and system is much smaller than the system energy scales, we introduce a two-level approximation for the USC system. (see Appendix D for details, and a careful analysis of the validity of this approximation). In our simulations, the initial state is , and the interaction Hamiltonian mainly induces a transition from to the first excited state . Since the transition matrix elements of the interaction Hamiltonian to the other excited states are negligible, we can approximate the low-energy states of the ultra-strongly-coupled system as a two-level system. In this case, and can be written aswhereand

In Fig. 7, we plot , corresponding to and , with this two-level system approximation. To check the validity of this simplified model, we plot with this model and using the full model in Fig. 7. These results show an excellent agreement. Also, we confirm that the behavior of the function of the nonlinear resonator for the two-level approximation (b) agrees well with the full Hamiltonian case (a), as shown in Fig. 6.

Figure 7

Numerical simulations of the expected values of and after the state is projected into or depending on the measurement results. We use the same parameters as those in Fig. 4.

The role of the AC Stark shift to induce correlations

With this two-level approximation, we can show that the large correlation between the nonlinear resonator and the ultra-strongly-coupled system originates from the combination of an AC Stark shift and an adiabatic transition. It is easy to see that the large number of photons in the nonlinear resonator induces an energy shift (AC Stark shift) of the USC two-level system. Since the photon number of the high-amplitude state is different from that of the low-amplitude state, the size of the AC Stark shift strongly depends on the state of the nonlinear resonator. As long as the timescale of the change in the nonlinear resonator photons is much smaller than , the state of the two-level system remains in the ground state of the following effective Hamiltonianwhere () is the average photon number of the high (low) amplitude state.

When the nonlinear measurement resonator becomes a mixed state of the low- and high-amplitude states, we expect that the AC Stark shift (whose amplitude depends on the nonlinear resonator state) induces an adiabatic change of the ground state of the two-level system. This leads to a large correlation between the USC system and the measurement resonator. To show the validity of this interpretation, we analytically calculate the of the ground state of the Hamiltonian in Eq. (3) where we substitute the numerically calculated photon numbers of the high (low) amplitude state for (). In Fig. 8, we compare these results with the numerical simulations[65,66] where the master equation with the simplified Hamiltonian is solved. We plot the result from to , because from to , the high amplitude state is not generated. There is a good agreement between these two results, leading us to conclude that the correlation between the two-level system and the nonlinear resonator is induced by the aforementioned adiabatic changes due to the AC Stark shift, whose amplitude depends on the nonlinear resonator state. For a more detailed explanation of the AC Stark shift, please refer to Appendix F. Also, to quantify such a correlation, in Appendix G, we discuss the evolution of the quantum discord, which also shows a strong correlation even when the coupling strength is much smaller than the system energy scale.

Figure 8

Numerical results and analytical solutions of the expected values of and after the nonlinear resonator is projected into a high-amplitude state or a low-amplitude state. In the analytical calculations, we use the simplified Hamiltonian described in Eq. (3). Note that we do not show the time evolution from to , because the high-amplitude state is not generated until approximately .

Moreover, we increase the ratio to check how the effect of the AC Stark shift will change. In Fig. 9(a), we plot , , and the function at where the effective energy is 10% of what was used in Fig. 4. Interestingly, although we increase , the backaction of the measurement shown in Fig. 9(a) becomes smaller than that shown in Fig. 4, which is also attributed to the combination of an AC Stark shift and an adiabatic transition. From Fig. 9(a), the system converges into an eigenstate of after the interaction, regardless of the measurement results of the nonlinear resonator. This can be understood by considering that the AC Stark effect becomes much larger than the effective energy , so that the state of the ultra-strongly-coupled system becomes an eigenstate of for both the high amplitude state and low amplitude state. Furthermore, it is worth mentioning that, from Fig. 9(b), the nonlinear resonator before the measurement almost becomes a high-amplitude state. For an ideal quantum projective measurements on the ground state of the ultra-strongly coupled system, the population in the low-amplitude state should be the same as that of the high-amplitude state, and so this result shows that the effective energy is still too large to realize a full projective measurement in the persistent current basis.

Figure 9

(a,b) Dynamics of the ultra-strongly-coupled system and the nonlinear resonator when the effective energy is 10 times smaller than those in Fig. 4. (a) The expected value of after the nonlinear resonator is projected into a high-amplitude state or a low-amplitude state. (b) The Q function of the nonlinear resonator at time . (c,d) Dynamics of the ultra-strongly-coupled system and the nonlinear resonator when the effective energy is 100 times smaller than those in Fig. 4. (c) The expected value of after the nonlinear resonator is projected into a high-amplitude state or a low-amplitude state. (d) The function of the nonlinear resonator at time . Except for the effective energy of the ultra-strongly coupled system, we use the same parameters as those in Fig. 4.

We also consider a case when the effective energy is 1% of that used in Fig. 4. In that case, becomes much larger than , and this cannot be explained just by the AC Stark shift. Moreover, from Fig. 9(d), the population of the high-amplitude state becomes comparable with that of the low-amplitude state. Therefore, in this regime, a strong projection of the ground state of the ultra-strongly-coupled system in the non-energy eigenbasis seems to be realized, which can be quantified by calculating the entanglement between the system and measurement apparatus, as shown in the next section.

Negativity

As a criteria of entanglement, and to understand how correlations between a nonlinear resonator and a USC system develop, we now consider the negativity. Suppose there is a Hilbert space of two systems, with a state . The definition of negativity ishere, is the partial transpose of the state taken over a subsystem A, and is the trace norm[67]. In our case, the subsystem A corresponds to the two-level system approximation of the USC system, and B to the nonlinear resonator.

In Fig. 10 we plot the negativity to quantify the entanglement between the ultra-strongly-coupled system and the nonlinear resonator. As we increase the ratio , the negativity also increases. These results show that a reasonably large entanglement between the ultra-strongly-coupled system and the nonlinear resonator is generated in the regime where we realize a projective measurement on the non-energy eigenbasis. However, due to the decoherence of the nonlinear resonator, the entanglement quickly degrades, and a classical correlation remains in these systems just before the measurement on the nonlinear resonator.

Figure 10

Entanglement between the ultra-strongly coupled system and the nonlinear resonator. We use the negativity as a measure of entanglement. From the top, we plot results with effective energies , , and . Except for the effective energy of the ultra-strongly coupled system, here we use the same parameters as those in Fig. 4.

Conclusions

In conclusion, we investigated non-energy-eigenbasis measurements using a nonlinear resonator measurement apparatus, for the purpose of manipulating an ultra-strong-coupling light-matter system. Here we considered a specific example of a circuit QED system, but the results give us intuitive insights into how non-energy-eigenbasis measurements work in more general scenarios. Interestingly, we found that, even if the effective coupling strength with the measurement device is one order of magnitude smaller than the typical energy scale of the ultra-strongly-coupled system, we can still observe a strong correlation with the measurement device. While this correlation is imperfect, it appears in a difficult parameter regime where traditional linear measurement devices fail completely, and thus maybe practically useful for partial measurement and control of USC systems.

Also, we confirmed that, as one intuitively expects, by increasing the coupling strength with the measurement device, strong entanglement between the system and measurement device can be generated, and we can realize strong projective measurements on the ground state of the ultra-strongly-coupled system.

Nonlinear measurement devices are attractive for their fast and strong response to the system, but they are also difficult to analyze theoretically because of their nonlinear nature. Our results help illuminate the mechanism of how an ultra-strongly coupled system interacts with a nonlinear measurement device in a non-trivial parameter regime, where QND condition is not satisfied

Appendix

A. Non-energy-eigenbasis measurements

Here we explain the reason why the non-energy eigenbasis measurement is difficult to realize. Naive calculations indicate that the non-energy eigenbasis measurements would require a violation of the rotating wave approximation, which needs a strong coupling between the system and apparatus. This seems to suggest that, unless the coupling between the system and measurement apparatus is as large as the resonant frequency of the system and measurement apparatus, it would be difficult to implement the non-energy basis measurements. However, our results show that this naive picture is actually wrong if we use the nonlinear resonator as a measurement apparatus.

We can explain these points more quantitatively as follows. Suppose the Hamiltonian which expresses the coupling between a qubit and a linear resonator is as follows,

In the rotating frame defined by a unitary operator  , we obtain

In the limit of a large , we can use a rotating wave approximation and we obtainin the rotating frame.

More generally, we have a Hamiltonianwhere and , (the superindex (S) denotes the system and the superindex (E) denotes the measurement apparatus). In a rotating frame defined by  , we havewhere . If the system and measurement apparatus are well detuned, we obtainwhere we used the rotating wave approximation. So the terms that commute with survive. This clearly shows that we can measure only an observable that commutes with if the rotating wave approximation is valid. This also means that we need a violation of the rotating wave approximation for the non-energy eigenbasis measurements.

B. Derivation of the interaction Hamiltonian between the nonlinear resonator and the qubit

In this work we rely on an interaction between a superconducting flux qubit coupled with a frequency tunable resonator. This is not a dispersive approximation to a dipolar coupling. In more detail, the flux qubit is described aswhere denotes the energy bias and Δ denotes the tunneling energy. The Pauli matrix denotes the population of a persistent current basis such as , where () denotes a left-sided (right-handed) persistent current.

The frequency tunable resonator is described aswhere denotes the frequency of the resonator. We assume that the resonator contains a SQUID structure, and we can tune the frequency of the resonator by changing an applied flux penetrating the SQUID structure. (For example, see ref. [68]).

We can derive the interaction between the flux qubit and the resonator as follows. The persistent current states of the flux qubit induce magnetic fields due to the Biot-Savart law, and this changes the penetrating magnetic flux of the SQUID in the resonator. So the frequency of the resonator depends on the state of the flux qubit. Suppose that (−) denotes the magnetic flux from the () state, and the resonator frequency will be approximately shifted by . This provides us with the following Hamiltonian.where . A similar Hamiltonian has been derived in[69] to represent the coupling between an NV center and a flux qubit. We now assume a large detuning between the flux qubit and resonator. In this case, the dipolar coupling is negligible.

C. Derivation of the coarse-graining measurement

In the case where there is noise in the measurement apparatus, when we have a position measurement, even if the result of the measurement apparatus is , the real value is not necessarily . To model such situations, we define a measurement operator as followswhere describes the strength of the noise. Also, satisfies the normalization condition

Here we consider a composite system which comprises of a system which we hope to readout (ultra-strongly coupled system) and its probe (nonlinear resonator). Also, the measurement result is divided into and . When we have a measurement on a composite system , the post-measurement state when the result is becomeswhere

By tracing out the probe system, we have the post-measurement state of the system we hope to readout as

Substituting , can be rewritten aswhere is a complementary error function, and is defined as

In the limit of , we havewhich is a noiseless measurement. Also, in the limit of , we obtainwhich shows that we cannot extract any information from the system.

D. Validity of the two-level approximation

D.1 Adiabatic approximation to the Rabi Hamiltonian

We now explain the adiabatic approximation to the Rabi Hamiltonian, that has also been used in previous works[17,19-21]. We will show that, within the framework of the adiabatic approximation, the ultra-strongly coupled system can be treated as a two-level system. The conventional Rabi Hamiltonian can be written as

The adiabatic approximation can be done when and the Rabi Hamiltonian can be diagonalized using the baseswhere , and are eigenstates of , is the eigenstates of , and is a displacement operator. The states and are degenerate in energy and their energy is . Then, considering that the term couples these terms, and only the transitions between the states of the same are taken into account in the adiabatic approximation, the Rabi hamiltonian can be rewritten aswherewhose eigenvalues are

Also, it can be easily shown that

So, as long as we apply the adiabatic approximation, the transition due to the term is between and . Since the interaction between the ultra-strongly coupled system and the nonlinear resonator can be expressed as , it is possible for us to consider that the ultra-strongly coupled system is driven only by the operator. Also, if the initial state is and the perturbation term is proportional only to , the dynamics is limited to . Therefore, as long as the adiabatic approximation is valid, we can consider our subsystem of the ultra-strongly coupled system as a two-level system.

D.2 Estimation of the deviation from the two-level approximation

By calculating the deviation from the two-level approximation, we show a quantitative analysis of how accurate the two-level system approximation is in our parameter regime. We consider a fidelity between the true ground state (the first excited state ) and (). It is possible to estimate the accuracy of our two-level approximation from this fidelity, and we derive a condition of the fidelity to be close to unity. Now, we defineand

Here, is the one defined in Eq. (4). In this way, we regard as the non-perturbative Hamiltonian and the perturbative Hamiltonian. By performing a perturbative calculation up to the lowest order, we obtainwhere is a normalization factor. Then by using perturbation theory, we have

It can be easily shown thatand

Then, we have

By assuming , we haveand, we obtainwhere

Also, we set .

Similarly, with regard to the first excited state, we can obtain

(Note that in Eqs. (5) and (6) are the same.) The fidelity and are calculated as

Then, we define

For , we have , and so we can consider f as an infidelity.

We plot  f for three regimes . Here, we fix . From Fig. 11, we can see that in these regimes the infidelity f is sufficiently small.

Figure 11

The infidelity versus qubit frequency for three regimes (red solid), (blue dashed), (green dashed). Here varies from to , where .

Also, we plot the numerically calculated

(j = 2, 3, 4) in Fig. 12 in the same regime, where , , are the second, third and fourth excited states, respectively. This shows the leakage from to unwanted states.

Figure 12

A measure of the leakage to excited states versus qubit frequency . Again varies from to . Here we set (a) , (b) , where .

From Figs. 11 and 12, we confirm that, as long as is satisfied, transitions from to the unwanted states such as , , are small, so that the two-level approximation should be valid in this regime.

E. Losses in the ultra-strongly coupled system

So far we did not include a full analysis of the USC losses because we assumed that the time scale of such losses would be much longer than the readout time. For completeness, here we present a brief analysis of the influence of such losses. Because including bath-induced transitions between all eigenstates in the full space is complex, here we restrict ourselves to the two-level approximation. We justify this approximation, in our relevant parameter regime, in the previous sections.

The interaction Hamiltonian between our ultra-strongly coupled system and its environment is described bywhere denotes the position operator and () denotes the environmental operator coupled with the resonator (qubit). Also, we incorporate the effect of a dephasing bath classically modeled aswhere is a time-dependent random variable and the ensemble average of is zero. In this case, it is well know that the Born-Markov secular Lindblad master equation can be written in the form[20]where is the system Hamiltonian and

Here, and are the eigenstates of the system Hamiltonian and and are the rates corresponding to the noise spectra of the qubit and resonator, respectively. Also,andwhere denotes the spectral density of the qubit dephasing at frequency . Here, we ignore the term , as this term is negligible when we operate at the “sweet spot” of the qubit. Owing to the two-level approximation, we consider only the lowest first two levels and , and defining , and , we obtainwhere . In Fig. 13(a), we plot the expectation value of and without the noise in the ultra-strongly coupled system. The two-level approximation shows an excellent agreement with the full Hamiltonian model. Also, in Fig. 13(b), we plot and including the noise in the ultra-strongly-coupled system with parameters[70] that are realized in recent experiments[71]. From these results, we can conclude that the noise in the ultra-strongly coupled system is almost negligible and does not have significance on the time scales in which we are interested.

Figure 13

The expectation values (using the full Hamiltonian in Eq. (1)) and (using the approximate Hamiltonian in Eq. (3)) (a) without and (b) with noise in the ultra strong system, after the coarse-graining measurements that projects the state into or , depending on the measurement results. We set the coarse-graining value as . The noise rate . For the other parameters, we use the same as those in Fig. 2 in the main text.

F. Detailed explanation of the AC Stark shift

We assume the regime where the time scale of the dynamics of photons in the nonlinear resonator is much smaller than . In this case, we can approximate the interaction term , where is the expectation value of the photon number in the nonlinear resonator.

Therefore, the approximated two-level system can be described with the effective Hamiltonian,

Now, as the initial state is , and , the initial state is also the ground state of . Due to the assumption that slowly changes compared with the evolution induced by , the state stays at the ground state of .

Then, the nonlinear resonator can produce the mixed state due to the bifurcation aswhere and are the high and low amplitude states at time t. Now, the approximated two-level system evolves according to and , which is a AC Stark shift. More concretely, the approximated two-level system birfurcates into two states: the ground states of , where is the expectation value of the photon number corresponding to the high and low amplitude states. Consequently, the dynamics of the composite system of the approximated two-level system and the nonlinear resonator can be described aswhere is the ground state of for high and low amplitude states. Considering the case that the number of photons in the high-amplitude states is extremely high, with the low-amplitude state having approximately zero photon, , and .

Therefore, the post measured state corresponding to high and low amplitude states are and . Thus, the results in Fig. 4 became asymmetric. In addition, the dynamics of the birfurcation crucially depends on the pumping strength (also the detuning, and the strength of the nonlinearlity). In some regime, a bifurcation does not appear, and the dynamics described above completely changes.

G Quantum discord

To elucidate the previous results further, we consider the quantum discord (QD), which is defined as follows. Two possible definitions of the mutual information of the state where is the von Neumann entropy for a state , is the reduced density operator for , and is the quantum generalization of the conditional entropy. In the purely classical case, one can show that these two definitions of the mutual information are equivalent. However, in the nonclassical case, these definitions do not necessarily coincide. Also, is dependent on the measurement basis for . Therefore, the QD is defined aswhere

Here is a projector when the result is , and the QD is basis independent and reflects only nonclassical correlations[72,73]. In our case, system corresponds to the approximated two-level system and system the nonlinear resonator. We set the measurement basis on the approximated two-level system as ,, where and are the eigenstates of . Given these definitions we find the which realize .

We plot the QD in Fig. 14. Interestingly, in contrast to the negativity, the QD, at , becomes larger as is decreased. This can be explained in the following way: if is sufficiently large, the state becomes a highly entangled state well approximated by the formwhich decays, due the measurement of the nonlinear cavity, to the mixturewhere and are high and low amplitude states of the nonlinear resonator. Since and are orthogonal to each other, is a classically correlated state without any superposition, implying the vanishing QD. On the other hand, when is small, the dynamics can be explained by an AC Stark shift and the state can be expressed aswhere is the probability that the nonlinear resonator is in the low (or high) amplitude state. The state is the ground state of . Here, and are not always orthogonal to each other, and, as such, the correlation in the mixture of the two could have a non-classical nature. Hence, the QD, in the long-time limit, tends to have a finite value when J/ is small.

Figure 14

Quantum discord between the ultra-strongly coupled system and the nonlinear resonator. From the top, we plot results with effective energies , , , and . Except for the effective energy of the ultra-strongly coupled system, here we use the same parameters as those in Fig. 4.