Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: our progress in understanding the large scale make - up of the universe has been propelled rapidly in the last 15 years by extremely detailed measurements of the cosmic microwave background ( cmb ) . current measurements of the cmb temperature anisotropies have allowed us to place stringent constraints on some of our universe s key parameters and are proving to be a driving force in fundamental physics ( see e.g.  @xcite ) . cosmic - variance places limits on the amount of information which can be extracted from the temperature anisotropies alone , however . the _ planck _  satellite , for example , is expected to reach this limit over the entire primary temperature anisotropy power spectrum @xcite . fortunately , the polarization of the cmb anisotropies provides us with a wealth of additional cosmological information which we have only just begun to probe . the benefits of cmb polarization data are numerous . perhaps most excitingly , polarization data provide a chance to detect primordial gravitational waves , through the measurement of b - mode polarization @xcite . the rapid improvements in cmb polarization data provided by experiments like pique @xcite , dasi @xcite , capmap @xcite , boomerang @xcite , cbi @xcite , maxipol @xcite , _ wmap _ @xcite , and upcoming experiments like bicep , clover @xcite , quad , and the _ planck _  satellite , ensure that polarization will become one of the most important tools we have for furthering our understanding of the universe in the years to come .    unfortunately , the polarization signal is much more difficult to measure than the temperature anisotropies . the magnitude of the cmb polarization anisotropies is at most 10 per cent that of the temperature anisotropies , making their detection alone a challenging task . this difficulty is being gradually overcome with improvements in detector technology . another problem with the cmb polarization , however , is that its signal is dominated over most of the sky by foregrounds . in particular , over the wavelength range where the cmb is typically measured , synchrotron and non - spherical dust emission from within our galaxy contaminate uncleaned polarization maps . it is not yet clear how well the problem of foreground contamination can be solved . its resolution will require both a detailed knowledge of the structure and magnitude of the foregrounds , and an understanding of how to remove them from the polarization data without introducing unwanted systematics . understandably , the issue of polarization foregrounds has received much attention in the literature . several papers have been written on the construction of templates for the emission @xcite , and a number of methods have been proposed for performing the foreground removal @xcite . testing the efficacy of these methods is a crucial step in the analysis pipeline , to compare the abilities of the various removal techniques and ultimately to provide confidence in the accuracy of cmb polarization data . tests for the existence of residual polarized foreground contamination have received less attention in the literature , however . the test used in the _ wmap _  team s analysis , where foregrounds are cleaned using a template method in pixel space , is a simple @xmath1 statistic @xcite . unfortunately , this statistic can not be used to make comparisons between different cleaning methods . the only currently proposed ` general ' statistic which is capable of this application is based on the bipolar power spectrum ( bips ) and has been developed by @xcite .    in this paper , we introduce and analyse another general method of testing for residual foreground contamination in cmb polarization maps , the @xmath0 statistic of @xcite . it possesses the virtues of having a simple interpretation as a measure of directionality in a cmb map and of being extremely rapid to compute . we argue that this statistic is well suited to the detection of foreground contamination , as both the synchrotron and dust foreground components have strong overall directionality on large scales . we then apply the @xmath0 statistic to the _ wmap _ 3-yr polarization maps to demonstrate its usefulness for real data . the @xmath2 statistic was first presented by @xcite ; the @xmath2 value for a sky map pixelized into @xmath3 pixels is defined as @xmath4 where @xmath5 runs over the unit sphere , and with the function @xmath6 given by @xmath7 here @xmath8 are weights chosen to compensate for noise and the effects of a cut sky , and @xmath9 is a local directionality vector calculated for the pixel @xmath10 . the @xmath6 value is a measure of how much ( or how little ) the @xmath9 of a map tend to point in a given direction . the @xmath2 statistic tells us how extreme these values are . once @xmath2 has been calculated for real sky data , we can compare its value to that found for simulated realizations of the cmb . calculating @xmath2 for a large number of these simulations gives us a histogram of expected values , and excess directionality in a cmb data - set appears as a value of @xmath2 which is an outlier of this distribution .    for sky maps of temperature , a reasonable quantity to use for @xmath9 is the local temperature gradient @xmath11 . for polarization maps , on the other hand , the polarization psuedo - vectors themselves are obvious candidates . we define @xmath9 to be a vector whose magnitude and direction are those of the polarization at point @xmath10 . to express @xmath9 explicitly in terms of the stokes parameters @xmath12 and @xmath13 , recall that the polarization magnitude is @xmath14 while the polarization direction lies in the tangent plane to the celestial sphere at @xmath10 and makes an angle with the meridian of @xmath15 measuring the angle from the meridian is a _ convention @xcite which we adopt here , although it is not universally followed . we define the vector @xmath9 to have magnitude @xmath16 and direction given by @xmath17 . because the polarization is a spin-2 quantity , it is represented by headless ` pseudo - vectors ' rather than vectors , so @xmath17 can always be rotated by @xmath18 without altering the polarization . however , due to the quadratic definition of @xmath6 , this ambiguity does not affect @xmath0 . moreover , because the function @xmath19 smoothly averages together polarization information over the whole sky , @xmath2 is only sensitive to large scale ( @xmath20 ) directionality .    the weights @xmath21 must be chosen to satisfy the null hypothesis that for isotropically distributed @xmath9 no preferred direction will be found on average , even when parts of the sky have been removed . we can see how these weights may be obtained by casting @xmath6 in matrix form as @xmath22 where @xmath23 is a @xmath24 matrix defined by @xmath25    in order for no preferred direction to be found , we must have @xmath6 independent of @xmath5 ( on average ) . from equation  [ eq.fmatrixdef ] it can be seen that this constraint may be realized by choosing suitably normalized @xmath21 such that @xmath26 is equal to the identity matrix for isotropically distributed @xmath9 . this criterion ensures that @xmath6 will be approximately constant over the sky when the underlying data lack a preferred direction , although it does not guarantee that fluctuations of @xmath19 about its ensemble average will be isotropic . as a result , @xmath2 may still pick out preferred directions arising from the noise structure or sky coverage , even when the underlying signal is isotropic . we find that such departures from anisotropy are weak in practice ; this point is discussed further in section  [ section.wmappolarizationanalysis ] .    because @xmath23 is symmetric , requiring @xmath27 gives only 6 independent equations on the weights . to completely determine them for a large number of pixels , the additional constraint is imposed that @xmath28 be minimized , where @xmath29 is the expected value of @xmath30 given by pixel noise and cosmological signal . for a large signal - to - noise ratio on an isotropic sky this is equivalent to minimizing the variance of the weights . when noise is not negligible , minimizing @xmath28 rather than @xmath31 has the favourable effect of down - weighting noisy pixels . more explicit details on the calculation of the weights are given in @xcite . once the weights @xmath21 have been calculated , the simple quadratic definition of @xmath6 makes it possible to calculate @xmath2 for a given data - set in @xmath32 operations , where @xmath3 is the number of pixels . this can be done by noting that the values of @xmath5 which extremize @xmath6 are given by the eigenvectors of @xmath23 . @xmath0 follows immediately by taking the ratio of the largest and smallest eigenvalues . @xmath2 is one of the simplest statistics which can be proposed for the identification of statistical anisotropy in a cmb map . to compare it to the bips method of @xcite , we note that it is ultimately a less powerful statistic , as bips can be used to probe both large and small scales , whereas @xmath2 is only sensitive to large scale anisotropies . the advantage of @xmath0 , however , is that it can be calculated in @xmath32 operations , an unbeatable scaling . this speed makes @xmath2 an excellent statistic for ` first look ' determination of anisotropy in cmb maps , and the concomitant possible foreground contamination . the two main sources of foreground contamination in cmb polarization maps at the @xmath33ghz frequencies where the cmb is usually probed are galactic synchrotron and thermal dust emission . polarized synchrotron emission results from the acceleration of cosmic - ray electrons as they orbit in the galactic magnetic field . polarized dust emission , on the other hand , is the result of non - spherical dust grains which tend to align their long axes perpendicular to the field ; these dust grains preferentially emit thermal radiation polarized along their long axis @xcite . the polarization of both the synchrotron and thermal dust emission has its origins in the anisotropy introduced by the galactic magnetic field , and at cmb wavelengths ( where galactic faraday rotation is negligible ) these two components should both be polarized preferentially in the same direction . thus , we expect that they have a common and distinctive directionality on the sky , dictated by the large scale structure of the galactic magnetic field . to illustrate the signature of this directionality under application of the @xmath0 statistic , we analyse the low resolution k - band polarization map at @xmath34ghz and the dust polarization template of the _ wmap _  team , masking both with the standard polarization mask ( p06 ) to mimic the situation when analysing maps for foreground contamination . for the k - band map we find a preferred direction at a galactic latitude @xmath35 , within @xmath36 of the poles . for the dust emission template we find that the preferred direction lies at @xmath37 , @xmath38 shy of the poles . thus , the significant presence of synchrotron and dust foregrounds in a polarization map is expected to result in an outlying value of @xmath0 , with a preferred direction in the vicinity of the galactic poles . the release of the _ wmap _ 3-year data has provided us with a glimpse at the structure of cmb polarization . here we analyse the @xmath12 and @xmath13 polarization maps for foreground contamination using the @xmath2 statistic . we have chosen to examine the _ wmap _  low resolution 3rd year ( v2 ) polarization maps at healpix resolution 4 , corresponding to 3072 pixels on the entire sphere . these are the maps used for low @xmath39 analysis of the _ wmap _  data , corresponding to the angular scales which @xmath2 effectively probes . at this resolution it is also feasible to accurately simulate realizations of the _ wmap _  noise structure using the full noise covariance matrix , which becomes necessary for large angle analyses .    in all of our simulations , we generate cosmological signals given by the ` _ wmap _  + all @xmath40cdm ' fit parameters @xcite . the details of this choice have little effect on the outcome of our simulations . all of the polarization maps are strongly noise - dominated , and we find that histograms of @xmath2 for simulated data with and without cosmological signal are indistinguishable for all of the _ wmap _  bands ; this is illustrated in fig . [ fig.dhistforvpolwwocosmo ] for simulated v - band maps . thus , we believe that , with an accurate treatment of the maps noise properties , any excess directionality which we detect in the _ wmap _  data can be attributed to residual foregrounds , and not a cosmological signal .     is the value of @xmath2 calculated for the v - band data themselves . ]    for all of our simulations , we add correlated gaussian noise to each unmasked map pixel based on the _ full qu covariance matrices _ for the corresponding _ wmap _  frequency band , as provided by the _ wmap _  team . we generate this noise using the cholesky decomposition of the covariance matrix . this method accounts for qq , qu and uu correlations both within each pixel and with other pixels across the sky .    for uncleaned maps , the @xmath2 statistic produces the expected results . for each band , we generate and analyze @xmath41 simulated skies . an example of this is shown in fig . [ fig.dhistforvpolwwocosmo ] . we can use the number of standard deviations ( at which the real data lie from the histogram mean ) as an approximate measure of the significance with which we detect foreground contamination . this is illustrated in fig . [ fig.significancesforintrapixel ] , which shows the expected dependence on radiometer band as the synchrotron emission decreases with frequency and the dust contribution increases , leading to an apparent minimum of foreground contamination near the v - band . we can see that use of the standard p06 polar mask results in a greater level of contamination than when we use the more stringent p02 mask . for the p06 mask , all of the preferred directions which are found for the real data lie within @xmath42 of the galactic poles , and for the p02 mask the preferred directions all lie within @xmath43 of the poles , further demonstrating the foreground origin of the directionality excess . with p06 masking , we find that @xmath44 of our v - band simulations have values of @xmath2 less than the value calculated for the _ wmap _  data . with p02 masking , this result drops to @xmath45 . the corresponding significance figures of merit are @xmath46 and @xmath47 respectively .    ) symbols are for the standard p06 mask ( 27 per cent sky cut ) , while the ( @xmath48 ) symbols are for the more stringent p02 masking ( 50 per cent sky cut for our downgrading scheme ) . the curve roughly follows that expected from the frequency dependence of contamination by synchrotron and dust foreground components . ] calculated for the q- and v - band data themselves . ]    in order to assess the efficacy of the foreground removal performed by the _ wmap _ team we also analyse the low resolution q- and v - band foreground reduced maps . these have had dust and synchrotron templates projected out to remove foreground contamination . the results of these analyses are shown in fig . [ fig.qvforeredhists ] . from this figure it can be seen that the @xmath12 and @xmath49 maps have similar noise properties under analysis with @xmath2 , as the simulation histograms almost completely overlap . we find no evidence for residual foreground contamination in these maps , their @xmath2 statistics both lying well within the range of probable values , and with the maximal @xmath2 directions for the q- and v - band data being @xmath50 and @xmath51 from the pole , respectively . our significance figure of merit is @xmath52 for the q - band and @xmath53 for the v - band maps . @xmath54 of the simulated @xmath2 values lie below the map value for the q - band , and @xmath55 of the simulated data lie below the v - band value . thus , we find no evidence for foreground contamination in these cleaned maps . it is interesting to note that if we had generated noise which only accounted for the intra - pixel or diagonal components of the noise covariance matrix then we would have found strong evidence for residual contamination . in the case of the v - band foreground reduced data , for example , we would have found some evidence for foreground contamination , with over @xmath44 of simulated @xmath2 values falling below that calculated for real data in the cleaned v - band maps ( corresponding to a value of @xmath46 for our significance figure ) . accounting for _ inter_-pixel noise correlations broadens the distribution of simulated @xmath2 values . indeed , the detailed noise structure of the _ wmap _  data appears to be such that it introduces a degree of anisotropy into the maps . this can best be seen by looking at the histograms of preferred directions found in a series of 30,000 v - band map simulations , given in fig . [ fig.diagintrafullcovhists ] . differing methods of noise generation were utilized for three sets of simulations . ` diagonal ' noise was generated from only the @xmath56 and @xmath57 elements of the covariance matrix , intra - pixel noise included @xmath58 terms , and the full covariance matrix treatment simulated inter - pixel correlations as well . unlike the isotropic distribution of preferred directions which is found when accounting only for diagonal noise correlations , there is a tendency towards anisotropy for the intra - pixel and full covariance simulations , with a definite slight preference for directions in the vicinity of @xmath59 for the full covariance simulations . the statistics of these two dimensional histograms support this visual finding , with the diagonal , intrapixel , and full covariance simulation methods having variance / mean values of @xmath60 , @xmath61 and @xmath62 respectively . for an isotropic sky , we expect the histogram values to follow a multinomial distribution , and to have a variance / mean of @xmath63 . we believe that the excess variance which is seen in the diagonal and intra - pixel simulation histograms can be attributed to the slight inherent bias possessed by @xmath2 , as discussed in section  [ section.thedstatistic ] . the much higher excess variance seen in the inter - pixel simulations can be attributed to the sum of this bias and another , stronger bias given by the noise covariance structure of the map . this bias is independent of the specifics of the @xmath2 statistic , and most likely can be detected in the future by other polarization statistics sensitive to anisotropy . the dependence of our results on the precise treatment of noise underscores the general need for care when dealing with large angle polarization data . we have investigated the @xmath2 statistic of @xcite in the context of cmb polarization data , finding that it is an excellent general tool for quickly assessing the magnitude of foreground contamination in polarization maps . under analysis with @xmath2 , imperfectly removed foregrounds leave a tell - tale signature of excess directionality toward the galactic poles . an analysis of the _ wmap _  3-year polarization data supports these statements . we find the expected dependence of foreground contamination on channel frequency in all of the _ wmap _  low resolution maps , with the magnitude of this excess disappearing in the cleaned , ` foreground reduced ' maps .    in the future , a better understanding of the polarization foregrounds will allow us to develop more sophisticated statistics to precisely evaluate the possibility of foreground contamination in cleaned maps . given the current state of understanding , however , we feel that calculation of @xmath2 is a useful technique for this purpose , as it is one of the simplest possible statistics which captures the directional signature of the polarization foregrounds . it has the advantages of being simple and rapidly computable , and can be considered complementary to other statistics such as the bips @xcite . hence it should be one of the diagnostic tests ready to assess foreground contamination in data from the _ planck _  satellite , as well as other large - angle polarization experiments . the work carried out in this paper made use of the healpix @xcite package for pixelization . this research was supported by the natural sciences and engineering research council of canada and the canadian space agency . we also acknowledge the contributions of patrick plettner for some related earlier work . baccigalupi c. , 2003 , new astron . , 47 , 1127 barkats d. , et al . , 2005 , apj , 619 , l127 basak s. , hajian a. , souradeep t. , 2006 , phys . d , 74 , 021301 bouchet f.r . , prunet s. , sethi s.k . , 1999 , mnras , 302 , 663 bunn e.f . , scott d. , 2000 , mnras , 313 , 331 davis , l.j . , greenstein , j.l . 1951 , apj , 114 , 206 de oliveira - costa a. , tegmark m. , odell c. , keating b. , timbie p. , efstathiou g. , smoot g. , 2003 , phys . d68 , 083003 grski k.m . , hivon e. , banday a.j . , wandelt b.d . , hansen f.k . , reinecke m. , bartelman m. , 2005 , apj , 622 , 759 hansen f.k . , banday a.j . , eriksen h.k . , grski k.m . , lilje p.b . , 2006 , apj , 648 , 784 hedman m. , barkats d. , gundersen j.o . , staggs s.t . , winstein b. , 2001 , apj , 548 , l111 hu , w. , white m. , 1997 , newa , 2 , 323 kamionkowski m. , kosowsky a. , 1998 , phys . d , 57 , 685 kamionkowski m. , kosowsky a. , stebbins , a. , 1997 , phys . d , 55 , 7368 kovac j. , leitch e.m . , pryke c. , carlstrom j.e . , halverson n.w . , holzapfel w.l . , 2002 , nature , 420 , 772 masi s. , et al . , 2006 , a&a , 458 , 687 page l. , et al . , 2007 , apjs , 170 , 335 _ planck _ bluebook _the scientific programme_ , esa - sci(2006)1 ; astro - ph/0604069v1 readhead a.c.s . , et al . , 2004 , science , 306 , 836 scott d. , smoot , g.f . , 2006 , in the review of particle physics , yao w .- m . , et al . , jphysg , 33 , 1 seljak , u. , zaldarriaga , m. , 1997 , phys . lett . , 78 , 2054 slosar , a. , seljak , u. , makarov , a. , 2004 , phys . d , 69 , 123003 spergel d.n . , et al . , 2006 , apjs , 170 , 377 taylor , a. , 2006 , newa rev . , 50 , 993 wu j.h.p . , et al . , 2007 , apj , 665 , 55
Abstract:

