- Research
- Open access
- Published:
Controlling the dynamic range of a Josephson parametric amplifier
EPJ Quantum Technology volume 1, Article number: 2 (2014)
Abstract
One of the central challenges in the development of parametric amplifiers is the control of the dynamic range relative to its gain and bandwidth, which typically limits quantum limited amplification to signals which contain only a few photons per inverse bandwidth. Here, we discuss the control of the dynamic range of Josephson parametric amplifiers by using Josephson junction arrays. We discuss gain, bandwidth, noise, and dynamic range properties of both a transmission line and a lumped element based parametric amplifier. Based on these investigations we derive useful design criteria, which may find broad application in the development of practical parametric amplifiers.
1 Introduction
Due to the rapidly evolving field of quantum optics and information processing with superconducting circuits the interest in low-noise amplifiers has dramatically increased in the past five years and has lead to a body of dedicated research on Josephson junction based amplifiers [1–15]. The most successful quantum limited detectors which have so far been realized in the microwave frequency range are based on the principle of parametric amplification [16–19]. Josephson parametric amplifiers (JPAs) have not only been used to generate squeezed radiation [6, 13, 20–22], but moreover enabled the realization of quantum feedback and post-selection based experiments [23–26], the efficient displacement measurement of nanomechanical oscillators [27] and the exploration of higher order photon field correlations [28, 29].
While JPAs have been demonstrated to operate close to the quantum limit, their performance is to date mostly limited by their relatively small dynamic range, i.e. the saturation of the gain for large input signals. Here, we discuss the control of the dynamic range by making use of Josephson junctions arrays in the parametric amplifier circuit, which we have already employed in recent experiments [29, 30]. After reviewing the principles of parametric amplification we discuss bandwidth and noise constraints in dependence on the circuit design, based on which we derive simple strategies for optimized circuit design.
2 Principles of parametric amplification
2.1 Parametric processes at microwave frequencies
In quantum optics the word parametric is used for processes in which a nonlinear refractive medium is employed for mixing different frequency components of light. Such processes are parametric in the sense that a coherent pump field, applied to a nonlinear medium, modulates its refractive index, which appears as a parameter in a semi-classical treatement. This time-varying parameter is affecting modes with frequencies detuned from the frequency of the pump field and can stimulate their population with photons. The energy for creating these photons is provided by the pump field.
The refractive index in optics is equivalent to the impedance of electrical circuits. In order to realize parametric processes at microwave frequencies we therefore modulate an effective impedance. This is achieved by varying the parameters of either a capacitive or an inductive element in time. Although there have been early proposals for fast time-varying capacitances [31], it now is considered to be more convenient to make use of dissipationless Josephson junctions for this purpose. In a regime in which the current I flowing through a Josephson junction is much smaller than its critical current its associated inductance is approximately . Applying an AC current through the junction using appropriate microwave drive fields therefore leads to the desired time-varying impedance. Because of the proportionality of the inductance L to the square of the current , such a drive results in a four-wave mixing process [32].
The effective impedance can alternatively be modulated by varying the magnetic flux threading a superconducting quantum interference device (SQUID) loop [5] such that the effective inductance is approximately modulated proportionally to the AC current flowing in the loop, . The quantity in this expression depends on the DC flux bias point of the SQUID loop. Since the relation between current and inductance is in this case linear, the magnetic flux drive results in a three-wave mixing process [33].
In order to enhance parametric amplification in a well-controlled frequency band while suppressing it for frequencies out of this band, the modulated Josephson inductance is frequently integrated into a microwave frequency resonator. This is the simplest way to control the band in which parametric amplification occurs. A number of variations of this basic idea are now explored. The circuit design has recently been modified to achieve a spatial separation of signal and idler modes [9, 21, 34–36] and to build traveling wave amplifiers, in which a field is amplified while propagating in forward direction coaxially with a pump field [37, 38]. Various drive mechanisms ranging from single and double pumps [8] to magnetic flux drives [5, 39, 40] have been explored. Being aware of this variety of possible approaches, we focus here on a single mode (degenerate) parametric amplifier driven with one pump tone close to its resonance frequency.
2.2 Circuit QED implementation of a parametric amplifier
The JPA essentially is a weakly nonlinear oscillator, in which the nonlinearity is provided by Josephon tunnel junctions. In practice, this is typically realized either as a transmission line resonator shunted by a SQUID [7, 12, 13], see Figure 1(a), or as a lumped element nonlinear oscillator [11]. The use of a SQUID instead of single tunnel junction guarantees tunability of the resonance frequency. Since resonator-based parametric amplifiers provide amplification in a narrow band only, tunability is highly desirable to match the band of amplification with the frequency of the signal to be amplified.
The relevant part of the Hamiltonian which describes the parametric amplifier considered here can be written as
where A labels the annihilation operator of the intra-resonator field. Expressions for the resonance frequency and the effective Kerr nonlinearity K are derived in Section 4 based on the full circuit model. In the following section we analytically study the dynamics of this system using the input-output formalism. Before presenting the mathematical derivations, we qualitatively describe different dynamical regimes of this nonlinear oscillator and explain the mechanism which leads to amplification.
If we assume for the moment that the JPA has no internal losses, all the incident power is reflected from the resonator and the classical response (i.e. reflection coefficient) is completely specified by the phase φ of the reflected field. In contrast to a linear system, where φ only depends on the frequency , it also depends on the power of the probe field in the case of a nonlinear oscillator. In Figure 1(b), the theoretically expected value of φ is plotted as a function of the probe amplitude for two characteristic drive frequencies. While the phase is constant for low drive powers (quasi-linear response), the phase changes significantly for increased drive power. Depending on the probe frequency we either find a bistable regime where two stable solutions exist [41, 42] or a regime where the phase has a unique solution (red and blue data sets in Figure 1(b)). In both cases the phase significantly depends on the input power. The bistable response can for example be used to realize a bifurcation amplifier [43, 44] and for nonlinear dispersive readout [45], which has been intensely studied in the context of circuit QED.
Since we are particularly interested in linear amplification the following discussion is focused on the regime, in which the response has a unique solution (blue data set). The mechanism of amplification can be understood qualitatively in the following way. If we imagine that the device is constantly driven at a frequency and power at which the reflected phase φ depends sensitively on power (see gray circle in Figure 1(b)), the system will strongly react to small perturbations. Such perturbations, which could be caused by an additional small signal field for example, are therefore translated into a large change of the output field.
We illustrate this process leading to amplification by plotting the resonator response for input fields with slightly varying amplitude and phase. In Figure 1(c) we indicate the input fields by a blue circle around the mean value (arrow). The small differences in amplitude of the input field translate into large changes in φ of the output field (red ellipse). If we interpret the arrow in Figure 1(c) as a constant pump field and its difference to the points on the blue circle as an additional signal, the signal is either amplified or deamplified depending on its phase relative to the pump.
The mechanism of amplification can thus be understood intuitively by considering the nonlinear response to a monochromatic drive field. In order to characterize the exact behavior of input fields with finite bandwidth we analyze the response in more detail below.
3 Input-output relations for the parametric amplifier
3.1 Classical nonlinear response
Here, we employ the input-output formalism [46, 47] to calculate the nonlinear resonator response discussed qualitatively in the previous section. The derivation presented here is inspired by Ref. [1]. A schematic of the input-output model is shown in Figure 2. The nonlinear resonator is coupled with rate κ to a transmission line, through which the pump and signal fields propagate. Based on this model and the Hamiltonian in Eq. (1) we obtain the following equation of motion for the intra-resonator field
In addition to the coupling to transmission line modes with rate κ we account for potential radiation loss mechanisms by introducing the coupling to modes with loss rate γ, compare Figure 2(a). A boundary condition equivalent to
also holds for the loss modes. When operating the device as a parametric amplifier, the input field is typically a sum of a strong coherent pump field and an additional weak signal field. Since this signal carries at least the vacuum noise, it is treated as a quantum field. In this formalism this particular situation is accounted for by decomposing each field mode into a sum of a classical part and a quantum part
where α, , represent the classical parts of the field which are associated with the pump, while a, , account for the quantum signal fields. Since all α’s are complex numbers the modes a satisfy the same bosonic commutation relations as modes A do. By multiplying the field modes defined in Eq. (4) with the additional exponential factor , one works in a frame rotating at the pump frequency . The strategy is to first solve the classical response for the pump field α exactly and then linearize the equation of motion for the weak quantum field a in the presence of the pump. Finally, we derive a scattering relation between input modes and reflected modes .
The steady state solution for the coherent pump field is determined by
which follows immediately by substituting Eq. (4) into Eq. (2) and collecting only the c-number terms. By multiplying both sides with their complex conjugate we get to the equation
which determines the average number of pump photons in the resonator. Eq. (6) reduces to
by defining the scale invariant quantities
δ is the detuning between pump and resonator frequency in units of the total resonator linewidth, is the dimensionless drive amplitude, and ξ is the product of drive power and nonlinearity, also expressed in dimensionless units. Finally, n is the mean number of pump photons in the resonator relative to the incident pump power. As an important consequence, we notice from Eq. (8) that only the product of drive power and nonlinearity determines the dynamics but not each quantity itself. Therefore, a small nonlinearity can at least in principle be compensated by increasing the drive power. Properties such as the gain-bandwidth product are therefore independent of the strength of the nonlinearity as long as the pump power is much larger than the power of amplified fluctuations. Furthermore, the solutions of Eq. (7) for negative ξ values are identical to those for positive ξ up to a sign change in δ. Since ξ is negative for the Josephson parametric amplifier, we focus on this particular case.
Equation (7) is a cubic equation in n and can therefore be solved analytically. We do not present the lengthy solutions here explicitly, but assume in the following that we have an explicit analytical expression for n in terms of δ and ξ. In Figure 2(b) we plot n for various parameters ξ as a function of δ. At the critical value the derivative diverges and thus the response of the parametric amplifier becomes extremely sensitive to small changes. For even stronger effective drive powers the cubic Eq. (7) has three real solutions. The solutions for the high and low photon numbers are stable, while the intermediate one is unstable. The system bifurcates in this regime as mentioned earlier. The critical detuning below which the system becomes bistable is . The critical point is the one at which both and vanish. In scale invariant units the maximal value of n is 4, which is reached at the detuning .
Experimentally, the system parameters are characterized by measuring the complex reflection coefficient . Based on the input-output relation and Eq. (5) we evaluate this reflection coefficient as
In Figure 2(c) we plot the absolute value of the reflection coefficient at for various loss rates γ. For vanishing losses all the incident drive power is reflected from the device and . Note that also in this case the resonance is clearly visible in the phase of the reflected signal (not shown here). When the loss rate γ becomes similar to the external coupling rate κ part of the radiation is dissipated into the loss modes. In the case of critical coupling all the coherent power is transmitted into the loss modes at resonance. This is equivalent to the case of a symmetrically coupled resonator, for which the transmission coefficient is one at resonance [48].
3.2 Linearized response for weak (quantum) signal fields
Under the assumption that the photon flux associated with the signal is much smaller than the photon flux of the pump field , we can drop terms such as , because they are small compared to the leading terms and . By neglecting these terms we obtain a linearized equation of motion for a in the presence of the pump field. In order to preserve the validity of this approximation even for larger input signals, the amplitude α of the pump field needs to be increased. Experimentally, this can be achieved by reducing the strength of the nonlinearity K as discussed in Section 3.3 in more detail. Substituting Eq. (4) into Eq. (2) and keeping only terms which are linear in a one finds
Since Eq. (10) is linear, we can solve it by decomposing all modes into their Fourier components
and equivalently for and . Note that the detuning Δ between signal frequencies and the pump frequency, is expressed here in units of the linewidth . Substituting the Fourier decompositions into Eq. (10) and comparing the coefficients of different harmonics, results in
where is the sum of all field modes incident on the resonator. Furthermore, in Eq. (12) ϕ is the phase of the intra-resonator pump field, defined by . The fact that Eq. (12) couples modes and can be interpreted as a wave mixing process. In order to express in terms of the input fields , Eq. (12) is rewritten as a matrix equation
By inverting the matrix on the right hand side, the quantum part of the intra-resonator field is expressed in terms of the incoming field
with . Using Eq. (3), the final transformation between input and output modes is
with
and
Eq. (15b) is the central result of this calculation. The output field at detuning Δ from the pump frequency is a sum of the input fields at frequencies Δ and −Δ multiplied with the signal gain factor and the idler gain factor , respectively. The additional noise contributions introduced via the loss modes vanish in the limit . In the ideal case , the coefficients and satisfy the relation
and Eq. (15b) is identical to a two-mode squeezing transformation [19, 49] with gain . The two-mode squeezing transformation describes a linear amplifier in its minimal form (compare Ref. [50]), of which we discuss characteristic properties in the following section.
3.3 Gain, bandwidth, noise and dynamic range
For simplicity we consider the case of no losses , for which the parametric amplifier response is described by Eq. (15b). An incoming signal at detuning Δ is thus amplified by the power gain and mixed with the frequency components at the opposite detuning from the pump. Characteristic properties of the parametric amplifier, such as the maximal gain and the bandwidth, are thus encoded in the quantity as a function of pump-resonator detuning δ, effective drive strength ξ and detuning between signal and pump Δ.
In Figure 3(a) we plot the gain for zero signal detuning as a function of δ and ξ. We find that the maximal gain increases with increasing drive strength ξ while the optimal value for δ at which this gain is reached, shifts approximately linearly with increasing ξ. The optimal values for δ are indicated as a dashed white line in Figure 3(a). Mathematically, the gain diverges when ξ approaches the critical value . In practice, the gain is limited to finite values due to the breakdown of the stiff pump approximation (see discussion below).
By changing the pump parameters ξ and δ we can adjust the gain to a desirable value, which is typically about 20 dB. Note that the gain can take values smaller than one, in the presence of finite internal losses . Once the pump parameters are fixed we characterize the bandwidth of the amplifier by analyzing the gain as a function of the signal detuning Δ. In Figure 3(b) we plot the gain as a function of Δ for the indicated values of and the corresponding optimal pump detunings δ (compare dashed white line in (a)). The gain curves are well approximated by Lorentzian lines as indicated by the dashed black lines in Figure 3(b). When the gain is increased, the band of amplification becomes narrower. This is quantitatively expressed by the gain-bandwidth relation , where B is the detuning Δ for which the gain reaches half of its maximal value. This gain-bandwidth relation follows from the Lorentzian approximation of the gain curves shown as black dashed lines in Figure 3(b)) and holds for gain values above a few dB. Remember that Δ is defined in units of the resonator linewidth , which means that the amplifier bandwidth equals approximately the resonator linewidth divided by the square root of the gain.
When operating the JPA, we also have to understand its behavior in terms of added noise. In the ideal case with zero loss rate (), the input-output relation of the parametric amplifier in Eq. (15b) has the minimal form of a scattering mode amplifier [19]. The amplification process reaches the vacuum limit as long as the input modes are cooled into the vacuum. In practice, however, the device may have finite loss γ which increases the effectively added noise by a factor of . This is due to the additional amplified noise, which originates from the modes and contributes to the output field (compare Eq. (15a)). Another potential source of noise is related to the stability of the resonance frequency of the parametric amplifier. Magnetic flux noise in the SQUID loop may lead to a fluctuating resonance frequency and thus a fluctuating effective gain.
In the derivation made in the previous sections we have assumed that the solution of the classical drive field is unaffected by the presence of additional signal and quantum fluctuations at the input. This is known as the stiff pump approximation [8], which assumes that the pump power at the output is equal to the pump power at the input of the JPA. This is of course an approximation, since the pump field provides the energy which is necessary for amplifying the input signal. The stiff pump approximation is valid as long as the pump power is significantly larger than the total output power of all amplified (quantum) signal and vacuum fields. In order to quantitatively analyze the pump depletion due to the presence of amplified fields we add the terms and to the left hand side of Eq. (5) and solve Eq. (5) and Eq. (10) self-consistently. This mean-field approach is similar to the one used in Ref. [8]. For our calculation we model the incoming signal field, which is to be amplified, as white noise with average photon number per unit time and bandwidth. Based on this model we find that the gain decreases when the signal strength exceeds a certain value, see Figure 4(a). The number of input photons at which this happens becomes smaller with decreasing ratio . This is expected because the pump power close to the bifurcation point is proportional to and provides the energy required for amplification. For small values the gain is reduced even for due to the amplification of vacuum fluctuations (see blue data points in Figure 4(a)). As a measure of the dynamic range we specify the 1 dB compression point of the JPA, i.e. the value of input photons at which the JPA gain decreases by 1 dB compared to the stiff pump approximated gain value. As shown in Figure 4(b), the 1 dB compression point increases proportionally to in the limit of . The presence of constant vacuum fluctuations leads to a gain compression by more than 1 dB even for when becomes smaller. We can qualitatively explain this behavior by comparing the pump power with the power of amplified signal and vacuum fields
Making use of the gain-bandwidth relation we find the following scaling of the ratio between the two powers
The results shown in Figure 4(b) together with Eq. (20) indicate that the validity of the stiff pump approximation is essentially determined by the ratio . The calculations furthermore show that the dynamic range can be increased by reducing the ratio of the JPA, which seems to be the case also for flux driven parametric amplifiers [40]. In Section 5 we discuss how to achieve small nonlinearities by making use of multiple SQUIDs connected in series.
4 Effective system parameters from distributed circuit model
In the previous section we have analyzed the model of a nonlinear resonator with resonance frequency , Kerr nonlinearity K and decay rate κ. Here, we explicitly derive this effective Hamiltonian from the full circuit model of a - transmission line resonator, which is terminated by a SQUID loop at the short-circuited end and coupled capacitively to a transmission line, see Figure 1(a). These calculations allow us to determine , K, κ from the distributed circuit parameters and give insight into potential limitations of the effective model. We also compare the obtained parameter relations with those of a lumped element parametric amplifier.
4.1 Resonator mode structure in the linear regime
In order to find the normal mode structure of the system, we first neglect its capacitive coupling to the transmission line as indicated in Figure 5. The derivation is similar to the ones in Refs. [40, 51, 52]. Dissipation effects due the environment are discussed in Section 4.3.
The total Lagrangian of the system in the magnetic flux field has a transmission line part and a term which describes the SQUID at position (Figure 5).
with the reduced flux quantum . Since we work in a limit in which the plasma frequency of the SQUID is much larger than the resonance frequencies of interest, we neglect the self-capacitance of the SQUID. Note that Ref. [52] provides a detailed study of the effect of the self-capacitance in various parameter regimes. We furthermore describe the SQUID as a single junction with tunable effective Josephson energy .
We first investigate the linear regime of the system, in which the cosine potential of the SQUID is approximated as a quadratic potential.
Due to the spatial derivative in the Lagrangian in Eq. (21) all local fields in the chain are coupled to their nearest neighbors and the normal mode structure is found by solving the Euler-Lagrange equation of the transmission line resonator. This results in the wave equation
with the phase velocity , of which the general solution can be written as a sum of normal modes
The valid wavevectors are determined by the boundary conditions at the two ends of the transmission line. The open end at requires that the current vanishes, which is implicitly satisfied by choosing the cosine ansatz in Eq. (24). On the shorted end the boundary condition is modified by the presence of the Josephson junction. In order to determine this boundary condition, we evaluate the Euler-Lagrange equation at position . For this purpose it is convenient to write the Lagrangian in a discretized form, see Figure 5(c) and compare Ref. [51]:
where and . Evaluating leads to the equation
Substituting the ansatz (24) into Eq. (26) and comparing the resulting coefficients of the independent variables , results in the transcendental equation
Here, we have defined the Josephson inductance . The infinite set of solutions of this equation determines the normal modes structure of the system in the linear regime. In the limit in which the SQUID inductance vanishes, Eq. (27) is solved by the poles of , and we recover the normal modes of the resonator
As a first order correction to this result in the limit of , we expand Eq. (27) to first order in and find or equivalently
For the fundamental mode with this linearized approximation is typically accurate even for inductance ratios up to , whereas for the higher harmonic modes the linearized equation breaks down for much smaller values of . A comparison between the exact solution based on Eq. (27) and the approximate solution in Eq. (29) is shown in Figure 6 for the first three resonant modes. When higher harmonics are expected to be relevant one should solve Eq. (27) numerically in order to determine the exact wave numbers .
4.2 Kerr nonlinear terms and effective Hamiltonian
Using the normal mode decomposition in Eq. (24) we reexpress the Lagrangian in Eq. (21) as a sum of oscillators which are only coupled via the boundary condition imposed by the SQUID. For the purposes of parametric amplification the phase drop across the junction is desired to be small, , i.e. the current flowing through the Josephson junction is small compared to its critical current. We can therefore expand the SQUID cosine potential and take into account only the first non-quadratic correction
In Section 5 we discuss under which circumstances such an approximation may break down. Substituting the normal mode decomposition Eq. (24) into the Taylor expansion of the Lagrangian results in
with the effective capacitances and inductances [51]
and the nonlinearity coefficients
As expected the linear part of the Lagrangian is diagonal in the normal mode basis. It describes a set of uncoupled LC oscillators for which the effective resonance frequencies coincide with the product of phase velocity and wave vector .
Based on the Lagrange function (31) we derive the Hamiltonian by introducing the conjugate charge variables . Performing a Legendre transformation and taking only self-interactions and two-mode interactions into account, results in the Hamiltonian
In a quantum regime and are operators which satisfy the commutation relation and it is convenient to write the Hamiltonian in terms of normal mode annihilation and creation operators [53]
with and . The abbreviation zpf stands for zero point fluctuations. Performing a rotating wave approximation (i.e. removing all terms with an unequal number of creation and annihilation operators), and neglecting the small photon number independent frequency shifts due to the nonlinear terms (i.e. Lamb shifts) we arrive at
with
The quantity is the Kerr nonlinearity of the fundamental mode, which is used for the parametric amplification process. The terms proportional to with unequal are cross Kerr interaction terms which couple different modes to each other. Such an interaction can for example be used for counting the number of photons in one mode by probing another one with a coherent field [54–61], similarly to a dispersive qubit measurement. Note that the values resulting from Eq. (37) are divided by the square of the number of SQUIDs, if an array is used instead of a single SQUID, as discussed in the following section.
4.3 Decay rate and resonance frequency correction for low Q resonators
Since the parametric amplifier bandwidth is proportional to the decay rate κ, typical devices are designed to have a low external quality factor, which is achieved by increasing the coupling capacitance between transmission line and resonator (Figure 1). The coupling of an oscillator to the environment shifts its resonance frequency [48], which can be significant if the coupling rate is large. When designing parametric amplifier devices, it is therefore necessary to take these shifts into account. Based on the effective inductance and capacitances calculated in Eq. (32) we find
for resonance frequency and decay rate of the j th mode of the parametric amplifier device. The external quality factor is given by .
4.4 Lumped element JPA
As already mentioned in the introduction, a JPA can also be realized as a lumped element resonator by shunting a SQUID with a large capacitance [10, 11]. In this case the resonator is described by the transmon Hamiltonian [62], which in the deep transmon limit takes the form of Eq. (1) with anharmonicity and resonance frequency . Also for this type of resonators the coupling rate κ to the transmission line can be designed independently of and by designing an appropriate capacitive network. Similarly as for the transmission line JPA, the description in terms of the effective Hamiltonian Eq. (1) is based on the assumption that for relevant resonator fields the phase drop across the Josephson junctions is small (compare Eq. (30)). In the following section we study the validity of this approximation when the resonator is driven close to the bifurcation point where we expect parametric amplification to occur and analyze its implications for realizing a parametric amplifier with large bandwidth and dynamic range.
5 Bandwidth and dynamic range constraints
5.1 Validity of the quartic approximation
In Section 3.3 we have shown that the dynamic range of the JPA scales with . Also the bandwidth becomes larger with increasing κ, which indicates that a large κ is desirable for JPAs. However, there are limitations on the maximal possible value for κ as discussed in the following.
For deriving the Hamiltonian in Eq. (1), or more generally Eq. (36), we have expanded the SQUID cosine potential to quartic order in the dimensionless flux variable , where is the phase drop across the SQUID. To guarantee that this approximation holds when we operate the device in the parametric amplification regime, we have to make sure that is small even when it is driven close to the bifurcation point. This is equivalent to keeping the current flowing through the SQUID small compared to the critical current.
To characterize the validity of the low order expansion of the cosine potential, we define the maximal coherent field inside the resonator as the one for which . This is the coherent amplitude, at which the current flowing through the SQUID equals its critical current. According to Eq. (35) and Eq. (24) a coherent field α in mode j leads to a maximal amplitude of across the tunnel junctions, based on which we define the critical amplitude as
The low order expansion of the cosine potential is only valid if the field inside the resonator α is much smaller than this maximal amplitude . In Section 3.1 we have found that the photon number in a resonator mode at the bifurcation point is . The ratio between and the maximal coherent photon number is therefore given by
In order to minimize the effect of higher order nonlinearities we want to keep this ratio small. It has the following two interesting limits
This is an important result, which sets clear constraints on both the maximally achievable bandwidth and the dynamic range of the JPA. If we want to be small, κ needs to be sufficiently small as well. One may also want to increase the dynamic range by reducing the Kerr nonlinearity . While this can in principle be achieved by choosing a small ratio between Josephson and geometric inductance, care has to be taken when using this approach because additional geometric inductance leads to a larger ratio as illustrated in Figure 7.
Interestingly, we find that in the lumped element case the Josephson inductance , and with it the Kerr nonlinearity K, can in principle be made smaller without affecting . However, in practice a small Josephson inductance has to be compensated by a large lumped element capacitor to retain the desired resonance frequency, which is challenging to realize without introducing additional parasitic geometric inductances. It therefore seems difficult to build a parametric amplifier with large bandwidth and high dynamic range at the same time using a single SQUID only. In the following we show how one can keep constant while decreasing the nonlinearity and thus increasing the dynamic range of the amplifier, by replacing the single SQUID with a serial array of M SQUIDs of M-times larger Josephson energy per SQUID (Figure 8).
5.2 Josephson junction arrays
For simplicity we assume that all SQUIDs in the array have the same effective Josephson energy . Since the spatial extent of the junction is still small compared to typical resonance wavelengths, we can treat the array as a lumped element. To derive the nonlinearity of the oscillator for this situation we investigate how the different terms in the Lagrangian scale with M.
Assuming that the phase drop from the flux node at the end of the transmission line resonator to the ground is homogeneously distributed over the array, we have the same phase drop across each SQUID, see Figure 8. As a result the quadratic term in the Lagrangian scales as
and thus remains constant. This agrees with our expectation, since the total linear Josephson inductance has not been changed. However, the quartic term scales like
which leads to a quadratic decrease in the effective Kerr nonlinearity and thus a quadratic increase in . Furthermore, the maximal photon number also scales as since the critical current of each junction is larger by a factor of M. In other words, the ratio only depends on the total Josephson inductance whereas the bifurcation power increases quadratically in M. We thus conclude that the dynamic range of a JPA can be increased without affecting the amplifier bandwidth, by using an array of SQUIDs instead of a single SQUID. This conclusion is valid for both the transmission line JPA and the lumped element JPA.
In practice, the Josephson energies in the array are not all equal due to inhomogeneous coupling to the external magnetic flux and scatter in the critical current of Josephson junctions due to unavoidable variations in fabrication. A quantitative analysis of the influence of such variations of Josephson energies on the parametric amplifier characteristics could be an interesting task for future studies. This would help to quantify limitations in the accessible tuning range of the parametric amplifier and a realistic understanding of the breakdown of the low order expansion of the cosine potential. For such an approach the methods used in Ref. [63] could turn out to be useful.
6 Conclusion
In summary, we have presented a detailed analysis of Josephson junction based parametric amplifiers, including a discussion of bandwidth, noise and dynamic range. By establishing relations between basic JPA properties and designable circuit parameters we have been able to derive two simple design strategies to achieve optimized JPA performance. On the one hand the contribution of the Josephson inductance to the total effective inductance of the resonant circuit has to be chosen sufficiently large. On the other hand the use of SQUID arrays instead of single SQUIDs provides the possibility to enhance the strength of the pump field close at the bifuraction point and with it the dynamic range of the JPA.
References
Yurke B, Buks E: J Lightwave Technol. 2006, 24: 5054.
Castellanos-Beltran MA, Lehnert KW: Appl Phys Lett. 2007., 91: Article ID 083509 Article ID 083509
Tholén E, Ergül A, Doherty E, Weber F, Grégis F, Haviland D: Appl Phys Lett. 2007., 90: Article ID 253509 Article ID 253509
Kinion D, Clarke J: Appl Phys Lett. 2008., 92: Article ID 172503 Article ID 172503
Yamamoto T, Inomata K, Watanabe M, Matsuba K, Miyazaki T, Oliver WD, Nakamura Y, Tsai JS: Appl Phys Lett. 2008., 93: Article ID 042510 Article ID 042510
Castellanos-Beltran MA, Irwin KD, Hilton GC, Vale LR, Lehnert KW: Nat Phys. 2008, 4: 929. 10.1038/nphys1090
Palacios-Laloy A, Nguyen F, Mallet F, Bertet P, Vion D, Esteve D: J Low Temp Phys. 2008, 151: 1034. 10.1007/s10909-008-9774-x
Kamal A, Marblestone A, Devoret M: Phys Rev B. 2009., 79: Article ID 184301 Article ID 184301
Bergeal N, Schackert F, Metcalfe M, Vijay R, Manucharyan VE, Frunzio L, Prober DE, Schoelkopf RJ, Girvin SM, Devoret MH: Nature. 2010, 465: 64. 10.1038/nature09035
Hatridge M, Vijay R, Slichter DH, Clarke J, Siddiqi I: Phys Rev B. 2011., 83: Article ID 134501 Article ID 134501
Vijay R, Slichter DH, Siddiqi I: Phys Rev Lett. 2011., 106: Article ID 110502 Article ID 110502
Gao J, Vale LR, Mates JAB, Schmidt DR, Hilton GC, Irwin KD, Mallet F, Castellanos-Beltran MA, Lehnert KW, Zmuidzinas J, Leduc HG: Appl Phys Lett. 2011., 98: Article ID 232508 Article ID 232508
Eichler C, Bozyigit D, Lang C, Baur M, Steffen L, Fink JM, Filipp S, Wallraff A: Phys Rev Lett. 2011., 107: Article ID 113601 Article ID 113601
Abdo B, Kamal A, Devoret M: Phys Rev B. 2013., 87: Article ID 014508 Article ID 014508
Kamal A, Clarke J, Devoret MH: Phys Rev B. 2012., 86: Article ID 144510 Article ID 144510
Louisell WH, Yariv A, Siegman AE: Phys Rev. 1961, 124: 1646. 10.1103/PhysRev.124.1646
Gordon JP, Louisell WH, Walker LR: Phys Rev. 1963, 129: 481. 10.1103/PhysRev.129.481
Mollow BR, Glauber RJ: Phys Rev. 1967, 160: 1076. 10.1103/PhysRev.160.1076
Clerk AA, Devoret MH, Girvin SM, Marquardt F, Schoelkopf RJ: Rev Mod Phys. 2010, 82: 1155. 10.1103/RevModPhys.82.1155
Mallet F, Castellanos-Beltran MA, Ku HS, Glancy S, Knill E, Irwin KD, Hilton GC, Vale LR, Lehnert KW: Phys Rev Lett. 2011., 106: Article ID 220502 Article ID 220502
Flurin E, Roch N, Mallet F, Devoret MH, Huard B: Phys Rev Lett. 2012., 109: Article ID 183901 Article ID 183901
Menzel EP, Di Candia R, Deppe F, Eder P, Zhong L, Ihmig M, Haeberlein M, Baust A, Hoffmann E, Ballester D, Inomata K, Yamamoto T, Nakamura Y, Solano E, Marx A, Gross R: Phys Rev Lett. 2012., 109: Article ID 250502 Article ID 250502
Vijay R, Macklin C, Slichter DH, Weber SJ, Murch KW, Naik R, Korotkov AN, Siddiqi I: Nature. 2012, 490: 77. 10.1038/nature11505
Johnson JE, Macklin C, Slichter DH, Vijay R, Weingarten EB, Clarke J, Siddiqi I: Phys Rev Lett. 2012., 109: Article ID 050506 Article ID 050506
Ristè D, van Leeuwen JG, Ku H-S, Lehnert KW, DiCarlo L: Phys Rev Lett. 2012., 109: Article ID 050507 Article ID 050507
Campagne-Ibarcq P, Flurin E, Roch N, Darson D, Morfin P, Mirrahimi M, Devoret MH, Mallet F, Huard B: Phys Rev X. 2013., 3: Article ID 021008 Article ID 021008
Teufel JD, Donner T, Li D, Harlow JW, Allman MS, Cicak K, Sirois AJ, Whittaker JD, Lehnert KW, Simmonds RW: Nature. 2011, 475: 359. 10.1038/nature10261
Eichler C, Bozyigit D, Wallraff A: Phys Rev A. 2012., 86: Article ID 032106 Article ID 032106
Eichler C, Lang C, Fink JM, Govenius J, Filipp S, Wallraff A: Phys Rev Lett. 2012., 109: Article ID 240501 Article ID 240501
Steffen L, Salathe Y, Oppliger M, Kurpiers P, Baur M, Lang C, Eichler C, Puebla-Hellmann G, Fedorov A, Wallraff A: Nature. 2013, 500: 319. 10.1038/nature12422
Louisell WH: Coupled Mode and Parametric Electronics. John Wiley, New York; 1960.
Slusher RE, Hollberg LW, Yurke B, Mertz JC, Valley JF: Phys Rev Lett. 1985, 55: 2409. 10.1103/PhysRevLett.55.2409
Burnham DC, Weinberg DL: Phys Rev Lett. 1970, 25: 84. 10.1103/PhysRevLett.25.84
Bergeal N, Vijay R, Manucharyan VE, Siddiqi I, Schoelkopf RJ, Girvin SM, Devoret MH: Nat Phys. 2010, 6: 296. 10.1038/nphys1516
Bergeal N, Schackert F, Frunzio L, Devoret MH: Phys Rev Lett. 2012., 108: Article ID 123902 Article ID 123902
Roch N, Flurin E, Nguyen F, Morfin P, Campagne-Ibarcq P, Devoret MH, Huard B: Phys Rev Lett. 2012., 108: Article ID 147701 Article ID 147701
Ho Eom B, Day PK, LeDuc HG, Zmuidzinas J: Nat Phys. 2012, 8: 623. 10.1038/nphys2356
Yaakobi O, Friedland L, Macklin C, Siddiqi I: Phys Rev B. 2013., 87: Article ID 144301 Article ID 144301
Wilson CM, Johansson G, Pourkabirian A, Simoen M, Johansson JR, Duty T, Nori F, Delsing P: Nature. 2011, 479: 376. 10.1038/nature10561
Wustmann W, Shumeiko V: Phys Rev B. 2013., 87: Article ID 184501 Article ID 184501
Dykman M, Krivoglaz M: Phys A, Stat Mech Appl. 1980, 104: 480. 10.1016/0378-4371(80)90010-2
Marthaler M, Dykman MI: Phys Rev A. 2006., 73: Article ID 042108 Article ID 042108
Siddiqi I, Vijay R, Pierre F, Wilson CM, Metcalfe M, Rigetti C, Frunzio L, Devoret MH: Phys Rev Lett. 2004., 93: Article ID 207002 Article ID 207002
Vijay R, Devoret MH, Siddiqi I: Rev Sci Instrum. 2009., 80: Article ID 111101 Article ID 111101
Mallet F, Ong FR, Palacios-Laloy A, Nguyen F, Bertet P, Vion D, Esteve D: Nat Phys. 2009, 5: 791. 10.1038/nphys1400
Gardiner CW, Collett MJ: Phys Rev A. 1985, 31: 3761. 10.1103/PhysRevA.31.3761
Walls DF, Milburn GJ: Quantum Optics. Springer, Berlin; 1994.
Göppl M, Fragner A, Baur M, Bianchetti R, Filipp S, Fink JM, Leek PJ, Puebla G, Steffen L, Wallraff A: J Appl Phys. 2008., 104: Article ID 113904 Article ID 113904
Braunstein SL, van Loock P: Rev Mod Phys. 2005, 77: 513. 10.1103/RevModPhys.77.513
Caves CM: Phys Rev D. 1982, 26: 1817. 10.1103/PhysRevD.26.1817
Wallquist M, Shumeiko VS, Wendin G: Phys Rev B. 2006., 74: Article ID 224506 Article ID 224506
Bourassa J, Beaudoin F, Gambetta JM, Blais A: Phys Rev A. 2012., 86: Article ID 013814 Article ID 013814
Girvin SM: Lectures delivered at Ecole d’Eté Les Houches. To be published by Oxford University Press, 2011. Girvin SM: Lectures delivered at Ecole d’Eté Les Houches. To be published by Oxford University Press, 2011.
Imoto N, Haus HA, Yamamoto Y: Phys Rev A. 1985, 32: 2287. 10.1103/PhysRevA.32.2287
Sanders BC, Milburn GJ: Phys Rev A. 1989, 39: 694. 10.1103/PhysRevA.39.694
Santamore DH, Goan H-S, Milburn GJ, Roukes ML: Phys Rev A. 2004., 70: Article ID 052105 Article ID 052105
Buks E, Yurke B: Phys Rev A. 2006., 73: Article ID 023815 Article ID 023815
Helmer F, Mariantoni M, Solano E, Marquardt F: Phys Rev A. 2009., 79: Article ID 052115 Article ID 052115
Johnson BR, Reed MD, Houck AA, Schuster DI, Bishop LS, Ginossar E, Gambetta JM, DiCarlo L, Frunzio L, Girvin SM, Schoelkopf RJ: Nat Phys. 2010, 6: 663. 10.1038/nphys1710
Suchoi O, Abdo B, Segev E, Shtempluck O, Blencowe MP, Buks E: Phys Rev B. 2010., 81: Article ID 174525 Article ID 174525
Kirchmair G, Vlastakis B, Leghtas Z, Nigg SE, Paik H, Ginossar E, Mirrahimi M, Frunzio L, Girvin SM, Schoelkopf RJ: Nature. 2013, 495: 205. 10.1038/nature11902
Koch J, Yu TM, Gambetta J, Houck AA, Schuster DI, Majer J, Blais A, Devoret MH, Girvin SM, Schoelkopf RJ: Phys Rev A. 2007., 76: Article ID 042319 Article ID 042319
Ferguson DG, Houck AA, Koch J: Phys Rev X. 2013., 3: Article ID 011003 Article ID 011003
Acknowledgements
The authors would like to acknowledge helpful discussions with Alexandre Blais, Vitaly Shumeiko and Sebastian Schmidt. This work was supported by the European Research Council (ERC) through a Starting Grant and by ETHZ.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare they have no competing interests.
Authors’ contributions
All authors contributed equally to the work in this article.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Eichler, C., Wallraff, A. Controlling the dynamic range of a Josephson parametric amplifier. EPJ Quantum Technol. 1, 2 (2014). https://doi.org/10.1140/epjqt2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjqt2