the leading model of gamma - ray bursts ( grbs ) involves a relativistic ( @xmath4 , where @xmath5 is the jet lorentz factor ) fireball jet , where the observed gamma rays are produced by radiation from fermi - accelerated electrons in optically thin shocks ( for reviews see e.g.  @xcite ) . in the most popular model for the more common long - duration grbs , the so - called collapsar model , the core of a massive star collapses to a black hole or neutron star , driving a highly relativistic jet which breaks out of the star  @xcite . within this scenario , high - energy neutrinos from photomeson interactions of accelerated protons
have been studied  @xcite . with @xmath6 size
neutrino detectors such as icecube  @xcite and km3net  @xcite being constructed , we are entering an age when we can test various predictions for high - energy neutrinos .
waxman and bahcall predicted @xmath7 neutrinos from external reverse shocks  @xcite , and similar predictions have been made in refs .
@xcite for the forward shock .
in addition , internal shocks can occur while the relativistic jet is still in the star .
mszros and waxman predicted a @xmath8 neutrino precursor burst from such a set up  @xcite .
it is expected that a precursor neutrino burst is present even if the jet does not successfully break out of the star , producing a choked jet that is dark in gamma rays and bright in neutrinos .
inferring from the observed rate of grbs , highly relativistic ( @xmath9 ) jets occur in perhaps @xmath10 of core - collapse supernovae  @xcite .
however , the jet signature may be more common , and a significantly higher portion of supernovae may be endowed with slower , mildly relativistic jets ; according to late - time radio observations , perhaps as high as a few @xmath11  @xcite .
the recent detection of low - luminosity grb 060218 also suggests from detection rates that low - luminosity grbs , with mildly relativistic ejecta , are more common than conventional high - luminosity grbs by factors as large as @xmath12  @xcite . the occurrence rate is high enough that despite their lower luminosity , their contribution to the neutrino background is comparable to or larger than conventional grbs  @xcite .
in addition to being more frequent , mildly relativistic jets are more baryon rich , which work positively for neutrino emission .
indeed , the detectability of neutrinos from internal shocks of mildly relativistic jets has been studied by razzaque , mszros , and waxman  @xcite and extended by ando and beacom  @xcite ( see also  @xcite ) . in these scenarios ,
internal shocks occur at @xmath13@xmath14 , smaller than the typical radius of type ib progenitors . accelerated protons interact via proton - photon ( @xmath15 ) and proton - proton ( @xmath16 ) interactions to produce mesons , which consequently produce neutrinos on decay .    in this paper
we investigate neutrinos arising from protons that are accelerated in the reverse shock , while the relativistic jet is still propagating inside the star .
we adopt a type ib progenitor as our model star , but we discuss how our results apply to type ic and ii progenitors also . the reverse shock is strongly motivated because it is a generic prediction of the initial stages of any jet moving through the star , regardless of its lorentz factor , and whether it breaks out of the star or not .
therefore , we consider two types of supernovae : one is the perhaps more numerous supernova containing a mildly relativistic ( @xmath17 ) jet , and the other is a supernova containing an ultrarelativistic ( @xmath18 ) jet that can ultimately cause the observed grbs . in both cases we consider choked as well as successful jets , depending on its duration .
we show that protons can be accelerated to @xmath0@xmath1 gev in the reverse shock .
the low and high - energy protons produce mesons efficiently by @xmath16 and @xmath15 interactions respectively . for energies in between , cooling mechanisms are faster , and result in a characteristic suppression feature in the meson spectrum .
mesons that are produced experience significant cooling due to the high proton and photon densities in the shocked jet . however , we still expect a reasonable number of neutrino events .
for instance , for a choked mildly relativistic jet with isotropic energy @xmath19 erg at 10 mpc , we expect @xmath20 neutrino events at icecube from the reverse shock .
interestingly , assuming the same jet lorentz factor , we find that a lower isotropic energy of @xmath21 erg produces more neutrino events ( @xmath22 events ) .
this is due to meson cooling .
specifically , we find that neutrinos are favored by high - lorentz factor , low - luminosity , and long - duration jets , as meson cooling is less significant for such jets .
the rate of core - collapse supernova within 10 mpc is expected to be @xmath313 yr@xmath23  @xcite , implying particle acceleration in reverse shocks can potentially be tested in the next decades .
neutrino emission from the reverse shock is potentially important for short - duration jets that are completely shocked and do not break out of the star , i.e. , choked jets .
we argue that while it is possible that internal shocks occur prior to jet shocking ( in particular for low @xmath5 jets ) , the allowed range of radii for internal shocks to occur is constrained from above by the reverse shock , and below by the fact that the stellar material is opaque to high - energy neutrinos .
this suggests that neutrino emission from choked jets are of similar magnitude to our estimates .
this paper is organized as follows . in sec .
ii we discuss the evolution of the jet as it propagates through the star . in sec .
iii we treat proton acceleration , and suppression of neutrino emission due to proton and meson cooling . in sec .
iv we utilize the results of the previous sections to investigate the detectability of high - energy neutrinos from reverse shocks . finally , we finish with discussions in sec .
v and conclusions in sec .
we first discuss the evolution of a relativistic @xmath24 jet as it strikes the stellar matter .
the aim is to obtain the evolution of the jet head and reverse shock velocities .
there are several numerical simulations on the propagation of an initially hot jet in the literature  @xcite , as well as analytic treatments of initially cold jets  @xcite .
we follow and extend the latter .
generally , two shocks form when a jet strikes external matter : a reverse shock that decelerates the head of the jet to a lorentz factor @xmath25 , and a forward shock that accelerates the external material to @xmath25 .
the jet and stellar material moving at @xmath25 are separated by a contact discontinuity , and are in pressure balance .
equations governing the evolution of the shocks are ( for @xmath24 )  @xcite @xmath26 where @xmath27 is the proton mass and @xmath28 is the speed of light .
thermodynamic quantities @xmath29 and @xmath30 ( particle number density and internal energy ) are measured in the fluids rest frames , and we have used the following notations : @xmath31 ( jet , unshocked ) , @xmath32 ( jet head , shocked ) , @xmath33 ( stellar , shocked ) , and @xmath34 ( stellar , unshocked ) .
lorentz factors are measured in the lab frame , except for @xmath35 which is the lorentz factor of the jet head measured in the frame of the unshocked jet , @xmath36 where @xmath37 is the velocity as a fraction of @xmath28 .
the jet particle density is given by ( for a jet with constant opening angle ) @xmath38 where @xmath39 is the isotropic jet luminosity .
we obtain @xmath25 by equating the pressure of the shocked jet head @xmath40 and the pressure of the shocked stellar matter @xmath41  @xcite . from eq .
( [ shockequations ] ) , this requires us to solve @xmath42 by substitution of eqs .
( [ gamma3bar ] ) and ( [ n_j ] ) .
in what follows we use @xmath43 from the progenitor model e20 from ref .
@xcite , an initially @xmath44 rotating star that is reduced to @xmath45 by mass - loss .
it has a @xmath46 fe core of radius @xmath47 , a @xmath48 he core that extends to @xmath49 , and a surrounding h envelope .
however , the results that we present here apply qualitatively to other models of different masses , with or without h envelopes . as an illustration of this point
we include other progenitor models in fig .
[ velocity ] .
jet head velocity @xmath50 as a function of radius . e15 and e20 are two progenitor models of mass 15@xmath51 and 20@xmath51 respectively , as detailed in the text .
l51l53 denotes the isotropic equivalent luminosity of the jet , @xmath52@xmath53 . for the l53 case ,
two jet lorentz factors are plotted : mildly relativistic ( @xmath54 , solid ) and ultrarelativistic ( @xmath55 , dashed ) . for l52 and l51 ,
mildly relativistic jets are shown.,width=326 ]    due to rotation , the mass along the rotational axis becomes substantially lower than the equatorial , and the rotational axis is an easy escape route for a jet .
we inject jets at a radius @xmath56 along the rotational axis , which corresponds to an enclosed mass of @xmath57 .
since the explosion and jet mechanisms are still under debate , they are treated parametrically with the following parameters : the isotropic equivalent energy @xmath58 and duration @xmath59 ( which yield the isotropic equivalent luminosity through @xmath60 ) , lorentz factor @xmath5 , and variability @xmath61 .
we consider two types of jets , ( i ) mildly relativistic jets with @xmath54 , and ( ii ) ultrarelativistic jets with @xmath55 .    in fig .
[ velocity ] we show the jet head velocity @xmath62 multiplied by @xmath25 as a function of radius , obtained from solving eq .
( [ v_h ] ) . for @xmath63 ,
we show our two jet types @xmath54 and @xmath55 , for two progenitor models .
the first progenitor model is e20 as described above , and the second is e15 , a similarly rotating star but with initial mass @xmath64  @xcite .
we see that the evolution of @xmath62 is independent of @xmath5 , and that the initial jet head velocity is always subrelativistic .
the latter is because the external density is very much larger than the jet density , @xmath65 , for conceivable set of jet parameters .
as the jet continues to propagate , the stellar density falls faster with radius ( @xmath66 ) than the jet density ( @xmath67 ) , and the jet head effectively accelerates , reaching mildly relativistic values at @xmath68 a few @xmath69 .
although the initial value of @xmath62 scales somewhat with @xmath39 , the radius at which the velocity reaches @xmath28 is only weakly dependent on jet parameters .
this is because the jet head accelerates to relativistic velocities when the external stellar density suddenly drops , at the edge of the he core .
in fact , rewriting eq .
( [ v_h ] ) in the relativistic limit ( @xmath70 ) yields @xmath71  @xcite .
note that the kinks in fig .
[ velocity ] corresponding to composition changes in the progenitor model are model dependent , e.g. , an abundance of oxygen resides at @xmath72 .
however , the general features of @xmath62 that concern us , including its initially @xmath73 value and its reaching @xmath74 at @xmath75 , are model independent .
we therefore consider only the e20 progenitor model in the rest of our paper .
[ t ]    .[table : times]summary of choked jet characteristic times . shown
are the core crossing time @xmath76 , and the jet duration time @xmath77 required for reverse shock crossing to occur at @xmath78 cm .
these values apply to both mildly relativistic and ultrarelativistic jets . [ cols="^,^,^,^,^",options="header " , ]     to estimate @xmath79 for a choked jet , we note that @xmath80 during propagation in the core .
the reverse shock thus crosses the entire length of the jet in a time @xmath81 , where @xmath82 is the jet length in the lab frame .
we obtain @xmath79 as the distance the jet head has traveled in time @xmath59 . as an example , adopting a jet luminosity @xmath83 , we obtain @xmath84 s ; a jet with duration @xmath85 is therefore choked , and we find that @xmath86 for both the mildly relativistic as well as ultrarelativistic jet cases . for calculation purposes we set @xmath87 .
this requires a jet duration of @xmath88@xmath89 , depending on the value of @xmath39 .
we summarize duration times in table [ table : times ] .
we note that we could have in principle chosen a smaller @xmath79 due to smaller @xmath77 .
however , as we discuss in sec .
v , the stellar matter becomes opaque to high - energy neutrinos for radii @xmath90 cm , and small radii are uninteresting .
we also note here that rotational effects decrease the matter density along the rotational axis , increasing the jet head velocity at small radii and thereby shortening the core crossing time @xmath76 .
a decrease in density also decreases neutrino opacity .
however , the effect is not expected to be large , because density modifications typically occur at radii smaller than or equal to @xmath91  @xcite .
next we discuss jets that successfully break out of the star , with @xmath92 .
these jets are progressively shocked as they move through the progenitor , but the reverse shock does not cross the entire jet . therefore , for successful jets , only the front portion of the jet is shocked .
we assume @xmath79 just inside of the he core radius . we stress that for successful jets , @xmath79 is _ not _ the point of reverse shock crossing , but the last point of shocking by the reverse shock .
a few words about the jet models must be noted .
first , we have assumed a conical jet geometry in eq .
( [ n_j ] ) .
this is a good approximation for jets propagating in a sufficiently sparse medium , such as the outer regions of stars . from analytic studies , the jet opening angle @xmath93 for @xmath94 and @xmath95 constant for @xmath96  @xcite .
we have therefore injected jets at a radius @xmath97 with profile eq .
( [ n_j ] ) . on the other hand , under conditions such as propagation in the dense iron core and/or when strong magnetic fields are associated with the launch of the jet , the high external matter density and/or magnetic hoop stress work to collimate the jet .
numerical studies find that jets remain tightly collimated in the iron core  @xcite , and further beyond the stronger the magnetic field  @xcite . for jets that remain strongly collimated ,
the reverse shock becomes increasingly less important due to the more penetrating nature of the jet .
further , it is less likely for these jets to form choked jets . for illustration , consider an extreme case where a jet of luminosity @xmath98 erg @xmath99 remains tightly cylindrical all the way until break out , with constant cylindrical radius @xmath100 cm  @xcite .
for such a jet , we find that @xmath101 6 s and @xmath102 s.    second , we have assumed a type ib progenitor .
a jet with @xmath103 erg s@xmath23 propagating through a type ic progenitor would clear the core in @xmath104 seconds , and jets are more likely to break out .
third , we have defined successful jets as those which break out of the star .
however , these do not necessarily always produce observable grbs  for small @xmath5 and @xmath61 , internal shocks will occur while the jet is optically thick to photons .
for example , @xmath105 cm ( we define @xmath106 for a quantity @xmath107 in cgs units ) , while the jet typically becomes optically thin to photons at @xmath108^{1/2 } \sim 10^{14}l_{\rm{iso},52}^{1/2}\gamma_{j,1}^{-1/2}$ ] cm , where @xmath109 is the thomson scattering cross section .
the jet is still transparent to high - energy neutrinos .
we summarize our various jet definitions in table  [ table : jets ] .
to conclude , we use the following crossing radii : ( i ) @xmath110 for choked jets and ( ii ) @xmath111 for successful jets .
in this section we discuss proton acceleration and neutrino production through decays of pion and kaon mesons .
the reverse shock is likely to be collisionless , as we discuss in section iiia .
if first - order fermi acceleration is realized , a fraction of the jet protons are accelerated , with a spectrum @xmath112 where the power - law index @xmath113 in the ultrarelativistic shock limit  @xcite , and @xmath114 for a nonrelativistic shock . in this paper , we optimistically adopt @xmath114 for all of our figures , but we present numerical results for @xmath115 2.1 , 2.3 , and 2.5 also .
the minimum proton energy would be @xmath116 , but otherwise it is not well known ; we assume @xmath117 , but its exact value does not affect our results .
the maximum proton acceleration energy is determined by comparing the proton s acceleration time scale to its energy loss time scales .
ando and beacom have shown that the photopion process is the most competitive cooling mechanism in the proton energy range of interest  @xcite .
we therefore focus on photopion production , and refer the reader to appendix  a for treatments of other cooling mechanisms such as synchrotron cooling .
we will however for completeness display all cooling mechanisms in fig .
[ cooling ] .
note that the bethe - heitler process turns out to be an important cooling process for @xmath16 interactions , as we will discuss in the next subsection .
the energy density in the shocked jet head is @xmath118 .
we have specifically used primes ( @xmath119 ) to remind that these quantities are in the jet head comoving frame .
we assume that a fraction @xmath120 of the shocked plasma internal energy is converted to relativistic electrons and magnetic fields , in analogy to grbs .
this gives the comoving magnetic field strength as @xmath121^{1/2 } ( \bar{\gamma}_h/\gamma_j ) \nonumber \\ & \approx & 1 \times 10^{10 } \ , \left ( \frac{\varepsilon_{b,-1}^{1/2}l_{\rm{iso},52}^{1/2 } } { r_{\times,10.7 } } \right ) \left ( \frac{\bar{\gamma}_{h,1}}{\gamma_{j,1 } } \right ) \ , \mathrm{g},\end{aligned}\ ] ] where we have assumed @xmath80 since @xmath122 for the jet scenarios studied here .
one can obtain from eq .
( [ gamma3bar ] ) that for example , @xmath123 yields @xmath124 .
the relativistic electrons immediately lose energy by synchrotron radiation , and the photons thermalize due to the large opacity @xmath125 , where @xmath126 and @xmath127 is the proper thickness of the shocked jet plasma . here
we adopt @xmath128 , where @xmath129  @xcite .
the black - body temperature is then @xmath130^{1/4 }   ( \bar{\gamma}_h/\gamma_j)^{1/2 } \nonumber \\ & \approx & 13 \ , \left ( \frac{\varepsilon_{e,-1}^{1/4}l_{\rm{iso},52}^{1/4 } } { r_{\times,10.7}^{1/2 } } \right )   \left ( \frac{\bar{\gamma}_{h,1}}{\gamma_{j,1 } } \right)^{1/2 } \,\mathrm{kev},\end{aligned}\ ] ] and the average photon energy is @xmath131 .
the average density is @xmath132^{3/4 } ( \bar{\gamma}_h/\gamma_j)^{3/2 } \nonumber \\ & \approx & 7 \times 10^{25 } \ , \left ( \frac{\varepsilon_{e,-1}^{3/4}l_{\rm{iso},52}^{3/4 } } { r_{\times,10.7}^{3/2 } } \right )   \left(\frac{\bar{\gamma}_{h,1}}{\gamma_{j,1 } } \right)^{3/2 } \,\mathrm{cm^{-3}}.\end{aligned}\ ] ]    the total proton cooling time is determined by the inverse sum @xmath133 , where @xmath134 runs through the various cooling mechanisms . as discussed in the opening paragraph ,
we focus on cooling by photopion production .
protons interact with the dense photon field eq .
( [ photons2 ] ) , and for sufficient proton - photon centre of inertia energies mesons are produced .
we determine the cooling time scale by @xmath135 .
we assume the conventional inelasticity @xmath136/2 $ ] where @xmath33 is the invariant mass of the system .
we fit @xmath137 to  @xcite , and use a blackbody photon spectrum of temperature eq .
( [ photons ] ) .
the time scale is plotted in fig .
[ cooling ] .    the acceleration time scale in the jet head frame
is given by @xmath138 , where @xmath139 is a constant if we assume the diffusion coefficient is proportional to the bohm diffusion coefficient , and @xmath140 is the larmor radius .
@xmath141 is a fairly conservative value , and @xmath142 can be achieved for mildly relativistic shocks  @xcite .
we use @xmath143 which gives agreement with previous studies  @xcite . using the magnetic field in eq .
( [ bfield ] ) , we obtain @xmath144 where @xmath145 is measured in gev .
inverse of proton cooling and acceleration time scales , in the shocked jet head frame , as functions of the proton energy .
the case of reverse shock crossing at @xmath146 is plotted for supernova jet model b , with other parameters @xmath54 , @xmath147 , and @xmath148 .
cooling mechanisms shown are : synchrotron ( blue , thin dashed ) , inverse - compton ( green , dot - dashed ) , bethe - heitler ( green , dotted ) , proton - proton ( black , long dashed ) , and photopion ( black , thick dashed ) .
also labeled are threshold energies @xmath149 and @xmath150 ( see sec .
photopion limits the proton energy to @xmath151.,width=326 ]    in fig .
[ cooling ] , we plot the proton acceleration and cooling time scales as functions of the proton energy , both in the jet head comoving frame . acceleration is stronger than cooling at low proton energies , but photopion cooling becomes stronger at higher energies .
it is evident that photopion cooling is the strongest cooling mechanism . equating @xmath152 , we obtain the maximum proton energy as @xmath153 .
the parameter dependency is @xmath154 in the flat region of @xmath155 .
we note that the reverse shock is likely to be collisionless , because the proton plasma frequency @xmath156 is larger than the typical radiation and particle collision frequencies . for example , the frequency of proton - photon collision is @xmath157 , while the frequency of proton - proton collision is @xmath158 .
depending on its energy , protons produce mesons via @xmath16 and/or @xmath15 interactions .
we consider the pion and kaon mesons .
the multiplicity in each @xmath16 and @xmath15 interaction is taken to be @xmath159 for pions and @xmath160 for kaons ; this matches the required ratio in the relevant energies  @xcite .
generally , kaons produce higher energy neutrinos  @xcite , but we additionally find that they can dominate those from pions because of less efficient kaon cooling .    for low energy protons , @xmath16 interactions dominate .
meson production by @xmath16 interaction is efficient since the opacity is very high , @xmath161 , where @xmath162 and @xmath163  @xcite .
proton cooling by pair production ( bethe - heitler , hereafter bh , see appendix for treatment ) overtakes @xmath16 interactions at higher energies ( see fig .
[ cooling ] ) .
we define this transition as @xmath164 . at even higher energies @xmath15 dominates over bh cooling , and we define this transition energy as @xmath165 .
the @xmath15 interaction is efficient at producing mesons since the opacity is @xmath166 , where we use @xmath167  @xcite . for model b in table i , @xmath168 gev and @xmath169 gev .
we assume mesons are produced with 20% of the parent proton energy . in the absence of competing cooling mechanisms , the meson spectrum follows the initial proton spectrum . on the other hand , when other cooling processes are faster than meson - yielding @xmath16 and @xmath15 interactions , the meson spectrum is suppressed by an amount equal to the ratio of time scales .
we find that the bethe - heitler process is a competitive cooling mechanism for @xmath16 interactions , and the resulting meson spectrum will be suppressed by a factor given as @xmath170 we here approximate @xmath171 as a power - law @xmath172 for @xmath173 , where @xmath174 is a constant .
we find that the value of @xmath174 is fairly model independent , with values @xmath175 ( 1.2 ) for @xmath5 = 10 ( 100 ) for model b.      in this section we discuss meson cooling and the neutrino spectrum .
mesons cool similarly to protons , by radiative and hadronic ( collisions with protons ) processes .
the meson cooling time scales are analogous to those for the proton and can be summarized as @xmath176 $ ] and @xmath177 , where @xmath178 and @xmath179 are the meson mass and energy , and @xmath180 and @xmath181 are the photon and magnetic field energy densities ( we assume inverse - compton process is in the thompson regime ; synchrotron dominates at high energies ) . here , @xmath182
@xcite is the cross section for meson - proton collisions and @xmath183  @xcite is the energy lost by the meson per collision .
the total cooling time is @xmath184 , and the time scales are @xmath185 and @xmath186 for pions . for kaons , the radiative cooling time scale is longer , because @xmath187 .
hadronic cooling dominates at low proton energies , and radiative cooling becomes important at higher energies .    charged pions and kaons decay into neutrinos through @xmath188 , with the muon neutrino taking 1/4 of the meson energy .
we neglect secondary neutrinos from muon decays since muons immediately lose energy by radiative cooling .
we define the break energy @xmath189 by equating @xmath190 , where @xmath191 and @xmath192 are the meson lorentz factor and proper lifetime .
below @xmath189 , neutrinos have a typical energy @xmath193 and a flat power per decade neutrino spectrum .
above @xmath189 , the spectrum is suppressed by a factor @xmath194 .
next , we define @xmath195 when radiative cooling begins to dominate , by equating @xmath196 . due to relativistic effects
, the neutrino energy in the observer frame is related to the parent meson energy in the jet head frame as @xmath197 and @xmath198 .
we define the suppression function @xmath199 as @xmath200 for some supernova jet parameters , it is possible that meson goes from decay dominated straight to radiation cooling dominated . for these cases
we define the break energy @xmath201 by equating @xmath202 , and the suppression function @xmath199 as @xmath203    the break energies for neutrinos from pion decay , in a choked mildly relativistic jet , are @xmath204 the parameter dependencies are @xmath205 and @xmath206 for the first and second break energies , respectively . on the other hand ,
break energies for neutrinos from kaon decay are higher due to the higher kaon mass and faster decay time ; @xmath207 the parameter dependencies are the same as for pions .
the break energies are smaller than those of internal shocks  @xcite with the same jet parameters , by as much as @xmath32 orders of magnitude .
this is a combination of the higher proton and photon densities in the shocked jet causing more effective cooling , and the fact that the neutrino energies in the shocked jet head receive no lorentz boosting .
next , the break energies for a choked ultrarelativistic jet are @xmath208 which can also be obtained from scaling eqs .
( [ break12])([break34 ] ) by their dependencies on @xmath5 importantly for neutrino emission , suppression due to radiative cooling is independent of @xmath5 , since the product of @xmath209 is unchanged .
also , note that @xmath210 , indicating that radiative cooling dominates the cooling process for pions . in this case , @xmath211    we show in fig .
[ suppression1 ] the resulting suppression factors for neutrinos from kaons , for our mildly relativistic and ultrarelativistic jets .
suppression is less intense for relativistic jets .
we also shade in for reference the energy range where the meson spectrum is suppressed by bethe - heitler cooling , from 5140 gev . note that for pions , radiative cooling dominates for energies of interest .
suppression factor as a function of neutrino energy , for neutrinos from kaon decays . shown are the choked mildly relativistic jet ( solid lines ) and choked ultrarelativistic jet ( dashed lines ) .
the kaon experiences less cooling for an ultrarelativistic jet . note that pions are always in the radiation cooling regime for neutrino energies shown .
we also shade in for reference the energies where bethe - heitler cooling is competitive and suppresses the meson spectrum.,width=326 ]
cumulative muon event number from choked jets with model b parameters .
thin - black ( thick - green ) lines denote pion ( kaon ) contributions .
the majority originates from kaons .
solid lines denote a mildly relativistic jet and a supernova distance of @xmath212 .
dashed lines denote an ultrarelativistic jet and two supernova distances - @xmath212 for illustrative purposes and @xmath213 for event calculations , as labeled . due to weaker kaon cooling ,
an ultrarelativistic jet yields more neutrinos than a mildly relativistic one .
the spectrum becomes almost flat around @xmath3 100 gev , due to bethe - heitler suppression of the meson spectrum .
the total number of events expected above the detector threshold energy is about 6 for the mildly relativistic case . for the ultrarelativistic case with distance
@xmath213 we expect 0.02 events .
also shown is the atmospheric neutrino background over 1 day in a @xmath214 circle.,width=326 ]    in this section we estimate the emission and detectability of high - energy neutrinos from reverse shock accelerated protons .
we consider a supernovae at a distance of @xmath215 , and use the code anis ( all neutrino interaction generator )  @xcite to calculate the neutrino induced muon spectrum at a neutrino detector with an effective area of @xmath216 , i.e.  an icecube class detector utilizing upgoing muons .
first , we estimate the fluence of neutrinos as  @xcite @xmath217 where @xmath218 is the meson multiplicity ( @xmath159 for pions and @xmath160 for kaons ) , @xmath219 is the branching ratio of meson decay into neutrinos ( @xmath159 for pions and @xmath220 for kaons ) , and @xmath221 is the fraction of proton energy carried by neutrinos in the absence of energy loses ; @xmath222 for pions and @xmath223 for kaons , since neutral and charged mesons are produced with roughly equal probability , and muon neutrinos carry roughly @xmath223 ( @xmath224 ) of the pion ( kaon ) energy in meson decay .
the functions @xmath225 and @xmath226 are the suppression factors due to proton and meson coolings , respectively , and the @xmath227 factor normalizes the proton spectrum to the jet energy ( for @xmath228 ) .    in fig .
[ muon ] we show the cumulative spectrum of neutrino induced muons , from a supernova with a choked jet .
we show three cases - a supernova distance of @xmath229 possessing a mildly relativistic @xmath54 jet ( solid lines ) , a supernova distance of @xmath229 possessing an ultrarelativistic @xmath55 jet ( dashed lines ) , and a similarly ultrarelativistic jet but supernova distance @xmath230 ( dashed lines labeled 400 mpc ) .
we use model b with parameters @xmath147 and @xmath231 .
we take into account the muon range , which effectively increases the detector volume .
we evaluate the muon energy @xmath232 when it enters the detector if it is produced outside , or at the production point if it is inside . considering a detection threshold muon energy of @xmath233
, we expect about 6 muon events from a 10 mpc supernova possessing a choked mildly relativistic jet .
similarly , we expect about 28 from an ultrarelativistic jet .
we have plotted the ultrarelativistic jet case at an optimistic distance of 10 mpc to illustrate the effect of meson cooling .
although the pion contributions are similar between the two jets , the kaon contribution is larger in the ultrarelativistic case due to the weaker suppression . in the remainder
we assume the distance 400 mpc for ultrarelativistic jet calculations . for model
b , this yields 0.02 events .
note that since the muon spectrum is steep , a lower detector energy threshold will yield much higher events .
additionally , one may observe the suppression due to bethe - heitler cooling , which is expected to manifest below @xmath234 gev .
for a steeper initial proton spectrum , motivated from studies of ultrarelativistic shocks , the expected event is 4 ( 2 , 1 ) for @xmath235 ( 2.3 , 2.5 ) . in all jets ,
the majority of neutrinos come from kaons , because kaons are more massive than pions and suffer less cooling ( the radiative cooling time scale is @xmath236 ) . also , the kaon decay time is slightly shorter than the pion .
the events discussed here arrive in a @xmath237 34 second time bin and @xmath3 3@xmath238 angular bin , allowing very strong rejection of atmospheric neutrino backgrounds .
we show in the figure the atmospheric neutrino background within @xmath159 day and a @xmath214 bin for comparison .
the expected neutrino signal from a choked mildly relativistic jet easily exceeds the atmospheric neutrino background .
same as fig .
[ muon ] but for model c parameters .
approximately 17 ( 0.04 ) events are expected above the detector threshold for a single mildly relativistic ( ultrarelativistic ) jet and supernova distance @xmath212 ( @xmath213).,width=326 ]    same as fig .
[ muon ] but for jets with successful jets . the supernova and jet parameters are @xmath239 , @xmath240 , and @xmath241 .
about 36 ( 0.1 ) events are expected above the detector threshold muon energy from a mildly relativistic ( ultrarelativistic ) jet supernova at a distance @xmath212 ( @xmath213 ) away.,width=326 ]    next we discuss a scan of supernova and jet parameters .
first , as shown above , a larger @xmath5 reduces kaon suppression and increases neutrino emission .
next , we consider jet model a of table  [ table : times ] , i.e. , with larger @xmath39 but smaller @xmath59 , and an overall larger @xmath58 . since the neutrino fluence scales with @xmath58 , one naively expects more neutrinos .
however , the expected muon event rate is in fact smaller , @xmath242 ( @xmath243 ) for a mildly relativistic ( ultrarelativistic ) jet .
this is because of stronger meson cooling : the proton and photon densities in the shocked jet scale as @xmath244 and @xmath245 , so a larger @xmath39 results in stronger meson cooling .
we conclude that neutrino emission from choked jets are more favorable from low - luminosity and long - duration jets , even if @xmath58 is smaller ( assuming the same jet lorentz factor ) . indeed
, model c predicts 17 ( 0.04 ) events at icecube .
we show the cumulative muon event spectrum for jet model c in fig .
[ muonc ] .    in fig .
[ muon2 ] we show the same as fig .
[ muon ] but for successful jets .
we show a mildly relativistic @xmath54 jet at @xmath212 ( solid lines ) , and an ultrarelativistic @xmath55 jet at @xmath213 ( dashed lines ) .
the other parameters used are @xmath246 , @xmath240 , @xmath241 , @xmath147 , and @xmath247 . with these parameters , @xmath84 s , and the cold jet
does indeed break out of the star .
we introduce an additional factor of 1/2 to the neutrino fluence for successful jets , because only about @xmath248 of the jet length is shocked by the reverse shock .
note that our chosen parameters have the equivalent @xmath58 to jet model a. however , we expect @xmath249 times more neutrino events than jet model a , because @xmath39 is a factor 10 smaller and @xmath79 is a factor 2 larger . indeed , evaluating the muon event rate yields about 36 ( 0.1 ) events per mildly relativistic ( ultrarelativistic ) jet . as in the case of choked jets ,
most of the contribution comes from kaons .
these events will show temporal and angular clustering , and in addition , these will appear as precursor signals to grb photons , allowing very strong rejection of background .
our result is that neutrinos from the reverse shock are strongly suppressed due to proton and meson cooling .
however , we still expect a reasonable number of events at neutrino detectors .
one reason for this is that proton cooling affects only energies near and below current detector threshold energies .
this produces a suppression feature in the spectrum but does not decrease event counts at icecube . for our models a , b , and c in table  [ table : times ] , we expect 2 , 6 , and 17 muon events respectively at a icecube class detector . here , we have calculated for a supernova @xmath212 away , possessing a @xmath54 mildly relativistic jet , and a @xmath228 initial proton spectrum .
if we assume a steeper index , @xmath115 2.1 ( 2.3 , 2.5 ) , we expect 4 ( 2 , 1 ) events for jet model b. in all cases , the neutrino flux is dominated by decay of kaons , since cooling is less efficient for kaons than pions
.    since neutrino emission from the reverse shock is suppressed by meson cooling , we find that neutrino emission is favored by parameter sets for which meson cooling is weaker  high - lorentz factor , low - luminosity , and long - duration jets .
this fact is readily observable from the neutrino event predictions of models a , b , and c. as further illustration , if we assume a @xmath55 ultrarelativistic jet we expect 11 , 28 , and 55 events from jet models a , b , and c respectively .
note that supernova and jet parameters are expected to show considerable scatter .
hence , given neutrino emission is cooling suppressed , the event count need not necessarily correlate with the total burst energy .
neutrino emission from the reverse shock is particularly important for choked jets , which are completely shock decelerated by the reverse shock .
is it possible for choked jets to produce more neutrino emission than our above estimates ?
if the duration is long enough such that @xmath250 , it is possible that internal shocks accelerate protons in the unshocked part of the jet , while the jet head is being shocked .
the internal shock accelerated protons do not travel far , producing mesons efficiently via @xmath16 interactions with the large target jet proton density  @xcite . for high - lorentz factor jets , the protons can travel and interact with termination shock photons  @xcite . in both cases meson cooling
is less severe than in the reverse shock , and we can expect more neutrino emission .
however , internal shocks must occur at large enough radii that the stellar material is not opaque to high - energy neutrinos .
the neutrino opacity is @xmath251 , where @xmath252 is the column number density of the star and @xmath253 is the sum of neutrino charged - current and neutral - current cross sections , which is approximately @xmath254 in the energy range in question . for a @xmath255 neutrino emitted at @xmath256 , @xmath257 and @xmath258  @xcite , and we find that @xmath259 . for a @xmath260 neutrino , @xmath261 .
thus , only neutrinos with energies @xmath262 can leave the star from @xmath263 . clearly , these results depend on the stellar density profile assumed , and rotational effects have been neglected in our simple estimate .
however , it seems clear that the stellar matter is opaque to high - energy neutrinos for small radii , constraining the range of allowed radii for internal shocks .    the case for the reverse shock in successful jets , which break out of the star , seem better .
this however is not due to any intrinsic property of successful jets .
it is merely because for a given total jet energy , successful jets have longer durations and lower luminosities compared to choked jets .
this results in less meson cooling ; for example , we expect 2 events from a choked jet ( of energy @xmath264 ) , while we expect 36 from a successful jet of the same total energy .
however , it should be noted that we expect internal shocks to be important for successful jets .
considering our successful jet with @xmath246 , @xmath240 , and @xmath265 ( so that internal shocks occur at @xmath266 ) , we find that the total muon event number from the internal shock is @xmath3100 times larger than our estimates from the reverse shock  @xcite .
this is due to less meson cooling in the cold jet , and we compare in fig .
[ suppression2 ] the suppression factor for mesons in the cold jet and hot jet .
it is clear that suppression in the cold jet is much weaker .
suppression factor for neutrinos from the reverse shock ( solid lines ) and the internal shock ( dot - dashed lines ) .
we adopt our successful mildly relativisitc jet ( @xmath246 , @xmath240 , @xmath265 ) , and show the suppression factors for pions ( thin - black lines ) and kaons ( thick - green lines ) .
the vertical dotted lines indicate break energies.,width=326 ]    despite being overwhelmed in terms of numbers , reverse shock neutrinos will arrive prior to neutrinos from internal shocks .
mszros and waxman  @xcite studied precursor signals arising from internal shock accelerated protons interacting with photons at @xmath267 , while razzaque , mszros , and waxman  @xcite extended this to include @xmath268 .
neutrinos from reverse shocks will precede even these precursor signals , provided the radius is large enough that the stellar material is thin to neutrinos .
the temporal profile of neutrinos therefore offers the possibility to distinguish locations of particle acceleration for exceptionally close supernovae .
the expected time difference between components is of order a few to tens seconds , which is within the time resolution of neutrino detectors .
we classified jets by whether they break out of the star or not , depending on their durations .
choked jets are bright in neutrinos and dark in gamma rays , while successful jets are bright in neutrinos but may or may not be bright in gamma rays , depending on whether internal shocks are optically thick or not
. a successful yet gamma dark jet may be typical for mildly relativistic jets , where the smaller @xmath5 implies internal shocks typically occur at smaller radii  @xcite .
these could be distinguished from choked jets by the presence of an afterglow like radio signature , for example .    in the current work
we have assumed a type ib progenitor .
however , if the progenitor is a type ic , the situation is slightly different .
having no he core , the core crossing time is shorter , @xmath269 s , for @xmath103 erg s@xmath23 , and choked jets will become rarer .
this implies the relative importance of internal shocks over reverse shocks . for successful ultrarelativistic jets in type ic progenitors
, the neutrino precursor signal can be as strong as those discussed previously  @xcite .
our results will however apply well to type ii progenitors .
the core - collapse supernova rate within 10 mpc is @xmath313 yr@xmath23  @xcite , and therefore , one could test the possibility of particle acceleration in the reverse shocks in the next decade . in particular ,
correlating with the optical signal of the supernova explosions ( not necessarily the jet - like signature )  @xcite or mev neutrino detections  @xcite would greatly help enhance the significance of the signal .
examples of recent close supernovae include sn  2007gr in ngc  1058 ( distance @xmath310 mpc )  @xcite , and a further list can be found in ref .
the actual rate could in fact be higher , as a systematic survey of all nearby galaxies has not been performed to date .
we showed that neutrino emission from reverse shock accelerated protons are strongly suppressed due to meson cooling in the hot jet and hot stellar material surrounding the reverse shock . despite this suppression
, we still expect a fairly promising number of events at icecube , ranging from @xmath242 to @xmath270 for a close supernova of @xmath212 hosting a choked mildly relativistic jet of total isotropic energy @xmath21@xmath271 .
we find that for a supernova of given energy @xmath58 , neutrino emission is favored by high - lorentz factor , low - luminosity , and long - duration jets , for which meson cooling is less severe .
our result sets the scale for neutrino emission from choked jets . for choked jets to provide more neutrinos ,
internal shocks are required to occur in the unshocked region of the jet . however ,
this is geometrically constrained from above by the reverse shock , and from below by the fact that the stellar matter is opaque to high - energy neutrinos for small radii .
regarding neutrino precursor signals which precede grb photons , we conclude that the reverse shock neutrinos likely plays a small role . however
, they will precede those neutrinos from internal shocks , allowing the temporal profile of precursor neutrinos from a close supernova to provide an identification of proton acceleration at the reverse shock .
we are grateful to john beacom for discussions and comments , and peter mszros for reading the script .
we also thank ehud nakar and yudai suwa for early discussions , and hylke koers for email discussions after initial submission .
sa was supported by the sherman fairchild foundation .
we discuss cooling mechanisms in the shocked jet head , which act on the proton ( section iiia , shown in fig .
[ cooling ] ) and mesons ( section iiic ) .
we explicitly denote by primes ( @xmath119 ) the jet head comoving frame .
all numerical values are derived for the proton .
first , the same magnetic field responsible for electron synchrotron radiation also causes energy losses in protons and mesons , given by @xmath272 where @xmath273 is the electron mass , and @xmath274 and @xmath275 are the particle mass and energy .
the numerical value is derived for the proton . due to the strong dependency on the mass , @xmath236 ,
synchrotron cooling is most effective for the pion meson and least effective for the proton .
second , protons and mesons lose energy by inverse compton ( ic ) scattering .
the most natural target photons are due to radiation cooling of ultrarelativistic electrons .
although a full treatment using a black - body spectrum is simple , a simplistic approximation using eq .
( [ photons ] ) and eq .
( [ photons2 ] ) is sufficient to demonstrate that ic cooling is subdominant . below a transition energy @xmath276 ( for the proton ) , inverse compton scattering proceeds in the thomson regime , and the cooling time scale is given @xmath277 again , note the strong dependency on the particle mass . above the transition energy ,
inverse compton scattering proceeds in the klein - nishina limit , where the cooling time scale becomes @xmath278    third , we consider cooling by the bethe - heitler ( bh ) interaction ( @xmath279 ) , which is important for protons .
our treatment is the same as that described in ref .
the bh cross section rises logarithmically with energy , @xmath280\}$ ]  @xcite .
the @xmath281 produced are at rest in the centre of mass frame of the proton and photon , and so the energy lost by the proton in each interaction is given by @xmath282 , where @xmath283 is the lorentz factor of the centre of inertia as seen in the comoving frame .
the energy loss rate is @xmath284 , and thus the cooling time is @xmath285 .
lastly we discuss cooling due to collisions with protons . assuming @xmath286 of the proton energy is lost in each collisions with a proton , the cooling time scale is @xmath287 for the proton , where @xmath163  @xcite .
analogous equations hold for mesons , e.g.  eq .
( [ mesonhadronic ] ) for the pion .
s.  ando , j.  f.  beacom and h.  yuksel , phys .
* 95 * , 171101 ( 2005 ) [ arxiv : astro - ph/0503321 ] . s.  ando , jcap * 0406 * , 007 ( 2004 ) [ arxiv : astro - ph/0405411 ] .
e.  nakar , phys .
rept .   *
442 * , 166 ( 2007 ) [ arxiv : astro - ph/0701748 ] .
m.  a.  aloy , e.  mueller , j.  m.  ibanez , j.  m.  marti and a.  macfadyen , astrophys .
j.   * 531 * , l119 ( 2000 ) [ arxiv : astro - ph/9911098 ] .
a.  burrows , l.  dessart , e.  livne , c.  d.  ott and j.  murphy , arxiv : astro - ph/0702539 .
t.  takiwaki , k.  kotake , s.  nagataki and k.  sato , astrophys .
j.   * 616 * , 1086 ( 2004 ) [ arxiv : astro - ph/0408388 ] .
a.  achterberg , y.  a.  gallant , j.  g.  kirk and a.  w.  guthmann , mon .  not .
* 328 * , 393 ( 2001 ) [ arxiv : astro - ph/0107530 ] .