Answer: polarization is the next frontier of cmb analysis , but its signal is dominated over much of the sky by foregrounds which must be carefully removed . to determine the efficacy of this cleaning it is necessary to have sensitive tests for residual foreground contamination in polarization sky maps . the dominant galactic foregrounds introduce a large - scale anisotropy on to the sky , so it makes sense to use a statistic sensitive to overall directionality for this purpose . here we adapt the rapidly computable @xmath0 statistic of bunn and scott to polarization data , and demonstrate its utility as a foreground monitor by applying it to the low resolution _ wmap _  3-yr sky maps . with a thorough simulation of the maps noise properties , we find no evidence for contamination in the foreground cleaned sky maps . [ firstpage ]    cosmic microwave background methods : numerical  cosmology : observations  cosmology : theory  large - scale structure


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the correspondence between classical and quantum dynamics of anharmonic systems has gained significant attention in the past few years @xcite . a short laser pulse impinging on an atom or a molecule excites a superposition of several stationary states , and the resulting wave packet follows the orbit of the corresponding classical particle in the initial stage of the time evolution . however , the nonequidistant nature of the involved energy spectra causes peculiar quantum effects , broadening of the initially well localized wave packets , revivals and partial revivals @xcite . partial revivals are in close connection with the formation of schrdinger - cat states , which , in this context , are coherent superpositions of two spatially separated , well localized wave packets @xcite . phase space description of vibrational schrdinger - cat state formation using animated wigner functions can be found in @xcite . however , these highly nonclassical states are expected to be particularly sensitive to decoherence @xcite . the aim of this paper is to analyze the process of decoherence for the spontaneously formed schrdinger - cat states in the anharmonic potential . we consider the decoherence model which relies on the interaction of the quantum system ( _ s _ ) with its environment ( _ e _ ) that has many degrees of freedom . the dynamics of the environment together with the quantum system under investigation is unitary , i.e. , the density operator @xmath0 of the coupled @xmath1 system always represents a pure state . starting from an initially uncorrelated density operator , the interaction builds up _ s__e _ entanglement and the reduced density operator @xmath2 where @xmath3 means a trace over the environmental degrees of freedom , turns into a mixture . assuming negligibly short relaxation times in the environment , one obtains a markovian master equation which is a useful tool for the dynamical investigation of the environment induced decoherence . it allows for the determination of pointer states @xcite , that is , states that are favored by decoherence @xcite , and also for the calculation of the characteristic time of decoherence for different initial states @xcite .    in the following we introduce a master equation that takes into account the fact that in a general anharmonic system the relaxation rate of each energy eigenstate is different this master equation is applied to the case of wave packet motion in the morse potential that is often used to describe a vibrating diatomic molecule . considering the phase - space description of decoherence , we show how the phase portrait of the system reflects the damping of revivals in the expectation values of the position and momentum operators due to the effect of the environment . we also calculate and plot the time evolution of the wigner function corresponding to the reduced density operator of the morse system . the wigner function picture visualizes the fact that although our master equation reduces to the amplitude damping equation @xcite in the harmonic limit , the anharmonic effects lead to a decoherence scheme which is similar to the phase relaxation @xcite of the harmonic oscillator ( ho ) . it is found that the reduced density operator that arises due to the decoherence can be identified with a mixture of states that are well - localized in the phase space and equally distributed along the orbit of the corresponding classical particle . we illustrate the generality of this decoherence scheme by presenting the time evolution of an energy eigenstate as well . we also calculate the decoherence time for various initial wave packets . we show that decoherence is faster for wave packets that correspond to a classical particle with a phase space orbit of larger diameter . the total hamiltonian of an anharmonic system and its environment has the form @xmath4 the self - hamiltonian of the system is written as : @xmath5 where the spectrum @xmath6 is assumed to be nondegenerate and discrete , but not necessarily equidistant . @xmath7 denotes the ground state energy of the system , and @xmath8 whenever @xmath9 . the environment is represented by a set of harmonic oscillators @xmath10 we assume the following interaction hamiltonian @xmath11 where @xmath12 is an appropriate system operator that transforms each eigenstate of @xmath13 into a superposition of different eigenstates corresponding to _ lower _ energy values . this is the application of the rotating wave approximation ( rwa ) to an anharmonic , multilevel system . the operator @xmath14 is the hermitian conjugate of @xmath12 , and for the sake of simplicity the coupling constants @xmath15 were taken to be real . in the case of a vibrating diatomic molecule @xmath16 describes the modes of the free electromagnetic field , while the terms @xmath12 and @xmath17 in the interaction hamiltonian are related to the molecular dipole moment operator . we assume that the environment is in thermal equilibrium at a given temperature @xmath18 , and the environmental oscillators have continuous distribution with frequency dependent density of states , @xmath19 . starting from the von neumann equation for the total density operator , standard techniques @xcite lead to the following markovian master equation in the schrdinger picture @xmath20 -\big ( \mathcal{x}^{\dagger}\mathcal{x}_{e}\rho_{s } + \mathcal{x}\mathcal{x}_{a}^{\dagger}\rho_{s}\nonumber   \\ & -&\mathcal{x}_{a}^{\dagger}\rho_{s}\mathcal{x } -\mathcal{x}_{e}\rho_{s}\mathcal{x}^{\dagger } + \mathrm{h.c . \label{me}\end{aligned}\ ] ] the matrix elements of the operators appearing in the nonunitary terms are given by @xmath21 where @xmath22 , @xmath23 denotes the average number of quanta in the corresponding mode of the environment , and the subscripts @xmath24 and @xmath25 refer to emission and absorption , respectively .    as we can see , the matrix elements ( [ mexelements ] ) of the operators that induce the transitions depend on the bohr frequency of the involved transition , which is a genuine anharmonic feature . in the special case of the ho , when @xmath13 has equidistant spectrum , and @xmath12 is identified with the usual annihilation operator @xmath25 , both @xmath26 and @xmath27 are proportional to @xmath28 , and eq . ( [ me ] ) reduces to the amplitude damping master equation @xcite at a finite temperature .    in certain cases one can further simplify eq . ( [ me ] ) . when the environment induced relaxation rates are much lower than the relevant bohr frequencies , the system hamiltonian induces oscillations that are very fast even on the time scale of decoherence and vanish on the average . ignoring these fast oscillations we arrive at the interaction picture master equation @xmath29 that has already been obtained in refs . @xcite in order to treat the spontaneous emission of a multilevel atom . in eq . ( [ me2 ] ) , @xmath30 denotes a relaxation rate , that is the probability of the @xmath31 transition per unit time , while @xmath32 , where @xmath33    however , due to the elimination of the fast oscillations related to @xmath13 , eq . ( [ me2 ] ) is not suitable for investigating the wave packet motion and decoherence simultaneously , therefore we propose to use eq . ( [ me ] ) . on the other hand we note that eq . ( [ me2 ] ) radically reduces the computational costs of calculating the time evolution for long times , which might be necessary when the system - environment coupling is very weak .    supposing that our knowledge is limited to the populations @xmath34 , both eq . ( [ me ] ) and eq . ( [ me2 ] ) leads to the pauli type equation @xmath35 requiring the condition of detailed balance @xcite in eq . ( [ pauli ] ) leads to the steady - state thermal distribution at the temperature of the environment .    when a diatomic molecule is considered in the environment of the free electromagnetic field , the operators @xmath12 and @xmath17 in eq . ( [ meinterh ] ) gain a clear interpretation : in the eigenbasis of @xmath13 they are the upper and lower triangular parts of the molecular dipole moment operator , @xmath36 . we will assume that @xmath36 is linear @xcite , that is , proportional to the displacement @xmath37 of the center of mass of the diatomic system from the equilibrium position . in this case @xmath38 , that is , @xmath39 where matrix elements of @xmath37 can be calculated using the algebraic method summarized in @xcite , and @xmath40 is an overall , frequency independent coupling constant .    however , in order to get insight into the interplay between wave packet motion and decoherence , it is worth considering a stronger molecule - environment interaction than the electromagnetic field modes can provide . keeping the structure of eqs . ( [ morsex ] ) , this can be done by increasing the value of @xmath41 . here we present calculations with two different coupling constants , @xmath42 and @xmath43 which are chosen so that at zero temperature @xmath44 and @xmath45 for @xmath46 and @xmath43 , respectively . this model allows for the numerical integration of the master equation ( [ me ] ) ( that provides more details of the dynamics than eq . ( [ me2 ] ) ) in a time interval that is long enough to identify the effects of decoherence . these effects can be summarized in a decoherence scheme ( see sec .  [ wignersec ] ) that has a clear physical interpretation , and which is valid also in the weak molecule - environment interaction , when ( [ me2 ] ) is more efficient to calculate the time evolution .    in the following we shall apply our master equation ( [ me ] ) to the case of the morse hamiltonian @xcite that is often used to describe a vibrating diatomic molecule . this hamiltonian has the dimensionless form @xmath47 , \label{ham}\ ] ] where the shape parameter @xmath48 is related to the dissociation energy @xmath49 , the reduced mass of the molecule @xmath50 , and the range parameter of the potential @xmath51 via @xmath52 the dimensionless displacement and momentum operators obey the canonical commutation relation @xmath53=i$ ] . the hamiltonian ( [ ham ] ) sustains @xmath54 + 1 $ ] normalizable eigenstates ( bound states ) , corresponding to the eigenvalues @xmath55 @xmath56 $ ] , where @xmath54 $ ] denotes the largest integer that is smaller than @xmath48 . the continuum above @xmath57}$ ] corresponds to the dissociated molecule with positive energies . for the sake of definiteness we have chosen the no molecule as our model , where @xmath58 . the initial wave packets of our analysis will be morse coherent states @xcite , @xmath59 , which are localized on the phase space around the point @xmath60 . the wigner function of a representative morse coherent state is shown in fig . [ wig1 ] a ) . these states can be prepared by an appropriate electromagnetic pulse that drives the vibrational state of the molecule starting from the ground state into an approximate coherent state . an example can be found in @xcite , where the effect of an external sinusoidal field is considered . let us note that the construction given in @xcite would allow us to use arbitrary initial states , but for our current purpose it suffices to consider states @xmath59 with negligible dissociation probability , i.e. , coherent states that practically can be expanded in terms of the bound states @xmath61 , @xmath62 $ ] . starting from @xmath63 as initial states , the qualitative behavior of the expectation value @xmath64 draws the limit of _ small _ oscillations . in the absence of environmental coupling ( i.e. , @xmath65 ) , for @xmath66 , @xmath67 ( as well as @xmath68 ) exhibits sinusoidal oscillations . [ tl][tl]@xmath69 , with @xmath70 . the parameters are @xmath46 , @xmath71 , and @xmath72 is the period of the small oscillations in the potential . the initial point @xmath73 , @xmath74 together with the starting direction is also indicated . [ portrait ] , title="fig:",width=321 ]    for larger initial displacements from the equilibrium position , the anharmonic effects become apparent : the amplitude of the oscillations decreases almost to zero , then faster oscillations with small amplitude appear but later we re - obtain almost exactly @xmath75 ( and @xmath76 as well ) , and the whole process starts again @xcite . without the influence of the environment the main attributes of @xmath67 and @xmath68 can be explained by referring to the various bohr frequencies that determine their time dependence : dephasing of these frequencies leads to the collapse of the expectation value , and we observe revival when they rephase again .    for the initial state of @xmath77 , with @xmath70 , the original phase of the eigenstates is restored @xcite around the full revival time @xmath78 , where @xmath72 is the period of the small oscillations in the potential . at @xmath79 and @xmath80 half and quarter revivals @xcite can be observed . [ portrait ] shows the damping of the revivals both in @xmath67 and @xmath68 when interaction with the environment is turned on . note that the phase portrait of the corresponding classical particle would be a helix with monotonically decreasing diameter , revivals are of quantum nature . however , fig . [ portrait ] does not provide a complete description of the time evolution in the phase space , this can be given by using wigner functions , see sec . [ wignersec ] . our master equation ( [ me ] ) describes decoherence as well as dissipation . however , the time scale of these processes is generally very different , and one can distinguish the stages of the time evolution that are dominated either by decoherence or dissipation @xcite . [ tl][tl][0.8]@xmath81 [ bl][bl][0.8]@xmath82 $ ] [ br][bl]@xmath83 ) . the coupling parameter ( see sec .  [ mastersec ] ) is @xmath46 and @xmath84 . the initial state was @xmath77 , with @xmath85 . [ dtdef ] , title="fig:",width=321 ]    in fig . [ dtdef ] an example is depicted showing how the method of time scale separation works . we have calculated the entropy @xmath86,\ ] ] as well as the quantity @xmath87 $ ] , which measures the purity of the reduced density operator . note that the @xmath88 operation without subscript refers to the trace in the system s hilbert space . decoherence time @xmath89 is defined as the time instant that divides the time axis into two parts where the character of the physical process is clearly different . initially both @xmath81 and @xmath90 $ ] change rapidly but having passed @xmath89 ( emphasized by a vertical line in fig . [ dtdef ] ) , the moduli of their derivative significantly decrease . after @xmath89 the entropy and the purity change on the time scale which is characteristic of the dissipation of the system s energy during the whole process @xcite . in other words , decoherence dominated time evolution turns into dissipation dominated dynamics around @xmath89 . in the next section we shall determine the density operators into which the process of decoherence drives the system . in connection with these results we have verified that the states around the decoherence time do not change appreciably in a time interval that covers the possible errors in determining @xmath89 . an interesting question is the dependence of the decoherence time on the initial state of the time evolution . we calculated @xmath89 as a function of the initial displacement for the case of displaced ground states ( that is , coherent states with zero momentum , @xmath91 ) as initial states . it was found that for all values of @xmath41 and @xmath18 , the decoherence time is longer for smaller initial displacements . additionally , for fixed @xmath41 and @xmath18 the function @xmath92 can be well approximated by an exponential curve @xmath93 . e. g. , for @xmath46 , @xmath84 and @xmath94 the parameters take the values @xmath95 and @xmath96 . it is known @xcite that quarter revivals in an anharmonic potential lead to the formation of schrdinger - cat states , i.e. , states that are superpositions of two distinct states localized in space @xcite as well as in momentum @xcite . on the other hand , smaller initial displacements correspond to classical phase space orbits with smaller diameter . consequently the quantum interference related to nonclassical states that are formed during the course of time cover a smaller area in the phase - space in this case . this means that our result is a manifestation of the general feature of decoherence that increasing the `` parameter of nonclassicality '' , which is the diameter of the corresponding classical orbit in our case , causes faster decoherence @xcite . a similar result was found in @xcite for the case of decoherence in a system of two - level atoms @xcite . in order to visualize the time evolution of the reduced density matrix of the morse system we have chosen the wigner function picture , which represents @xmath97 as a function over the classical phase space @xmath98 this description allows us to investigate the correspondence between classical and quantum dynamics .    first we consider the ideal case without environment . then , in the initial stage of the time evolution , the positive hill corresponding to the wave packet @xmath59 follows the orbit of the classical particle that has started from @xmath99 at @xmath100 . however , due to the uncertainty relation , the wigner function as a quasiprobability distribution has a finite width , and this fact combined with the form of the morse potential implies the stretching of the wigner function along the classical orbit in the course of time . ( see ref . @xcite for similar results with the husimi @xmath101 function . ) after a certain time the increasingly broadened wave packet becomes able to interfere with itself , and around the quarter revival time one can observe two positive hills chasing each other at the opposite sides of the classical orbit . the strong oscillations of @xmath102 between the hills represent the quantum correlation of the constituents of this molecular schrdinger - cat state @xcite . later on the initial wigner function is restored almost exactly and schrdinger - cat state formation starts again . detailed wigner function description of these processes that are related to the free time evolution can be found in @xcite .    in the case when environmental effects are present , we found that decoherence follows a general scheme . a representative series of wigner functions is shown in fig .  [ wig1 ] the snapshots correspond to the initial state and time instants when the first and third schrdinger - cat state formation would occur in the absence of the environment . consequently , the wigner function in fig .  [ wig1 ] b ) corresponds almost to a schrdinger - cat state , but this state is already a mixture . however , there are still negative parts of the function in between the positive hills centered at @xmath103 and @xmath104 . the `` ridge '' that connects these hills along the classical orbit is absent in a pure schrdinger - cat state . later on this ridge becomes more and more pronounced and at the decoherence time we arrive at the positive wigner function of fig . [ wig1 ] c ) and d ) . according to the contour plot fig . [ wig1 ] d ) , the highest values of this function trace out the phase space orbit of the corresponding classical particle . that is , @xmath105 , the reduced density matrix that arises as a result of decoherence , can be interpreted as a mixture of localized states that are equally distributed along the orbit of the corresponding classical phase space orbit . it is worth comparing this result with the case of the ho , when the master equation ( [ me ] ) reduces to the amplitude damping equation @xcite , see sec . [ mastersec ] . it is known that harmonic oscillator coherent states are robust against the decoherence described by the amplitude damping master equation ( as well as against the caldeira - leggett @xcite master equation @xcite ) , the initial superposition of coherent states turns into the statistical mixture of essentially the same states . this is a consequence of the facts that these states are eigenstates of the destruction operator @xmath25 , and the operators in the nonunitary terms of eq . ( [ me ] ) are proportional to @xmath25 and @xmath106 in the harmonic case . none of these statements can be transferred to the anharmonic system , where the morse coherent states do not remain localized during the course of time , even without environment . therefore the scheme of the decoherence is qualitatively different for the harmonic and anharmonic oscillators : our results in the anharmonic system are similar to the phase relaxation in the harmonic case @xcite , where the energy of the system remains unchanged , but the phase information is completely destroyed . we note that a similar result was obtained in ref . @xcite , where the rotational degrees of freedom were considered as a reservoir for the harmonic vibration of hot alkaline dimers . our decoherence scheme is universal to a large extent . in the investigated domain of the coupling constants @xmath107 and temperatures ranging from @xmath108 to @xmath109 , it is found to be valid for all initial states , not only for coherent states . [ wig2 ] shows an example when the initial state is not a wave packet , it is the fifth bound state , corresponding to @xmath110 , which is very close to @xmath111 , so direct comparison with fig . [ wig1 ] is possible . as we can see , although the two wigner functions are initially obviously very different , they follow different routes ( that takes different times ) to the _ same _ state : fig . [ wig1 ] c ) and fig . [ wig2 ] c ) are practically identical . the final plot in fig .  [ wig2 ] indicates how the wigner function represents the long way to thermal equilibrium with the environment : the distribution becomes wider and the hole in the middle disappears . it is expected that the loss of phase information has observable consequences . according to the franck - condon principle , the absorption spectrum of a molecule around the frequency corresponding to an electronic transition between two electronic surfaces depends on the vibrational state . the time dependence of the spectrum should exhibit the differences between the pure state of an oscillating wave packet and the state @xmath105 and the thermal state . more sophisticated experimental methods based on the detection of fluorescence @xcite or fluorescence intensity fluctuations @xcite , surely have the capacity of observing the dephasing phenomenon considered in this paper . we investigated the decoherence of wave packets in the morse potential . the decoherence time for various initial states was calculated and it was found that the larger is the diameter of the phase space orbit described by a wave packet , the faster is the decoherence . we obtained a general decoherence scheme , which has a clear physical interpretation : the reduced density operator that is the result of the decoherence is a mixture of states localized along the corresponding classical phase space orbit .
Abstract:

Answer: the time evolution of anharmonic molecular wave packets is investigated under the influence of the environment consisting of harmonic oscillators . these oscillators represent photon or phonon modes and assumed to be in thermal equilibrium . our model explicitly incorporates the fact that in the case of a nonequidistant spectrum the rates of the environment induced transitions are different for each transition . the nonunitary time evolution is visualized by the aid of the wigner function related to the vibrational state of the molecule . the time scale of decoherence is much shorter than that of dissipation , and gives rise to states which are mixtures of localized states along the phase space orbit of the corresponding classical particle . this behavior is to a large extent independent of the coupling strength , the temperature of the environment and also of the initial state .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: experimental advances in generating strongly interacting photons opened the opportunity to explore various many - body quantum phenomena in a new context  @xcite . photonic systems are unique in the sense that it is naturally an open quantum system where the adding and destroying of photons can be accomplished by drive and dissipation in a controlled way . thus it is an ideal system to study many outstanding questions on open systems dynamics , dissipative phase transitions  @xcite and the effects of interactions in a dissipative environment . the simplest model to study strongly interacting bosons on a lattice is the celebrated bose - hubbard model where atoms have on - site interactions and can hop across lattices  @xcite . there have been several recent proposals  @xcite on achieving this model with photons and open systems , such as with photons in coupled cavity arrays , superconducting circuit qed and polaritons . with added drive and dissipation , the bose - hubbard model does not exhibit superfluid or mott insulator phases but gives rise to mixed state nonequilibrium steady states and phases  @xcite , the nature and generation of which are not fully understood .    a starting point for understanding the complex dynamics of driven dissipative photonic cavity arrays can be a two - site model . such a two - site driven dissipative nonlinear model may be realized with systems such as two coupled photon cavities  @xcite , circuit qed systems  @xcite and two coupled micropillars  @xcite . this system is also referred to as a photonic molecule  @xcite , dimer or a double well . for the closed system , the physics of double well has been studied in great detail  @xcite , giving rise to phenomena such as the josephson effect , matter wave interference and self - trapped and symmetry breaking states , among others . for open systems , studies of two coupled cavities have appeared in several contexts in both theory  @xcite and experiment  @xcite . for an end - driven cavity , studies focused on the topics of unconventional photon blockade  @xcite and multi - stability  @xcite , among others . for symmetrically driven cavities , studies appeared on quantum correlations  @xcite , classical to quantum phase transitions for a jaynes - cummings dimer  @xcite , and symmetry breaking for incoherent drives  @xcite , among others .    in this article , we explore the physics of two coupled nonlinear cavities in a dissipative setting where both sites are driven coherently and symmetrically . we perform semiclassical and quantum analysis of the system investigating the complex interplay of many competing terms such as hopping , interaction , drive , dissipation and detuning . in a semiclassical treatment , we show that the nonequilibrium steady states have asymmetric number density in the two cavities in addition to the expected symmetry preserving states . these states are the driven - dissipative analog of the closed system double well symmetry breaking or self - trapped states  @xcite with a fundamental difference that these are not minimum energy states but long - time steady states resulting from the competition in drive and dissipation . these can be understood from the bistability of a single driven cavity ; when two cavities are coupled with small but nonzero tunneling , the low density and high density bistable branches of a single cavity are hybridized to form two symmetry breaking steady states with unequal photon occupations in the two cavities . beyond a critical coupling and drive , the symmetry breaking states do not form . we analyze the occurrence and stability of the semiclassical solutions , finding that there can be up to nine solutions , a maximum of four of which are stable containing pairs of symmetry breaking and symmetry preserving states in a multistable region . we delineate a phase diagram for the symmetry broken states in the parameter space of drive and tunneling . and inter - cavity mode overlap contributes to tunneling with strength @xmath0 . the system is coupled to a markovian bath and photons can decay or leak out with rate @xmath1 . coherent pumping ( @xmath2 ) replenishes the photons . we treat the case when drive and dissipation are the same for both cavities . two coupled cavities are the simplest setting for an array of coupled cavities where the interplay of dissipation , drive and interaction plays a crucial role that are being investigated in a wider context.,scaledwidth=35.0% ]    we then study the system by solving the full quantum mechanical master equation and using the method of quantum trajectories  @xcite . in a full quantum treatment , when quantum fluctuations are taken into account , the symmetry breaking states are no longer seen , similar to the case of single - cavity bistable states  @xcite . however , in quantum trajectory simulations of the dynamics , we show that quantum jump statistics of number differences reveal the presence of underlying semiclassical bistability indicating the presence of symmetry breaking states . finally , we present a set of analytical expressions for the steady state correlation functions using the complex p - representation  @xcite expressing the master equation in the form of fokker - planck equation . the solutions work for small tunnel couplings ; we discuss the regime of validity of the solutions in the parameter space . the article is organized as follows . in section i , we introduce the model . in section ii , we present semiclassical analysis and analyze the driven - dissipative symmetry breaking states . in section iii , we present quantum trajectory analysis and full quantum treatment of the master equation . in section iv , we derive analytical solutions for the steady state correlation functions , and summarize our results in section v. we consider two cavities with kerr nonlinearity  @xcite that are coupled by tunneling  @xcite . the cavities are coherently driven in a dissipative setting , with both drive and dissipation acting equally on both sites . [ fig : schematic ] shows a schematic of the set up . the system can be described by the following hamiltonian @xmath3 where @xmath4 and @xmath5 are the annihilation operators for the two cavities labelled @xmath6 and @xmath7 . here @xmath0 is the inter - cavity tunneling strength , @xmath8 is the anharmonicity or the nonlinear kerr - type interaction strength , @xmath2 is the driving field for cavity 1 and 2 , and @xmath9 is the detuning , after the rotating wave approximation for the drive term . @xmath10 and @xmath11 are the cavity resonance frequency and driving frequency , respectively .    . the four steady states are reached with four different initial values of @xmath12.,scaledwidth=48.0% ]    for our open system , there can be photon losses induced by spontaneous decay or cavity leakage . in the approximation that the system is weakly coupled to a markovian bath , the dynamics of the density matrix can be modeled by a quantum master equation in the lindblad form @xmath13+\gamma [ 2 ( { \hat a}_1 { \hat \rho } { \hat a}_1^{\dagger}+{\hat a}_2 { \hat \rho } { \hat a}_2^{\dagger } ) \nonumber \\ & & -({\hat a}_1^{\dagger } { \hat a}_1+{\hat a}_2^{\dagger}{\hat a}_2 ) { \hat \rho } - { \hat \rho}({\hat a}_1^{\dagger } { \hat a}_1+{\hat a}_2^{\dagger } { \hat a}_2)],\end{aligned}\ ] ] where @xmath1 is the dissipation rate in each cavity . we set @xmath14 throughout . ignoring quantum fluctuations , the mean field amplitudes for the field operators are @xmath15 . we assume @xmath16 and @xmath17 , where @xmath18 are occupation numbers and @xmath19 are phases for the two cavities , and @xmath20 . in the semiclassical approximation the correlation functions factorize , and the equations of motion become @xmath21 here @xmath22 and @xmath23 . we find the steady states of these equations by solving the four coupled differential equations eq . [ eq : eqofmotion ] , examining the long - time dynamics for different sets of initial conditions . for most values of parameters we get a single steady state . however , for the value of parameters @xmath24 , @xmath25 , we get four different steady states as depicted in fig .  [ fig : dynamics ] showing their long time dynamics . we take @xmath25 for all our calculations in this article , essentially giving other parameters in units of dissipation . two of the steady states have equal number of photon occupations in the two cavities as in figs .  [ fig : dynamics](a ) and ( b ) , one with low occupations and the other with higher occupations . in addition to these , we get steady states where the photon occupations are different in the two cavities as shown in figs .  [ fig : dynamics](c ) and ( d ) . this pair of states are asymmetric , mirror images of each other , and localized more in one of the cavities . these states are the driven - dissipative analog of the symmetry breaking states of nonlinear josephson coupled junctions as observed in bec and with photons  @xcite in a closed system . we call them _ driven - dissipative symmetry breaking states_. unlike the single cavity semiclassical solutions , the phases here are important as the relative phase between the two cavities can not be gauged away . for the symmetry preserving states we find @xmath26 , and for the symmetry breaking states @xmath27 . another method we employ to find the steady state solutions is to set @xmath28 and find roots of the four coupled equations : @xmath29 here we have four equations with four unknowns @xmath30 and @xmath31 that are real . we assume @xmath32 and solve for the variables @xmath33 which is then transformed back into numbers and phases of the two cavities as @xmath34 . the solution set contains both real and complex solutions . there is only one real solution for most of the parameter values ; this is the regime where there is no bistability . for certain parameter values , there can be more than one real solutions , up to a maximum of nine and a maximum of four of which are stable . the stable solutions that we obtain with this method match with what we obtain from the previous method as steady states of long - term dynamics starting with different initial conditions .    ) plane for @xmath35 and different values of @xmath36 , and @xmath37 respectively . the area labelled 1 has one symmetry preserving steady state , 2 has two symmetry preserving states in a bistable region while 4 has two symmetry breaking and two symmetry preserving states . ( b ) shows @xmath38 and @xmath39 as a function of @xmath0 for a fixed @xmath40 . ( c ) shows @xmath38 and @xmath39 as a function of @xmath2 for a fixed @xmath41 . these depict how the symmetry breaking states appear and disappear as a function of tunnelling and drive . we see that there is a critical @xmath0 value beyond which the symmetry breaking states do not appear.,scaledwidth=42.0% ]    we find stability by introducing small fluctuations to the steady state , obtaining linearized equations for the fluctuations and examining their eigenvalues . introducing small fluctuations about the steady state @xmath42 , @xmath43 , @xmath44 here the stability matrix @xmath45 is following : @xmath46 where @xmath47 . when the real part of the eigenvalues of @xmath45 are positive the solution is stable . furthermore , we look at the determinant and trace , as done in ref . @xcite , for which the hurwitz criterion for stability requires that the trace and determinant of this matrix is nonzero and positive for stable eigenvalues .        we can further analyze the steady state equations eq . [ eq : eqsteadystate ] to get two state equations @xmath48 eqs . [ eq : eqdrive ] capture the parameter dependence for both the symmetry preserving and symmetry breaking semiclassical steady states . for symmetry preserving states , when @xmath49 , @xmath26 and @xmath50 , this simplifies to one equation @xmath51 further taking the limit @xmath52 reproduces the single - cavity result of drummond and walls  @xcite @xmath53 eq . [ eq : singlecavity ] is the bistability state equation for a single cavity . the extension of this to a coupled cavity gives us the state equations eqs . [ eq : eqdrive ] and [ eq : coupledbistability ] , expressing the dependence of @xmath0 in going from a single - cavity bistability to the coupled cavity bistability ( eq .  [ eq : coupledbistability ] ) and multistability ( eq .  [ eq : eqdrive ] ) .    the symmetry breaking states are borne out of the bistable behavior of a single cavity . we can understand this by introducing an infinitesimal coupling to the two cavities as illustrated in fig  [ fig : sbstates ] . as soon as we turn on @xmath0 , two bistable branches split into four steady states whose number occupations suggest that they are made up of the @xmath54 and @xmath55 density bistable states of individual cavities in this fashion - @xmath56 , @xmath57 , @xmath58 , @xmath59 . the @xmath58 and @xmath59 states are the symmetry breaking states . for a single cavity , there are three solutions to the semiclassical equations ; two are stable and one is unstable . coupling the two cavities gives rise to a maximum of nine solutions in some regime ; out of which a maximum of four are stable . there are regions of one , two or four stable solutions . the regime of four stable solutions is the multistable region where two are symmetry preserving states and two are symmetry breaking states . [ fig : sbstates](c ) shows both stable and unstable branch of the solutions marked in blue and red respectively . a pair of symmetry breaking states never appear in isolation but always come with a pair of symmetric states .    in fig . [ fig : phasediag ] we analyze the steady state number and phase differences between the two cavities and delineate a phase diagram in the @xmath60 phase . [ fig : phasediag](b ) shows the number difference @xmath38 as a function of @xmath0 for a fixed @xmath2 , and we see that the self - trapped states appear as soon as we turn on @xmath0 and disappears for a critical @xmath0 . here the positive and negative values of @xmath38 represent the pair of symmetry breaking states while the symmetric states are on zero of the @xmath61-axis . similarly , fig . [ fig : phasediag](c ) shows the dependence on @xmath2 for a fixed @xmath0 . with many of such slices for fixed @xmath2 and @xmath0 , we can now draw a phase diagram for the presence of symmetry broken states in the tunneling - drive ( @xmath60 ) plane for @xmath35 and different interactions ( a ) @xmath36 and ( d ) @xmath62 . the shades differentiate three regions with one , two and four stable steady states . the region labelled @xmath6 has one stable symmetry preserving steady state . the region labelled @xmath7 is a bistable regime with two symmetry preserving states which continues to grow as @xmath0 and @xmath2 increases . the boundaries of the bistability region can be obtained from the analytic expression eq . [ eq : coupledbistability ] which basically gives the turning points for the coupled cavity bistability in a similar way as can be derived for a single cavity  @xcite . the boundaries of bistability region are shown in figs . [ fig : phasediag](a ) , ( d ) , having a near linear relationship in @xmath2 and @xmath0 . the region labelled 4 and shaded dark blue is the region of multistability with two symmetry breaking states in addition to two symmetry preserving states . the blue region shrinks and straightens upward as we increase the nonlinearity @xmath8 as shown in fig . [ fig : phasediag](d ) . the characteristics of this region are intimately connected to the whole parameter space of detuning , nonlinearity and that of bistability . we note that the symmetry breaking region always stays inside the bistability region , and therefore can not go beyond the bistable phase boundaries to the right , where the slope a higher slope as @xmath8 increases . as a function of @xmath0 , the phase diagram shows a reentrant behavior as we can see in fig . [ fig : phasediag](a ) near @xmath63 .    in our model we do not break any symmetry of the system externally . yet we get steady states which break the symmetry of the system and a stability diagram where there are multiple states for a single value of drive . for coherent pumping as we treat here , these driven dissipative symmetry breaking has not yet been observed . we note that the recent observation of spontaneous mirror symmetry breaking in coupled photonic - crystal nanolasers  @xcite used incoherent pumping for which the mechanism and results are different in that the symmetry breaking states are not related to bistability as in our case . in terms of the bifurcation properties , unlike the supercritical pitchfork bifurfactions in experiment  @xcite , we get here subcritical bifurcation  @xcite as associated with bistable behavior . pumping just one end cavity excites the system modes in a different way and gives rise to multistability  @xcite , which has been recently observed in photonic microcavities  @xcite . creating an external phase asymmetry such as in drive ( driving two sites with @xmath2 and @xmath64 ) or tunnelling ( using @xmath0 instead of @xmath65 in the bose - hubbard hamiltonian ) , also gives rise to symmetry broken states , whose criticality and entanglement has been investigated recently  @xcite . analysis from the perspectives of semiclassical discrete nonlinear schrdinger equation ( dnls ) for lattices and soliton physics have been reported in ref . to understand more features of the driven dissipative bose - hubbard dimer , here we analyze the problem quantum mechanically taking into account the quantum fluctuations and using two methods  first , we examine the dynamics by numerically solving the lindblad master equation , and second , we do a quantum trajectory or monte carlo wavefunction analysis  @xcite .    in fig . [ fig : trajectory](b ) we show the semiclassical multistability diagram which exhibits up to four steady states for a single value of the drive . we integrate the lindblad master equation eq .  2 to find the quantum steady states which we find to be unique for each value of the drive . occupation ( @xmath66 ) and normalized second - order correlator ( @xmath67 ) are shown in the fig . [ fig : trajectory](b ) as a function of drive , where @xmath68 . in the multistability region , @xmath67 exhibits a peak which is due to the presence of enhanced quantum fluctuations that cause the underlying semiclassical multistability . this is similar to the case of single cavity bistability  @xcite . however , the @xmath67 signatures only indicate that this is a region of bistability or multistability without giving us any clue about the presence of sysmmetry breaking states . to reveal this feature , we perform quantum trajectory analysis below . quantum trajectory method provides exact results for physical observables under the ensemble averaging of trajectories of wavefunctions . we specifically would like to find out how the quantum jumps occur in the multi - stable region and whether they reveal any signature of the underlying semiclassical symmetry braking solutions . first we pick a specific driving field @xmath40 to run the quantum trajectory simulations . simulation results are analyzed in fig . [ fig : trajectory ] . panel ( a ) shows a histogram of the quantum jumps as they happened as a function of the occupation numbers of one cavity plotted vertically along the semiclassical multistability diagram . the two peaks of the histogram coincide with the lower and upper branches of the multistability diagram as expected for bistability . now if we investigate the dynamics of @xmath69 and @xmath70 separately , a typical such evolution for a single trajectory is given in fig . [ fig : trajectory](c ) . this shows clearly that the photons spend most of the time fluctuating near the two bistable branches . panel ( d ) shows a magnified region where the two cavities have unequal populations during the evolution . a histogram of jumps in the variable of photon population difference ( @xmath71 ) for @xmath40 is shown in panel ( f ) . the distribution in ( f ) not only shows a single peak at @xmath72 , but also broad side peaks at @xmath73 . for comparison , [ fig : trajectory](e ) and ( g ) show the statistics for quantum jumps outside the multistable region at @xmath74 and @xmath75 respectively , showing single lorentzian peaks at @xmath72 . the highest peaks centered at @xmath72 correspond to the jumps related to symmetry preserving states , whereas , we interpret the broad side peaks of the distribution originating from the symmetry breaking states . in our model , the symmetry breaking states always coexist with the symmetric states . so the side peaks are always overshadowed by a prominent central peak . as for a quantum signature of driven - dissipative symmetry breaking states , we show here that a typical quantum trajectory reveals this feature . recent experiments  @xcite have used homodyne detection of photocurrents to detect the quantum trajectory underlying a bistability , and hence , the observation of quantum features of dimer symmetry breaking states should be within experimental reach . analytic solutions for steady states of driven - dissipative systems are very rare . for a single cavity with kerr nonlinearity , analytic solutions for the steady - state was presented by drummond and walls in ref . since then , there have been a number of instances where analytic solutions have been obtained for various scenarios in a single - cavity  @xcite . drummond and walls  @xcite used the complex p - representation to obtain a fokker - planck equation and present closed - form analytic expressions for the correlation functions . we apply a similar technique here for the two - coupled nonlinear cavities and present analytic solutions for the correlation functions expressed as a series . the solutions are applicable for small @xmath0 and deviate slowly from the exact numerical results as @xmath0 is increased . considering the fact that no analytic solutions are known for coupled cavities , we present our investigations here which may find applications in the study of weak coupling limit and in finding an improved solution . a short derivation is provided below .    in comparison to the single cavity case  @xcite , the additional term in the dimer hamiltonian is the hopping term , @xmath76 . after a rotating wave approximation , the contribution of this term to the master equation dynamics becomes : @xmath77 $ ] . now the full master equation with kerr nonlinearity is @xmath78-i\frac{u}{2}[{\hat a}_{1}^{\dagger 2}{\hat a}_{1}^{2}+{\hat a}_{2}^{\dagger 2 } { \hat a}_{2}^{2},{\hat \rho}]+[f { \hat a}_{1}^{\dagger}+f^ { * } { \hat a}_{1},{\hat \rho}]+[f { \hat a}_{2}^{\dagger}+f^ { * } { \hat a}_{2},{\hat \rho}]\nonumber\\   &   & + ij[{\hat a}_{1}^{\dagger}{\hat a}_{2}+{\hat a}_{2}^{\dagger}{\hat a}_{1},{\hat \rho}]+\gamma[2{\hat a}_{1}{\hat \rho } { \hat a}_{1}^{\dagger}-{\hat \rho } { \hat a}_{1}^{\dagger}{\hat a}_{1}-{\hat a}_{1}^{\dagger}{\hat a}_{1}{\hat \rho}]+\gamma [ 2{\hat a}_{2}{\hat \rho } { \hat a}_{2}^{\dagger}-{\hat \rho } { \hat a}_{2}^{\dagger}{\hat a}_{2}-{\hat a}_{2}^{\dagger}{\hat a}_{2}{\hat \rho}].\end{aligned}\ ] ] the contribution to the @xmath79 function dynamics from the hopping term is @xmath80 where @xmath81 . now the fokker - planck equation for the @xmath79 function is given as @xmath82 where the drift matrix and the diffusion vector respectively are @xmath83 in order to have a specific form of analytic solutions for the steady states of the fokker - planck equation , the potential conditions must be satisfied  @xcite . now let us check whether the potential conditions , @xmath84 , are satisfied for this multidimensional fokker - planck equation . here @xmath85 denoting @xmath86 as @xmath87 respectively , we find @xmath88 , and @xmath89 . these four potential conditions are satisfied . however , two potential conditions are not satisfied : @xmath90 and @xmath91 . here @xmath92 if @xmath93 and @xmath94 if @xmath95 . they are also approximately satisfied in the limit @xmath96 is small . we present the analytic solutions under this restricted condition . the steady - state @xmath79 function is @xmath97.\end{aligned}\ ] ] we can rewrite the steady state @xmath79 function @xmath98,\end{aligned}\ ] ] where @xmath99 , and @xmath100 . we now get the analytic expression for any order of correlations . the zeroth order correlator or the normalization integral is @xmath101 first order on - site correlation function ( occupation ) is @xmath102 and the second order on - site correlation function is @xmath103    , @xmath35 , @xmath41 . for these set of parameters , our analytic expressions give a very good match with the exact numerics for the first and second order correlation functions @xmath66 and @xmath104 shown in ( a ) . for small drives @xmath2 , the normalized second order correlation @xmath67 however deviates . please see the discussions in the texts for explanations . ( b ) shows comparisons for a fixed @xmath63 and varying coupling @xmath0 . we find that for small values of @xmath0 the solutions are comparable but deviates slowly as we go to higher @xmath0.,scaledwidth=42.0% ]    a comparison of the results of the analytic solutions and numerical solutions of the master equation is shown in fig . [ fig : compare ] . in fig . [ fig : compare](a ) , for the set of parameters used , it is clear that the numeric and analytic results match very well for occupations @xmath66 and second - order correlation @xmath105 . for small values of drive , however , the normalized second - order correlator @xmath106 shows a deviation . this is due to the fact that for this specific value of parameters , slight deviations in @xmath105 in the numerator and @xmath107 in the denominator enhances the discrepancies when @xmath66 is very small . the mismatch could also be due to the potential condition not satisfied in that region , as discussed earlier . in fig . [ fig : compare](b ) we show comparisons for increasing @xmath0 for a fixed @xmath2 , and find deviations that come from the potential conditions not being satisfied . for increasing coupling the solutions deviate more and more , being comparable for small values of @xmath0 .    unlike the single cavity case , the series here do not sum to a closed form hypergeometric function . instead the computation involves a series involving six variables leading to difficulties in the convergence of the series . due to the properties of the gamma function , sometimes the series converges fast . more specifically , @xmath108 increases very fast with @xmath109 when @xmath109 is positive . thus if we have the the detuning positive ( i.e. the real part of @xmath110 and @xmath111 positive ) , the summation can converge for indices ( @xmath112 ) up to around 10 . if we have the detuning negative , we may need to compute up to a bigger index to make the series converge . we have checked the results for convergence for indices up to 30 . for a smaller ratio of @xmath96 as in our example fig . [ fig : compare ] , convergence is faster . it is conceivable that we can use monte carlo sampling technique to perform this multidimensional sum . a more fundamental issue is the regime of validity of the analytic expressions originating from the potential conditions as discussed earlier : @xmath90 and @xmath91 . the potential conditions are approximately satisfied when @xmath96 is small . from our investigation we find that when @xmath96 is small , as shown in fig . [ fig : compare ] , the analytical results are comparable to the numerical . the potential conditions are also satisfied when @xmath93 , which involves both the occupation and phase in the two cavities to be equal . in the regime where quantum fluctuations are important such as in the bistability and the symmetry breaking region , @xmath113 and the potential conditions are more likely to be violated . analytical results are rare for a driven - dissipative system especially when going beyond a single cavity . the results presented here , although with a restricted regime of validity , are therefore important in our view . the analytic solutions presented here can be extended to multiple cavities . we investigated the physics of two coupled nonlinear cavities in a lossy setting where both sites are driven coherently and equally . we performed semiclassical , quantum and analytical analyses of the system . in a semiclassical treatment , we find that the nonequilibrium steady states can have asymmetric number density in the two cavities which appear in addition to the symmetry preserving states . these states are the driven - dissipative analog of the double well self - trapped or symmetry broken states . their appearance can be understood from the bistability of a single driven cavity ; when two cavities are coupled , the low density and high density branches of single cavity bistable states hybridize to form two symmetry breaking states with unequal photon occupations in the two cavities . we examined the properties and stability of the semiclassical solutions , finding that there can be up to nine solutions of which a maximum of four are stable , giving rise to a pair of symmetry preserving and a pair of symmetry breaking states . we presented a phase diagram for these states in the tunneling - drive space . we further studied the system using the method of quantum trajectories and by solving the full quantum mechanical master equation . in a full quantum treatment , when quantum fluctuations are taken into account , the coupled cavity bistable self - trapped states no longer appear , a case similar to that of single - cavity bistability . however , in a quantum trajectory analysis of the dynamics , we found that a histogram of quantum jumps in number differences reveal the presence of semiclassical bistability with strong indication of symmetry breaking states . finally , we presented analytical solutions for the steady state correlation functions using the complex p - representation and forming a fokker - planck equation . we pointed out the regime of validity and limitations of this analytic solution .    coupled cavity arrays are an exciting system to explore a host of important phenomena such as nonequilibrium dynamics , open system physics , strongly interacting photons and quantum many - body physics . we took the simplest coupled cavity model of a dimer and analyzed it using several different methods . the physics explored here and our predictions are within experimental reach in coupled cavity dimers , in terms of multistable states and single trajectory measurements . besides fundamental physics , bistability and dimers have applications in optical memories and quantum correlation devices . the insights we gained on semiclassical and quantum nature of photons for two coupled cavities can also be useful for an array . we thank howard carmichael , ryan wilson and christian ciuti for helpful discussions . we knowledge the support from afosr - muri , onr - yip , aro - muri , nsf - pfc at the jqi , and the sloan foundation . mh thanks the kavli institute of theoretical physics ( kitp ) at santa barbara for hospitality . kitp is supported by nsf phy11 - 25915 .
Abstract:

Answer: we investigate two coupled nonlinear cavities that are coherently driven in a dissipative environment . we perform semiclassical , numerical and analytical quantum studies of this dimer model when both cavities are symmetrically driven . in the semiclassical analysis , we find steady - state solutions with different photon occupations in two cavities . such states can be considered analogs of the closed system double well symmetry breaking states . we analyze the occurrence and properties of these localized states in the system parameter space and examine how the symmetry breaking states , in form of a bistable pair , are associated to the single cavity bistable behavior . in a full quantum calculation of the master equation dynamics that includes quantum fluctuations , the symmetry breaking states and bistability disappear due to the quantum fluctuations . in quantum trajectory picture , we observe enhanced quantum jumps and switching which indicate the presence of the underlying semiclassical symmetry breaking states . finally , we present a set of analytical solutions for the steady state correlation functions using the complex p - representation and discuss its regime of validity .