superconducting qubits are solid - state artificial atoms , comprising lithographically defined josephson tunnel junctions and superconducting interconnects . when cooled to milli - kelvin temperatures , these qubits exhibit quantized states of charge , flux , or junction phase depending on the circuit design parameters  @xcite . associated with these quantized states
is a spectrum of energy levels , tunable via an external control parameter , e.g. , an applied electric or magnetic field .
although generally only the lowest two energy eigenstates are utilized for quantum information science applications , the energy spectrum indeed extends to higher - energy levels corresponding to higher - excited states of the circuit .
the separation between pairs of energy levels typically falls in the radio frequency and microwave regimes , and resonantly driving the artificial atom with a harmonic field can couple and induce quantum - state transitions .    due to their relatively large size , superconducting artificial atoms can be strongly coupled to their external control fields .
it is this feature , along with their quantum coherence , that we utilize in the present article .
a large - amplitude harmonic control field can drive an artificial atom throughout its energy - level spectrum . when driven through an avoided level crossing
, a landau - zener - stckelberg ( lzs ) transition occurs .
this is a coherent process akin to a beamsplitter for photons , taking an input state of the atom and outputting a superposition of states .
repeated passages through an avoided crossing act as an atom interferometer , causing the atomic superposition state to interfere quantum mechanically with itself .
since the weighting of the superposition state depends sensitively on the size of the avoided crossing and the velocity ( change in relative energy between levels per unit time ) with which it is traversed , the quantum interference reflects the energy spectrum of the artificial atom . in turn
, the quantum interference can be leveraged to facilitate non - adiabatic quantum control .
we begin this article with an introduction to the superconducting persistent - current qubit  @xcite and an overview of landau - zener - stckelberg ( lzs ) transitions .
we then present three experimental works that utilize lzs transitions in a strongly - driven persistent - current qubit .
the first is mach - zehnder - type interferometry between repeated lzs transitions  @xcite , for which we observed quantum interference fringes in n - photon transition rates  @xcite , with n=1  50 .
the second is microwave - induced cooling  @xcite , by which we achieved effective qubit temperatures less than 3 mk , a factor 10 - 100 times lower than the environmental temperature .
the third is amplitude spectroscopy  @xcite , a spectroscopy technique that monitors the system response to amplitude rather than frequency .
amplitude spectroscopy allowed us to probe the energy spectra of our artificial atom from 0.01 - 120 ghz , while driving it at a fixed frequency near 0.2 ghz .
finally , we conclude by considering the application of lzs transitions to quantum information science and technology .
superconducting artificial atoms exhibit a high degree of quantum coherence , and there have been numerous proposals and demonstrations of quantum phenomena in these systems , many derived from the fields of atomic physics and quantum optics .
a few examples include : coherent superpositions of macroscopic states  @xcite , rabi oscillations  @xcite , landau - zener transitions  @xcite , stckelberg oscillations  @xcite , microwave cooling  @xcite , electromagnetically induced transparency  @xcite , geometrical phase  @xcite , and cavity quantum electrodynamics  @xcite .
significant progress has also been made toward their application to quantum information science  @xcite , including state initialization  @xcite , tunable  @xcite and long - distance  @xcite coupling , quantum control  @xcite , quantum state  @xcite and process  @xcite tomography , and measurement  @xcite .
for a recent review of superconducting qubits , see ref .  .        here
, we use a superconducting persistent - current qubit to implement an artificial atom  @xcite .
the persistent - current qubit is a superconducting loop interrupted by three josephson junctions ( fig .
[ fig:1]a ) .
when biased with a static magnetic flux @xmath0 , where @xmath1 is the superconducting flux quantum , the system assumes a double - well potential profile ( fig .
[ fig:1]b ) .
the diabatic ground state of the left ( right ) well corresponds to a persistent current @xmath2 with clockwise ( counterclockwise ) circulation .
these two diabatic energy levels ( fig .
[ fig:2]a ) have a separation @xmath3 linear in the flux detuning @xmath4 .
higher - excited states of the double - well potential ( see fig .  [ fig:5]c ) will be considered in sections  [ sec : cooling ] and  [ sec : as ] .
the two - level hamiltonian for the lowest two diabatic states is shown in eq .
[ eq : hamiltonian ] . at detuning @xmath5 ,
the double - well potential is symmetric and the diabatic - state energies are degenerate . at this `` degeneracy point , ''
resonant tunneling between the diabatic states opens an avoided level crossing of energy @xmath6 . here
, the qubit states are @xmath7 eigenstates , corresponding to symmetric and anti - symmetric combinations of diabatic circulating - current states .
detuning the flux away from this point tilts the double well , allowing us to tune the eigenstates and eigenenergies of the artificial atom .
far from the degeneracy point the qubit states are approximately @xmath8 eigenstates , the diabatic states with well - defined circulating current .
the qubit is read out using a hysteretic dc squid ( superconducting quantum interference device ) , a sensitive magnetometer that can distinguish the flux generated by circulating current states .
in addition to the static flux biases , the artificial atom is controlled and read out using the pulses illustrated in fig .
[ fig:1]c .
as we describe below , the qubit is first prepared in its ground state using a harmonic cooling pulse with amplitude @xmath9 and frequency @xmath10 .
quantum - state transitions are then driven using a harmonic rf pulse with amplitude @xmath11 and frequency @xmath12 .
these fields are mutually coupled to the qubit through a small antenna .
this is followed by a squid readout current pulse using the `` sample and hold '' technique  @xcite .
if the sample current exceeds the squid switching current , a voltage pulse will appear at the output during the hold phase .
threshold detection looks for the presence or absence of a squid voltage , and this constitutes a digital measurement of the qubit state .
alternatively , although not used in these experiments , we have incorporated the squid into a resonant circuit and realized qubit readout via the shift in resonance frequency and phase for both the linear and non - linear resonance regimes  @xcite .
the `` qubit step '' is shown in fig .
[ fig:1]d as a function of the squid sample current and the flux detuning .
the diabatic states @xmath13 and @xmath14 correspond to different levels of sample current ( dashed lines ) located symmetrically about the degeneracy point @xmath15=0 .
this plot constitutes a cumulative switching current distribution of the squid for each @xmath15 value .
additionally , a 3-@xmath16s rf pulse at 1.2 ghz is applied to the qubit , and resonant transitions can be observed as fingers extending down ( up ) from state @xmath13 ( @xmath14 ) when @xmath171.2 ghz becomes resonant with the energy - level separation .
a best - estimator ( dashed - dotted line ) can be determined to provide the best statistical discrimination between states @xmath13 and @xmath14 .
the resulting qubit step with its saturated n - photon resonances along the best estimator line is shown in fig .
[ fig:1]e .
the device used in this work was fabricated at lincoln laboratory using a fully - planarized niobium trilayer process and optical lithography .
it has a critical current density @xmath18 , and the characteristic josephson and charging energies are @xmath19 @xmath20 respectively .
the ratio of the qubit josephson junction ( jj ) areas is @xmath21 , and @xmath22 .
although dependent on the flux detuning , the approximate values for coherence times are : interwell relaxation time @xmath23 , intrawell relaxation time @xmath24 , homogeneous transverse decay time @xmath25 , inhomogeneous transverse decay time @xmath26 .
the experiments were performed in a dilution refrigerator at a base temperature of 20 mk .
the device was magnetically shielded , and all electrical lines were carefully filtered and attenuated to reduce noise ( see ref .   for details ) .
landau - zener - stckelberg ( lzs ) transitions occur when a quantum system is driven through an energy - level avoided crossing  @xcite .
the resulting quantum dynamics of the lzs mechanism   in driven systems  @xcite have been developed  @xcite within a two - level coherent scattering formalism  @xcite with potential application to quantum information science and technology  @xcite .
we model a lzs transition at a single crossing using the two - level hamiltonian @xmath27 in which @xmath28 is the energy difference between diabatic - state energy levels ( dashed lines in fig .
[ fig:2]a ) , and @xmath6 is the size of the avoided crossing and corresponds to the coupling strength between diabatic states @xmath29 and @xmath30 .
the system is prepared in state @xmath31 at a static value @xmath32 far from the avoided crossing ( blue dot , fig .
[ fig:2]a ) and driven longitudinally from that point by making @xmath33 time dependent .
we distinguish here between `` weak driving , '' in which the resulting excursion is not large enough to reach the avoided crossing , and `` strong driving '' for which the avoided crossing is traversed . under strong - driving conditions , the asymptotic probability @xmath34 of a transition between diabatic states ,
@xmath35 is governed by the relative - energy sweep rate @xmath36 @xmath37 evaluated at the time @xmath38 at which the system is swept through the avoided crossing . in the original lzs formulation , the system was driven with a fixed sweep rate , whereas we will later consider a harmonically driven system .
we also note that we have elected to monitor the probability @xmath34 of a transition between the diabatic states @xmath39 and @xmath31 , because our readout detector is sensitive to changes in diabatic state .
this can be written with no loss of generality in terms of the probability @xmath40 of a transition between the system eigenstates .
there are two strong - driving limits characterized by the relative sizes of the sweep velocity @xmath41 and the avoided crossing @xmath6 . in the adiabatic limit ( fig .
[ fig:2]a ) , the sweep velocity is small ( @xmath42 ) and the probability of a transition from state @xmath31 to state @xmath39 approaches unity , @xmath43 . in this case , the system dynamics are slow enough that the system adiabatically follows the ground eigenstate through the avoided crossing . in the nonadiabatic limit ( fig .
[ fig:2]b ) , the sweep velocity is large ( @xmath44 ) and the probability of a transition approaches zero @xmath45 . in this case , the dynamics are too fast for the system to follow ; the system remains in diabatic state @xmath31 and , thereby , effectively jumps the energy gap at the avoided crossing .
more generally , a superposition state @xmath46 results from an excursion through the avoided crossing , as illustrated in fig .
[ fig:2]c .
following an idea discussed by shytov _
et al . _
@xcite , the lzs transition acts as a beamsplitter for the atomic state .
the amplitudes @xmath47 and @xmath48 are determined by a unitary transformation @xmath49 , effectively a @xmath50 `` beamsplitter matrix '' comprising complex reflection @xmath51 and transmission @xmath52 coefficients related to the adiabaticity parameter  @xcite @xmath53 present in eq .
[ eq : p_lzs ] such that @xmath54 , @xmath55 , and @xmath56 .
note that we have defined @xmath57 and @xmath58 from the perspective of a beamsplitter , which `` transmits '' ( `` reflects '' ) an input state to the same ( opposite ) diabatic state , respectively .
after the transition , a relative phase @xmath59 accrues due to the energy difference between the states @xmath39 and @xmath31 . if the drive @xmath33 then returns the system through the avoided crossing a second time , the atomic state collides and quantum mechanically interferes with itself during the second lzs transition ( fig .
[ fig:2]d ) .
the cumulative result is an atom - state interferometer whose output state depends on the lzs transition amplitudes and the interference phase .
an analogy can be made to an optical mach - zehnder interferometer : the atomic states play the role of the photon modes , the lzs transitions play the role of the photon beamsplitters , and the energy - level splitting , which determines the interference phase , plays the role of optical path length difference .
the quantum interference is robust provided the evolution time of the state through the interferometer is short compared with the atom s coherence times .
in addition to superconducting artificial atoms  @xcite , this concept is generally applicable to other solid - state artificial atoms  @xcite and generalized spin systems  @xcite ( e.g. , molecular magnets  @xcite , natural atoms  @xcite , and molecules  @xcite ) that exhibit avoided level crossings , and it is extensible to multiple energy levels as we demonstrate in section  [ fig:5 ] .
the structure of the n - photon spectroscopy peaks seen in fig .
[ fig:1]e consists of regularly - spaced resonances positioned according to the condition @xmath171.2 ghz being resonant with the energy level separation .
notably , however , for this particular value of driving amplitude , the n=1 , 3 , and 6 photon peaks are missing . as we describe in this section ,
the presence and absence of these peaks arise from mach - zehnder - type quantum interference at a level crossing .
the interference phase is tunable via the driving amplitude , leading to a `` bessel staircase '' interference pattern in the observed spectroscopy .
the interference oscillations are known as stckelberg oscillations  @xcite , and they have been observed in both artificial  @xcite and natural  @xcite atomic systems .        in a conventional mach - zehnder interferometer
, an optical signal is passed through two beamsplitters .
the first beamsplitter coherently divides the signal into two output paths , which may have different effective optical path lengths .
these paths are then recombined at the second beamsplitter , where the signal waves interfere and exit the interferometer through the two output ports .
an intensity measurement at either output port exhibits interference fringes depending on the relative path length .    here , we instead utilize lzs transitions at level crossings as beamsplitters for the atomic state  @xcite .
we drive the persistent - current qubit with a harmonic driving field of the form @xmath60 with @xmath61 the driving frequency and @xmath62 the field amplitude ( in units of energy ) , which is proportional to the rf voltage at the source . as illustrated in fig .
[ fig:2]d , the qubit state undergoes two lzs transitions during one period of the driving field . the first lzs transition at time @xmath63 splits the qubit into a superposition of excited and ground states .
a relative phase @xmath59 accumulates @xmath64 until the second lzs transition at time @xmath65 , at which point the qubit state collides with itself and interferes .
interference fringes appear at half - integer and integer values of @xmath66 , which is tunable via the drive amplitude .
it is clear , however , that a second phase must also play a role in this problem , since the qubit state continues to evolve for the remainder of the driving period .
it is physically meaningful to consider the total phase @xmath67 accumulated over a single period : @xmath68 which is _ independent _ of the driving amplitude . over many periods of the driving field ,
the cascaded pairs of lzs transitions ( cascaded interferometers ) will constructively interfere provided @xmath69 .
one can view this as a `` time - domain '' formulation for the familiar n - photon resonance condition @xmath70 where n is the number of photons involved in the transition .
it is only when the product @xmath71 equals the energy separation @xmath72 that the cascaded lzs transitions lead to a net buildup of state population and , as a function of @xmath59 , the observed interference fringes .
these oscillations are related to photoassisted transport  @xcite and rabi oscillations  @xcite in the multiphoton regime .
mach - zehnder - type interference in the discrete - resonance limit driven towards saturation is shown in fig .
[ fig:3]a for driving frequency @xmath73 mhz .
this frequency is larger than the resonance linewidth ( fig .
[ fig:3]b ) , and so individual n - photon resonances can be resolved .
there are two main features observable in this plot .
the first is the presence of equally - spaced n - photon transitions as a function of flux detuning , symmetrically located about the qubit step at @xmath5 . as one might expect , the onset of the higher - photon transitions requires larger driving amplitudes . remarkably , we observe up to 50-photon transitions in this scan .
the second main feature is that for each n - photon resonance , the spectroscopy appears and disappears as a function of amplitude , which sweeps the interference phase @xmath59 .
for example , 14 oscillation lobes are visible for the 1-photon transition . a vertical slice of the spectroscopy ( dashed white line in fig .  [ fig:3]a )
is plotted in fig .
[ fig:3]c ; at this particular amplitude , one can see the enhancement and suppression of the spectroscopy peaks as was observed in fig .
[ fig:1]e .
the mach - zehnder interference for an n - photon transition yields a modified amplitude - dependent matrix element  @xcite @xmath74 where @xmath75 is the nth - order bessel - function , and its argument @xmath76 is the dimensionless driving amplitude .
intuitively , the bessel - function dependence arises , because the qubit is driven harmonically through energy levels that are linear in flux detuning and , as a result , the interference phase accumulates with a harmonic time dependence . the transition rate and , therefore , the population transfer , at each n - photon resonance is related to the matrix element squared , @xmath77 .
at specific amplitudes , despite driving the system resonantly , no net transition occurs due to a complete destructive mach - zehnder interference corresponding to the zeros of @xmath75 ; this is the coherent destruction of tunneling  @xcite condition for driven n - photon transitions  @xcite .
more generally , as a function of the driving amplitude , we have observed the continuous evolution between regions of enhanced and suppressed tunneling as dictated by @xmath75 over several bessel - function lobes ( fig .
[ fig:3]a ) .    in the discrete resonance limit ,
n - photon resonances are distinguishable , because the coherence time of the qubit is sufficiently longer than the driving period . by reducing the frequency to @xmath78 mhz ( fig .
[ fig:3]a , right panel ) , we effectively made the drive period comparable with the linewidth , @xmath79 ( fig .  [ fig:3]b ) . in this spectroscopy plot ,
the individual n - photon resonances are no longer resolvable .
however , the mach - zehnder interference fringes ( vertical slice , fig .  [ fig:3]c ) are clearly visible , because the qubit remains coherent during the critical fraction of the drive period during which the phase @xmath59 accumulates and the mach - zehnder interference occurs . intuitively , provided @xmath80 , the output populations of each quantum interferometer are then preserved ( frozen ) until the subsequent interferometer is reached in the following period ; therefore , although the resonance condition is lost , the mach - zehnder quantum interference remains .
this behavior can be contrasted with driven rabi oscillations , for which there would be no signature of coherence in the limit @xmath79 .
the previous discussion involved driving transitions in the lowest two energy levels in the double - well potential of our artificial atom ( fig .
[ fig:4]a ) , which constitute the two - level qubit subsystem of a more complex energy level diagram ( fig .  [ fig:5]c ) .
when higher - excited states are accessed , the driven system behavior can be markedly different from the population saturation observed when only two levels are involved .
for example , at least three levels are required to achieve population inversion , and such a multi - level artificial atom coupled to a microwave cavity has been used to demonstrate masing ( microwave lasing )  @xcite . in that work ,
josephson quasi - particle states were driven to achieve inversion .
alternatively , population inversion can be established by accessing an ancilliary excited state via direct or lzs transitions .
this will be briefly discussed in the next section .    here , by reversing the cycle that leads to population inversion ,
we show that one can pump population from the qubit excited state @xmath81 to the qubit ground state @xmath82 ( fig .
[ fig:4]a ) via an ancillary energy level @xmath83 . in the case
where the population in @xmath81 results from thermal excitation , the transfer of population to @xmath82 effectively cools the qubit by lowering its effective temperature .
this kind of active cooling represents a means to initialize and reset qubits with high fidelity , key elements for quantum information science and technology .
alternatively , the pumping mechanism can be used to refrigerate environmental degrees of freedom  @xcite , or to cool neighboring quantum systems  @xcite .
more explicitly , for a qubit in equilibrium with its environment , the population in @xmath84 that is thermally excited from @xmath82 follows the boltzmann relation @xmath85 ,   \label{boltzmann}\ ] ] where @xmath86 are the qubit populations for energy levels @xmath87 , @xmath88 , @xmath89 is the boltzmann constant , and @xmath90 is the bath temperature . to cool the qubit subsystem below @xmath90 ,
a microwave magnetic flux of amplitude @xmath91 and frequency @xmath12 targets the @xmath92 transition , driving the state-@xmath84 thermal population to state @xmath83 , from which it quickly relaxes to the ground state @xmath82 .
efficient cooling only occurs when the driving - induced population transfer to @xmath82 is faster than the thermal repopulation of @xmath84 .
the hierarchy of relaxation and absorption rates required , @xmath93 , is achieved in our system owing to a relatively weak tunneling between wells , which inhibits the interwell relaxation and absorption processes @xmath94 and @xmath95 , compared with the relatively strong intrawell relaxation process @xmath96 .
figure [ fig:4]b shows the qubit step at @xmath90 = 150 mk in equilibrium with the bath ( top ) and after a 3-@xmath16s cooling pulse at 5 mhz ( bottom ) . under equilibrium conditions ,
the average level populations exhibit a thermally - broadened qubit step about @xmath97 , with a width proportional to @xmath90 .
the presence of microwave excitation targeting the @xmath92 transition , followed by relaxation , acts to increase the ground - state population and , thereby , sharpens the qubit step .
cooling can thus be quantified in terms of an effective temperature @xmath98 , a signature that is evident from the narrowing of the qubit steps in fig .
[ fig:4]b after cooling .
more precisely , using the notation from fig .
[ fig:2 ] , the effective qubit temperature is obtained by fitting an effective temperature that would have been required in equilibrium to achieve the observed qubit population @xmath99 , @xmath100.\ ] ]    universal cooling ( cooling that is independent of flux detuning ) occurs near an optimal driving amplitude @xmath101 ( fig .  [ fig:4]c ) .
this is demonstrated in fig .
[ fig:4]d where we present the @xmath82 state population @xmath102 measured as a function of the microwave amplitude @xmath91 and flux detuning @xmath103 for @xmath104 . cooling and
the diamond feature can be understood in terms of the energy level diagram ( fig .
[ fig:5]c ) . as the amplitude of the microwave pulse
is increased from @xmath105 , population transfer first occurs when the @xmath106 avoided crossing is reached , _
i.e. _ @xmath107 ; this defines the front side of the observed diamond , symmetric about the qubit step ( see also fig .  [
fig:3]a ) . for amplitudes
@xmath108 , the @xmath109 ( @xmath110 ) side avoided crossing dominates the dynamics , resulting in a second pair of thresholds @xmath111 , which define the back side of the diamond . in the region outside of the diamond s backside ,
the qubit is cooled .
as the diamond narrows to the point @xmath112 , the narrowest qubit step is observed .
this is the universal cooling condition : only one of the two side avoided crossings ( @xmath109 or @xmath110 ) is reached and , thereby , strong transitions with relaxation to the ground state result for all @xmath103 . in contrast , for @xmath113 , both side avoided crossings ( @xmath109 and @xmath110 ) are reached simultaneously for @xmath114 , leading once again to a large population transfer between @xmath82 and @xmath84 , and opening the second diamond feature ( see fig .  [ fig:5 ] ) .
the cooling exhibits a rich structure as a function of driving frequency and detuning , resulting from the manner in which state @xmath83 is accessed ( fig . [
fig:4]c ) .
transitions occur via a ( multiphoton ) resonant or adiabatic passage process when the driving frequency is high or low enough , respectively  @xcite . at high frequencies ( 800 and 400 mhz in fig .
[ fig:4]e ) well - resolved resonances of n - photon transitions are observed and cooling is thus maximized near resonances . at intermediate frequencies ( 400 and 200 mhz ) ,
mach - zehnder interference at the side avoided crossing @xmath115 becomes more prominent and modulates the intensity of the n - photon resonances  @xcite .
below @xmath116 , individual resonances are no longer discernible , but as in fig .
[ fig:3]c , the modulation envelope persists  @xcite . at the lowest frequencies ( @xmath117 ) ,
state @xmath83 is reached via adiabatic passage through the @xmath115 crossing ( fig .
[ fig:4]c ) , and the population transfer and cooling become conveniently independent of detuning ( see @xmath118 in fig .  [
fig:4]e ) . as shown in fig .
[ fig:4]e , we achieve an effective qubit temperature @xmath119 = 3 mk , even for @xmath90 = 300 mk . in our qubit ,
our determination of @xmath119 was limited primarily by decoherence ( linewidth ) , which limited the resolution with which we could distinguish the states @xmath82 and @xmath84 near degeneracy .
nonetheless , we can estimate the ideally resolvable cooling factor @xmath120 for this type of cooling process using eq .
[ boltzmann ] , @xmath121 where @xmath122  ghz is the energy separation where the relaxation @xmath96 occurs and @xmath123  ghz for our qubit , yielding a cooling factor @xmath124 . for a bath temperature
@xmath125  mk , this would correspond to an effective temperature @xmath126  @xmath16k in our qubit . cooling a qubit in equilibrium with the bath
requires a characteristic cooling time . in turn
, a cooled qubit will thermalize to the environmental bath temperature over a characteristic equilibration time .
the relationship between these two times determines if it is possible to drive the qubit while it is still cold .
we found in this qubit that equilibration times are at least one order of magnitude larger than cooling times at @xmath127 mk and up to three order of magnitude larger at @xmath128 mk  @xcite .
this allowed us ample time to drive the qubit after cooling it .
the implementation of an active cooling pulse prior to a generic driving pulse is highly advantageous . on the one hand , it sensibly shortens measurement times , enabling us to acquire data at repetition rates that far exceed the intrinsic equilibration rate due to interwell relaxation after each measurement trial . by adopting active cooling
, we gained a factor 50 in data acquisition speed , limited by the bandwidth of our readout circuit . on the other hand , by analogy to the cooling of trapped ions and atoms , active cooling
greatly reduces thermal smearing , allowing us to analyze features in the data that would have been hidden otherwise .
this type of active cooling protocol was required to obtain the detailed , clean data in fig .
[ fig:3 ] , where we could clearly resolve resonances separated by only 270 mhz ( 13 mk ) .
in fact , the necessity of active cooling becomes even more evident in the next section where , without the cooling pulse , the observed level of detail could not have been resolved over such a large parameter space in practical acquisition times .
frequency - dependent absorption and emission spectroscopy has long played a fundamental role in the characterization of quantum systems .
the development of coherent microwave ( maser ) and optical ( laser ) sources , high - intensity radiation with tunable , narrow spectral linewidth , has further enabled targeted absorption spectroscopy of atoms and molecules with high frequency resolution  @xcite . however , the application of broadband frequency spectroscopy is not universally straightforward .
this is particularly relevant for artificial atoms engineered from solid - state materials which , when cooled to cryogenic temperatures , assume quantized energy levels that extend into microwave and millimeter wave regimes ( 10 - 300 ghz ) .
although certainly not an impossible task , a broadband frequency - based spectroscopic study of these atoms in excess of around 50 ghz becomes extremely challenging and expensive to implement due to numerous frequency - dependent effects ( _ e.g. _ , frequency dispersion and the requisite tolerances to control impedance ) .
amplitude spectroscopy is a technique that allows broadband spectroscopic characterization of a quantum system . with amplitude spectroscopy
, spectroscopic information is obtained from the system response to driving - field amplitude at a _ fixed _ frequency .
the resulting spectroscopic interference patterns , `` spectroscopy diamonds , '' are mediated by multilevel lzs transitions and mach - zehnder - type interferometry , and they serve as a fingerprint of the artificial atom s multilevel energy spectrum ( figs .  [ fig:5]a and  [ fig:5]b ) . the energy spectrum
is then determined by analyzing the atomic fingerprint . in this way , the amplitude spectroscopy technique complements frequency spectroscopy : although a less direct approach , it allows one to probe the energy level structure of a quantum system over extraordinarily large ( even practically prohibitive ) bandwidths by circumventing many of the challenges associated with a frequency - based approach .    in general ,
the spectroscopy diamonds arise due to an interplay between the static flux detuning @xmath129 and driving amplitude @xmath11 .
as described in sections  [ sec : lzs ] and  [ sec : mz ] , transitions occur when an avoided crossing is reached for a particular set of values ( @xmath129,@xmath11 ) .
for example , at a flux detuning @xmath130 , the diamond boundaries occur at @xmath131 ( fig .
[ fig:5 ] a and c ) .
the mach - zehnder interference due to a phase accumulation @xmath132 at a given avoided crossing @xmath133 can be modulated by varying both @xmath129 and @xmath11 .
there are two important contributions to the diamond spectroscopy patterns : lzs - mediated transitions and intrawell relaxation . in fig .
[ fig:5]d , we show a subset of the diamond interference pattern in fig .
[ fig:5]a .
arrows indicate lines of constant phase accumulation @xmath134 in ( @xmath129,@xmath11)-space that leads to lzs transitions at each of the three listed avoided crossings , @xmath106 , @xmath109 , and @xmath110 . where these lines cross , competition ( coordination ) between avoided crossings act to suppress ( enhance ) the net transition rate between pairs of energy levels .
the arrangement of these crossing lines leads to the checkerboard patterns observed inside and outside the diamonds .
the second contribution , intrawell relaxation , provides another means to connect energy levels and results in both cooling and population inversion . in section  [ sec : cooling ] ,
[ fig:4]c , @xmath109 mediated the transition @xmath92 , and intrawell relaxation then mediated the transition @xmath96 ; the net result was cooling , since the flux @xmath129 was positive , making @xmath82 the ground state .
however , in fig .
[ fig:5]d @xmath129 is negative .
furthermore , in the upper - right corner of fig .
[ fig:5]d , both crossings @xmath109 and @xmath110 are accessed . in the bright red regions
, @xmath109 still tends to cause transitions @xmath92 , and relaxation puts that population in state @xmath82 .
however , the interference condition at @xmath110 on the other side of the energy level diagram ( fig .
[ fig:5]c ) tends to keep the population in @xmath82 .
therefore , population builds in state @xmath82 , the first excited state for negative flux detuning , resulting in strong population inversion . varying the interference conditions at @xmath109 and @xmath110 by changing ( @xmath129,@xmath11 ) causes the observed modulation between population inversion and cooling .
we have developed several techniques for extracting information about the energy levels from the spectroscopy interference patterns  @xcite .
the key metrics are the positions of the avoided crossings in flux , the values of the splittings @xmath133 , and the slopes of the energy levels . with this information
, one can reconstruct a large portion of the energy level diagram .
the positions of the avoided crossings can be determined precisely from the diamond boundaries , because the onset of each diamond is associated with lzs transitions at a particular level crossing .
the splitting of each avoided crossing can be obtained essentially by fitting the lzs formula in eq .
[ eq : p_lzs ] to the mach - zehnder interference patterns .
alternatively , one can study the dynamical population transfer between states using the pulsed , short - time implementation of amplitude spectroscopy  @xcite .
the energy - level slopes can be determined by two means .
the first is by relating the separation between mach - zehnder interference nodes to the expected phase accumulation @xmath59 , which depends sensitively on the energy - level slope .
alternatively , we show with rudner _
_ in ref .   that the two - dimensional fourier transform of the diamond patterns yields a family of _ one - dimensional _ sinusoids in fourier space ; the periods of these sinusoids are related to the energy - level slopes .
the intuition is that the fourier transform inverts the energy domain of the spectroscopy to the time domain ( scaled by @xmath135 ) .
this means that the sinusoids in fourier space image the time - dependence of the quantum phase of the qubit , which , in turn , depends sensitively on the energy - level slopes .    using amplitude spectroscopy
, we were able to scan the energy level diagram continuously beyond the fourth energy - level avoided crossing ( @xmath136 , @xmath137 ) with splitting values ranging from @xmath138  ghz to @xmath139  ghz  @xcite .
the equivalent information obtained using frequency spectroscopy would have required scanning frequencies from 0.01 ghz out to beyond 100 ghz ( in this device , each avoided crossing is separated by approximately 25 ghz , and so each subsequent crossing in fig .
[ fig:5 ] is raised an additional 25 ghz above the ground state ) .
remarkably , with amplitude spectroscopy , the entire scan performed in fig .
[ fig:5]a was performed at a fixed frequency @xmath140  mhz .
the scan in fig .
[ fig:5]b shows the amplitude spectroscopy of the same system for a fixed frequency @xmath141  mhz , clearly in the discrete resonance limit .
the resonance condition adds another constraint , making a more complex and interesting checkerboard pattern .
strongly driving a superconducting artificial atom through an avoided level crossing results in a landau - zener - stckelberg transition , which , in general , creates a superposition of atomic states whose weighting depends on the size of the avoided crossing and the velocity with which it is traversed . in this sense , as we discussed in section  [ sec : lzs ] , the lzs mechanism acts as a beamsplitter for artificial atoms .    in section  [ sec : mz ] , we described how harmonically driving the system cascades two lzs transitions per driving period , resulting in an atomic analog to a mach - zehnder interferometer .
the relative phase acquired between lzs transitions is the interference phase of the interferometer , and it is tunable by the driving amplitude .
the buildup of population over many driving periods exhibits stckelberg oscillations as a function of the driving amplitude ( interference phase ) due to the cascaded mach - zehnder - type interference effect . to observe these oscillations ,
the coherence time must only be longer than the small portion of the drive period during which the interference phase accrues , and the relaxation time must be long enough to maintain the population until readout .    the total phase accumulated over one period , in contrast , is amplitude independent . for coherence times
longer than the drive period , cascaded interferometers yield a net population change when this round - trip phase accumulation is @xmath142 per driving period , a condition which can be viewed as the `` time - domain '' counterpart to the @xmath143-photon resonance condition @xmath144 . by making the drive period commensurate with the coherence time
, we showed that we could still observe the stckelberg oscillations , even though the discrete resonances were no longer discernable .
we utilized strong driving and the lzs mechanism with higher - energy levels to achieve both cooling and population inversion in our artificial atom . in section
[ sec : cooling ] , we described using a microwave pumping scheme to cool the atomic degrees of freedom a factor 10 - 100 times colder than the ambient dilution refrigerator temperature .
the scheme involved pumping unwanted thermal population to an ancillary excited state , from which it relaxed to the ground state . in section  [ sec : as ] , we showed that by reversing the order , we could pump population through an ancillary state to achieve inversion .    the energy level structure can be probed over extraordinarily large bandwidth using the amplitude spectroscopy approach presented in section  [ sec : as ] . since the lzs mechanism and mach - zehnder interference are sensitive to the defining features in the energy level diagram ( energy band slopes , level splittings , and their locations ) , the interference patterns that result from sweeping the amplitude are a `` fingerprint '' of the artificial atom s energy spectrum . using amplitude spectroscopy at a fixed driving frequency of only 185 mhz , we could access continuously multiple energy levels from about 10 mhz out to beyond 120 ghz .
large - amplitude driving and the lzs mechanism have application to quantum information science and technology .
active cooling has utility in state initialization and refreshing ancilla qubits in quantum error correction protocols .
amplitude spectroscopy provides a means to ascertain over large bandwidth the energy level structure of a qubit system beyond the lowest two levels .
the mach - zehnder - type interference can facilitate nonadiabatic control schemes , in which the quantum interference at an avoided crossing is used to achieve state transitions that approach the intrinsic coupling rate @xmath6 . in cold atoms and molecules ,
this kind of non - adiabatic control has been used to drive transitions that would otherwise be challenging to realize in a weak - driving limit .
we thank t.p .
orlando , d.m .
berns , l.s .
levitov , m.s .
rudner , a.v .
shytov , and k.k .
berggren with whom we collaborated closely on this work ; j. bylander , s. gustavsson , b. turek , and j. sage for carefully reading the manuscript ; v. bolkhovsky , g. fitch , t.r jordan , e. macedo , p. murphy , k. parrillo , r. slattery , and t. weir at mit lincoln laboratory for technical assistance .
this work was supported by afosr and lps ( f49620 - 01 - 1 - 0457 ) under the durint program , and by the u.s .
the work at lincoln laboratory was sponsored by the us dod under air force contract no .
fa8721 - 05-c-0002 .                        c.  h. van der wal , a. c. j. ter haar , f. k. wilhelm , r. n. schouten , c. j. p.
m. harmans , t. p. orlando , s. lloyd , j. e. mooij , quantum superposition of macroscopic persistent - current states , _ science _ * 290 * , 773 - 777 ( 2000 ) .
a.  j. berkley , h. xu , r. c. ramos , m. a. gubrud , f. w. strauch , p. r. johnson , j. r. anderson , a. j. dragt , c. j. lobb , f. c. wellstood , entangled macroscopic quantum states in two superconducting qubits , _ science _ * 300 * , 1548 - 1550 ( 2003 ) .
a. izmalkov , m. grajcar , e. ilichev , n. oukhanski , th .
wagner , h .-
meyer , w. krech , m. h. s. amin , a. maassen van den brink and a. m. zagoskin , observation of macroscopic landau - zener transitions in a superconducting device , _ europhys .
_ 65 , 844 - 849 ( 2004 ) .
a. izmalkov , s.j.w .
van der ploeg , s.n .
shevchenko , m. grajcar , e. ilichev , u. hbner , a.n .
omelyanchouk , and h .- g.meyer , consistency of ground state and spectroscopic measurements on flux qubits , _ phys .
_ 101 , 017003 ( 2008 ) .
m. grajcar , s.h.w .
van der ploeg , a. izmalkov , e. ilichev , h .-
meyer , a. fedorov , a. shnirman , and g. schn , sisyphus cooling and amplification by a superconducting qubit , _ nature phys .
_ 4 , 612 - 616 ( 2008 ) .
murali , z. dutton , w.d .
oliver , d.s .
crankshaw , and t.p .
orlando , probing decoherence with electromagnetically induced transparency in superconductive quantum circuits , _ phys .
* 93 * , 087003 ( 2004 ) .
z. dutton , k.v.r.m .
murali , w.d .
oliver , and t.p .
orlando , electromagnetically induced transparency in superconducting quantum circuits : effects of decoherence , tunneling , and multilevel crosstalk , _ phys .
b _ * 73 * , 104516 ( 2006 ) .
leek , j.m .
fink , a. blais , r. bianchetti , m. goppl , j.m .
gambetta , d.i .
schuster , l. frunzio , r.j .
schoelkopf , and a. wallraff , observation of berry s phase in a solid - state qubit , _ science _ * 318 * , 1889 - 1892 ( 2007 ) .
a. wallraff , d. i. schuster , a. blais , l. frunzio , r .- s .
huang , j. majer , s. kumar , s. m. girvin and r. j. schoelkopf , strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics , _ nature _ * 431 * , 162 - 167 ( 2004 ) .
j.  johansson , s. saito , t. meno , h. nakano , m. ueda , h. tanaka , k. semba , and h. takayanagi , vacuum rabi oscillations in a macroscopic superconducting qubit lc oscillator system , _ phys .
* 96 * , 127006 ( 2006 ) .
f. deppe , m. mariantoni , e.p .
menzel , a. marx , s. saito , k. kakuyanagi , h. tanaka , t. menon , k. semba , h. takayanagi , e. solano , and r. gross two - photon probe of the jaynes - cummings model and controlled symmetry breaking in circuit qed , _ nature phys .
_ * 4 * , 686 - 691 ( 2008 ) .
fink , m. goppl , m. baur , r. bianchetti , p.j .
leek , a. blais , and a. wallraff , climbing the jaynes - cummings ladder and observing its root n nonlinearity in a cavity qed system , _ nature _ * 454 * , 315 - 318 ( 2008 ) .
a. fragner , m. goppl , j.m .
fink , m. baur , r. bianchetti , p.j .
leek , a. blais , and a. wallraff , resolving vacuum fluctuations in an electrical circuit by measuring the lamb shift , _ science _ * 322 * , 1357 - 1360 ( 2008 ) .
m. hofheinz , e.m .
weig , m. ansmann , r.c .
bialczak , e. lucero , m. neeley , a.d .
oconnell , h. wang , j.m .
martinis , and a.n .
cleland , generation of fock states in a superconducting quantum circuit , _ nature _ * 454 * , 310 - 314 ( 2008 ) .
van der ploeg , a. izmalkov , a. maassen van den brink , u. hbner , m. grajcar , e. ilichev , h .-
meyer , and a.m. zagoskin , controllable coupling of superconducting flux qubits , _ phys .
* 98 * , 057004 ( 2007 ) .
j. majer , j.m .
chow , j.m .
gambetta , j. koch , b.r .
johnson , j.a .
schreier , l. frunzio , d.i .
schuster , a.a .
houck , a. wallraff , a. blais , m.h .
devoret , s.m .
girvin , and r.j .
schoelkopf coupling superconducting qubits via a cavity bus , _ nature _ * 449 * , 443 - 447 ( 2007 ) .
r.  mcdermott , r. w. simmonds , m. steffen , k. b. cooper , k. cicak , k. d. osborn , s. oh , d. p. pappas , j. m. martinis , simultaneous state measurement of coupled josephson phase qubits , _ science _ * 307 * , 1299 - 1302 ( 2005 ) .
m. steffen , a. ansmann , r.c .
bialczak , n. katz , e. lucero , r.  mcdermott , m. neeley , e.m .
weig , a.n .
cleland , and j. m. martinis , measurement of the entanglement of two superconducting qubits via state tomography , _ science _ * 313 * , 1423 - 1425 ( 2006 ) .
m. steffen , a. ansmann , r.  mcdermott , n. katz , r.c .
bialczak , e. lucero , m. neeley , e.m .
weig , a.n .
cleland , and j. m. martinis , state tomography of capacitively shunted phase qubits with high fidelity , _ phys .
lett . _ * 97 * , 050502 ( 2006 ) .
m. neeley , m. ansmann , r.c .
bialczak , m. hofheinz , n. katz , e. lucero , a. oconnell , h. wang , a.n .
cleland , and j. m. martinis , process tomography of quantum memory in a josephson - phase qubit coupled to a two - level state , _ nature phys . _ * 4 * , 523 - 526 ( 2008 ) .
n. katz , m. ansmann , r.c .
bialczak , e. lucero , r. mcdermott , m. neeley , m. steffen , e.m .
weig , a.n .
cleland , j.m .
martinis , a.n .
korotkov coherent state evolution in a superconducting qubit from partial - collapse measurement , _ science _ * 312 * , 1498 - 1500 ( 2006 )                                                                s. yoakum , l. sirko , p. m. koch , stckelberg oscillations in the multiphoton excitation of helium rydberg atoms : observation with a pulse of coherent , field and suppression by additive noise _ phys .
* 69 * , 1919 - 1922 ( 1992 ) .
f. lang , p.v.d .
straten , b. brandsttter , g. thalhammer , k. winkler , p.s .
julienne , r. grimm , and j. hecker denschlag , cruising through molecular bound - state manifolds with radiofrequency , _ nature physics _ * 4 * , 223 - 226 ( 2008 ) .
p.  k. tien and j. p. gordon , multiphoton process observed in the interaction of microwave fields with the tunneling between superconductor films _ phys .
* 129 * , 647 ( 1963 ) .